Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock capturing”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Simulation of Acoustics for Ares I Scale Model Acoustic Tests

The Ares I Scale Model Acoustics Test (ASMAT) is a series of live-fire tests of scaled rocket motors meant to simulate the conditions of the Ares I launch configuration. These tests have provided a well documented set of high fidelity acoustic measurements useful for validation including data taken over a range of test conditions and containing phenomena like Ignition Over-Pressure and water suppression of acoustics. To take advantage of this data, a digital representation of the ASMAT test setup has been constructed and test firings of the motor have been simulated using the Loci/CHEM computational fluid dynamics software. Results from ASMAT simulations with the rocket in both held down and elevated configurations, as well as with and without water suppression have been compared to acoustic data collected from similar live-fire tests. Results of acoustic comparisons have shown good correlation with the amplitude and temporal shape of pressure features and reasonable spectral accuracy up to approximately 1000 Hz. Major plume and acoustic features have been well captured including the plume shock structure, the igniter pulse transient, and the ignition overpressure.

Putnam, Gabriel↗
NASA NTRS2011

Validation and Simulation of ARES I Scale Model Acoustic Test -1- Pathfinder Development

The Ares I Scale Model Acoustics Test (ASMAT) is a series of live-fire tests of scaled rocket motors meant to simulate the conditions of the Ares I launch configuration. These tests have provided a well documented set of high fidelity measurements useful for validation including data taken over a range of test conditions and containing phenomena like Ignition Over-Pressure and water suppression of acoustics. To take advantage of this data, a digital representation of the ASMAT test setup has been constructed and test firings of the motor have been simulated using the Loci/CHEM computational fluid dynamics software. Within this first of a series of papers, results from ASMAT simulations with the rocket in a held down configuration and without water suppression have then been compared to acoustic data collected from similar live-fire tests to assess the accuracy of the simulations. Detailed evaluations of the mesh features, mesh length scales relative to acoustic signals, Courant-Friedrichs-Lewy numbers, and spatial residual sources have been performed to support this assessment. Results of acoustic comparisons have shown good correlation with the amplitude and temporal shape of pressure features and reasonable spectral accuracy up to approximately 1000 Hz. Major plume and acoustic features have been well captured including the plume shock structure, the igniter pulse transient, and the ignition overpressure. Finally, acoustic propagation patterns illustrated a previously unconsidered issue of tower placement inline with the high intensity overpressure propagation path.

Putnam, G. C.↗
NASA NTRS1992

Eno-Osher schemes for Euler equations

The combination of the Osher approximate Riemann solver for the Euler equations and various ENO schemes is discussed for one-dimensional flow. The three basic approaches, viz. the ENO scheme using primitive variable reconstruction, either with Cauchy-Kowalewski procedure for time integration or the TVD Runge-Kutta scheme, and the flux-ENO method are tested on different shock tube cases. The shock tube cases were chosen to present a serious challenge to the ENO schemes in order to test their ability to capture flow discontinuities, such as shocks. Also the effect of the ordering of the eigen values, viz. natural or reversed ordering, in the Osher scheme is investigated. The ENO schemes are tested up to fifth order accuracy in space and time. The ENO-Osher scheme using the Cauchy-Kowalewski procedure for time integration is found to be the most accurate and robust compared with the other methods and is also computationally efficient. The tests showed that the ENO schemes perform reasonably well, but have problems in cases where two discontinuities are close together. In that case there are not enough points in the smooth part of the flow to create a non-oscillatory interpolation.

Vandervegt, Jacobus J.↗
NASA NTRS1986

Hypervelocity particle capture: Some considerations regarding suitable target media

Hypervelocity particles colliding with passive capture media will be traversed by shock waves; depending on the stress amplitude, the particle may remain solid or it may melt or vaporize. Any capture mechanism considered for cosmic dust collection in low Earth-orbit must be designed such that sample alteration and hence loss of scientific information is minimized. Capture of pristine particles is fundamentally difficult, because the specific heat of melting and even vaporization is exceeded upon impact at typical, geocentric encounter velocities. From the results of calculated and observed melting behaviors it is concluded that shock stresses in excess of 50 GPA should be avoided during hypervelocity particle capture on board Space Station and that stresses 20 GPa, even at 15 km/s collision velocities, should constitute desirable instrument design goals. Some principal characteristics of the capture medium that may satisfy these requirements are identified.

