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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Solar cycle variation of interplanetary disturbances observed as Doppler scintillation transients

Interplanetary disturbances characterized by plasma that is more turbulence and/or moves faster than the background solar wind are readily defected as transients in Doppler scintillation measurements of the near-Sun solar wind. Systematic analysis of over 23,000 hours of Pioneer Venus Orbiter Doppler measurements obtained inside 0.5 AU during 1979-1987 have made it possible for the first time to investigate the frequency of occurrence of Doppler scintillation transients under solar minimum conditions and to determine its dependence on solar cycle. On the basis of a total of 142 transients, Doppler scintillation transient rates vary from a high of 0.22 in 1979 (one every 4.6 days) to a low of 0.077 transients/d in 1986 (one every 13 days), a decrease by almost a factor of 3 from solar maximum to solar minimum. This solar cycle variation, the strongest yet of any solar wind Doppler scintillation property, is highly correlated with both solar activity characterized by sunspot number and the coronal mass ejection rates deduced from Solswind and Solar Maximum Mission (SMM) coronagraph observations. These results indicate that coronal mass ejections and Doppler scintillation transients are closely related not just during solar maximum, as occasional individual comparisons have shown in the past, but throughout the entire solar cycle, and strengthen the notation that the Doppler scintillation and optical transients are different manifestations of the same physical phenomenon. The magnitudes of the transients, as described by the ratio of peak to pretransient scintillation levels (EF for enhancement factor), and their distribution iwth heliocentric distance also vary with solar cycle. While EF tends to diminish with increasing heliocentric distance during high solar activity, it is more evenly distributed during low solar activity. EF is also lower during solar minimum, as 13% of the transients during solar maximum have values exceeding 23, the highest EF observed during solar minimum. These results are consistent with the fact that occasional major fast-moving interplanetary shocks that are observed during solar maximum are very rate during solar minimum.

Woo, Richard↗
NASA NTRS2009

Implementation of Radiation, Ablation, and Free Energy Minimization Modules for Coupled Simulations of Hypersonic Flow

A description of models and boundary conditions required for coupling radiation and ablation physics to a hypersonic flow simulation is provided. Chemical equilibrium routines for varying elemental mass fraction are required in the flow solver to integrate with the equilibrium chemistry assumption employed in the ablation models. The capability also enables an equilibrium catalytic wall boundary condition in the non-ablating case. The paper focuses on numerical implementation issues using FIRE II, Mars return, and Apollo 4 applications to provide context for discussion. Variable relaxation factors applied to the Jacobian elements of partial equilibrium relations required for convergence are defined. Challenges of strong radiation coupling in a shock capturing algorithm are addressed. Results are presented to show how the current suite of models responds to a wide variety of conditions involving coupled radiation and ablation.

Gnoffo, Peter A.↗
NASA NTRS2011

Flow Simulation of Supersonic Inlet with Bypass Annular Duct

A relaxed isentropic compression supersonic inlet is a new concept that produces smaller cowl drag than a conventional inlet, but incurs lower total pressure recovery and increased flow distortion in the (radially) outer flowpath. A supersonic inlet comprising a bypass annulus to the relaxed isentropic compression inlet dumps out airflow of low quality through the bypass duct. A reliable computational fluid dynamics solution can provide considerable useful information to ascertain quantitatively relative merits of the concept, and further provide a basis for optimizing the design. For a fast and reliable performance evaluation of the inlet performance, an equivalent axisymmetric model whose area changes accounts for geometric and physical (blockage) effects resulting from the original complex three-dimensional configuration is proposed. In addition, full three-dimensional calculations are conducted for studying flow phenomena and verifying the validity of the equivalent model. The inlet-engine coupling is carried out by embedding numerical propulsion system simulation engine data into the flow solver for interactive boundary conditions at the engine fan face and exhaust plane. It was found that the blockage resulting from complex three-dimensional geometries in the bypass duct causes significant degradation of inlet performance by pushing the terminal normal shock upstream.

Kim, HyoungJin↗
NASA NTRS1986

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.↗
DOE OSTI2024

Double cylinder implosion experiments at the National Ignition Facility

Here, cylindrical implosion experiments are used to directly measure instability growth in a convergent geometry, providing a wealth of data for model validation. Double cylinders are a natural extension of the platform and enable measurements at a classically unstable interface, the outer surface of the inner cylinder, which experiences no ablative stabilization from the laser drive. However, the utility of this platform relies upon maintaining adequate axial uniformity of the inner cylinder during the implosion. Although previous smaller-scale double cylinder experiments exhibited acceptable levels of axial uniformity, radiation-hydrodynamics simulations of larger-scale double cylinders predict more axial non-uniformity induced by the impedance mismatch as the shock wraps around the axial ends of the inner cylinder. A mechanism to reduce axial non-uniformity in these larger double cylinder implosions is presented, and preliminary experimental data confirms the efficacy of the selected mitigation approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Re-Examination of the N 2 O + O Reaction

