Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Strong reconnection electric fields in shock-driven turbulence

Turbulent magnetic reconnection in a quasi-parallel shock under parameters relevant to the Earth's bow shock is investigated by means of a two-dimensional particle-in-cell simulation. The addressed aspects include the reconnection electric field, the reconnection rate, and the electron and the ion outflow speeds. In the shock transition region, many current sheets are generated in shock-driven turbulence, and electron-only reconnection and reconnection where both ions and electrons are involved can occur in those current sheets. The electron outflow speed in electron-only reconnection shows a positive correlation with the theoretical speed, which is close to the local electron Alfvén speed, and a strong convection electric field is generated by the large electron outflow. As a result, the reconnection electric field becomes much larger than those in the standard magnetopause or magnetotail reconnection. In shock-driven reconnection that involves ion dynamics, both electron outflows and ion outflows can reach of the order of 10 times the Alfvén speed in the X-line rest frame, leading to a reconnection electric field the same order as that in electron-only reconnection. An electron-only reconnection event observed by the magnetospheric multiscale mission downstream of a quasi-parallel shock is qualitatively similar to those in the simulation and shows that the outflow speed reaches approximately half the local electron Alfvén speed, supporting the simulation prediction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

New insights on wind turbine wakes from large‐eddy simulation: Wake contraction, dual nature, and temperature effects

Abstract Large‐eddy simulation (LES) has been adopted to study wind turbine wakes because it can capture fine‐scale details of turbulent wind flows and interactions with wind turbines. Here, we use the LES version of the Weather Research and Forecasting (WRF) model with an actuator disk model to gain insights on several wake effects that have been traditionally difficult to measure. The first finding is that the wake has a “dual nature,” meaning that the wind speed deficit behaves differently from the added turbulent kinetic energy (TKE) and the two are not co‐located in space. For example, the wind speed deficit peaks at hub height and reaches the ground within 8D (D is the rotor diameter), but added TKE peaks near the rotor tip and generally remains aloft. Second, temperature changes near the ground are driven by the added TKE in the rotor area and by atmospheric stability. The combination of these two factors determines the sign and intensity of the vertical heat flux divergence below the rotor, with convergence and warming associated with stable conditions and weak divergence and modest cooling with unstable conditions. Third, wakes do not expand indefinitely, as suggested by similarity theory applied to the wind speed deficit, but eventually stop expanding and actually contract, at different rates depending on atmospheric stability. The implication of these findings is that, in order to study wakes, it is not sufficient to focus on wind speed deficit alone, because TKE is also important and yet behaves differently from the wind speed deficit.

17 WIND ENERGY↗
DOE OSTI2025

Autonomous sensor suite for evaluating fish-turbine interactions and environmental impacts in marine renewable energy and hydropower

Marine renewable energy (MRE) harnesses ocean-based resources such as waves, tides, currents, and thermal or salinity gradients for sustainable power generation. It has the potential to complement existing renewable resources, support remote communities, and contribute to decarbonization efforts. However, understanding the hydrodynamic forces created by MRE devices and their impacts on marine life is critical for responsible deployment. Here, to address these concerns, advanced sensor devices, including the Marine Sensor Fish (MSF), Sensor Fish Mini (SF Mini), and Flexible Sensor Fish (FSF), were developed to measure interactions between aquatic organisms and MRE systems. This paper details the design, manufacturing, calibration, and field deployment of these sensor suites, highlighting their ability to capture key physical stressors such as shear forces, pressure changes, and collision impacts. The MSF successfully evaluated turbine interactions at a tidal turbine in the Salish Sea, capturing data on turbulence, collision impact, and pressure gradients. The SF Mini validated hydrodynamic conditions in scaled hydraulic models, supporting computational fluid dynamics simulations. The FSF, with its flexible silicone body, measured species-specific impacts in turbulent environments. This research demonstrates the potential of Sensor Fish technology to advance sustainable marine energy systems by reducing biological impacts and informing environmentally sustainable designs.

Ecological impacts↗
DOE OSTI2020

