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 1,207 records · Page 67

The solar wind - Advances in our knowledge through two solar cycles

As the Pioneer and Voyager spacecraft have moved outward they have gradually unfolded a view of distant regions of the heliosphere. Information on the solar wind velocity, density and temperature as a function of distance out to more than 40 AU has been gathered. Meanwhile the description of the solar wind has evolved. Long-standing questions on the sources of the wind causing geomagnetic activity were clarified by the discovery of coronal holes and coronal mass ejections. The propagation of the resultant solar wind disturbances through the heliosphere has been studied using both observations and models. Plasma physical processes have been studied. This review focuses on the development of the concepts that have been used to describe the solar wind in the three dimensional heliosphere over the last two solar cycles. Collisionless shocks, transient disturbances in space, disturbance propagation and the distant solar wind are discussed.

Feynman, Joan↗

Molecular line emission models of Herbig-Haro objects. II - HCO(+) emission

We present time-dependent models of the chemistry and temperature of interstellar molecular gas clumps that are exposed to the radiation from propagating stellar-jet shocks. The X-ray, EUV, and FUV radiation from the shock initiates ion chemistry and also heats the gas in the clumps. Using representative parameters, we show that, on the shock transit time between the clumps, the abundances of the ionized molecular species that are produced in the clumps can exceed the values determined from steady state models by several orders of magnitude. Collisional excitation by the heated gas can lead to measurable line emission from several ionized species; as in previous investigations of X-ray-irradiated molecular gas, we find that electron impacts contribute significantly to this process. We apply these results to the interpretation of the HCO(+) line emission that has already been detected in several Herbig-Haro objects. We demonstrate that this picture provides a natural explanation of the fact that the line intensity typically peaks ahead of the associated shock, as well as of the reported low line-center velocities and narrow line widths. We tabulate several diagnostic line intensities of HCO(+) and other molecular species that may be used to infer the physical conditions in the emitting gas.

Wolfire, Mark G.↗

The method of space-time conservation element and solution element-applications to one-dimensional and two-dimensional time-marching flow problems

A nontraditional numerical method for solving conservation laws is being developed. The new method is designed from a physicist's perspective, i.e., its development is based more on physics than numerics. Even though it uses only the simplest approximation techniques, a 2D time-marching Euler solver developed recently using the new method is capable of generating nearly perfect solutions for a 2D shock reflection problem used by Helen Yee and others. Moreover, a recent application of this solver to computational aeroacoustics (CAA) problems reveals that: (1) accuracy of its results is comparable to that of a 6th order compact difference scheme even though nominally the current solver is only of 2nd-order accuracy; (2) generally, the non-reflecting boundary condition can be implemented in a simple way without involving characteristic variables; and (3) most importantly, the current solver is capable of handling both continuous and discontinuous flows very well and thus provides a unique numerical tool for solving those flow problems where the interactions between sound waves and shocks are important, such as the noise field around a supersonic over- or under-expansion jet.

Chang, Sin-Chung↗

Miniature laser ignited bellows motor

A miniature optically ignited actuation device has been demonstrated using a laser diode as an ignition source. This pyrotechnic driven motor provides between 4 and 6 lbs of linear force across a 0.090 inch diameter surface. The physical envelope of the device is 1/2 inch long and 1/8 inch diameter. This unique application of optical energy can be used as a mechanical link in optical arming systems or other applications where low shock actuation is desired and space is limited. An analysis was performed to determine pyrotechnic materials suitable to actuate a bellows device constructed of aluminum or stainless steel. The aluminum bellows was chosen for further development and several candidate pyrotechnics were evaluated. The velocity profile and delivered force were quantified using an non-intrusive optical motion sensor.

Renfro, Steven L.↗

Computational Modeling and Validation for Hypersonic Inlets

Hypersonic inlet research activity at NASA is reviewed. The basis for the paper is the experimental tests performed with three inlets: the NASA Lewis Research Center Mach 5, the McDonnell Douglas Mach 12, and the NASA Langley Mach 18. Both three-dimensional PNS and NS codes have been used to compute the flow within the three inlets. Modeling assumptions in the codes involve the turbulence model, the nature of the boundary layer, shock wave-boundary layer interaction, and the flow spilled to the outside of the inlet. Use of the codes and the experimental data are helping to develop a clearer understanding of the inlet flow physics and to focus on the modeling improvements required in order to arrive at validated codes.

Povinelli, Louis A.↗

