Global Sensitivity Analysis for Patient-Specific Aortic Simulations: The Role of Geometry, Boundary Condition and Large Eddy Simulation Modeling Parameters

2020 ◽  
Vol 143 (2) ◽  
Author(s):  
Huijuan Xu ◽  
Davide Baroli ◽  
Alessandro Veneziani
Keyword(s):  
Stochastic Process ◽  
Sensitivity Analysis ◽  
Large Eddy Simulation ◽  
Numerical Simulations ◽  
Global Sensitivity Analysis ◽  
Patient Specific ◽  
Global Sensitivity ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Inflow Conditions

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.

Large eddy simulation of multi-regime burner: a reaction mechanism sensitivity analysis

2022 ◽  
Author(s):  
Lorenzo Angelilli ◽  
Pietro Paolo Ciottoli ◽  
Francisco E. Hernandez Perez ◽  
Mauro Valorani ◽  
Hong G. Im ◽  
...  
Keyword(s):  
Sensitivity Analysis ◽  
Reaction Mechanism ◽  
Large Eddy Simulation ◽  
Eddy Simulation ◽  
Large Eddy

scholarly journals Developments in turbulence research: a review based on the 1999 Programme of the Isaac Newton Institute, Cambridge

2001 ◽  
Vol 436 ◽  
pp. 353-391 ◽  
Author(s):  
J. C. R. HUNT ◽  
N. D. SANDHAM ◽  
J. C. VASSILICOS ◽  
B. E. LAUNDER ◽  
P. A. MONKEWITZ ◽  
...  
Keyword(s):  
Reynolds Number ◽  
Large Eddy Simulation ◽  
Numerical Simulations ◽  
Turbulent Flows ◽  
Small Scale ◽  
Isaac Newton ◽  
Eddy Simulation ◽  
Isaac Newton Institute ◽  
Large Eddy ◽  
Eddy Structures

Recent research is making progress in framing more precisely the basic dynamical and statistical questions about turbulence and in answering them. It is helping both to define the likely limits to current methods for modelling industrial and environmental turbulent flows, and to suggest new approaches to overcome these limitations. Our selective review is based on the themes and new results that emerged from more than 300 presentations during the Programme held in 1999 at the Isaac Newton Institute, Cambridge, UK, and on research reported elsewhere. A general conclusion is that, although turbulence is not a universal state of nature, there are certain statistical measures and kinematic features of the small-scale flow field that occur in most turbulent flows, while the large-scale eddy motions have qualitative similarities within particular types of turbulence defined by the mean flow, initial or boundary conditions, and in some cases, the range of Reynolds numbers involved. The forced transition to turbulence of laminar flows caused by strong external disturbances was shown to be highly dependent on their amplitude, location, and the type of flow. Global and elliptical instabilities explain much of the three-dimensional and sudden nature of the transition phenomena. A review of experimental results shows how the structure of turbulence, especially in shear flows, continues to change as the Reynolds number of the turbulence increases well above about 104 in ways that current numerical simulations cannot reproduce. Studies of the dynamics of small eddy structures and their mutual interactions indicate that there is a set of characteristic mechanisms in which vortices develop (vortex stretching, roll-up of instability sheets, formation of vortex tubes) and another set in which they break up (through instabilities and self- destructive interactions). Numerical simulations and theoretical arguments suggest that these often occur sequentially in randomly occurring cycles. The factors that determine the overall spectrum of turbulence were reviewed. For a narrow distribution of eddy scales, the form of the spectrum can be defined by characteristic forms of individual eddies. However, if the distribution covers a wide range of scales (as in elongated eddies in the ‘wall’ layer of turbulent boundary layers), they collectively determine the spectra (as assumed in classical theory). Mathematical analyses of the Navier–Stokes and Euler equations applied to eddy structures lead to certain limits being defined regarding the tendencies of the vorticity field to become infinitely large locally. Approximate solutions for eigen modes and Fourier components reveal striking features of the temporal, near-wall structure such as bursting, and of the very elongated, spatial spectra of sheared inhomogeneous turbulence; but other kinds of eddy concepts are needed in less structured parts of the turbulence. Renormalized perturbation methods can now calculate consistently, and in good agreement with experiment, the evolution of second- and third-order spectra of homogeneous and isotropic turbulence. The fact that these calculations do not explicitly include high-order moments and extreme events, suggests that they may play a minor role in the basic dynamics. New methods of approximate numerical simulations of the larger scales of turbulence or ‘very large eddy simulation’ (VLES) based on using statistical models for the smaller scales (as is common in meteorological modelling) enable some turbulent flows with a non-local and non-equilibrium structure, such as impinging or convective flows, to be calculated more efficiently than by using large eddy simulation (LES), and more accurately than by using ‘engineering’ models for statistics at a single point. Generally it is shown that where the turbulence in a fluid volume is changing rapidly and is very inhomogeneous there are flows where even the most complex ‘engineering’ Reynolds stress transport models are only satisfactory with some special adaptation; this may entail the use of transport equations for the third moments or non-universal modelling methods designed explicitly for particular types of flow. LES methods may also need flow-specific corrections for accurate modelling of different types of very high Reynolds number turbulent flow including those near rigid surfaces.This paper is dedicated to the memory of George Batchelor who was the inspiration of so much research in turbulence and who died on 30th March 2000. These results were presented at the last fluid mechanics seminar in DAMTP Cambridge that he attended in November 1999.

