SEARCH · Engineering Papers

Results for “dissipation rate”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

DNS of ignition and flame stabilization in a simplified gas turbine premixer

With the increasing need for fuel flexibility, mitigation of auto-ignition (AI) inside gas turbine (GT) premixers becomes crucial. They must be designed to yield a sufficiently homogeneous fuel-air mixture to achieve low emissions while at the same time avoiding the occurrence of AI and subsequent flame stabilization. This challenge requires a detailed understanding of turbulent mixing and chemistry interactions. In the present work, a direct numerical simulation (DNS) of an array of jets in crossflow (JICF), representative of an industrial GT premixer, is reported to shed light on these complex phenomena. It is found that AI kernels form in the aft part of the premixer and coalesce into a flame front that then propagates upstream, mainly through the boundary layer, and successively engulfs the jets. This, therefore, suggests a significant role of the jet array pattern on the flame stabilization. It is noted that AI kernels continue to form independently during the whole time of the simulation. To clarify the contribution of AI and diffusion in the ignition kernels and the main flame, chemical explosive mode analysis (CEMA) is employed jointly with a kernel tracking algorithm. It is found that during the initial formation of the flame, many ignition kernels form in mixtures with low scalar dissipation rate and large contribution from AI mode. As they quickly grow, they merge into a single flame front that becomes increasingly more diffusion-assisted over time, balancing the AI mode. Turbulence is shown to have a significant enhancing effect in lean premixed flames, but further analysis is required to fully characterize it. These findings are relevant for the industrial premixer studied, and also for novel micromix concepts that may be used in the next generation of GT combustion systems.

ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Characterizing in-stream turbulent flow for tidal energy converter siting in Cook Inlet, Alaska

Cook Inlet in Alaska is the most promising location for tidal energy development in the U.S. due to its significant tidal range of approximately 10 meters and high volume flux. The inlet's unique geometry and flow characteristics make it the most energetic tidal stream in the nation, with GW-scale potential energy capacity. With the growing interest in tidal energy converter (TEC) deployment in this area, we implemented a regional-scale, 3D hydrodynamic modeling framework to predict tidal current and turbulence characteristics that can assist TEC designers and project managers. We validated the model results extensively using various datasets collected with bottom-mounted acoustic Doppler current profilers and velocimeters. The comparison between the model outputs and observational data highlighted the effectiveness of the 3D FVCOM model and the Mellor-Yamada Level 2.5 Turbulence Model in accurately assessing macro-scale kinetic energy, turbulence intensity, and the production and dissipation rates at a prospective TEC site. Using two months of model simulation data, we examined the channel cross-section for TEC deployment, focusing on undisturbed power density and macro-scale turbulent properties. Further, our findings indicate that understanding the turbulence characteristics and flow properties can enhance Stage I/II resource characterization by identifying optimal locations for TECs and their layouts within the channel. Furthermore, we demonstrated that TEC designers can utilize macro-scale turbulence data from 3D coastal models as boundary conditions for other turbulence models, allowing for a more detailed resolution of the turbulence structure at TEC siting locations. Ultimately, this work emphasizes the importance of estimating flow and turbulence conditions in energetic systems to understand turbulent sites better and improve resource characterization.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗

DOE OSTI2022

Low-Reynolds-number oscillating boundary layers on adiabatic slopes

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

42 ENGINEERING↗

DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗

DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗

DOE OSTI2025

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗

DOE OSTI2025

A Computational Framework for Simulations of Dissipative Nonadiabatic Dynamics on Hybrid Oscillator-Qubit Quantum Devices

Here, we introduce a computational framework for simulating nonadiabatic vibronic dynamics on circuit quantum electrodynamics (cQED) platforms. Our approach leverages hybrid oscillator-qubit quantum hardware with midcircuit measurements and resets, enabling the incorporation of environmental effects such as dissipation and dephasing. To demonstrate its capabilities, we simulate energy transfer dynamics in a triad model of photosynthetic chromophores inspired by natural antenna systems. We specifically investigate the role of dissipation during the relaxation dynamics following photoexcitation, where electronic transitions are coupled to the evolution of quantum vibrational modes. Our results indicate that hybrid oscillator-qubit devices, operating with noise levels below the intrinsic dissipation rates of typical molecular antenna systems, can achieve the simulation fidelity required for practical computations on near-term and early fault-tolerant quantum computing platforms.

Hamiltonians↗

DOE OSTI2024

