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,603 records · Page 89

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser-plasma instabilities in long scale-length plasmas relevant to shock-ignition

We present 2D particle-in-cell simulations of laser-plasma instabilities (LPI) performed using conditions relevant to experiments on the OMEGA and NIF laser facilities. The laser intensity used is 2×10 15 W cm -2 , which allows comparison with previous experiments and simulations investigating shock ignition, though this is a relatively low intensity for the scheme. We find that the large convective gain of the instabilities leads to dynamics that are largely controlled by pump depletion and to scattering or absorption by LPIs predominantly occurring away from the quarter-critical density. The two-plasmon decay (TPD) instability is dominant for OMEGA-type conditions and remains important at the NIF-scale. In both cases, most absorption from TPD occurs near its Landau cutoff density, and we examine its nonlinear dynamics in this region, focusing on the ion-acoustic waves responsible for saturation. For NIF-type conditions, stimulated Raman scattering (SRS) and stimulated Brillouin scattering divert significant fractions of laser energy, and we find that SRS backscatter can occur with high gain at low densities due to kinetic effects. An in-depth analysis of hot-electron production is performed to identify the main sources and their characteristic electron temperature. For both simulations presented, the overall hot-electron temperature lies between 30 and 35 keV. This is relatively low and will likely increase shock pressure; however, the distributions also contain a significant number of high-energy (E k >100 keV) electrons that would likely cause unacceptable preheat. We, therefore, suggest some strategies that may be utilized to minimize this high-energy component. I. INTROD

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GnarlyX: Eulerian multi-material hydrodynamics coupled to equation of state and hyperelastic, plastic constitutive models

GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.

42 ENGINEERING↗

The near-earth cross-tail current sheet - Detailed ISEE 1 and 2 case studies

Three near geomagnetic tail current sheet crossings of the ISEE 1 and 2 satellites, on April 5, 1979, are examined in detail. All are assoicated with the passage of an interplanetary shock and the region of variable solar wind pressure behind it. The general geometry of field reversing current sheets is discussed, and this geometry is examined for the cases studied, by using the ISEE 1 and 2 coorbiting satellite data sets. A new technique is employed which removes the effects of a variable sheet normal velocity for the first time. This allows us to calculate firm upper bounds on the current sheet thicknesses, and by utilizing certain physically motivated assumptions, determine the most probable actual sheet thicknesses, and inclinations of the field lines within these sheets. Current density profiles derived with this technique show the main cross-tail current sheet to be a structure that is many thermal ion-gyroradii thick and which is sometimes imbedded in a region that is three or more times thicker and contains much smaller current densities. These profiles also exhibit a considerable amount of fine structure in the sheet which appears as narrow peaks in the current density distributions. Possible explanations for these structures, and for the overall sheet structure itself, are examined.

Mccomas, D. J.↗

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

Lunar and martian meteorite delivery services

Launch mechanisms for lunar and martian meteorites have been investigated, by integrating physical modeling constraints, geochemical cosmic-ray exposure (CRE) constraints, and petrologic constraints. The potential source region for lunar meteorites is remarkably small compared to the final crater volume. CRE constraints indicate that most launches start at depths of less than or equal to 3.2 m, and cratering theory implies derivation of suitably accelerated objects from a subvolume with diameter only about 0.3 x the final crater diameter. The shallow depth provenance is probably related to shock-wave interference, enhanced by the lunar regolith's extremely low compressional wave velocity. CRE constraints alone imply that four to five separate launch events are represented among the eight well-studied lunar meteorites. Most of the lunar meteorites are regolith breccias, which tend to show only limited compositional diversity within any kilometer-scale region of the Moon. Several others are polymict breccias, which also show relatively subdued compositional diversity, compared to igneous rocks. The observed diversity among these samples in terms of abundances of mare basalt and KREEP, and in Mg/(Mg + Fe) ratio, implies that among eight well-studied lunar meteorites only two potential source craters pairings are plausible: between Asuka-881757 + Y-793169 (most probable) and between Y-793274 + EET875721. Altogether, these eight lunar meteorites apparently represent at least six separate source craters, including three in the past 10(exp 5) years and five in the past 10(exp 6) years. CRE constraints imply that SNC meteorites are launched from systematically greater than lunar meteorites. SNCs are also systematically bigger, and all nine well-studied SNCs are uncommonly young (by martian standards) mafic igneous rocks. Comparison between Viking and Apollo results reveals that rocks the size of common meteorites are remarkably scarce in the martian regolith, probably due to pervasive weathering. A plausible explanation for these trends is that most old, small, and shallow rocks on Mars have been weakened by a two-stage process of brecciation followed by pervasive weathering, to the point where they seldom survive the stresses of spallation off the planet. The scarcity of source-crater pairing among the lunar meteorites implies that these objects can be launched from craters much smaller than previously estimated and tends to support the suggestion of Rabinowitz (1993) that present-day cratering rates for the Earth-Moon region may be higher than previously estimated.

