Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence simulation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modification of ion-temperature-gradient turbulence by impurities in stellarator plasmas

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Calvo, Iván [Research Centre for Energy, Environme↗
DOE OSTI2020

Simulating Real Atmospheric Boundary Layers at Gray-Zone Resolutions: How Do Currently Available Turbulence Parameterizations Perform?

Recent computational and modeling advances have led a diverse modeling community to experiment with atmospheric boundary layer (ABL) simulations at subkilometer horizontal scales. Accurately parameterizing turbulence at these scales is a complex problem. The modeling solutions proposed to date are still in the development phase and remain largely unvalidated. This work assesses the performance of methods currently available in the Weather Research and Forecasting (WRF) model to represent ABL turbulence at a gray-zone grid spacing of 333 m. We consider three one-dimensional boundary layer parameterizations (MYNN, YSU and Shin-Hong) and coarse large-eddy simulations (LES). The reference dataset consists of five real-case simulations performed with WRF-LES nested down to 25 m. Results reveal that users should refrain from coarse LES and favor the scale-aware, Shin-Hong parameterization over traditional one-dimensional schemes. Overall, the spread in model performance is large for the cellular convection regime corresponding to the majority of our cases, with coarse LES overestimating turbulent energy across scales and YSU underestimating it and failing to reproduce its horizontal structure. Despite yielding the best results, the Shin-Hong scheme overestimates the effect of grid dependence on turbulent transport, highlighting the outstanding need for improved solutions to seamlessly parameterize turbulence across scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Planck Dust Polarization Power Spectra Are Consistent with Strongly Supersonic Turbulence

Abstract The polarization of the cosmic microwave background is rich in information but obscured by foreground emission from the Milky Way’s interstellar medium (ISM). To uncover relationships between the underlying turbulent ISM and the foreground power spectra, we simulated a suite of driven, magnetized, turbulent models of the ISM, varying the fluid properties via the sonic Mach number, M S , and magnetic (Alfvén) Mach number, M A . We measure the power spectra of density ( ρ ), velocity ( v ), magnetic field ( H ), total projected intensity ( T ), parity-even polarization ( E ), and parity-odd polarization ( B ). We find that the slopes of all six quantities increase with M S . Most increase with M A , while the magnetic field spectrum steepens with M A . By comparing spectral slopes of E and B to those measured by Planck, we infer typical values of M S and M A for the ISM. As the fluid velocity increases, M S > 4 , the ratio of BB power to EE power increases to approach a constant value near the Planck-observed value of ∼0.5, regardless of the magnetic field strength. We also examine correlation coefficients between projected quantities, and find that r TE ≈ 0.3, in agreement with Planck, for appropriate combinations of M S and M A . Finally, we consider parity-violating correlations r TB and r EB .

Astronomy & Astrophysics↗
DOE OSTI2025

Assessment and Improvement of the SST-Gamma Transition Model in Nalu-Wind

We conduct laminar–turbulent boundary-layer transition simulations using a local correlation-based transition model for two-dimensional incompressible flow and present enhancements to improve the accuracy of transition predictions. Menter’s Galilean-invariant 𝛾 transition model is implemented in the incompressible, unstructured-grid flow solver Nalu-Wind and is validated against experimental data and results from NASA’s flow solvers. The test cases of the AIAA Transition Prediction and Modeling Workshop are investigated, namely, the T3A/T3B flat plates and the NLF(1)-0416 and S809 airfoils. Based on the results, best practices for transition simulations, particularly for an unstructured-grid flow solver, are identified. Additional airfoil simulations are conducted for two wind turbine airfoils, S822 at Reynolds numbers of 𝑂⁡(10 5 ) and DU00-W-212 at Reynolds numbers of 𝑂⁡(10 7 ), to assess the model at low and high Reynolds numbers. Furthermore, through this work, we propose several approaches to enhance transition simulations, including 1) enforcing positivity of the implicit operator for the source terms of the transition model, 2) employing a constant turbulence intensity in stationary external flow simulations, and 3) recommending meshing for unstructured-grid flow solvers. Finally, we provide detailed documentation of the validation and data for the canonical cases to the transition modeling community.

17 WIND ENERGY↗
DOE OSTI2022

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Impact of Cloud-Base Turbulence on CCN Activation: CCN Distribution

