skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Global Sensitivity Analysis for Patient-Specific Aortic Simulations: The Role of Geometry, Boundary Condition and Large Eddy Simulation Modeling Parameters
Abstract Numerical simulations for computational hemodynamics in clinical settings require a combination of many ingredients, mathematical models, solvers and patient-specific data. The sensitivity of the solutions to these factors may be critical, particularly when we have a partial or noisy knowledge of data. Uncertainty quantification is crucial to assess the reliability of the results. We present here an extensive sensitivity analysis in aortic flow simulations, to quantify the dependence of clinically relevant quantities to the patient-specific geometry and the inflow boundary conditions. Geometry and inflow conditions are generally believed to have a major impact on numerical simulations. We resort to a global sensitivity analysis, (i.e., not restricted to a linearization around a working point), based on polynomial chaos expansion (PCE) and the associated Sobol' indices. We regard the geometry and the inflow conditions as the realization of a parametric stochastic process. To construct a physically consistent stochastic process for the geometry, we use a set of longitudinal-in-time images of a patient with an abdominal aortic aneurysm (AAA) to parametrize geometrical variations. Aortic flow is highly disturbed during systole. This leads to high computational costs, even amplified in a sensitivity analysis -when many simulations are needed. To mitigate this, we consider here a large Eddy simulation (LES) model. Our model depends in particular on a user-defined parameter called filter radius. We borrowed the tools of the global sensitivity analysis to assess the sensitivity of the solution to this parameter too. The targeted quantities of interest (QoI) include: the total kinetic energy (TKE), the time-average wall shear stress (TAWSS), and the oscillatory shear index (OSI). The results show that these indexes are mostly sensitive to the geometry. Also, we find that the sensitivity may be different during different instants of the heartbeat and in different regions of the domain of interest. This analysis helps to assess the reliability of in silico tools for clinical applications.  more » « less
Award ID(s):
1620406
PAR ID:
10232750
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Biomechanical Engineering
Volume:
143
Issue:
2
ISSN:
0148-0731
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Transcatheter aortic valve replacement (TAVR) first received FDA approval for high-risk surgical patients in 2011 and has been approved for low-risk surgical patients since 2019. It is now the most common type of aortic valve replacement, and its use continues to accelerate. Computer modeling and simulation (CM&S) is a tool to aid in TAVR device design, regulatory approval, and indication in patient-specific care. This study introduces a computational fluid-structure interaction (FSI) model of TAVR with Medtronic’s CoreValve Evolut R device using the immersed finite element-difference (IFED) method. We perform dynamic simulations of crimping and deployment of the Evolut R, as well as device behavior across the cardiac cycle in a patient-specific aortic root anatomy reconstructed from computed tomography (CT) image data. These IFED simulations, which incorporate biomechanics models fit to experimental tensile test data, automatically capture the contact within the device and between the self-expanding stent and native anatomy. Further, we apply realistic driving and loading conditions based on clinical measurements of human ventricular and aortic pressures and flow rates to demonstrate that our Evolut R model supports a physiological diastolic pressure load and provides informative clinical performance predictions. 
    more » « less
  2. Heart valve disease (HVD), a significant cardiovascular complication, is one of the leading global causes of morbidity and mortality. Treatment for HVD often involves medical devices such as bioprosthetic valves. However, the design and optimization of these devices require a thorough understanding of their biomechanical and hemodynamic interactions with patient-specific anatomical structures. Parametric procedural geometry has become a powerful tool in enhancing the efficiency and accuracy of design optimization for such devices, allowing researchers to systematically explore a wide range of possible configurations. In this work, we present a robust framework for parametric and procedural modeling of stented bioprosthetic heart valves and patient-specific aortic geometries. The framework employs non-uniform rational B-splines (NURBS)-based geometric parameterization, enabling precise control over key anatomical and design variables. By enabling a modular and expandable workflow, the framework supports iterative optimization of valve designs to achieve improved hemodynamic performance and durability. We demonstrate its applicability through simulations on bioprosthetic aortic valves, highlighting the impact of geometric parameters on valve function and their potential for personalized device design. By coupling parametric geometry with computational tools, this framework offers researchers and engineers a streamlined pathway toward innovative and patient-specific cardiovascular solutions. 
    more » « less
  3. Carotid artery diseases, such as atherosclerosis, are a major cause of death in the United States. Wall shear stresses are known to prompt plaque formation, but there is limited understanding of the complex flow structures underlying these stresses and how they differ in a pre-disposed high-risk patient cohort. A ‘healthy’ and a novel ‘pre-disposed’ carotid artery bifurcation model was determined based on patient-averaged clinical data, where the ‘pre-disposed’ model represents a pathological anatomy. Computational fluid dynamic simulations were performed using a physiological flow based on healthy human subjects. A main hairpin vortical structure in the internal carotid artery sinus was observed, which locally increased instantaneous wall shear stress. In the pre-disposed geometry, this vortical structure starts at an earlier instance in the cardiac flow cycle and persists over a much shorter period, where the second half of the cardiac cycle is dominated by perturbed secondary flow structures and vortices. This coincides with weaker favorable axial pressure gradient peaks over the sinus for the ‘pre-disposed’ geometry. The findings reveal a strong correlation between vortical structures and wall shear stress and imply that an intact internal carotid artery sinus hairpin vortical structure has a physiologically beneficial role by increasing local wall shear stresses. The deterioration of this beneficial vortical structure is expected to play a significant role in atherosclerotic plaque formation. 
    more » « less
  4. Pamies, Pep (Ed.)
    Preclinical models of aortic stenosis can induce left ventricular pressure overload and coarsely control the severity of aortic constriction. However, they do not recapitulate the haemodynamics and flow patterns associated with the disease. Here we report the development of a customizable soft robotic aortic sleeve that can mimic the haemodynamics and biomechanics of aortic stenosis. By allowing for the adjustment of actuation patterns and blood-flow dynamics, the robotic sleeve recapitulates clinically relevant haemodynamics in a porcine model of aortic stenosis, as we show via in vivo echocardiography and catheterization studies, and a combination of in vitro and computational analyses. Using in vivo and in vitro magnetic resonance imaging, we also quantified the four-dimensional blood-flow velocity profiles associated with the disease and with bicommissural and unicommissural defects re-created by the robotic sleeve. The design of the sleeve, which can be adjusted on the basis of computed tomography data, allows for the design of patient-specific devices that may guide clinical decisions and improve the management and treatment of patients with aortic stenosis. 
    more » « less
  5. Abstract This study evaluates a hypothesis for the role of vertical wind shear in deep convection initiation (DCI) that was introduced in Part I by examining behavior of a series of numerical simulations. The hypothesis states, “Initial moist updrafts that exceed a width and shear threshold will ‘root’ within a progressively deeper steering current with time, increase their low-level cloud-relative flow and inflow, widen, and subsequently reduce their susceptibility to entrainment-driven dilution, evolving toward a quasi-steady self-sustaining state.” A theoretical model that embodied key elements of the hypothesis was developed in Part I, and the behavior of this model was explored within a multidimensional environmental parameter space. Remarkably similar behavior is evident in the simulations studied here to that of the theoretical model, both in terms of the temporal evolution of DCI and in the sensitivity of DCI to environmental parameters. Notably, both the simulations and theoretical model experience a bifurcation in outcomes, whereby nascent clouds that are narrower than a given initial radiusR0threshold quickly decay and those above theR0threshold undergo DCI. An important assumption in the theoretical model, which states that the cloud-relative flow of the background environmentVCRdetermines cloud radiusR, is scrutinized in the simulations. It is shown that storm-induced inflow is small relative toVCRbeyond a few kilometers from the updraft edge, andVCRtherefore plays a predominant role in transporting conditionally unstable air to the updraft. Thus, the critical role ofVCRin determiningRis validated. 
    more » « less