Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “moving shock”

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

The effects of mixture preburning on detonation wave propagation

Pressure gain combustion in the form of continuous detonations can provide a significant increase in the efficiency of a variety of propulsion and energy conversion devices. In this regard, rotating detonation engines (RDEs) that utilize an azimuthally-moving detonation wave in annular systems are increasingly seen as a viable approach to realizing pressure gain combustion. However, practical RDEs that employ non- premixed fuel and oxidizer injection need to minimize losses through a number of mechanisms, including turbulence-induced shock-front variations, incomplete fuel-air mixing, and premature deflagration. In this study, a canonical stratified detonation configuration is used to understand the impact of preburning on detonation efficiency. It was found that heat release ahead of the detonation wave leads to weaker shock fronts, delayed combustion of partially-oxidized fuel-air mixture, and non-compact heat release. Furthermore, large variations in wave speeds were observed, which is consistent with wave behavior in full-scale RDEs. Peak pressures in the compression region or near triple points were considerably lower than the theoretically-predicted values for ideal detonations. Analysis of the detonation structure indicates that this deflagration process is parasitic in nature, reducing the detonation efficiency but also leading to heat release far behind the wave that cannot directly strengthen the shock wave. As a result, this parasitic combustion leads to commensal combustion (heat release far downstream of the wave), indicating that it is the root cause of combustion efficiency losses.

42 ENGINEERING↗
DOE OSTI2021

Femtosecond Reflectance Spectroscopy for Energetic Material Diagnostics

Understanding the fundamental mechanisms underpinning shock initiation is critical to predicting energetic material (EM) safety and performance. Currently, the timescales and pathways by which shock-excited lattice modes transfer energy into specific chemical bonds remains an open question. Towards understanding these mechanisms, our group has previously measured the vibrational energy transfer (VET) pathways in several energetic thin films using broadband, femtosecond transient absorption spectroscopy. However, new technologies are needed to move beyond these thin film surrogates and measure broadband VET pathways in realistic EM morphologies. Herein, we describe a new broadband, femtosecond, attenuated total reflectance spectroscopy apparatus. Performance of the system is benchmarked against published data and the first VET results from a pressed EM pellet are presented. This technology enables fundamental studies of VET dynamics across sample configurations and environments (pressure, temperature, etc .) and supports the potential use of VET studies in the non-destructive surveillance of EM components.

42 ENGINEERING↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Explosive dispersal of particles in high speed environments

In this paper, we present the results of the explosive dispersal of particles in high-speed environments. We carry out Euler–Lagrange numerical simulations of a source at quiescent ambient conditions as well as moving at Mach numbers of 3 and 6. Particle volume fractions of 0%, 1%, and 4.5% are presented. The detonation profile is computed with the Jones–Wilkins–Lee equation of state using a reactive burn model. Non-static cases provide a framework to consider the effect of a bow shock and pre-existing high-speed flow conditions on the dispersal process. We also compute averages of both static and dynamic pressures, as well as impulse density histories on virtual probe planes to characterize the momentum of the flow and particles that would deposit on a target. Results suggest that the presence of the particles can have a substantial effect on the pressure average of the virtual target planes.

Durant, B. A.↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Interaction of Inner Heliosheath Ions with Electromagnetic Ion Cyclotron Waves

Abstract Inner heliosheath (IHS) ions are expected to be exposed to various waves, shocks, and turbulence, which can affect ion distributions and thus their charge-exchange rates with interstellar neutral atoms. This work addresses the potential significance of electromagnetic ion cyclotron (EMIC) waves under expected IHS conditions. From a kinetic dispersion relation, we find the possibility of frequent triggering of EMIC instability in the IHS. The threshold anisotropy of proton temperatures required for the instability is small, ( T ⊥ – T ∣∣ )/ T ∣∣ ≈ 0.1 or less, mainly due to high plasma β (ratio of the plasma to magnetic pressures). Numerical calculations on the scattering of ions (protons, He + , and He 2+ with energy of 0.01–50 keV) based on two models for EMIC waves with a moderate intensity indicate significant scattering in the pitch angle (mostly a few tens of degrees) and energy (mostly a few tens of percent) although details depend on the energy and pitch angle of each species and adopted EMIC wave models. This occurs on a short timescale (<100 times the gyro-period of each ion species). Resonant scattering in a few to a few tens of keV (corresponding to the expected pickup ion energy in the IHS) is easily expected unless the wavenumber is too large. The scattering effect is distinguished among different species such that for lower gyrofrequency ions (He + versus He 2+ and He 2+ versus protons), the main scattering effect moves toward a lower energy domain. All these results imply continuous disturbance of ion distributions by possibly prevailing EMIC waves in the IHS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Reflected shock-initiated ignition probed via simultaneous lateral and endwall high-speed imaging with a transparent, cylindrical test-section

