Lid driven cavity

edit 

Related papers

G. Kosec; A local numerical solution of a fluid-flow problem on an irregular domain, Advances in engineering software, vol. 120, 2018 [DOI: 10.1016/j.advengsoft.2016.05.010]

Click here to return back to Fluid Mechanics

Intro

Let's try MLSM and described solution procedures first on a lid driven cavity problem that stands for a standard benchmark test for validation of the fluid flow solvers. It has been proposed in 1982 and since then solved by many researchers with wide spectra of different numerical methods. The test is still widely studied and used for validation of novel methods and numerical principles. Case is schematically presented in Figure 1 followed by few basic analyses, just to show that things work as they should. Details about comparison can be found in paper

Figure 1: Scheme of a lid driven cavity case

Code

The snippet of the MLSM code for an explicit ACM method with CBS looks like: (full examples, including implicit versions, can be found under the examples in the code repository).

 1     for (int step = 0; step <=O.t_steps; ++step) {
 2         // Navier-Stokes
 3         for (auto c:interior) 
 4             // Navier-Stokes
 5             v2[c] = v1[c] + O.dt * ( - 1/O.rho    * op.grad(P1,c)
 6                                      + O.mu/O.rho * op.lap(v1, c)
 7                                      -              op.grad(v1,c)*v1[c]);
 8         }
 9         // Mass continuity
10         for (auto i:interior) { 
11             P2[i] = P1[i] - O.dl * O.dt * O.rho * op.div(v2,i) + O.dl2 * O.dl * O.dt * O.dt * op.lap(P1,i);
12         }
13         v1.swap(v2);
14         P1.swap(P2);
15     }

Results

Figure 2: Cross section velocity profiles.