Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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,153 records · Page 64

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗

Improved design of a high response slotted plate overboard bypass valve for supersonic inlets

The electrohydraulically actuated slotted-plate bypass valve used to control the position of the normal shock during wind-tunnel investigations of supersonic inlets also has proven to be a valuable device for determining inlet dynamics and creating airflow disturbances. Operation of previous valves at high frequencies (to 100 Hz) for extended testing has resulted in numerous failures. An improved bypass-valve design is presented which increases the cyclic tolerance of the device considerably over past designs. The use of dynamic limit criteria to obtain an optimum actuator-piston size results in a frequency response which is flat within + or - 3 decibels to 120 Hz for a peak-to-peak variation of 20 percent of full area.

Webb, J. A., Jr.↗

The University of Iowa Helios solar wind plasma wave experiment /E 5a/

This document describes the University of Iowa solar wind plasma wave experiment for the Helios missions (Experiment 5a). The objective of this experiment is the investigation of naturally occurring plasma instabilities and electromagnetic waves in the solar wind. To carry out this investigation, the experiment consists primarily of a 16-channel spectrum analyzer connected to the electric field antennas. The spectrum analyzer covers the frequency range from 20 Hz to 200 kHz and has an amplitude dynamic range which extends from .3 microvolts/m to 30 mV/m per channel. This spectrum analyzer, the antenna potential measurements, the shock alarm system and the supporting electronics are discussed in detail.

Gurnett, D. A.↗

Solar flare effects on the zodiacal light

Results are reported for an observational and theoretical study of possible solar-flare effects on the zodiacal light. A seven-year search for flare-induced brightness variations in the light is described which discovered no significant changes in either the level of brightness from night to night or the shape of the main zodiacal-light cone at elongations greater than 25 deg from the sun. Several possible dynamical effects of solar flares on the zodiacal dust cloud are investigated theoretically, including increased ion drag from a passing shock wave, possible rotational bursting of dust particles as a result of ion impacts, and fluorescence of the dust particles. It is found that the kinetic energy carried by shock-front ions is not sufficient to cause the dust to fluoresce, that increased ion drag may be significant only at heliocentric distances less than 0.1 AU, and that particles larger than 10 microns cannot attain bursting speeds at the appropriate distances. The question is considered of whether changes in the brightness of the polarized component as a result of variations in electron number density can be observed from earth or the Helios solar probe.

Misconi, N. Y.↗

Atmospheric effects on inlets for supersonic cruise aircraft

Mixed-compression inlet dynamic behavior in the vicinity of unstart, was simulated and analyzed to investigate time response of an inlet's normal shock to independent disturbances in ambient temperature and pressure and relative velocity (longitudinal gust), with and without inlet controls active. The results indicate that atmospheric disturbances may be more important than internal disturbances in setting inlet controls requirements because they are usually not anticipated and because normal shock response to rapid atmospheric disturbances is not attenuated by the inlet, as it is for engine induced disturbances. However, before inlet control requirements can be fully assessed, more statistics on extreme atmospheric disturbances are needed.

Cole, G. L.↗

Grain motions in the solar nebula

Isotopic analyses of meteorites suggest the possibility that some interaction between supernova ejecta and grains occurred in the solar nebula. In particular, the dynamics of grain motions in the solar nebula can explain the observed mixing of nucleosynthetic components. The effect of a shock wave on the motions of grains is examined. On the basis of calculations, it is estimated that if grains carried the isotopic anomalies investigated by Lee, Papanastassoiu, and Wasserburg (1978), then those grains could be no bigger than 0.0002 cm in size. A scenario is suggested in which the sluggishness of grains provides a natural way to concentrate and mix the nucleosynthetic components carried by grains in the ejecta and in the solar nebula.

Margolis, S. H.↗

A shock capturing technique for hypersonic, chemically relaxing flows

