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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Limited Dynamic Triggering in the Utah Region, USA

The state of Utah, USA, experiences around 3800 catalogued earthquakes per year, highlighting that the region is seismically active and susceptible to earthquakes. Following the 2002 Denali Fault (M7.9) earthquake in Alaska, the region showed an elevated seismicity rate for 3 weeks following the passage of high amplitude surface waves, suggesting that the region may be particularly susceptible to dynamic triggering. With over 23,396 faults and each fault presenting a potential fault for triggering, we systematically search for dynamic triggering throughout the state of Utah caused by large, global earthquakes with M ≥ 7. Specifically, we analyse earthquake catalogues and all available waveform data to determine statistical increases of seismicity rate following the passage of seismic arrivals. While we find instances of dynamic triggering, our results show that these events occur sparsely in the region. In total, less than 20 per cent of the 273 main shocks that occur from 2000 to the end of 2017 show a statistical indication of dynamic triggering throughout the Utah region, highlighting that dynamic triggering is limited for stresses created by transient signals from global M ≥ 7 earthquakes, with the exception being the Denali Fault (M7.9), Alaska earthquake (i.e. an instance of significant triggering).

58 GEOSCIENCES↗
DOE OSTI2021

In Situ X-ray Diffraction of High Density Polyethylene During Dynamic Compression: Polymer Chain Compression and Phase Transformations

Polymers are widely implemented materials often used for safety and security in our daily lives. However, we do not have a thorough understanding of the phase transformations of these materials when exposed to extreme condition, high strain rate, events. To study polymers in these severe environments, we performed plate impact experiments that impart a one-dimensional shock wave, at high pressures and temperatures, to observe the materials response. Generally, velocimetry is the experimental diagnostic used to capture the movement of high velocity interfaces to then calculate the mechanical variables associated with shock compression. However, when a material undergoes a phase transition, such as between two crystalline phases or melting to decomposition during shock compression, velocimetry does not explicitly indicate what structural change occurred. At the Dynamic Compression Sector at the Advanced Photon Source gas gun driven compression has been paired with x-ray diffraction, providing the ability to study both the shock parameters and structural changes concurrently. We explored high density polyethylene (HDPE) both on its reactants Hugoniot (<10.12 GPa) and within its transition region (24.3 GPa) where HDPE is hypothesized to decompose from reactants to products. With x-ray diffraction we observed both polymer chain compression and crystal structure reorientation from orthorhombic to monoclinic above 6 GPa, as well as either a possible phase change or chemical decomposition at 24.3 GPa. With these experiments we were able to determine the crystalline and amorphous Hugoniot loci at pressures below the polymer transition region.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Wave and particle evolution downstream of quasi-perpendicular shocks

Distributions of ions heated in quasi-perpendicular bow shocks have large perpendicular temperature anisotropies that provide free energy for the growth of Alfven ion cyclotron (AIC) and mirror waves. These modes are often obsreved in the Earth's magnetosheath. Using two-dimensional hybrid simulations, we show that these waves are produced near the shock front and convected downstream rather than being produced locally downstream. The wave activity reduces the proton anisotropy to magnetosheath levels within a few tens of gyroradii of the shock but takes significantly longer to reduce the anisotropy of He(++) ions. The waves are primarily driven by proton anisotropy and the dynamics of the helium ions is controlled by the proton waves. Downstream of high Mach number shocks, mirror waves compete effectively with AIC waves. Downstream of low Mach number shocks, AIC waves dominate.

Mckean, M. E.↗
NASA NTRS1972

Feasibility study of inlet shock stability system of YF-12

The feasibility of self actuating bleed valves as a shock stabilization system in the inlet of the YF-12 is considered for vortex valves, slide valves, and poppet valves. Analytical estimation of valve performance indicates that only the slide and poppet valves located in the inlet cowl can meet the desired steady state stabilizing flows, and of the two the poppet valve is substantially faster in response to dynamic disturbances. The poppet valve is, therefore, selected as the best shock stability system for the YF-12 inlet.

Blausey, G. C.↗
NASA NTRS1995

Merged interaction regions and large-scale fluctuations observed by the Voyager 2 in the distant heliosphere

