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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

The Link Between Shocks, Turbulence, and Magnetic Reconnection in Collisionless Plasmas

Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence and magnetic reconnection. In particular collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These 'wavefronts' maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, 'stir' the downstream region, creating large scale disturbances such as vortices, sunward flows, and can trigger flux ropes along the magnetopause. This phenomenology closes the loop between shocks, turbulence and magnetic reconnection in ways previously unrealized. These interconnections appear generic for the collisionless plasmas typical of space, and are expected even at planar shocks, although they will also occur at curved shocks as occur at planets or around ejecta.

MAGNETOSPHERE↗
NASA NTRSDate not supplied

Numerical Simulations of an Axisymmetric Shock-Wave/Boundary-Layer Interaction

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed of two Mach 2.5 axisymmetric shock-wave/boundary-layer interactions (SWBLIs) based off of the NASA Turbulent Computational Fluid Dynamics Validation Experiment (TCFDVE). The goals of the simulations were to (1) determine how well RANS can predict the flow physics of an axisymmetric SWBLI and (2) generate a new validation case for the NASA Turbulence Modeling Resource (TMR). The simulations were run with two different turbulence models and three different computational fluid dynamics (CFD) codes while looking at two different conical shock generator angles, one of which was shown to cause flow separation in the experiment. The results in this paper show that, as expected, the CFD simulations agreed well with each other and the experiment when SWBLI-induced separation was not present. When separation was present, differences in the predicted separation characteristics were noted. These included a stronger interaction and subsequently longer separation extent for the Menter Shear-Stress Transport k-ω (SST) turbulence model simulations compared to both the Spalart-Allmaras (SA) turbulence model simulations and the experiment.

Computational Fluid Dynamics↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗
DOE OSTI2020

Detectability of standing accretion shock instabilities activity in supernova neutrino signals

We introduce a novel methodology for establishing the presence of standing accretion shock instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice- based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test- statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda et al. [Astrophys. J. 851, 62 (2017).] as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to ~6 kpc for IceCube and up to ~3 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1962

ENVIRONMENTAL TESTING

Dynamic environments and environmental testing - shock and vibration testing

ENVIRONMENTAL TESTING↗
NASA NTRS1972

Review of delta wing space shuttle vehicle dynamics

The unsteady aerodynamics of the delta planform, high cross range, shuttle orbiter were investigated. It has been found that these vehicles are subject to five unsteady flow phenomena that could compromise the flight dynamics. They are: (1) leeside shock induced separation, (2) sudden leading edge stall, (3) vortex burst, (4) bow shock-flap shock interaction, (5) forebody vorticity. Trajectory shaping is seen as the most powerful means of avoiding the detrimental effects of the stall phenomena. However, stall must be fixed or controlled when traversing the stall region. The other phenomena may be controlled by carefully programmed control deflections and some configuration modification. Ways to alter the occurrence of the various flow conditions are explored.

Reding, J. P.↗
NASA NTRS1990

Supersonic convection

Numerical simulations with high spatial resolution are used to study that the combined effects of stratification, pressure gradients, and nonadiabatic processes can lead to the formation of regions of supersonic motions near the upper thermal boundary layer. Within these regions, the dynamics is dominated by nonstationary shock structures. These form near the downflow sites and propagate upstream along the boundary layer to the upflow regions where they weaken and eventually disappear. The shock cycle, consisting of the formation, propagation, and disappearance of shock structures, has a time scale comparable to the sound crossing time over a portion of the convective cell, giving rise to vigorous time dependence in the convection.

Cattaneo, Fausto↗
NASA NTRS1994

Radiative properties of molecular nitrogen ions produced by helium Penning ionization and argon effects

The development of hypersonic aerospace vehicles requires a better understanding on the thermal and chemical nonequilibrium kinetics of participating species in shock layers. The computational fluid dynamic (CFD) codes developed for such flowfields overestimate the radiation in the spectral region of 300 - 600 nm. A speculation for this overestimation is that inclusion of Ar, CO2, and H2O at the upper atmosphere flight region makes a significant impact on radiative kinetics of molecular nitrogen ions. To define the effects of minority species on the radiative kinetics of N2(+), an experimental setup was made by using the helium Penning ionization. The vibrational and rotational temperature were measured by mapping the vibrational and rotational distributions of N2(+) emission with high spectroscopic resolution and absolute intensity measurements. Measured vibrational temperatures were in the range from 18,000 to 36,000 K, and rotational temperatures were in the range from 300 to 370 K. The irradiance of 391.44 nm line and rotational and vibrational temperatures were analyzed to define argon and CO2 effects on the N2(+) emission. When Ar or CO2 is injected with N2, the rotational temperature did not change. The irradiances were reduced by 34 percent and 78 percent for the 50 percent of mixture of Ar and CO2, respectively. The vibrational temperatures were increased by 24.1 percent and 82.9 percent for the 50 percent of mixture of Ar and CO2, respectively. It appears that there are no significant effects from small concentrations of Ar and CO2 at the upper atmosphere flight region.