1831 Large-Eddy Simulation of Atmospheric Flow by MSSG-A : Parameter Sensitivity Analysis

2008 ◽  
Vol 2008.2 (0) ◽  
pp. 181-182
Author(s):  
Shingo KANAZAWA ◽  
Takeshi KAMIO ◽  
Takeshi SUGIMURA ◽  
Makoto IIDA ◽  
Keiko TAKAHASHI ◽  
...  
Keyword(s):  
Sensitivity Analysis ◽  
Large Eddy Simulation ◽  
Parameter Sensitivity ◽  
Parameter Sensitivity Analysis ◽  
Atmospheric Flow ◽  
Eddy Simulation ◽  
Large Eddy

Developing large eddy simulation for turbomachinery applications

2009 ◽  
Vol 367 (1899) ◽  
pp. 2999-3013 ◽  
Author(s):  
Simon J. Eastwood ◽  
Paul G. Tucker ◽  
Hao Xia ◽  
Christian Klostermeier
Keyword(s):  
Large Eddy Simulation ◽  
Leading Edge ◽  
Problem Definition ◽  
Navier Stokes ◽  
Complex Geometries ◽  
Subgrid Scale ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Chevron Nozzle ◽  
Inflow Conditions

For jets, large eddy resolving simulations are compared for a range of numerical schemes with no subgrid scale (SGS) model and for a range of SGS models with the same scheme. There is little variation in results for the different SGS models, and it is shown that, for schemes which tend towards having dissipative elements, the SGS model can be abandoned, giving what can be termed numerical large eddy simulation (NLES). More complex geometries are investigated, including coaxial and chevron nozzle jets. A near-wall Reynolds-averaged Navier–Stokes (RANS) model is used to cover over streak-like structures that cannot be resolved. Compressor and turbine flows are also successfully computed using a similar NLES–RANS strategy. Upstream of the compressor leading edge, the RANS layer is helpful in preventing premature separation. Capturing the correct flow over the turbine is particularly challenging, but nonetheless the RANS layer is helpful. In relation to the SGS model, for the flows considered, evidence suggests issues such as inflow conditions, problem definition and transition are more influential.

Partially Averaged Navier-Stokes Turbulence Modeling of Flow in 5x5 PWR Fuel Assembly With Spacer Grid

2018 ◽  
Author(s):  
Giacomo Busco ◽  
Yassin A. Hassan
Keyword(s):  
Large Eddy Simulation ◽  
Numerical Simulations ◽  
Stokes Equations ◽  
Direct Numerical Simulations ◽  
Navier Stokes ◽  
Spacer Grid ◽  
Navier Stokes Equations ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Reynolds Averaged Navier Stokes

The highly turbulent flow inside a pressurized water reactor makes unpractical the use of scale resolving simulations, due to the large number of space and time turbulent structures. The high computational cost associated with typical large eddies simulations or direct numerical simulations techniques is unsuitable due to the large spatiotemporal resolution required. Partially averaged Navier-Stokes turbulence model is presented as bridging model between Reynolds averaged Navier-Stokes equations and direct numerical simulations. As filtered representation of the Navier-Stokes equations, the model is able to continuously shift its energy-based filter, inside the turbulence spectrum, being able to resolve the turbulent scales of interest. The choice of energy based cut-off filters gives the chance to directly impose the degree of needed resolution, where the most important large scales unsteadiness are resolved at minimal computational expenses. The partially averaged Navier-Stokes modelling approach has been tested for a Reynolds number of 14,000, inside a 5 × 5 fuel bundle, with a single spacer grid and split-type mixing vanes. Four different filters have been tested, whose resolution ranged from Reynolds averaged Navier-Stokes and large eddy simulation. A comparison with large eddy simulation will be presented. The results show that the partially averaged Navier-Stokes modeling produces results comparable to those of large eddy simulation when the appropriate cut-off energy filter is chosen. The turbulence models results will be compared with the available particle image velocimetry experimental data.

Sign in / Sign up

Export Citation Format

Share Document