Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock chemistry”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Time-dependent simulation of reflected-shock/boundary layer interaction

An initial numerical/experimental investigation has been done to better understand multi-dimensional flow phenomena inside pulse facilities. Time-dependent quasi-one-dimensional and axisymmetric numerical simulations of complete shock tube flow are compared with experimental pressure traces recorded at the NASA Ames electric-arc driven shock tube facility (from cold driver shots). Of particular interest is the interaction between the reflected shock wave and the boundary layer. Evidence of the shock bifurcation caused by this interaction is clearly seen in the present experimental data. The axisymmetric simulations reproduce this phenomenon and demonstrate how this interaction can provide a mechanism for driver gas to contaminate the stagnation region. The simulations incorporate finite-rate chemistry, a moving mesh and laminar viscosity.

Wilson, Gregory J.↗
NASA NTRS1995

Computational Study of Flow Establishment in Hypersonic Pulse Facilities

This paper presents a study of the temporal evolution of the combustion flowfield established by the interaction of ram-accelerator-type projectiles with an explosive gas mixture accelerated to hypersonic speeds in an expansion tube. The Navier-Stokes equations for a chemically reacting gas are solved in a fully coupled manner using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order, total variation diminishing (TVD) scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of expansion tube flowfields are then presented for two different configurations. In particular, the development of the shock-induced combustion process is followed. In one case, designed to ensure ignition only in the boundary layer, the lateral extent of the combustion front during the initial transient phase was surprisingly large. The time histories of the calculated thrust and drag forces on the ram accelerator projectile are also presented.

Yungster, S.↗
NASA NTRS2011

Study of Unsteady, Sphere-Driven, Shock-Induced Combustion for Application to Hypervelocity Airbreathing Propulsion

A premixed, shock-induced combustion engine has been proposed in the past as a viable option for operating in the Mach 10 to 15 range in a single stage to orbit vehicle. In this approach, a shock is used to initiate combustion in a premixed fuel/air mixture. Apparent advantages over a conventional scramjet engine include a shorter combustor that, in turn, results in reduced weight and heating loads. There are a number of technical challenges that must be understood and resolved for a practical system: premixing of fuel and air upstream of the combustor without premature combustion, understanding and control of instabilities of the shock-induced combustion front, ability to produce sufficient thrust, and the ability to operate over a range of Mach numbers. This study evaluated the stability of the shock-induced combustion front in a model problem of a sphere traveling in a fuel/air mixture at high Mach numbers. A new, rapid analysis method was developed and applied to study such flows. In this method the axisymmetric, body-centric Navier-Stokes equations were expanded about the stagnation streamline of a sphere using the local similarity hypothesis in order to reduce the axisymmetric equations to a quasi-1D set of equations. These reduced sets of equations were solved in the stagnation region for a number of flow conditions in a premixed, hydrogen/air mixture. Predictions from the quasi-1D analysis showed very similar stable or unstable behavior of the shock-induced combustion front as compared to experimental studies and higher-fidelity computational results. This rapid analysis tool could be used in parametric studies to investigate effects of fuel rich/lean mixtures, non-uniformity in mixing, contaminants in the mixture, and different chemistry models.

Axdahl, Erik↗
DOE OSTI2020

Nucleation and growth of voids in shock loaded copper bicrystals

Understanding the evolution of damage and deformation due to spall at grain boundaries can provide a basis for connecting micro- to macroscale failure behavior in metals under extreme conditions. Copper bicrystal samples were shock loaded using flyer-plate impacts in a light gas gun with shock stresses ranging from 3 to 6 GPa. Pulse duration as well as crystal orientation along the shock direction were varied for a fixed boundary misorientation to determine their effects on void nucleation and coalescence. Samples were soft recovered and cross-sectioned to characterize damage using electron backscattering diffraction and scanning electron microscopy to gather information on damage characteristics at and around the GB, with emphasis on growth of boundary and bulk voids. Chemistry and composition analysis were also performed on samples to determine if trace elements present in a sample affected the threshold for void nucleation. Results show that the kinetics of damage growth at the boundary are strongly affected by stress level and impurities. It was found that the boundary selected had a similar or even lower tendency to show damage than the bulk at low pulse durations and amplitudes. As pulse duration and amplitude increased damage localized at the boundary, which was found to consist of many small voids, indicating that the boundary experienced rapid void nucleation and coalescence. Furthermore, the presence of impurities correlated strongly with scatter on damage evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental and fuel-surrogates modeling study of the high-pressure pyrolysis of specialty cetane number fuels: implications for fall-off in ethylene unimolecular dissociation

