skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Patient‐specific CFD modelling in the thoracic aorta with PC‐MRI–based boundary conditions: A least‐square three‐element Windkessel approach
Title: Patient‐specific CFD modelling in the thoracic aorta with PC‐MRI–based boundary conditions: A least‐square three‐element Windkessel approach
Abstract

The increasing use of computational fluid dynamics for simulating blood flow in clinics demands the identification of appropriate patient‐specific boundary conditions for the customization of the mathematical models. These conditions should ideally be retrieved from measurements. However, finite resolution of devices as well as other practical/ethical reasons prevent the construction of complete data sets necessary to make the mathematical problems well posed. Available data need to be completed by modelling assumptions, whose impact on the final solution has to be carefully addressed. Focusing on aortic vascular districts and related pathologies, we present here a method for efficiently and robustly prescribing phase contrast MRI–based patient‐specific data as boundary conditions at the domain of interest. In particular, for the outlets, the basic idea is to obtain pressure conditions from an appropriate elaboration of available flow rates on the basis of a 3D/0D dimensionally heterogeneous modelling. The key point is that the parameters are obtained by a constrained optimization procedure. The rationale is that pressure conditions have a reduced impact on the numerical solution compared with velocity conditions, yielding a simulation framework less exposed to noise and inconsistency of the data, as well as to the arbitrariness of the underlying modelling assumptions. Numerical results confirm the reliability of the approach in comparison with other patient‐specific approaches adopted in the literature.

  more » « less