Gas Transfer Across Air‐Water Interfaces in Inland Waters: From Micro‐Eddies to Super‐Statistics

In inland water covering lakes, reservoirs, and ponds, the gas exchange of slightly soluble gases such as carbon dioxide, dimethyl sulfide, methane, or oxygen across a clean and nearly flat air‐water interface is routinely described using a water‐side mean gas transfer velocity $\overline{k_{L}}$, where overline indicates time or ensemble averaging. The micro‐eddy surface renewal model predicts $\overline{k_{L}}$ = α o Sc -1/2 ($v\bar{ϵ}$) 1/4 , where Sc is the molecular Schmidt number, $v$ is the water kinematic viscosity, and $\bar{ϵ}$ is the waterside mean turbulent kinetic energy dissipation rate at or near the interface. While α o = 0.39 - 0.46 has been reported across a number of data sets, others report large scatter or variability around this value range. It is shown here that this scatter can be partly explained by high temporal variability in instantaneous ϵ around $\bar{ϵ}$, a mechanism that was not previously considered. As the coefficient of variation (CV e ) in ϵ increases, α o must be adjusted by a multiplier (1 = CV e 2 ) -3/32 that was derived from a log‐normal model for the probability density function of ϵ. Reported variations in α o with a macro‐scale Reynolds number can also be partly attributed to intermittency effects in ϵ. Such intermittency is characterized by the long‐range (i.e., power‐law decay) spatial auto‐correlation function of ϵ. That α o varies with a macro‐scale Reynolds number does not necessarily violate the micro‐eddy model. Instead, it points to a coordination between the macro‐ and micro‐scales arising from the transfer of energy across scales in the energy cascade.

Batchelor scale↗

DOE OSTI2025

The Turbulent Pressure Spectrum Within the Roughness Sublayer of a Subarctic Forest Canopy

The turbulent static pressure spectrum E pp (k x ) as a function of longitudinal wavenumber k x in the roughness sublayer of forested canopies is of interest to a plethora of problems such as pressure transport in the turbulent kinetic energy budget, pressure pumping from snow or forest floor, and coupling between flow within and above canopies. Long term static pressure measurements above a sub-arctic forested canopy for near-neutral conditions during the winter and spring were collected and analyzed for three snow cover conditions: trees and ground covered with snow, trees are snow free but the ground is covered with snow, and snow free cover. In all three cases, it is shown that E pp (k x ) obeys the attached eddy hypothesis at low wavenumbers (1/δ < k x < 1/z) —with E pp (k x ) ∝ u$^4_*$k$^{-1}_x$ and Kolmogorov scaling in the inertial subrange at higher wavenumbers—with E pp (k x ) ∝ ϵ 4/3 k x -7/3 , where u * is the friction velocity at the canopy top, is the mean turbulent kinetic energy dissipation rate, z is the distance from the snow top, and is the boundary layer depth. The implications of these two scaling laws to the normalized root-mean squared pressure C p = o p /u$^2_*$ and its newly proposed logarithmic scaling with normalized wall-normal distance z/δ are discussed for snow covered and snow free vegetation conditions. The work here also shows that k$^{-1}_x$ the in the E pp (k x ) appears more extensive and robust than its longitudinal velocity counterpart.

