Computational Fluid Dynamics as good as it gets.

M. Meyer, A. Devesa, S. Hickel, X.Y. Hu, N.A. Adams (2010)
Journal of Computational Physics 229: 6300-6317. doi: 10.1016/j.jcp.2010.04.040

We propose a conservative, second-order accurate immersed interface method for representing incompressible fluid flows over complex three dimensional solid obstacles on a staggered Cartesian grid. The method is based on a finite-volume discretization of the incompressible Navier–Stokes equations which is modified locally in cells that are cut by the interface in such a way that accuracy and conservativity are maintained.

A level-set technique is used for description and tracking of the interface geometry, so that an extension of the method to moving boundaries and flexible walls is straightforward. Numerical stability is ensured for small cells by a conservative mixing procedure. Discrete conservation and sharp representation of the fluid–solid interface render the method particularly suitable for Large-Eddy Simulations of high-Reynolds number flows. Accuracy, second-order grid convergence and robustness of the method is demonstrated for several test cases: inclined channel flow at Re = 20, flow over a square cylinder at Re = 100, flow over a circular cylinder at Re = 40, Re = 100 and Re = 3900, as well as turbulent channel flow with periodic constrictions at Re = 10595.

 

Control volume for the velocity components u – – and v  · – ·, and control volume for the pressure p for a staggered grid in 2D; arrows indicate locations of velocity components, circles locations of pressure.

 

Fluid volume fraction αi,j,, intersection of (Vi,j,k ∩ Ω) and interface with cell Γi,j,k, for the cut cell (i, j, k).

 

Linear approximation of the velocity gradient of a velocity cell for the calculation of the viscous interface interaction term.

 

Mixing procedure for the small fluid cell (i, j, k).

 

AMR mesh for the LES of the flow over a circular at Re = 3900.
Results for LES of the flow over a circular at Re = 3900: (a) streamwise vorticity-isosurface ±1.5 s−1 ; (b) pressure coefficient Cp along the surface of a circular cylinder: θ = 0 in the stagnation point; □ experimental data Norberg extracted from Kravchenko and Moin (2000); – – simulation of Kravchenko and Moin (2000);  — present method.