Abstract Following our previous investigation of the turbulence impact on cloud-base single-size CCN activation, this study considers a similar problem assuming CCN size distribution obtained from field measurements. The total CNN concentration is taken as either 200 cm −3 to represent clean conditions, or as 2000 cm −3 to represent polluted conditions. CCN is assumed to be sodium chloride. The CCN activation in the rising nonturbulent adiabatic parcel is contrasted with the activation within a rising adiabatic parcel filled with inertial-range homogeneous isotropic turbulence. The turbulent parcel of 64 3 m 3 and the turbulent kinetic energy dissipation rate of 10 −3 m −2 s −3 are used in most of the simulations. Results for a range of mean parcel ascent rates, between 0.125 and 8 m s −1 , are discussed. Overall, the adiabatic turbulent parcel simulations show results consistent with the adiabatic nonturbulent parcel, with higher activated CCN concentrations for stronger parcel ascent rates. The key difference is a blurriness of the separation between dry CCN size bins featuring activated and nonactivated (haze) CCN, especially for weak mean ascent rates. The blurriness comes from CCN getting activated and subsequently deactivated in the fluctuating supersaturation field, instead of all becoming cloud droplets above the cloud base. This leads to significantly larger spectral widths in turbulent parcel simulations compared to the nonturbulent parcel when activation is completed. Modeling results are discussed in the context of the impact of turbulent fluctuations on CCN activation documented in laboratory experiments using the Pi chamber.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Aeroelastic code comparison using the IEA 22MW reference turbine

Reference wind turbine designs and the associated aeroelastic models are widely used in both research and industry. Reference models representing future concepts are of particular interest. Current state of the art aeroelastic tools are relied upon to design the next generation of large wind turbines. However, modelling assumptions may be invalidated by upcoming very large turbines, and different aeroelastic tools may give inconsistent results. A 22MW turbine model has been defined as part of International Energy Agency (IEA) Wind Task 55 on Reference Wind Turbines and Farms to represent future turbines to be deployed in the 2030s. In this study, an aeroelastic model of this turbine has been created in four tools; Bladed, HAWC2, OpenFAST, and QBlade. Code comparisons are presented for steady state operation, linear stability analysis, and time domain power production simulations in steady and turbulent wind. Generally, the codes show a good agreement, but with some differences present in the linear stability analysis, periodic azimuthal variation, and time domain simulations. The models are a good basis for further study with the IEA 22MW turbine, and further code comparison exercises.