A fully coupled, shock capturing technique is presented for chemically reacting flows at high Mach numbers. The technique makes use of a total variation diminishing (TVD) dissipation operator which results in sharp, crisp shocks. The eigenvalues and eigenvectors of the fully coupled system, which includes species conversion equations in addition to the gas dynamics equations, are analytically derived for a general reacting gas. Species production terms for a model dissociating gas are introduced and are included in the algorithm. The convective terms are solved using a first-order TVD scheme while the source terms are solved using a fourth-order Runge-Kutta scheme to enhance stability. Results from one-dimensional numerical experiments are shown for a two species and a three species gas.

Eberhardt, S.↗

Impact Cratering Calculations

Many Martian craters are surrounded by ejecta blankets which appear to have been fluidized forming lobate and layered deposits terminated by one or more continuous distal scarps, or ramparts. One of the first hypotheses for the formation of so-called rampart ejecta features was shock-melting of subsurface ice, entrainment of liquid water into the ejecta blanket, and subsequent fluidized flow. Our work quantifies this concept. Rampart ejecta found on all but the youngest volcanic and polar regions, and the different rampart ejecta morphologies are correlated with crater size and terrain. In addition, the minimum diameter of craters with rampart features decreases with increasing latitude indicating that ice laden crust resides closer to the surface as one goes poleward on Mars. Our second goal in was to determine what strength model(s) reproduce the faults and complex features found in large scale gravity driven craters. Collapse features found in large scale craters require that the rock strength weaken as a result of the shock processing of rock and the later cratering shear flows. In addition to the presence of molten silicate in the intensely shocked region, the presence of water, either ambient, or the result of shock melting of ice weakens rock. There are several other mechanisms for the reduction of strength in geologic materials including dynamic tensile and shear induced fracturing. Fracturing is a mechanism for large reductions in strength. We found that by incorporating damage into the models that we could in a single integrated impact calculation, starting in the atmosphere produce final crater profiles having the major features found in the field measurements (central uplifts, inner ring, terracing and faulting). This was accomplished with undamaged surface strengths (0.1 GPa) and in depth strengths (1.0 GPa).

Ahrens, Thomas J.↗

Rocket Engine Nozzle Side Load Transient Analysis Methodology: A Practical Approach

At the sea level, a phenomenon common with all rocket engines, especially for a highly over-expanded nozzle, during ignition and shutdown is that of flow separation as the plume fills and empties the nozzle, Since the flow will be separated randomly. it will generate side loads, i.e. non-axial forces. Since rocket engines are designed to produce axial thrust to power the vehicles, it is not desirable to be excited by non-axial input forcing functions, In the past, several engine failures were attributed to side loads. During the development stage, in order to design/size the rocket engine components and to reduce the risks, the local dynamic environments as well as dynamic interface loads have to be defined. The methodology developed here is the way to determine the peak loads and shock environments for new engine components. In the past it is not feasible to predict the shock environments, e.g. shock response spectra, from one engine to the other, because it is not scaleable. Therefore, the problem has been resolved and the shock environments can be defined in the early stage of new engine development. Additional information is included in the original extended abstract.

Shi, John J.↗

Overset Techniques for Hypersonic Multibody Configurations with the DPLR Solver

Three unit problems in shock-shock/shock-boundary layer interactions are considered in the evaluation overset techniques with the Data Parallel Line Relaxation (DPLR) computational fluid dynamics solver, a three dimensional Navier-Stokes solver . The unit problems considered are those of two stacked hemispherical cylinders (of different diameters and lengths, and at various orientations relative to each other or relative to the nozzle axis) tested in a hypersonic wind tunnel. These problems are taken as representative of a Two-Stage-To-Orbit design. The objective of the present presentation would be to discuss the techniques used to develop suitable overset grid systems and then evaluate their respective solutions by comparing to corresponding point matched grid solutions and experimental data. Both successful and unsuccessful techniques would be discussed. All solutions would be calculated using the DPLR solver and SUGGAR will be used to develop the domain connectivity information.

Hyatt, Andrew James↗

Best Practices for Unstructured Grid Shock Fitting