NSF-PAR ID:
10073841
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical Methods in Biomedical Engineering
Volume:
34
Issue:
11
ISSN:
2040-7939
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Journal of Biomechanical Engineering)
    null (Ed.)
    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
  2. ; ; ; ; ( , Geophysical Journal International)
    SUMMARY

    Tsunami generation by offshore earthquakes is a problem of scientific interest and practical relevance, and one that requires numerical modelling for data interpretation and hazard assessment. Most numerical models utilize two-step methods with one-way coupling between separate earthquake and tsunami models, based on approximations that might limit the applicability and accuracy of the resulting solution. In particular, standard methods focus exclusively on tsunami wave modelling, neglecting larger amplitude ocean acoustic and seismic waves that are superimposed on tsunami waves in the source region. In this study, we compare four earthquake-tsunami modelling methods. We identify dimensionless parameters to quantitatively approximate dominant wave modes in the earthquake-tsunami source region, highlighting how the method assumptions affect the results and discuss which methods are appropriate for various applications such as interpretation of data from offshore instruments in the source region. Most methods couple a 3-D solid earth model, which provides the seismic wavefield or at least the static elastic displacements, with a 2-D depth-averaged shallow water tsunami model. Assuming the ocean is incompressible and tsunami propagation is negligible over the earthquake duration leads to the instantaneous source method, which equates the static earthquake seafloor uplift with the initial tsunami sea surface height. For longer duration earthquakes, it is appropriate to follow the time-dependent source method, which uses time-dependent earthquake seafloor velocity as a forcing term in the tsunami mass balance. Neither method captures ocean acoustic or seismic waves, motivating more advanced methods that capture the full wavefield. The superposition method of Saito et al. solves the 3-D elastic and acoustic equations to model the seismic wavefield and response of a compressible ocean without gravity. Then, changes in sea surface height from the zero-gravity solution are used as a forcing term in a separate tsunami simulation, typically run with a shallow water solver. A superposition of the earthquake and tsunami solutions provides an approximation to the complete wavefield. This method is algorithmically a two-step method. The complete wavefield is captured in the fully coupled method, which utilizes a coupled solid Earth and compressible ocean model with gravity. The fully coupled method, recently incorporated into the 3-D open-source code SeisSol, simultaneously solves earthquake rupture, seismic waves and ocean response (including gravity). We show that the superposition method emerges as an approximation to the fully coupled method subject to often well-justified assumptions. Furthermore, using the fully coupled method, we examine how the source spectrum and ocean depth influence the expression of oceanic Rayleigh waves. Understanding the range of validity of each method, as well as its computational expense, facilitates the selection of modelling methods for the accurate assessment of earthquake and tsunami hazards and the interpretation of data from offshore instruments.

     
    more » « less
  3. ; ; ; ; ( , Journal of Fluid Mechanics)
    null (Ed.)
    This paper investigates the steady axisymmetric structure of the cold boundary-layer flow surrounding fire whirls developing over localized fuel sources lying on a horizontal surface. The inviscid swirling motion found outside the boundary layer, driven by the entrainment of the buoyant turbulent plume of hot combustion products that develops above the fire, is described by an irrotational solution, obtained by combining Taylor's self-similar solution for the motion in the axial plane with the azimuthal motion induced by a line vortex of circulation $2 {\rm \pi}\Gamma$ . The development of the boundary layer from a prescribed radial location is determined by numerical integration for different swirl levels, measured by the value of the radial-to-azimuthal velocity ratio $\sigma$ at the initial radial location. As in the case $\sigma =0$ , treated in the seminal boundary-layer analysis of Burggraf et al. ( Phys. Fluids , vol. 14, 1971, pp. 1821–1833), the pressure gradient associated with the centripetal acceleration of the inviscid flow is seen to generate a pronounced radial inflow. Specific attention is given to the terminal shape of the boundary-layer velocity near the axis, which displays a three-layered structure that is described by matched asymptotic expansions. The resulting composite expansion, dependent on the level of ambient swirl through the parameter $\sigma$ , is employed as boundary condition to describe the deflection of the boundary-layer flow near the axis to form a vertical swirl jet. Numerical solutions of the resulting non-slender collision region for different values of $\sigma$ are presented both for inviscid flow and for viscous flow with moderately large values of the controlling Reynolds number $\Gamma /\nu$ . The velocity description provided is useful in mathematical formulations of localized fire-whirl flows, providing consistent boundary conditions accounting for the ambient swirl level. 
    more » « less
  4. ; ( , ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik)
    Abstract

    A flow vessel with an elastic wall can deform significantly due to viscous fluid flow within it, even at vanishing Reynolds number (no fluid inertia). Deformation leads to an enhancement of throughput due to the change in cross‐sectional area. The latter gives rise to a non‐constant pressure gradient in the flow‐wise direction and, hence, to a nonlinear flow rate–pressure drop relation (unlike the Hagen–Poiseuille law for a rigid tube). Many biofluids are non‐Newtonian, and are well approximated by generalized Newtonian (say, power‐law) rheological models. Consequently, we analyze the problem of steady low Reynolds number flow of a generalized Newtonian fluid through a slender elastic tube by coupling fluid lubrication theory to a structural problem posed in terms of Donnell shell theory. A perturbative approach (in the slenderness parameter) yields analytical solutions for both the flow and the deformation. Using matched asymptotics, we obtain a uniformly valid solution for the tube's radial displacement, which features both a boundary layer and a corner layer caused by localized bending near the clamped ends. In doing so, we obtain a “generalized Hagen–Poiseuille law” for soft microtubes. We benchmark the mathematical predictions against three‐dimensional two‐way coupled direct numerical simulations (DNS) of flow and deformation performed using the commercial computational engineering platform by ANSYS. The simulations show good agreement and establish the range of validity of the theory. Finally, we discuss the implications of the theory on the problem of the flow‐induced deformation of a blood vessel, which is featured in some textbooks.

     
    more » « less
  5. ; ( , Monthly Weather Review)
    Abstract

    We present a new ensemble of 36 numerical experiments aimed at comprehensively gauging the sensitivity of nested large-eddy simulations (LES) driven by large-scale dynamics. Specifically, we explore 36 multiscale configurations of the Weather Research and Forecasting (WRF) Model to simulate the boundary layer flow over the complex topography at the Perdigão field site, with five nested domains discretized at horizontal resolutions ranging from 11.25 km to 30 m. Each ensemble member has a unique combination of the following input factors: (i) large-scale initial and boundary conditions, (ii) subgrid turbulence modeling in thegray zoneof turbulence, (iii) subgrid-scale (SGS) models in LES, and (iv) topography and land-cover datasets. We probe their relative importance for LES calculations of velocity, temperature, and moisture fields. Variance decomposition analysis unravels large sensitivities to topography and land-use datasets and very weak sensitivity to the LES SGS model. Discrepancies within ensemble members can be as large as 2.5 m s−1for the time-averaged near-surface wind speed on the ridge and as large as 10 m s−1without time averaging. At specific time points, a large fraction of this sensitivity can be explained by the different turbulence models in the gray zone domains. We implement a horizontal momentum and moisture budget routine in WRF to further elucidate the mechanisms behind the observed sensitivity, paving the way for an increased understanding of the tangible effects of the gray zone of turbulence problem.

    Significance Statement

    Several science and engineering applications, including wind turbine siting and operations, weather prediction, and downscaling of climate projections, call for high-resolution numerical simulations of the lowest part of the atmosphere. Recent studies have highlighted that such high-resolution simulations, coupled with large-scale models, are challenging and require several important assumptions. With a new set of numerical experiments, we evaluate and compare the significance of different assumptions and outstanding challenges in multiscale modeling (i.e., coupling large-scale models and high-resolution atmospheric simulations). The ultimate goal of this analysis is to put each individual assumption into the wider perspective of a realistic problem and quantify its relative importance compared to other important modeling choices.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email