Horz, F.↗
NASA NTRS1998

Physics of Interplanetary Dust Collection with Aerogel

This report presents the results of research undertaken to study various problems associated with hypervelocity capture of dust particles in aerogel. The primary topics investigated were the properties of shocked aerogel and the requirements for reliable capture of particles on the STARDUST mission. In particular, the viscosity of shocked aerogel has been an open question. The results presented here suggest that the viscosity of aerogel at high impact velocities is negligible, although there remains some uncertainty about lower velocities. The model adopted for viscosity treats the mixture of polymeric silica and decomposition products and finds that, for particle velocities of 6-7 km/s, the viscosity is similar to that typical of light gasses at STP. Expressions for the Hugoniot of aerogel as a function of density were also obtained from the available data. All aerogels of interest for cosmic dust collectors have very similar shock velocity-particle velocity Hugoniot curves. The strength behavior of aerogel for low-speed penetration was measured, but further work is needed to study the proper way to apply this to the issue of terminal deceleration of a dust particle. Preliminary calculations designed to maximize the penetration depths were performed to determine the required density of aerogel to reliably stop a particle in a 3 cm thickness of aerogel (the path length expected for a normal impact into the STARDUST collector). In order to stop a particle of density rho(sub p) and diameter d(sub p), the mean density of the aerogel collector should be no less than that given by the expression bar rho(sub 0) = 1.085 X 10(exp -4 )rho(sub p)d(sub p), for densities measured in g/ cu cm and the particle diameter measured in micrometers.

Anderson, William W.↗
NASA NTRS2015

Cart3D Analysis of Plume and Shock Interaction Effects on Sonic Boom

A plume and shock interaction study was developed to collect data and perform CFD on a configuration where a nozzle plume passed through the shock generated from the wing or tail of a supersonic vehicle. The wing or tail was simulated with a wedge-shaped shock generator. Three configurations were analyzed consisting of two strut mounted wedges and one propulsion pod with an aft deck from a low boom vehicle concept. Research efforts at NASA were intended to enable future supersonic flight over land in the United States. Two of these efforts provided data for regulatory change and enabled design of low boom aircraft. Research has determined that sonic boom is a function of aircraft lift and volume distribution. Through careful tailoring of these variables, the sonic boom of concept vehicles has been reduced. One aspect of vehicle tailoring involved how the aircraft engine exhaust interacted with aft surfaces on a supersonic aircraft, such as the tail and wing trailing edges. In this work, results from Euler CFD simulations are compared to experimental data collected on sub-scale components in a wind tunnel. Three configurations are studied to simulate the nozzle plume interaction with representative wing and tail surfaces. Results demonstrate how the plume and tail shock structure moves with increasing nozzle pressure ratio. The CFD captures the main features of the plume and shock interaction. Differences are observed in the plume and deck shock structure that warrant further research and investigation.

supersonic wind tunnels↗
NASA NTRS1993

Collisions and close encounters between massive main-sequence stars

Smooth-particle hydrodynamics is used to investigate collisions and close encounters between two massive main-sequence stars. Approximate scaling relations and fitting formulas are given for the amount of mass loss and for the critical impact parameters for capture or merging. The more massive stars, which have smaller ratios of specific binding energy to the square of escape velocity, are more easily disrupted in collisions. As the relative velocity increases, the orbital energy of the colliding stars can only be dissipated by shock heating, and the critical capture radius decreases much faster than predicted by linear theory. The average fractional mass loss per collision in a cluster is found not to depend sensitively on the stellar velocity dispersion.

Lai, Dong↗
NASA NTRS2021

ARTEMIS Observations of Plasma Waves in Laminar and Perturbed Interplanetary Shocks

