Search NASA⌕ Search

SEARCH · Search NASA

Results for “turbulent length scale”

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

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 OSTI2023

Analysis of real-fluid thermodynamic effects on turbulent statistics in transcritical channel flows

Wall-bounded turbulence at high-pressure transcritical conditions with intense density fluctuations are encountered in many technical applications. In this study, we analyze the turbulent energy transport in transcritical channel flows specifically focusing on dissipation rate, turbulent kinetic energy budgets, heat fluxes, and momentum-fluctuation statistics; results from this analysis are used to guide the development of turbulent scaling laws. Furthermore, we find that the dissipation rate of turbulent kinetic energy is dominated by the enstrophy in the logarithmic layer, and the fluctuating viscosity results in the reduced tilting of the vortical structures and the attenuation of streamwise vorticity in the near-wall layer; the fluctuating viscosity attenuates the dissipation rate by reducing the shear strain and the enstrophy production. Local equilibrium of the turbulent kinetic energy exists in the logarithmic layer. We show that the real-fluid thermodynamic effects significantly change the turbulent heat flux correlated with the sweep and the ejection events; the density changes alter the turbulent transport and result in noticeable magnitudes of density-fluctuation-related momentum-fluctuation statistics. From these results, scaling laws for the turbulent length scales and turbulent kinetic energy budgets are proposed, thereby contributing to improvement of the wall models in large-eddy simulations and Reynolds-averaged Navier-Stokes (RANS) simulations.

42 ENGINEERING↗
DOE OSTI2025

On Master-Length Scale Formulations for Stable Conditions in Turbulence Closure Models

Three formulations of the turbulence-length scales used in numerical modelling of atmospheric flows are compared. The comparison is made using the Mellor–Yamada–Nakanishi–Niino turbulence closure model within the stable boundary layer local similarity framework. With an appropriate choice of model constants, the model predictions are barely discernible and compare well with the empirical data obtained from the SHEBA campaign.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
DOE OSTI2023

Effect of the integral length scales of turbulent inflows on wind turbine loads

As wind turbines become larger, the fluctuations in the inflow become increasingly influential in the turbine structural loading. These fluctuations are characterized by the integral length scale, which represents the average size of the largest energy-containing turbulent eddies. Current design standards neglect the varying integral length scales that characterize inflows of wind turbines in operation. Using large-eddy simulations, we generate turbulent inflows of varying integral length scales and quantify the loads of the IEA 15-MW reference wind turbine. Results illustrate that the impact of turbulence on rotor and tower loads is up to 10 times greater than the impact of the mean shear profile. Increasing integral length scales from 0.3x to 0.5x the rotor diameter reduces blade root flapwise moments and rotor and tower loads. Increasing integral length scales from 0.5x to 1.7x the rotor diameter increases the rotor aerodynamic thrust force and the blade root flapwise shear loads and decreases the tower base torsional moment and the tilting and yawing rotor aerodynamic moments. Additionally, turbulence intensity has a greater impact on wind turbine loads than integral length scales. In conclusion, findings indicate that design standards should consider varying integral length scales for accurate wind turbine loading characterization in turbulent inflow conditions.

17 WIND ENERGY↗
DOE OSTI2022

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Visualization of post-detonation fireball flowfields and comparison to CFD modeling

Visualization of flow structures within post-detonation fireballs has been performed for benchmark validation of numerical simulations. Custom pressed PETN explosives with a 12-mm diameter hemispherical form factor were used to produce a spherically symmetric post-detonation flow with low soot yield. Hydroxyl-radical planar laser induce fluorescence (OH-PLIF) was employed to visualize the structure ranging from approximately 10 μs to 35 μs after shock breakout from the explosive pellet. Fireball simulations were performed using the HyBurn Computational Fluid Dynamics (CFD) package. Experimental OH-PLIF results were compared to synthetic OH-PLIF from post-processing of CFD simulations. From the comparison of experimental and synthetic OH-PLIF images, CFD is shown to replicate much of the flow structure observed in the experiments, revealing potential differences in turbulent length scales and OH kinetics. Furthermore, results provide significant advancement in experimental resolution of these harsh turbulent combustion environments and validate physical models thereof.

