Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pressure gradient”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Importance of $\delta B_{\|}$ on ETG stability, turbulence, and transport in NSTX

This study employs electron-scale gyrokinetic simulations to investigate the electron temperature gradient (ETG) driven instabilities, turbulence, and transport in the pedestal region of the National Spherical Torus Experiment, comparing non-lithiated (narrow pedestal) and lithiated (wide pedestal) scenarios. Our findings reveal that, in the non-lithiated case, a branch of strongly unstable ETG modes exhibiting finite parallel magnetic field fluctuations ($\delta B_{\parallel} \neq 0$) emerges at the pedestal top and upper density pedestal region. This branch is uncovered only when $\delta B_{\parallel}$ is retained in the simulations and is associated with substantial electrostatic electron heat flux. This region of strong ETG transport corresponds to the only region in the plasma where the pressure gradient is far below the critical gradient for kinetic ballooning modes. We investigated the origin of this finite $\delta B_{\parallel}$ ETG branch by analyzing the gyrokinetic field equations. Nonlinear saturation is also analyzed and contrasted for simulations with and without $\delta B_{\parallel}$. In contrast with the nonlithiated case, ETG modes in the lithiated case produce substantial transport in the steep gradient region, but are negligible at the pedestal top.

ETG↗
DOE OSTI2025

Algorithmic Advancements for High-Order Self-Gravitating Hydrodynamics

Self-gravity plays a key role in the formation and evolution of many astronomical objects. Though gravity is often dominant at large scales, other forces (e.g., gas pressure gradients, radiation, and/or magnetic fields) often compete. It is therefore essential for numerical simulations to evaluate their interplay accurately and robustly. Hanawa & Mullen derived a 4th-order accurate finite volume scheme to solve the equations of self-gravitating hydrodynamics on a uniform Cartesian grid. In this work, we supply improvements to the algorithm that (1) mitigate spurious gravitational circulation and (2) greatly simplify the evaluation of the high order corrections. The proposed algorithm provides the gravitational acceleration (ρg) and the gravitational energy release (ρv · g) as source terms for the hydrodynamic equations, all while preserving conservation of linear momentum. Spurious heating and/or cooling associated with truncation error in the numerical evaluation of the gravitational energy release decreases in proportion to the fourth power of the cell width. We demonstrate fourth order convergence on smooth problems (e.g., 3D inclined sound wave propagation and 3D equilibria). An application test tracks the spherical collapse of a polytrope by an imposed, sudden decrease of the central gas pressure; a bounce and second collapse (associated with a spherical accretion shock) are robustly captured by the high order algorithm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
DOE OSTI2025

Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids

Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.

42 ENGINEERING↗
DOE OSTI2021

On the interplay between magnetic field and anisotropy in holographic QCD

We investigate the combined effects of anisotropy and a magnetic field in strongly interacting gauge theories by the gauge/gravity correspondence. Our main motivation is the quark-gluon plasma produced in off-central heavy-ion collisions which exhibits large anisotropy in pressure gradients as well as large external magnetic fields. We explore two different configurations, with the anisotropy either parallel or perpendicular to the magnetic field, focusing on the competition and interplay between the two. A detailed study of the RG flow in the ground state reveals a rich structure where depending on which of the two, anisotropy or magnetic field, is stronger, intermediate geometries with approximate AdS 4 × $\mathbb{R}$ and AdS 3 × $\mathbb{R}$ 2 factors arise. This competition is also manifest in the phase structure at finite temperature, specifically in the dependence of the chiral transition temperature on anisotropy and magnetic field, from which we infer the presence of inverse magnetic and anisotropic catalyses of the chiral condensate. Finally, we consider other salient observables in the theory, including the quark-antiquark potential, shear viscosity, entanglement entropy and the butterfly velocity. We demonstrate that they serve as good probes of the theory, in particular, distinguishing between the effects of the magnetic field and anisotropy in the ground and plasma states. We also find that the butterfly velocity, which codifies how fast information propagates in the plasma, exhibits a rich structure as a function of temperature, anisotropy and magnetic field, exceeding the conformal value in certain regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Coupled Thermo-Hydro-Mechanical Processes in Fractured Rocks: Some Past Scientific Highlights and Future Research Directions

Abstract Coupled thermo-hydro-mechanical (THM) processes in fractured rocks have been a topic of intense scientific research for more than 30 years. The present paper takes a look into the past and highlights some scientific advances which are of an unusual “out-of-the-box” nature, and then looks forward and discusses possible directions of future research in this interesting field of study. Concerning future research directions, we see a trend from a focus on coupled THM processes in single fractures or a few interacting fractures, to the study of coupled THM behavior in complex fracture network systems where the fractures act collectively giving rise to local stress concentration points and points of large pressure gradients. Three examples of future research directions are presented. First is an effort towards identifying characterizing parameters of a fracture network that play a direct controlling role in major coupled THM phenomena (such as induced seismicity and flow channeling), rather than parameters of stochastic distributions of fractures in the network. The second example of research direction is accounting for the heterogeneity and hierarchy of fractures in a fault or fracture zone which has been associated with major THM events in a number of geo-energy projects. The third example is at the opposite end of the first; here it is recognized that in some cases, the coupled THM processes in fractured rocks may be controlled dominantly by only a few key bridges. Identification, characterization, and evaluation of these key bridges should be one of the important research directions in the coming days.

