skip to main content


Title: Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm
Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems but need improvements when the problems become harder when (a) the geometry becomes complex (eg, a few branches to a full pulmonary artery), (b) the model becomes more complex (eg, fluid-only to coupled fluid-structure interaction), (c) both the fluid and wall models become highly nonlinear, and (d) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from computed tomography (CT) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations, and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know, this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10 000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.  more » « less
Award ID(s):
1720366
NSF-PAR ID:
10155913
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
International journal for numerical methods in biomedical engineering
ISSN:
2040-7939
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient‐specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems but need improvements when the problems become harder when (a) the geometry becomes complex (eg, a few branches to a full pulmonary artery), (b) the model becomes more complex (eg, fluid‐only to coupled fluid‐structure interaction), (c) both the fluid and wall models become highly nonlinear, and (d) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid‐structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient‐specific, full size pulmonary artery obtained from computed tomography (CT) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier‐Stokes equations, and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know, this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10 000 processor cores on a supercomputer for solving the fluid‐structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.

     
    more » « less
  2. Abstract. We consider a nonlinear, moving boundary, fluid-structure interaction problem between a time dependent incompressible, viscous fluid flow, and an elastic structure composed of a cylindrical shell supported by a mesh of elastic rods. The fluid flow is modeled by the time-dependent Navier- Stokes equations in a three-dimensional cylindrical domain, while the lateral wall of the cylinder is modeled by the two-dimensional linearly elastic Koiter shell equations coupled to a one-dimensional system of conservation laws defined on a graph domain, describing a mesh of curved rods. The mesh supported shell allows displacements in all three spatial directions. Two-way coupling based on kinematic and dynamic coupling conditions is assumed between the fluid and composite structure, and between the mesh of curved rods and Koiter shell. Problems of this type arise in many ap- plications, including blood flow through arteries treated with vascular prostheses called stents. We prove the existence of a weak solution to this nonlinear, moving boundary problem by using the time discretization via Lie operator splitting method combined with an Arbitrary Lagrangian-Eulerian approach, and a non-trivial extension of the Aubin-Lions-Simon compactness result to problems on moving domains. 
    more » « less
  3. Simulation of blood flows in the pulmonary artery provides some insight into certain diseases by examining the relationship between some continuum metrics, e.g., the wall shear stress acting on the vascular endothelium, which responds to flow-induced mechanical forces by releasing vasodilators/constrictors. V. Kheyfets, in his previous work, studies numerically a patient-specific pulmonary circulation to show that decreasing wall shear stress is correlated with increasing pulmonary vascular impedance. In this paper, we develop a scalable parallel algorithm based on domain decomposition methods to investigate an unsteady model with patient-specific pulsatile waveforms as the inlet boundary condition. 
    more » « less
  4. Summary

    We present a spatially varying Robin interface condition for solving fluid‐structure interaction problems involving incompressible fluid flows and nonuniform flexible structures. Recent studies have shown that for uniform structures with constant material and geometric properties, a constant one‐parameter Robin interface condition can improve the stability and accuracy of partitioned numerical solution procedures. In this work, we generalize the parameter to a spatially varying function that depends on the structure's local material and geometric properties, without varying the exact solution of the coupled fluid‐structure system. We present an algorithm to implement the Robin interface condition in an embedded boundary method for coupling a projection‐based incompressible viscous flow solver with a nonlinear finite element structural solver. We demonstrate the numerical effects of the spatially varying Robin interface condition using two example problems: a simplified model problem featuring a nonuniform Euler‐Bernoulli beam interacting with an inviscid flow and a generalized Turek‐Hron problem featuring a nonuniform, highly flexible beam interacting with a viscous laminar flow. Both cases show that a spatially varying Robin interface condition can clearly improve numerical accuracy (by up to two orders of magnitude in one instance) for the same computational cost. Using the second example problem, we also demonstrate and compare two models for determining the local value of the combination function in the Robin interface condition.

     
    more » « less
  5. Abstract

    Spiral/helical forms of blood flow have been observed in large arteries of the cardiovascular system, but their benefits remain underappreciated. Spiral flow has been postulated to improve near‐wall washout, promoting anti‐atherothrombotic conditions. This research aims to study the washout characteristics of spiral flow, specifically, its ability to increase velocity and wall shear stress (WSS) in atherothrombotic‐prone regions. Using 1.2 cm diameter angled test‐conduits (45°, 90°, 135°) with known recirculation/stasis regions at the bend corners, spiral flow washout potential was evaluated in terms of low velocity and low WSS. Two sub‐studies were conducted: the first utilized a spiral flow‐inducing device to enable qualitative analysis of washout‐potential in both computational fluid dynamic (CFD) simulations and benchtop ultrasound visualization; the second used CFD to study the impact of several induced helical wavelengths on the conduit‐dependent recirculation/stasis zones. Physical models of the angled conduits and spiral flow‐inducer were 3D‐printed to facilitate ultrasound visualization. Compared to straight flow, spiral flow generated by the flow‐inducer significantly cleared the recirculation/stasis zones at the corners of the angled conduits. CFD simulations demonstrated that past a geometry‐dependent threshold, increased helical content improved washout, denoted by decreased regions of low velocity and low WSS. Overall, spiral flow markedly improved washout in difficult to reach areas in the angled conduits. This has several important clinical implications: spiral flow shows great promise in reducing blood‐transport‐related complications and can be used to enhance the performance of future medical devices (eg grafts, mechanical circulatory support devices, hemodialysis access ports).

     
    more » « less