Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Recovery of a high-pressure phase formed under laser-driven compression

The recovery of metastable structures formed at high pressure has been a long-standing goal in the field of condensed matter physics. While laser-driven compression has been used as a method to generate novel structures at high pressure, to date no high-pressure phases have been quenched to ambient conditions. Here we demonstrate, using in situ x-ray diffraction and recovery methods, the successful quench of a high-pressure phase which was formed under laser-driven shock compression. We show that tailoring the pressure release path from a shock-compressed state to eliminate sample spall, and therefore excess heating, increases the recovery yield of the high-pressure ω phase of zirconium from 0% to 48%. Our results have important implications for the quenchability of novel phases of matter demonstrated to occur at extreme pressures using nanosecond laser-driven compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Core-collapse supernova neutrino emission and detection informed by state-of-the-art three-dimensional numerical models

ABSTRACT Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, namely Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral standing accretion shock instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate (${\gtrsim} 50 {{\ \rm per\ cent}}$), but the time-integrated signal has a relatively modest asymmetry (${\lesssim} 20 {{\ \rm per\ cent}}$). Both features are associated with the lepton-number emission self-sustained asymmetry and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavours of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of ∼20 per cent if the distance to the CCSN is ≲6 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gravitational-wave signal of a core-collapse supernova explosion of a 15M ⊙ star

