Authors: John Singleton and Rajesh Bhaskaran, Cornell University
Problem Specification
1. Pre-Analysis & Start-Up
2. Geometry
3. Mesh
4. Physics Setup
5. Numerical Solution
6. Numerical Results
7. Verification & Validation
Exercises
Comments
Verification & Validation
- Perform additional runs with timesteps of 0.02 and 1s. Comment on the stability of the scheme relating to the timestep by comparing the residuals.
- Compare the residuals and the vorticity figures for a mesh of 128x64 or another mesh of your choosing. Is there any noticeable differences between the results using different mesh sizes?
Effect of Uniform Initial Condition
In our simulation, we patched the cells behind the cylinder to have a starting velocity of .2 m/s so that we would get the vortex shedding earlier. This also allowed us to better see the oscilliation of the lift coefficient as flow time increased.
Effect of Timestep on Solution
In order to effectively capture the vortex shedding, a small enough timestep was needed. It was recommeded to have roughly 20 - 25 timesteps per cycle hence we used a timestep of .2. The resulting plot in the numerical results showing how the lift coefficient oscillated due to the vortex shedding is shown below. To compare the effect of the timestep size, we continue to plot the lift coefficient after the oscillations have been established with a timestep of .02 seconds.
Timestep = .2 s
Timestep = .2 -> .02s
Strouhal Number Comparison