The ‘Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun’ (ARTEMIS) mission provides a unique opportunity to study the structure of interplanetary shocks and the associated generation of plasma waves with frequencies between ~50-8000 Hz due to its long duration electric and magnetic field burst waveform captures. We compare wave properties and occurrence rates at 11 quasi-perpendicular interplanetary shocks with burst data within 10 minutes (~3200 proton gyroradii upstream, ~1900 downstream) of the shock ramp. A perturbed shock is defined as possessing a large amplitude whistler precursor in the quasi-static magnetic field with an amplitude greater than ⅓ the difference between the upstream and downstream average magnetic field magnitudes; laminar shocks lack these large precursors and have a smooth, step function-like transition. In addition to wave modes previously observed, including ion acoustic, whistler, and electrostatic solitary waves, waves in the ion acoustic frequency range that show rapid temporal frequency change are common. Three shocks had burst captures in the ramp; of these, the two laminar shocks with burst data in the ramp contained a wide range of large amplitude wave modes in the ramp whereas the one perturbed shock contained no such waves. Thus, energy dissipation through wave-particle interactions is more prominent in these two laminar shocks than in the perturbed shock. Based on observations from all 11 shocks, the The wave occurrence rates for laminar shocks are higher in the transition region, especially the ramp, than downstream. In contrast, perturbed shocks have approximately 2-3 times the wave occurrence rate downstream than laminar shocks.

L. A. Davis↗
NASA NTRS1989

Vortex paradigm for shock-accelerated density-stratified interfaces

Numerical simulations are reported that capture the observations of Haas and Sturtevant (1985) when a planar shock accelerates an inclined planar interface between two media. A vortex paradigm is presented for interpreting the evolution of shock-accelerated density-stratified interfaces beyond early times. The simulations illustrate the paradigm and are in good agreement with the results of Haas and Sturtevant. Color images of space-time diagrams of the one-space integrated vorticity function are given to show the primary and secondary features of the emerging vortex structures.

Hawley, John F.↗
NASA NTRS2014

USM3D Predictions of Supersonic Nozzle Flow

This study focused on the NASA Tetrahedral Unstructured Software System CFD code (USM3D) capability to predict supersonic plume flow. Previous studies, published in 2004 and 2009, investigated USM3D's results versus historical experimental data. This current study continued that comparison however focusing on the use of the volume souring to capture the shear layers and internal shock structure of the plume. This study was conducted using two benchmark axisymmetric supersonic jet experimental data sets. The study showed that with the use of volume sourcing, USM3D was able to capture and model a jet plume's shear layer and internal shock structure.

Carter, Melissa B.↗
NASA NTRS2013

STEREO and Wind Observations of Intense Cyclotron Harmonic Waves at the Earth's Bow Shock and Inside the Magnetosheath

We present the first observations of electron cyclotron harmonic waves at the Earth's bow shock from STEREO and Wind burst waveform captures. These waves are observed at magnetic field gradients at a variety of shock geometries ranging from quasi-parallel to nearly perpendicular along with whistler mode waves, ion acoustic waves, and electrostatic solitary waves. Large amplitude cyclotron harmonic waveforms are also observed in the magnetosheath in association with magnetic field gradients convected past the bow shock. Amplitudes of the cyclotron harmonic waves range from a few tens to more than 500 millivolts/meter peak-peak. A comparison between the short (15 meters) and long (100 meters) Wind spin plane antennas shows a similar response at low harmonics and a stronger response on the short antenna at higher harmonics. This indicates that wavelengths are not significantly larger than 100 meters, consistent with the electron cyclotron radius. Waveforms are broadband and polarizations are distinctively comma-shaped with significant power both perpendicular and parallel to the magnetic field. Harmonics tend to be more prominent in the perpendicular directions. These observations indicate that the waves consist of a combination of perpendicular Bernstein waves and field-aligned waves without harmonics. A likely source is the electron cyclotron drift instability which is a coupling between Bernstein and ion acoustic waves. These waves are the most common type of high-frequency wave seen by STEREO during bow shock crossings and magnetosheath traversals and our observations suggest that they are an important component of the high-frequency turbulent spectrum in these regions.

bow shock↗
NASA NTRS2008

Transient Three-Dimensional Startup Side Load Analysis of a Regeneratively Cooled Nozzle

