Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak turbulence”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Staircase resiliency in a fluctuating cellular array

Inhomogeneous mixing by stationary convective cells set in a fixed array is a particularly simple route to layering. Layered profile structures, or staircases, have been observed in many systems, including drift-wave turbulence in magnetic confinement devices. The simplest type of staircase occurs in passive-scalar advection, due to the existence and interplay of two disparate timescales, the cell turn-over (τ H ), and the cell diffusion (τ D ) time. In this simple system, we study the resiliency of the staircase structure in the presence of global transverse shear and weak vortex scattering. The fixed cellular array is then generalized to a fluctuating vortex array in a series of numerical experiments. The focus is on regimes of low-modest effective Reynolds numbers, as found in magnetic fusion devices. By systematically perturbing the elements of the vortex array, we learn that staircases form and are resilient (although steps become less regular, due to cell mergers) over a broad range of Reynolds numbers. The criteria for resiliency are (a) τ D >>τ H and (b) a sufficiently high profile curvature (κ ≥ 1.5). We learn that scalar concentration travels along regions of shear, thus staircase barriers form first, and scalar concentration "homogenizes" in vortices later. The scattering of vortices induces a lower effective speed of scalar concentration front propagation. The paths are those of the least time. We observe that if background diffusion is kept fixed, the cell geometric properties can be used to derive an approximation for the effective diffusivity of the scalar. Furthermore, the effective diffusivity of the fluctuating vortex array does not deviate significantly from that of the fixed cellular array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

ECEI characterization of pedestal fluctuations in quiescent H-mode plasmas in DIII-D

Here, electron cyclotron emission imaging (ECEI) is employed to characterize the magneto hydraulics dynamics (MHD) fluctuations at the quiescent H-mode pedestals in DIII-D. Pedestal MHD fluctuations cause ECE radiation temperature fluctuations $\delta {T_{e{\text{,rad}}}}$ in both the pedestal and scrape-off-layer (SOL). A synthetic ECE platform is utilized for detailed interpretation of the ECEI signals in the SOL and pedestal regions. It is observed that the ECE radiation $\delta {T_{e{\text{,rad}}}}$, which is located in the SOL region according to the cold and optically thick plasma resonance assumption, is extremely sensitive to MHD radial displacements near the separatrix, exhibiting radiation inversion to $\delta {T_{e{\text{,rad}}}}$ at the pedestal. Here, the radiation inversion refers to the opposite phase between the radiation fluctuation at the pedestal and the radiation fluctuation at the SOL. Consequently, the quasi-coherent MHD (QCM), which displays a radiation inversion, is found to be consistent with an MHD radial structure that has a strong displacement near the separatrix. In contrast, the edge harmonic oscillation (EHO), which displays weak or no inversion, is found to be consistent with an MHD radial displacement structure peaking at the pedestal top. The ECEI data, interpreted with synthetic ECE, are in qualitative agreement with beam emission spectroscopy measurements on DIII-D for the relative radial extent and localization of the EHO and QCM. The high sensitivity of ECE radiation to separatrix displacements can be used to detect turbulence or MHD fluctuations near the separatrix, which may affect the transport across the separatrix and the wetted area in the divertor. The $\delta {T_{e{\text{,rad}}}}$ inversion measured with an ECE or ECEI system potentially provides important information on the magnetic field $\left| {{B_{{\text{sep}}}}} \right|$ at the separatrix, which helps constrain the pedestal equilibrium reconstruction and achieve an unambiguous mapping of the ECE/ECEI system with respect to the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Accuracy of volumetric flow rate inflow/outflow measurement by integrating PIV velocity fields

In this paper, the usefulness of planar Particle Image Velocimetry (PIV) measurements to determine the volumetric flow rate from a rectangular or round jet exit is assessed. Both two-component PIV (2C–PIV) and through-plane three-component (3C) stereo PIV (SPIV) data sets were acquired at the exit of a turbulent, rectangular nozzle and a round nozzle with Reynolds numbers between 10 000 and 100 000. The PIV data sets were processed using a variety of algorithms. The time-averaged results were then spatially integrated across the jet exits and compared to a calibrated flow meter. Recommendations for each method are developed and discussed with potential drawbacks. The accuracy of the measurement was found to be a weak function of the Reynolds number of the flow. Furthermore, two-component PIV was found to underestimate the volumetric flow rate by 1%–4% depending on the integration scheme and SPIV underestimated volumetric flow rate by 2%.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Multimodel simulation of vertical gas transfer in a temperate lake

