Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock simulation test”

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 145 records · Page 8

A Simulator for Neyer Tests of Explosives

Explosives and explosive devices such as detonators are typically tested by applying a range of stimuli such as voltage or mechanical shock, and recording binary “detonated/did not detonate” responses. These are analyzed using maximum likelihood or generalized linear models to provide estimates of quantities such as the all-fire and no-fire points. Given that the true threshold for detonation is unknown a priori , sequential design methods are typically used to optimize the set of test points. One popular method, implemented in commercial software, is Neyer’s algorithm. To support simulation and experimental design, we have developed code in the R programming language to duplicate the functions of the Neyer software. We provide code for the simulator along with a description and examples of usage.

42 ENGINEERING↗

Application of Laser Induced Fluorescence for flow diagnostics in arc jets

Laser Induced Fluorescence (LIF) studies are carried out on copper atoms, nitric oxide and oxygen molecules in the arc jet flows at the NASA Johnson Space Center Reentry Simulation Testing Facility. Measurements are taken in the free stream and from a blunt body shock layer. Tests are performed under different flow conditions to determine the feasibility and sensitivity of the LIF technique for various species. The LIF technique is developed as a part of high enthalpy flow diagnostics and will be useful to elucidate the rotational and vibrational temperatures. Adequate sensitivity for the detection of O2 and NO is demonstrated. Preliminary experiments for temperature and velocity measurements using copper as a tracer material are shown to be a useful new application of the laser resonant Doppler scattering technique.

Arepalli, Sivaram↗

Simulation of transient flow in a shock tunnel and a high Mach number nozzle

A finite volume Navier-Stokes code was used to simulate the shock reflection and nozzle starting processes in an axisymmetric shock tube and a high Mach number nozzle. The simulated nozzle starting processes were found to match the classical quasi-1-D theory and some features of the experimental measurements. The shock reflection simulation illustrated a new mechanism for the driver gas contamination of the stagnated test gas.

Jacobs, P. A.↗

Application of the space-time conservation element and solution element method to shock-tube problem

An Euler solver based on the method of space-time conservation element and solution element is in this paper to simulate shock-tube flows involving shock waves, contact discontinuities, expansion waves and their intersections. Seven test problems are considered to examine the capability of this method. The numerical results, when compared with exact solutions and/or numerical solutions by other methods, indicate that the present method can accurately resolve strong shock and contact discontinuities without using any ad hoc techniques which are used only at the neighborhood of a discontinuity.

Wang, Xiao-Yen↗

Boundary Layer Protuberance Simulations in Channel Nozzle Arc-Jet

Two protuberance designs were modeled in the channel nozzle of the NASA Johnson Space Center Atmospheric Reentry Materials and Structures Facility with the Data-Parallel Line Relaxation computational fluid dynamics code. The heating on the protuberance was compared to nominal baseline heating at a single fixed arc-jet condition in order to obtain heating augmentation factors for flight traceability in the Boundary Layer Transition Flight Experiment on Space Shuttle Orbiter flights STS-119 and STS-128. The arc-jet simulations were performed in conjunction with the actual ground tests performed on the protuberances. The arc-jet simulations included non-uniform inflow conditions based on the current best practices methodology and used variable enthalpy and constant mass flow rate across the throat. Channel walls were modeled as fully catalytic isothermal surfaces, while the test section (consisting of Reaction Cured Glass tiles) was modeled as a partially catalytic radiative equilibrium wall. The results of the protuberance and baseline simulations were compared to the applicable ground test results, and the effects of the protuberance shock on the opposite channel wall were investigated.

Marichalar, J. J.↗

Asymmetric Uncertainty Expression for High Gradient Aerodynamics

When the physics of the flow around an aircraft changes very abruptly either in time or space (e.g., flow separation/reattachment, boundary layer transition, unsteadiness, shocks, etc), the measurements that are performed in a simulated environment like a wind tunnel test or a computational simulation will most likely incorrectly predict the exact location of where (or when) the change in physics happens. There are many reasons for this, includ- ing the error introduced by simulating a real system at a smaller scale and at non-ideal conditions, or the error due to turbulence models in a computational simulation. The un- certainty analysis principles that have been developed and are being implemented today do not fully account for uncertainty in the knowledge of the location of abrupt physics changes or sharp gradients, leading to a potentially underestimated uncertainty in those areas. To address this problem, a new asymmetric aerodynamic uncertainty expression containing an extra term to account for a phase-uncertainty, the magnitude of which is emphasized in the high-gradient aerodynamic regions is proposed in this paper. Additionally, based on previous work, a method for dispersing aerodynamic data within asymmetric uncer- tainty bounds in a more realistic way has been developed for use within Monte Carlo-type analyses.

Pinier, Jeremy T↗

Formation, Survival and Evolution of Organics After Impact Processing of Asteroidal and Cometary Simulants

This work examines the formation and evolution of amino acids during impact events on the prebiotic Earth Laboratory experiments, which simulate impact processes, provide direct tests of questions regarding the conditions (concentration of amino acid, shock stress or impact velocity, matrix material) under which amino acid synthesis and polymerization occur.