Ion versus Electron Heating in Compressively Driven Astrophysical Gyrokinetic Turbulence

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio Q i /Q e as a function of the compressive-to-Alfvénic driving power ratio P compr /P AW , of the ratio of ion thermal pressure to magnetic pressure β i , and of the ratio of ion-to-electron background temperatures T i /T e . It is shown that Q i /Q e is an increasing function of P compr /P AW . When the compressive driving is sufficiently large, Q i /Q e approaches ≃P compr /P AW . This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high β i , meaning that the compressive driving is directly connected to the ion-entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low βi and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Turbulent transport and reactions of plant-emitted hydrocarbons in an Amazonian rain forest

The processes governing the temporal and spatial patterns of isoprene and monoterpenes emitted by a rainforest in the central Amazon region of Brazil is investigated using a combination of field experiments and numerical simulations. Specifically, Large Eddy Simulations (LES) were used to resolve emissions of isoprene and monoterpenes, turbulent transport, and air chemistry. The coupled chemistry-transport LES included the effects of isoprene and monoterpenes reactivity due to reactions with hydroxyl radical and ozone. The LES results are used to compute vertically resolved budgets of isoprene and monoterpenes in the rainforest canopy in response to emissions, turbulent transport, surface deposition, and air chemistry. Our results indicated that emission and dispersion dominated the isoprene budget as the gases were transported out of the canopy space. In a region limited by nitrogen oxides (with prevailing nitric oxide levels of < 0.5 parts per billion), the in-canopy chemical destruction removed approximately 10% of locally emitted monoterpenes. Hydroxyl radical production rates from the ozonolysis of monoterpenes amounted to ≈ 2 x 10 6 radicals cm -3 s -1 and had similar magnitude to the light-dependent hydroxyl radical formation. One key conclusion was that the Amazonia rainforest abundantly emitted monoterpenes whose in-canopy ozonolysis yielded hydroxyl radicals in amounts similar to the magnitude of light-dependent formation. Reactions of monoterpenes and isoprene with hydroxyl radical and ozone were necessary for the maintenance of the Amazon rainforest canopy as a photochemically active environment suitable to generate oxidants and secondary organic aerosols.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2024

Formulation, Implementation and Validation of a 1D Boundary Layer Inflow Scheme for the QUIC Modeling System

Recent studies have highlighted the importance of accurate meteorological conditions for urban transport and dispersion calculations. In this work, we present a novel scheme to compute the meteorological input in the Quick Urban & Industrial Complex () diagnostic urban wind solver to improve the characterization of upstream wind veer and shear in the Atmospheric Boundary Layer (ABL). The new formulation is based on a coupled set of Ordinary Differential Equations (ODEs) derived from the Reynolds Averaged Navier–Stokes (RANS) equations, and is fast to compute. Building upon recent progress in modeling the idealized ABL, we include effects from surface roughness, turbulent stress, Coriolis force, buoyancy and baroclinicity. We verify the performance of the new scheme with canonical Large Eddy Simulation (LES) tests with the GPU-accelerated FastEddy"Equation missing" solver in neutral, stable, unstable and baroclinic conditions with different surface roughness. Furthermore, we evaluate QUIC calculations with and without the new inflow scheme with real data from the Urban Threat Dispersion (UTD) field experiment, which includes Lidar-based wind measurements as well as concentration observations from multiple outdoor releases of a non-reactive tracer in downtown New York City. Compared to previous inflow capabilities that were limited to a constant wind direction with height, we show that the new scheme can model wind veer in the ABL and enhance the prediction of the surface cross-isobaric angle, improving evaluation statistics of simulated concentrations paired in time and space with UTD measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Distributions of Particles Accelerated by Strong Alfvénic Turbulence

This work presents a model for generating nonthermal power-law tails of particles’ energy probability density functions in turbulent collisionless plasmas, applicable to both nonrelativistic and relativistic scenarios. We propose that strong Alfvénic turbulence energizes plasma particles through curvature acceleration, particularly for particles with Larmor radii comparable to the scales of turbulence. When the energy density of the energized particles increases, the efficiency of the energy exchange process diminishes. As a result, the acceleration process saturates, leading to power-law distributions of particle momentum and energy. In the nonrelativistic case, the momentum probability density function scales as f(p)dp ∝ p −3 dp, while in the ultrarelativistic case, the energy probability density function scales as f(γ)dγ ∝ γ −3 dγ, where γ is the Lorentz factor. This model provides a unified framework for understanding particle acceleration in both energy regimes, complementing existing analytical approaches. The predicted scalings are consistent with available observations of energetic ion distributions in the heliosphere and with the findings from numerical simulations of ultrarelativistic particle acceleration in magnetically dominated plasma turbulence.

