Search NASASearch

SEARCH · Search NASA

Results for “shock waves”

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.

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (approx. 9 inches from the source) dominated by direct wave propagation, mid-field environment (approx. 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This report documents the outcome of the assessment.

Steven J Gentz

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (~9 inches from the source) dominated by direct wave propagation, mid-field environment (approximately 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This document contains appendices to the Volume I report.

Steven J Gentz

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading: Appendices

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (approx. 9 inches from the source) dominated by direct wave propagation, mid-field environment (approx. 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This document contains appendices to the Volume I report.

Gentz, Steven J

Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading: Appendices - Volume 2, Part 1

The NASA Engineering and Safety Center (NESC) received a request to develop an analysis model based on both frequency response and wave propagation analyses for predicting shock response spectrum (SRS) on composite materials subjected to pyroshock loading. The model would account for near-field environment (approximately 9 inches from the source) dominated by direct wave propagation, mid-field environment (approximately 2 feet from the source) characterized by wave propagation and structural resonances, and far-field environment dominated by lower frequency bending waves in the structure. This document contains appendices to the Volume I report.

Steven J Gentz

Fluid-Thermal-Structural Interactions Induced by an Asymmetric Shock-Wave/Boundary-Layer Interaction in a Mach-6 Compression Corner

An experimental study is conducted of the fluid-thermal-structural interaction of a clamped compliant panel exposed to a three dimensional shock-wave/boundary-layer interaction (SWBLI) induced by a Mach-6 compression ramp with a spanwise nonuniform incoming boundary layer. The nonuniform boundary layer was produced by placing trips on one side of the upstream flat plate, resulting in largely turbulent flow on the tripped side and transitional flow on the untripped side. Measurements of the flowfield confirmed that the tripped boundary layer contained elevated levels of unsteadiness, and the SWBLI was observed to vary from attached to fully separated as the ramp angle was increased from 10◦ to 38◦; the separation region on the tripped side of the panel was noticeably smaller, showing the elevated turbulence levels of the tripped-side flow to remain relatively localized rather than diffusing across the whole model. Full-field, time-resolved panel deformations were measured using high-speed photogrammetry and the vibrational response at each compression angle was characterized. Although the measured modes conformed largely to those from classical clamped-plate theory, some skewing of the mode shapes was observed. IR thermography highlighted regions of the compliant region where elevated temperatures were likely to promote thermal softening effects to the transient panel response. The quasi-static deformation and stress field was used to characterize the internal stress factor of each mode and showed a meaningful relationship between transient panel response and stress contained within each mode: modes with antinodes lying in high-stress areas of the plate tended to exhibit increases in vibrational frequency and decreases in vibrational power, whereas the opposite was true for modes with antinodes in low-stress areas.

Spectral Proper Orthogonal Decomposition

Boundary Layer Analysis of Shock Tube Flows

Shock tubes offer a controlled environment to reproduce kinetic and radiative phenomena characteristic of atmospheric entry flows under ground-test conditions. The boundary layer developing behind the incident shock wave determines the available test time, influences particle residence times important for similarity scaling, and can affect radiative energy transport. In this study, we couple a quasi-1D space marcher with the compressible boundary-layer equations to numerically compute the post-shock flow in a shock frame of reference for various test gas mixtures representative of different planetary atmospheres. The solvers are individually verified against CFD simulations and analytical correlations available in the literature. Coupled solutions are computed for a finite-rate chemistry in the boundary layer and a non-catalytic isothermal wall. Results yield refined estimates of the maximum separation distance, as well as insights into concentration profiles of relevant species within the boundary layer.

Andrea Fagnani

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

Time-resolved, high-speed self-aligned focusing schlieren images were acquired in the NASA Glenn Research Center Transonic Turbine Blade Cascade facility to help understand the aerodynamic behavior of high-pressure, thick trailing edge turbine blades. The trailing edge thickness of 9% of axial chord tested represents simulated ceramic matrix composite fabrication constraints, which was verified previously to possess a high-loss flow regime at high inlet turbulence conditions over a narrow range of Reynolds numbers and at a fixed design exit Mach number of 0.74. Our high-speed images, which were acquired at 10 distinct Reynolds numbers, show a significant increase in energy from flow oscillations due to transonic vortex shedding at Reynolds numbers corresponding to the high loss conditions. For those conditions, strong acoustic waves turn into shock waves. Spectral Proper Orthogonal Decomposition of the high-speed images shows acoustic waves from trailing edge vortex shedding at all conditions, with increased spectral energy at the high-loss conditions and slightly increasing frequency (about 6%) as a function of Reynolds number. Analysis of potential feedback timing is performed using velocity fields from a previous LES simulation, considering different feedback mechanisms. Most noteworthy is the acoustic/shock-boundary layer interaction mechanism on the suction surface at the blade geometric throat, which likely plays an important role in realistic curved blade passages.

Trailing Edge

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

In this work, we will discuss observations from images acquired from a time-resolved, high-speed self-aligned focusing Schlieren campaign that was performed at the CW-22 facility at NASA Glenn Research Center to understand the dynamic behavior of thick trailing-edge high-pressure turbine blades simulating a ceramic matrix composite (CMC) construction at high inlet turbulence conditions. For the CMC-9 blade, which has a trailing edge thickness of 9% of the axial chord, we identified a regime where an excessive total pressure loss (loss anomaly) is observed only for a narrow range of Reynolds numbers at a fixed exit Mach number of 0.74. The loss anomaly is qualitatively explained by our images, which were taken at 10 distinct blade Reynolds numbers spanning a factor of 6. The images show a significant increase in energy related to the oscillations due to transonic vortex shedding at the Reynolds numbers related to the high loss conditions. This increased energy leads to the formation of strong acoustic waves that turn into shock waves at the highest loss conditions. From our observations stemming from Spectral POD analysis of the high-speed images, we see the acoustic waves produced by the trailing edge vortex shedding exist in all conditions tested; but the shedding frequency has a very slight trend upwards as the Reynolds number is increased, varying about 6% in the range tested. Considering this variation of shedding frequency as a function of Reynolds number, which is well-established for other bluff-body flows, we stipulate there may be a potential feedback mechanism involving an acoustic information transfer path across neighboring blades that may explain why only a narrow range of Reynolds numbers displays strong, shock-forming vortex shedding. We consider a few feedback paths and examine the timing based on the mean flow field from a high-resolution LES simulation. It appears that all feedback mechanisms are viable, presenting an integer number of delay cycles with respect to disturbances generated at the trailing edge. Most noteworthy, however, is the acoustic/shock-boundary layer interaction mechanism at the blade geometric throat. Based on our analysis, this mechanism likely plays an important role in realistic, curved turbine blade passages.

Aerodynamics

Bosch Carbon Catalyst Regeneration in Water III: Experiments for and Related to Thermal Shock

Bosch technology has the potential to close the atmospheric revitalization loop by providing theoretically 100% regeneration of metabolic carbon dioxide into oxygen. The main obstacle to using Bosch technology is over-pressurization issues in the carbon formation reactor. Prior work on Bosch catalyst regeneration established that agitating carbon coated beads in water broke up pressure causing clogs and regenerated the catalyst surface. This process proved to be time consuming with experiments taking at least one hour to regenerate ten grams of carbon coated beads. In an attempt to shorten the processing time, thermally shocking the carbon coated beads prior to mechanical agitation was tested. Agitating the carbon coated beads after exposure to a saturated steam environment was also tested to ensure that the steam generated by thermally shocking the carbon coated beads would not cause adverse effects. The rationale behind thermally shocking the carbon coated beads is that the bead would pull away from the carbon shell due to a difference in the thermal expansion coefficients. The results of this project show that the processing time can be cut in half when compared to previous work, and that saturated steam has little to no effect on regeneration.

Bosch Process

Updates and Modernization of the Chemical Equilibrium with Applications (CEA) Code

NASA’s Chemical Equilibrium with Applications (CEA) code is a foundational tool for propulsion system analysis. It provides equilibrium chemistry, rocket performance, shock, and detonation calculations used across NASA and the broader aerospace community. NASA Engineering and Safety Center (NESC) Activity TI-22-01730 modernized the legacy CEA2 Fortran code into CEA v3, a Fortran 2008, object-oriented software package with expanded interface support, updated thermochemical data, improved maintainability, and substantially improved workflow integration. The modernized code preserves backward compatibility with legacy CEA input workflows while enabling direct use from modern analysis environments, including Python, C, MATLAB, and automated design studies.

Combustion

Experimental Demonstration of a Supersonic Spike Ejector

Spike diffusers are presented as a means to improve supersonic ejector efficiency via reduction in shock losses. Design, fabrication, and testing of an articulated prototype is described. Experimental results are compared to empirical performance curves representing the loci of optimal one- and two-stage conventional ejectors driven by dry air and discharging to atmosphere. Operated at a compression ratio of ~10 with all inputs equal, the spike ejector is shown to have achieved up to 2.19X and 1.83X the secondary mass throughput of the one- and two-stage references, respectively. This suggests the potential of spike diffusers to reduce ejector operating costs by ~50%.

wind tunnel

Bosch Carbon Catalyst Regeneration in Water I: Bosch Carbon Removal from Reacted Beads Utilizing Mechanical Agitation

The ability to recover oxygen from metabolic carbon dioxide is crucial for crewed long duration missions in which the atmosphere must be recycled. Two physiochemical methods are primarily employed for this purpose: the Bosch process and the Sabatier reaction. This paper focuses on the Bosch process. It utilizes hydrogen in the presence of a catalyst to form water and solid carbon out of metabolic carbon dioxide. The solid carbon product deposits onto the catalyst and eventually causes reactor over pressurization issues. Multiple methods are currently being investigated to regenerate the catalyst’s surface and break up the pressure causing carbon clogs. This paper explores whether water and mechanical abrasion would work to remove Bosch carbon by testing multiple washing conditions. The different washing conditions involved agitation intensity, heating, duration, and volume of water. It was found that it is possible to break up the clumps of carbon formed in the reactor through agitation with room temperature water in one hour and shows promise for completely cleaning the carbon off the beads’ surface.

Agitation

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics

Remote sensing of hurricane waves

A program was undertaken to obtain data on wind generated ocean waves in hurricanes by remote sensing techniques. The sensor, a synthetic aperture radar, has collected the first data ever on the directional wave climate throughout a hurricane. This information was found to have inconsistencies with present hurricane wave generation models.

David B King

Wave Rotor Enhanced Nuclear (WREN) Propulsion: NASA Innovative Advanced Concepts (NIAC) - Phase I Final Report

Nuclear Thermal Propulsion (NTP) is identified as one of the preferred propulsion technologies for manned missions throughout the solar system (NASA MSFC).[1, 2] The state-ofthe-art NTP cycle is based on a solid core Nuclear Engine for Rocket Vehicle Application (NERVA)[3] class technology (Fig. 1) that is envisioned to provide a specific impulse of 900 seconds doubling chemical rocket performance (450 seconds). Even with this impressive increase, the NTP NERVA designs still have issues providing adequate initial to final mass fractions for high ΔV missions.[4] Nuclear Electric Propulsion (NEP) can provide extremely high Isp (2,000 to over 10,000 seconds) but with only low thrust and limits on mass to power ratios. The need for an electric power source also adds the issue of heat rejection in space where thermal energy conversion is at best 30-40% under ideal conditions. NASA Space Technology Mission Directorate (STMD) has recently expressed interest in finding advanced nuclear propulsion technology through the NASA Go:Thrust RFI.[5, 6] A novel Wave Rotor (WR) topping cycle has been proposed for our NASA NIAC concept. It promises to deliver similar thrust as NERVA class NTP propulsion, but with Isp in the 1,200-2,000 second range. Coupled with an NEP cycle, the duty cycle Isp can further be increased (1,800-4,000 seconds) with minimal addition of dry mass. This bimodal design enables fast transit trajectories for manned missions to Mars and revolutionizes the deep space exploration of our solar system.

Nuclear Thermal Propulsion