A numerical method for solving the 3D unsteady incompressible Navier–Stokes equations in curvilinear domains with complex immersed boundaries

Abstract

A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [1]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow through a mechanical, bileaflet heart valve mounted in a model straight aorta with an anatomical-like triple sinus.

Numerical Study on Hydrodynamic Performance of Bionic Caudal Fin

“…However, a large number of problems in biofluid dynamics involve interactions between deformable elastic bodies and incompressible viscous fluids. These fluid-structure interaction (FSI) problems involve the complicated interplay between a viscous fluid, deformable body, and free-moving boundary, making them difficult to discern [12][13][14][15][16][17][18][19][20][21][22][23].…”

Numerical Study on the Effect of Air–Sea–Land Interaction on the Atmospheric Boundary Layer in Coastal Area

“…To address this issue, we use an LES tool called CURVIB-LES [13,14]. The CURVIB-LES is based on the IB method [15].…”

Emerging Trends in Heart Valve Engineering: Part IV. Computational Modeling and Experimental Studies

“…Accordingly, well-resolved three-dimensional simulations require high spatial resolution that presently necessitates the use of high-performance computing resources. The IB method was introduced to model heart valves, 124,125 and both fictitious domain [37][38][39][40]147 and the IB methods 50,54,[59][60][61][62]100,101,110,111,121,127 have been used substantially in modeling heart valve dynamics.…”