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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗
DOE OSTI2022

Meso-Scale Electrodynamic Coupling of the Earth Magnetosphere-Ionosphere System

Within the fully integrated magnetosphere-ionosphere system, many electrodynamic processes interact with each other. We review recent advances in understanding three major meso-scale coupling processes within the system: the transient field-aligned currents (FACs), mid-latitude plasma convection, and auroral particle precipitation. (1) Transient FACs arise due to disturbances from either dayside or nightside magnetosphere. As the interplanetary shocks suddenly compress the dayside magnetosphere, short-lived FACs are induced at high latitudes with their polarity successively changing. Magnetotail dynamics, such as substorm injections, can also disturb the current structures, leading to the formation of substorm current wedges and ring current disruption. (2) The mid-latitude plasma convection is closely associated with electric fields in the system. Recent studies have unraveled some important features and mechanisms of subauroral fast flows. (3) Charged particles, while drifting around the Earth, often experience precipitating loss down to the upper atmosphere, enhancing the auroral conductivity. Recent studies have been devoted to developing more self-consistent geospace circulation models by including a better representation of the auroral conductance. It is expected that including these new advances in geospace circulation models could promisingly strengthen their forecasting capability in space weather applications. The remaining challenges especially in the global modeling of the circulation system are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Detectability of standing accretion shock instabilities activity in supernova neutrino signals

We introduce a novel methodology for establishing the presence of standing accretion shock instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice- based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test- statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda et al. [Astrophys. J. 851, 62 (2017).] as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to ~6 kpc for IceCube and up to ~3 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Signatures of an energetic charged body streaming in a plasma

A charged body moving in a plasma can excite a variety of linear and nonlinear waves in the form of trailing wakes, fore-wake shocks, and precursor solitons. These structures can further interact with the background plasma to create secondary effects that can serve as signatures of the passage of the charged body. Using particle-in-cell simulations, we carry out a basic investigation of the dynamics of a plasma system that is being traversed by an energetic charged body. Using two different shapes of this charged source, namely, an idealized infinite length planar source and a two-dimensional thin rectangular source, we examine the differences in the nature of the excited wave structures and their consequent impact on the background plasma. Our simulations reveal interesting features such as the dependence of the precursor speeds on the total charge on the driving source, local particle trapping, and energization of the trapped particles in various regions along the traversal path leading to the formation of energetic charged beamlets. So our basic findings could find practical applications such as in analyzing the trajectories of charged objects like space debris orbiting in the ionosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An experimental and detailed kinetic modeling study of the pyrolysis and oxidation of allene and propyne over a wide range of conditions

Allene and propyne are important intermediates in the pyrolysis and oxidation of higher hydrocarbon fuels, and they are also a major source of propargyl radical formation, which can recombine into different C 6 H 6 isomers and finally produce soot. In a prior work (Panigrahy et al., “A comprehensive experimental and improved kinetic modeling study on the pyrolysis and oxidation of propyne”, Proc. Combust. Inst 38 (2021)), the pyrolysis, ignition , and laminar flame speed of propyne were investigated. To understand the kinetic features of initial fuel breakdown and oxidation of the two C 3 H 4 isomers, new measurements for allene pyrolysis and oxidation are conducted in the present paper at the same operating conditions as those studied previously for propyne. Ignition delay times of allene are measured using a high-pressure shock tube and a heated twin-opposed piston rapid compression machine in the temperature range 690– 1450 K at equivalence ratios of 0.5, 1.0 and 2.0 in ‘air’, and at pressures of 10 and 30 bar. Pyrolysis species measurements of allene and propyne are also performed using a gas chromatography integrated singlepulse shock tube in the temperature range 1000–1700 K at pressure of 2 and 5 bar. Furthermore, laminar flame speeds of allene are measured at elevated gas temperatures of 373 K at pressures of 1 and 2 bar for a wide range of equivalence ratios from 0.6 to 1.5. A newly updated kinetic mechanism developed for this study is the first model that can well reproduce all of the experimental results for both allene and propyne. It is observed that in the pyrolysis process, allene dissociates faster than propyne. Both isomers exhibit similar ignition delay times at high temperatures (>1000 K), while, at intermediate temperatures (770–1000 K) propyne is the faster to ignite, and at lower temperatures (< 770 K) allene becomes more reactive. Furthermore, laminar flame speeds for propyne are found to be slightly faster than those for allene under the conditions studied in this work

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

Shocks propagate in a 2D dusty plasma with less attenuation than due to gas friction alone

