An immersed interface method for discrete surfaces

JOURNAL OF COMPUTATIONAL PHYSICS

卷:400

Article

Kolahdouz, Ebrahim M.^1,2  Bhalla, Amneet Pal Singh³  Craven, Brent A.²  Griffith, Boyce E.^1,4,5,6,7,8 

[1] Univ N Carolina, Dept Math, Chapel Hill, NC 27515 USA

[2] US FDA, Div Appl Mech, Off Sci & Engn Labs, Ctr Devices & Radiol Hlth, Silver Spring, MD USA

[3] San Diego State Univ, Dept Mech Engn, San Diego, CA 92182 USA

[4] Univ N Carolina, Dept Appl Phys Sci, Chapel Hill, NC 27515 USA

[5] Univ N Carolina, Dept Biomed Engn, Chapel Hill, NC 27515 USA

[6] Univ N Carolina, Carolina Ctr Interdisciplinary Appl Math, Chapel Hill, NC 27515 USA

[7] Univ N Carolina, Computat Med Program, Sch Med, Chapel Hill, NC 27515 USA

[8] Univ N Carolina, McAllister Heart Inst, Sch Med, Chapel Hill, NC 27515 USA

关键词: Immersed boundary method;    Immersed interface method;    Finite element;    Fluid-structure interaction;    Jump conditions;    Complex geometries;   

DOI  :  10.1016/j.jcp.2019.07.052

来源: Elsevier

【 摘 要 】

Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary (IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction (FSI) in such systems, but a difficulty of the IB formulation of these problems is that the pressure and viscous stress are generally discontinuous at fluid-structure interfaces. The conventional IB method regularizes these discontinuities, which typically yields low-order accuracy at these interfaces. The immersed interface method (IIM) is an IB-like approach to FSI that sharply imposes stress jump conditions, enabling higher-order accuracy, but prior applications of the IIM have been largely restricted to numerical methods that rely on smooth representations of the interface geometry. This paper introduces an immersed interface formulation that uses only a C-0 representation of the immersed interface, such as those provided by standard nodal Lagrangian finite element methods. Verification examples for models with prescribed interface motion demonstrate that the method sharply resolves stress discontinuities along immersed boundaries while avoiding the need for analytic information about the interface geometry. Our results also demonstrate that only the lowest-order jump conditions for the pressure and velocity gradient are required to realize global second-order accuracy. Specifically, we demonstrate second-order global convergence rates along with nearly second-order local convergence in the Eulerian velocity field, and between first- and second-order global convergence rates along with approximately first-order local convergence for the Eulerian pressure field. We also demonstrate approximately second-order local convergence in the interfacial displacement and velocity along with first-order local convergence in the fluid traction along the interface. As a demonstration of the method's ability to tackle more complex geometries, the present approach is also used to simulate flow in a patient-averaged anatomical model of the inferior vena cava, which is the large vein that carries deoxygenated blood from the lower and middle body back to the heart. Comparisons of the general hemodynamics and wall shear stress obtained by the present IIM and a body-fitted discretization approach show that the present method yields results that are in good agreement with those obtained by the body-fitted approach. (C) 2019 Elsevier Inc. All rights reserved.

【 授权许可】

Free