A TRAnsparent cylindrical (TRAC) shock tube test section has been developed to explore the ignition structure behind reflected shockwaves via simultaneous lateral and endwall high-speed imaging using three fuel/oxidizer mixtures: 0.1% n-C 7 H 16 /1.10% O 2 /4.14% N 2 /94.66% Ar, 0.1% iC 8 H 18 /1.25% O 2 /4.69% N 2 /93.96% Ar, and 5% CH 4 /10% O 2 /85% CO 2 . During mild ignition, the first exothermic centers were located outside of the observable test section. It is shown that with increasing temperature the far-wall ignition events move closer to the endwall, until a strong ignition is achieved. The temperature gradients that promote mild ignition were quantified and were found to have a slight dependence on fuel with the RMS temperature gradients being 0.80% for n-heptane and 1.07% for isooctane. The impact of temperature gradients on metrics used to improve kinetic models, such as species time histories, is noteworthy. From this temperature gradient, an experimental correlation is provided to estimate the transition temperature from mild to strong ignition for a mixture. The effect of diagnostic location in non-homogeneous reactions is quantified in reference to ignition delay. In conclusion, the unique and insightful perspective on the ignition process under heavy shock bifurcation is also discussed.

02 PETROLEUM↗
DOE OSTI2024

Metals strengthen with increasing temperature at extreme strain rates

The strength of materials depends on the rate at which they are tested, as defects, for example dislocations, that move in response to applied strains have intrinsic kinetic limitations. As the deformation strain rate increases, more strengthening mechanisms become active and increase the strength. However, the regime in which this transition happens has been difficult to access with traditional micromechanical strength measurements. Here, with microballistic impact testing at strain rates greater than 10 6 s −1 , and without shock conflation, we show that the strength of copper increases by about 30% for a 157 °C increase in temperature, an effect also observed in pure titanium and gold. This effect is counterintuitive, as almost all materials soften when heated under normal conditions. This anomalous thermal strengthening across several pure metals is the result of a change in the controlling deformation mechanism from thermally activated strengthening to ballistic transport of dislocations, which experience drag through phonon interactions. These results point to a pathway to better model and predict materials properties under various extreme strain rate conditions, from high-speed manufacturing operations to hypersonic transport.

36 MATERIALS SCIENCE↗
DOE OSTI2023

NucD-24 Nuclear Diagnostics Development: Time resolved 12 C-rhoR on GCD and NIS technology development for polar LOS and data quality [Slides]

The Neutron Imaging team is readying for the NIS2 polar LOS development and future improvements to data quality on all LOS given consistently higher yields at NIF. Next-gen nano guide scintillators for NIS2 and CMOS cameras to replace outdated CCDs will be studied. Digital image plates will be tested for better data quality on passive imaging systems. The transition into the ignition regime is expected to show a significant signature in the ablator areal density. In the marginal ignition regime, a decrease of 12C-rhoR (ablator areal density) is expected but once the burn propagation produces a shock and compresses the ablator the trend reverses. GCD with the PD-PMT can resolve this effect in the ablator areal density by a temporal shift in the carbon signal. To give quantitative results a calibration of the detector for the Hohlraum/TMP signal is needed. Shots using Al, Si, and Au pucks will yield the needed information to move forward with this project.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Understanding the evolution of liquid and solid microjets from grooved Sn and Cu samples using radiography

Experiments were performed on grooved Sn and Cu samples to study the temporal evolution of microjets. Jets were generated by the impact of gun-launched flyer plates against the back of grooved targets made from either Cu or Sn (groove depth of ~250 μm). The Hugoniot states in the various Sn targets encompassed conditions where solid phases are maintained throughout (7 and 16 GPa) and also conditions where melting occurs upon the release of compression (25 and 34 GPa); the transition occurs near a Hugoniot pressure of 23 GPa. Cu targets at 27 and 56 GPa provide comparisons in which the jets move at similar speeds but remain solid. In all cases, the spatial distribution of mass within the microjets was measured using high-speed synchrotron radiography. The result is a time history of the jet thickness profile from which quantities like total jet mass and jet velocity can be derived. In both the solid and liquid states, we generally observe that an increase in the shock strength leads to an increase in jet mass. However, this trend breaks down for Hugoniot states near the transition from continuously solid to melted-on-release. This is evidenced by the observation that there was no difference in the rate of mass flow in Sn jets at 16 and 25 GPa, while similar pressure jumps on either side of this range caused substantial changes in the jet mass. This contrasts with the behavior of smaller polishing defects that were present on the same samples (~1 μm deep). From these, no ejecta mass was detected below the melt boundary, but obvious microjets were generated once melting occurred. This indicates that crossing the bulk melt-on-release threshold can alternately promote or inhibit the flow of mass into microjets based on the amplitude of the initial perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Regional seismic velocity changes following the 2019 M w 7.1 Ridgecrest, California earthquake from autocorrelations and P / S converted waves

