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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Vibrational and rotational analysis of the emission spectra of the arc jet flow

Applying atomic and molecular physics to the analysis of the radiation emitted from the arc jet flow provides a means for determining the species and excitation temperature of the constituents of the flow. A rotational and vibrational analysis of the spectra obtained from the radiation emitted in the shock layer and in the free stream of the jet flow was performed, specifically, in the shock layer bands of the First Negative Group of ionized molecular nitrogen and in the free stream of the gamma-system of nitric oxide.

Blount, Charles E.↗
NASA NTRS1990

A near-infrared spectroscopic study of the starburst core of M82

Near-IR spectroscopy of the M82 starburst core is presented, including complete J, H, and K band spectra with a resolution of 0.0035-micron for the inner 60 pc of the galaxy. Also, spatial profiles along the starburst ridge are presented for Br-gamma, molecular hydrogen, and forbidden Fe II line fluxes. Emission from shocked molecular hydrogen is detected from the core of M82. The distribution of features across the starburst disk are mapped to study the relationships between spectral diagnostics. The observations are used to test the appropriateness of single-beam, aggregate models for studying the physical conditions in starbursts.

Lester, D. F.↗
NASA NTRS1995

Organic molecules in the gas phase of dense interstellar clouds

Since a previous Committee on Space Research (COSPAR) review on this subject, the number of molecular species identified by astronomers in dense interstellar clouds or in the envelopes expelled by evolved stars has grown from about eighty to approximately one hundred. Recent detections in stellar envelopes include the radical CP, the second phosphorus-containing astronomical molecule; SiN, the first astronomical molecule with a Si-N bond; and the HCCN radical. In the dense interstellar clouds recent detections or verifications of previous possible identifications include the H3O(+) ion, which is a critical intermediary in the production of H2O and O2; the CCO radical, which is isoelectronic with HCCN; the SO(+) ion, which appears to be diagnostic of shock chemistry; two new isomers of cyanoacetylene, HCCNC and CCCNH; and the two cumulenes H2C3 and H2C4. Some recent work is also described on the mapping of interstellar clouds in multiple molecular transitions in order to separate variations in chemical abundance from gradients in physical parameters.

Irvine, W. M.↗
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↗
NASA NTRS1983

Nearfield observations of tones generated from supersonic jet flows

The physics of discrete tone generation (screech) from improperly expanded supersonic jets have been investigated experimentally. Both nearfield acoustic survey and phase-averaged schlieren flow visualization were conducted. The dominant mode of flow oscillation that occurred during intense screech and its coupling with the internal shock structure and external sound field were established. Deductions made from the observed coupling mechanism revealed the importance of Helmholtz number and shock cell spacing in the maintenance of screech. Measured amplitude variation of screech are compared with computed amplitude variation of the dominant instability. The results suggest that jet instability plays a direct role in the generation of screech.

Yu, J. C.↗
NASA NTRS2003

Effects of Pre-Existing Target Structure on the Formation of Large Craters

The shapes of large-scale craters and the mechanics responsible for melt generation are influenced by broad and small-scale structures present in a target prior to impact. For example, well-developed systems of fractures often create craters that appear square in outline, good examples being Meteor Crater, AZ and the square craters of 433 Eros. Pre-broken target material also affects melt generation. Kieffer has shown how the shock wave generated in Coconino sandstone at Meteor crater created reverberations which, in combination with the natural target heterogeneity present, created peaks and troughs in pressure and compressed density as individual grains collided to produce a range of shock mineralogies and melts within neighboring samples. In this study, we further explore how pre-existing target structure influences various aspects of the cratering process. We combine experimental and numerical techniques to explore the connection between the scales of the impact generated shock wave and the pre-existing target structure. We focus on the propagation of shock waves in coarse, granular media, emphasizing its consequences on excavation, crater growth, ejecta production, cratering efficiency, melt generation, and crater shape. As a baseline, we present a first series of results for idealized targets where the particles are all identical in size and possess the same shock impedance. We will also present a few results, whereby we increase the complexities of the target properties by varying the grain size, strength, impedance and frictional properties. In addition, we investigate the origin and implications of reverberations that are created by the presence of physical and chemical heterogeneity in a target.

Barnouin-Jha, O. S.↗
NASA NTRS2010

Instability of a Supersonic Boundary-Layer with Localized Roughness

A localized 3-D roughness causes boundary-layer separation and (weak) shocks. Most importantly, streamwise vortices occur which induce streamwise (low U, high T) streaks. Immersed boundary method (volume force) suitable to represent roughness element in DNS. Favorable comparison between bi-global stability theory and DNS for a "y-mode" Outlook: Understand the flow physics (investigate "z-modes" in DNS through sinuous spanwise forcing, study origin of the beat in DNS).

Marxen, Olaf↗
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↗