In a dusty plasma, an impulsively generated shock, i.e., blast wave, was observed to decay less than would be expected due to gas friction alone. In the experiment, a single layer of microparticles was levitated in a radio frequency glow-discharge plasma. In this layer, the microparticles were self-organized as a 2D solid-like strongly coupled plasma, which was perturbed by the piston-like mechanical movement of a wire. To excite a blast wave, the wire's motion was abruptly stopped, so that the input of mechanical energy ceased at a known time. It was seen that, as it propagated across the layer, the blast wave's amplitude persisted with little decay. This result extends similar findings, in previous experiments with 3D microparticle clouds, to the case of 2D clouds. In our cloud, out-of-plane displacements were observed, lending support to the possibility that an instability, driven by wakes in the ion flow, provides energy that sustains the blast wave's amplitude despite the presence of gas damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Synthesis of Si-Fe Chondrule-like Dust Analogues in RF Discharge Plasmas

Chondrules are tiny particles that occur in stony meteorites and are considered as the building blocks of early asteroids and planets. It is believed that they were formed by the fast heating of the dust in the solar nebula. To date, there is no lab-scale experimental study of the formation of chondrules from the initial gas phase precursors following fast heating and crystallisation. The motivation of this work is a pre-trial study of the formation of chnodrule-like particles. The formation of meteorites in the space environment is associated with the aggregation of small particles or molecular clouds under the influence of shock waves or high-energy gas discharges in the solar nebula. In this work, the properties of product formation at the nanoscale-level were investigated using different feedstock materials which are the dominant elements in the meteorite. The structural and morphological properties of the synthesised Si-Fe nanomaterials were analysed by scanning/transmission electron microscopy (SEM/TEM), and chemical composition was analysed by X-ray energy-dispersive spectroscopy (EDS). The identification of crystalline phases was carried out by X-ray diffraction (XRD), whereas the presence of an Fe-Si system in the synthesised particles was demonstrated by Mössbauer spectroscopy. The obtained materials were exposed to the relatively high-energy pulsed plasma beam on the substrate with the aim to emulate the possible fast heating and melting of the formed nanoparticles. The formation steps of growing synthetic (engineered) chondro-like particles and nanostructures in laboratory conditions is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Thermomechanical Degradation of Sintered Copper Under High-Temperature Thermal Shock

The need for reliable bonded interface materials is critical to realize the performance benefits of wide-bandgap devices in power electronics modules, especially in operating temperatures greater than 150 degrees C. In this paper, we investigate the thermomechanical performance of sintered copper (Cu) as a large-area attachment, bonded between Cu baseplates and active-metal-bonded substrates, under accelerated thermal shock (-40 degrees C to 200 degrees C) conditions. In the fabrication phase of the samples, we used different stencil patterns and found that the grid and stripe patterns resulted in a better outgassing of the residual organics during the sintering process, thereby ensuring a substantially improved bond quality compared to a full-area print. The paste consisted of Cu microflakes, and we performed sintering using a Budatec SP300 sintering press at 275 degrees C with 15 MPa of bonding pressure for 5 minutes in a nitrogen atmosphere. Under accelerated experiments, we monitored the degradation of the sintered Cu bond in the samples through C-mode scanning acoustic microscope (C-SAM) images. To quantify the defect percentage in C-SAM images, we investigated image denoising techniques to exclude the pattern prints. Finally, we cross-sectioned a sample and obtained scanning electron microscope images, which revealed adhesive fracture as the dominant failure mechanism.

ADVANCED PROPULSION SYSTEMS,INORGANIC, ORGANIC, PH↗
DOE OSTI2020

Preshot Report for the NIF Discovery Science DDCyl Campaign: D_Hyd_DDCyl_DRT_DDD/EEE/FFF

The first Discovery Science shot day with the NIF DDCyl campaign demonstrated that Rayleigh-Taylor instability growth during the deceleration phase is scale-invariant for targets that varied in radial dimension and implosion time-scale by a factor of three. Smaller cylindrical implosions at the OMEGA laser facility were compared to the larger cylinders fielded in our first NIF outing. Both sets of targets had an initial m = 20 sinusoidal perturbation machined on the inner surface of an embedded aluminum marker layer, with the amplitude larger by the same factor of three for the NIF-scale targets. We observed growth factors of 14 at late time (t=λ = 7:7 ns, where λ = 1 for OMEGA and λ = 3 for NIF). These used a high density (300 mg/cm³) CH foam core, providing a very long period of deceleration before the rebounding shock hit the incoming spike tips, but the final convergence was only CR≈ 2:25. For our second shot day with this campaign, we are attempting to increase the convergence ratio to CR≥ 5 before reshock, while still measuring instability growth during the convergence/deceleration phase. As discussed in previous reports, the best way to increase convergence before reshock is to lower the density of the central foam. The preshot designs called for the use of a 30 mg/cm³ CH foam, though the as-built foams are closer to 40-42 mg/cm³, which is predicted to result in a final CR≈ 5:1 before reshock. Other than the lower foam density, the target design is identical to those used previously.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

