Abstract:We present a new model and a novel loosely coupled partitioned numerical scheme modeling fluid-structure interaction (FSI) in blood flow allowing non-zero longitudinal displacement. Arterial walls are modeled by a {linearly viscoelastic, cylindrical Koiter shell model capturing both radial and longitudinal displacement}. Fluid flow is modeled by the Navier-Stokes equations for an incompressible, viscous fluid. The two are fully coupled via kinematic and dynamic coupling conditions. Our numerical scheme is based on a new modified Lie operator splitting that decouples the fluid and structure sub-problems in a way that leads to a loosely coupled scheme which is {unconditionally} stable. This was achieved by a clever use of the kinematic coupling condition at the fluid and structure sub-problems, leading to an implicit coupling between the fluid and structure velocities. The proposed scheme is a modification of the recently introduced "kinematically coupled scheme" for which the newly proposed modified Lie splitting significantly increases the accuracy. The performance and accuracy of the scheme were studied on a couple of instructive examples including a comparison with a monolithic scheme. It was shown that the accuracy of our scheme was comparable to that of the monolithic scheme, while our scheme retains all the main advantages of partitioned schemes, such as modularity, simple implementation, and low computational costs.

What problem does this paper attempt to address?

This paper aims to solve the fluid - structure interaction (FSI) problem in blood flow, especially the modeling and numerical solution when non - zero longitudinal displacement occurs in the arterial wall. Specifically, the paper attempts to: 1. **Improve the existing fluid - structure interaction models**: Traditional fluid - structure interaction models usually only consider radial displacement and ignore the longitudinal displacement of the arterial wall under the action of blood flow. In this paper, by introducing a linear viscoelastic cylindrical shell model (Koiter shell model), the radial and longitudinal displacements of the arterial wall are captured simultaneously, thus describing the movement of the arterial wall more accurately. 2. **Develop a new numerical solution scheme**: In order to solve the complex fluid - structure interaction problem efficiently and stably, this paper proposes a new loosely - coupled partitioned numerical scheme. This scheme is based on the improved Lie operator splitting method. By skillfully using the dynamic and kinematic coupling conditions, the implicit coupling of the fluid and structure sub - problems is achieved, ensuring the unconditional stability of the numerical scheme and improving the computational accuracy. 3. **Verify the effectiveness of the new model and numerical scheme**: By comparing with the results of the monolithic scheme, it is proved that the proposed loosely - coupled scheme not only has comparable accuracy, but also retains the advantages of the partitioned solver, such as modularity, easy implementation and low computational cost. ### Main contributions - **Model innovation**: For the first time, a linear viscoelastic cylindrical shell model that can capture both radial and longitudinal displacements is introduced, which is more in line with the actual physiological situation. - **Improvement of numerical methods**: A new loosely - coupled numerical scheme based on the improved Lie splitting is proposed, which solves the stability problem in the traditional loosely - coupled method and significantly improves the computational accuracy. - **Application verification**: The effectiveness and accuracy of the new method are verified through multiple examples, especially in simulating the movement of the arterial wall. ### Key formulas - **Energy expression of the cylindrical shell model**: \[ E(\eta)=\frac{h}{2}\int_{0}^{L}(A\gamma(\eta)+B\gamma(\dot{\eta})):\gamma(\xi)R \, dz+\frac{h^{3}}{24}\int_{0}^{L}(A\rho(\eta)+B\rho(\dot{\eta})):\rho(\xi)R \, dz \] where $\gamma(\eta)$ and $\rho(\eta)$ represent the changes of the strain tensor and the curvature tensor respectively. - **Fluid dynamic equations**: \[ \rho_{f}\left(\frac{\partial u}{\partial t}+u\cdot\nabla u\right)=\nabla\cdot\sigma \quad \text{in } \Omega(t) \] \[ \nabla\cdot u = 0 \quad \text{in } \Omega(t) \] where $\sigma=-pI + 2\mu D(u)$ is the fluid stress tensor, and $D(u)=\frac{1}{2}(\nabla u+(\nabla u)^{T})$ is the strain - rate tensor. - **Coupling conditions**: - Kinematic coupling condition (velocity continuity): \[ u(\hat{z}+\eta_{z}(\hat{z},t),R+\eta_{r}(\hat{z},t),t)=\frac{\partial\eta}{\partial t}(\hat{z},t) \] - Dynamic coupling condition (contact force balance): \[ f_{z}=-J\hat{\sigma}n|_{\Gamma(t)}\cdot e_{z} \] \[ f_{r}=-J

Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement

Existence of a weak solution to a nonlinear fluid-structure interaction problem modeling the flow of an incompressible, viscous fluid in a cylinder with deformable walls

Efficient split-step schemes for fluid-structure interaction involving incompressible generalised Newtonian flows

Existence and Regularity Results for a Nonlinear Fluid-Structure Interaction Problem with Three-Dimensional Structural Displacement

On the Impact of Fluid Structure Interaction in Blood Flow Simulations: Stenotic Coronary Artery Benchmark

A monolithic fluid–structure interaction approach using mixed LSFEM with high-order time integration

Vascular fluid-structure interaction: unified continuum formulation, image-based mesh generation pipeline, and scalable fully implicit solver technology

Stability and convergence analysis of the kinematically coupled scheme for fluid-structure interaction

Numerical modeling of fluid-structure interaction in arteries with anisotropic polyconvex hyperelastic and anisotropic viscoelastic material models at finite strains

A reduced unified continuum formulation for vascular fluid-structure interaction

A Fully Parallelizable Loosely Coupled Scheme for Fluid-Poroelastic Structure Interaction Problems

Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm

An immersed boundary method for the fluid-structure interaction of slender flexible structures in viscous fluid

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

A sharp MLS penalty immersed finite element method for fluid-structure interaction of highly deformable slender body in turbulent flow

Modeling blood flow in viscoelastic vessels: the 1D augmented fluid-structure interaction system

A stable and explicit fluid–structure interaction solver based on lattice-Boltzmann and immersed boundary methods

Monolithic Approach to Thermal Fluid-Structure Interaction with Nonconforming Interfaces

An immersed boundary-lattice Boltzmann method with hybrid multiple relaxation times for viscoplastic fluid-structure interaction problems

Fluid Structure Interaction Using Modal Superposition and Lagrangian CFD

A monolithic fluid-structure interaction formulation for solid and liquid membranes including free-surface contact