Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK WAVE PROFILE”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for a range of high speed Earth entry experiments conducted in the Electric Arc Shock Tube at NASA Ames. The nonequilibrium radiative transport equations are solved in NEQAIR using flow field variables from the full facility CFD simulations of the EAST shock tube performed by US3D ow solver. These physics-based flow calculations lead to a significantly different post-shock gas state and associated radiation field as compared to that based on a simplified but computationally inexpensive calculation for flow over a blunt-body with appropriate initial conditions. The radiation spectra and radiance profiles are computed for an extensive range of wavelengths, from deep VUV to IR, which are pertinent to the emission characteristics of high enthalpy shock waves in air. The radiation properties of the shocked gas are calculated both in the nonequilibrium region at the shock, and in the equilibrium region behind the shock. Numerical predictions are found to be consistent with the experimental observations.

Simulations↗
NASA NTRS2013

An Early and Comprehensive Millimetre and Centimetre Wave and X-Ray Study of SN 2011dh: a Non-Equipartition Blast Wave Expanding into a Massive Stellar Wind

Only a handful of supernovae (SNe) have been studied in multiwavelengths from the radio to X-rays, starting a few days after the explosion. The early detection and classification of the nearby Type IIb SN 2011dh/PTF 11eon in M51 provides a unique opportunity to conduct such observations. We present detailed data obtained at one of the youngest phase ever of a core-collapse SN (days 3-12 after the explosion) in the radio, millimetre and X-rays; when combined with optical data, this allows us to explore the early evolution of the SN blast wave and its surroundings. Our analysis shows that the expanding SN shock wave does not exhibit equipartition (epsilon(sub e)/epsilon(sub B) approx. 1000), and is expanding into circumstellar material that is consistent with a density profile falling like R(exp −2). Within modelling uncertainties we find an average velocity of the fast parts of the ejecta of 15 000 +/- 1800 km/s, contrary to previous analysis. This velocity places SN 2011dh in an intermediate blast wave regime between the previously defined compact and extended SN Type IIb subtypes. Our results highlight the importance of early (approx.1 d) high-frequency observations of future events. Moreover, we show the importance of combined radio/X-ray observations for determining the microphysics ratio epsilon(sub e)/epsilon(sub B).

gamma rays↗
NASA NTRS2018

Nozzle Plume/Shock Interaction Experimental and Computational Sonic Boom Analyses from the NASA Ames 9- by 7-Foot Supersonic Wind Tunnel

A wind tunnel test and a computational study were conducted to investigate the complex interactions between a supersonic nozzle plume and shock waves of differing strengths generated from various aft surfaces typical of supersonic aircraft. These analytically-defined aft surfaces were representative of horizontal tails of various sizes, and an aft deck. CFD simulations of many proposed model configurations allowed for assessments of the detailed flow interactions of components in close proximity to the nozzle, as well as assessments of the nozzle jet flow itself. The evaluation of the computational results for many candidate configurations guided the design of model components. The interactions of the waveforms from these surfaces with the jet exhaust plume can have significant adverse effects on the loudness of the sonic boom if the surfaces are not carefully integrated into an aircraft design. The greatest discrepancy in estimating sonic boom loudness for low-boom flight vehicles is currently in predicting the signatures from the aft part of an aircraft, including the interactions with the plume flow. The objectives of this test were to gain a better understanding of these interactions, and to provide a detailed experimental database from multiple sources for use as validation cases for CFD tool development. The subject test was run in the NASA Ames 9- by 7-Ft Supersonic Wind Tunnel in February 2016 at Mach numbers of 1.6 and 2.0, and was funded by the NASA Commercial Supersonics Technology (CST) Project. The nozzle flow was provided by high-pressure air (HPA) pumped through the model, and pressure signature data were acquired with the NASA 14-inch sonic boom pressure rail. The rail measured the locations of the shocks and expansions at various distances and off-track angles from the model. This enabled the impact of the nozzle plume/shock interactions on the near- and mid-field sonic boom pressure waveforms to be quantified. Schlieren images of the flow field around and behind the model were obtained with an RBOS (Retroreflective Background-Oriented Schlieren) technique to determine the origins of the shock and expansion waves, to identify the shape and boundaries of the plume, and to determine the changes in incoming and exiting waveforms within the plume. A total pressure rake was positioned closely behind the model nozzle in order to measure the total pressure profiles of the flow above, within, and below the nozzle exhaust. Model angles and positions in the tunnel were measured by photogrammetry using two cameras since the lack of a model force balance prevented the measurement of model deflections under load.Navier-Stokes computations using two different CFD codes were compared to the experimental sonic boom pressure signature data, and the rake total pressure data in the plume. A computational schlieren technique was used to compare the computed flow field with the RBOS images. The computational results were also used to complement the test data with flow field quantities that could not be measured, such as Mach number and pressure distributions to distinguish shock waves and expansion waves.