58 GEOSCIENCES↗
DOE OSTI2026

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Condensation heat transfer and pressure drop of low-global warming potential refrigerants in smooth aluminum tubes

In this work we report condensation heat transfer coefficient and pressure drop data for R-454B, R-32, and R-454C. These refrigerants are low Global Warming Potential (GWP) replacements for R-410A. We obtained these data in a smooth 7.2 mm outer diameter (5.94 mm inner diameter) aluminum tube. In these experiments, the key variables were condensation temperature (35 °C ≤ T cond ≤ 50 °C) and refrigerant mass flux (151.1 ≤ G ≤ 347.3 kg m –2 s –1 ). In general, for similar nominal operating conditions, we measured the highest heat transfer coefficients for R-32 and the lowest for R-454C. Among the correlations we used to predict the heat transfer data for all three refrigerants, the one developed by Thome et al. (2003) predicted the data with the greatest accuracy. Deviations between experimental data and predictions of this correlation were 21.7% for R-454B (with SBG correction), 22.0% for R-32, and 18.8% for R-454C (with SBG correction). Experimental pressure gradient data were predicted by the Friedel correlation Friedel (1979) with deviations of 13.8% for R-454B, 9.8% for R-32, and 10.2% for R-454C. In this work, we also compare the predictions of the Cavallini et al. correlation Cavallini et al. (2006) for R-454B and R-32 in smooth copper versus smooth aluminum tubes. To improve this correlation’s predictive capability for aluminum tubes, we have proposed some modifications. After incorporating these modifications, the deviations between the correlation’s predictions and our data were 10% for R-454B (with SBG correction), 8.1% for R-32, and 5.2% for R-454C (with SBG correction). Furthermore, these data and findings will assist the HVAC and refrigeration industry to implement low-GWP refrigerants in condensers using smooth aluminum tubes.

Aluminum tube↗
DOE OSTI2021

Carbon aerogels with integrated engineered macroporous architectures for improved mass transport

Carbon aerogels (CAs) combine unique properties including ultra-high surface area, high electrical conductivity, corrosion resistance, and robust mechanical properties making them ideal materials for electrochemical applications. Traditional CA synthesis results in isotropic, random nanoporous networks that work well for applications relying on diffusional mass transport. However, many applications would benefit from integration of engineered macroporous network structures that enable directed pressure-gradient-driven mass transport. In this study, we report on using 3D-printed sacrificial polymeric templates to generate templated CAs (t-CAs) with integrated engineered nonrandom macroporous network structures. Specifically, we used projection micro-stereo-lithography (PμSL) and two-photon polymerization direct laser writing (2PP-DLW) to fabricate millimeter-to-centimeter-sized 3D sacrificial polymeric templates with features ranging from tens of microns (PμSL) to 100s of nanometers (2PP-DLW). T-CAs were fabricated by infiltrating the templates with resorcinol-formaldehyde (RF) precursor solution, followed by carbonization at 1050 °C to simultaneously convert the RF gel to a CA and decompose the 3D-printed template, leaving an embedded templated macroporous network structure behind. X-ray computer tomography confirms integration of the macroporous architecture defined by the template. The templated macroporous architecture improves mass transport in t-CAs compared to traditional bulk CA as demonstrated by more uniform activation and their response in electrochemical cyclic voltammetry and galvanostatic charge-discharge tests.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of initial water saturation on the performance of fracturing fluids with and without polyallylamine under simulated EGS conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO 2 - reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO 2 , can be used as a potentially less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO 2 fracturing fluids outperformed other fluids such as water, CO 2 , and CO 2 /water in generating large fractures with less fluid consumed. The objective of this work is to investigate the effect of initial water saturation on the performance of StimuFrac by conducting hydraulic fracturing tests with ½ foot cubic rock samples held under representative EGS stress/temperature conditions and by using cyclic injection strategies (under constant injection rate). The resulting fracture hydraulic conductivities, breakdown pressures, and volumes of fluids required are compared. Methods/Procedures/Process: To simulate geothermal reservoir conditions, in all tests, the rock sample was held under triaxial confinement and at 200 °C, and different volumes of water were initially injected into the rock sample before any fracturing processes were initiated. For the single-cycle PAA (or water) alternating CO 2 (PAG or WAG) injection fracturing experiments, one complete cycle consisted of two steps: (1) injecting a PAA slug (or water slug) followed by (2) injecting CO 2 to initiate and propagate the fracture. For experiments involving multiple injection cycles, the CO 2 injection pressure is increased until it peaks and begins to decline (indicating fracture initiation at this moment), and then continued being injected for another 30 seconds to propagate the fracture. Then, these two-step cycles [injection of PAA (or water) followed by CO 2 injection (up to 2-4 mL/min)] are repeated. Applications/Significance/Novelty: The results of this study suggest that water saturation significantly affects the fracturing fluid transmission into the rock pore space, thus affecting the fracture initiation and propagation. In this study, fracturing tests via a single injection cycle or multiple injection cycles were performed. Splitting the rock samples in half after testing reveals that fracture propagation is significantly limited under high water saturation conditions (three-day initial water injection) compared to stimulation experiments performed in hot dry rock. The fractures propagate less than 1/3 of the distance from the wellbore to the outer rock surface, and in some cases, no fracture is generated. This may be caused by leak-off dominating the fracturing process and the fluid injection rate is insufficient to overcome leak-off, even under high injection rate conditions. Additionally, CO 2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate sufficiently high-pressure gradients away from the near-wellbore region. Under low water saturation conditions (dry rock or after 1-day initial water injection), PAA/CO 2 consistently generated significantly larger fractures compared with the other fluids. CO 2 generated large fractures only in the hot dry rock and only when using high injection rates, though data variability is high.

