Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

An ignition criterion for inertial fusion boosted by microturbulence

Turbulence on fine spatial scales enhances fusion reactivity, enabling ignition at lower temperature. A modified Lawson-like ignition criterion is derived for inertially confined plasmas harboring turbulent kinetic energy. For some turbulent energy spectra, hot spots ignite at lower energy density and smaller volume. While detrimental mixing effects typically accompany turbulence and obscure these advantages, targets might be engineered to drive flow in regions where it is beneficial. The optimal length scale for this driving is identified, typically lying in the micrometer range.

Alpha particles↗
DOE OSTI2021

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗
DOE OSTI2025

Nonlinear simulation of under-resolved flows with shocks

Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

100 kHz OH PLIF of High-Pressure H2/NH3 and H2/CH4 Flames

The reaction zones in a multi-element, multi-stage, micromix injector are investigated using planar laser-induced fluorescence (PLIF) imaging of the hydroxyl (OH) radical at a 100 kHz repetition rate. The combustor is operated at a pressure of ∼10 bar. Multiple fuel blends of hydrogen, methane, and ammonia are considered to study the effect of changing fuel composition on the flame structure. Compact regions of high OH signal are observed in each flame for fuel mixtures with high hydrogen content, indicating rapid reactant consumption in the shear layers. Instantaneous images exhibit wrinkled flame fronts and the formation of unreacted pockets due to highly turbulent burning. The addition of methane is observed to reduce the reactivity and flame speed, resulting in more diffuse regions of OH between each flame. The shape of the reactant core is altered, and the length scale of the flame surface wrinkling is increased. Ammonia addition has a more severe effect on the distribution of OH and the local flame structure than methane addition. A time series of images shows the presence of OH vortical structures and significant turbulent fluctuations in the flames near the wall, suggesting the presence of a recirculation zone that could drive the observed disparity in flame lengths.

Hodge, Alexander J. [Purdue University]↗
DOE OSTI2026

Wind Conditions for Solar Energy Facility Design Based on Field Measurements at Nevada Solar One

This study characterises near-surface wind conditions at the Nevada Solar One (NSO) concentrated solar power plant using two years of measurements from a 15-m mast located 30 m west of the solar arrays. The surrounding terrain is fairly flat within 3 km, but nearby solar photovoltaic (PV) and concentrated solar power (CSP) installations introduce strong directional wind inhomogeneity, which motivates a sector-based analysis that distinguishes undisturbed inflow from waked flow. The analysis uses data from three sonic anemometers to assess whether turbulence in the surface layer can be represented as Gaussian and whether standard surface-layer scaling remains valid below 10 m height. The roughness length and displacement height are estimated for twelve directional sectors under near-neutral wind conditions, and all three velocity components show departures from Gaussian behaviour. Spectral analysis indicates that the turbulence contains low-frequency fluctuations with amplitudes higher than those predicted by standard models such as the Kaimal model or the Simiu and Scanlan model. The study introduces a new turbulence spectral model based on the eddy-surface-layer framework to better represent near-surface wind turbulence. The analysis further quantifies the influence of the surrounding power plants on mean wind profiles, one-point spectra, and vertical coherence for all velocity components. The solar collector arrays located 30 m east of the mast significantly modify the turbulence structure, whereas PV installations situated a few hundred metres away have only a limited influence on the turbulence measured at the mast. These findings highlight the limitations of conventional turbulence models for near-surface wind over solar infrastructure and provide improved parameterisations to support the structural design. The implications for the structural response of parabolic troughs and PV modules to such turbulence remain uncertain and require further investigation.

14 SOLAR ENERGY↗
DOE OSTI2026

Pixel-Resolved Long-Context Learning for Turbulence at Exascale: Resolving Small-scale Eddies Toward the Viscous Limit