sonic boom↗
DOE OSTI2024

Multiwavelength Afterglow Analysis of GRB 221009A: Unveiling the Evolution of a Cooling Break in a Wind-like Medium

Gamma-ray bursts (GRBs) are the most energetic explosions in the Universe, and their afterglow emission provides an opportunity to probe the physics of relativistic shock waves in an extreme environment. Several key pieces for completing the picture of GRB afterglow physics are still missing, including jet properties, the emission mechanism, and particle acceleration. Here, we present a study of the afterglow emission of GRB 221009A, the most energetic GRB ever observed. Using optical, X-ray, and gamma-ray data up to approximately 2 days after the trigger, we trace the evolution of the multiwavelength spectrum and the physical parameters behind the emission process. The broadband spectrum is consistent with the synchrotron emission emitted by relativistic electrons with its index of p = 2.29 ± 0.02. We identify a break energy at keV and an exponential cutoff at GeV in the observed multiwavelength spectrum. The break energy increases in time from $16.0^{+7.1}_{-4.9}$ keV at 0.65 days to $46.8^{+25.0}_{-15.5}$ keV at 1.68 days, favoring a stellar-wind-like profile of the circumburst medium with k = 2.4 ± 0.1 as in ρ(r) ∝ r -k . The high-energy attenuation at around 0.4 to 4 GeV is attributed to the maximum of the particle acceleration in the relativistic shock wave. This study confirms that the synchrotron process can explain the multiwavelength afterglow emission and its evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A Curated Dataset of Regional Meteor Events with Simultaneous Optical and Infrasound Observations (2006–2011)

We present a curated, openly accessible dataset of 71 regional meteor events simultaneously recorded by optical and infrasound instrumentation between 2006 and 2011. These events were captured during an observational campaign using the all-sky cameras of the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. Each entry provides optical trajectory measurements, infrasound waveforms, and atmospheric specification profiles. The integration of optical and acoustic data enables robust linkage between observed acoustic signals and specific points along meteor trajectories, offering new opportunities to examine shock wave generation, propagation, and energy deposition processes. This release fills a critical observational gap by providing the first validated, openly accessible archive of simultaneous optical–infrasound meteor observations that supports trajectory reconstruction, acoustic propagation modeling, and energy deposition analyses. By making these data openly available in a structured format, this work establishes a durable reference resource that advances reproducibility, fosters cross-disciplinary research, and underpins future developments in meteor physics, atmospheric acoustics, and planetary defense.

astrometry↗
NASA NTRS2010

Determination of the Rate Coefficients of the SO2 plus O plus M yields SO3 plus M Reaction

Rate coefficients of the title reaction R(sub 31) (SO2 +O+M yields SO3 +M) and R(sub 56) (SO2 + HO2 yields SO3 +OH), important in the conversion of S(IV) to S(VI),were obtained at T =970-1150 K and rho (sub ave) = 16.2 micro mol/cubic cm behind reflected shock waves by a perturbation method. Shock-heated H2/ O2/Ar mixtures were perturbed by adding small amounts of SO2 (1%, 2%, and 3%) and the OH temporal profiles were then measured using laser absorption spectroscopy. Reaction rate coefficients were elucidated by matching the characteristic reaction times acquired from the individual experimental absorption profiles via simultaneous optimization of k(sub 31) and k(sub 56) values in the reaction modeling (for satisfactory matches to the observed characteristic times, it was necessary to take into account R(sub 56)). In the experimental conditions of this study, R(sub 31) is in the low-pressure limit. The rate coefficient expressions fitted using the combined data of this study and the previous experimental results are k(sub 31,0)/[Ar] = 2.9 10(exp 35) T(exp ?6.0) exp(?4780 K/T ) + 6.1 10(exp 24) T(exp ?3.0) exp(?1980 K/T ) cm(sup 6) mol(exp ?2)/ s at T = 300-2500 K; k(sub 56) = 1.36 10(exp 11) exp(?3420 K/T ) cm(exp 3)/mol/s at T = 970-1150 K. Computer simulations of typical aircraft engine environments, using the reaction mechanism with the above k(sub 31,0) and k(sub 56) expressions, gave the maximum S(IV) to S(VI) conversion yield of ca. 3.5% and 2.5% for the constant density and constant pressure flow condition, respectively. Moreover, maximum conversions occur at rather higher temperatures (?1200 K) than that where the maximum k(sub 31,0) value is located (approximately 800 K). This is because the conversion yield is dependent upon not only the k(sup 31,0) and k(sup 56) values (production flux) but also the availability of H, O, and HO2 in the system (consumption flux).