58 GEOSCIENCES↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Integration of high confinement, high poloidal beta plasma with dual radiated power and detachment controls for divertor protection and ELM suppression

Divertor detachment without serious core confinement quality loss in DIII-D’s high poloidal β scenario has been combined with impurity-induced ELM mitigation without disruption. Use of Ne previously granted access to a detached, non-ELMing regime that retained high confinement quality due to stimulation of Internal Transport Barrier (ITB) growth, but suffered from on-going core fuel dilution and high disruptivity. Excess Ne accumulation in the core plasma has now been avoided by feeding back core radiated power ($P_{rad}$) measurements to control Ne seeding, rather than using attachment fraction ($A_{frac}$) control with Ne; this also reduces disruptivity. At the same time, N 2 seeding is used in a feedback loop with $A_{frac}$ measurements, which previously posed low disruption risk. In this way, the effect of Ne in the core is managed while avoiding excess seeding, and N 2 acts to correct for any excess heat exhaust that might interfere with detachment. The average Ne flow rate was 38% of what was used in pure Ne $A_{frac}$ control, plus average N 2 flow that was 43% of pure N 2 $A_{frac}$ control, even while meeting an even deeper detachment target. Meanwhile, a steep pressure gradient in the core plasma reduces sensitivity to impurity-induced degradation of the pedestal and reasonable confinement quality was maintained despite operational challenges that blocked formation of an ITB in these experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Advances and challenges in CO2 foam technologies for enhanced oil recovery in carbonate reservoirs

Utilization of surfactants for generating foam has the potential for fluid-mobility control during CO2 flooding processes, leading to improved oil recovery. This manuscript is a review of the most important aspects for the design of CO2-foams with mobility control in carbonate reservoirs, including recent advances in novel surfactants and analytical techniques for analysis of surfactant adsorption on carbonate minerals and their thermal stability. Several key parameters and properties regarding foam transport in porous media are reviewed such as, the minimum pressure gradient required for foam generation, the effects of the partition coefficient of a CO2-soluble surfactant on foam transport in rock cores, visualization of foam flow by either microfluidic or PET/CT imaging processes in heterogeneous or fractured systems and the effect of oil on foam transport under reservoir conditions. Also are discussed the most recent advances of two major foam modeling methods, including the semi-empirical STARS foam model and the population balance model, as well as CO2 foam pilot tests and the factors governing their successes.

Jian, Guoqing↗
DOE OSTI2023

Decoding the flow evolution in Au+Au reactions at 1.23 A GeV using hadron flow correlations and dileptons

We investigate the development of the directed, v 1 , and elliptic flow, v 2 , in heavy ion collisions in mid-central Au+Au reactions at E lab = 1.23A GeV. We demonstrate that the elliptic flow of hot and dense matter is initially positive (v 2 >0) due to the early pressure gradient. This positive v 2 transfers its momentum to the spectators, which leads to the creation of the directed flow v 1 . In turn, the spectator shadowing of the in-plane expansion leads to a preferred decoupling of hadrons in the out-of-plane direction and results in a negative v 2 for the observable final state hadrons. We propose a measurement of v 1 −v 2 flow correlations and of the elliptic flow of dileptons as methods to pin down this evolution pattern. The elliptic flow of the dileptons allows then to determine the early-state EoS more precisely, because it avoids the strong modifications of the momentum distribution due to shadowing seen in the protons. This opens the unique opportunity for the HADES and CBM collaborations to measure the Equation-of-State directly at 2-3 times nuclear saturation density.

Reichert, Tom↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