The merged interaction regions (MIRs) and large-scale fluctuations of the heliospheric magnetic field play a major role in the dynamics of the solar wind, the position and motion of the termination shock and heliopause, the triggering of radio emissions, and the modulation of cosmic rays. The structure of MIRs and large-scale fluctuations varies with distance from the sun and with solar activity. Here we compare Voyager 2 observations near the maximum of solar activity (1989 through 1991) with those during the declining phase of solar activity (1992 thorough 1994). Global MIRs with strong magnetic fields, preceded by a strong shock, were observed near solar maximum. During the declining phase of the solar cycle, the MIRs had significantly weaker magnetic fields. In both cases the pickup protons, identified by an analysis of pressure balanced structures, play a major role in the dynamical evolution of the MIRs beyond 30 AU. The large-scale magnetic field fluctuations have significantly greater amplitudes near solar maximum than during the declining phase of the solar cycle.

Burlaga, L. F.↗
DOE OSTI2025

An Internal Digital Image Correlation Technique for High-Strain Rate Dynamic Experiments

Full-field, quantitative visualization techniques, such as digital image correlation (DIC), have unlocked vast opportunities for experimental mechanics. However, DIC has traditionally been a surface measurement technique, and has not been extended to perform measurements on the interior of specimens for dynamic, full-scale laboratory experiments. This limitation restricts the scope of physics which can be investigated through DIC measurements, especially in the context of heterogeneous materials. The focus of this study is to develop a method for performing internal DIC measurements in dynamic experiments. The aim is to demonstrate its feasibility and accuracy across a range of stresses (up to 650 MPa), strain rates (10 3 - 10 6 s -1 ), and high-strain rate loading conditions (e.g., ramped and shock wave loading). Internal DIC is developed based on the concept of applying a speckle pattern at an inner-plane of a transparent specimen. The high-speed imaging configuration is coupled to the traditional dynamic experimental setups, and is focused on the internal speckle pattern. During the experiment, while the sample deforms dynamically, in-plane, two-dimensional deformations are measured via correlation of the internal speckle pattern. In this study, the viability and accuracy of the internal DIC technique is demonstrated for split-Hopkinson (Kolsky) pressure bar (SHPB) and plate impact experiments. The internal DIC experimental technique is successfully demonstrated in both the SHPB and plate impact experiments. In the SHPB setting, the accuracy of the technique is excellent throughout the deformation regime, with measurement noise of approximately 0.2% strain. In the case of plate impact experiments, the technique performs well, with error and measurement noise of 1% strain. The internal DIC technique has been developed and demonstrated to work well for full-scale dynamic high-strain rate and shock laboratory experiments, and the accuracy is quantified. Here, the technique can aid in investigating the physics and mechanics of the dynamic behavior of materials, including local deformation fields around dynamically loaded material heterogeneities.

36 MATERIALS SCIENCE↗
NASA NTRS2020

X-ray Evolution of the Nova V959 Mon Suggests a Delayed Ejection and a Non-Radiative Shock

X-ray observations of shocked gas in novae can provide a useful probe of the dynamics of the ejecta. Here we report on X-ray observations of the nova V959 Mon, which was also detected in GeV gamma-rays with the Fermi satellite. We find that the X-ray spectra are consistent with a two-temperature plasma model with non-solar abundances. We interpret the X-rays as due to shock interaction between the slow equatorial torus and the fast polar outflow that were inferred from radio observations ofV959 Mon. We further propose that the hotter component, responsible for most of the flux, is from the reverse shock driven into the fast outflow. We find a systematic drop in the column density of the absorber between days 60 and 140, consistent with the expectations for such a picture. We present intriguing evidence for a delay of around 40 d in the expulsion of the ejecta from the central binary. Moreover, we infer a relatively small (a few times 10−6Msun) ejecta mass ahead of the shock, considerably lower than the mass of 104K gas inferred from radio observations. Finally, we infer that the dominant X-ray shock was likely notradiative at the time of our observations, and that the shock power was considerably higher than the observed X-ray luminosity. It is unclear why high X-ray luminosity, closer to the inferred shock power, is never seen in novae at early times, when the shock is expected to have high enough density to be radiative.

Thomas Nelson↗
NASA NTRS1995

Collapse and fragmentation of molecular cloud cores. 2: Collapse induced by stellar shock waves