Turbulence plays a crucial role in multiphysics applications, including aerodynamics, fusion, and combustion. Accurately capturing turbulence's multiscale characteristics is essential for reliable predictions of multiphysics interactions, but remains a grand challenge even for exascale supercomputers and advanced deep learning models. The extreme-resolution data required to represent turbulence, ranging from billions to trillions of grid points, pose prohibitive computational costs for models based on architectures like vision transformers. To address this challenge, we introduce a multiscale hierarchical Turbulence Transformer that reduces sequence length from billions to a few millions and a novel RingX sequence parallelism approach that enables scalable long-context learning. We perform scaling and science runs on the Frontier supercomputer. Our approach demonstrates excellent performance up to 1.1 EFLOPS on 32,768 AMD GPUs, with a scaling efficiency of 94\%. To our knowledge, this is the first AI model for turbulence that can capture small-scale eddies down to the dissipative range in three-dimensional turbulence at high Reynolds numbers.

Yin, Junqi [ORNL] (ORCID:0000000338435520)↗
DOE OSTI2023

Skyrmion‐Excited Spin‐Wave Fractal Networks

Magnetic skyrmions exhibit unique, technologically relevant pseudo-particle behaviors which arise from their topological protection, including well-defined, 3D dynamic modes that occur at microwave frequencies. During dynamic excitation, spin waves are ejected into the interstitial regions between skyrmions, creating the magnetic equivalent of a turbulent sea. However, since the spin waves in these systems have a well-defined length scale, and the skyrmions are on an ordered lattice, ordered structures from spin-wave interference can precipitate from the chaos. This work uses small-angle neutron scattering (SANS) to capture the dynamics in hybrid skyrmions and investigate the spin-wave structure. Performing simultaneous ferromagnetic resonance and SANS, the diffraction pattern shows a large increase in low-angle scattering intensity, which is present only in the resonance condition. This scattering pattern is best fit using a mass fractal model, which suggests the spin waves form a long-range fractal network. The fractal structure is constructed of fundamental units with a size that encodes the spin-wave emissions and are constrained by the skyrmion lattice. Furthermore, these results offer critical insights into the nanoscale dynamics of skyrmions, identify a new dynamic spin-wave fractal structure, and demonstrate SANS as a unique tool to probe high-speed dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal fluctuations in the dissipation range of homogeneous isotropic turbulence

In this work, using fluctuating hydrodynamics we investigate the effect of thermal fluctuations in the dissipation range of homogeneous isotropic turbulence. Simulations confirm theoretical predictions that the energy spectrum is dominated by these fluctuations at length scales comparable to the Kolmogorov length. We also find that the extreme intermittency in the far-dissipation range predicted by Kraichnan is replaced by Gaussian thermal equipartition.

42 ENGINEERING↗
DOE OSTI2020

Evaluating turbulence models at high energy densities

The RESHOCK campaign at Lawrence Livermore National Laboratory has been working on increasing our understanding of the evolution of turbulent, unstable plasma interfaces which are applicable in various HED applications including ICF implosions. A common approach to model such interfaces involves the use of Renolds-Averaged-Navier-Stokes (RANS) models whose parameters have been constrained by theory and experiments. These parameter sets are not unique and the applicability of one set of tuned model parameters to systems with differing shock strengths or material densities has not been well tested. Our new NIF experiments, along with related studies of turbulent shear flows, provide data and model validation for plasma interfaces in the HED regime. Our experiments measure mixing-layer width at an unstable interface in the plasma regime and are specifically designed to challenge mix models of the Reynolds-Averaged-Navier-Stokes (RANS) type where turbulence is assumed to be fully developed, i.e. the flow exhibits a broad spectrum of length scales without memory of the initial condition.4 We utilize precise control of the initial interface conditions (densities on both sides of the interface, an initial rippled perturbation pattern consisting of equally weighted wavelengths between 10 and 20 microns, and initial material compositions) along with a repeatable drive history

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hierarchical-embedding autoencoder with a predictor as efficient architecture for learning time-evolution in multi-scale turbulent flows