In this paper we report on the gravitational-wave signal computed in the context of a three-dimensional simulation of a core-collapse supernova explosion of a 15 M ⊙ star. The simulation was performed with our neutrino hydrodynamics code chimera . We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (less than 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the standing accretion shock instability (SASI), and high-frequency emission (greater than 600 Hz) emanates from the proto–neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto–neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical–acoustic origin, to it. We compare the results of our three-dimensional model analysis with those obtained from the model’s two-dimensional counterpart and find a significant reduction in the strain amplitudes in the former case, as well as significant reductions in all related quantities. Our dominant proto–neutron star gravitational wave emission is not well characterized by emission from surface g modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto–neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. This simple analysis indicates that the spectrum of gravitational wave emission between approximately 20 Hz and approximately 1 kHz, stemming from neutrino-driven convection, the SASI, accretion onto the proto–neutron star, and proto–neutron star convection, will be accessible for a Galactic event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

R&D Program for HEP High-Power Targets at Fermilab

A high-power target system is a key beam element to complete future High Energy Physics (HEP) experiments. In the recent past, major accelerator facilities have been limited in beam power not by their accelerators, but by the beam intercepting device survivability. The target must then endure high power pulsed beam, leading to high cycle thermal stresses/pressures and thermal shocks. The increased beam power will also create significant challenges such as corrosion and radiation damage that can cause harmful effects on the material and degrade their mechanical and thermal properties during irradiation. This can eventually lead to the failure of the material and drastically reduce the lifetime of targets and beam intercepting devices. In order to operate reliable beam-intercepting devices in the framework of energy and intensity increase projects of the future, it is essential to develop a strong R&D program and have synergy with various expertise. After presenting the high power targetry challenges facing next generation multi-MW accelerators, we will give an overview of Fermilab’s R&D program in support of High Power Targetry development. The RaDIATE collaboration (Radiation Damage In Accelerator Target Environment), managed by Fermilab, also draws on existing expertise in related fields to execute a coordinated strategy for high power targetry R&D between the 14 international member institutions.

Pellemoine, Frederique [Fermilab]↗
DOE OSTI2024

Creating and studying a scaled interplanetary coronal mass ejection

The Sun, being an active star, undergoes eruptions of magnetized plasma that reach the Earth and cause the aurorae near the poles. These eruptions, called coronal mass ejections (CMEs), send plasma and magnetic fields out into space. CMEs that reach planetary orbits are called interplanetary coronal mass ejections (ICMEs) and are a source of geomagnetic storms, which can cause major damage to our modern electrical systems with limited warning. To study ICME propagation, we devised a scaled experiment using the Big Red Ball (BRB) plasma containment device at the Wisconsin Plasma Physics Laboratory. These experiments inject a compact torus of plasma as an ICME through an ambient plasma inside the BRB, which acts as the interplanetary medium. Magnetic and temperature probes provide three-dimensional magnetic field information in time and space, as well as temperature and density as a function of time. Using this information, we can identify features in the compact torus that are consistent with those in real ICMEs. We also identify the shock, sheath, and ejecta similar to the structure of an ICME event. This experiment acts as a first step to providing information that can inform predictive models, which can give us time to shield our satellites and large electrical systems in the event that a powerful ICME were to strike.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗
DOE OSTI2025

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗
DOE OSTI2021

Model of ramp compression of diamond from ab initio simulations

Ramp compression experiments characterize high-pressure states of matter at temperatures well below those present in shock compression. However, because temperature is typically not directly measured during ramp compression, it is uncertain how much heating occurs under these shock-free conditions. Here, we performed a series of ab initio simulations on carbon in order to match the density-stress measurements of Smith et al. [Smith et al., Nature (London) 511, 330 (2014)]. We considered isotropically as well as uniaxially compressed solid carbon in the diamond and BC8 phases, with and without defects, as well as liquid carbon. Our idealized model ascribes heating during ramp compression to an initially uniaxially compressed cell transforming isochorically into an isotropically (hydrostatic equivalent) compressed state having lower internal energy, hence higher temperature so as to conserve energy. Multiple such heating events can occur during a single ramp experiment, leading to higher temperatures than with isentropic compression. Comparison with experiments shows that heating alone does not explain the equation of state measurements on diamond, instead implying that a significant uniaxial stress component remains present at high compression. The temperature predictions of our ramp compression model remain to be verified with future laboratory measurements.

58 GEOSCIENCES↗
DOE OSTI2020

Microbial electrolysis cell recovery after inducing operational failure conditions

Microbial Electrolysis Cells (MECs) are often documented for their ability to produce hydrogen through new sources and new configurations. However, few studies attempt to document the recovery of the bioanode after being exposed to operational failure conditions. This study attempted to compare the behavior and recovery of MECs after being exposed to four operational failure conditions: acidification, aeration, osmotic shock, and voltage reversal. For each failure condition, three time points were tested: before failure, immediately after failure, and after recovery efforts. Of the modes tested, acidification caused the largest loss in performance and operational efficiency, while osmotic shock caused the second largest loss in performance. Aeration and voltage reversal caused negligible losses in performance immediately after failure conditions and after the recovery period. Aeration and voltage reversal had a minimal effect on the removal of individual compounds. However, acetate accumulated after acidification, and propionic acid accumulated after osmotic shock. Furthermore these findings could be useful for determining if failure will require more significant repair efforts in commercially deployed devices. While acidification created permanent losses to performance, if device failures like aeration, osmotic shock, and voltage reversal are caught early (less than 15 min), full recovery without bioanode replacement is likely.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Phase transformation in boron under shock compression

Using first-principles molecular dynamics, we calculated the equation of state and shock Hugoniot of various boron phases. Here, we find a large mismatch between Hugoniots based on existing knowledge of the equilibrium phase diagram and those measured by shock experiments, which could be reconciled if the α-B 12 /β phases are significantly over-pressurized in boron under shock compression. Our results also indicate that there exist an anomaly and negative Clapeyron slope along the melting curve of boron at 100 GPa and 1500–3000 K. These results enable an in-depth understanding of matter under shock compression, in particular the significance of compression-rate dependence of phase transitions and kinetic effects in experimental measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Testing the thermal Sunyaev-Zel’dovich power spectrum of a halo model using hydrodynamical simulations

Statistical properties of large-scale cosmological structures serve as powerful tools for constraining the cosmological properties of our Universe. Tracing the gas pressure, the thermal Sunyaev-Zel’dovich (tSZ) effect is a biased probe of mass distribution and, hence, can be used to test the physics of feedback or cosmological models. Therefore, it is crucial to develop robust modelling of hot gas pressure for applications to tSZ surveys. Since gas collapses into bound structures, it is expected that most of the tSZ signal is within halos produced by cosmic accretion shocks. Hence, simple empirical halo models can be used to predict the tSZ power spectra. In this study, we employed the HMx halo model to compare the tSZ power spectra with those of several hydrodynamical simulations: the Horizon suite and the Magneticum simulation. We examine various contributions to the tSZ power spectrum across different redshifts, including the one- and two-halo term decomposition, the amount of bound gas, the importance of different masses, and the electron pressure profiles. Our comparison of the tSZ power spectrum reveals discrepancies between the halo model and cosmological simulations that increase with redshift. We find a 20% to 50% difference between the measured and predicted tSZ angular power spectrum over the multipole range ℓ = 10 3 − 10 4 . Our analysis reveals that these differences are driven by the excess of power in the predicted two-halo term at low k and in the one-halo term at high k . At higher redshifts ( z ∼ 3), simulations indicate that more power comes from outside the virial radius than from inside, suggesting a limitation in the applicability of the halo model. We also observe differences in the pressure profiles, despite the fair level of agreement on the tSZ power spectrum at low redshift with the default calibration of the halo model. In conclusion, our study suggests that the properties of the halo model need to be carefully controlled against real or mock data to be proven useful for cosmological purposes.

Ayçoberry, Emma (ORCID:0000000292351195)↗
DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
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 OSTI2022

MeV–GeV Gamma-Ray Emission from SNR G327.1–1.1 Discovered by the Fermi-LAT

Abstract We report the discovery of MeV–GeV γ -ray emission by the Fermi-LAT positionally coincident with the TeV pulsar wind nebula (PWN) HESS J1554–550 within the host supernova remnant (SNR) G327.1–1.1. The γ -ray emission is point-like and faint but significant (>4 σ ) in the 300 MeV–2 TeV energy range. We report here the Fermi-LAT analysis of the observed γ -ray emission followed by a detailed multiwavelength investigation to understand the nature of the emission. The central pulsar powering the PWN within G327.1–1.1 has not been detected in any wave band; however, it is likely embedded within the X-ray nebula, which is displaced from the center of the radio nebula. The γ -ray emission is faint and therefore a pulsation search to determine if the pulsar may be contributing is not feasible. Prior detailed multiwavelength reports revealed an SNR system that is old, τ ∼ 18,000 yr, where the interaction of the reverse shock with the PWN is underway or has recently occurred. We find that the γ -ray emission agrees remarkably well with a detailed broadband model constructed in a prior report based on independent hydrodynamical and semianalytic simulations of an evolved PWN. We further investigate the physical implications of the model for the PWN evolutionary stage incorporating the new Fermi-LAT data and attempt to model the distinct particle components based on a spatial separation analysis of the displaced PWN counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Proton acceleration in an overdense hydrogen plasma by intense CO 2 laser pulses with nonlinear propagation effects in the underdense pre-plasma

We report on proton acceleration from intense CO 2 laser-irradiated hydrogen plasmas at near-critical densities, with the density gradient steepened by Nd:YAG laser ablation-driven hydrodynamic shocks. While the experimental results, such as the quasi-monoenergetic proton spectra and their scaling with respect to the laser energy, are generally in agreement with the simulations, certain laser shots produced significantly higher proton energies than anticipated during the experiment. The increased proton energy may be linked to nonlinear propagation effects in the steepened plasma density ramp before the critical surface, including relativistic self-focusing and, for the case of temporally-structured laser pulses observed in the experiment, focusing of the trailing pulse through the plasma channel formed by the leading pulse 25 ps ahead. The occurrence of channel focusing in the underdense hydrogen plasma is supported by a subsequent pump-probe experiment with a dark-field imaging technique, where the formation of ion channels was observed after the passage of an intense CO 2 laser pulse

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