Aslan, Toprak [Finnish Meteorological Inst., Helsi↗

DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Scaling Exponents of Turbulent Static Pressure Structure Function in the Inertial Subrange

The measured variations in the turbulent static pressure structure function D pp (r) with scale $r$ in the roughness sublayer above a subarctic forest are empirically shown to exhibit exponents that are smaller than $r$ 4/3 predicted for the inertial subrange (ISR). Three hypotheses are offered to explain these deviations. The first is based on conventional intermittency correction to the averaged turbulent kinetic energy dissipation rate, the second is based on shearing introducing deviations from locally isotropic state that must be sensed by both velocity and pressure structure functions, and the third is based on large and inertial scale pressure interactions that persist at values of $r$ within the resolvable ISR. The third hypothesis is shown to yield superior results, which allows a new formulation for D pp (r) to be derived that accommodates such finite interactions.

structure function↗

DOE OSTI2023

Hybrid magnonics in hybrid perovskite antiferromagnets

Abstract Hybrid magnonic systems are a newcomer for pursuing coherent information processing owing to their rich quantum engineering functionalities. One prototypical example is hybrid magnonics in antiferromagnets with an easy-plane anisotropy that resembles a quantum-mechanically mixed two-level spin system through the coupling of acoustic and optical magnons. Generally, the coupling between these orthogonal modes is forbidden due to their opposite parity. Here we show that the Dzyaloshinskii–Moriya-Interaction (DMI), a chiral antisymmetric interaction that occurs in magnetic systems with low symmetry, can lift this restriction. We report that layered hybrid perovskite antiferromagnets with an interlayer DMI can lead to a strong intrinsic magnon-magnon coupling strength up to 0.24 GHz, which is four times greater than the dissipation rates of the acoustic/optical modes. Our work shows that the DMI in these hybrid antiferromagnets holds promise for leveraging magnon-magnon coupling by harnessing symmetry breaking in a highly tunable, solution-processable layered magnetic platform.

14 SOLAR ENERGY↗

DOE OSTI2020

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth’s bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (d i ), and the other has a wavelength less than 1 d i . Here, electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Mechanism for sequestering magnetic energy at large scales in shear-flow turbulence

Straining of magnetic fields by large-scale shear flow, which is generally assumed to lead to intensification and generation of small scales, is reexamined in light of the persistent observation of large-scale magnetic fields in astrophysics. It is shown that, in magnetohydrodynamic turbulence, unstable shear flows have the unexpected effect of sequestering magnetic energy at large scales due to counteracting straining motion of nonlinearly excited large-scale stable eigenmodes. This effect is quantified via dissipation rates, energy transfer rates, and visualizations of magnetic field evolution by artificially removing the stable modes. In conclusion, these analyses show that predictions based upon physics of the linear instability alone miss substantial dynamics, including those of magnetic fluctuations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Three-dimensional shear-flow instability saturation via stable modes

Turbulence in three dimensions (3D) supports vortex stretching that has long been known to accomplish energy transfer to small scales. Moreover, net energy transfer from large-scale, forced, unstable flow-gradients to smaller scales is achieved by gradient-flattening instability. Despite such enforcement of energy transfer to small scales, it is shown here that the shear-flow-instability-supplied 3D-fluctuation energy is largely inverse-transferred from the fluctuation to the mean-flow gradient, and such inverse transfer is more efficient for turbulent fluctuations in 3D than in two dimensions (2D). The transfer is due to linearly stable eigenmodes that are excited nonlinearly. The stable modes, thus, reduce both the nonlinear energy cascade to small scales and the viscous dissipation rate. The vortex-tube stretching is also suppressed. Up-gradient momentum transport by the stable modes counters the instability-driven down-gradient transport, which also is more effective in 3D than in 2D (≈70% vs ≈50%). From unstable modes, these stable modes nonlinearly receive energy via zero-frequency fluctuations that vary only in the direction orthogonal to the plane of 2D shear flow. The more widely occurring 3D turbulence is thus inherently different from the commonly studied 2D turbulence, despite both saturating via stable modes.

Mechanics↗

DOE OSTI2024

Adaptive nonequilibrium design of actin-based metamaterials: Fundamental and practical limits of control

The adaptive and surprising emergent properties of biological materials self-assembled in far-from-equilibrium environments serve as an inspiration for efforts to design nanomaterials. In particular, controlling the conditions of self-assembly can modulate material properties, but there is no systematic understanding of either how to parameterize external control or how controllable a given material can be. Here, we demonstrate that branched actin networks can be encoded with metamaterial properties by dynamically controlling the applied force under which they grow and that the protocols can be selected using multi-task reinforcement learning. These actin networks have tunable responses over a large dynamic range depending on the chosen external protocol, providing a pathway to encoding “memory” within these structures. Interestingly, we obtain a bound that relates the dissipation rate and the rate of “encoding” that gives insight into the constraints on control—both physical and information theoretical. Taken together, these results emphasize the utility and necessity of nonequilibrium control for designing self-assembled nanostructures.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Quasioptical propagation and absorption of electron cyclotron waves: simulations and experiment

Electron-cyclotron resonance heating is an essential tool for manipulating fusion plasmas and has been diligently studied for decades through numerical simulations and experiments. Here, we report the first comparison of simulations using the quasioptical code PARADE (PAraxial RAy DEscription) with experiments using a target-plate system on the Large Helical Device. We find that the numerical and experimental results agree qualitatively, which makes PARADE a promising code for modeling of electron-cyclotron wave beams in future experiments. The signature advantage of PARADE is its ability to capture both mode conversion and inhomogeneity of the dissipation rate across the beam cross section. Accounting for these effects, which can noticeably affect the power deposition, significantly improves the agreement between numerical results and experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