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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

arcjetCV: A New Machine Learning Application for Extracting Time-Resolved Recession Measurements From Arc Jet Test Videos

Arc jet Computer Vision (ArcjetCV) is a software application built to automate analysis of arc jet ground test video footage. This includes tracking material recession and the shock-material standoff distance. This provides a new capability to resolve and validate new physics associated with non-linear processes. This is an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcjetCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

machine learning↗
DOE OSTI2020

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗
DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Stability and Characteristics of Lower-hybrid Drift Waves: Dependence on Electron Beta and Cross-field Relative Drift

Lower-hybrid drift waves (LHDWs) are frequently observed microinstabilities in both space and laboratory plasmas. Despite decades of study, the relationship between electrostatic (ES-LHDW) and electromagnetic (EM-LHDW) variants and the plasma parameters controlling their stability remains unclear. Here, we systematically examine LHDW behavior by solving the local linear dispersion relation over a wide range of plasma and field conditions. Our results demonstrate that ES-LHDWs and EM-LHDWs are not distinct modes but rather two different regimes of the same drift wave whose character evolves smoothly with electron beta (β e ) and the cross-field electron drift velocity relative to ions, normalized to the ion sound speed (u 0x /C s ). The nature of the waves changes from electrostatic to electromagnetic when β e increases. Growth rates increase with u 0x /C s but decrease with β e , while the most unstable wavelength remains nearly universal, with kρ e ∼ 0.8 (k is the magnitude of the wave vector and ρ e is the electron gyroradius). We further present quasi-linear estimates of nonlinear saturation properties, including energy partition among electric fields, magnetic fields, and particle kinetic responses. We show that ES-LHDWs reach higher electric-field saturation amplitudes, whereas EM-LHDWs generate strong magnetic perturbations and parallel electric fields that may enable efficient particle heating. Comparisons with the classical model reveal that retaining electromagnetic effects is essential for accurate predictions of frequency, growth rate, and the propagation angle. These findings provide a unified framework for understanding LHDWs across diverse collisionless plasma environments, including current sheets of magnetic reconnection, shear layers, collisionless shocks, and boundary regions.

Solar coronal waves↗
NASA NTRS1966

Discontinuities in the solar wind.

Discontinuities in plasma parameters of solar wind satisfy conditions of Mhd and can be treated as shock waves and stationary contact surfaces

SHOCK WAVE↗
NASA NTRS1972

On Greenstadt's binary index criterion

The Greenstadt-type criterion for distinguishing between laminar and pulsating (i.e., irregular) bow shock crossings can be expressed as a critical angle between the shock normal and the magnetic field. When the field is more tangent to the shock plane than this critical angle, the shock is laminar. More importantly, the requirement of Galilean invariance of the physics underlying such a criterion reveals that the Greenstadt p index must vary inversely as the solar wind velcoity.

Northrop, T. G.↗
NASA NTRS1973

On Greenstadt's binary index criterion.

The Greenstadt-type criterion for distinguishing between laminar and 'pulsating' (i.e., irregular) bow shock crossings can be expressed as a critical angle between the shock normal and the magnetic field. When the field is more tangent to the shock plane than this critical angle, the shock is laminar. More importantly, the requirement of Galilean invariance of the physics underlying such a criterion reveals that the Greenstadt p index must vary inversely as the solar wind velocity.

Northrop, T. G.↗
NASA NTRS1991

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗
NASA NTRS1995

Physics of magnetospheric boundary layers

This final report was concerned with the ideas that: (1) magnetospheric boundary layers link disparate regions of the magnetosphere-solar wind system together; and (2) global behavior of the magnetosphere can be understood only by understanding its internal linking mechanisms and those with the solar wind. The research project involved simultaneous research on the global-, meso-, and micro-scale physics of the magnetosphere and its boundary layers, which included the bow shock, the magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical, and simulation projects were performed on these subjects, as well as comparisons of theoretical results with observational data. Other related activity included in the research included: (1) prediction of geomagnetic activity; (2) global MHD (magnetohydrodynamic) simulations; (3) Alfven resonance heating; and (4) Critical Ionization Velocity (CIV) effect. In the appendixes are list of personnel involved, list of papers published; and reprints or photocopies of papers produced for this report.

Cairns, Iver H.↗
NASA NTRS1995

Experimental and Computational Investigation of the Tip Clearance Flow in a Transonic Axial Compressor Rotor

Experimental and computational techniques are used to investigate tip clearance flows in a transonic axial compressor rotor at design and part speed conditions. Laser anemometer data acquired in the endwall region are presented for operating conditions near peak efficiency and near stall at 100% design speed and at near peak efficiency at 60% design speed. The role of the passage shock/leakage vortex interaction in generating endwall blockage is discussed. As a result of the shock/vortex interaction at design speed, the radial influence of the tip clearance flow extends to 20 times the physical tip clearance height. At part speed, in the absence of the shock, the radial extent is only 5 times the tip clearance height. Both measurements and analysis indicate that under part-speed operating conditions a second vortex, which does not originate from the tip leakage flow, forms in the endwall region within the blade passage and exits the passage near midpitch. Mixing of the leakage vortex with primary flow downstream of the rotor at both design and part speed conditions is also discussed.

