A divergence‐free interpolation scheme for the immersed boundary method

Abstract

SUMMARYThe immersed boundary approach for the modeling of complex geometries in incompressible flows is examined critically from the perspective of satisfying boundary conditions and mass conservation. It is shown that the system of discretized equations for mass and momentum can be inconsistent, if the velocity is used in defining the force density to satisfy the boundary conditions. As a result, the velocity is generally not divergence free and the pressure at locations in the vicinity of the immersed boundary is not physical. However, the use of the pseudo-velocities in defining the force density, as frequently done when the governing equations are solved using a fractional step or projection method, combined with the use of the specified velocity on the immersed boundary, is shown to result in a consistent set of equations which allows a divergence-free velocity but, depending on the time step, is shown to have the undesirable effects of inaccurately satisfying the boundary conditions and allowing a significant permeability of the immersed boundary. If the time step is reduced sufficiently, the boundary conditions on the immersed boundary can be satisfied. However, this entails an unacceptable increase in computational expense. Two new methods that satisfy the boundary conditions and allow a divergence-free velocity while avoiding the increased computational expense are presented and shown to be second-order accurate in space. The first new method is based on local time step reduction. This method is suitable for problems where the immersed boundary does not move. For these problems, the first new method is shown to be closely related to the second new method. The second new method uses an optimization scheme to minimize the deviation from the interpolation stencil used to represent the immersed boundary while ensuring a divergence-free velocity. This method performs well for all problems, including those where the immersed boundary moves relative to the grid. Additional results include showing that the force density that is added to satisfy the boundary conditions at the immersed boundary is unbounded as the time step is reduced and that the pressure in the vicinity of the immersed boundary is unphysical, being strongly a function of the time step. A method of computing the total force on an immersed boundary which takes into account

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

“…Thus, strict conservation of quantities such as mass, momentum or kinetic energy is not observed near the irregular boundary. The most severe manifestations of these shortcomings is the occurrence of non-divergence free velocities or unphysical oscillations of the pressure in the vicinity of the immersed boundary [43,30]. Numerous revisions of these interpolations are still proposed for improving the accuracy and consistency of this class of IB methods [30,3,49,43].…”

A second-order cut-cell method for the numerical simulation of 2D flows past obstacles

“…The coupling of momentum forcing methods with projection schemes introduces difficulties in imposing at the same time level the contuinity equation and the boundary conditions on the immersed interface (see [9] for example). The incompressibility of the flow may even be violated in the vicinity of the immersed boundary [10]. This may result in serious diffulties when computing boundary layers at high Reynolds number.…”

Accurate Cartesian-grid simulations of near-body flows at intermediate Reynolds numbers

“…A well known issue with this class of algorithms is their tendency to introduce large non-physical pressure oscillations (see [8] for example). Muldoon [9] even showed that the pressure could locally increase without bound as the time step goes to zero. These oscillations are caused by the lack of smoothness of the velocity across the boundary before the projection step [10].…”