In recent decades, several lake models of varying complexity have been developed and incorporated in numerical weather prediction systems and climate models. To foster enhanced forecasting ability and verification, improvement of these lake models remains essential. This especially applies to the limited simulation capabilities of biogeochemical processes in lakes and greenhouse gas exchanges with the atmosphere. Here we present multi-model simulations of physical variables and dissolved gas dynamics in a temperate lake (Harp Lake, Canada). The five models (ALBM, FLake, LAKE, LAKEoneD, MTCR-1) considered within this most recent round of the Lake Model Intercomparison Project (LakeMIP) all captured the seasonal temperature variability well. In contrast, none of the models is able to reproduce the exact dates of ice-cover and ice-off, leading to considerable errors in the simulation of eddy diffusivity around those dates.We then conducted an additional modeling experiment with a diffusing passive tracer to isolate the effect of the eddy diffusivity on gas concentration. Remarkably, sophisticated k -" models do not demonstrate a significant difference in the vertical diffusion of a passive tracer compared to models with much simpler turbulence closures. All models simulate less intensive spring overturn compared to autumn. Reduced mixing in the models consequently leads to the accumulation of the passive tracer distribution in the water column. The lake models with a comprehensive biogeochemical module, such as ALBM and LAKE, predict dissolved oxygen dynamics adequate to the observed data. However, for the surface carbon dioxide concentration the correlation between modeled (ALBM, LAKE) and observed data is weak (~0.3). Overall our results indicate the need to improve the representation of physical and biogeochemical processes in lake models, thereby contributing to enhanced weather prediction and climate projection capabilities.

Guseva, Sofya↗
DOE OSTI2026

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗
DOE OSTI2026

Influence of Operating Conditions on Ethanol Passive Prechamber Cold-start Operation

Cold-start operation in spark ignition engines is characterized by many drawbacks such as poor fuel vaporization, substantial wall film formation, and weak ignition which can result in unstable combustion and high concentrations of carbon monoxide (CO), nitrogen oxides (NOx), and unburnt hydrocarbons (UHC) prior to catalyst light-off. These challenges are particularly important for ethanol as the fuel’s high latent heat of vaporization and low vapor pressure can hinder reliable ignition during cold-start. In this context, passive prechamber ignition systems offer a viable pathway to improve ignitability by producing hot turbulent jets that can promote combustion in the main chamber. In this study, the influence of operating conditions on jet formation and combustion characteristics in an ethanol-fueled passive prechamber spark ignition engine is investigated. Two intake pressures (40 kPa and 60 kPa) and two engine speeds (450 and 600 rpm) are studied to represent various stages of the cold-start ramp-up. Results indicate that lower intake pressures lead to delayed and asymmetric turbulent jet formation, and lower peak heat release rates. Lower engine speeds, on the other hand, produced higher peak heat release and faster combustion due to greater residence time in the engine. This greater residence time also increases the likelihood of autoignition in the main chamber, potentially promoting quicker heat release for lower engine speeds. These results provide insight into the role of operating conditions on cold-start ramp-up of passive prechamber engines.

Banagiri, Shrikar [ORNL] (ORCID:0000000239745099)↗
DOE OSTI2021

Gravitomagnetic vorticity generation in black hole accretion discs: a potential spatial constraint on plasma flow stability

We calculate the vorticity generation rate in the accretion disc near a slowly rotating black hole in the low-velocity, weak-field limit of general relativity. Specifically, we find that the frame-dragging effect due to the black hole’s rotation – manifested through the gravitomagnetic field – can generate vorticity in a moving plasma in the accretion disc. The mechanism remains operational as long as the accretion disc has non-negligible vertical height and is independent of the exact thermodynamical profile of the disc. The enstrophy density generation rate, as a measure of turbulence and dissipation, is presented, which indicates that the frame-dragging effect can disrupt the stability of the disc away from the z = 0 plane.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

First-principles-integrated Study of Blazar Synchrotron Radiation and Polarization Signatures from Magnetic Turbulence

