Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TURBULENCE EFFECT”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗
DOE OSTI2025

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2024

The Effects of Turbulence on Heliosheath Ions and Implications for Energetic Neutral Atoms

The distribution of ions in the heliosheath—the region between the heliospheric termination shock and the heliopause—is important for understanding remote observations of energetic neutral atoms (ENAs). The ion distributions were estimated previously based on hybrid simulations of the heating and evolution of solar wind and interstellar pickup ions across the solar wind termination shock, but these estimates only provide the distributions near the shock. In this work, we use self-consistent hybrid kinetic simulations to investigate the effects of turbulence on ion distributions in the heliosheath. The simulations are compared against Voyager observations, constraining the feasible amplitude and compressibility of turbulence. We find that the heating due to turbulent dissipation can lead to a significant increase in the temperature of thermal solar wind ions. Both turbulent velocity fluctuations and the heating of solar wind ions increase the charge-exchange source for ENAs at low energies (around 100 eV), where current ENA models underpredict observations by more than an order of magnitude. However, the effects of turbulence are likely not strong enough to fully explain these discrepancies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Magnetic geometry effects on turbulent density fluctuations in Wendelstein 7-X

Line-integrated density fluctuation measurements from the core of the optimised stellarator Wendelstein 7-X (W7-X) are studied in the context of magnetic field geometry changes. The amplitude as well as the spectral distribution of fluctuations is affected by varying the magnetic field configuration. Weaker turbulent fluctuations are observed in configurations with higher rotational transform and lower magnetic mirror depth. Simulations with the gyrokinetic codes stella and GENE-3D are used to investigate the role of magnetic geometry for the experimental observations. The stabilising effect of the high rational transform configuration is reproduced by simulations and can be associated with stronger shaping of the magnetic field. This leads to increased finite Larmor radius and local shear effects, which limit the toroidal extent of ion temperature gradient (ITG) modes. For W7-X configurations with higher magnetic mirror depth, gyrokinetic theory and simulations predict weaker fluctuations, since they are closer to having the maximum-J property, which has a stabilising effect on trapped-electron mode and ITG. The opposite trend is observed experimentally, which can therefore not be explained by geometry effects on turbulence alone.

Wendelstein-7X↗
DOE OSTI2024

Temperature effect on turbulent burning velocity of lean premixed hydrogen/air flames

Hydrogen has drawn great attention in recent years as a carbon-free fuel. The turbulent burning velocity (S T ) is an important parameter for the design and modeling of hydrogen-fueled engines given the high propagation speed of hydrogen flames. It has been well documented that S T of hydrogen flames can be dramatically increased by thermo-diffusive effects which are sensitive to thermodynamic conditions. Previous studies have mainly focused on the pressure effect on S T of lean hydrogen flames, while the temperature effect has been largely ignored. Here, in the present study, the turbulent burning velocity for a lean hydrogen/air mixture over a wide range of temperatures (300–641 K) and pressures (1–15 atm) is investigated through direct numerical simulations of statistically planar turbulent premixed flames. Results show that the variation of normalized turbulent burning velocity (S T /S L , where S L is the laminar flame speed) with temperature and pressure is mainly controlled by the variation of the stretching factor I 0 . While S T /S L is only marginally dependent on temperature at the atmospheric pressure, it exhibits a decreasing trend with temperature at an elevated pressure (10 atm). This is associated with different temperature dependencies of flame surface area enlargement at the two different pressures, despite the monotonically decreasing trends of I 0 with temperature at both pressures. In addition, under engine-relevant conditions where the temperature and pressure increase simultaneously, the promotion effect of pressure is found to be largely canceled out by the suppression effect of temperature, leading to only a slight increase in I 0 and S T /S L . The observed trends are further explained through detailed flame dynamic analysis. Furthermore, I 0 at different temperatures and pressures is found to correlate very well with the enhancement of fuel consumption rate in the critically strained laminar flames. The present study elucidates the strong impact of temperature on S T of lean premixed hydrogen/air flames at elevated pressures and provides new insights into the modeling of S T , especially for engine-relevant conditions.