The standard scenario for low-mass star formation involves 'inside-out' collapse of a dense molecular cloud core following loss of magnetic field support through ambipolar diffusion. However, isotopic anomalies in presolar grains and meteoritical inclusions imply that the collapse of the presolar cloud may have been triggered by a stellar shock wave. This paper explores 'outside-in' collapse, that is, protostellar collapse initiated directly by the compression of quiescent dense cloud cores impacted by relatively slow stellar shock waves. A second-order accurate, gravitational hydrodynamics code has been used to study both the spherically symmetrical and three-dimensional evolution of initially centrally condensed, isothermal, self-gravitating, solar-mass cloud cores that are struck by stellar shock waves with velocities up to 25 km/s and postshock temperatures of 10 to 10,000 K. The models show that such mild shock waves do not completely shred and destroy the cloud, and that the dynamical ram pressure can compress the cloud to the verge of self-gravitational collapse. However, compression caused by a high postshock temperature is a considerably more effective means of inducing collapse. Shock-induced collapse produces high initial mass accretion rates (greater than 10(exp -4) solar mass/yr in a solar-mass cloud) that decline rapidly to much lower values, depending on the presence (approximately 10(exp -6) solar mass/yr) or absence (approximately 10(exp -8) to 10(exp -7) solar mass/yr) of an infinite reservoir of mass. Stellar mass accretion rates approximately 10(exp -7) solar mass/yr have been previously inferred from the luminosities of T Tauri stars; balanced mass accretion (stellar rate = envelope rate) at approximately 10(exp -7) solar mass/yr could then be possible if accretion occurs from a finite mass reservoir. Fluid tracers are used to determine what fraction of the stellar shock material is incorporated into the resulting protostellar object and disk; roughly half the impinging material is injected into the collapsing cloud core when there is a high postshock temperature. The models are consistent with a scenario where an AGB star wind triggered the collapse of the presolar cloud while injecting about 0.01 solar mass of matter derived from the AGB star envelope, as has been separately inferred on the basis of nucleosynthesis calculations.

Boss, Alan P.↗
NASA NTRS1994

Interferometric investigations of compressible dynamic stall over a transiently pitching airfoil

The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Chandrasekhara, M. S.↗
NASA NTRS1985

Coupling of newborn ions to the solar wind by electromagnetic instabilities and their interaction with the bow shock

The process by which the solar wind assimilates newly ionized atoms is important for understanding the presence of planetary or interstellar helium in the solar wind, the dynamics of the Active Magnetospheric Particle Tracer Explorers (AMPTE) lithium releases in front of the earth's bow shock, and the formation of cometary tails. In this paper is examined how newborn ions can be coupled to the solar wind in the direction parallel to the magnetic field by means of electromagnetic instabilities driven by the distribution of newborn ions. The linear properties of three instabilities are analyzed and compared with numerical solutions of the linear dispersion equation, while their nonlinear behavior is followed by means of computer simulation to obtain the characteristic time for the pickup process. With a primary emphasis on the AMPTE lithiuim releases, various degrees of realism are introduced into the calculations to model the upstream conditions and the intersection of the lithium with the bow shock. It is shown that a time-dependent shock model is needed to correctly reproduce the amount of lithium which is transmitted through the shock and that the resulting lithium ion distribution is still likely to be subject to the same type of instabilities in the magnetosheath. Applications of these results to comets, in particular the artificial comet expected to be generated by the AMPTE barium release in the magnetosheath, is also briefly discussed.

Winske, D.↗
NASA NTRS1973

Design and performance of the Ames electric-arc shock tunnel

A high enthalpy shock tunnel using arc-heated helium as the driver gas was designed for gas dynamic research at total stream energies from 7,000 to 35,000 j/g. The arc driver was found to be a relatively efficient energy converter. Tailored shock Mach numbers from 11.5 to 14.6 in air were demonstrated. A nozzle calibration with a total stream enthalpy of 18,600 j/g of air gave test times of 1.5 to 2.0 m sec at flow Mach numbers from 16 to 24.

Reller, J. O., Jr.↗
NASA NTRS1992

Development of a 3-D upwind PNS code for chemically reacting hypersonic flowfields

Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.

Tannehill, J. C.↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evaluation of a novel mid-infrared pyrometer for dynamic material experiments

