A simple and efficient second-order immersed-boundary method for the incompressible Navier–Stokes equations

An immersed-boundary method for the incompressible Navier–Stokes equations is presented. It employs discrete forcing for a sharp discrimination of the solid-fluid interface, and achieves second-order accuracy, demonstrated in examples with highly complex three-dimensional geometries. The method is implicit, meaning that the point in the solid which is nearest to the interface is accounted for implicitly, which benefits stability and convergence properties; the correction is also implicit in time (without requiring a matrix inversion), although the temporal integration scheme is fully explicit. The method stands out for its simplicity and efficiency: when integrated with second-order finite differences, only the weight of the center point of the Laplacian stencil in the momentum equation is modified, and no corrections for the continuity equation and the pressure are required. The immersed-boundary method, its performance and its accuracy are first verified on simple problems, and then put to test on a simple laminar, two-dimensional flow and on two more complex examples: the turbulent flow in a channel with a sinusoidal wall, and the flow in a human nasal cavity, whose extreme anatomical complexity mandates an accurate treatment of the boundary.