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 775 records · Page 43

Evidence for Dissociation and Ionization in Shock Compressed Nitrogen to 800 GPa

Triple bonding in the nitrogen molecule (N 2 ) is among the strongest chemical bonds with a dissociation enthalpy of 9.8 eV/molecule. Nitrogen is therefore an excellent test bed for theoretical and numerical methods aimed at understanding how bonding evolves under the influence of the extreme pressures and temperatures of the warm dense matter regime. Here, we report laser-driven shock experiments on fluid molecular nitrogen up to 800 GPa and 4.0 g/cm 3 . Line-imaging velocimetry measurements and impedance matching method with a quartz reference yield shock equation of state data of initially precompressed nitrogen. Comparison with numerical simulations using path integral Monte Carlo and density functional theory molecular dynamics reveals clear signatures of chemical dissociation and the onset of L-shell ionization. Furthermore, combining data along multiple shock Hugoniot curves starting from densities between 0.76 and 1.29 g/cm 3 , our study documents how pressure and density affect these changes in chemical bonding and provides benchmarks for future theoretical developments in this regime, with applications for planetary interior modeling, high energy density science, and inertial confinement fusion research.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pioneer 10 observations of the solar wind and its interaction with Jupiter - Plasma electron results

Typical plasma electron spectra obtained during transversal of the Jovian magnetosheath are discussed. Data from the outer magnetosphere show an enhanced energy tail between approximately 50 eV and 200 eV, indicating heating of the solar wind electrons across the Jovian bow shock. Transversal of the inner magnetosphere produced data implying even higher electron temperatures. Results are consistent with ion observations of a high Mach number bow shock and magnetosheath flow field. They are consistent with a contraction of the Jovian magnetosphere resulting from an increase in the solar wind dynamic pressure similar to, but on a larger scale than that observed in the case of the earth's bow shock.

Intriligator, D. S.↗

Graphics and flow visualization in computational fluid dynamics

Techniques for displaying two- and three-dimensional flowfield solutions are described. Several methods of illustrating flow structure are addressed including particle tracing, simulated oil flow, and shock finding. These are incorporated into an interactive graphics program for CFD flowfields, called PLOT3D. Emphasis is made on the difficulty in visualizing three-dimensional flow features, and the importance of color, fast 3D image manipulation, and dynamic movie play-back in displaying such flows. The need for advanced algorithms to identify shock waves, vortices, and separation lines is pointed out. It is likely that the supercomputer will be needed for this process because of the size of 3D and/or unsteady CFD databases.

Buning, P. G.↗

On Theoretical Broadband Shock-Associated Noise Near-Field Cross-Spectra

The cross-spectral acoustic analogy is used to predict auto-spectra and cross-spectra of broadband shock-associated noise in the near-field and far-field from a range of heated and unheated supersonic off-design jets. A single equivalent source model is proposed for the near-field, mid-field, and far-field terms, that contains flow-field statistics of the shock wave shear layer interactions. Flow-field statistics are modeled based upon experimental observation and computational fluid dynamics solutions. An axisymmetric assumption is used to reduce the model to a closed-form equation involving a double summation over the equivalent source at each shock wave shear layer interaction. Predictions are compared with a wide variety of measurements at numerous jet Mach numbers and temperature ratios from multiple facilities. Auto-spectral predictions of broadband shock-associated noise in the near-field and far-field capture trends observed in measurement and other prediction theories. Predictions of spatial coherence of broadband shock-associated noise accurately capture the peak coherent intensity, frequency, and spectral width.

Miller, Steven A. E.↗

Slot Nozzle Effects for Reduced Sonic Boom on a Generic Supersonic Wing Section