SUMMARY We examine regional transient changes of seismic velocities generated by the 2019 Mw 7.1 Ridgecrest earthquake in California, using autocorrelations of moving time windows in continuous waveforms recorded at regional stations. We focus on traveltime differences in a prominent phase generated by an interface around 2 km depth, associated with transmitted Pp waves and converted Ps waves from the ongoing microseismicity. Synthetic tests demonstrate the feasibility of the method for monitoring seismic velocity changes. Taking advantage of the numerous aftershocks in the early period following the main shock, we obtain a temporal resolution of velocity changes up to 20 min in the early post-main-shock period. The results reveal regional coseismic velocity drops in the top 1–3 km with an average value of ∼2 per cent over distances up to 100 km from the Ridgecrest event. These average velocity drops are likely dominated by larger changes in the shallow materials and are followed by rapid recoveries on timescales of days. Around the north end of the Ridgecrest rupture and the nearby Coso geothermal region, the observed coseismic velocity drops are up to ∼8 per cent. The method allows monitoring temporal changes of seismic velocities with high temporal resolution, fast computation and precise spatial mapping of changes. The results suggest that significant temporal changes of seismic velocities of shallow materials are commonly generated on a regional scale by large events.

58 GEOSCIENCES↗
DOE OSTI2022

An Investigation into the Surface Skidding Response of PBX 9501 and PBX 9502

While polymer-bonded explosives, PBX, have reduced sensitivity to ignition from mechanical shock or heating compared to conventional explosives, the investigation of the mechanical ignition mechanisms for PBX remains vital to assessing the safety during machining and general handling. Under frictional heating, hot spot generation due to high melting point contaminants in the interface is a suspected source for increased probability of detonation. Here, to investigate such frictional heating, samples of PBX 9501 and PBX 9502 were impacted and skidded against a tangentially moving target and thermally imaged. Temporally resolved temperature and forces were simultaneously measured with and without 300 µm silica grit at the frictional interface. A trend of increasing temperature was observed as the speed and tangential force on the sample increased. Grit particles in the frictional interface were found to act as frictional heat concentration spots owing to the grit's protrusion from the surface and lesser susceptibility to melting compared to the PBX. The coefficient of friction between PBX 9501 and the skidding surface was observed to be largely constant for forces and speeds at the lower end of the test range. In contrast, at high tangential speeds, the coefficient was found to be significantly lower.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy

A significant number of hypervelocity stars with velocities between 1500 − 2500 km s −1 have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new AREPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M ⊙ remains bound, moving at a velocity of ∼2800 km s −1 . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.

Astronomy and AstroPhysics↗
DOE OSTI2021

Comments on the SURFplus reactive burn model

The SURFplus reactive burn model is intended for heterogeneous solid HE that produce an excess amount of carbon. It utilizes a fast hotspot reaction (SURF model) plus a second slow reaction for the energy release from carbon clustering. The fast reaction dominates shock initiation, though the fast rate parameters are slightly affected by the late energy release from the slow reaction. The slow reaction has a significant affect on propagating detonation waves. Propagating detonation waves are characterized by the curvature effect; i.e., detonation speed $D_n$ as function of the local front curvature κ. The $D_n(κ)$ curves have a qualitative different shape. This results from the sensitivity of the slope of the $D_n(κ)$ curve to the reaction-zone width up to the sonic point, and the sonic point moving as κ increases from the end of the slow reaction to the neighborhood of the end of the fast reaction. For PBX 9502, the steep slope for small curvature is determined by the slow reaction, while the lower slope for larger curvature is determined by the fast reaction rate in the high detonation pressure regime. The remainder of this paper is organized as follows. In section 2 we review the assumptions of the Shaw and Johnson carbon cluster model used by the SURFplus model for the slow reaction. Section 3 describes in detail the mixture model for a partly burned HE. The form of the slow rate used by the SURFplus model is described in section 4. We conclude with a brief discussion of applying a burn model with 2 rates to either a metalized or a composite PBX.

42 ENGINEERING↗
DOE OSTI2021

Shear Strength and Release of Large Cryogenic Pellets from the Barrel of a Shattered Pellet Injector for Disruption Mitigation

