Fluid-structure interaction modeling of blood flow in the pulmonary arteries using the unified continuum and variational multiscale formulation

Ju, Liu; Yang, Weiguang; Lan, Ingrid S.; Marsden, Alison L.

doi:10.1016/j.mechrescom.2020.103556

Cited by 36 publications

(33 citation statements)

References 34 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…We also note that ∞ = 0.5 is commonly adopted as a balanced choice for CFD and FSI investigations. 6,20,32,35,36 Remark 2. One may introduce an additional set of parameters p m , p f , and p for the pressure variable (see p. 89 of Reference 11).…”

Section: The Generalized-schemes Under Comparisonmentioning

confidence: 99%

See 1 more Smart Citation

A note on the accuracy of the generalized‐α scheme for the incompressible Navier‐Stokes equations

Lan

Tikenogullari

et al. 2020

Numerical Meth Engineering

Self Cite

View full text Add to dashboard Cite

We investigate the temporal accuracy of two generalized-schemes for the incompressible Navier-Stokes equations. In a widely-adopted approach, the pressure is collocated at the time step t n + 1 while the remainder of the Navier-Stokes equations is discretized following the generalized-scheme. That scheme has been claimed to be second-order accurate in time. We developed a suite of numerical code using inf-sup stable higher-order non-uniform rational B-spline (NURBS) elements for spatial discretization. In doing so, we are able to achieve high spatial accuracy and to investigate asymptotic temporal convergence behavior. Numerical evidence suggests that only first-order accuracy is achieved, at least for the pressure, in this aforesaid temporal discretization approach. On the other hand, evaluating the pressure at the intermediate time step t n+ f recovers second-order accuracy, and the numerical implementation is simplified. We recommend this second approach as the generalized-scheme of choice when integrating the incompressible Navier-Stokes equations.

show abstract

Section: The Generalized-schemes Under Comparisonmentioning

confidence: 99%

“…It sometimes yields favorable numerical performance in flow simulations 3 and presents the only scenario in which Scheme‐1 and Scheme‐2 become identical. We also note that ϱ ∞ = 0.5 is commonly adopted as a balanced choice for CFD and FSI investigations 6,20,32,35,36 …”

Section: Governing Equations and The Spatiotemporal Discretizationmentioning

confidence: 99%

A note on the accuracy of the generalized‐α scheme for the incompressible Navier‐Stokes equations

Lan

Tikenogullari

et al. 2020

Numerical Meth Engineering

Self Cite

View full text Add to dashboard Cite

show abstract

“…This limitation was resolved by coupling the hemodynamics with enclosing wall motion using fluid-structure interaction (FSI) simulations. 16,28,29,38,39 One FSI method involved a coupledmomentum model with the assumption of linearized stiffness for the arterial wall, 35,38,39 while Gutierrez et al 9 used small deformation theory. Su et al 28,29 employed idealized geometries.…”

Section: Introductionmentioning

confidence: 99%

“…A vast majority of modeling studies 5,10,11,16,21,38,39 have dealt with pulmonary arterial hypertension (PAH, WHO group I). However, left heart disease (WHO group II) and lung disease (WHO group III) are causative factors for PH in the majority of diagnosed cases.…”

Section: Introductionmentioning

confidence: 99%

Patient-Specific Computational Analysis of Hemodynamics in Adult Pulmonary Hypertension

et al. 2021

View full text Add to dashboard Cite

