Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Supersaturation Variability from Scalar Mixing: Evaluation of a New Subgrid-Scale Model Using Direct Numerical Simulations of Turbulent Rayleigh–Bénard Convection

We report that supersaturation fluctuations in the atmosphere are critical for cloud processes. A nonlinear dependence on two scalars - water vapor and temperature - leads to different behavior than single scalars in turbulent convection. For modeling such multiscalar processes at subgrid scales (SGS) in large-eddy simulations (LES) or convection-permitting models, a new SGS scheme is implemented in CM1 that solves equations for SGS water vapor and temperature fluctuations and their covariance. The SGS model is evaluated using benchmark direct-numerical simulations (DNS) of turbulent Rayleigh–Bénard convection with water vapor as in the Michigan Tech Pi Cloud Chamber. This idealized setup allows thorough evaluation of the SGS model without complications from other atmospheric processes. DNS results compare favorably with measurements from the chamber. Results from LES using the new SGS model compare well with DNS, including profiles of water vapor and temperature variances, their covariance, and supersaturation variance. SGS supersaturation fluctuations scale appropriately with changes to the LES grid spacing, with the magnitude of SGS fluctuations decreasing relative to those at the resolved scale as the grid spacing is decreased. Sensitivities of covariance and supersaturation statistics to changes in water vapor flux relative to thermal flux are also investigated by modifying the sidewall conditions. Relative changes in water vapor flux substantially decrease the covariance and increase supersaturation fluctuations even away from boundaries.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

High-fidelity Experiments of Turbulent Buoyant Jets from Simultaneous Particle Image Velocimetry and Laser Induced Fluorescence

Accurate models of turbulent buoyant flows are essential for the design of nuclear reactors thermal hydraulics and passive safety systems. However, available models fail to fully capture the physics of turbulent mixing when buoyancy becomes predominant with respect to momentum. Therefore, high-fidelity experiments of well-controlled fundamental flows are needed to develop and validate more accurate models. We performed experiments of positive and negative turbulent buoyant jets, both in uniform and stratified environments, with the aim of understanding the thermal hydraulics of turbulent mixing with variable density and providing high-fidelity data for the development and validation of turbulence models. Non-intrusive, simultaneous Particle Image Velocimetry and Laser Induced Fluorescence measurements were carried out to acquire instantaneous velocity and concentration fields on a vertical section parallel to the axis of a jet in the self-similar region. The Refractive Index Matching method was applied to measure high-resolution buoyant jets with up to 8.6% density difference. These data are free of the typical errors that characterize optical measurements of buoyancy driven flows (e.g., natural and mixed convection) where the refractive index of the fluid is inhomogeneous throughout the measurement domain. Turbulent statistics and entrainment of buoyant jets in uniform and stratified environments are presented. These data are compared with non-buoyant jets in uniform environment, as a reference to investigate the effects of buoyancy and stratification on turbulent mixing. The results will be used for the assessment of current turbulence models and as basis for the development of a new one that captures turbulent mixing.

42 ENGINEERING↗
DOE OSTI2022

How Generalizable is a Machine-Learning Approach for Modeling Hub-Height Turbulence Intensity?

Hub-height turbulence intensity is essential for a variety of wind energy applications. However, simulating it is a challenging task. Simple analytical models have been proposed in the literature, but they all come with significant limitations. Even state-of-the-art numerical weather prediction models, such as the Weather Research and Forecasting model, currently struggle to predict hub-height turbulence intensity. Here, we propose a machine-learning-based approach to predict hub-height turbulence intensity from other hub-height and ground-level atmospheric measurements, using observations from the Perdigao field campaign and the Southern Great Plains atmospheric observatory. We consider a random forest regression model, which we validate first at the site used for training and then under a more robust round-robin approach, and compare its performance to a multivariate linear regression. The random forest successfully outperforms the linear regression in modeling hub-height turbulence intensity, with a normalized root-mean-square error as low as 0.014 when using 30-minute average data. In order to achieve such low root-mean-square error values, the knowledge of hub-height turbulence kinetic energy (which can instead be modeled in the Weather Research and Forecasting model) is needed. Interestingly, we find that the performance of the random forest generalizes well when considering a round-robin validation (i.e., when the algorithm is trained at one site such as Perdigao or Southern Great Plains) and then applied to model hub-height turbulence intensity at the other location.

17 WIND ENERGY↗
DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Physics of turbulence spreading and explicit nonlocality

In this paper, we systematically derive a model for turbulence spreading from the basic kinetic equation. The model contains explicit nonlocal nonlinear diffusion and nonlocal growth. When the nonlocality scale parameter δ b (banana width) vanishes, this model reduces to the usual turbulence spreading model. We elucidate the mechanisms of nonlinear saturation and nonlocal growth. Results show that nonlocal effects, especially the nonlocal growth, thicken the turbulence spreading front and increase the speed of front propagation. More turbulence intensity penetrates the stable region when δ b increases. The penetration depth Δ p is proportional to $ \delta_b/L_T $, therefore the fraction of turbulence in the unstable region scales as 1 – δ b* . Lastly, the transport scales the same way.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Incorporating long-range dependence and fractal features in turbulence spectra

