-
Machine-learning-assisted Blending of Data-Driven Turbulence Models
Authors:
Mourad Oulghelou,
Soufiane Cherroud,
Xavier Merle,
Paola Cinnella
Abstract:
A machine learning-based methodology for blending data-driven turbulent closures for the Reynolds-Averaged Navier-Stokes (RANS) equations is proposed to improve the generalizability across different flow scenarios. Data-driven models based on sparse Bayesian learning and symbolic regression are pre-trained for specific flow classes, including turbulent channel flows, separated flows, and axisymmet…
▽ More
A machine learning-based methodology for blending data-driven turbulent closures for the Reynolds-Averaged Navier-Stokes (RANS) equations is proposed to improve the generalizability across different flow scenarios. Data-driven models based on sparse Bayesian learning and symbolic regression are pre-trained for specific flow classes, including turbulent channel flows, separated flows, and axisymmetric jets. The specialized models (called "\textit{experts}") are then blended using a set of weighting functions that reflect the local likelihood of a model to capture the reference data. The weighting functions are expressed as functions of a set of local flow features, and Random Forest regressors (RFR) are trained alongside the expert models to learn the mapping between the features and the weights. The training dataset includes additional data for flows ranging from shear flows at equilibrium to separated flows. The resulting blended model is then used to predict unseen flows, including a turbulent zero-pressure-gradient flat plate and a wall-mounted hump involving attached and separated boundary layer regions. The results show that the blending strategy adapts to the local flow conditions and effectively leverages the strengths of individual models, promoting the best-performing models locally, in a manner consistent with the training scenarios. The blended model provides a robust approximation for unseen flows that exhibit similar underlying physical characteristics to the training cases, while also successfully extrapolating to flows with different geometries and physical configurations.
△ Less
Submitted 18 October, 2024;
originally announced October 2024.
-
NeurIPS 2024 ML4CFD Competition: Harnessing Machine Learning for Computational Fluid Dynamics in Airfoil Design
Authors:
Mouadh Yagoubi,
David Danan,
Milad Leyli-abadi,
Jean-Patrick Brunet,
Jocelyn Ahmed Mazari,
Florent Bonnet,
maroua gmati,
Asma Farjallah,
Paola Cinnella,
Patrick Gallinari,
Marc Schoenauer
Abstract:
The integration of machine learning (ML) techniques for addressing intricate physics problems is increasingly recognized as a promising avenue for expediting simulations. However, assessing ML-derived physical models poses a significant challenge for their adoption within industrial contexts. This competition is designed to promote the development of innovative ML approaches for tackling physical…
▽ More
The integration of machine learning (ML) techniques for addressing intricate physics problems is increasingly recognized as a promising avenue for expediting simulations. However, assessing ML-derived physical models poses a significant challenge for their adoption within industrial contexts. This competition is designed to promote the development of innovative ML approaches for tackling physical challenges, leveraging our recently introduced unified evaluation framework known as Learning Industrial Physical Simulations (LIPS). Building upon the preliminary edition held from November 2023 to March 2024, this iteration centers on a task fundamental to a well-established physical application: airfoil design simulation, utilizing our proposed AirfRANS dataset. The competition evaluates solutions based on various criteria encompassing ML accuracy, computational efficiency, Out-Of-Distribution performance, and adherence to physical principles. Notably, this competition represents a pioneering effort in exploring ML-driven surrogate methods aimed at optimizing the trade-off between computational efficiency and accuracy in physical simulations. Hosted on the Codabench platform, the competition offers online training and evaluation for all participating solutions.
△ Less
Submitted 30 June, 2024;
originally announced July 2024.
-
Opportunities for machine learning in scientific discovery
Authors:
Ricardo Vinuesa,
Jean Rabault,
Hossein Azizpour,
Stefan Bauer,
Bingni W. Brunton,
Arne Elofsson,
Elias Jarlebring,
Hedvig Kjellstrom,
Stefano Markidis,
David Marlevi,
Paola Cinnella,
Steven L. Brunton
Abstract:
Technological advancements have substantially increased computational power and data availability, enabling the application of powerful machine-learning (ML) techniques across various fields. However, our ability to leverage ML methods for scientific discovery, {\it i.e.} to obtain fundamental and formalized knowledge about natural processes, is still in its infancy. In this review, we explore how…
▽ More
Technological advancements have substantially increased computational power and data availability, enabling the application of powerful machine-learning (ML) techniques across various fields. However, our ability to leverage ML methods for scientific discovery, {\it i.e.} to obtain fundamental and formalized knowledge about natural processes, is still in its infancy. In this review, we explore how the scientific community can increasingly leverage ML techniques to achieve scientific discoveries. We observe that the applicability and opportunity of ML depends strongly on the nature of the problem domain, and whether we have full ({\it e.g.}, turbulence), partial ({\it e.g.}, computational biochemistry), or no ({\it e.g.}, neuroscience) {\it a-priori} knowledge about the governing equations and physical properties of the system. Although challenges remain, principled use of ML is opening up new avenues for fundamental scientific discoveries. Throughout these diverse fields, there is a theme that ML is enabling researchers to embrace complexity in observational data that was previously intractable to classic analysis and numerical investigations.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Data-driven turbulence modeling
Authors:
Paola Cinnella
Abstract:
This chapter provides an introduction to data-driven techniques for the development and calibration of closure models for the Reynolds-Averaged Navier--Stokes (RANS) equations. RANS models are the workhorse for engineering applications of computational fluid dynamics (CFD) and are expected to play an important role for decades to come. However, RANS model inadequacies for complex, non-equilibrium…
▽ More
This chapter provides an introduction to data-driven techniques for the development and calibration of closure models for the Reynolds-Averaged Navier--Stokes (RANS) equations. RANS models are the workhorse for engineering applications of computational fluid dynamics (CFD) and are expected to play an important role for decades to come. However, RANS model inadequacies for complex, non-equilibrium flows and uncertainties in modeling assumptions and calibration data are still a major obstacle to the predictive capability of RANS simulations. In the following, we briefly recall the origin and limitations of RANS models, and then review their shortcomings and uncertainties. Then, we provide an introduction to data-driven approaches to RANS turbulence modeling. The latter can range from simple model parameter inference to sophisticated machine learning techniques. We conclude with some perspectives on current and future research trends.
△ Less
Submitted 13 April, 2024;
originally announced April 2024.
-
Enhancing non-intrusive Reduced Order Models with space-dependent aggregation methods
Authors:
Anna Ivagnes,
Niccolò Tonicello,
Paola Cinnella,
Gianluigi Rozza
Abstract:
In this manuscript, we combine non-intrusive reduced order models (ROMs) with space-dependent aggregation techniques to build a mixed-ROM. The prediction of the mixed formulation is given by a convex linear combination of the predictions of some previously-trained ROMs, where we assign to each model a space-dependent weight. The ROMs taken into account to build the mixed model exploit different re…
▽ More
In this manuscript, we combine non-intrusive reduced order models (ROMs) with space-dependent aggregation techniques to build a mixed-ROM. The prediction of the mixed formulation is given by a convex linear combination of the predictions of some previously-trained ROMs, where we assign to each model a space-dependent weight. The ROMs taken into account to build the mixed model exploit different reduction techniques, such as Proper Orthogonal Decomposition (POD) and AutoEncoders (AE), and/or different approximation techniques, namely a Radial Basis Function Interpolation (RBF), a Gaussian Process Regression (GPR) or a feed-forward Artificial Neural Network (ANN). The contribution of each model is retained with higher weights in the regions where the model performs best, and, vice versa, with smaller weights where the model has a lower accuracy with respect to the other models. Finally, a regression technique, namely a Random Forest, is exploited to evaluate the weights for unseen conditions. The performance of the aggregated model is evaluated on two different test cases: the 2D flow past a NACA 4412 airfoil, with an angle of attack of 5 degrees, having as parameter the Reynolds number varying between 1e5 and 1e6 and a transonic flow over a NACA 0012 airfoil, considering as parameter the angle of attack. In both cases, the mixed-ROM has provided improved accuracy with respect to each individual ROM technique.
△ Less
Submitted 8 March, 2024;
originally announced March 2024.
-
A priori tests of turbulence models for compressible flows
Authors:
Sciacovelli L.,
Cannici A.,
Passiatore D.,
Cinnella P
Abstract:
A priori tests of turbulence models for the compressible Reynolds-Averaged Navier--Stokes (RANS) are performed by using Direct Numerical Simulations (DNS) data of zero-pressure-gradient flat-plate turbulent boundary layers. The DNS database covers a wide range of operating conditions, ranging from supersonic ($M_\infty = 2.25$) up to the high-enthalpy hypersonic regime ($M_\infty= 12.48$). Several…
▽ More
A priori tests of turbulence models for the compressible Reynolds-Averaged Navier--Stokes (RANS) are performed by using Direct Numerical Simulations (DNS) data of zero-pressure-gradient flat-plate turbulent boundary layers. The DNS database covers a wide range of operating conditions, ranging from supersonic ($M_\infty = 2.25$) up to the high-enthalpy hypersonic regime ($M_\infty= 12.48$). Several RANS closures and compressibility corrections in the literature are assessed against the exact terms from the DNS. Particular attention is paid to closure models for the turbulent heat fluxes and the dilatational dissipation, as well as to the analysis of turbulence/chemistry and turbulence/vibrational relaxation interactions for the high-enthalpy simulations.
△ Less
Submitted 15 October, 2023;
originally announced October 2023.
-
Space-dependent Aggregation of Stochastic Data-driven Turbulence Models
Authors:
Soufiane Cherroud,
Xavier Merle,
Paola Cinnella,
Xavier Gloerfelt
Abstract:
A stochastic Machine-Learning approach is developed for data-driven Reynolds-Averaged Navier-Stokes (RANS) predictions of turbulent flows, with quantified model uncertainty. This is done by combining a Bayesian symbolic identification methodology for learning stochastic RANS model corrections for selected classes of flows (expert models), and a Mixture-of-Experts methodology that aggregates their…
▽ More
A stochastic Machine-Learning approach is developed for data-driven Reynolds-Averaged Navier-Stokes (RANS) predictions of turbulent flows, with quantified model uncertainty. This is done by combining a Bayesian symbolic identification methodology for learning stochastic RANS model corrections for selected classes of flows (expert models), and a Mixture-of-Experts methodology that aggregates their predictions. The expert models are learned using the recently proposed SBL-SpaRTA algorithm, which generates sparse analytical expressions of the corrective terms with model parameters described by probability distributions. They outperform the baseline RANS model for flows similar to those used for training, but their generalization to different flows is not warranted. With the aim of quantifying the predictive uncertainty associated with the data-driven models while improving predictive accuracy and generalization capabilities, a space-dependent model aggregation technique (XMA) is then adopted. A gating function, which assigns each model a performance score (weight) based on a vector of local flow features, is trained alongside the expert models. The weights can be interpreted as the probability that a candidate model will outperform its competitors given the flow behavior at a given location. Predictions of unseen flows are formulated as a locally weighted average of the stochastic solutions of the expert models. A prediction uncertainty estimate is obtained by propagating the models' posterior parameter distributions and by evaluating the inter-model prediction variance. The expectancy of the XMA prediction is significantly more accurate than the baseline deterministic solution and the individual solutions of the experts for well-documented benchmark flows not included in the training set, while providing consistent estimates of the predictive variance.
△ Less
Submitted 3 May, 2024; v1 submitted 29 June, 2023;
originally announced June 2023.
-
Space-dependent turbulence model aggregation using machine learning
Authors:
Maximilien de Zordo-Banliat,
Grégory Dergham,
Xavier Merle,
Paola Cinnella
Abstract:
In this article, we propose a data-driven methodology for combining the solutions of a set of competing turbulence models. The individual model predictions are linearly combined for providing an ensemble solution accompanied by estimates of predictive uncertainty due to the turbulence model choice. First, for a set of training flow configurations we assign to component models high weights in the r…
▽ More
In this article, we propose a data-driven methodology for combining the solutions of a set of competing turbulence models. The individual model predictions are linearly combined for providing an ensemble solution accompanied by estimates of predictive uncertainty due to the turbulence model choice. First, for a set of training flow configurations we assign to component models high weights in the regions where they best perform, and vice versa, by introducing a measure of distance between high-fidelity data and individual model predictions. The model weights are then mapped into a space of features, representative of local flow physics, and regressed by a Random Forests (RF) algorithm. The RF regressor is finally employed to infer spatial distributions of the model weights for unseen configurations. Predictions of new cases are constructed as a convex linear combination of the underlying models solutions, while the between model variance provides information about regions of high model uncertainty. The method is demonstrated for a class of flows through the compressor cascade NACA65 V103 at Re~3e5. The results show that the aggregated solution outperforms the accuracy of individual models for the quantity used to inform the RF regressor, and performs well for other quantities well-correlated to the preceding one. The estimated uncertainty intervals are generally consistent with the target high-fidelity data. The present approach then represents a viable methodology for a more objective selection and combination of alternative turbulence models in configurations of interest for engineering practice
△ Less
Submitted 21 January, 2023;
originally announced January 2023.
-
AirfRANS: High Fidelity Computational Fluid Dynamics Dataset for Approximating Reynolds-Averaged Navier-Stokes Solutions
Authors:
Florent Bonnet,
Ahmed Jocelyn Mazari,
Paola Cinnella,
Patrick Gallinari
Abstract:
Surrogate models are necessary to optimize meaningful quantities in physical dynamics as their recursive numerical resolutions are often prohibitively expensive. It is mainly the case for fluid dynamics and the resolution of Navier-Stokes equations. However, despite the fast-growing field of data-driven models for physical systems, reference datasets representing real-world phenomena are lacking.…
▽ More
Surrogate models are necessary to optimize meaningful quantities in physical dynamics as their recursive numerical resolutions are often prohibitively expensive. It is mainly the case for fluid dynamics and the resolution of Navier-Stokes equations. However, despite the fast-growing field of data-driven models for physical systems, reference datasets representing real-world phenomena are lacking. In this work, we develop AirfRANS, a dataset for studying the two-dimensional incompressible steady-state Reynolds-Averaged Navier-Stokes equations over airfoils at a subsonic regime and for different angles of attacks. We also introduce metrics on the stress forces at the surface of geometries and visualization of boundary layers to assess the capabilities of models to accurately predict the meaningful information of the problem. Finally, we propose deep learning baselines on four machine learning tasks to study AirfRANS under different constraints for generalization considerations: big and scarce data regime, Reynolds number, and angle of attack extrapolation.
△ Less
Submitted 1 June, 2023; v1 submitted 14 December, 2022;
originally announced December 2022.
-
Shock impingement on a transitional hypersonic high-enthalpy boundary layer
Authors:
Donatella Passiatore,
Luca Sciacovelli,
Paola Cinnella,
Giuseppe Pascazio
Abstract:
The dynamics of a shock wave impinging on a transitional high-enthalpy boundary layer out of thermochemical equilibrium is investigated for the first time by means of a direct numerical simulation. The freestream Mach number is equal to 9 and the oblique shock impinges with a cooled flat-plate boundary layer with an angle of 10°, generating a reversal flow region. In conjunction with freestream di…
▽ More
The dynamics of a shock wave impinging on a transitional high-enthalpy boundary layer out of thermochemical equilibrium is investigated for the first time by means of a direct numerical simulation. The freestream Mach number is equal to 9 and the oblique shock impinges with a cooled flat-plate boundary layer with an angle of 10°, generating a reversal flow region. In conjunction with freestream disturbances, the shock impingement triggers a transition to a fully turbulent regime shortly downstream of the interaction region. Accordingly, wall properties emphasize the presence of a laminar region, a recirculation bubble, a transitional zone and fully turbulent region. The breakdown to turbulence is characterized by an anomalous increase of skin friction and wall heat flux, due to the particular shock pattern. At the considered thermodynamic conditions the flow is found to be in a state of thermal non-equilibrium throughout, with non-equilibrium effects enhanced and sustained by the shock-induced laminar/turbulent transition, while chemical activity is almost negligible due to wall cooling. In the interaction region, relaxation towards thermal equilibrium is delayed and the fluctuating values of the rototranslational and the vibrational temperatures strongly differ, despite the strong wall-cooling. The fully turbulent portion exhibits evolutions of streamwise velocity, Reynolds stresses and turbulent Mach number in good accordance with previous results for highly-compressible cooled-wall boundary layers in thermal nonequilibrium, with turbulent motions sustaining thermal nonequilibrium. Nevertheless, the vibrational energy is found to contribute little to the total wall heat flux.
△ Less
Submitted 13 December, 2022;
originally announced December 2022.
-
CFD-driven Symbolic Identification of Algebraic Reynolds-Stress Models
Authors:
I. Ben Hassan SaÏdi,
M. Schmelzer,
P. Cinnella,
F. Grasso
Abstract:
A CFD-driven deterministic symbolic identification algorithm for learning explicit algebraic Reynolds-stress models (EARSM) from high-fidelity data is developed building on
the frozen-training SpaRTA algorithm of [1].
Corrections for the Reynolds stress tensor and the production of transported turbulent quantities of a baseline linear eddy viscosity model (LEVM) are expressed as functions of t…
▽ More
A CFD-driven deterministic symbolic identification algorithm for learning explicit algebraic Reynolds-stress models (EARSM) from high-fidelity data is developed building on
the frozen-training SpaRTA algorithm of [1].
Corrections for the Reynolds stress tensor and the production of transported turbulent quantities of a baseline linear eddy viscosity model (LEVM) are expressed as functions of tensor polynomials selected from a library of candidate functions. The CFD-driven training consists in solving a blackbox optimization problem in which the fitness of candidate EARSM models is evaluated by running RANS simulations. Unlike the frozen-training approach, the proposed methodology is not restricted to data sets for which full fields of high-fidelity data are available. However, the solution of a high-dimensional expensive blackbox function optimization problem is required. Several steps are then undertaken to reduce the associated computational burden. First, a sensitivity analysis is used to identify the most influential terms and to reduce the dimensionality of the search space. Afterwards, the Constrained Optimization using Response Surface (CORS) algorithm, which approximates the black-box cost function using a response surface constructed from a limited number of CFD solves, is used to find the optimal model parameters. Model discovery and cross-validation is performed for three configurations of 2D turbulent separated flows in channels of variable section using different sets of training data to show the flexibility of the method. The discovered models are then applied to the prediction of an unseen 2D separated flow with higher Reynolds number and different geometry. The predictions for the new case are shown to be not only more accurate than the baseline LEVM, but also of a multi-purpose EARSM model derived from purely physical arguments.
△ Less
Submitted 2 December, 2021; v1 submitted 19 April, 2021;
originally announced April 2021.
-
Finite-rate chemistry effects in turbulent hypersonic boundary layers: a direct numerical simulation study
Authors:
Donatella Passiatore,
Luca Sciacovelli,
Paola Cinnella,
Giuseppe Pascazio
Abstract:
The influence of high-enthalpy effects on hypersonic turbulent boundary layers is investigated by means of direct numerical simulations (DNS). A quasi-adiabatic flat-plate air flow at free-stream Mach number equal to 10 is simulated up to fully-developed turbulent conditions using a five-species, chemically-reacting model. A companion DNS based on a frozen-chemistry assumption is also carried out,…
▽ More
The influence of high-enthalpy effects on hypersonic turbulent boundary layers is investigated by means of direct numerical simulations (DNS). A quasi-adiabatic flat-plate air flow at free-stream Mach number equal to 10 is simulated up to fully-developed turbulent conditions using a five-species, chemically-reacting model. A companion DNS based on a frozen-chemistry assumption is also carried out, in order to isolate the effect of finite-rate chemical reactions and assess their influence on turbulent quantities. In order to reduce uncertainties associated with turbulence generation at the inlet of the computational domain, both simulations are initiated in the laminar flow region and the flow is let to evolve up to the fully turbulent regime. Modal forcing by means of localized suction and blowing is used to trigger laminar-to-turbulent transition. The high temperatures reached in the near wall region including the viscous and buffer sublayers activate significant dissociation of both oxygen and nitrogen. This modifies in turn the thermodynamic and transport properties of the reacting mixture, affecting the first-order statistics of thermodynamic quantities. Due to the endothermic nature of the chemical reactions in the forward direction, temperature and density fluctuations in the reacting layer are smaller than in the frozen-chemistry flow. However, the first- and second-order statistics of the velocity field are found to be little affected by the chemical reactions under a scaling that accounts for the modified fluid properties. We also observed that the Strong Reynolds Analogy (SRA) remains well respected despite the severe hypersonic conditions and that the computed skin friction coefficient distributions match well the results of the Renard-Deck decomposition extended to compressible flows.
△ Less
Submitted 30 March, 2021;
originally announced March 2021.
-
Assessment of a high-order shock-capturing central-difference scheme for hypersonic turbulent flow simulations
Authors:
Luca Sciacovelli,
Donatella Passiatore,
Paola Cinnella,
Giuseppe Pascazio
Abstract:
High-speed turbulent flows are encountered in most space-related applications (including exploration, tourism and defense fields) and represent a subject of growing interest in the last decades. A major challenge in performing high-fidelity simulations of such flows resides in the stringent requirements for the numerical schemes to be used. These must be robust enough to handle strong, unsteady di…
▽ More
High-speed turbulent flows are encountered in most space-related applications (including exploration, tourism and defense fields) and represent a subject of growing interest in the last decades. A major challenge in performing high-fidelity simulations of such flows resides in the stringent requirements for the numerical schemes to be used. These must be robust enough to handle strong, unsteady discontinuities, while ensuring low amounts of intrinsic dissipation in smooth flow regions. Furthermore, the wide range of temporal and spatial active scales leads to concurrent needs for numerical stabilization and accurate representation of the smallest resolved flow scales in cases of under-resolved configurations. In this paper, we present a finite-difference high-order shock-capturing technique based on Jameson's artificial diffusivity methodology. The resulting scheme is ninth-order-accurate far from discontinuities and relies on the addition of artificial dissipation close to large gradients. The shock detector is slightly revised to enhance its selectivity and avoid spurious activations of the shock-capturing term. A suite of test cases ranging from 1D to 3D configurations (namely, shock tubes, Shu-Osher problem, isentropic vortex advection, under-expanded jet, compressible Taylor-Green Vortex, supersonic and hypersonic turbulent boundary layers) is analysed in order to test the capability of the proposed numerical strategy to handle a large variety of problems, ranging from calorically-perfect air to multi-species reactive flows. Results obtained on under-resolved grids are also considered to test the applicability of the proposed strategy in the context of implicit Large-Eddy Simulations.
△ Less
Submitted 30 March, 2021;
originally announced March 2021.
-
Discovery of Algebraic Reynolds-Stress Models Using Sparse Symbolic Regression
Authors:
Martin Schmelzer,
Richard P. Dwight,
Paola Cinnella
Abstract:
** This article is published (open-access). ** A novel deterministic symbolic regression method SpaRTA (Sparse Regression of Turbulent Stress Anisotropy) is introduced to infer algebraic stress models for the closure of RANS equations directly from high-fidelity LES or DNS data. The models are written as tensor polynomials and are built from a library of candidate functions. The machine-learning m…
▽ More
** This article is published (open-access). ** A novel deterministic symbolic regression method SpaRTA (Sparse Regression of Turbulent Stress Anisotropy) is introduced to infer algebraic stress models for the closure of RANS equations directly from high-fidelity LES or DNS data. The models are written as tensor polynomials and are built from a library of candidate functions. The machine-learning method is based on elastic net regularisation which promotes sparsity of the inferred models. By being data-driven the method relaxes assumptions commonly made in the process of model development. Model-discovery and cross-validation is performed for three cases of separating flows, i.e. periodic hills ($Re$=10595), converging-diverging channel ($Re$=12600) and curved backward-facing step ($Re$=13700). The predictions of the discovered models are significantly improved over the $k$-$ω$ SST also for a true prediction of the flow over periodic hills at $Re$=37000. This study shows a systematic assessment of SpaRTA for rapid machine-learning of robust corrections for standard RANS turbulence models.
△ Less
Submitted 28 February, 2020; v1 submitted 17 May, 2019;
originally announced May 2019.
-
Quantification of Model Uncertainty in RANS Simulations: A Review
Authors:
Heng Xiao,
Paola Cinnella
Abstract:
In computational fluid dynamics simulations of industrial flows, models based on the Reynolds-averaged Navier--Stokes (RANS) equations are expected to play an important role in decades to come. However, model uncertainties are still a major obstacle for the predictive capability of RANS simulations. This review examines both the parametric and structural uncertainties in turbulence models. We revi…
▽ More
In computational fluid dynamics simulations of industrial flows, models based on the Reynolds-averaged Navier--Stokes (RANS) equations are expected to play an important role in decades to come. However, model uncertainties are still a major obstacle for the predictive capability of RANS simulations. This review examines both the parametric and structural uncertainties in turbulence models. We review recent literature on data-free (uncertainty propagation) and data-driven (statistical inference) approaches for quantifying and reducing model uncertainties in RANS simulations. Moreover, the fundamentals of uncertainty propagation and Bayesian inference are introduced in the context of RANS model uncertainty quantification. Finally, the literature on uncertainties in scale-resolving simulations is briefly reviewed with particular emphasis on large eddy simulations.
△ Less
Submitted 27 September, 2018; v1 submitted 27 June, 2018;
originally announced June 2018.
-
A return to eddy viscosity model for epistemic UQ in RANS closures
Authors:
W. N. Edeling,
G. Iaccarino,
P. Cinnella
Abstract:
For the purpose of Uncertainty Quantification (UQ) of Reynolds-Averaged Navier-Stokes closures, we introduce a framework in which perturbations in the eigenvalues of the anisotropy tensor are made in order to bound a Quantity-of-Interest based on limiting states of turbulence. To make the perturbations representative of local flow features, we introduce two additional transport equations for linea…
▽ More
For the purpose of Uncertainty Quantification (UQ) of Reynolds-Averaged Navier-Stokes closures, we introduce a framework in which perturbations in the eigenvalues of the anisotropy tensor are made in order to bound a Quantity-of-Interest based on limiting states of turbulence. To make the perturbations representative of local flow features, we introduce two additional transport equations for linear combinations of these aforementioned eigenvalues. The location, magnitude and direction of the eigenvalue perturbations are now governed by the model transport equations. The general behavior of our discrepancy model is determined by two coefficients, resulting in a low-dimensional UQ problem. We will furthermore show that the behavior of the model is intuitive and rooted in the physical interpretation of misalignment between the mean strain and Reynolds stresses.
△ Less
Submitted 15 May, 2017;
originally announced May 2017.