Lid-driven cavity

 

he lid-driven cavity case, described also in the tutorial cases, may be used as a validation case of the flow solver inside AeroSolved. The well-known setting of lid-driven cavity flow has been studied by many authors. In Ghia et al. (1982) data are presented of the velocity profiles across the vertical and horizontal center lines of a two-dimensional cavity. The lid-driven cavity case, available in the AeroSolved code repository under the path: cases/validation/lidDrivenCavity, contains a configuration to simulate a cavity at a Reynolds number of 10000. The Reynolds number is controlled by the viscosity, which is set to 10-4 kg/(m s). All other quantities are set to unity, i.e., density ρ = 1 kg/m3, lid velocity U = 1 m/s and length L = 1 m. The adopted mesh is 300 x 300 x 1 cells. The simulation runs until 200 s on 9 processor cores in parallel, and is then sampled along the horizontal and vertical center lines, to yield the two relevant velocity profiles. The validation case can be executed by calling the following scripts for preparation, running and post processing: prep.sh, run.sh (will take a long time) and post.sh. The post.sh script requires Python with Matplotlib to be installed, and creates two .pdf files (plotU.pdf and plotV.pdf) that contain the comparison of the simulated velocity profiles and the reference velocity profiles (Ghia et al.). These two plots are shown below, where u,v are dimension less velocities scaled by U and x,y are the dimension less spatial coordinates scaled by L.

plotUplotV