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 in the case when (1) the geometry becomes complex (from a few branches to a full pulmonary artery), (2) the model becomes more complex (from fluid-only calculation to coupled fluid-structure interaction calculation), (3) both the fluid and wall models become highly nonlinear, and (4) 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 CT (Computed Tomography) 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.

Stabilized finite element methods for a blood flow model of arteriosclerosis

A Stabilized Finite Element Formulation of Non-Newtonian Fluid Model of Blood Flow in A Bifurcated Channel with Overlapping Stenosis

Mathematical and numerical analysis for PDE systems modeling intravascular drug release from arterial stents and transport in arterial tissue

Numerical and Lyapunov-Based Investigation of the Effect of Stenosis on Blood Transport Stability Using a Control-Theoretic PDE Model of Cardiovascular Flow

A non-Newtonian fluid model for blood flow through arteries under stenotic conditions

Numerical Analysis of Blood Flow Through Stenosed Microvessels Using a Multi-Phase Model

Towards a Mathematical Model for the Solidification and Rupture of Blood in Stenosed Arteries

Unsteady viscous flow model on moving the domain through a stenotic artery

Analysis of physiological pulsating flow of fractional Maxwell fluid in a locally narrow artery

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

Multifaceted simulation: Finite volume and finite element modeling of blood flow in multiple stenosed arteries

Simulation of a Pulsatile Non-Newtonian Flow Past a Stenosed 2D Artery with Atherosclerosis

Multiscale constitutive framework of 1D blood flow modeling: Asymptotic limits and numerical methods

Numerical Analysis of the Non-Newtonian Blood Flow in the Non-Planar Artery with Bifurcation

Extended one-dimensional reduced model for blood flow within a stenotic artery

Non-Newtonian and poroelastic effects in simulations of arterial flows

Simulating the Blood Flow for the Aorta with a Stenosis

Blood-Flow Modelling Along and Trough a Braided Multi-Layer Metallic Stent

Probabilistic 3D Object Recognition Based on Multiple Interpretations Generation

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