The objective of this effort is to develop a computational methodology to capture the startup side load physics and to anchor the computed aerodynamic side loads with the available data from a regeneratively cooled, high-aspect-ratio nozzle, hot-fired at sea level. The computational methodology is based on an unstructured-grid, pressure-based, reacting flow computational fluid dynamics and heat transfer formulation, a transient 5 s inlet history based on an engine system simulation, and a wall temperature distribution to reflect the effect of regenerative cooling. To understand the effect of regenerative wall cooling, two transient computations were performed using the boundary conditions of adiabatic and cooled walls, respectively. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with the pulsation of shocks across the lip, although the combustion wave is commonly eliminated with the sparklers during actual test. The test measured two side load events: a secondary and lower side load, followed by a primary and peak side load. Results from both wall boundary conditions captured the free-shock separation to restricted-shock separation transition with computed side loads matching the measured secondary side load. For the primary side load, the cooled wall transient produced restricted-shock pulsation across the nozzle lip with peak side load matching that of the test, while the adiabatic wall transient captured shock transitions and free-shock pulsation across the lip with computed peak side load 50% lower than that of the measurement. The computed dominant pulsation frequency of the cooled wall nozzle agrees with that of a separate test, while that of the adiabatic wall nozzle is more than 50% lower than that of the measurement. The computed teepee-like formation and the tangential motion of the shocks during lip pulsation also qualitatively agree with those of test observations. Moreover, a third transient computation was performed with a proportionately shortened 1 s sequence, and lower side loads were obtained with the higher ramp rate.

Wang, Ten-See↗
NASA NTRS2019

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗
NASA NTRS2015

On Theoretical Broadband Shock-Associated Noise Near-Field Cross-Spectra

The cross-spectral acoustic analogy is used to predict auto-spectra and cross-spectra of broadband shock-associated noise in the near-field and far-field from a range of heated and unheated supersonic off-design jets. A single equivalent source model is proposed for the near-field, mid-field, and far-field terms, that contains flow-field statistics of the shock wave shear layer interactions. Flow-field statistics are modeled based upon experimental observation and computational fluid dynamics solutions. An axisymmetric assumption is used to reduce the model to a closed-form equation involving a double summation over the equivalent source at each shock wave shear layer interaction. Predictions are compared with a wide variety of measurements at numerous jet Mach numbers and temperature ratios from multiple facilities. Auto-spectral predictions of broadband shock-associated noise in the near-field and far-field capture trends observed in measurement and other prediction theories. Predictions of spatial coherence of broadband shock-associated noise accurately capture the peak coherent intensity, frequency, and spectral width.

Miller, Steven A. E.↗
NASA NTRS2004

Experimental, Theoretical, and Computational Investigation of Separated Nozzle Flows

A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.

Hunter, Craig A.↗
NASA NTRS1985

Flow structure capturing on overset patched meshes

The present paper describes the application of a simple, robust stable implicit approach to solutions of the conservative equations of gasdynamics on either composite or overset meshes to capture flow structures such as shocks with high resolution in complex geometric domains. Without requiring special flux conservative operators, but rather, interpolating conservative variable data at mesh boundaries, the well posed upwind method provides accurate and oscillation free solutions, even when shocks cross the interior patch boundaries. In three problems with flow complexities that require grid refinement, the paper demonstrates the capability to conveniently carry out for gasdynamics the adaptive refined meshing strategy in overset patches proposed by Berger and Oliger (1984), and it extends this technique to rapidly convergent implicit methods for the Euler and Navier-Stokes equations. The numerical experiments show concretely, in a realistic aerodynamic problem, the savings in mesh points (about an order of magnitude here in two dimensions) for similar accuracy that flow structure aligned adaptive patched meshing affords compared to uniform grid refinement.

Venkatapathy, E.↗
NASA NTRS1986

An upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method does not require the addition of user specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming scheme in terms of accuracy, stability, computer time and storage, and programming effort. The new algorithm has been validated by applying it to three laminar test cases including flat plate boundary-layer flow, hypersonic flow past a 15 deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with the results obtained using the conventional Beam-Warming algorithm.

Lawrence, S. L.↗