Wang, Yiqing [Argonne National Laboratory (ANL), A↗
DOE OSTI2025

Effect of Surface Roughness on Dynamic Stall in Pitching Motion

Dynamic stall plays a critical role in determining the performance and stability of a wide range of fluid-dynamics systems in various engineering applications. This unsteady aerodynamic phenomenon is particularly significant for maneuvering aircraft wings, jet aircraft subjected to gust encounters, helicopter rotor blades, and wind turbine blades [1–3]. The prediction of the dynamic stall vortex (DSV) is challenging due to factors such as unsteady aerodynamics, three-dimensional (3-D) effects, turbulence and flow separation, incoming gust, and surface impact effects [4–6]. Hence, advanced computational techniques and modeling approaches in computational fluid dynamics (CFD) would be required to enhance the accuracy and reliability of DSV predictions in dynamic stall scenarios. In the past, Batther and Lee [7] employed delayed detached eddy simulations (DDES) to understand the flow physics associated with the onset of dynamic stall. Their approach demonstrated that DDES achieves results comparable to those obtained from large-eddy simulations at a reduced computational cost. In another study, Khalifa et al. [8] examined the 3-D aspects of dynamic stall on a NACA 0012 airfoil using DES solvers. In conclusion, the findings underlined the superiority of 3-D simulations over two-dimensional approaches, particularly in predicting the lift coefficient values and capturing dynamic stall stages more precisely.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A generative artificial intelligence framework for long-time plasma turbulence simulations

Generative deep learning techniques are employed in a novel framework for the construction of surrogate models capturing the spatiotemporal dynamics of 2D plasma turbulence. The proposed Generative Artificial Intelligence Turbulence (GAIT) framework enables the acceleration of turbulence simulations for long-time transport studies. GAIT leverages a convolutional variational auto-encoder and a recurrent neural network to generate new turbulence data from existing simulations, extending the time horizon of transport studies with minimal computational cost. The application of the GAIT framework to plasma turbulence using the Hasegawa–Wakatani (HW) model is presented, evaluating its performance via various analyses. Very good agreement is found between the GAIT and the HW models in the spatiotemporal Fourier and Proper Orthogonal Decomposition spectra, the flow topology characterized by the Okubo–Weiss parameter, and the time autocorrelation function of turbulent fluctuations. Excellent agreement has also been obtained in the probability distribution function of particle displacements and the effective turbulent diffusivity. In-depth analyses of the latent space of turbulent states, choice of hyperparameters and alternative deep learning models for the time prediction are presented. Our results highlight the potential of Artificial Intelligence-based surrogate models to overcome the computational challenges in turbulence simulation, which can be extended to other situations such as geophysical fluid dynamics.

Artificial intelligence↗
DOE OSTI2025

Sensitivity analysis of a parameterization of coagulation in an aqueous transport model

A three-dimensional aqueous transport model, ALGE, was developed at the Savannah River National Laboratory (SRNL) as a tool for emergency response. Recently, coagulation and break-up of sediment and particulate (contaminant adhered to sediment) matter has been added as a new model capability, representing both processes as functions of salinity and turbulence. A sensitivity analysis was performed by comparing time series of sediment and particulate concentrations at various locations when the model’s parameters (particle size and density) are altered. The results show significant differences between experiments in both freshwater and saltwater systems parameters vary (P-values < 0.05). From this sensitivity study, we have concluded the importance of incorporating both salinity and turbulent effects, as they have a large role in the physical process of coagulation and breakup of sediment. By accounting for these physical processes, we have made the ALGE model a more robust tool for emergency response.

Aqueous fate and transport↗
DOE OSTI2025

Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows

We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.

compressible flows↗
DOE OSTI2026

Effects of density turbulence on helicon wave propagation in the core plasmas

Radio frequency wave propagation can be significantly affected by density irregularities, such as filaments in the scrape-off layer or instabilities in the core plasma. In this study, we examine the impact of edge turbulence on helicon wave propagation using the Petra-M simulation tool. To analyze the effect of edge turbulence, we utilize a realistic background plasma derived from XGC simulations, which includes spatial density fluctuations at the edge. This focus is particularly relevant for the DIII-D configuration, characterized by edge density turbulence of the core plasmas in a wide pedestal QH-mode. We focus on helicon wave propagation in the core plasma since the slow mode cannot propagate into the core plasma when the density is higher than the lower hybrid resonance, where a wide pedestal QH-mode occurs. The simulation results indicate that edge density fluctuations have a substantial impact on helicon wave coupling due to scattering in the core. We specifically demonstrate that the toroidal mode number and the level of density fluctuations are significant factors influencing scattering. While we have minimized the excitation of slow modes from the antenna, we still observe mode-converted slow modes resulting from the incoming helicon waves. The insights derived from these simulations will inform upcoming tokamak experiments regarding helicon antenna coupling in long pulse scenarios.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2025

Role of perturbed parallel magnetic field effects in predicting turbulent transport in NSTX

This study presents analysis of gyrokinetic simulations on the National Spherical Torus Experiment (NSTX) to investigate the effects of electromagnetic fields on plasma turbulence and transport. The simulations, performed with varying levels of fidelity using the gyrokinetic CGYRO code, include electrostatic (ES), single-field electromagnetic (EM1), and two-field electromagnetic (EM2) models. A detailed comparison across the simulation database reveals that electromagnetic effects increase both predicted growth rates and quasilinear fluxes, with EM2 simulations producing stronger turbulence than ES and EM1 cases. Quasilinear modeling using QLGYRO demonstrates that while the perturbed parallel magnetic field (δB ∥ ) does not drastically affect the total flux at experimental gradients, it leads to a shift in the dominant instability, altering mode structures from microtearing to kinetic ballooning modes (KBMs). The proximity of the plasma profiles to the KBM threshold is explored, with the experimental conditions being near the onset of KBM-driven transport. The KBM, with its large growth rates, is identified as a potential driver of electron temperature flattening, as it can rapidly transport heat across flux surfaces. Performing stability analysis shows core-localized unstable a low- mode that could contribute to the flattening at the early times of the discharge. TGYRO predictive modeling, incorporating both TGLF and QLGYRO, indicates that the inclusion of δB ∥ significantly improves the accuracy of temperature profile predictions in NSTX high-beta plasmas, although challenges remain in modeling the sharp flux discontinuities caused by KBM-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗
DOE OSTI2025

Real-time steerable frequency-stepped Doppler backscattering (DBS) system for local helicon wave electric field measurements on the DIII-D tokamak

A new frequency-stepped Doppler backscattering (DBS) system has been integrated into a real-time steerable electron cyclotron heating launcher system to simultaneously probe local background turbulence (f < 10 MHz) and high-frequency (20–550 MHz) density fluctuations in the DIII-D tokamak. The launcher allows for 2D steering (horizontally and vertically) over wide angular ranges to optimize probe location and wavenumber response. The vertical steering can be optimized during a discharge in real time. The new DBS system employs a programmable frequency synthesizer with adjustable dwell time as a source to launch either O or X-mode polarized millimeter waves. This system can step in real-time over the entire E-band frequency range (60–90 GHz). This combination of capabilities allows for the diagnosis of the complex internal spatial structure of high power (>200 kW) helicon waves (476 MHz) injected from an external antenna during helicon current drive experiments in DIII-D. Broadband density fluctuations around the helicon frequency are observed during real-time scans of measurement location and wavenumber during these experiments. Analysis indicates that these broadband high-frequency fluctuations are a result of backscattering of the DBS millimeter-wave probe beam from plasma turbulence modulated by the helicon wave. It is observed that background turbulence is effectively locally “tagged” with the helicon wave electric field, forming images of the turbulent spectrum in the overall density fluctuation spectrum that appear as high-frequency sidebands of the turbulence. These observations of background turbulence and high-frequency fluctuations open up the possibility of monitoring local helicon wave amplitude by comparing the high-frequency signal amplitude to the simultaneously measured background turbulence. In combination with the real-time measurement location and wavenumber scanning capabilities (offered by real-time frequency-stepping and steering), this allows rapid determination of the spatial distribution of the helicon wave power during steady-state plasma operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dampening of the Precipitation Response to Aerosol Pollution From Turbulence in Cumulus Clouds

Using aircraft observations and a cloud model with particle‐based microphysics, it is shown that the enhancement of drop collision‐coalescence from turbulence in clouds not only leads to earlier onset of rain in warm cumulus clouds as past studies have suggested, but also significantly dampens the precipitation susceptibility to aerosol loading. Enhanced drop coalescence from turbulence substantially increases the production of drizzle embryos just above cloud base, which in turn act as seeds that accelerate rain drop growth at mid and upper cloud levels even in highly polluted conditions. In contrast, pollution aerosols strongly inhibit rainfall when the commonly assumed gravitational‐only collision kernel is used and turbulent coalescence is neglected. Overall, turbulence‐enhanced drop coalescence strongly influences the response of warm cumulus clouds and precipitation to aerosol loading, suggesting that the effects of turbulent coalescence should be included in Earth system model representations of aerosol‐cloud‐precipitation interactions and aerosol indirect radiative forcing.

Chandrakar, Kamal Kant [NSF National Center for At↗