exobiology↗

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗

USM3D Simulations for Third Sonic Boom Workshop

The NASA USM3D flow solver was used to compute test cases for the Third AIAA Sonic Boom Prediction Workshop (SBPW3). The test cases include an axisymmetric equivalent area body, a near field biconvex shock-plume interaction wind tunnel model, and the C608 Low Boom Flight Demonstrator. Numerical simulations were conducted on the mixed element grids and the tetrahedral grids provided by the workshop committee, as well as a family of grids generated by an in-house approach for sonic boom analyses known as Boom Grid. The near-field pressure signatures were extracted, propagated to the ground and the perceived loudness levels on the ground was computed. The USM3D near-field pressure signatures,corresponding ground signatures, and loudness levels on the ground are compared with mean values from other workshop participants.

Supersonics↗

Advanced Shock Position Control for Mode Transition in a Turbine Based Combined Cycle Engine Inlet Model

A dual flow-path inlet system is being tested to evaluate methodologies for a Turbine Based Combined Cycle (TBCC) propulsion system to perform a controlled inlet mode transition. Prior to experimental testing, simulation models are used to test, debug, and validate potential control algorithms. One simulation package being used for testing is the High Mach Transient Engine Cycle Code simulation, known as HiTECC. This paper discusses the closed loop control system, which utilizes a shock location sensor to improve inlet performance and operability. Even though the shock location feedback has a coarse resolution, the feedback allows for a reduction in steady state error and, in some cases, better performance than with previous proposed pressure ratio based methods. This paper demonstrates the design and benefit with the implementation of a proportional-integral controller, an H-Infinity based controller, and a disturbance observer based controller.

Csank, Jeffrey T.↗

Experimental program for real gas flow code validation at NASA Ames Research Center

The experimental program for validating real gas hypersonic flow codes at NASA Ames Rsearch Center is described. Ground-based test facilities used include ballistic ranges, shock tubes and shock tunnels, arc jet facilities and heated-air hypersonic wind tunnels. Also included are large-scale computer systems for kinetic theory simulations and benchmark code solutions. Flight tests consist of the Aeroassist Flight Experiment, the Space Shuttle, Project Fire 2, and planetary probes such as Galileo, Pioneer Venus, and PAET.

Deiwert, George S.↗

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗

Design and subscale tests of a diffuser system for a Mach 4 scramjet test facility

A fixed-geometry diffuser system was designed for use in Mach 4 free-jet tests of a hydrogen-burning modular scramjet engine. The scramjet engine has a rectangular cross section that blocks up to 33 percent of the existing tunnel nozzle exit area and swept leading edges of the scramjet produce an asymmetric downward flow that sharply increases when the scramjet inlet unstarts. Effects of these conditions on the operation of the diffuser system are not accurately predictable. An experimental investigation has therefore been conducted in unheated air using a subscale model of the tunnel-scramjet-diffuser system to substantiate the design. Test results showed that the preliminary design had to be modified to obtain an acceptable configuration of the nozzle exit, test cabin, and diffuser that would provide shock-free flow at the scramjet inlet for simulated Mach 4 flight at altitudes from 16,764 to 20,422 meters.

Andrews, E. H., Jr.↗

Numerical models of two complex hypersonic flowfields

Numerical solutions of the compressible Navier-Stokes equations assuming a perfect gas and laminar flow are applied to two hypersonic experiments conducted at NASA Langley Research Center. The first is cooling by injection of a gas jet through the nose of an ogive-cone, and the second is the aerothermal environment in the gap formed by the wing and elevon section of a test model of the Space Shuttle. The simulations demonstrate that the model predicts accurate pressures and shock stand-off for the gas jet and accurate pressures for the wing-elevon cove. Experimental heating rates for the gas jet indicate turbulent mixing to be occurring, an effect not modeled in the numerical solution.

Macaraeg, M. G.↗

Characterization of Titan Entry Radiative Heating in the Low Density Electric Arc Shock Tube

Prediction of aftbody radiative heating for Dragonfly’s entry at Titan requires accurate models to inform the design of its thermal protection systems. Previously reported investigations of shock layer aftbody radiation in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center focused on the lower altitude/higher density part of the trajectory. This paper presents the first optical emission spectroscopy measurements of an ongoing test series in the new aluminum Low Density Shock Tube (LDST), recently installed at EAST, focusing on the radiative heating characterization through the higher altitude portion of Dragonfly’s entry path. Experimental data collected at a fill pressure of 0.1 Torr and shock velocities between 5.2 and 6.2 km/s show favorable agreement to similar conditions previously tested in the 4" High-Velocity Shock Tube, as well as improvements in the signal to noise ratio and available test time. Simulations employing DPLR and NEQAIR underpredict the experimental radiance, with increasing discrepancy from the NIR towards the UV wavelength range, thus motivating further efforts in improving the current non-Boltzmann model for the CN molecule.

A Fagnani↗

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