Blazar emission is dominated by nonthermal radiation processes that are highly variable across the entire electromagnetic spectrum. Turbulence, which can be a major source of nonthermal particle acceleration, can widely exist in the blazar emission region. The Turbulent Extreme Multi-Zone (TEMZ) model has been used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) simulations have also revealed the stochastic nature of the turbulent emission region and particle acceleration therein. However, radiation signatures have not been systematically studied via first-principles-integrated simulations. In this paper, we perform combined PIC and polarized radiative transfer simulations to study synchrotron emission from magnetic turbulence in the blazar emission region. We find that the multiwavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability timescale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle swings with arbitrary amplitudes and duration, in either direction, that are not associated with changes in flux or PD. Surprisingly, these swings, which are stochastic in nature, can appear either bumpy or smooth, although large-amplitude swings (>180°) are very rare, as expected. Our radiation and polarization signatures from first-principles-integrated simulations are consistent with the TEMZ model, except that in the latter, there is a weak correlation, with zero lag, between flux and degree of polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Effect of triangularity on ion-temperature-gradient-driven turbulence

The linear and nonlinear properties of ion-temperature-gradient-driven turbulence with adiabatic electrons are modeled for axisymmetric configurations for a broad range of triangularities δ, both negative and positive. Peak linear growth rates decrease with negative δ but increase and shift toward a finite radial wavenumber kx with positive δ. The growth-rate spectrum broadens as a function of kx with negative δ and significantly narrows with positive δ. Here, the effect of triangularity on linear instability properties can be explained through its impact on magnetic polarization and curvature. Nonlinear heat flux is weakly dependent on triangularity for |δ|≤0.5, decreasing significantly with extreme δ, regardless of sign. Zonal modes play an important role in nonlinear saturation in the configurations studied, and artificially suppressing zonal modes increased nonlinear heat flux by a factor of about four for negative δ, increasing with positive δ by almost a factor of 20. Proxies for zonal-flow damping and drive suggest that zonal flows are enhanced with increasing positive δ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observed Covariations in Boundary Layer and Cumulus Cloud Layer Processes

In this paper we examine variations in boundary-layer processes spanning the shallow to deep cumulus transition. This is accomplished by differentiating boundary layer properties on the basis of convective outcomes, ranging from shallow to deep, as observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma, USA. Doppler lidar, radar, and radiosonde data are combined to determine statistical differences in boundary layer and cloud layer properties using a large sample (236) of days with a range of convective outcomes: shallow, congestus, and deep convection. In these analyses, the radar characterizes diurnal cloud depth, the lidar quantifies updraft and downdraft properties in the subcloud layer, and daily radiosonde data provides the convective inhibition (CIN). Combined, these data are used to test the hypothesis that deep convection occurs when the strength of the boundary layer turbulence (i.e., TKE) exceeds the strength of the energy barrier (i.e., CIN) at the top of the CBL. Results show that days with deep convective clouds have significantly lower vertical velocity variance and weaker updrafts within the subcloud layer. However, CIN values are also found to be significantly lower on deep convective days, allowing for these weaker updrafts to penetrate the energy barrier and reach the level of free convection (LFC). In contrast, shallow convective outcomes occur when the updrafts are strong in an absolute sense, but are weak when compared to the strength of the energy barrier. These findings support the use of the CIN/TKE framework in parameterizing convection in coarse resolution models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Scalable Real-Time Data Assimilation Framework for Predicting Turbulent Atmosphere Dynamics

AI-based foundation models like FourCastNet, GraphCast are revolutionizing weather and climate predictions but are not yet ready for operational use. Their limitation lies in the absence of a data assimilation system to incorporate real-time Earth system observations, crucial for accurately forecasting events like tropical cyclones. To overcome these obstacles, we introduce a generic real-time data assimilation framework and demonstrate its end-to-end performance on the Frontier supercomputer. This framework comprises two primary modules: an ensemble score filter (EnSF), which significantly outperforms the state-of-the-art data assimilation method, and a vision transformer-based surrogate capable of real-time adaptation through the integration of observational data. We demonstrate both the strong and weak scaling of our framework up to 1024 GPUs on the Exascale supercomputer, Frontier. Our results not only illustrate the framework's exceptional scalability on high-performance computing systems, but also demonstrate the importance of supercomputers in real-time data assimilation for weather and climate predictions.

Lu, Dan↗
DOE OSTI2024

Improving GCM Predictability of Mixed-Phase Clouds and Aerosol Interactions at High Latitudes with ARM Observations