Single pulse shock tube experiments were conducted at 50 atm nominal pressure and 4 ms nominal reaction time over a temperature range of 900–1800 K, to study the pyrolysis speciation of a multi-component jet fuel, F-24, and six cetane number (CN) specialty fuels - CN30, CN35, CN40, CN45, CN50, and CN55. Gas chromatography (GC) was used to qualitatively and quantitatively analyze the post shock gases. The relationship between the formation of key pyrolysis species and the chemically controlled combustion propensity as reflected by the cetane number of each fuel was examined. A surrogate-based mechanism from the CRECK Modelling Group and chemical-functional group based optimized surrogates (CFGO) were used to simulate the pyrolysis speciation results. The model was able to capture the chemistry of most species except two important pyrolysis intermediates – ethylene and acetylene. Chemical kinetic analyses were performed to identify the important reactions which affect the chemistry of these species; however, the rate parameters of critical reactions were found to be unsuitable for simulating the present high-pressure studies. Here, to address this unsuitability, a theory-based fall-off analysis for three reactions representing the decomposition of ethylene and subsequent formation of acetylene was performed, and these are included in an updated version of the CRECK mechanism. This update resolves discrepancies between the experimental results and simulations for ethylene and acetylene. Reaction flux analyses using the updated surrogate model were also performed to identify the important reaction pathways responsible for the formation of crucial species and to provide an analysis of the chemistry of complex multi-component fuel systems. The fundamental reactions responsible for driving pyrolysis chemistry were greatly influenced by the chemical functional groups present in these fuels. In addition to updating the rate parameters of specific reactions to improve modeling, this study also emphasizes the effectiveness of the fuel-surrogate approach, where surrogates representing the chemical functional group composition of the parent fuel serve as a valuable tool for predicting the combustion chemistry of novel fuels.

Chemical Kinetics↗
NASA NTRS1978

Studies of Brazilian meteorites. XIII - Mineralogy, petrology, and chemistry of the Putinga, Rio Grande do Sul, chondrite

The Putinga, Rio Grande do Sul chondrite is described and classified as an L6. The mineral composition and some significant ratios of elements are reported, and the reasons for assignment to the L group and to petrologic type 6 are explained. The analysis suggests that maskelynite of oligoclase composition was formed by solid-state shock transformation of previously existing well-crystallized plagioclase at estimated shock pressures of about 250-350 kbar. This finding indicates that recrystallization (formation of well-crystallized oligoclase) preceded shock transformation formation of the maskelynite.

Keil, K.↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Comparisons of thermochemical nonequilibrium viscous flowfield predictions for AFE vehicle

Numerical solutions for the Aeroassist Flight Experiment vehicle were obtained from three methods at a trajectory point corresponding to the maximum aerodynamic heating. The flow regime and vehicle's speed require a viscous model of laminar flow and finite-rate chemical and thermal modeling. The computational domain covers both forebody and base such that the shock layer and near wake flowfield are included. Because of differences in computational grids, methods of solution, and models of rate equations, the results are generally in poor agreement. Temperature and species concentrations are strongly affected by the physical model equations and associated parameters. The chemistry model based on 11 species is found to yield lower translational temperature profile near stagnation than those from a seven-species model. The vibrational temperature varies according to the modeling details. All solutions indicate that strong neutral and/or molecular dissociation and weak ionization take place at the forebody and vibrational freezing is present in the afterbody expansion region where the vibrational temperature is higher than the translational temperature. Some forms of shear layer emanating from the aerobrake skirt coalesce in the region of reversed flow behind the vehicle.

Tam, Luen T.↗
NASA NTRS2002

Oceanic Impact: Mechanisms and Environmental Perturbations

The contents include the following: Oceanic impacts-a growing field of fundamental geoscience. Shock metamorphism on the ocean floor (numerical simulations). Numerical modeling of impact-induced modifications of the deep-sea floor. Computer modelling of the water resurge at a marine impact: the Lockne crater, Sweden. Experimental investigation of the role of water in impact vaporization chemistry. Calcareous plankton stratigraphy around the Pliocene Eltanin asteroid impact area (SE Pacific): documentation and application for geological and paleoceanographic reconstruction. Composition of impact melt debris from the Eltanin impact strewn field, Bellingshausen Sea. Iridium concentrations and abundances of meteoritic ejecta from the Eltanin impact in sediment cores from Polarstern expedition ANT XII/4. Unmelted meteoritic debris collected from Eltanin ejecta in Polarstern cores from expedition ANT XII/4. Impact tsunami-Eltanin. Ancient impact structures on modern continental shelves: The Chesapeake Bay, Montagnais, and Toms Canyon craters, Atlantic margin of North America. The Mjolnir marine impact crater porosity anomaly. Kardla (Hiiu-maa Island, Estonia) - the buried and well-preserved Ordovician marine impact structure. Long-term effect of the Kardla crater (Hiiu-maa, Estonia) on Late Ordovician carbonate sedimentation. The middle Devonian Kaluga impact crater (Russia): new interpretation of marine setting.

Gersonde, Rainer↗
DOE OSTI2023

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

An improved PNS scheme for predicting complex three-dimensional hypersonic flows

