Modification of the LS-STAG Immersed Boundary Method for Simulating Turbulent Flows

14 / 16

I.K. Marchevsky, V.V. Puzikova
32 ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2017. № 5
Flow simulation at high Reynolds numbers 5 7 (Re =10 10 ). The flow was
simulated at the Reynolds numbers 5 7 Re =10 10 (on non-uniform meshes
240 296 with 4 = 5 10 t    and 480 592 with 4 =10 t   ). Results obtained on mesh
480 592 are very close to experimental data [18], see Fig. 6. At 5 Re = 2 10 , the
boundary layer on the cylinder surface undergoes a transition from laminar to
turbulent [19]. This transition leads to a delay of the separation of flow from the
cylinder surface causing a substantial reduction in the drag force. This is often
referred to as ''drag crisis''. This phenomenon was simulated by using modified
LS-STAG immersed boundary method (Fig. 6).
Conclusions. The key points of the LS-STAG method [2, 21] extension for
RANS/LES/DES turbulence models were described. For the shear Reynolds stresses
and for the eddy viscosity an additional mesh (xy-mesh) is introduced. The general
approach to the construction of the LS-STAG discretization for differential equations
of the EVM models on the additional xy-mesh shown. The Smagorinsky, Spalart —
Allmaras, , k   k  and k  SST turbulence models are considered. To validate
modified LS-STAG immersed boundary method the flow past a circular airfoil at
2 7 Re =10 10 was simulated. Computational results are in good agreement with
established results from the literature. Also, the so-called ''drag crisis'' phenomenon of
circular cylinder at 5 6 Re =10 10 was simulated.
REFERENCES
[1] Mittal R., Iaccarino G. Immersed boundary methods. Ann. Rev. Fluid Mech., 2005, vol. 37,
pp. 239–261. DOI: 10.1146/annurev.fluid.37.061903.175743
Available at:
[2] Cheny Y., Botella O. The LS-STAG method: A new immersed boundary/level-set method for
the computation of incompressible viscous flows in complex moving geometries with good
conservation properties. J. Comput. Phys., 2010, vol. 229, iss. 4, pp. 1043–1076.
DOI: 10.1016/j.jcp.2009.10.007
[3] Spalart P.R. Strategies for turbulence modelling and simulations. Int. J. Heat and Fluid Flow.
2000, vol. 21, iss. 3, pp. 252–263. DOI: 10.1016/S0142-727X(00)00007-2
[4] Spalart P.R., Allmaras S.R. A one-equation turbulence model for aerodynamic flows. Recherche
Aerospatiale, 1994, no. 1, pp. 5–21.
[5] Jones W.P., Launder B.E. The prediction of laminarization with a two-equation model of
turbulence. Int. J. Heat Mass Transfer, 1972, vol. 15, iss. 2, pp. 301–314.
DOI: 10.1016/0017-9310(72)90076-2
[6] Wilcox D.C. Reassessment of the scale-determining equation for advanced turbulence models.
AIAA Journal, 1988, vol. 26, no. 11, pp. 1299–1310. DOI: 10.2514/3.10041
Available at:
[7] Menter F.R. Two-equation eddy-viscosity turbulence models for engineering applications.
AIAA Journal, 1994, vol. 32, no. 8, pp. 1598–1605. DOI: 10.2514/3.12149
Available at: