High-order fluid–structure interaction in 2D and 3D application to blood flow in arteries

Abstract

a b s t r a c tThis paper addresses the numerical approximation of fluid-structure interaction (FSI) problems through the arbitrary Lagrangian Eulerian (ALE) framework, high-order methods and a Dirichlet-Newmann approach for the coupling. The paper is divided into two main parts. The first part concerns the discretization method for the FSI problem. We introduce an improved ALE map, capable of handling curved geometries in 2D and 3D in a unified manner, that is based on a local differential operator. We also propose a minimal continuous interior penalty (CIP) stabilization term for the fluid discretization that accounts for a smaller computational effort, while stabilizing the flow regime. The second part is dedicated to validating our numerical strategy through a benchmark and some applications to blood flow in arteries.

Simulation of two-fluid flows using a finite element/level set method. Application to bubbles and vesicle dynamics

“…The previous section described the strategy we used to track the interface. We couple it now to the Navier Stokes equation solver described in [17]. In the current section, we present a validation of this two-fluid flow framework.…”

A hybridizable discontinuous Galerkin method for modeling fluid–structure interaction

“…The majority of FSI research has dealt with low-order (second or less) modeling; however, at least two groups are researching high-order Galerkin Finite Element Method based FSI, including Persson et al [37,38] and Pena et al [39,40]. Perrson et al use a mixed spatial discretization, with DG for the fluid and CG for the solid, while Pena et al have developed a framework for high-order FSI capable of either CG or DG formulations.…”

Hemodynamic simulations in the cerebral venous network: A study on the influence of different modeling assumptions

“…Before investigating the impact of several modeling assumptions, our objective in the present section is to assess the impact of the numerical strategy in use. We ensure that our sensitivity analysis is performed in a controlled numerical environment within the open-source developed and validated software library Feel++, see section 3.2, and in [5,6,9,8,4]. Our numerical strategy is in line with [24] as we execute high resolution simulations and not normal resolution simulations: both the spacial and the temporal discretization methods are second-order schemes and a small time step (∆t = 10 −3 s) was chosen, in order to adequately resolve the complex flow features.…”