The Impact of the 8–10 March 2012 Geomagnetic Storm on Inner Zone Protons as Measured by Van Allen Probes

The Relativistic Electron Proton Telescope (REPT) instrument on the Van Allen Probes observed a double-peaked inner zone proton population throughout the 7 year lifetime of the mission. M. Hudson et al. (2023) showed that a strong SEP event accompanied by a CME-shock in early March 2012 provided the Solar Energetic Proton (SEP) source for the higher L trapped proton population, which then diffused radially inward to be observed by REPT at L = 2. The study followed trajectories of SEP protons launched isotropically from a sphere at 7 Re for 2.5 hr in fields calculated by the LFM-RCM global MHD model, which includes electric fields needed to model the transport and trapping of the protons by the shock, and then a radial diffusion simulation was run for 2 years using the result from the test-particle simulation as the initial condition. The simulation result was compared with REPT measurement in November 2013 and showed reasonable agreement. However, the simulation overestimated the Phase Space Density by a factor of four due to lack of field line curvature scattering during the storm in the model. In this study, a test-particle simulation is performed for 2 days following the injection and trapping of protons in March 2012 using TS05 fields to simulate the field line curvature scattering of the trapped SEP due to the buildup of the ring current during the geomagnetic storm. The resulting sample distribution was then weighted using the flux at the end of the two-hour MHD-test particle simulation. A radial diffusion simulation is then run using the initial profile that included the loss effect, with improved comparison with REPT measurements after 2 years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Growth of concomitant laser-driven collisionless and resistive electron filamentation instabilities over large spatiotemporal scales

Collective processes in plasmas often induce microinstabilities that play an important role in many space or laboratory plasma environments. Particularly of note is the Weibel-type current filamentation instability, which is believed to drive the creation of collisionless shocks in weakly magnetized astrophysical plasmas. In this paper, this instability class is studied through interactions of ultraintense and short laser pulses with solid foils, leading to localized generation of megaelectronvolt electrons. Proton radiographic measurements of both low- and high-resistivity targets show two distinct, superimposed electromagnetic field patterns arising from the interpenetration of the megaelectronvolt electrons and the background plasma. Particle-in-cell simulations and theoretical estimates suggest that the collisionless Weibel instability building up in the dilute expanding plasmas formed at the target surfaces causes the observed azimuthally symmetric electromagnetic filaments. For a sufficiently high resistivity of the target foil, an additional resistive instability is triggered in the bulk target, giving rise to radially elongated filaments. The data reveal the growth of both filamentation instabilities over large temporal (tens of picoseconds) and spatial (hundreds of micrometres) scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Flat Spectra of Energetic Particles in Interplanetary Shock Precursors

The observed energy spectra of accelerated particles at interplanetary shocks often do not match the diffusive shock acceleration (DSA) theory predictions. In some cases, the particle flux forms a plateau over a wide range of energies, extending upstream of the shock for up to seven flux e-folds before submerging into the background spectrum. Remarkably, at and downstream of the shock we have studied in detail, the flux falls off in energy as ϵ -1 , consistent with the DSA prediction for a strong shock. The upstream plateau suggests a particle transport mechanism different from those traditionally employed in DSA models. We show that a standard (linear) DSA solution based on a widely accepted diffusive particle transport with an underlying resonant wave–particle interaction is inconsistent with the plateau in the particle flux. To resolve this contradiction, we modify the DSA theory in two ways. First, we include a dependence of the particle diffusivity κ on the particle flux F (nonlinear particle transport). Second, we invoke short-scale magnetic perturbations that are self-consistently generated by, but not resonant with, accelerated particles. They lead to the particle diffusivity increasing with the particle energy as ∝ϵ 3/2 that simultaneously decreases with the particle flux as 1/F. The combination of these two trends results in the flat spectrum upstream. We speculate that nonmonotonic spatial variations of the upstream spectrum, apart from being time-dependent, may also result from non-DSA acceleration mechanisms at work upstream, such as stochastic Fermi or magnetic pumping acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Shock compression of diamond single crystals to 120 GPa: Refractive index and nonlinear photoelasticity