Shattered pellet injection (SPI) has been chosen as the baseline disruption mitigation system on ITER due to its ability to rapidly inject material deep into the plasma to greatly increase the plasma density and radiate the thermal energy. SPI utilizes a mechanical punch or high-pressure gas to release and accelerate a pellet that has been cryogenically desublimated in the barrel of a pipe gun. Various material injection combinations could possibly be implemented during different phases of a disruption event to radiate plasma energy, reduce electromagnetic loads on machine components, avoid the formation of runaway electrons, or to dissipate runaway electrons that form. Each injection phase could possibly utilize combinations of deuterium, neon, or argon. In this paper we outline experimental measurements of pellet material shear strength at SPI operating temperatures to understand the force needed to release SPI pellets. Deuterium, neon, argon, and deuterium-neon mixture pellets with diameters of 8.5, 12.5, and 15.7 mm are formed at a range of relevant gas pressures and temperatures and dislodged from the cold zone with a slow-moving piston driven by a motor. The slow-moving piston is kept above the triple point temperature of the material while the pellet is forming, then cooled to below the triple point temperature before contacting the pellet to minimize any thermal conduction to the pellet. The piston incorporates a load cell to measure the force applied when the pellet breaks away from the cold zone in the barrel.The ability of the gas and punch methods to exceed the shear strength of the studied pellet materials for release has been analyzed. High pressure gas delivered by fast opening valves produce pressure shock to the pellet due to supersonic expansion of the propellant gas. Pressure (and therefore, force) oscillations are present due to transverse density propagation throughout the breech volume. Mechanical punches deliver an impact force through a high-kinetic energy impact. The effect of the mechanical shock on the pellet has been explored and will be presented in this paper. Scaling to larger ITER size SPI pellets will be described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of the Mach Stem Triple Point

Among the more complex characteristics of an explosion’s shock front is the creation of the Mach stem and its triple point. The scientific footage of the US Atmospheric Nuclear Weapons Testing Program has a film that may show the locus of points formed by the triple point, allowing us to visually investigate the trajectory of these shockwaves. During a nuclear blast, once the X-ray diffusion process ends, an incident shockwave forms in all directions. When the incident shockwave contacts the ground, a reflected shockwave is created. The overall magnitude and impact of these two shockwaves are proportional to the increasing total pressure and temperature, allowing the disturbances to move faster than the local speed of sound.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Parallel Implicit Hydrodynamics for High Explosive Burn Calculations

High explosives in hostile environments will require calculational capabilities that model processes, which evolve on timescales from minutes to nanoseconds. Eventually the HE will begin to move metal. To handle this temporal evolution an implicit hydrodynamics coupled to the chemical release of the HE energy is required. In addition, the use of chemical kinetics to model the transition from, the initially, slow heating of a confined high explosive through to deflagration and on to detonation requires many computational zones to model high explosive engineering systems. This requirement means that a fully parallel implicit hydrodynamics is essential. In this paper we present the calculation of a nonlinear matrix equation for the advanced particle pressure that has been made parallel and implemented in our AMR code, BABBO. This new parallel implicit hydrodynamics has been applied to a cookoff problem, as well as, one and two dimensional shock problems. Results are presented and discussed.

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

CO Laser Absorption Measurements During Syngas Combustion at High Pressure

Abstract Syngas is a desirable fuel for combustion in the Allam-Fetvedt cycle, which involves combustion under supercritical-CO2 conditions. While some work has been conducted in collecting ignition delay times (IDT) at the extreme pressures required by these systems, significant model deficiencies remain. Additionally, considerable barriers in terms of nonideal gas dynamic effects have been shown for these experiments in shock tubes. Further investigation into the fundamental combustion kinetics of H2/CO/CO2 mixtures is required. Time-resolved speciation measurements for target species have been shown to better aid in improving the understanding of underlying chemical kinetics than global ignition delay time measurements. Therefore, laser absorption measurements of CO were measured behind reflected shock waves during combustion of syngas at 5 and 10 bar and temperatures between 1080 and 2100 K. The mixtures investigated utilized H2-to-CO ratios of 1:1 and 1:4, respectively, each at stoichiometric conditions, allowing for discussions of the effect of initial fuel composition. A ratio of fuel to CO2 of 1:2 was also utilized to represent commercially available syngas. The mixtures were diluted in helium and argon (20% He, 76.5% Ar) to minimize thermal effects and to expedite CO thermal relaxation during the experiment. The resulting CO time histories were then compared to modern chemical kinetics mechanisms, and disagreement is seen for this system, which is assumed to be fairly well known. This study elucidates particular chemistry that needs improvement in moving toward a better understanding of syngas combustion at elevated pressures.

Engineering↗