Improved immersed boundary/wall modeling method for RANS solver coupled with wall functions: application to Cartesian grid systems

A two-dimensional (2D) Reynolds-averaged Navier–Stokes (RANS) equations solver with k–ω turbulence closure is developed, employing immersed boundary (IB) technique on Cartesian grids. Generalized wall functions are introduced to enhance computational efficiency for problems with high Reynolds numbers. To address existing challenges in applying wall functions within IB methods, a novel, effective and easy-to-implement strategy is proposed. Another distinguishing feature of this turbulent-flow solver is that it employs the highly accurate immersed-boundary generalized harmonic polynomial cell (IB-GHPC) method to solve the Poisson equation for fluid pressure. The new solver is firstly validated by simulating channel flows on both hydraulically smooth and rough walls, achieving excellent agreement with benchmark experimental and numerical studies for various flow parameters including velocity, turbulent kinetic energy and shear stress. For channel flow simulations, our implementation of generalized wall functions using the proposed strategy results in a remarkable reduction of grid nodes by over 80%. Moreover, the solver is applied to simulate flow around both smooth and rough cylinders, producing promising results for drag, lift, and pressure coefficients. Our analysis demonstrates a robust performance of the developed solver in modeling turbulent flows based on Cartesian grids, offering a substantial improvement in computational efficiency for tackling problems involving large Reynolds numbers.