NASA has conducted research programs to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas. Restrictions are due to the disturbance from the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Results from two-dimensional computational fluid dynamic (CFD) analyses (performed on a baseline Mach 2.0 nozzle in a simulated Mach 2.2 flow) indicate that over-expanded and under-expanded operation of the nozzle has an effect on the N-wave boom signature. Analyses demonstrate the feasibility of reducing the magnitude of the sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock structure. This work was extended to study the impact of integrating a high aspect ratio exhaust nozzle or long slot nozzle on the trailing edge of a supersonic wing. The nozzle is operated in a highly under-expanded condition, creating a large exhaust plume and a shock at the trailing edge of the wing. This shock interacts with and suppresses the expansion wave caused by the wing, a major contributor to the sonic boom signature. The goal was to reduce the near field pressures caused by the expansion using a slot nozzle located at the wing trailing edge. Results from CFD analysis on a simulated wing cross-section and a slot nozzle indicate potential reductions in sonic boom signature compared to a baseline wing with no propulsion or trailing edge exhaust. Future studies could investigate if this effect could be useful on a supersonic aircraft for main propulsion, auxiliary propulsion, or flow control.

Castner, Raymond S.↗

Slot Nozzle Effects for Reduced Sonic Boom on a Generic Supersonic Wing Section

NASA has conducted research programs to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas. Restrictions are due to the disturbance from the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Results from two-dimensional computational fluid dynamic (CFD) analyses (performed on a baseline Mach 2.0 nozzle in a simulated Mach 2.2 flow) indicate that over-expanded and under-expanded operation of the nozzle has an effect on the N-wave boom signature. Analyses demonstrate the feasibility of reducing the magnitude of the sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock structure. This work was extended to study the impact of integrating a high aspect ratio exhaust nozzle or long slot nozzle on the trailing edge of a supersonic wing. The nozzle is operated in a highly under-expanded condition, creating a large exhaust plume and a shock at the trailing edge of the wing. This shock interacts with and suppresses the expansion wave caused by the wing, a major contributor to the sonic boom signature. The goal was to reduce the near field pressures caused by the expansion using a slot nozzle located at the wing trailing edge. Results from CFD analysis on a simulated wing cross-section and a slot nozzle indicate potential reductions in sonic boom signature compared to a baseline wing with no propulsion or trailing edge exhaust. Future studies could investigate if this effect could be useful on a supersonic aircraft for main propulsion, auxiliary propulsion, or flow control.

Caster, Raymond S.↗

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗

Dynamic experiments to study the α - ε phase transition in cerium

The ability to understand and predict the response of matter at extreme conditions requires knowledge of a material’s equation-of-state including the location of phase boundaries, transition kinetics, and the evolution of material strength. Cerium is a material with a complex phase diagram that continues to attract significant scientific interest. Recent dynamic experiments have provided information on the low-pressure γ–α phase transition, sound speed, and Hugoniot data for the higher-pressure α phase, as well as the incipient shock melt transition. Despite these efforts, there are still regions of the phase diagram that are largely unexplored dynamically, including the high-pressure region below the melt boundary. Along a room temperature isotherm, diamond anvil cell data report a transition to the ϵ phase between 13 and 17 GPa. At higher temperatures, similar diamond anvil cell data show significant disagreement regarding the existence, location, and slope of the ϵ-phase boundary. In this work, double-shock loading was used to access the α–ϵ region of the phase diagram to obtain equation-of-state information and to determine the location of the ϵ-phase boundary for shock loading.

74 ATOMIC AND MOLECULAR PHYSICS↗

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Material effects on dynamics in triple-nozzle gas-puff Z pinches

