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

  • Authors:
  • Martina Bukač;SunčIca Anić;Roland Glowinski;Josip TambačA;Annalisa Quaini

  • Affiliations:
  • Department of Mathematics, University of Houston, 4800 Calhoun Rd, Houston, TX 77204, USA;Department of Mathematics, University of Houston, 4800 Calhoun Rd, Houston, TX 77204, USA;Department of Mathematics, University of Houston, 4800 Calhoun Rd, Houston, TX 77204, USA and Laboratoire Jacques-Louis Lions, Université P. et M. Curie, 4 Place Jussieu, 75005 Paris, France;Department of Mathematics, University of Zagreb, Bijenička 30, 10000 Zagreb, Croatia;Department of Mathematics, University of Houston, 4800 Calhoun Rd, Houston, TX 77204, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2013

Quantified Score

Hi-index 31.45

Visualization

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.