Hwang, S. M.↗
NASA NTRS2003

Particle Acceleration in SN1006 Shock Waves

This grant is for the analysis of FUSE observations of particle acceleration in supernova remnant SN1006 shock waves. We have performed quick look analysis of the data, but because the source is faint and because the O VI emission lines on SN1006 are extremely broad, extreme care is needed for background subtraction and profile fitting. Moreover, the bulk of the analysis in will consist of model calculations. The Ly beta and O VI lines are clearly detected at the position in the NW filament of SN1006, but not in the NE position where non-thermal X-rays are strong. The lack of O VI emission in the NE places an upper limit on the pre-shock density there.

Raymond, John C.↗
DOE OSTI2020

Why Are There so Few Reports of High-Energy Electron Drift Resonances? Role of Radial Phase Space Density Gradients

Models of monochromatic Pc5 (2–7 mHz) ultralow frequency (ULF) wave interactions with high energy (greater than ~1 MeV) electrons predict drift resonant interactions that can cause rapid radial transport and acceleration. There are few reports of electron drift resonance at energies greater than ~1 MeV, in contrast to lower energies; moreover, all previous reports occur in the aftermath of interplanetary shocks. These two facts are difficult to reconcile with theory and numerical simulations predicting that greater than ~1 MeV drift resonances should occur more often and in a wider variety of driving conditions. In this study, we show that a combination of observational sampling biases and nominal radial phase space density gradients is one explanation for this discrepancy between theory and observations. In particular, we examine electron dynamics in two case studies with very similar satellite coverage, solar wind conditions, and Pc5 wave properties, yet with different radial phase space density profiles. Using global wave and particle observations, we show that the events have vastly different particle responses despite having similar wave properties. Placing these results in context with past studies, we further show here that nominal radial PSD gradients near geostationary orbit can mask the expected drift resonance particle response and explain (1) the small number of past greater than ~1 MeV drift resonance reports and (2) the restriction of these reports to interplanetary shock events. We argue that future observational studies characterizing radial transport via drift resonance should examine global particle dynamics, including observations of the radial phase space density profile.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Experimental and numerical investigation of an oblique shock wave/turbulent boundary layer interaction with continuous suction

An numerical and experimental investigation has been conducted into the interaction of an incident oblique shock wave with a turbulent boundary layer, for the cases of a rough plate and a porous plate with suction, at a nominal Mach number of 2.5 and flow deflection angles of 0, 4, 6, and 8 deg. Attention is given to the pitot pressure profiles, wall static pressures, and porous plate local bleed distributions measured for the two plates. Suction is found to increase the strength of the incident shock required to separate the boundary layer; for all shock strengths tested, separation is completely eliminated.

Benhachmi, Driss↗
NASA NTRS1949

Gas Motion in a Local Supersonic Region and Conditions of Potential-Flow Breakdown

For a certain Mach number of the oncoming flow, the local velocity first reaches the value of the local velocity of sound (M = 1) at some point on the surface of the body located within the flow. This Mach number is designated the critical Mach number M(sub cr). By increasing the flow velocity, a supersonic local region is formed bounded by the body contour and the line of transition from subsonic to supersonic velocity. As is shown by observations with the Toepler apparatus, at a certain flow Mach number M > M(sub cr) a shock wave is formed near the body that closes the local supersonic region from behind. The formation of the shock wave is associated with the appearance of an additional resistance defined as the wave drag. In this paper, certain features are described of the flow in the local supersonic region, which is bounded by the contour of the body and the transition line, and conditions are sought for which the potential flow with the local supersonic region becomes impossible and a shock wave occurs. In the first part of the paper, the general properties of the potential flow in the local supersonic region, bounded by the contour of the profile and the transition line, are established. It is found that at the transition line, if it is not a line of discontinuity, the law of monotonic variation of the angle of inclination of the velocity vector holds (monotonic law). An approximation is given for the change in velocity at the contour of the body. The flow about a contour having a straight part is studied. In the second part of the paper, an approximation is given of the magnitudes of the accelerations at the interior points of the supersonic region. With the aid of these approximations, it is shown that for profiles convex to the flow the breakdown of the potential flow,associated with an increase of the Mach number of the oncoming flow, cannot be due to the formation of an envelope of the characteristics within the supersonic region. On the basis of the monotonic law, the transitional Mach number M is found, beyond which the potential flow with local supersonic region becomes impossible.

Nikolskii, A. A.↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
NASA NTRS1952

Transonic flow past a wedge profile with detached bow wave