Here, the gas-puff Z pinch has a long history with myriad applications as an efficient neutron or x-ray source. Its simplicity as a load configuration makes it suitable for studying fundamental plasma physics phenomena such as instabilities and energy transport. For example, the implosion of cylindrical shells onto a fusion fuel are inherently susceptible to instability growth on their external surfaces; if such instabilities are unmitigated, then the consequences in terms of degraded performance can be substantial. Similarly, mitigating heat transport from a hot fuel to its colder surrounding container can make fusion conditions more easily achievable. Here we have conducted a systematic study of triple-nozzle (outer liner, inner liner, fuel) gas puffs using two-dimensional (2D) magnetohydrodynamic simulations to investigate the effect of load material on the relevant dynamics. Analogous to past studies on spherical blast waves and converging shock waves, a trend emerges linking increased radiative cooling, lower adiabatic index, and increased magneto-Rayleigh-Taylor instability growth. Notably, our results suggest that, for the present configuration, Ar radiates less than both Ne and Kr during the early stages of the implosion while mass is being swept up and perturbations begin to seed instability growth. Consequently, pinches with Ar on the outer surface exhibit more stable 2D behavior. Here we also present a parameter scan of thermonuclear neutron yield, Y, as a function of peak current, I pk and dopant concentration with Ne or Ar, depending on the inner liner material. Above 6 MA, our results suggest Y ∝ $I$$^{5}_{pk}$ and even substantial mixing (10% by volume) of Ne into the fuel does not drastically reduce yield, suggesting an Ar/Ne/fuel configuration may reliably achieve DD thermonuclear yields of 10 13 – 10 14 /cm in the 10–20 MA range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

2024 Milestone Report: Site 300 mPDV Optical Fiducials

In the early 2000s Photon Doppler Velocimetry (PDV) replaced the Fabry-Perot many beam system1 and has since become ubiquitous across dynamic experimental platforms to measure velocity and event times such as shock breakout, key variables for high pressure physics research. The advent of optically multiplexed oscilloscope channels to create multiplexed PDV (mPDV) increased portability and reduced price per data point. However, these advantages came with the cost of additional fielding complexity and thermally induced timing drift, which directly affects the ability to use PDV for high precision time measurement. Controlling the temperature or incorporating optical fiducials allows researchers to characterize, reduce, and correct this thermal drift in analysis. Our implemented optical fiducial or “timing marker” allows us to reduce thermal drift uncertainty from the 10s of nanoseconds down to the 100s of picoseconds with minimal added complexity to existing systems. The implementation of this timing marker normalizes uncertainty across optical delays, brings our facility cross timing into the sub nanosecond regime and lets us identify/study anomalies in our data. This improvement increases experiment reliability and quality enabling a new class of high precision experiments at S300.

42 ENGINEERING↗

Energy deposition in the earth's atmosphere due to impact of solar activity-generated disturbances

Energy deposition in and dynamic responses of the terrestrial atmosphere to solar flare-generated shocks and other physical processes - such as particle precipitation and local heating - are investigated self-consistently in the context of hydrodynamics, the problem being treated as an initial boundary-value problem. It is extremely difficult to construct a general model for the line solar activity-magnetosphere-atmosphere; however, a limited model for this link is possible. The paper describes such a model, and presents some results on energy deposition into the earth's atmosphere due to solar activity-generated disturbances. Results from the present calculations are presented and discussed.

Wu, S. T.↗

Flow downstream of the heliospheric terminal shock. I - Irrotational flow

Recent reports of remote detection of the heliospheric terminal shock place it near 50 AU. These conflict with standard models which, when combined with current data on the local interstellar medium, place the shock beyod 100 AU. Resolution of this discrepancy has led to hypotheses that invoke cosmic ray pressure, momentum exchange with interstellar neutrals, and magnetic field effects between the shock and the contact discontinuity dividing the solar wind from interstellar plasma. These hypotheses depend not only on properties of the interstellar medium, but also on the downstream three-dimensional flow between the shock and the contact discontinuity, in the region called the 'heliosheath'. The downstream flow field in the absence of magnetic fields is examined here under the assumptions that the flow everywhere outside the shock can be approximated as irrotational and incompressible. It is found, in particular, that the distance between the terminal shock and the contact discontinuity is less than the heliocentric distance to the terminal shock, effectively eliminating magnetic field effects in the heliosheath as being dynamically important.

Suess, Steven T.↗

Vibrational relaxation in hypersonic flow fields

Mathematical formulations of vibrational relaxation are derived from first principles for application to fluid dynamic computations of hypersonic flow fields. Relaxation within and immediately behind shock waves is shown to be substantially faster than that described in current numerical codes. The result should be a significant reduction in nonequilibrium radiation overshoot in shock layers and in radiative heating of hypersonic vehicles; these results are precisely the trends needed to bring theoretical predictions more in line with flight data. Errors in existing formulations are identified and qualitative comparisons are made.