The reaction of N 2 O with O is a key step in consumption of nitrous oxide in thermal processes. It has two product channels, NO + NO (R2) and N 2 + O 2 (R3). The rate constant for R2 has been measured both in the forward and the reverse direction at elevated temperature and is well established. However, the rate constant for the N 2 + O 2 channel (R3) has been difficult to quantify and has significant error limits. Here, the direct reaction on the triplet surface has a barrier of around 40 kcal mol –1 , and it is too slow for the N 2 + O 2 channel to have any practical significance. Recently, Pham and Lin (2022) suggested an alternative low activation energy reaction path that involves intersystem crossing and reaction on the singlet surface. In the present work, we re-examined a wide range of experiments relevant for the N 2 O + O reaction through kinetic modeling, paying attention to the impact of artifacts such as impurities and surface reactions. Experimental results from shock tubes and batch reactors on the final NO yield in N 2 O decomposition, covering temperatures of 973–2200 K and pressures of 0.013–11.5 atm, support k 3 ~ 0, consistent with the high activation energy for reaction on the triplet surface and a low probability of ISC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

X-ray imaging of liquid-liquid Mg cavitation

The coupling of dynamic loading platforms with synchrotron light sources have enabled experiments examining void formation in optically opaque materials. Much like spallation in a solid, cavitation occurs when rarefaction waves collide within a liquid. The resulting deformation leads to a complex flow field with numerous voids. It has been suggested that compaction of such a structure may be a significant source of ejecta (i.e. shallow bubble collapse). Experiments imaging supported and unsupported shocks interacting with a free surface in shock melted Mg with synchrotron radiation are presented. Here, the results are used to examine the structure of a cavitation plane that forms in the unsupported case and provide an estimate of the porosity within this region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ultrafast x-ray detection of low-spin iron in molten silicate under deep planetary interior conditions

The spin state of Fe can alter the key physical properties of silicate melts, affecting the early differentiation and the dynamic stability of the melts in the deep rocky planets. The low-spin state of Fe can increase the affinity of Fe for the melt over the solid phases and the electrical conductivity of melt at high pressures. However, the spin state of Fe has never been measured in dense silicate melts due to experimental challenges. We report detection of dominantly low-spin Fe in dynamically compressed olivine melt at 150 to 256 gigapascals and 3000 to 6000 kelvin using laser-driven shock wave compression combined with femtosecond x-ray diffraction and x-ray emission spectroscopy using an x-ray free electron laser. The observation of dominantly low-spin Fe supports gravitationally stable melt in the deep mantle and generation of a dynamo from the silicate melt portion of rocky planets.

58 GEOSCIENCES↗
DOE OSTI2020

FICH: Feature instability cavitation history [Slides]

We seek to improve our ability to predict and understand the Richtmyer-Meshkov (RM) unstable dynamics of explosively shocked isolated features interfaced with vacuum or gases, and an HE products cavity. Explosively shocked features form as 2D or 3D RM phenomena. Liquid cavitation cause cavitation bubbles when shockwaves collide, and cavitation waves from wave interactions with isolated features form that source mass into the RM unstable sheets or spikes, and the cavitation waves may slice through bulk material into the HE products cavity.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

SZ–X-Ray Surface Brightness Fluctuations in the SPT-XMM Clusters

The hot plasma in galaxy clusters, the intracluster medium, is expected to be shaped by subsonic turbulent motions, which are key for heating, cooling, and transport mechanisms. The turbulent motions contribute to the nonthermal pressure, which, if not accounted for, consequently imparts a hydrostatic mass bias. Accessing information about turbulent motions is thus of major astrophysical and cosmological interest. Characteristics of turbulent motions can be indirectly accessed through surface brightness fluctuations. This study expands on our pilot investigations of surface brightness fluctuations in the Sunyaev–Zel’dovich and in X-ray data by examining, for the first time, a large sample of 60 clusters using both SPT-SZ and XMM-Newton data and spans the redshift range 0.2 < z < 1.5, thus constraining the respective pressure and density fluctuations within 0.6 R $_{500}$. We deem density fluctuations to be of sufficient quality for 32 clusters, finding mild correlations between the peak of the amplitude spectra of density fluctuations and various dynamical parameters. We infer turbulent velocities from density fluctuations with an average Mach number $\mathcal{M}$ 3D = 0.52 ± 0.14, in agreement with numerical simulations. For clusters with inferred turbulent Mach numbers from fluctuations in both pressure, $\mathcal{M}$ P , and density, $\mathcal{M}$ ρ , we find broad agreement between $\mathcal{M}$ P and $\mathcal{M}$ ρ . Our results suggest either a bimodal or a skewed unimodal Mach number distribution, with the majority of clusters being turbulence-dominated (subsonic) while the remainder are shock-dominated (supersonic).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Physical properties of glasses exposed to Earth-facing and trailing-side environments on LDEF