We introduce an advanced turbulence spectrum model developed from mathematical foundations from a covariance function class and empirically validated using extensive field data. This model captures the complex dynamics of long-range dependence, and fractal characteristics prevalent in riverine and atmospheric boundary layer (ABL) flows that are ignored by classical spectrum models, such as IEC (International Electrotechnical Commission) von Kármán and Kaimal model. The model delineates scaling behaviors across distinct frequency bands and offers substantial flexibility through five well-defined parameters each characterizing a distinct physical aspect of the velocity time series. A detailed procedure for obtaining each parameter from time series data is outlined. The comprehensive validations with field data from tidal currents and ABL flows substantiate the model’s fidelity in accurately replicating observed phenomena. This validation establishes the reliability of the proposed model and, when incorporated into stochastic full-field simulators such as TurbSim, demonstrates its potential to advance the predictive modeling and analysis of turbulent flows in environmental science and engineering contexts.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗
DOE OSTI2020

Numerical Investigation of Fuel Property Effects on Mixed-Mode Combustion in a Spark-Ignition Engine

In this research, lean mixed-mode combustion is numerically investigated using computational fluid dynamics (CFD) in a spark-ignition engine. A new E30 fuel surrogate is developed using a neural network model with matched octane numbers. A skeletal mechanism is also developed by automated mechanism reduction and by incorporating a NO x submechanism. A hybrid approach that couples the G-equation model and the well-stirred reactor model is employed for turbulent combustion modeling. The developed CFD model is shown to well predict pressure and apparent heat release rate (AHRR) traces compared with experiment. Two types of combustion cycles (deflagration-only and mixed-mode cycles) are observed. The mixed-mode cycles feature early flame propagation and subsequent end-gas auto-ignition, leading to two distinctive AHRR peaks. The validated CFD model is then employed to investigate the effects of NO x chemistry. The NO x chemistry is found to promote auto-ignition through the residual gas, while the deflagration phase remains largely unaffected. Sensitivity analysis is finally performed to understand effects of fuel properties, including heat of vaporization (HoV) and laminar flame speed (S L ). An increased HoV tends to suppress auto-ignition through charge cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased S L is found to significantly promote both flame propagation and end-gas auto-ignition. The promoting effect of S L on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end-gas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Advanced surrogate model for electron-scale turbulence in tokamak pedestals

We derive an advanced surrogate model for predicting turbulent transport at the edge of tokamaks driven by electron temperature gradient (ETG) modes. Our derivation is based on a recently developed sensitivity-driven sparse grid interpolation approach for uncertainty quantification and sensitivity analysis at scale, which informs the set of parameters that define the surrogate model as a scaling law. Our model reveals that ETG-driven electron heat flux is influenced by the safety factor q, electron beta β e and normalized electron Debye length λ D , in addition to well-established parameters such as the electron temperature and density gradients. To assess the trustworthiness of our model's predictions beyond training, we compute prediction intervals using bootstrapping. The surrogate model's predictive power is tested across a wide range of parameter values, including within-distribution testing parameters (to verify our model) as well as out-of-bounds and out-of-distribution testing (to validate the proposed model). Overall, validation efforts show that our model competes well with, or can even outperform, existing scaling laws in predicting ETG-driven transport.

fusion plasma↗
DOE OSTI2025

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗
DOE OSTI2023

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Learning turbulent flows with generative models for super resolution and sparse flow reconstruction

Neural operators are promising surrogates for dynamical systems but when trained with standard L 2 losses they tend to oversmooth fine-scale turbulent structures. Here, we show that combining operator learning with generative modeling overcomes this limitation. We consider three practical turbulent-flow challenges where conventional neural operators fail: spatio-temporal super-resolution, forecasting, and sparse flow reconstruction. For Schlieren jet super-resolution, an adversarially trained neural operator (adv-NO) reduces the energy-spectrum error by 15 × while preserving sharp gradients at neural operator-like inference cost. For 3D homogeneous isotropic turbulence, adv-NO trained on only 160 timesteps from a single trajectory forecasts accurately for five eddy-turnover times and offers 114 × wall-clock speed-up at inference than the baseline diffusion-based forecasters, enabling near-real-time rollouts. For reconstructing cylinder wake flows from highly sparse Particle Tracking Velocimetry-like inputs, a conditional generative model infers full 3D velocity and pressure fields with correct phase alignment and statistics. These advances enable accurate reconstruction and forecasting at low compute cost, bringing near-real-time analysis and control within reach in experimental and computational fluid mechanics.

Fluid dynamics↗
DOE OSTI2021

Self-consistent, high-order spatial profiles in a model for two-fluid turbulent mixing

