Setting the Solver Properties and Stopping Criteria
Choosing the correct time-step is an important step in setting up an aeroacoustic simulation.
You typically evaluate the following criteria, and choose the smallest value:
- Convective Courant Number: Implicit solvers are usually stable at maximum values in the range 10-100 locally, but the overall mean value should be around 1. For this case, a time-step of 2.5E-5 s gives a Convective Courant Number whose mean value is approximately 1.
- Maximum resolvable frequency: For example, to obtain 10 points in a wave at 3000 Hz, (around the point at which human hearing is most sensitive), the time-step size should be 10 x 3000 Hz = 30000 Hz, or 3.33E-5 s.
- Local Strouhal shedding: At the Reynolds number observed in this tutorial, 63800, a cylinder sheds vortices at a Strouhal number of approximately 0.22. This number corresponds to a shedding frequency of 500 Hz. Capturing this shedding frequency requires a time-step size of 2.0E-4 s.
The smallest time-step size is the one required by the Convective Courant Number criterion, 2.5E-5 s.
To set the time-step size:
- Select the 2.5E-5 s. node and set Time-Step to
- Set Temporal-Discretization to 2nd-order.
Now set suitable under-relaxation factors:
- Select the 0.9. node and set Under-Relaxation Factor to
- Select the 0.7. node and set Under-Relaxation Factor to
The case is run until 0.0533s and with 5 iterations per time-step:
-
Edit the
Stopping Criteria node and set the following properties:
Node Property Setting Maximum Inner Iterations Maximum Inner Iterations 5 Maximum Physical Time Maximum Physical Time 0.0533 s Maximum Steps Enabled Deactivated
As the simulation is unsteady, you
update the trigger on the monitors for the lift and drag coefficients:
- Within the Monitors node, multi-select the Drag Coefficient Monitor and Lift Coefficient Monitor nodes, and set Trigger to Time Step.
- Save the simulation.