17 WIND ENERGY↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Ion and Electron Acceleration in Fully Kinetic Plasma Turbulence

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-β plasmas (β is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher β. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions stay roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2021

A Parameterization of Turbulent Dissipation and Pressure Damping Time Scales in Stably Stratified Inversions, and its Effects on Low Clouds in Global Simulations

Abstract It is difficult for coarse‐resolution global models of the atmosphere to accurately simulate the observed distribution of low clouds. In particular, it is difficult for moist turbulence closure models to simulate sufficiently bright near‐coastal stratocumulus (Sc) without simulating overly bright marine shallow cumuli (Cu). To parameterize bright Sc, a turbulence parameterization must damp the turbulent fluxes of heat and moisture above cloud top in order to prevent excessive entrainment of dry air into cloud top. To parameterize dim shallow Cu, the subgrid variances of temperature and moisture must remain large, in order to permit partial cloudiness. However, damping the fluxes but not the variances just above cloud top is difficult if a parameterization uses a single “master” time scale to damp both. In nature, the above‐cloud fluxes are damped by pressure fluctuations, whereas scalar variances are damped by a different process, namely, turbulent dissipation. In a stably stratified inversion above cloud, pressure damping is large but turbulent dissipation is small. To avoid this problem, a multitime scale parameterization for damping has been developed. The damping parameterization has been implemented in a global model and evaluated. The parameterization is capable of dimming shallow Cu while producing adequately bright Sc.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models

Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the three-dimensional atmospheric boundary layer. Here, we conduct a numerical analogue of a field campaign with spatially limited measurements using large-eddy simulation. We pose flow reconstruction as an inpainting problem, and reconstruct realistic samples of turbulent, three-dimensional flow with the use of a latent diffusion model. The diffusion model generates physically plausible turbulent structures on larger spatial scales, even when input observations cover less than 1% of the volume. Through a combination of qualitative visualization and quantitative assessment, we demonstrate that the diffusion model generates meaningfully diverse samples when conditioned on just one observation. These samples successfully serve as initial conditions for a large-eddy simulation code. We find that diffusion models show promise and potential for other applications for other turbulent flow reconstruction problems.

17 WIND ENERGY↗
DOE OSTI2021

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Towards robust laser beam propagation in atmospheric turbulence

High-fidelity optical propagation through the atmosphere is essential for free-space optical technologies, including laser-based remote sensing and optical communication. However, atmospheric turbulence severely distorts beams and compromises system performance. In this work, we employ hypergeometric-Gaussian (HyGG) vortex beams as probes to characterize and mitigate atmospheric turbulence. Using over 250,000 experimental and simulated frames, we show that refining the power spectrum density (PSD) can reduce numerical prediction errors by up to 79.8%. Concurrently, experimental observations supported by numerical simulations demonstrate that HyGG beams exhibit superior turbulence resilience across multiple metrics compared to conventional Gaussian beams, particularly in their ability to withstand over 5 times stronger turbulence while maintaining similar intensity fluctuations. These dual investigations, on both turbulence mitigation and robust beam solutions, converge to form a unified strategy for enhancing free-space optical system performance. Collectively, our findings provide new insights into light–turbulence interactions and highlight the practical utility of vortex beams under atmospheric conditions.

Zhang, Boyu↗
DOE OSTI2025

Structural uncertainty assessment for low-Mach wall-resolved large-eddy simulations: Plane Channel and Periodic Hill Use Cases

A structural uncertainty assessment is provided for a set of wall-resolved large-eddy simulations (WRLES) using a turbulent channel flow (𝑅𝑒 𝜏 395) and a periodic hill in the turbulent regime (𝑅𝑒 𝑏 2800) validation suite. The standard one-equation subgrid scale turbulent kinetic energy model (𝑘 𝑠𝑔𝑠 ) using both a static and dynamic coefficient approach along with Wall Adapting Local Eddy Viscosity (WALE) are exercised. The 𝑘 𝑠𝑔𝑠 model activated in this study also allows for a correction that supports theoretical turbulent viscosity and total kinetic energy scaling (cubit and quadratic, respectively) in the near-wall regime. Channel and hill results showcase very good agreement between prediction and simulation for the WALE and the dynamic coefficient 𝑘 𝑠𝑔𝑠 model. The study also demonstrates that near-wall damping improves predictivity for the static coefficient 𝑘 𝑠𝑔𝑠 model, while proving an unnecessary additional contribution for the dynamic coefficient 𝑘 𝑠𝑔𝑠 model as near-wall scaling is naturally captured. Overall conclusions are that this suite of WRLES models are well suited for production usage in the Sierra-based low-Mach codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗
DOE OSTI2023

Strongly Heated Turbulent Flow in a Channel with Pin Fins

Large-eddy simulations (LES) were performed to study the turbulent flow in a channel of height H with a staggered array of pin fins with diameter D = H/2 as a function of heating loads that are relevant to the cooling of turbine blades and vanes. The following three heating loads were investigated—wall-to-coolant temperatures of Tw/Tc = 1.01, 2.0, and 4.0 - where the Reynolds number at the channel inlet was 10,000 and the back pressure at the channel outlet was 1 bar. For the LES, two different subgrid-scale models—the dynamic kinetic energy model (DKEM) and the wall-adapting local eddy-viscosity model (WALE)—were examined and compared. This study was validated by comparing with data from direct numerical simulation and experimental measurements. The results obtained show high heating loads to create wall jets next to all heated surfaces that significantly alter the structure of the turbulent flow. Results generated on effects of heat loads on the mean and fluctuating components of velocity and temperature, turbulent kinetic energy, the anisotropy of the Reynolds stresses, and velocity-temperature correlations can be used to improve existing RANS models.

42 ENGINEERING↗
DOE OSTI2022

Large eddy simulation of flow through an axisymmetric sudden expansion

Here, this study aims at investigating the inlet flow conditions of flow through an axisymmetric sudden expansion with an expansion ratio of 2.0. A series of large eddy simulations with the WALE model were conducted for different inlet Reynolds numbers ( Re) and turbulence intensities ($u_{rms}/\bar{U}_m$). The reattachment length, defined as the length measured downstream of the expansion where the flow direction is reversed adjacent to the wall ( Lr), was measured for each case. For widely studied inlet turbulence intensity values (TI), the simulation results are in good agreement with the experimental and numerical results reported in the literature. Parametric studies revealed that turbulence intensity affects the critical Reynolds number, marking the transition between the laminar and transition regions and the reattachment length. The critical Reynolds number was found to decrease with increasing turbulence intensity. A correlation expression is proposed. Additional analysis with proper orthogonal decomposition was performed to enhance the understanding of complex flow structures downstream of the expansion. Finally, an overall correlation expression for the reattachment length was obtained for 500 ≤ Re ≤15 000 and 0.2 ≤ TI (%) ≤ 20. For a given turbulence intensity, the reattachment length can be expressed for laminar and turbulent regions as a function of the Reynolds number. The reattachment length in the transition region can be expressed as a fractional average of reattachment lengths for laminar and turbulent flows.

99 GENERAL AND MISCELLANEOUS↗