Warren, Paul H.↗

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗

Shock initiation of cyclotol (75/25) at both ambient temperature and 70 °C

Here, six shock initiation experiments have been carried out to study the shock sensitivity and develop the unreacted Hugoniot for the Cyclotol composition 75/25 %wt RDX/TNT. The experiments were carried out both at ambient temperature, and elevated to 70 °C, close to the TNT melt temperature of 80 °C. Two sets of three experiments were fired at nominal pressures of 4.1, 6.0 and 7.7 GPa. Comparison has been made to the previously studied (80/20 %wt RDX/TNT) composition and analysis shows that the five percent decrease in RDX content has resulted in a significant increase in shock sensitivity. The increase in initial temperature has yielded an unexpectedly small increase in shock sensitivity, particularly relative to the effects of the compositional change. Wave profiles will be presented, along with a discussion of the unreacted Hugoniot and relative sensitivities presented in the Pop plot.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bayesian inferences of electrical current delivered to shocked transmission lines

Small radius Bdot measurements of electrical current delivered by pulsed power devices are routinely compromised by electrode/convolute plasmas endemic to multiterawatt transmission lines. Inferences of delivered current unaffected by these issues have recently been obtained by numerically optimizing consistency between model predictions and local experimental velocimetry data, but these are only unique for shockless velocity profiles. In the paper we describe a more general Bayesian method capable of inferring current despite the presence of shocks. Additionally, we describe uncertainty estimates and use of the technique on experimental data. This technique is the first to provide uncertainty estimates on the full current trace delivered to an inertial confinement fusion target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface flow visualization requirements for testing in NTF

Surface flow visualization is required to be correlated at all the facilities to allow the following: (1) an understanding of the chordwise and spanwise extent of laminar flow; (2) the change in shock location for various trip configurations; (3) the effectiveness of the chosen boundary layer trip (since its specifications will change as a function of Reynolds number); and (4) a comparison of the separation patterns at the buffet conditions. A requirement to tunnel-to-tunnel correlation testing involving the NTF is pointed out. A specification for a surface flow visualization system to be used in the NTF is defined. Recognition of the special limitations in the NTF including physical and visual accessibility, high operating cost, flow contamination requirements, as well as the need for on-line review of the results to develop the final trip configurations in a timely manner, leads to the listed requirements.

Bengelink, R. L.↗

A review of high-speed, convective, heat-transfer computation methods

The objective of this report is to provide useful engineering formulations and to instill a modest degree of physical understanding of the phenomena governing convective aerodynamic heating at high flight speeds. Some physical insight is not only essential to the application of the information presented here, but also to the effective use of computer codes which may be available to the reader. A discussion is given of cold-wall, laminar boundary layer heating. A brief presentation of the complex boundary layer transition phenomenon follows. Next, cold-wall turbulent boundary layer heating is discussed. This topic is followed by a brief coverage of separated flow-region and shock-interaction heating. A review of heat protection methods follows, including the influence of mass addition on laminar and turbulent boundary layers. Also discussed are a discussion of finite-difference computer codes and a comparison of some results from these codes. An extensive list of references is also provided from sources such as the various AIAA journals and NASA reports which are available in the open literature.

Tauber, Michael E.↗

Chemical environments of submarine hydrothermal systems

The paper synthesizes diverse information about the inorganic geochemistry of submarine hydrothermal systems, provides a description of the fundamental physical and chemical properties of these systems, and examines the implications of high-temperature, fluid-driven processes for organic synthesis. Emphasis is on a few general features, i.e., pressure, temperature, oxidation states, fluid composition, and mineral alteration, because these features will control whether organic synthesis can occur in hydrothermal systems.

Shock, Everett L.↗

GRB 120422A: a Low-Luminosity Gamma-Ray Burst Driven by a Central Engine