Upwinding is incorporated into a numerical technique for predicting hypersonic viscous flows over lifting configurations at moderate angles of attack. A general real-gas flux-vector-splitting technique based on Van Leer's (1982) approach is employed to model upwinding, and three techniques are examined for flux-vector differencing. The three methods are evaluated by applying them to an axisymmetric configuration with a 10-deg afterbody flare. The results indicate that an oscillation-free shock front can be described by using first-order full upwinding across the embedded shock and central-differencing for the other zones. This combined approach is found to be highly convergent for the near-wall region, and its performance is examined for predicting a Mach 15 flow over a finned missile. Attention is given to the effects of gas chemistry which can significantly affect the flows over the missile configurations.

Bhutta, Bilal A.↗
DOE OSTI2020

Structural Transformation and Chemical Stability of a Shock-Compressed Insensitive High Explosive Single Crystal: Time-Resolved Raman Spectroscopy

Despite the considerable interest in insensitive high explosives (IHE) as a safer alternative to conventional high explosives, a good understanding of the low sensitivity of IHEs to shock initiation is lacking. In particular, real-time measurements to directly probe the molecular-level response of shock compressed IHE single crystals constitute an important need. To address this need, plate impact experiments were conducted to determine time-resolved changes in the Raman spectra of 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shock compressed up to 20 GPa longitudinal stress. The Raman measurements examined vibrational frequencies from 800 cm-1 to 1500 cm-1 with 15 nanosecond time resolution and were conducted at several peak stresses. At 4 – 6 GPa, two new Raman peaks appeared, in addition to the original peaks, consistent with onset of the α' – ε structural transformation reported previously in static compression work. The measured spectra indicated completion of the transformation at 10 GPa. Raman data to 20 GPa showed neither additional transformations nor any indication of chemical decomposition. This finding, though consistent with recent continuum measurements, is in marked contrast to the chemical decomposition observed at lower stresses in shock compressed conventional high explosive single crystals. Our Raman results support the previous suggestion that strengthening of intra- and inter-molecular bonds, due to the α' – ε structural transformation, plays a significant role in the insensitivity of FOX-7 single crystals to shock initiation. In conclusion, the present work, in conjunction with previous static compression studies, provides the first experimental insight into the molecular-level response of a shock compressed IHE single crystal and can serve as a bench mark for theoretical studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗
DOE OSTI2022

Fast imaging of laser ablation of cerium: Dependence of plasma plume dynamics and structure on native oxide growth

We have imaged the evolution of plasma following nanosecond pulsed laser ablation of cerium in air with sub-10 ns resolution. We find the ablated plasma plume expansion and shape to be strongly dependent on the time after mechanical polishing of the sample surface. After an initial shock front common to both new and aged samples (velocities up to 16 km s ₋1 ), we observe a second front that is relatively localized to the surface vicinity of newly polished samples, but in the case of older samples it moves outward at velocities of up to 4 km s ₋1 . We attribute this behavior to the formation and growth of a native oxide layer on the order of hours after polishing. These results demonstrate that plasma imaging can be used to diagnose the presence of, and evaluate the extent of, thin surface oxides or other compounds that form soon after exposure to some atmosphere. These findings also highlight the need to explicitly consider the presence of such overlayers when modeling laser/metal interaction, when making measurements of vapor phase chemistry, and when analyzing the composition, phase, and morphology of solid particulates formed after the ablation of reactive metals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Forebody and base region real gas flow in severe planetary entry by a factored implicit numerical method. II - Equilibrium reactive gas

The factored-implicit, gas-dynamic algorithm has been adapted to the numerical simulation of equilibrium reactive flows. Changes required in the perfect gas version of the algorithm are developed, and the method of coupling gas-dynamic and chemistry variables is discussed. A flow-field solution that approximates a Jovian entry case was obtained by this method and compared with the same solution obtained by HYVIS, a computer program much used for the study of planetary entry. Comparison of surface pressure distribution and stagnation line shock-layer profiles indicates that the two solutions agree well.

Davy, W. C.↗
NASA NTRS1992

Scramjet combustor and nozzle computations

The SHIP3D PNS code is used to conduct a 3D parametric study of a scramjet combustor and nozzle configuration at flight Mach numbers of 8, 10, 12, 16, and 20. The study demonstrated the powerful computational strategy of using a highly efficient, modified PNS code and a nonkinetic chemistry model to analyze test data and conduct parametric studies of scramjet combustors and nozzles. The mixing efficiency was found to decrease by 0.2 when the flight Mach number increased from 8 to 12 and showed little change at higher Mach numbers. The protruding ramp injector generated an oblique shock that impinged on the cowl. Ramp-slot injection produced better mixing, both because of the interaction of this shock with the slot fuel, and because of better fuel distribution in the duct. The former effect was enhanced when the ramp angle was increased. Mixing was also enhanced when the ramp fuel injection angle and sweep angle were increased.

Kamath, Pradeep S.↗
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↗