提出了一种模拟高雷诺数湍流的扩散界面浸入边界法（IBM）。该方法采用壁面模型来减少壁面附近的网格量。为了实施壁面模型和边界条件，在物面外部设置了2个辅助层（一系列拉格朗日点）：外侧的参考层用于实施壁面模型来得到壁面剪应力，内侧的强制层用于实施离壁的边界条件。在实施壁面模型时，将动量方程沿物面法向积分，从而把物面切向的速度修正量与由壁面模型得到的壁面剪应力联系起来，而速度的法向分量则按二次曲线分布来近似重构。在实施边界条件时，为了严格满足无穿透条件，应用了基于隐式速度修正的IBM。最后，利用绕平板湍流和绕NACA0012翼型湍流来对方法的可行性进行了验证。
A diffuse-interface Immersed Boundary Method (IBM) for simulating high Reynolds number turbulent flows is proposed. This method uses the wall model to reduce the number of grids near the wall. To implement the wall model and boundary conditions, two auxiliary layers (a series of Lagrangian points) are set outside the wall: the outer reference layer carries out the wall model to calculate the wall shear stress, and the inner enforced layer implements boundary conditions off the wall. When the wall model is implemented, the momentum equation is integrated along the normal direction of the wall, such that the tangential velocity correction is related to the wall shear stress obtained by the wall model, and the normal component of the velocity is approximately reconstructed according to the parabolic distribution. Furthermore, when implementing boundary conditions, we adopt the implicit velocity correction-based IBM to satisfy the non-penetration condition. Finally, the feasibility of this method is verified by numerical simulation of the turbulent flows around the flat plate and the NACA0012 airfoil.
[1] PESKIN C S. Numerical analysis of blood flow in the heart[J]. Journal of Computational Physics, 1977, 25(3): 220-252.
[2] CABOT W, MOIN P. Approximate wall boundary conditions in the large-eddy simulation of high Reynolds number flow[J]. Flow, Turbulence and Combustion, 1999, 63(1): 269-291.
[3] SHU C, LIU N, CHEW Y T. A novel immersed boundary velocity correction-lattice Boltzmann method and its application to simulate flow past a circular cylinder[J]. Journal of Computational Physics, 2007, 226(2): 1607-1622.
[4] WU J, SHU C. Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications[J]. Journal of Computational Physics, 2009, 228(6): 1963-1979.
[5] SPALDING D B. A single formula for the "law of the wall"[J]. Journal of Applied Mechanics, 1961, 28(3): 455-458.
[6] WANG M, MOIN P. Dynamic wall modeling for large-eddy simulation of complex turbulent flows[J]. Physics of Fluids, 2002, 14(7): 2043-2051.
[7] TESSICINI F, IACCARINO G, FATICA M, et al. Wall modeling for large-eddy simulation using an immersed boundary method[R].Center for Turbulence Research, 2002.
[8] KALITZIN G, IACCARINO G. Toward immersed boundary simulation of high reynolds number flows[R].Center for Turbulence Research, 2003.
[9] CRISTALLO A, VERZICCO R. Combined immersed boundary/large-eddy-simulations of incompressible three dimensional complex flows[J]. Flow, Turbulence and Combustion, 2006, 77(1-4): 3-26.
[10] LEE J D, RUFFIN S. Development of a turbulent wall-function based viscous Cartesian-grid methodology[C]//45th AIAA Aerospace Sciences Meeting and Exhibit. Reston: AIAA, 2007.
[11] CAPIZZANO F. Turbulent wall model for immersed boundary methods[J]. AIAA Journal, 2011, 49(11): 2367-2381.
[12] CAPIZZANO F. Coupling a wall diffusion model with an immersed boundary technique[J]. AIAA Journal, 2016, 54(2): 728-734.
[13] JI C, MUNJIZA A, WILLIAMS J J R. A novel iterative direct-forcing immersed boundary method and its finite volume applications[J]. Journal of Computational Physics, 2012, 231(4): 1797-1821.
[14] CHANG P H, LIAO C C, HSU H W, et al. Simulations of laminar and turbulent flows over periodic hills with immersed boundary method[J]. Computers & Fluids, 2014, 92: 233-243.
[15] ZHOU C H. RANS simulation of high-Re turbulent flows using an immersed boundary method in conjunction with wall modeling[J]. Computers & Fluids, 2017, 143: 73-89.
[16] PU T M, ZHOU C H. An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows[J]. International Journal for Numerical Methods in Fluids, 2018, 87(5): 217-238.
[17] MITTAL R, IACCARINO G. Immersed boundary methods[J]. Annual Review of Fluid Mechanics, 2005, 37(1): 239-261.
[18] SHI B J, YANG X L, JIN G D, et al. Wall-modeling for large-eddy simulation of flows around an axisymmetric body using the diffuse-interface immersed boundary method[J]. Applied Mathematics and Mechanics, 2019, 40(3): 305-320.
[19] WANG Y, SHU C, YANG L M. Boundary condition-enforced immersed boundary-lattice Boltzmann flux solver for thermal flows with Neumann boundary conditions[J]. Journal of Computational Physics, 2016, 306: 237-252.
[20] KALITZIN G, MEDIC G, IACCARINO G, et al. Near-wall behavior of RANS turbulence models and implications for wall functions[J]. Journal of Computational Physics, 2005, 204(1): 265-291.
[21] MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8): 1598-1605.
[22] KNOPP T, ALRUTZ T, SCHWAMBORN D. A grid and flow adaptive wall-function method for RANS turbulence modelling[J]. Journal of Computational Physics, 2006, 220(1): 19-40.
[23] SHU C, WANG Y, TEO C J, et al. Development of lattice boltzmann flux solver for simulation of incompressible flows[J]. Advances in Applied Mathematics and Mechanics, 2014, 6(4): 436-460.
[24] JESPERSEN D C, PULLIAM T H, CHILDS M L. OVERFLOW turbulence modeling resource validation results[R].Washington, D.C.: NASA Ames Research Center, 2016.
CJA