Temperature is a complicated thermodynamic parameter to measure in dynamic compression experiments. Optical pyrometry is a general-purpose “work-horse” technique for measuring temperature from a radiant surface on these experimental platforms. The optical pyrometry channels are commonly held to the visible or Near-Infrared spectrum, which provides high fidelity temperature measurement for shock temperature above ∼1200–1500 K. However, low temperature (T < 1200 K) dynamic material experiments, including low pressure or quasi-isentropic studies, as well as experiments with complex thermodynamic paths, require Mid-Infrared (Mid-IR) for high fidelity measurements. This article outlines the design, testing, and characterization of a novel Mid-IR pyrometer system that can be configured between 2.5 and 5.0 µm, suitable for lower temperature measurements and for increasing the fidelity and precision of higher temperature measurements. Experimental validation was done on two separate gas gun platforms, with two separate impact velocities, achieving temperatures between 450 and 1100 K.

Dutra, E. C. (ORCID:0000000167649195)↗
NASA NTRS1993

Some aspects of the aeroacoustics of high-speed jets

Some of the background to contemporary jet aeroacoustics is addressed. Then scaling laws for noise generation by low-Mach-number airflows and by turbulence convected at 'not so low' Mach number is reviewed. These laws take into account the influence of Doppler effects associated with the convection of aeroacoustic sources. Next, a uniformly valid Doppler-effect approximation exhibits the transition, with increasing Mach number of convection, from compact-source radiation at low Mach numbers to a statistical assemblage of conical shock waves radiated by eddies convected at supersonic speed. In jets, for example, supersonic eddy convection is typically found for jet exit speeds exceeding twice the atmospheric speed of sound. The Lecture continues by describing a new dynamical theory of the nonlinear propagation of such statistically random assemblages of conical shock waves. It is shown, both by a general theoretical analysis and by an illustrative computational study, how their propagation is dominated by a characteristic 'bunching' process. That process associated with a tendency for shock waves that have already formed unions with other shock waves to acquire an increased proneness to form further unions - acts so as to enhance the high-frequency part of the spectrum of noise emission from jets at these high exit speeds.

Lighthill, James↗
DOE OSTI2026

ECLEIRS: Exact conservation law embedded identification of reduced states for parameterized nonlinear conservation laws from sparse and noisy data

Multi-query applications such as parameter estimation, uncertainty quantification and design optimization for parameterized partial differential equation (PDE) systems are expensive. While reduced/latent state dynamics approaches for parameterized PDEs offer a viable alternative, these approaches rely on high-quality data and struggle with highly sparse spatiotemporal noisy measurements typically obtained from experiments. Furthermore, there is no guarantee that these models satisfy governing physical conservation laws. In this article, we propose a reduced state dynamics approach, referred to as ECLEIRS, that embeds exact conservation in the solution and flux representation by utilizing a space-time divergence-free neural network formulation. We compare ECLEIRS with other reduced state dynamics approaches, those that do not enforce any physical constraints and those with physics-informed loss functions, for three shock-propagation problems: 1-D advection, 1-D Burgers and 2-D Euler equations. In conclusion, the numerical experiments conducted in this study demonstrate that ECLEIRS provides the most accurate prediction of dynamics for unseen parameters even in the presence of highly sparse and noisy data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Compensating slice emittance growth in high brightness photoinjectors using sacrificial charge

Achieving maximum electron beam brightness in photoinjectors requires detailed control of the 3D bunch shape and precise tuning of the beam focusing. Even in state-of-the-art designs, slice emittance growth due to nonlinear space charge forces and partial nonlaminarity often remains non-negligible. In this work, we introduce a new means to linearize the transverse slice phase space: a sacrificial portion of the bunch’s own charge distribution, formed into a wavebroken shock front by highly nonlinear space charge forces within the gun, whose downstream purpose is to dynamically linearize the desired bunch core. We show that linearization of an appropriately prepared bunch can be achieved via strongly nonlaminar focusing of the sacrificial shock front, while the inner core focuses laminarly. This leads to a natural spatial separation of the two distributions: a dense core surrounded by a diffuse halo of sacrificial charge that can be collimated. Multiobjective genetic algorithm optimizations of the ultracompact x-ray free electron laser injector employ this concept, and we interpret it with an analytic model that agrees well with the simulations. In simulation, we demonstrate a final bunch charge of 100 pC, peak current ∼ 30 A, and a sacrificial charge of 150 pC (250 pC total emitted from cathode) with normalized emittance growth of < 20 nm rad due to space charge. This implies a maximum achievable brightness approximately an order of magnitude greater than existing free electron laser injector designs. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