Alfven waves↗
DOE OSTI2023

Evolution and hot electron generation of laser–plasma instabilities in direct-drive inertial confinement fusion

A series of 2D in-plane plane wave particle-in-cell simulations find distinctive paths of laser-plasma instability evolution in OMEGA-scale implosions, depending on the initial electron temperature. At low temperatures, two-plasmon decay (TPD) dominates in both initial growth and the steady state. At high temperatures, the initial dominant modes switch to stimulated Raman scattering, but TPD still dominates a steady state characterized by cavitation and Langmuir turbulence. A hot electron scaling is also obtained from the simulations that, when combined with laser/plasma conditions from hydro simulations, can predict hot electron generation in implosions that do not employ smoothing-by-spectral-dispersion (SSD). It also shows that under the same laser/plasma conditions, SSD can reduce hot electron generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Are Marine Low Cloud Droplet Concentrations Buffered by Entrained Aitken‐Mode Aerosol (Final Technical Report)

During the summertime, the high-latitude oceans come to life with green phytoplankton, which gain their energy from sunlight and are food for sea creatures small and large. When the phytoplankton are eaten or die, sulfur-rich gases are released under the ocean surface and mix into the air. Observations suggest that, over the Southern Ocean, frequent storms lift this air high into the atmosphere while raining out particulates like salt. As a result, the sulfur-rich air then spawns high concentrations of small ‘Aitken-mode’ aerosol particles. We hypothesize that these particles work their way down into the marine boundary layer, where they can replenish the supply of cloud-condensation nuclei scavenged by frequent precipitation. Further, this process maintains high concentrations of liquid cloud droplets in austral summer, promoting more sunlight to be reflected to space. We call this ‘Aitken buffering’. The primary objective of this project has been to document and test what role Aitken-mode aerosols play in clouds over the Southern Ocean and elsewhere. This effort has included three main components: 1) developing a computer model that realistically simulates the aerosol processes and the small-scale turbulent air motions that move aerosols around and create the clouds, 2) using that model to interpret and extend these observations for process understanding, by allowing different factors that contribute to the aerosol budget, such as surface wind speed, precipitation, surface gas exchange, etc. to be separated, and 3) studying Aitken-mode aerosol and its variability with a focus over the Southern Ocean and Antarctica, using data from Atmospheric Radiation Measurement (ARM) sites and other available observations. Initial computer studies in more idealized conditions found that elevated concentrations of Aitken-mode aerosols above the clouds could help prevent the breakup of those clouds by acting as cloud-condensation nuclei after they were entrained into the cloudy boundary layer. Simulations of a day during the ACE-ENA field campaign showed that Aitken-mode aerosols could also prevent cloud breakup under more realistic conditions. A new method that extracts information about Aitken-mode aerosols from measurements of aerosols onto which cloud droplets can form finds that Aitken-mode aerosols do vary seasonally over the Southern Ocean, with a peak in summertime, as described above. Other work during this project has focused on understanding how patterns of water vapor, clouds and precipitation are coupled within low-lying clouds over the oceans, and also on how cloud droplets cluster within clouds and how the distribution of cloud droplet sizes change as dry air is mixed into clouds, with the latter studies also using observations from ACE-ENA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2020

Collaborative Research: The Evolution of Magnetic Skeletons During 3D Reconnection. Final report