A theoretical study has been made of the aerodynamic characteristics at zero angle of attack of a thin, doubly symmetrical double-wedge profile in the range of supersonic flight speed in which the bow wave is detached. The analysis utilizes the equations of the transonic small-disturbance theory and involves no assumptions beyond those implicit in this theory. The mixed flow about the front half of the profile is calculated by relaxation solution of boundary conditions along the shock polar and sonic line. The purely subsonic flow about the rear of the profile is found by means of the method of characteristics specialized to the transonic small-disturbance theory. Complete calculations were made for four values of the transonic similarity parameter. These were found sufficient to bridge the gap between the previous results of Guderley and Yoshihara at a Mach number of 1 and the results which are readily obtained when the bow wave is attached and the flow is completely supersonic.

Vincenti, Walter G↗
DOE OSTI2026

Staged Z-pinch experiment at the double-EAGLE pulsed power facility

The staged Z-pinch is a potential high-energy gain fusion concept where a high atomic number liner implodes on a deuterium target using a pulsed multi-MA current source. Over the past several years this concept has been studied on 0.5–1.0 MA facilities using Ar and Kr gas puffs injected at R = 1.2 cm. Neutron yield up to (2.5 ± 0.34) x 10 10 was measured, Ruskov et al (2023 IEEE Trans. Plasma Sci. 51 3310–6). Here, in this study, we present experimental results from the 4 MA, 110 ns current rise time Double-EAGLE facility at L3Harris (currently, Fisica Inc.) where a larger radius nozzle created gas density profiles peaked at R = 2.5 cm. Modeling with the MACH2 and FLASH codes indicates that the larger radius allows stronger acceleration of the liner plasma and generation of shock waves that preheat a target plasma layer next to the liner to temperatures $T_i >$1 keV. The resulting counter thermal pressure on the liner plasma limits the growth of the Magneto-Rayleigh–Taylor (MRT) instability at the liner-vacuum boundary and the implosion proceeds in a relatively stable manner. Near bang time the enormous target plasma thermal pressure smoothens the MRT perturbations developed during the earlier implosion stages. Time integrated x-ray pinhole images with cutoff energy of 100 eV confirm that a long (∼3 cm), stable and uniform high energy density plasma column is formed in the final implosion stage. Consistent neutron yield in the 10 10 –10 11 range was measured for both Ar and Kr liners imploding on a deuterium target.

Physics - Plasma physics↗
DOE OSTI2023

Multi-shock design tools in double shell models

The essence of the project was to implement a program in python that demonstrated certain models put forth by fusion physicists. Across the time spent at Los Alamos National Laboratory, much was spent programming. Though fundamentally programmatical in nature, the task contained considerable difficulty hid within mathematically deriving and manipulating the models while maintaining data discretization. The result of solving these complications gave invaluable experience in applied mathematics. The final product, a piece of lengthy code, ended up having the full range of desired capabilities, with plans for more. Given the laser power history from a pulse shot on NIF, Omega, or any other Inertial Confinement Fusion facility, the code maps the terrain of radiation temperature inside the hohlraum throughout the period of the shot. Additionally, shock wave convergence prediction was implemented as additional capability. As the name of the project suggests, this tool acts as a compass of direction for researchers designing laser power profiles. Further development of this tool will undoubtedly assist in the research conducted on single and double-shell indirect ICF, consequently expanding the capabilities of operating stewardship over the national nuclear stockpile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Separated trailing-edge flow at a transonic Mach number

Results from an experimental investigation of a small-scale separated flow at a sharp trailing edge are presented. The separated flow was a result of sustained adverse pressure gradients. Experiments were carried out at a transonic Mach number, without any shock wave in the flow, and at a high Reynolds number. Measurements include surface pressure distributions, mean velocity, turbulent shear stress, turbulent kinetic energy, and some streamwise turbulence intensity profiles in the region from just upstream of separation to downstream into the near-wake, following wake-closure. The results show noticeable changes in the mean and turbulent quantities in the separated region although the wall adverse pressure gradients are relatively weak. Turbulent normal stress becomes important in the vicinity of separation and in the separated zone with regard to both mean-flow development and turbulence-energy production. The experimental data are used to assess the performance of two boundary-layer calculation methods.

Viswanath, P. R.↗
NASA NTRS1979

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. Part 1: Pressure distribution

Analytical solutions are derived which incorporate additional physical effects as higher order terms for the case when the sonic line is very close to the wall. The functional form used for the undisturbed velocity profile is described to indicate how various parameters will be calculated for later comparison with experiment. The basic solutions for the pressure distribution are derived. Corrections are added for flow along a wall having longitudinal curvature and for flow in a circular pipe, and comparisons with available experimental data are shown.

Messiter, A. F.↗