Observations of chemical processing in the circumstellar environment

High resolution interferometer and single-dish observations of young, deeply embedded stellar systems reveal a complex chemistry in the circumstellar environments of low to intermediate mass stars. Depletions of gas-phase molecules, grain mantle evaporation, and shock interactions actively drive chemical processes in different regions around young stars. We present results for two systems, IRAS 05338-0624 and NCG 1333 IRAS 4, to illustrate the behavior found and to examine the physical processes at work.

NASA Program Exobiology↗

Particle Acceleration by Dense Impulsive Structures Moving in Ambient Magnetospheric Plasma. 3-D Hybrid Kinetic Modeling and MMS Observations

High resolution observations of dense plasma impulsive structures moving through an ambient background magnetospheric flows were captured by the Magnetospheric Multiscale mission. The observations show particle heating and acceleration, shock-like wave formation, and whistler wave excitation inside the interface between the dense impulsive plasma structures and the ambient plasma. A multiscale hybrid kinetic simulation provides an explanation of the observed wave-particle interactions with the assumption that the dense plasma structures may be represented by plasma clouds which are formed at the magnetopause layer due to reconnection processes. Plain Language Summary: Dense, impulsive plasma structures moving through a background plasma were captured by the NASA Magnetospheric Multiscale mission. The observations show that the dense structures can generate strong perturbations in the electromagnetic field and shock-like waves. Interactions between these electromagnetic waves and the particles results in particle acceleration. 3-D hybrid kinetic modeling (particle description for ions and fluid description for electrons) was used to investigate the plasma physics of the observed structures. It was assumed that the plasma clouds were formed by magnetic field reconnection inside the magnetopause, which is the interface between the solar wind particles and the cold low-density magnetospheric plasma. The work helps us understand the plasma environment at the interface between the Earth and solar wind, near planetary moons, within astrophysical explosions, and possibly at the interface between the solar wind and local interstellar medium.

A. S. Lipatov↗

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗

Imaging the dynamics of initial laser-driven shocks and blowoff plasmas in polystyrene under laser-direct-drive fusion conditions

The dynamics of laser-driven shock propagation in solid ablator material, along with the concomitant processes of shinethrough and blowoff plasma plume expansion, at conditions relevant to laser-direct-drive fusion are measured using a tabletop system capable of supporting rapid dataset development with significantly higher resolution compared to existing methods. The experimental results are directly compared with current radiation-hydrodynamic simulations using the two-dimensional code . Discrepancies between experiments and simulations are evidenced and attributed to limitation of simulation to account for important processes involved in the solid-to-plasma transition. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock Tube Radiation Measurements in Nitrogen

Spectrally and spatially resolved radiance has been measured in the Electric Arc Shock Tube (EAST) facility, with the aim of improving fundamental understanding of high enthalpy flows in pure nitrogen. These tests provide data to inform models used for simulations of high speed flight in nitrogen rich atmospheres, such as Earth or Titan. The experiments presented in this paper cover conditions from approximately 6 km/s to 11 km/s at an initial pressure of 0.2 Torr. A wide range of physics, with different degrees of non-equilibrium and nitrogen dissociation, are covered. The EAST data are presented in different formats for analysis and comparisons. These formats include the spectral radiance at equilibrium (where appropriate), the spatial dependence of radiance over defined wavelength ranges and the mean non-equilibrium spectral radiance (the so-called "spectral non-equilibrium metric"). All the information needed to simulate each experimental trace, including free-stream conditions, shock time of arrival (i.e. x-t) relation, and the spectral and spatial resolution functions, are provided. Equilibrium radiation calculations are shown as a reference. It is the intention of this paper to motivate code comparisons benchmarked against this data set.

shock tube↗

Oxidation and pyrolysis of methyl propyl ether

The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical Investigation of Dual-Mode Scramjet Combustor with Large Upstream Interaction