We investigated how the plasmoid instability affects the evolution of topological features of a magnetic field and the process of magnetic reconnection with three sub-projects. (i) We extensively studied the criterion for plasmoid-mediated current sheet disruption and the onset of fast reconnection. Scaling relations with respect to plasma parameters have been established, using both direct numerical simulations and analytic theories, for the disruption time, the current sheet width, the instability linear growth rate, and the dominant mode wavenumber. (ii) We applied the theories of plasmoid-mediated current sheet disruption to magnetohydrodynamic turbulence and predicted a steepening of the turbulence energy spectrum. This steepening has been confirmed by direct numerical simulations. (iii) We performed numerical simulations of ultraviolet burst events in the solar atmosphere with a spine-dome magnetic topology. The plasmoid instability leads to unsteady reconnection, which is consistent with solar observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility

Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.

buoyant mixing↗
DOE OSTI2021

Millimeter-wave beam scattering and induced broadening by plasma turbulence in the TCV tokamak

The scattering of millimeter-wave beams from electron density fluctuations and the associated beam broadening are experimentally demonstrated. Using a dedicated setup, instantaneous deflection and (de-)focusing of the beam due to density blobs on the beam path are shown to agree with full-wave simulations. The detected time-averaged wave power transmitted through the turbulent plasma is reproduced by the radiative-transfer model implemented in the WKBeam code, which predicts a ~50% turbulence-induced broadening of the beam cross-section. Finally, the role of core turbulence for the considered geometry is highlighted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Global Hybrid Simulations of Interaction Between Interplanetary Rotational Discontinuity and Bow Shock/Magnetosphere: Can Ion-Scale Magnetic Reconnection be Driven by Rotational Discontinuity Downstream of Quasi-Parallel Shock?

Ion-scale magnetic reconnection has been observed downstream of the terrestrial quasi-parallel (Q-∥) shock. Whether it is driven by interplanetary discontinuities or turbulent Q-∥ shock, however, is unclear. Using three-dimensional global hybrid simulation, we investigate the generation of magnetic reconnection downstream of the Q-∥ shock, while an interplanetary rotational discontinuity (RD) is launched to the bow shock. Cases with various solar wind Alfvén Mach numbers, M A = 3.0 to 8, and propagation directions of the RD are presented. The propagation direction n is assumed to be in the GSE xz plane and pointing earthward, with n = (-sin(θ 12 /2),0,-cos(θ 12 /2)), where θ 12 is the angle between the upstream (B1) and downstream (B2) magnetic fields across the transmitted RD. It is found that magnetic reconnection occurs inside the RD downstream of the Q-∥ shock, forming flux ropes extending along the dawn-dusk direction about tens of ion inertial lengths. Large-amplitude low-frequency waves originated from the Q-∥ shock lead to the bending and squeezing of the field lines around the RD, which play an important role in triggering reconnection inside the RD. As the RD impacts the dayside magnetopause, magnetopause reconnection takes place between the field lines behind the RD and geomagnetic field lines. Nevertheless, no reconnection is found downstream of the Q-∥ shock itself or outside the RD in the magnetosheath. The existent and structure of reconnection in the magnetosheath are found to strongly depend on the parameters M A and n. Our simulation shows that ion-scale magnetic reconnection is driven by an external driver in the form of the compression of an RD around the bow shock and in the magnetosheath, rather than caused by the turbulent Q-∥ shock alone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Measuring scale-dependent shape anisotropy by coarse-graining: Application to inhomogeneous Rayleigh-Taylor turbulence

Here, we generalize the “filtering spectrum” to probe scales along different directions by spatial coarse-graining. This multidimensional filtering spectrum quantifies the spectral content of flows that are not necessarily homogeneous. From multidimensional spectral information, we propose a simple metric for shape anisotropy at various scales. The method is applied to simulations of two-dimensional (2D) and 3D Rayleigh-Taylor (RT) turbulence, which is inhomogeneous and anisotropic. We show that 3D RT has clear shape anisotropy at large scales with approximately 4:3 vertical to horizontal aspect ratio but tends toward isotropy at small scales as expected. In sharp contrast, we find that RT in 2D simulations, which are still the main modeling framework for many applications, is isotropic at large scales and its shape anisotropy increases at smaller scales where structures tend to be horizontally elongated. While this may be surprising, it is consistent with recent results that large-scale isotropy in 2D RT is due to the generation of a large-scale overturning circulation via an upscale cascade, while small-scale anisotropy is due to the stable stratification resultant from such overturning and the inefficient mixing in 2D.

42 ENGINEERING↗