Detonator↗
DOE OSTI2024

Dynamic Wind Loading on CSP Collectors Caused by Turbulent Wind Fluctuations: Insights from a 2-Year Field Campaign

Concentrating Solar Power (CSP) is a promising solar technology for electricity generation with thermal energy storage and with the additional benefit of industrial heat production. Wind loading on CSP collector structures, such as parabolic troughs or heliostats, is one of the primary drivers of their structural design costs. In particular, dynamic wind loading is a major source of uncertainty in the collector design process, which heavily relies on wind tunnel testing. In the field, the turbulent nature of the incoming wind creates fluctuating loads (support structure loads and resulting mirror deflections) on the collectors, with impacts on fatigue lifetime and optical performance. As is well known, wind tunnel tests cannot entirely reproduce the complex turbulent wind conditions typically observed at full-scale plants. To shed light on this topic, NREL initiated a field campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs. We use the published 2-year dataset of high-resolution combined wind and structural loads measurements [1] to characterize the dynamic structural wind response. For quantifying dynamic wind loading, we apply the concept of admittance functions, which are spectral transfer functions that couple the turbulent wind to resulting structural loads (aerodynamic admittance), and to the structural response (mechanical admittance). In practice, aerodynamic admittance describes which turbulent eddy sizes are effective in creating structural loads. The mechanical admittance describes in which frequency ranges these loads are reinforced or dampened by the structure. While these functions are an established concept in civil engineering, their recent application to a single full-scale heliostat [2] proved their broader applicability to CSP collectors. Here, we present a characterization of admittance functions for full-scale parabolic trough collectors and show how wind characteristics (mean wind speed and direction, turbulent kinetic energy, turbulent length scales), the sun-tracking trough angle, and row position alter the admittance functions. Further, we study to which extent the admittance functions are universal for a specific trough geometry and how our findings compare to reported heliostat results. References [1] https://data.openei.org/submissions/5938. [2] Blume, K., Roger, M., and Pitz-Paal, R. 2023b. "Simplified analytical model to describe wind loads and wind-induced tracking deviations of heliostats." Solar Energy, 256, 96-109. https://doi.org/10.1016/j.solener.2023.03.055.

admittance functions↗
DOE OSTI2023

In Situ Wind and Turbulence Measurements in a Field of Full-Size Parabolic Trough Solar Collectors

Concentrated Solar Power (CSP) is a promising method for using Solar power for electricity generation with thermal energy storage. One of the primary drivers of structural design costs of CSP collector structures is wind loading. To date, the design of these structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights on structural dynamic response, and generating a first-of-a-kind, comprehensive, high-resolution wind-loading dataset available for validating simulations of wind loading on collector structures. The measurements at NSO consist of Sonic anemometers on masts at different heights to characterize the incoming flow and conditions within and above the trough field, in combination with a Doppler Lidar scanning the horizontal plane above the troughs. The wind measurements at NSO have been continuously operating since October 2021 and provide a year-long dataset characterizing wind and turbulence conditions. The structural load measurements start in November 2022 and will complement the wind measurements. In this poster, we present first results of the wind measurement campaign by highlighting days with different atmospheric flow conditions. We identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. The highest loads are expected when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest loads and block the subsequent rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales change. These affected conditions produce unique load cases on the structures that will be captured by the load measurements. If the wind blows along the trough rows, wind and turbulence conditions are less influenced by the troughs.

atmospheric turbulence↗
DOE OSTI2023

In Situ Wind and Turbulence Measurements in a Field of Full-Size Parabolic Trough Collectors

