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 more » 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. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
International journal for numerical methods in biomedical engineering
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 basedmore »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.« less
  2. 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.
  3. This study presents the first 3D two-way coupled fluid structure interaction (FSI) simulation of a hybrid anechoic wind tunnel (HAWT) test section with modeling all important effects, such as turbulence, Kevlar wall porosity and deflection, and reveals for the first time the complete 3D flow structure associated with a lifting model placed into a HAWT. The Kevlar deflections are captured using finite element analysis (FEA) with shell elements operated under a membrane condition. Three-dimensional RANS CFD simulations are used to resolve the flow field. Aerodynamic experimental results are available and are compared against the FSI results. Quantitatively, the pressure coefficientsmore »on the airfoil are in good agreement with experimental results. The lift coefficient was slightly underpredicted while the drag was overpredicted by the CFD simulations. The flow structure downstream of the airfoil showed good agreement with the experiments, particularly over the wind tunnel walls where the Kevlar windows interact with the flow field. A discrepancy between previous experimental observations and juncture flow-induced vortices at the ends of the airfoil is found to stem from the limited ability of turbulence models. The qualitative behavior of the flow, including airfoil pressures and cross-sectional flow structure is well captured in the CFD. From the structural side, the behavior of the Kevlar windows and the flow developing over them is closely related to the aerodynamic pressure field induced by the airfoil. The Kevlar displacement and the transpiration velocity across the material is dominated by flow blockage effects, generated aerodynamic lift, and the wake of the airfoil. The airfoil wake increases the Kevlar window displacement, which was previously not resolved by two-dimensional panel-method simulations. The static pressure distribution over the Kevlar windows is symmetrical about the tunnel mid-height, confirming a dominantly two-dimensional flow field.« less
  4. Nonlinear fluid–structure interaction (FSI) problems on unstructured meshes in 3D appear in many applications in science and engineering, such as vibration analysis of aircrafts and patient-specific diagnosis of cardiovascular diseases. In this work, we develop a highly scalable, parallel algorithmic and software framework for FSI problems consisting of a nonlinear fluid system and a nonlinear solid system, that are coupled monolithically.
  5. Abstract Aerosol jet printing (AJP) is a direct-write additive manufacturing technique, which has emerged as a high-resolution method for the fabrication of a broad spectrum of electronic devices. Despite the advantages and critical applications of AJP in the printed-electronics industry, AJP process is intrinsically unstable, complex, and prone to unexpected gradual drifts, which adversely affect the morphology and consequently the functional performance of a printed electronic device. Therefore, in situ process monitoring and control in AJP is an inevitable need. In this respect, in addition to experimental characterization of the AJP process, physical models would be required to explain themore »underlying aerodynamic phenomena in AJP. The goal of this research work is to establish a physics-based computational platform for prediction of aerosol flow regimes and ultimately, physics-driven control of the AJP process. In pursuit of this goal, the objective is to forward a three-dimensional (3D) compressible, turbulent, multiphase computational fluid dynamics (CFD) model to investigate the aerodynamics behind: (i) aerosol generation, (ii) aerosol transport, and (iii) aerosol deposition on a moving free surface in the AJP process. The complex geometries of the deposition head as well as the pneumatic atomizer were modeled in the ansys-fluent environment, based on patented designs in addition to accurate measurements, obtained from 3D X-ray micro-computed tomography (μ-CT) imaging. The entire volume of the constructed geometries was subsequently meshed using a mixture of smooth and soft quadrilateral elements, with consideration of layers of inflation to obtain an accurate solution near the walls. A combined approach, based on the density-based and pressure-based Navier–Stokes formation, was adopted to obtain steady-state solutions and to bring the conservation imbalances below a specified linearization tolerance (i.e., 10−6). Turbulence was modeled using the realizable k-ε viscous model with scalable wall functions. A coupled two-phase flow model was, in addition, set up to track a large number of injected particles. The boundary conditions of the CFD model were defined based on experimental sensor data, recorded from the AJP control system. The accuracy of the model was validated using a factorial experiment, composed of AJ-deposition of a silver nanoparticle ink on a polyimide substrate. The outcomes of this study pave the way for the implementation of physics-driven in situ monitoring and control of AJP.« less