The goal of this project is to improve the predictability of mixed-phase clouds and aerosol interactions in the Community Atmosphere Model version 6 (CAM6) through comparison with the ARM observations. There are three main objectives of the proposed study: (1) Improve the representation of ice microphysical processes in mixed-phase clouds; (2) Test the performance of ice microphysics in CESM-CAM6 with the ARM observations in northern and southern high latitudes; and (3) Examine mixed-phase cloud microphysics-aerosol-turbulence-radiation interactions in CESM-CAM6. In this project, we have (1) Improved the representation of ice microphysical processes in mixed-phase clouds in CESM-CAM6 by implementing the marine organic aerosol (MOA) and treating the ice nucleating particles (INPs) from MOA and its impacts on mixed-phase clouds. We made a first attempt to represent different secondary ice production (SIP) mechanisms in a GCM (CESM2-CAM6). We found that misrepresentation of these ice formation processes in global climate models (GCMs) leads to too weak negative cloud feedback over the Southern Ocean (SO) and too high climate sensitivity in the models. In addition to ice formation processes, we improved the treatment of ice depositional growth through the Wegener–Bergeron–Findeisen (WBF) process by considering the subgrid heterogeneous distributions between liquid droplets and ice crystals in mixed-phase clouds; (2) Tested the performance of our improved representations of ice microphysics with the ARM observations at high latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3–3 × 10 4 M ⊙ yr −1 , typical of short-period compact-object binaries in common envelopes. Our models account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on the accretion rates, these systems can develop both proton- and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radiation trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominating during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that depend on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is conducive to creating a wide range of nuclear products through a variety of effects, including nuclear statistical equilibrium freeze-out and the r-, p-, and γ-processes. We also observe a novel multistep process in reheated trajectories, consisting of proton-capture and photodissociation reactions operating on r-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freeze-out timescales capable of making r-process elements with high overabundances through a disequilibrium effect between neutrons and α-particles that does not require an excess of neutrons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Environmental sensitivities of shallow-cumulus dilution – Part 1: Selected thermodynamic conditions

Cumulus entrainment, and its consequent dilution of buoyant cloud cores, strongly regulates the life cycle of shallow cumuli yet remains poorly understood. Herein, new insights into this problem are obtained through large-eddy simulations that systematically investigate the sensitivity of shallow-cumulus dilution to cloud-layer relative humidity (RH), cloud- and subcloud-layer depths, and continentality (i.e., the land–ocean contrast). The simulated cloud-core dilution is found to be strongly sensitive to continentality, with fractional dilution rates twice as large over the ocean as over land. Using a similarity theory based on the turbulent-kinetic-energy (TKE) budget, the reduced cloud-core dilution over land is attributed to larger cloud-base mass flux (mb), driven by stronger surface heating and subcloud turbulence. As mb increases, the fractional dilution rate must decrease to maintain energetic equilibrium. A positive sensitivity is also found to cloud-layer RH, with the core dilution increasing by 25 %–50 % for a 10 % enhancement in RH. This sensitivity is interpreted using the buoyancy-sorting hypothesis, in that mixtures of cloud and environmental air are more likely to become negatively buoyant and detrain (rather than diluting the cloud core) in drier cloud layers. By contrast, the sensitivities of (marine) shallow-cumulus dilution to cloud- and subcloud-layer depths are weak, with a 3 % decrease for a doubling for the former and a 4 % reduction in dilution for a 50 % deeper subcloud layer. These surprisingly weak sensitivities are readily explained by offsetting effects in the TKE similarity theory. Altogether, these experimental findings provide useful, though still incomplete, guidance for flow-dependent shallow-cumulus entrainment formulations in large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ~ 264 keV for 242 MeV calcium ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison with macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long, reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Scattering theory in noncanonical phase space: A Drift-Kinetic collision operator for weakly collisional plasmas

After developing a scattering theory for grazing collisions in general noncanonical phase spaces, we introduce a guiding center collision operator in five-dimensional phase space designed for plasma regimes characterized by long wavelengths (relative to the Larmor radius), low frequencies (relative to the cyclotron frequency), and weak collisionality (where repeated Coulomb collisions induce cumulatively small changes in particle magnetic moment). The collision operator is fully determined by the noncanonical Hamiltonian structure of guiding center dynamics and exhibits a metriplectic structure, ensuring the conservation of particle number, momentum, energy, and interior Casimir invariants. It also satisfies an H-theorem, allowing for deviations from an equilibrium Maxwellian distribution due to the nontrivial kernel of the noncanonical guiding center Poisson tensor, spanned by the magnetic moment. We propose that this collision operator and its underlying mathematical structure may offer valuable insight into the study of turbulence, transport, and self-organizing phenomena in both laboratory and astrophysical plasmas.

Hamiltonian mechanics↗
DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

intrinsic rotation↗