Concentrated Solar Power (CSP) is a promising method for using Solar power for electricity generation with thermal energy storage. One of the primary drivers of structural design costs of CSP collector structures is wind loading. To date, the design of these structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights on structural dynamic response, and generating a first-of-a-kind, comprehensive, high-resolution wind-loading dataset available for validating simulations of wind loading on collector structures. The measurements at NSO consist of Sonic anemometers on masts at different heights to characterize the incoming flow and conditions within and above the trough field, in combination with a Doppler Lidar scanning the horizontal plane above the troughs. The wind measurements at NSO have been continuously operating since October 2021 and provide a year-long dataset characterizing wind and turbulence conditions. The structural load measurements start in November 2022 and will complement the wind measurements. In this poster, we present first results of the wind measurement campaign by highlighting days with different atmospheric flow conditions. We identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. The highest loads are expected when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest loads and block the subsequent rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales change. These affected conditions produce unique load cases on the structures that will be captured by the load measurements. If the wind blows along the trough rows, wind and turbulence conditions are less influenced by the troughs.

concentrated solar power↗
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
DOE OSTI2024

DRDMannTurb: A Python package for scalable, data-driven synthetic turbulence

Synthetic turbulence models (STMs) are used in wind engineering to generate realistic flow fields and are employed as inputs to industrial wind simulations. Examples include prescribing inlet conditions in large eddy simulations that model loads on wind turbines and tall buildings. We are interested in STMs capable of generating fluctuations based on prescribed second-moment statistics since such models can simulate environmental conditions that closely resemble on-site observations. To this end, the widely used Mann model (see Mann, 1994, 1998) is the inspiration for DRDMannTurb. The Mann model is described by three physical parameters: a magnitude parameter influencing the global variance of the wind field and corresponding to the Kolmogorov constant multiplied by the rate of viscous dissipation of the turbulent kinetic energy to the two-thirds, αϵ 2/3 , a turbulence length scale parameter L, and a nondimensional parameter Γ related to the lifetime of the eddies. A number of studies, as well as international standards (e.g., those by the International Electrotechnical Commission (IEC)), include recommended values for these three parameters with the goal of standardizing wind simulations according to observed energy spectra. Yet, having only three parameters, the Mann model faces limitations in accurately representing the diversity of observable spectra. This Python package enables users to extend the Mann model and more accurately fit field measurements through flexible neural network models of the eddy lifetime function. Following Keith et al. (2021), we refer to this class of models as Deep Rapid Distortion (DRD) models. DRDMannTurb also includes a general module implementing an efficient method for synthetic turbulence generation based on a domain decomposition technique. This technique is also described in Keith et al. (2021).

17 WIND ENERGY↗
DOE OSTI2023

Wind Loading on Parabolic Trough Collectors: Wind and Structural Loads Measurements at an Operational Powerplant

Concentrated Solar Power (CSP) is a promising method for using solar power for electricity generation with thermal energy storage for industrial applications. Solar collectors constitute almost 1/3 of the total cost of the power plant. One of the primary drivers of unrealiability of these collectors is wind-driven loading of mirrors, support structures, and drives. To date, the design of the solar collector structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) parabolic trough powerplant. At this plant, parabolic trough solar collectors track the sun from east to west in the course of a day and face varying wind loads depending on the wind properties and the angle of the troughs. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights into structural dynamic response, and generating a comprehensive wind-loading dataset for validating simulations of wind loading on collector structures. The measurements at NSO consist of sonic anemometers on masts at different heights to characterize the incoming flow and conditions at four trough rows at the edge of the trough field. In addition, a Doppler Lidar scans the horizontal plane above the troughs. The wind measurements at NSO have been continuously collecting data since October 2021 and are combined with structural load measurements that started in November 2022. The load measurements were installed on the same four outermost trough rows and include support structure bending moments, drive torque moments, dynamic accelerations of the spaceframe, mirror displacement, and tilt angles. Measurements are planned to continue until May 2023, providing a first-of-a-kind, high-resolution multi-month dataset of combined wind and load measurements. In this presentation, we show first results of the measurement campaign. Based on the measurements, we identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. Most interactions between incoming wind and the trough field are observed when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest wind speed and block the downwind rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales are smaller after the first row. This leads to the highest static loads (bending and torque moments) at the first row but potentially increased dynamic loads within the field. We illustrate our findings with case studies focused on different wind conditions and will present fatigue analysis to highlight the impact of wind-driven loads on collector structures.

CSP collectors↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

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↗