An immersed boundary velocity correction method combined with virtual body-fitted grid for simulation of incompressible flows

In this work, a virtual body-fitted grid is introduced into the velocity correction-based immersed boundary method (IBM) to simulate incompressible flows. The impact of the immersed boundary is indirectly transmitted to the flow field via a virtual body-fitted grid. In this method, the fractional step technique consisting of the predictor and the corrector is adopted. The prediction step is executed on the Eulerian mesh, and the correction step is done on the virtual grid to fulfill the no-slip boundary condition. After the correction step, the corrected velocity field on the virtual grid is then assigned to that on the Eulerian mesh to update the flow field. Being able to adjust the grid spacing flexibly, the virtual body-fitted grid alleviates the shortcomings of the conventional IBM that uses the smooth Dirac delta function to associate Lagrangian points with their surrounding Eulerian points. As a result, the present method is easy to apply to non-uniform Cartesian grids, which is inapplicable to the conventional IBM with the smooth Dirac delta function. Numerical experiments concerning flow past a circular cylinder and a NACA0012 airfoil demonstrate the advantages of the present method, i.e., fewer Lagrangian points are required to avoid the streamline penetration of boundary and the range of “diffuse interface” can be narrowed by reducing the normal grid spacing of the virtual body-fitted grid to improve numerical results on a coarse mesh. In addition, an accuracy assessment on the decaying vortex problem reveals that the present IBM has a second-order accuracy.