Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple shocks”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Benchmarking impact hydrocodes in the strength regime: Implications for modeling deflection by a kinetic impactor

The Double Asteroid Redirection Test (DART) is a NASA-sponsored mission that will be the first direct test of the kinetic impactor technique for planetary defense. Here, the DART spacecraft will impact into Didymos-B, the moon of the binary system 65803 Didymos and the resulting period change will be measured from Earth. Impact simulations will be used to predict the crater size and momentum enhancement expected from the DART impact. Because the specific material properties (strength, porosity, internal structure) of the Didymos-B target are unknown, a wide variety of numerical simulations must be performed to better understand possible impact outcomes. This simulation campaign will involve a large parameter space being simulated using multiple different shock physics hydrocodes. In order to understand better the behaviors and properties of numerical simulation codes applicable to the DART impact, a benchmarking and validation program using different numerical codes to solve a set of standard problems was designed and implemented.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

A modular, multi-diagnostic, automated shock tube for gas-phase chemistry

A new shock tube has been constructed for investigations of high-temperature chemical kinetics with an emphasis on combustion chemistry. This instrument includes a diaphragmless driver and electrical control of valving. A diaphragmless design significantly improves repeatability of experimental conditions vs the use of diaphragms and leads to an approximate order of magnitude reduction in turnaround time between experiments. Electrical control of valves, combined with diaphragmless operation, also enables remote and automated operation of the shock tube. The design allows for both incident and reflected shock experiments with multiple diagnostics. The performance of the shock tube is demonstrated by reproducing previous literature measurements on the unimolecular decomposition of isobutyl nitrite and cyclohexene.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Empirical estimation of weather-driven yield shocks using biophysical characteristics for U.S. rainfed and irrigated maize, soybeans, and winter wheat

Agricultural yields are highly susceptible to changes in weather system patterns, including both annual and sub-annual changes in temperature and precipitation. Understanding the impacts of future changes in meteorological variables on crop yields have been widely studied in both empirical and process-based models. This work presents a structural econometric approach that combines historical weather data with information about the biophysical growth cycles of maize, winter wheat, and soy to predict the year-to-year yield responses for rainfed crops. Modeled soil moisture and temperature are taken as key predictors, which allows testing the fitted rainfed model’s ability to predict irrigated yields by assuming irrigation produces the numerically optimal level of soil moisture. This approach is grounded in known biophysical processes, which improves confidence in the model predictions. It produces estimates of interannual variability in yields, which enable study of the effects of extreme events on agricultural productivity. Finally, it enables prediction of the potential impacts of changing weather patterns on irrigated crops in areas that are currently primarily rainfed. We present the results of the empirical model, fitted with rainfed data; out-of-sample validation on irrigated crops; and projections of yield shocks under multiple future climate scenarios. Under a bias-corrected GFDL RCP8.5 scenario, this approach predicts average yield change, relative to 2006-2020 average yields, across U.S. counties in the 2040-2060 period of -16.9% and -13.4% for rainfed and irrigated maize, of -19.1% and -18.4% for rainfed and irrigated soy, and of -4.0% and -2.0% for rainfed and irrigated winter wheat. And in the 2070-2090 period of -30.3% and -28.8 % for rainfed and irrigated maize, of -33.1% and -32.8% for rainfed and irrigated soy, and of -6.5% and -5.2%for rainfed and irrigated winter wheat.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Evaluation of earthquake stress parameters and its scaling during the 2016-2017 Amatrice-Norcia-Visso sequence -- Part I