Miller, George, III↗
NASA NTRS1994

The CO/SiO radiative instability in cool star atmospheres revisited

We revisit the formation of radiative instabilities in cool star atmospheres and compare our results with those given by Muchmore, Nuth, & Stencel. We have considered the combined influence of CO and SiO molecules and have computed models for a grid of effective temperatures and geometrical dilution factors for the stellar radiation. Our results are based on the analysis of the energy balance of gas elements with prescribed thermodynamic properties. Our results show that radiative instabilities are most likely primarily caused by CO, whereas SiO is expected to play only a minor role, except when the CO density is reduced compared to Local Thermodynamic Equilibrium (LTE) values or the CO band can be assumed to be optically thick. The onset of radiative instabilities is expected to be strongly modified when dynamic phenomena such as stochastic shocks are present. Our results provide strong evidence that dust formation can most likely occur via a radiative instability alone. Therefore, we present a revised version of the Muchmore et al. dust formation paradigm, which also considers hydrodynamic cooling. The new paradigm is particularly relevant in cases where dust is formed relatively close to the stellar photosphere.

Cuntz, M.↗
NASA NTRS2007

Alfven Waves in the Solar Wind, Magnetosheath, and Outer Magnetosphere

Alfven waves Propagating outward from the Sun are ubiquitous in the solar wind and play a major role in the solar wind-magnetosphere interaction. The passage of the waves generally occurs in the form of a series of discrete steepened discontinuities, each of which results in an abrupt change in the interplanetary magnetic field direction. Some orientations of the magnetic field permit particles energized at the Earth's bow shock to gain access to the foreshock region immediately upstream from the Earth's bow shock. The thermal pressure associated with these particles can greatly perturb solar wind plasma and magnetic field parameters shortly prior to their interaction with the Earth's bow shock and magnetosphere. The corresponding dynamic pressure variations batter the magnetosphere, driving magnetopause motion and transient compressions of the magnetospheric magnetic field. Alfven waves transmit information concerning the dynamic pressure variations applied to the magnetosphere to the ionosphere, where they generate the traveling convection vortices (TCVs) seen in high-latitude ground magnetograms. Finally, the sense of Alfvenic perturbations transmitted into the magnetosheath reverses across local noon because magnetosheath magnetic field lines drape against the magnetopause. The corresponding change in velocity perturbations must apply a weak torque to the Earth's magnetosphere.

Sibeck, D. G.↗
NASA NTRS2017

Recent Advancements in the Infrared Flow Visualization System for the NASA Ames Unitary Plan Wind Tunnels

The following details recent efforts undertaken at the NASA Ames Unitary Plan wind tunnels to design and deploy an advanced, production-level infrared (IR) flow visualization data system. Highly sensitive IR cameras, coupled with in-line image processing, have enabled the visualization of wind tunnel model surface flow features as they develop in real-time. Boundary layer transition, shock impingement, junction flow, vortex dynamics, and buffet are routinely observed in both transonic and supersonic flow regimes all without the need of dedicated ramps in test section total temperature. Successful measurements have been performed on wing-body sting mounted test articles, semi-span floor mounted aircraft models, and sting mounted launch vehicle configurations. The unique requirements of imaging in production wind tunnel testing has led to advancements in the deployment of advanced IR cameras in a harsh test environment, robust data acquisition storage and workflow, real-time image processing algorithms, and evaluation of optimal surface treatments. The addition of a multi-camera IR flow visualization data system to the Ames UPWT has demonstrated itself to be a valuable analyses tool in the study of new and old aircraft/launch vehicle aerodynamics and has provided new insight for the evaluation of computational techniques.

Infrared Thermography↗
NASA NTRS2019

Cluster Analysis of Spectroscopic Line Profiles and EUV Emission in RMHD Simulations and Observations of the Solar Atmosphere

Spatially-resolved observations from the IRIS, SDO/AIA, and other space mission and ground-based telescopes, coupled with realistic 3D RMHD simulations, are a powerful tool for analysis of processes in the solar atmosphere. To better understand the dynamical and thermodynamic properties in the simulation data and their connection to observations, it is essential to determine similarities in the behaviors of the synthesized and observed emission. However, the complexity of observational data and physical processes makes comparison of observations and modeling results difficult. In this work, we show the initial results of application of K-Means clustering (unsupervised machine learning) algorithm to two different problems: 1) recognition of the typical spectroscopic line profiles observed by IRIS during solar flares and their typical dynamic behavior; 2) recognition of shocks and heating events in synthetic AIA emission data obtained from StellarBox quiet-Sun simulations. The average silhouette width technique for the KMeans algorithm is utilized in different ways to obtain optimal numbers of clusters. We discuss application of the emission clustering to visualizations of the computational volume, understanding its evolutionary trends and behavior patterns, and inversion (reconstruction) of physical properties of the solar atmosphere from synthesizes emission data.

Sadykov, Viacheslav↗
NASA NTRS2003

Magnetohydrodynamic Jump Conditions for Oblique Relativistic Shocks with Gyrotropic Pressure

