Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock-driven turbulence”

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

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

14 records
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↗
NASA NTRS2022

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.

magnetic reconnection↗
DOE OSTI2022

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal (rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multiphase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(\sigma _{\rho /\langle \rho \rangle }^{2\Gamma }-1)^{1/2}(1-\sigma _{\rho /\langle \rho \rangle }^{-2})^{-1/2}\mathcal {M}^{-1}\Gamma ^{-1/2}$, with the density dispersion σρ/<ρ>, the turbulent Mach number $\mathcal {M}$, and the polytropic exponent Γ. Purely solenoidal and purely compressive driving correspond to b ~ 1/3 and b ~ 1, respectively. As a result, using simulations in which a shock is driven into a multiphase medium with structures of different sizes and Γ < 1, we find b ~ 1 for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Small-scale fluctuating magnetic fields of order nG are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel 3D magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use 384 × 512 × 384 tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay (E kin ∝ t -1.15 ) without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short (~0.3 of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian–Lagrangian bridge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Kinetic characterization of 3D magnetic reconnection: A transformative step (Final Report)

We investigated magnetic reconnection in the Earth’s magnetopause, magnetotail, and the bow shock, by means of fully kinetic simulations and analysis of space data obtained by NASA’s Magnetospheric Multiscale mission. We elucidated signatures of particle energization in particle distribution functions in the vicinity of reconnection X-line, properties of reconnection in shock-driven turbulence, and wave excitations and nonlinear structures. We advanced the understanding of magnetic reconnection from the viewpoint of particle kinetics in the diffusion region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Rayleigh–Taylor Instability With Varying Periods of Zero Acceleration

We present our findings from a numerical investigation of the acceleration-driven Rayleigh–Taylor Instability, modulated by varying periods without an applied acceleration field. It is well known from studies on shock-driven Richtmyer–Meshkov instability that mixing without external forcing grows with a scaling exponent as ≈ t 0.20-0.28 When the Rayleigh–Taylor Instability is subjected to varying periods of “zero” acceleration, the structural changes to the mixing layer remain remarkably small. After the acceleration is re-applied, the mixing layer quickly resumes the profile of development it would have had if there had been no intermission. As a result, this behavior contrasts in particular with the strong sensitivity that is found to other variable acceleration profiles examined previously in the literature.

42 ENGINEERING↗
DOE OSTI2020

Scalar Power Spectra and Scalar Structure Function Evolution in the Richtmyer–Meshkov Instability Upon Reshock

The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