Dual-mode scramjet combustor configuration with significant upstream interaction is investigated numerically, The possibility of scaling the domain to accelerate the convergence and reduce the computational time is explored. The supersonic combustor configuration was selected to provide an understanding of key features of upstream interaction and to identify physical and numerical issues relating to modeling of dual-mode configurations. The numerical analysis was performed with vitiated air at freestream Math number of 2.5 using hydrogen as the sonic injectant. Results are presented for two-dimensional models and a three-dimensional jet-to-jet symmetric geometry. Comparisons are made with experimental results. Two-dimensional and three-dimensional results show substantial oblique shock train reaching upstream of the fuel injectors. Flow characteristics slow numerical convergence, while the upstream interaction slowly increases with further iterations. As the flow field develops, the symmetric assumption breaks down. A large separation zone develops and extends further upstream of the step. This asymmetric flow structure is not seen in the experimental data. Results obtained using a sub-scale domain (both two-dimensional and three-dimensional) qualitatively recover the flow physics obtained from full-scale simulations. All results show that numerical modeling using a scaled geometry provides good agreement with full-scale numerical results and experimental results for this configuration. This study supports the argument that numerical scaling is useful in simulating dual-mode scramjet combustor flowfields and could provide an excellent convergence acceleration technique for dual-mode simulations.

Mohieldin, T. O.↗

Relation Between Type II Bursts and CMEs Inferred from STEREO Observations

The inner coronagraph (COR1) of the Solar Terrestrial Relations Observatory (STEREO) mission has made it possible to observe coronal mass ejections (CMEs) a in the spatial domain overlapping with that of the metric type II radio bursts. The type II bursts were associated with generally weak flares (mostly B and C class soft X-ray flares), but the CMEs were quite energetic. Using CME data for a set of type II bursts during the declining phase of solar cycle 23, we determine the CME height when the type II bursts start, thus giving an estimate of the heliocentric distance at which CME-driven shocks form. This distance has been determined to be approximately 1.5Rs (solar radii), which coincides with the distance at which the Alfv?n speed profile has a minimum value. We also use type II radio observations from STEREO/WAVES and Wind/WAVES observations to show that CMEs with moderate speed drive either weak shocks or no shock at all when they attain a height where the Alfv?n speed peaks (?3Rs ? 4Rs). Thus the shocks seem to be most efficient in accelerating electrons in the heliocentric distance range of 1.5Rs to 4Rs. By combining the radial variation of the CME speed in the inner corona (CME speed increase) and interplanetary medium (speed decrease) we were able to correctly account for the deviations from the universal drift-rate spectrum of type II bursts, thus confirming the close physical connection between type II bursts and CMEs. The average height (approximately 1.5 Rs) of STEREO CMEs at the time of type II bursts is smaller than that (2.2 Rs) obtained for SOHO (Solar and Heliospheric Observatory) CMEs. We suggest that this may indicate, at least partly, the density reduction in the corona between the maximum and declining phases, so a given plasma level occurs closer to the Sun in the latter phase. In two cases, there was a diffuse shock-like feature ahead of the main body of the CME, indicating a standoff distance of 1Rs - 2Rs by the time the CME left the LASCO field of view.

Gopalswamy, N.↗

Relation Between Type II Bursts and CMEs Inferred from STEREO Observations

The inner coronagraph (COR1) of the Solar Terrestrial Relations Observatory (STEREO) mission has made it possible to observe CMEs in the spatial domain overlapping with that of the metric type II radio bursts. The type II bursts were associated with generally weak flares (mostly B and C class soft X-ray flares), but the CMEs were quite energetic. Using CME data for a set of type II bursts during the declining phase of solar cycle 23, we determine the CME height when the type II bursts start, thus giving an estimate of the heliocentric distance at which CME-driven shocks form. This distance has been determined to be approx. 1.5Rs (solar radii), which coincides with the distance at which the Alfven speed profile has a minimum value.We also use type II radio observations from STEREO/WAVES and Wind/WAVES observations to show that CMEs with moderate speed drive either weak shocks or no shock at all when they attain a height where the Alfven speed peaks (approx. 3Rs - 4Rs). Thus the shocks seem to be most efficient in accelerating electrons in the heliocentric distance range of 1.5Rs to 4Rs. By combining the radial variation of the CME speed in the inner corona (CME speed increase) and interplanetary medium (speed decrease) we were able to correctly account for the deviations from the universal drift-rate spectrum of type II bursts, thus confirming the close physical connection between type II bursts and CMEs. The average height (approx 1.5Rs) of STEREO CMEs at the time of type II bursts is smaller than that (2.2Rs) obtained for SOHO (Solar and Heliospheric Observatory) CMEs. We suggest that this may indicate, at least partly, the density reduction in the corona between the maximum and declining phases, so a given plasma level occurs closer to the Sun in the latter phase. In two cases, there was a diffuse shock-like feature ahead of the main body of the CME, indicating a standoff distance of 1Rs - 2Rs by the time the CME left the LASCO field of view.

Solar Energetic Particles↗