Shock jump conditions, i.e., the specification of the downstream parameters of the gas in terms of the upstream parameters, are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. This is done by combining the continuity of particle number flux and the electromagnetic boundary conditions at the shock with the magnetohydrodynamic conservation laws derived from the stress-energy tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. Examples assuming isotropic pressure illustrate how the shock compression ratio depends on the shock speed and obliquity. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically with- out additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., approximately 20%) produce significant changes in the shock compression ratio, r , at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks; such models include the modification of the fluid velocity profile due to the contribution of energetic particles to the momentum and energy fluxes.

Double, Glen P.↗
DOE OSTI2024

Shock Response of Polymethyl Methacrylate (PMMA) Under Explosive Loading

Abstract The dynamic response of polymethyl methacrylate (PMMA) is well understood for one-dimensional planar impact shocks, but limited research has been performed on the response of PMMA under spherical shock loading. In this work, the shock decay of an explosively-driven shock wave into PMMA was experimentally measured. PMMA cubes of various geometries were explosively loaded with an RP-80 detonator to produce the explosive shock wave. High-speed schlieren imaging was implemented to measure the explosively-driven shock wave velocity throughout the PMMA cubes. Photon Doppler velocimetry (PDV) was used to measure the particle velocity imparted by the shock wave at the surface of the cubes. The material shock response was studied at distances from 21.91 to 133.3 mm from the explosive source. The particle velocity history measured by PDV was compared to the wave profile visualized in the high-speed images. The shock wave pulse amplitude decreased with increased distance from the source. The conducted experiments extend the PMMA shock Hugoniot relating to the lower shock and particle velocity regime.

Torres, S. M. (ORCID:0009000287177206)↗
DOE OSTI2021

Visualization of shocked material instabilities using a fast-framing camera and XFEL four-pulse train

Many questions regarding dynamic materials could be answered by using time-resolved ultra-fast imaging techniques to characterize the physical and chemical behavior of materials in extreme conditions and their evolution on the nanosecond scale. In this work, we perform multi-frame phase-contrast imaging (PCI) of micro-voids in low density polymers under laser-driven shock compression. At the Matter in Extreme Conditions (MEC) Instrument at the Linac Coherent Light Source (LCLS), we used a train of four x-ray free electron laser (XFEL) pulses to probe the evolution of the samples. To visualize the void and shock wave interaction, here, we deployed the Icarus V2 detector to record up to four XFEL pulses, separated by 1-3 nanoseconds. In this work, we image elastic waves interacting with the micro-voids at a pressure of several GPa. Monitoring how the material’s heterogeneities, like micro-voids, dictate its response to a compressive wave is important for benchmarking the performances of inertial confinement fusion energy materials. For the first time in a single sample, we have combined an ultrafast x-ray framing camera and four XFEL pulse train to create an ultrafast movie of micro-void evolution under laser-driven shock compression. Eventually, we hope this technique will resolve the material density as it evolves dynamically under laser shock compression.

fusion↗
DOE OSTI2020

Phase Transitions in (Mg,Fe) 2 SiO 4 Olivine under Shock Compression

The response of forsterite, Mg 2 SiO 4 , under dynamic compression is of fundamental importance for understanding its phase transformations and high-pressure behavior. Here, we have carried out an in situ X-ray diffraction study of laser-shocked polycrystalline and single-crystal forsterite from 19 to 122 GPa using the Matter in Extreme Conditions end-station of the Linac Coherent Light Source. Under laser-based shock loading, forsterite does not transform to the high-pressure equilibrium assemblage of MgSiO 3 bridgmanite and MgO periclase, as has been suggested previously. Instead, we observe forsterite and forsterite III, a metastable polymorph of Mg 2 SiO 4 , coexisting in a mixed-phase region from 33 to 75 GPa for both polycrystalline and single-crystal samples. Densities inferred from X-ray diffraction are consistent with earlier gas-gun shock data. At higher stress, the response is sample-dependent. Polycrystalline samples undergo amorphization above 79 GPa. For [010]- and [001]-oriented crystals, a mixture of crystalline and amorphous material is observed to 108 GPa, whereas the [100]-oriented forsterite adopts an unknown phase at 122 GPa. The first two sharp diffraction peaks of amorphous Mg 2 SiO 4 show a similar trend with compression as those observed for MgSiO 3 in both recent static- and laser-driven shock experiments. This study provides new insight into the transformation of forsterite under nanosecond-duration shock loading. This work emphasizes the importance of formation of metastable phases along the Hugoniot and adds to evidence that the 300-K single-crystal diamond anvil cell experiments have relevance for understanding structures formed under shock compression. In particular, the metastable phase forsterite III has now been shown to form under dynamic compression from 10s to 100s of nanoseconds as well as under 300-K static compression. Upon compression to higher pressures, Mg 2 SiO 4 transforms to an amorphous phase. These results have broad relevance for understanding the behavior of silicates under dynamic compression.

36 MATERIALS SCIENCE↗