The Amatrice–Norcia–Visso sequence is distinguished by complex behaviour that is somewhat atypical of main-shock–aftershock sequences, as there were multiple large main shocks that continued for months. Here, we focus on the Amatrice sequence (main shock 2016 August 24, M w = 5.97) to evaluate the apparent stress values and magnitude-dependent scaling in order to improve our knowledge of processes that control small and large earthquakes within this active region of Italy. Apparent stress is proportional to the ratio of radiated seismic energy and seismic moment, and as such, these stress parameters play an important role in hazard prediction as they have a strong effect on the observed and predicted ground shaking. We analyse 83 events of the sequence from 2016 August 24 to October 16, within a radius of 20 km from the main shock and with an M w ranging between 5.97 and 2.72. Taking advantage of the averaging nature of coda waves, we analyse coda-envelope-based spectral ratios between neighbouring event pairs. We use equations proposed by Walter et al. to consider stable, low-frequency and high-frequency spectral ratio levels which provide measures of the corner frequency and apparent stress ratios of the events within the sequence. The results demonstrate non-self-similar behaviour within the sequence, indicating a change in dynamics between the largest events and the smaller aftershocks. The apparent stress and corner frequency estimates are compared to those obtained by Malagnini and Munafò who utilized hundreds of direct S-wave spectral ratio measurements to obtain their results. Although our analysis is based only on 83 events, our results are in very good agreement, demonstrating once more that the use of coda waves is very stable and provides lower variance measures than those using direct waves. A comparison with recent Central Apennines source scaling models derived from various seismic sequences (1997–1998 Colfiorito, 2002 San Giuliano di Puglia, 2009 L'Aquila) shows that the Amatrice sequence source scaling in this study is well represented by the models proposed by Pacor et al. and Malagnini and Mayeda.

58 GEOSCIENCES↗
DOE OSTI2015

HTK Lethality Analysis CET

Sandia National Laboratories (Sandia) and the Lockheed Martin Corporation (LMC) collaborated to establish performance of the hit-to-kill lethal mechanism compared to other interceptor kill mechanisms. Multiple eulerian shock physics and couple shock physics/structural simulations were performed to evaluate the performance of two interceptor kill mechanisms – hit to kill and blast fragment – against a conventional target. These simulations showed each kill mechanism’s advantages and disadvantages. Sandia/DOE benefitted by exercising the Sandia-developed computational codes with scenarios that had not been simulated before, providing insight and experience with highly relevant and timely intercept scenarios that will be representative of future work. LMCO benefitted by the experience using the coupled code which they acquired from SNL. They also were provided data that was relevant to interceptor systems being considered and in development. The public benefitted by providing information that will be useful in choosing an interceptor system that will be used for national defense at lower cost than relying solely on an extensive testing program.

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

High-throughput screening of Protein-DNA interactions in Sorghum

Drought is the primary abiotic factors negatively affecting agriculture, causing tremendous reductions in yield world-wide. Drought responsive adaptation strategies require complex molecular responses to including hormone signaling and production of osmo-protectants, as well as developmental responses including stomatal closure, changes to leaf morphology, and altered root system architecture. A common strategy to better understand these adaptation mechanisms is to compare plants with differing drought tolerances and has been successful in multiple crop plants including wheat, rice, and maize. To better understand drought adaptation strategies, we compared the cytosolic- and organelle-enriched protein profiles of leaves from two Sorghum bicolor genotypes with differing pre-flowering drought tolerances after 8-weeks of growth under water limitation. In agreement with previous literature findings, we observed significant drought-induced changes in the abundance of multiple heat shock proteins and dehydrins in both genotypes. We also identified unique proteins that might play a role in genotype dependent drought adaptation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Electron Heating in 2D Particle-in-cell Simulations of Quasi-perpendicular Low-beta Shocks

Abstract We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta ( β p = 0.25) collisionless ion–electron shocks with mass ratio m i / m e = 200, fast Mach number M ms = 1 –4, and upstream magnetic field angle θ Bn = 55°–85° from the shock normal n ˆ . It is known that shock electron heating is described by an ambipolar, B -parallel electric potential jump, Δ ϕ ∥ , that scales roughly linearly with the electron temperature jump. Our simulations have Δ ϕ ∥ / ( 0.5 m i u sh 2 ) ∼ 0.1 –0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measure ϕ ∥ , including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate of ϕ ∥ in our low- β p shocks. We further focus on two θ Bn = 65° shocks: a M s = 4 ( M A = 1.8 ) case with a long, 30 d i precursor of whistler waves along n ˆ , and a M s = 7 ( M A = 3.2 ) case with a shorter, 5 d i precursor of whistlers oblique to both n ˆ and B ; d i is the ion skin depth. Within the precursors, ϕ ∥ has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the M s = 4 , θ Bn = 65° case, ϕ ∥ shows a weak dependence on the electron plasma-to-cyclotron frequency ratio ω pe /Ω ce , and ϕ ∥ decreases by a factor of 2 as m i / m e is raised to the true proton–electron value of 1836.

Tran, Aaron (ORCID:0000000334834890)↗
NASA NTRS1982

Magnetic field studies by Voyager 2 - Preliminary results at Saturn

Results of Voyager 2 studies of the magnetosphere and planetary magnetic field of Saturn are presented. Magnetometer studies have confirmed the results obtained by Voyager 1, indicating the magnetic field to be that of a centered dipole of moment 0.21 gauss Saturn radii-cubed, tilted approximately 0.8 deg from the rotation axis and a maximum measured field intensity of 1187 nT at latitude 17.3 deg N just before periapsis. Voyager 2 observed multiple bow shock and magnetopause crossings during its inbound and outbound trajectories, which were complementary to those of Voyager 1, including magnetopause crossing at 18.5 Saturn radii on the inbound trajectory, and at 48.4-50.9 Saturn radii outbound indicative of magnetospheric expansion due to changing solar wind conditions. Throughout the outbound passage, the magnetospheric field was observed to be relatively steady and smooth, with no evidence for any azimuthal asymmetry or magnetic anomaly. Results thus are incapable of accounting for the observed periodic modulation of the Saturnian kilometric radio emissions.

Ness, N. F.↗
NASA NTRS1983

Multiple spacecraft observations of interplanetary shocks: Characteristics of the upstream ULF turbulence

All interplanetary shocks observed by ISEE-3 and either ISEE-1 or ISEE-2 or both in 1978 and 1979 are examined for evidence of upstream waves. In order to characterize the properties of these shocks it is necessary to determine accurate shock normals. An overdetermined set of equations were inverted to obtain shock normals, velocities and error estimates for all these shocks. Tests of the method indicate it is quite reliable. Using these normals the Mach number and angle were between the interplanetary magnetic field and the shock normal for each shock. The upstream waves were separated into two classes: whistler mode precursors which occur at low Mach numbers and upstream turbulence whose amplitude at Mach numbers greater than 1.5 is controlled by the angle of the field to the shock normal. The former waves are right hand circularly polarized and quite monochromatic. The latter waves are more linearly polarized and have a broadband featureless spectrum.

Russell, C. T.↗
NASA NTRS1983

Multiple spacecraft observations of interplanetary shocks ISEE three-dimensional plasma measurements

ISEE 1 and ISEE 3 three-dimensional solar wind plasma measurements are used together with magnetic field measurements across five previously studied interplanetary shocks to test the accuracy of the mixed-mode shock-normal determination technique and to test whether the shock properties are best approximated with a ratio of specific heats of 5/3 or 2. In the shocks examined, the assumption that the velocity jump was along the normal provided an estimate of the shock normal within 15 deg of the best fit normal 50 percent of the time and within 50 deg, 90 percent of the time. The mixed-mode normals lay within 12 deg of the best fit normal 50 percent of the time and within 36 deg, 90 percent of the time. Part of this deviation may be due to differences in the orientation of the local normal from that of the average normal. Finally, the jump in plasma and field across the shock is better predicted from the Rankine-Hugoniot equations using a ratio of specific heats of 5/3 rather than 2.

Russell, C. T.↗
NASA NTRS1988

An improved version of NCOREL: A computer program for 3-D nonlinear supersonic potential flow computations

A computer code called NCOREL (for Nonconical Relaxation) has been developed to solve for supersonic full potential flows over complex geometries. The method first solves for the conical at the apex and then marches downstream in a spherical coordinate system. Implicit relaxation techniques are used to numerically solve the full potential equation at each subsequent crossflow plane. Many improvements have been made to the original code including more reliable numerics for computing wing-body flows with multiple embedded shocks, inlet flow through simulation, wake model and entropy corrections. Line relaxation or approximate factorization schemes are optionally available. Improved internal grid generation using analytic conformal mappings, supported by a simple geometric Harris wave drag input that was originally developed for panel methods and internal geometry package are some of the new features.

Siclari, Michael J.↗
NASA NTRS1992

Multiple switch-off slow shock solutions

A new electromagnetic quasi-neutral time-dependent implicit hybrid model for simulating the slow-mode switch-off shock is developed, in which the electrons are treated as a finite-mass adiabatic fluid (gamma sub e = 5/3), whose motion along the shock normal provides local charge neutrality. The model is used to examine the effects of various boundary conditions applied during shock formation. It is shown that, for each upstream condition, two distinctly different slow-shock structures can be obtained: a low-shock structure with a trailing magnetic wave (TMW) and a slow-shock structure without a TMW, depending strongly on the boundary conditions for the electric and magnetic fields at the physical piston which generates the shock. The results indicate that the slow-shock structure with a TMW is nonlinearly unstable with respect to changes in the upstream magnetic field, suggesting that the downstream structure of slow shocks need not be associated with a TMW.

Vu, H. X.↗
NASA NTRS1993

Observations of entrainment and time variability in the HH 47 jet

We present new Fabry-Perot images of the HH 47 jet that show the first clear evidence for entrainment in a jet from a young star. The material in the jet moves faster down the axis of the flow and slower at the edges, similar to viscous flow in a pipe. The higher excitation lines occur along the edges of the jet, as expected if entrainment accelerates and heats the ambient material. We confirm previous observations of multiple bow shocks in this system. Together, time variability and entrainment produce much of the observed shock-excited gas in this object. Our data show that the 'wiggles' along the jet are not caused by jet material tied to a spiraling magnetic field, but instead result from time variability, variable ejection angles, or inhomogeneities in the flow. The gas entrained in the HH 47 jet may be atomic; our results do not provide direct evidence that stellar jets drive molecular outflows.

Hartigan, Patrick↗
NASA NTRS1959

The Total-Pressure Recovery and Drag Characteristics of Several Auxiliary Inlets at Transonic Speeds

Several flush and scoop-type auxiliary inlets have been tested for a range of Mach numbers from 0.55 to 1.3 to determine their transonic total-pressure recovery and drag characteristics. The inlet dimensions were comparable with the thickness of the boundary layer in which they were tested. Results indicate that flush inlets should be inclined at very shallow angles with respect to the surface for optimum total-pressure recovery and drag characteristics. Deep, narrow inlets have lower drag than wide shallow ones at Mach numbers greater than 0.9 but at lower Mach numbers the wider inlets proved superior. Inlets with a shallow approach ramp, 7 deg, and diverging ramp walls which incorporated boundary-layer bypass had lower drag than any other inlet tested for Mach numbers up to 1.2 and had the highest pressure recovery of all of the flush inlets. The scoop inlets, which operated in a higher velocity flow than the flush inlets, had higher drag coefficients. Several of these auxiliary inlets projected multiple, periodic shock waves into the stream when they were operated at low mass-flow ratios.

Dennard, John S.↗
NASA NTRS2011

A Chandra X-Ray Survey of Ejecta in the Cassiopeia A Supernova Remnant

We present a survey of the X-ray emitting ejecta in the Cassiopeia A supernova remnant based on an extensive analysis of over 6000 spectral regions extracted on 2.5-10" angular scales using the Chandra 1 Ms observation. We interpret these results in the context of hydrodynamical models for the evolution of the remnant. The distributions of fitted temperature and ionization age are highly peaked and suggest that the ejecta were subjected to multiple secondary shocks. Based on the fitted emission measure and element abundances, and an estimate of the emitting volume, we derive masses for the X-ray emitting ejecta as well as showing the distribution of the mass of various elements over the remnant. The total shocked Fe mass appears to be roughly 0.14 Solar Mass, which accounts for nearly all of the mass expected in Fe ejecta. We find two populations of Fe ejecta, that associated with normal Si-burning and that associated with alpha-rich freeze-out, with a mass ratio of approximately 2:1. Surprisingly, essentially all of this Fe (both components) is well outside the central regions of the SNR, presumably having been ejected by hydrodynamic instabilities during the explosion. We discuss this, and its implications for the neutron star kick.

Hwang, Una↗
NASA NTRSDate not supplied

Implementation of a Compact Self-Aligned Focusing Schlieren System for Measurements with the SBLT-II Flight Test Article

Previous flight testing of a supersonic natural laminar flow airfoil on the NASA F-15Baircraft showed evidence of multiple impinging shocks (shock train) of unknown origin, which generated sufficient crossflow on the airfoil to affect the primary measurement goal of transition location. A compact self-aligned focusing schlieren (SAFS) system is being developed to visualize this shock train on the forward surface of the Supersonic Boundary Layer Transition - PhaseII (SBLT-II) flight test article. Details of the design of the SAFS flight system, the constraints imposed by the imaging pod that will be mounted to the F-15B aircraft, preliminary laboratory results obtained with the system, and challenges associated with integrating the imaging system into the aircraft are presented. The resulting schlieren images obtained from this flight test will allow for a three-dimensional reconstruction of the unknown shock structure, validating the use of SAFS for flight testing.

Joshua M. Weisberger↗
DOE OSTI2024

Laser pulse-length dependent ablation and shock generation in silicon at 5 × 10 14 W/ cm 2 intensities

The effect of laser pulse duration on energy coupling into a planar silicon target is investigated in experiments at the OMEGA-EP facility by varying the laser pulse length τ —spanning 3 orders of magnitude from 100 ps to 10 ns—while maintaining a constant peak laser intensity, I 0 = 5 × 10 14 W/ cm 2 . In theoretical models, the ablation pressure primarily scales for a given material with laser intensity and wavelength, which are all fixed variables here, allowing us to explore the specific role of laser pulse duration. Two-dimensional radiation-hydrodynamics simulations benchmarked with optical probing of the expanding plasma show that the pulse duration is critical for the ablation pressure to reach a steady state. Moreover, the pulse duration impacts shock decay and multiple wave effects, which strongly dictate the evolving shock profile that propagates within the laser-shocked target as ultimately measured by rear-surface diagnostics. The shock velocities inferred from the theoretical model, after considering shock decay, impedance matching, and shock Hugoniot, are found to be in good agreement with velocimetry measurements. However, discrepancies are observed with simulations for the shorter (0.1 ns) and longer (10 ns) pulse durations, which are respectively attributed to unaccounted contributions of kinetic absorption mechanisms and instabilities in simulations. Published by the American Physical Society 2024

Bailly-Grandvaux, M. (ORCID:0000000175294013)↗