Meador, Willard E.↗

Theory of Type 3 and Type 2 Solar Radio Emissions

The main features of some current theories of type III and type II bursts are outlined. Among the most common solar radio bursts, type III bursts are produced at frequencies of 10 kHz to a few GHz when electron beams are ejected from solar active regions, entering the corona and solar wind at typical speeds of 0.1c. These beams provide energy to generate Langmuir waves via a streaming instability. In the current stochastic-growth theory, Langmuir waves grow in clumps associated with random low-frequency density fluctuations, leading to the observed spiky waves. Nonlinear wave-wave interactions then lead to secondary emission of observable radio waves near the fundamental and harmonic of the plasma frequency. Subsequent scattering processes modify the dynamic radio spectra, while back-reaction of Langmuir waves on the beam causes it to fluctuate about a state of marginal stability. Theories based on these ideas can account for the observed properties of type III bursts, including the in situ waves and the dynamic spectra of the radiation. Type 11 bursts are associated with shock waves propagating through the corona and interplanetary space and radiating from roughly 30 kHz to 1 GHz. Their basic emission mechanisms are believed to be similar to those of type III events and radiation from Earth's foreshock. However, several sub-classes of type II bursts may exist with different source regions and detailed characteristics. Theoretical models for type II bursts are briefly reviewed, focusing on a model with emission from a foreshock region upstream of the shock for which observational evidence has just been reported.

Robinson, P. A.↗

A Parametric Study of the Aerodynamic Characteristics of Nose-Cylinder-Flare Bodies at a Mach Number of 6.0

Force tests were conducted at a Mach number of 6.0 on nose-cylinder-flare bodies to determine the effect of nose shape, cylinder length, flare angle, and flare length on the longitudinal aerodynamic characteristics. A particular investigation was conducted to determine the effect of flare angle for constant flare length, surface area, and diameter. Results indicated that at a Reynolds number of approximately 0.92 x l0 (exp 6) (based on body diameter), the boundary-layer separation effects were significant only with respect to the slope of the normal-force and pitching-moment curve at low angles of attack. The variations of the aerodynamic characteristics with the various parameters were, in general, similar to those predicted by Newtonian theory below a flare angle of 30 degrees and a ratio of flare base diameter to cylinder diameter of less than approximately 2.2. The limiting diameter ratio is consistent with the extent of the low-constant dynamic-pressure region near the body caused by the bow-shock influences as predicted by axisymmetric characteristic theory. The effects of the various parameters for the flares that exceeded the limiting diameter ratio follow the trends predicted by the computed flow-field properties. The axial force for these flare configurations at zero angle of attack was, in general, computed within 10 percent by using these properties. For a constant flare length and surface area the flare effectiveness increased with increasing flare angle; however, for constant flare diameter only the axial-force coefficient was affected by flare angle.

Ashby, George C., Jr.↗

Identifying the Community Structure of the Food-Trade International Multi-Network

Achieving international food security requires improved understanding of how international trade networks connect countries around the world through the import-export flows of food commodities. The properties of international food trade networks are still poorly documented, especially from a multi-network perspective. In particular, nothing is known about the multi-network's community structure. Here we find that the individual crop-specific layers of the multi-network have densely connected trading groups, a consistent characteristic over the period 2001-2011. Further, the multi-network is characterized by low variability over this period but with substantial heterogeneity across layers in each year. In particular, the layers are mostly assortative: more-intensively connected countries tend to import from and export to countries that are themselves more connected. We also fit econometric models to identify social, economic and geographic factors explaining the probability that any two countries are co-present in the same community. Our estimates indicate that the probability of country pairs belonging to the same food trade community depends more on geopolitical and economic factors-such as geographical proximity and trade-agreement co-membership-than on country economic size and/or income. These community-structure findings of the multi-network are especially valuable for efforts to understand past and emerging dynamics in the global food system, especially those that examine potential 'shocks' to global food trade.

complex networks↗