We introduce a scale-aware, data-driven deep learning modeling framework for accurately predicting the time evolution of multi-scale turbulent plasma and liquid flows. The approach is motivated by the idea of scale separation. Structures of vastly different length scales emerge in these systems, and interactions between these structures occur only locally. To exploit this structure, the flow state is transformed by a hierarchical, fully convolutional autoencoder, not into a single embedding layer as in conventional convolutional surrogate models, but into a series of embedding layers. A stepwise training strategy ensures that fine-scale features are encoded on a high-resolution grid, while larger structures are represented on progressively coarser layers. The time evolution predictor advances all embedding layers in sync, capturing local interactions between features at the same scale as well as between all scales. This approach enables efficient modeling of multi-scale systems since negligible interactions between distant, small-scale structures do not need to be directly modeled. Our hierarchical-embedding autoencoder with a predictor framework is evaluated on canonical examples of multi-scale turbulence: two-dimensional Kolmogorov flow and Hasegawa–Wakatani plasma turbulence. In both cases, the proposed framework significantly improves predictive accuracy relative to conventional convolutional network architectures. A significant improvement in prediction accuracy was observed for crucial statistical characteristics of the Hasegawa–Wakatani plasma as well as for individual trajectories of the Kolmogorov flow turbulence. Importantly, the model's rollout for the Hasegawa–Wakatani problem demonstrates a four-order-of-magnitude speedup compared to traditional numerical solvers.