Pulmonary hypertension (PH) is a progressive disease characterized by elevated pressure and vascular resistance in the pulmonary arteries. Nearly 250,000 hospitalizations occur annually in the US with PH as the primary or secondary condition. A definitive diagnosis of PH requires right heart catheterization (RHC) in addition to a chest computed tomography, a walking test, and others. While RHC is the gold standard for diagnosing PH, it is invasive and posseses inherent risks and contraindications. In this work, we characterized the patient-specific pulmonary hemodynamics in silico for diverse PH WHO groups. We grouped patients on the basis of mean pulmonary arterial pressure (mPAP) into three disease severity groups: at-risk (18mmHg mPAP<25mmHg, denoted with A), mild (25mmHg mPAP<40mmHg, denoted with M), and severe (mPAP ! 40mmHg, denoted with S). The pulsatile flow hemodynamics was simulated by evaluating the three-dimensional Navier-Stokes system of equations using a flow solver developed by customizing OpenFOAM libraries (v5.0, The OpenFOAM Foundation). Quasi patient-specific boundary conditions were implemented using a Womersley inlet velocity profile and transient resistance outflow conditions. Hemodynamic indices such as spatially averaged wall shear stress (SAWSS), wall shear stress gradient (WSSG), timeaveraged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT), were evaluated along with the clinical metrics pulmonary vascular resistance (PVR), stroke volume (SV) and compliance (C), to assess possible spatiotemporal correlations. We observed statistically significant decreases in SAWSS, WSSG, and TAWSS, and increases in OSI and RRT with disease severity. PVR was moderately correlated with SAWSS and RRT at the mid-notch stage of the cardiac cycle when these indices were computed using the global pulmonary arterial geometry. These results are promising in the context of a long-term goal of identifying computational biomarkers that can serve as surrogates for invasive diagnostic protocols of PH.

show abstract

“…Moreover, it has been shown that FSI simulations seem to accurately capture flow fields over a deformable arterial wall with patient-specific geometries, whereas models based on rigid arterial walls seem to overestimate flow velocities and WSS values. 37,38 Therefore, an FSI model has been used to simulate a patient-specific pulmonary artery, 39 aortic dissection, 40 coronary artery, 37 and carotid artery. 41 Specifically, Liu et al used MRI (magnetic resonance imaging) data to construct the input geometry and captured the vessel wall deformation following hyper-elastic constitutive law.…”

Section: Introductionmentioning

confidence: 99%

Direct numerical simulation of a pulsatile flow in a stenotic channel using immersed boundary method

Mirfendereski

Park

2021

Engineering Reports

View full text Add to dashboard Cite

A three-dimensional direct numerical simulation model coupled with the immersed boundary method has been developed to simulate a pulsatile flow in a planar channel with single and double one-sided semicircular constrictions. For relevance to blood flow in large arteries, simulations have been performed at Reynolds numbers of 750 and 1000. Flow physics and resultant wall shear stress (WSS)-based hemodynamic parameters are presented. The instantaneous vortex dynamics, mean flow characteristics, and turbulent energy spectra are evaluated for flow physics. Subsequently, three WSS-based parameters, namely the time-averaged WSS, oscillatory shear index, and relative residence time, are calculated over the stenotic wall and correlated with flow physics to identify the regions prone to atherosclerotic plaque progression. Results show that the double stenotic channel leads to high-intensity and broadband turbulent characteristics downstream, promoting critical values of the WSS-based parameters in the post-stenotic areas. In addition, the inter-space area between two stenoses displays multiple strong recirculations, making this area highly prone to atherosclerosis progression. The effect of stenosis degree on the WSS-based parameters is studied up to 60% degree. As the degree of occlusion is increased, larger regions are involved with the nonphysiological ranges of the WSS-based parameters.

show abstract

Fluid-structure interaction modeling of blood flow in the pulmonary arteries using the unified continuum and variational multiscale formulation

Cited by 36 publications

References 34 publications

A note on the accuracy of the generalized‐α scheme for the incompressible Navier‐Stokes equations

A note on the accuracy of the generalized‐α scheme for the incompressible Navier‐Stokes equations

Patient-Specific Computational Analysis of Hemodynamics in Adult Pulmonary Hypertension

Direct numerical simulation of a pulsatile flow in a stenotic channel using immersed boundary method

Contact Info

Product

Resources

About