A Reynolds-averaged Navier-Stokes model is presented with the property that it admits self-consistent, high-order spatial profiles in simulations of two-fluid turbulent mixing layers. Whereas previous models have been limited by the assumption of a linear mixing profile, the present paper relaxes this assumption and, as a result, is shown to achieve much better agreement with experimental profiles. Similarity analysis is presented to derive constraints on model coefficients to enforce desired self-similar growth rates that are fully consistent with the high-order spatial profiles. Through this similarity analysis, it is shown that care must be taken in model construction, as it is possible to construct certain terms in such a way as to leave growth rates unconstrained. This model, termed the k-Φ-L-a-V model, is then applied in simulations of Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz mixing layers. Overall, these simulations confirm that the expected growth parameters are recovered and high-order spatial profiles are maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Difference in load predictions obtained with effective turbulence vs. a dynamic wake meandering modeling approach

Abstract. According to the international standard for wind turbine design, the effects of wind turbine wakes on structural loads can be considered in two ways: (1) by augmenting the ambient turbulence levels with the effective turbulence model (EFF) and then calculating the resulting loads and (2) by performing dynamic wake meandering (DWM) simulations, which compute wake effects and loads for all turbines on a farm at once. There is no definitive answer in scientific literature as to the consequences of choosing one model over the other, but the two approaches are unarguably very different. The work presented here expounds on these differences and investigates to what extent they affect the simulated structural loads. We consider an idealized 4×4 rectangular array of National Renewable Energy Laboratory 5 MW wind turbines with a spacing of 5 by 8 rotor diameters and three wind speed scenarios at high ambient turbulence. Load simulations are performed in OpenFAST with EFF and in FAST.Farm with the DWM implementation. We compare ambient turbulence, wind farm turbulence, and loads between both approaches. When omnidirectional results are compared, EFF wind farm turbulence intensity is consistently higher by 0.2 % (above-rated wind speed) to 2.7 % (below-rated wind speed). However, for certain wind directions, the EFF turbulence can be lower than FAST.Farm by almost 9 %. Wind speeds within the farm were found to differ by up to 3 m s−1 due to the lack of wake deficits in the EFF approach, leading to longer tails toward low values in the FAST.Farm mean load distributions. Consistent with the turbulence results, the median EFF load standard deviations are also consistently higher, by a maximum of 20 % and 17 % for blade-root out-of-plane and tower-base fore-aft moments, respectively.

17 WIND ENERGY↗
DOE OSTI2023

A data-driven method for modelling dissipation rates in stratified turbulence

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$ . Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$ , where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$ , and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$ ), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

42 ENGINEERING↗
DOE OSTI2022

Accurate and Efficient Parametric Model-Order Reduction for Turbulent Thermal Transport

This project produced new algorithms and software for low-cost exploration of turbulent thermal-fluids behavior under parametric variation by using reduced-order models (ROMs). The ROMs numerically solve the governing equations for fluid motion by using a small set of basis functions (typically, N=20-200 modes) to represent the solution. The base modes are computed as optimal combinations of solutions from expensive high-fidelity "anchor-point" solutions involving millions of unknowns, which are typically generated by solving the full Navier-Stokes equations on a supercomputer. The ROM solution is itself a combination of the base modes, where the basis coefficients are determined by evolving an NxN system of nonlinear equations. The overarching idea is to use the inexpensive ROM to predict solutions under conditions where the parameters differ from the anchor-point conditions. Several ingredients are required to make ROMs useful for thermal hydraulics analysis. These include: a stable and accurate ROM that is capable of reproducing the large-scale dynamics of turbulent flow, error indicators than can guide the choice of anchor points, and low-cost mechanisms for evaluating nonlinear terms in the reduced equations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Phase Contrast Imaging Diagnostic on DIII-D (1994-2007); Phase Contrast Imaging Diagnostic on DIII-D and Alcator C-Mod (2007-2010); Experimental Study of Waves, Turbulence and Transport, and Model Validation with Phase Contrast Imaging on DIII-D and Alcator C-Mod (2010-2012); Phase Contrast Imaging Diagnostic Upgrade on DIII-D for Measuring High-k Turbulence and Alfven modes (2021-2016); Phase Contrast Imaging Diagnostic Upgrade on DIII-D for Full Multi-Scale Measurements (2016-2020) (Final Report)

This report describes the first twenty-five years of development and use of the Phase Contrast Imaging (PCI) diagnostic for measuring plasma fluctuations on the DIII–D tokamak. One of the key challenges in moving tokamak technology forward to the stage of a nuclear reactor remains understanding how energy and particles are transported across the magnetic field by the action of turbulence and waves, and PCI provides a highly flexible and capable measurement of these fluctuations. The PCI diagnostic uses a laser beam to probe the plasma and optical filtering to produce an amplitude variation on a linear detector array which provides an image of the density waves in the plasma. The PCI covers an exceptional range in frequency and wavelength of the measured fluctuations, allowing observation of plasma-sized magnetohydrodynamic (MHD) modes and ion- and electron-scale turbulence. Modulation of the beam amplitude allows measurement of radio frequency (RF) waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