Unstructured grid solvers have well-known issues predicting surface heat fluxes when strong shocks are present. Various efforts have been made to address the underlying numerical issues that cause the erroneous predictions. The present work addresses some of the shortcomings of unstructured grid solvers, not by addressing the numerics, but by applying structured grid best practices to unstructured grids. A methodology for robust shock detection and shock fitting is outlined and applied to production relevant cases. Results achieved by using the Loci-CHEM Computational Fluid Dynamics solver are provided.

McCloud, Peter L.↗

NASA Supersonics Research — Leading Toward Quiet Supersonic Flight!!

- Brief history of supersonic flight - Shock waves and sonic booms - Sonic boom pressure signatures and loudness - Tools for studying quiet supersonic flight - CFD (Computational Fluid Dynamics) - Wind tunnels - Flight test - Quesst Mission - X-59 airplane - Mission profile - Flight simulations - Wind tunnel data validations - Video: Supersonic flight research leading to Quesst

sonic boom↗

NASA Supersonics Research and the X-59 — How They Will Benefit You!!

- Brief history of supersonic flight - Shock waves and sonic booms - Sonic boom pressure signatures and loudness - Tools for studying quiet supersonic flight - CFD (Computational Fluid Dynamics) - Wind tunnels - Flight test - Quesst Mission - X-59 airplane - Mission profile - Flight simulations - Wind tunnel data validations - Video: Supersonic flight research leading to Quesst

sonic boom↗

Nasa Supersonics Research and the X-59 —To Enable Future Supersonic Travel!

- Brief history of supersonic flight - Shock waves and sonic booms - Sonic boom pressure signatures and loudness - Tools for studying quiet supersonic flight - CFD (Computational Fluid Dynamics) - Wind tunnels - Flight test - Quesst Mission - X-59 airplane - Mission profile - Flight simulations - Wind tunnel data validations - Video: Supersonic flight research leading to Quesst

sonic boom↗

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

New Constraints on the Melting Temperature and Phase Stability of Shocked Iron up to 270 GPa Probed by Ultrafast X-Ray Absorption Spectroscopy

Studying the properties and phase diagram of iron at high-pressure and high-temperature conditions has relevant implications for Earth’s inner structure and dynamics and the temperature of the inner core boundary (ICB) at 330 GPa. Also, a hexagonal-closed packed to body-centered cubic (bcc) phase transition has been predicted by many theoretical works but observed only in a few experiments. The recent coupling of high-power laser with advanced x-ray sources from synchrotrons allows for novel approaches to address these issues. Here, we present a study on shock compressed iron up to 270 GPa and 5800 K probed by single-pulse (100 ps FWHM) x-ray absorption spectroscopy (XAS). Based on the analysis of the XAS spectra, we provide structural identification and bulk temperature measurements along the Hugoniot up to the melting. These results rule out the predicted transition to a high-temperature bcc phase and allow one to discriminate among existing equations of state models and melting curves. In particular, we report the first bulk temperature measurement in shock compressed iron on the melting plateau at 240(20) GPa and 5345(600) K. The melting curve resulting from our work extrapolates to a temperature of 6202(514) K at 330 GPa and represents a refined upper bound for the ICB temperature. Published by the American Physical Society 2024

Balugani, S.↗

The Solar Wind at 20-30 AU

Pioneer 10 sampled the interplanetary plasma over the range 20 to 30 astronomical units, during the period 1979-1983. The median flow speed is about 400 km/s, and at 20 AU the median density, proton temperature and dynamic pressure are, respectively, 0.025 cm-3, 10(4) K, and 6x10-11 dyne cm-2. It is shown that the average solar wind flow speed does not vary significantly with increasing heliocentric distance, and the density falls off as R-2, as predicted by simple solar wind models. The day-to-day variations in solar wind parameters are smaller at larger distance. Very large shocks however, were detected beyond 25 AU. Comparison of Pioneer 10 and 11 observations at similar distances but different phases of the solar activity cycle shows that solar wind dynamic pressure varies over a wider range during epochs of high solar activity. The variation near 20 AU is likely to be smaller at Voyager 2 Uranus encounter than observed by Pioneer 10 in the 1979-80 period.

Barnes, A.↗