Khrabry, Alexander I. [Princeton Univ., NJ (United↗
DOE OSTI2023

Adjustments to the law of the wall above an Amazon forest explained by a spectral link

Modification to the law of the wall represented by a dimensionless correction function ϕ RSL (z/h) is derived using atmospheric turbulence measurements collected at two sites in the Amazon in near-neutral stratification, where z is the distance from the forest floor and h is the mean canopy height. The sites are the Amazon Tall Tower Observatory for z/h∈ [1,2.3] and the Green Ocean Amazon (GoAmazon) site for z/h∈ [1,1.4]. Here, a link between the vertical velocity spectrum E ww (k) (k is the longitudinal wavenumber) and ϕ RSL is then established using a co-spectral budget (CSB) model interpreted by the moving-equilibrium hypothesis. The key finding is that ϕ RSL is determined by the ratio of two turbulent viscosities and is given as ν t,BL /ν t,RSL, where ν t,RSL = (1/A)∫$^{∞}_{0}$ τ(k)E ww (k)dk, ν t,BL = k v (z−d)u * , τ(k) is a scale-dependent decorrelation time scale between velocity components, A = C R /(1−C I ) = 4.5 is predicted from the Rotta constant C R = 1.8, and the isotropization of production constant C I = 3/5 given by rapid distortion theory, k v is the von Kármán constant, u * is the friction velocity at the canopy top, and d is the zero-plane displacement. Because the transfer of energy across scales is conserved in E ww (k) and is determined by the turbulent kinetic energy dissipation rate (ε), the CSB model also predicts that ϕ RSL scales with L BL /L d , where L BL is the length scale of attached eddies to z = d, and L d = u$_{*}^{3}$/ε is a macro-scale dissipation length.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Alignment of the magnetic field in star-forming regions and why it might be difficult to observe

ABSTRACT Magnetic fields are an important component of the interstellar medium (ISM) and exhibit strongly varying field strengths and a non-trivial correlation with the gas density. Its dynamical impact varies between individual regions of the ISM and correlates with the orientation of the field with respect to the gas structures. Using high-resolution magneto-hydrodynamical simulations of the ISM, we explore the connection between the orientation of the field and the dynamical state of the gas. We find that the onset of gravitational instability in molecular gas above a density of $\rho \sim 10^{-21}\, \mathrm{g\, cm}^{-3}$$(n\sim 400\, \mathrm{cm}^{-3})$ coincides with an alignment of the magnetic field lines and the gas flow. At this transition, the gradient of the density changes from mainly perpendicular to preferentially parallel to the field lines. A connection between the three-dimensional alignment and projected two-dimensional observables is non-trivial, because of a large dispersion of the magnetic field orientation along the line of sight. The turbulent correlation lengths can be small compared to the typical integration lengths. As a consequence, the small-scale signal of the orientation can sensitively depend on the line of sight or the dynamical state of the cloud can fluctuate stochastically or be completely averaged out. With higher spatial resolution more small-scale structures are resolved, which aggravates the link between magneto-hydrodynamical quantities and projected observables.

Girichidis, Philipp (ORCID:0000000293009914)↗
DOE OSTI2025

Invariant discovery of features across multiple length scales: Applications in microscopy and autonomous materials characterization

Physical imaging is a foundational characterization method in areas from condensed matter physics and chemistry to astronomy and spans length scales from atomic to universe. Images encapsulate crucial data regarding atomic bonding, materials microstructures, and dynamic phenomena such as microstructural evolution and turbulence, among other phenomena. The challenge lies in effectively extracting and interpreting this information. Variational Autoencoders (VAEs) have emerged as powerful tools for identifying the underlying factors of variation in image data, providing a systematic approach to distilling meaningful patterns from complex data sets. However, a significant hurdle in their application is the definition and selection of appropriate descriptors reflecting local structures. Here, we introduce the scale-invariant VAE approach (SI-VAE) based on the progressive training of the VAE with the descriptors sampled at different length scales. The SI-VAE allows the discovery of the length scale-dependent factors of variation in the system. Here, we illustrate this approach using the ferroelectric domain images and generalize it to the movies of the electron-beam induced phenomena in graphene and topography evolution across combinatorial libraries. This approach can further be used to initialize the decision making in automated experiments including structure–property discovery and can be applied across a broad range of imaging methods. This approach is universal and can be applied to any spatially resolved data including both experimental imaging studies and simulations, and can be particularly useful for exploration of phenomena such as turbulence and scale-invariant transformation fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Scaling Behavior of a Turbulent Kinetic Energy Closure Scheme for the Stably Stratified Atmosphere: A Steady-State Analysis

Abstract We present a turbulent kinetic energy (TKE) closure scheme for the stably stratified atmosphere in which the mixing lengths for momentum and heat are not parameterized in the same manner. The key difference is that, while the mixing length for heat tends toward the stability independent mixing length for momentum in neutrally stratified conditions, it tends toward one based on the Brunt–Väisälä time scale and square root of the TKE in the limit of large stability. This enables a unique steady-state solution for TKE to be obtained, which we demonstrate would otherwise be impossible if the mixing lengths were the same. Despite the model’s relative simplicity, it is shown to perform reasonably well with observational data from the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) using commonly employed model constants. Analyzing the scaling behavior of the nondimensional velocity and potential temperature gradients, or of the stability (correction) functions, reveals that for large stability the present model scales in the same manner as the first-order operational scheme of Viterbo et al. Alternatively, it appears as a blend of two cases of the TKE closure scheme of Baas et al. Critically, because a unique steady-state TKE can be obtained, the present model avoids the nonphysical behavior identified in one of the cases of Baas et al.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Unveiling the Significance of Correlations in K-Space and Configuration Space for Drift Wave Turbulence in Tokamaks

Turbulence and transport phenomena play a crucial role in the confinement and stability of tokamak plasmas. Turbulent fluctuations in certain physical quantities, such as density or temperature fluctuations, can have a wide range of spatial scales, and understanding their correlation length is important for predicting and controlling the behavior of the plasma. The correlation length in the radial direction is identified as the critical length in real space. The dynamics in real space are of significant interest because transport in configuration space is primarily focused on them. When investigating transport caused by the E×B drift, the correlation length in real space represents the size of E×B whirls. It was numerically discovered that in drift wave turbulence, this length is inversely proportional to the normalized mode number of the fastest growing mode relative to the drift frequency. Considerable time was required before a proper analytical derivation of this condition was accomplished. Therefore, a connection has been established between phenomena occurring in real space and those occurring in k-space. Although accompanied by a turbulent spectrum in k-space with a substantial width, transport in real space is uniquely determined by the correlation length, allowing for accurate transport calculations through the dynamics of a single mode. Naturally, the dynamics are subject to nonlinear effects, with resonance broadening in frequency being the most significant nonlinear effect. Thus, mode number space is once again involved. Resonance broadening leads to the detuning of waves from particles, permitting a fluid treatment. It should be emphasized that the consideration here involves the total electric field, including the induction part, which becomes particularly important at higher beta plasmas.

Weiland, Jan (ORCID:0000000253496025)↗
DOE OSTI2022

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Neoclassical transport and profile prediction in transport barriers

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Trinczek, Silvia [Princeton Plasma Physics Laborat↗