The exposure of 108 glass samples and 12 glass-ceramic samples to Earth-orbit environments permitted measurements which establish the effects of each environment. Examination of five glass types and one glass ceramic located on both the Earth-facing side and the trailing edge revealed no reduction in strength within experimental limits. Strength measurements subjected less than 5 percent of the sample surface area to stresses above 90 percent of the glass's failure strength. Seven micrometeorite or space debris impacts occurred on trailing edge samples. One of those impacts occurred in a location which was subjected to 50 percent of the applied stress at failure. Micrometeorite or space debris impacts were not observed on Earth-facing samples. The physical shape and structure of the impact sites were carefully examined using stereographic scanning electron microscopy. These impacts induce a stress concentration at the damaged region which influences mechanical strength. The flaw size produced by such damage was examined to determine the magnitude of strength degradation in micrometeorite or space-debris impacted glasses. Scanning electron microscopy revealed topographical details of impact sites which included central melt zones and glass fiber production. The overall crater structure is similar to much larger impacts of large meteorite on the Moon in that the melt crater is surrounded by shocked regions of material which fracture zones and spall areas. Residual stresses arising from shock compression and cooling of the fused zone cannot currently be included in fracture mechanics analyses based on simple flaw size examination.

Wiedlocher, David E.↗
NASA NTRS1994

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗
NASA NTRS2005

Warm Molecular Gas Traced with CO J = 7 --> 6 in the Galaxy's Central 2 Parsecs: Dynamical Heating of the Circumnuclear Disks

We present an 11" resolution map of the central 2 pc of the Galaxy in the CO J = 7 --> 6 rotational transition. The CO emission shows rotation about Sgr A* but also evidence for noncircular turbulent motion and a clumpy morphology. We combine our data set with available CO measurements to model the physical conditions in the disk. We find that the molecular gas in the region is both warm and dense, with T approx. 200-300 K and n(sub H2) approx. (5-7) x 10(exp 4) cm(exp -3). The mass of warm molecular gas we measure in the central 2 pc is at least 2000 M(solar), about 20 times the UV-excited atomic gas mass, ruling out a UV heating scenario for the molecular material. We compare the available spectral tracers with theoretical models and conclude that molecular gas is heated with magnetohydrodynamic shocks with v approx. 10-20 km s(exp -1) and B approx. 0.3- 0.5 mG. Using the conditions derived with the CO analysis, we include the other important coolants, neutral oxygen and molecular hydrogen, to estimate the total cooling budget of the molecular material. We derive a mass-to-luminosity ratio of approx. 2-3 M(solar)(L(solar)exp -1), which is consistent with the total power dissipated via turbulent decay in 0.1 pc cells with v(sub rms) approx. 15 kilometers per second. These size and velocity scales are comparable to the observed clumping scale and the velocity dispersion. At this rate, the material near Sgr A* is dissipating its orbital energy on an orbital timescale and cannot last for more than a few orbits. Our conclusions support a scenario in which the features near Sgr A* such as the circumnuclear disk and northern arm are generated by infalling clouds with low specific angular momentum.

galatic nuclei↗
DOE OSTI2025

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Self-Consistent Relativistic Electron Scattering using the Sherlock Scattering Model for X-ray Diagnostics

We present on a new, self-consistent, arbitrary-temperature Romberg integration scheme for modeling electron scattering in materials in a LANL Lagrangian Shock Hydro (LSH) code. Electron beam-target interactions are fundamental to a wide range of scientific and technological applications. When high-energy electron beams hit their target, they may scatter, deposit energy, or ionize the source. These processes govern the behavior and outcomes in nanotechnology manufacturing, electron microscopy, and modern X-ray diagnostics. Simulating these interactions is essential for interpreting experimental results, predicting material responses, and designing efficient tools and experiments. At Los Alamos, this is done using a LSH code, which is a multi-dimension, multi-material, massively parallel, multi-physics code used to simulate applications from asteroid impacts to electron beam interactions. By effectively and efficiently modeling the way that electrons scatter from the beam we can bolster these simulations and more accurately predict experimental outcomes. The model currently implemented in the LSH of interest is based on work by Papp and does not self-consistently preserve momentum in the slightly relativistic regime; here we adopt a model proposed by Braams and Karney and implement a Romberg integration scheme to compute the diffusion tensor. In this paper we will provide background on the Braams-Karney diffusion tensor as well as the Romberg integration scheme we employed to numerically solve for it. We will show that our integration scheme is accurate in solving for the set of scalar potentials used to re-express the diffusion tensor in differential form, and in solving for the diffusion coefficients in the larger LSH code. By using this diffusion tensor rather than the existing Papp one, and numerically integrating it with a Romberg method, we produce much more accurate, self-consistent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

The interstellar medium on the Gamma Cassiopeiae line of sight

The interstellar medium on the Gamma Cas line of sight is studied through the observations of O I, H2, H I, D I, and Ar I absorption features with the Copernicus satellite. By using the velocity structure of the line of sight previously determined through atomic nitrogen, it is demonstrated that these neutrals are located in the same physical regions and that the O I and Ar I abundances could be solar, on the average, relative to N I in all components. The 1066.660 A Ar I line profile is contaminated by an unidentified line, already detected in the spectrum of Zeta Pup. Furthermore, a peculiar behavior of the Ar I abundance in the different components is reported, possibly associated with high-velocity interstellar gas present on the line of sight. It is concluded that the Gamma Cas line of sight could intercept a dense cloud which has been engulfed and disrupted by a shock.

Ferlet, R.↗