GRB 120422A is a low-luminosity gamma-ray burst (GRB) associated with a bright supernova, which distinguishesitself by its relatively short T(sub 90) (approximately 5 s) and an energetic and steep-decaying X-ray tail. We analyze the Swift BurstAlert Telescope and X-ray Telescope data and discuss the physical implications. We show that the steep declineearly in the X-ray light curve can be interpreted as the curvature tail of a late emission episode around 58-86 s,with a curved instantaneous spectrum at the end of the emission episode. Together with the main activity in thefirst 20 s and the weak emission from 40 s to 60 s, the prompt emission is variable, which points to a centralengine origin in contrast to a shock-breakout origin, which is used to interpret some other nearby low-luminosity supernova GRBs. Both the curvature effect model and interpreting the early shallow decay as the coasting externalforward shock emission in a wind medium provide a constraint on the bulk Lorentz factor to be around several.Comparing the properties ofGRB 120422A and other supernova GRBs,we find that themain criterion to distinguish engine-driven GRBs from shock-breakout GRBs is the time-averaged -ray luminosity. Engine-driven GRBs likelyhave a luminosity above approximately 10(sup 48) erg s(sup -1).

engine↗

Report on the solar physics-plasma physics workshop

The paper summarizes discussions held between solar physicists and plasma physicists on the interface between solar and plasma physics, with emphasis placed on the question of what laboratory experiments, or computer experiments, could be pursued to test proposed mechanisms involved in solar phenomena. Major areas discussed include nonthermal plasma on the sun, spectroscopic data needed in solar plasma diagnostics, types of magnetic field structures in the sun's atmosphere, the possibility of MHD phenomena involved in solar eruptive phenomena, the role of non-MHD instabilities in energy release in solar flares, particle acceleration in solar flares, shock waves in the sun's atmosphere, and mechanisms of radio emission from the sun.

Sturrock, P. A.↗

Heliospheric magnetic fields and plasmas

A survey of the existing literature on heliospheric physics, covering the period 1972-1982, is presented. Attention is given to observations and theories germane to the examination of the heliosphere as a large-scale astrophysical system that is part of the earth's environment. The literature includes data and models for magnetic sectors and the large-scale magnetic field, the large-scale plasma structure, and models and observed variations in the solar wind. Consideration is also devoted to the transient and corotating streams and shocks, the composition of the solar wind, and to MHD turbulence, waves, and discontinuities. More intensive investigations of the region near 1 AU are recommended, particularly to characterize the coronal source of the solar wind. The solar polar mission will be the first to provide radial measurements for comparisons with previous exclusively ecliptic measurements of solar activities.

Burlaga, L. F.↗

A review of high-speed, convective, heat-transfer computation methods

The objective is to provide useful engineering formulations and to instill a modest degree of physical understanding of the phenomena governing convective aerodynamic heating at high flight speeds. Some physical insight is not only essential to the application of the information presented here, but also to the effective use of computer codes which may be available to the reader. Given first is a discussion of cold-wall, laminar boundary layer heating. A brief presentation of the complex boundary layer transition phenomenon follows. Next, cold-wall turbulent boundary layer heating is discussed. This topic is followed by a brief coverage of separated flow-region and shock-interaction heating. A review of heat protection methods follows, including the influence of mass addition on laminar and turbulent boundary layers. Next is a discussion of finite-difference computer codes and a comparison of some results from these codes. An extensive list of references is also provided from sources such as the various AIAA journals and NASA reports which are available in the open literature.

Tauber, Michael E.↗

Cratering in glasses impacted by debris or micrometeorites

Mechanical strength measurements on five glasses and one glass-ceramic exposed on LDEF revealed no damage exceeding experimental limits of error. The measurement technique subjected less than 5 percent of the sample surface area to stresses above 90 percent of the failure strength. Seven micrometeorite or space debris impacts occurred at locations which were not in that portion of the sample subjected to greater than 90 percent of the applied stress. As a result of this, the impact events on the sample were not detected in the mechanical strength measurements. The physical form and structure of the impact sites was carefully examined to determine the influence of those events upon stress concentration associated with the impact and the resulting mechanical strength. The size of the impact site, insofar as it determines flaw size for fracture purposes, was examined. Surface topography of the impacts reveals that six of the seven sites display impact melting. The classical melt crater structure is surrounded by a zone of fractured glass. Residual stresses arising from shock compression and from cooling of the fused zone cannot be included in the fracture mechanics analyses based on simple flaw size measurements. Strategies for refining estimates of mechanical strength degradation by impact events are presented.

Wiedlocher, David E.↗