Suder, Kenneth L.↗
NASA NTRS1995

Pair cascades in extragalactic jets. 1: Gamma rays

A model of the approximately 0.1-10 GeV gamma-ray jets observed by the EGRET instrument on the Compton Gamma Ray Observatory (CGRO) is developed. It is shown that the soft X-ray background in an active galactic nuclei (AGN) contributes an opacity to pair production and that a gamma-ray photosphere or 'gamma-sphere' can be defined whose radius increases with gamma-ray energy E(sub gamma). It is proposed that the observed gamma-ray emission is due to inverse Compton scattering of the ambient soft X-rays by relativistic pairs accelerated in situ by shock fronts in a relativistic jet. For a wide range of assumed physical conditions, the emission at a given E(sub gamma) originates from near the associated gamma-spheres; emission from below the gamma-sphere initiates a cascade down to the energy where the gamma-rays can escape freely. In this model, the slope of the emergent gamma-ray spectrum is determined by the scattered, soft X-ray spectrum and the variation of the particle acceleration rate with jet radius. In general it is expected that the variation in the gamma-ray flux will be either slower or later at higher energy. It is also shown that the efficiency of conversion of energy from injected high-energy pairs to 0.1-10 GeV gamma-rays is typically high so that the models are radiatively efficient. It is argued that the observed gamma-ray jets are likely to be particle-dominated, though magnetically confined. The gamma-ray spectrum should continue down to an energy approximately 5 MeV emitted from an annihilation radius within which the pair content of the jet is limited by annihilation. This is probably the site of the beamed hard X-ray emission. It is speculated that the relativistic jets associated with radio-loud AGNs are powered electromagnetically by a spinning black hole and that they are collimated by an encircling MHD wind leaving the accretion disk at a slower speed. Powerful FR2 radio sources are formed when the hole spins rapidly and the relativistic core accelerates the MHD sheath; low-power FR1 sources ensue when the opposite occurs. Finally, it is suggested that the key factor which determines whether or not a given active nucleus can form a jet and a radio to gamma-ray nonthermal continuum is the central density of mass-losing stars which, when large, precludes the formation of a super-Alfvenic, collimating wind.

Blandford, R. D.↗
NASA NTRS1999

Automated Fluid Feature Extraction from Transient Simulations

In the past, feature extraction and identification were interesting concepts, but not required to understand the underlying physics of a steady flow field. This is because the results of the more traditional tools like iso-surfaces, cuts and streamlines were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of much interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snap-shot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments like pV3). And methods must be developed to abstract the feature and display it in a manner that physically makes sense. The following is a list of the important physical phenomena found in transient (and steady-state) fluid flow: (1) Shocks, (2) Vortex cores, (3) Regions of recirculation, (4) Boundary layers, (5) Wakes. Three papers and an initial specification for the (The Fluid eXtraction tool kit) FX Programmer's guide were included. The papers, submitted to the AIAA Computational Fluid Dynamics Conference, are entitled : (1) Using Residence Time for the Extraction of Recirculation Regions, (2) Shock Detection from Computational Fluid Dynamics results and (3) On the Velocity Gradient Tensor and Fluid Feature Extraction.

Haimes, Robert↗
NASA NTRS2016

Flash Spectroscopy: Emission Lines From the Ionized Circumstellar Material Around 10-Day-Old Type II Supernovae

Supernovae (SNe) embedded in dense circumstellar material (CSM) may show prominent emission lines in their early-time spectra (< or =10 days after the explosion), owing to recombination of the CSM ionized by the shock-breakout flash. From such spectra ("flash spectroscopy"), we can measure various physical properties of the CSM, as well as the mass-loss rate of the progenitor during the year prior to its explosion. Searching through the Palomar Transient Factory (PTF and iPTF) SN spectroscopy databases from 2009 through 2014, we found 12 SNe II showing flash-ionized (FI) signatures in their first spectra. All are younger than 10 days. These events constitute 14% of all 84 SNe in our sample having a spectrum within 10 days from explosion, and 18% of SNe II observed at ages <5 days, thereby setting lower limits on the fraction of FI events. We classified as "blue/featureless" (BF) those events having a first spectrum that is similar to that of a blackbody, without any emission or absorption signatures. It is possible that some BF events had FI signatures at an earlier phase than observed, or that they lack dense CSM around the progenitor. Within 2 days after explosion, 8 out of 11 SNe in our sample are either BF events or show FI signatures. Interestingly, we found that 19 out of 21 SNe brighter than an absolute magnitude M(sub R) = −18.2 belong to the FI or BF groups, and that all FI events peaked above M(sub R) = −17.6 mag, significantly brighter than average SNe II.

supernovae: general↗
NASA NTRS2020

Electron Energy Partition Across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

Solar wind↗