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

TL;DR

This paper introduces a simple, efficient, and second-order accurate immersed-boundary method for solving the incompressible Navier--Stokes equations, capable of handling complex geometries with high accuracy and stability.

Contribution

It presents an implicit, second-order immersed-boundary method that simplifies implementation and improves stability without requiring matrix inversions or complex corrections.

Findings

Achieves second-order accuracy in complex geometries

Demonstrates stability and efficiency in turbulent and anatomical flows

No corrections needed for pressure or continuity equations

Abstract

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 implemented alongside 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.