The optical response of transparent solids at extreme conditions is important for both fundamental science and many applications. Strong transparent solids are of particular interest for use as optical windows in dynamic compression experiments. Due to diamond’s exceptional strength and optical properties, laser-driven shock experiments and plate impact experiments were carried out to examine the diamond optical response for shock wave compression along two different crystal orientations. Using laser interferometry at 532 nm and 1550 nm wavelengths, optical transparency was observed and refractive indices were determined for [100] diamond at stresses up to 119 GPa and for [111] diamond at stresses up to 87 GPa. From these results, the nonlinear photoelastic response for [100] and [111] diamond was determined, revealing significant dependence on both crystal orientation and laser wavelength. To enable [100] diamond as an interferometry window in dynamic compression experiments, the requisite window corrections were determined for 532 nm and 1550 nm wavelengths. Because of diamond’s excellent x-ray transparency, the present findings will be particularly useful for incorporating [100] diamonds as windows in dynamic compression experiments involving x-ray diffraction or other x-ray diagnostics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Understanding the transport and break up of reactive ejecta

This manuscript investigates reactive- versus hydrodynamic-breakup processes of ejecta. For this paper, the reactive metal is cerium (Ce) and the nonreactive metal is tin (Sn), the nonreactive gas is helium (He) and the reactive gas is deuterium (D 2 ) or hydrogen (H 2 ). Experiments were performed in vacuum and the reactive- and nonreactive-gases at various pressures, where we endeavored to match the post-shock gas densities to differentiate between reactive- versus hydrodynamic-breakup processes. Hydrodynamic breakup sensitively links to the Weber number (gas density, liquid fragment diameter, surface tension, and the square of the relative velocity between the fragment and the gas), whereas reactive breakup links to the reactive dynamics which includes two processes. In one case the reactive metal breaks up into smaller fragments as rapidly as the reaction rate, and in the other a crust grows on the liquid fragments as the reactions occur, a diffusion limited process. In the latter case, the particle diameters increase with time as the crust grows. In this process, which is indicated by the data, particles breakup as the CeD 2 loses strength with increasing temperature, leaving an exponentially increasing diameter.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mode-Selective Vibrational Energy Transfer Dynamics in 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX) Thin Films

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. In this work, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N–N and axial NO2 stretching modes with the long-lived, excited phonon bath.

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

Global heliospheric termination shock strength in the solar–interstellar interaction

A heliospheric termination shock (HTS) surrounds our Solar System at approximately 100 astronomical units from the Sun, where the expanding solar wind (SW) is compressed and heated before encountering the interstellar medium. HTS-accelerated particles govern the pressure balance with the interstellar medium, but little is known about the global properties of the HTS beyond in situ measurements from Voyager in only two directions of the sky. Here we fill this gap by extracting the HTS strength using particle-in-cell, test particle and magnetohydrodynamic simulations, constrained by Interstellar Boundary Explorer observations of energetic neutral atoms produced from HTS-accelerated particles. Our results reveal there is a higher compression near the poles during solar minimum compared with solar maximum due to the higher Mach number flow. North–south asymmetries arise from the disparate evolution of the polar coronal holes, while minimum compression near the flanks is probably due to SW slowing from mass loading over a greater distance to the HTS. The results imply a strong connection between the HTS strength and the SW and interstellar medium dynamics.

Astrophysical plasmas↗
DOE OSTI2025

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Measurements of enhanced performance in an indirect drive inertial confinement fusion experiment when reducing the contact area of the capsule support

Experimental results from indirectly driven inertial confinement fusion experiments testing the performance gained from using an alternate capsule tent support are reported. The polar tent describes an alternate geometry for the thin membrane used to support the Deuterium–Tritium (DT) filled capsule. In this work, the contact area is reduced by 23 times by locating the tent support close to the poles of the capsule. The polar tent experiments are repeats of previous 3 shock 1.63 MJ, 400 TW high foot experiments and use a 165 μm thick silicon doped carbon hydrogen plastic (CH) shell. Using the polar tent support, we report a DT neutron yield of 1.07 ×1016, 76% higher than the expected Y DT ∝V 7.7 scaling. This is, at the time of writing, the highest neutron yield to date from a CH shell implosion. Furthermore, we find that the inferred pressure when using the polar tent is significantly above the model based on analytic scaling even when accounting for tent effects. Analysis of x-ray and neutron images shows the reduction of lobes produced by nominal tent features. The reduction of these features in the polar tent experiments leads to decreased low mode (P2 and P4) asymmetry compared to the nominal tent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