SpectraSim Documentation
Parameter · 07 Run Settings

Maximum Courant number

simulation_control.transient.max_courant Step 07 · Run Settings

The adaptive time step of a transient run is chosen so that the largest Courant number in the mesh stays at this value. A larger value is faster per simulated second and less stable.

Last updated

Theory

The Courant number of a cell is the number of cells the flow crosses in one time step, equation (1). With adjustTimeStep yes, OpenFOAM computes the largest Courant number of the mesh after each step and scales the next step towards ComaxCo_{max}, growing by at most 20 % per step and never beyond Maximum time step, equation (2) (OpenFOAM user guide, 2026).

The largest Courant number occurs in the smallest cells, the boundary-layer cells next to fast walls, so these cells set the time step of the whole run: about 1.6×10−51.6\times10^{-5} s at Courant 1 on a 60 mm pipe.

PIMPLE with several outer correctors per step (PIMPLE outer correctors) tolerates Courant numbers above 1. A transient run always uses first-order upwind convection for kk and ω\omega (k and ω convection scheme).

Equations

Co=Δt2VP∑f∣ϕf∣Co = \frac{\Delta t}{2V_P}\sum_f \lvert\phi_f\rvert
(1)
Δtn+1=min⁡ ⁣(min⁡ ⁣(λ, 1+0.1λ, 1.2)Δtn, Δtmax),λ=ComaxComaxn\Delta t^{n+1} = \min\!\left(\min\!\left(\lambda,\ 1 + 0.1\lambda,\ 1.2\right)\Delta t^{n},\ \Delta t_{max}\right), \qquad \lambda = \frac{Co_{max}}{Co^{n}_{max}}
(2)
wheremeaningunit
CoCoCourant number of cell P-
Δt\Delta tTime steps
VPV_PCell volumem³
ϕf\phi_fVolume flux through face fm³/s
ComaxnCo^{n}_{max}Largest Courant number of the mesh at step nn-
Δtmax\Delta t_{max}Maximum time steps

Default and rule

Default
2Measured
Rule

2.

Source

Validation runs of the program (residual-sensitivity study):

  • Courant 5 with one outer corrector, an earlier default, diverged on every transient case tried, a plain verification pipe included: ω\omega rose to about 106510^{65} and the time step collapsed to about 10−1310^{-13} s within a few hundred steps.
  • Courant 2-3.5 with 2 outer correctors and first-order kk, ω\omega convection was stable on the pipe and on the mufflers.
  • On case V2 (a two-inlet single-chamber muffler with internal inlet and outlet tubes) at Courant 3.5, the run was stable through at least 0.15 s of simulated time, with steps of about 10−510^{-5} s.

Equation (2) is OpenFOAM's time-step adjustment as used by the transient incompressible solver.

Valid range > 0; validated from 2 to 3.5 with 2 outer correctors

Values

Above 0. With 1 outer corrector, keep it at or below 1.

Above 3.5

Not validated. Courant 5 diverged on every case of the study with the earlier settings.

References 1

  1. OpenCFD Ltd. (2026). OpenFOAM user guide. openfoam.com. Web page