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 415 records · Page 23

Wind-tunnel preflight test program for aeroassist flight experiment

A proposed Aeroassist Flight Experiment (AFE) will simulate return from geosynchronous orbit to low-earth orbit to provide fundamental information about the hypersonic flow conditions surrounding aeroassisted orbital transfer vehicles (AOTV's). Future AOTV design must rely heavily on computational fluid dynamic (CFD) computer codes now under development, since the entry conditions are beyond the simulation capabilities of ground-based test facilities. However, existing codes and best available ground-based test facilities must provide for design of the AFE. The AFE ground-based test program provides data for calibration of existing and developing pertinent CFD codes. Confidence in code prediction capability increases with ability to predict forces and moments, pressure and heat-transfer distributions, shock shapes, and surface streamline directions. The ground-based test program provides these data over a wide range of hypersonic test parameters. This paper describes the test program, models, facilities, and some representative results.

Wells, William L.↗

Two-equation turbulence modeling for 3-D hypersonic flows

An investigation to verify, incorporate and develop two-equation turbulence models for three-dimensional high speed flows is presented. The current design effort of hypersonic vehicles has led to an intensive study of turbulence models for compressible hypersonic flows. This research complements an extensive review of experimental data and the current development of 2D turbulence models. The review of experimental data on 2D and 3D flows includes complex hypersonic flows with pressure profiles, skin friction, wall heat transfer, and turbulence statistics data. In a parallel effort, turbulence models for high speed flows have been tested against flat plate boundary layers, and are being tested against the 2D database. In the present paper, we present the results of 3D Navier-Stokes numerical simulations with an improved k-omega two-equation turbulence model against experimental data and empirical correlations of an adiabatic flat plate boundary layer, a cold wall flat plate boundary layer, and a 3D database flow, the interaction of an oblique shock wave and a thick turbulent boundary layer with a free stream Mach number = 8.18 and Reynolds number = 5 x 10 to the 6th.

Bardina, J. E.↗

Numerical simulation by TVD schemes of complex shock reflections from airfoils at high angle of attack

The shock-capturing capability of total variation diminishing (TVD) schemes is demonstrated for a more realistic complex shock-diffraction problem for which the experimental data are available. Second-order explicit upwind and symmetric TVD schemes are used to solve the time-dependent Euler equations of gas dynamics for the interaction of a blast wave with an airfoil at high angle-of-attack. The test cases considered are a time-dependent moving curved-shock wave and a contant moving planar-shock wave impinging at an angle-of-attack 30 deg on a NACA 0018 airfoil. Good agreement is obtained between isopycnic contours computed by the TVD schemes and those from experimental interferograms. No drastic difference in flow-field structure is found between the curved- and planar-shock wave cases, except for a difference in density level near the lower surface of the airfoil. Computation for cases with higher shock Mach numbers is also possible. Numerical experiments show that the symmetric TVD scheme is less sensitive to the boundary conditions treatment than the upwind scheme.

Moon, Young J.↗

Sierra/SD – Verification Test Manual – 5.22

Verification and validation (V&V) of scientific computing programs are important at Sandia National Labs due to the expanding role of computational simulation in managing the United States nuclear stockpile. The complexities of structural response calculations used to analyze physical problems, the varieties of codes applied to the calculations, and the importance of accurate predictions when assessing field conditions demand confidence in the consistency and accuracy of computer codes. Confidence in the accuracy of the predictions arising from computer simulations must ultimately be gained through verification and validation. The Sierra salinas structural dynamics analysis code, Sierra/SD, is used at the DOE Laboratories, and in several DOD projects. The roles of Sierra/SD in the qualification of weapon systems and components for normal and hostile environments throughout the Stockpile-to-Target Sequence include to, • Redesign weapon components. • Certify weapon components and systems for target environments such as hypersonic vehicles. • Certify that components will survive the thermal mechanical shock loads associated with hostile environments. • Evaluate current stockpile issues, including issues associated with uncertainty quantification. • Address many other problems that are encountered in stockpile management. The Sierra/SD verification plan is described, and an evolving set of key verification tests are described in detail. The verification tests ensure the correctness of the mathematics and numerical algorithms associated with functionality describing engineering phenomena. Development is in accordance with a set of tailored Software Quality Engineering (SQE) practices. SQE practices guide the overall verification and validation effort.

97 MATHEMATICS AND COMPUTING↗

NEQAIRv14.0 Release Notes: Nonequilibrium and Equilibrium Radiative Transport Spectra Program

NEQAIR v14.0 is the first parallelized version of NEQAIR. Starting from the last version of the code that went through the internal software release process at NASA Ames (NEQAIR 2008), there have been significant updates to the physics in the code and the computational efficiency. NEQAIR v14.0 supersedes NEQAIR v13.2, v13.1 and the suite of NEQAIR2009 versions. These updates have predominantly been performed by Brett Cruden and Aaron Brandis from ERC Inc at NASA Ames Research Center in 2013 and 2014. A new naming convention is being adopted with this current release. The current and future versions of the code will be named NEQAIR vY.X. The Y will refer to a major release increment. Minor revisions and update releases will involve incrementing X. This is to keep NEQAIR more in line with common software release practices. NEQAIR v14.0 is a standalone software tool for line-by-line spectral computation of radiative intensities and/or radiative heat flux, with one-dimensional transport of radiation. In order to accomplish this, NEQAIR v14.0, as in previous versions, requires the specification of distances (in cm), temperatures (in K) and number densities (in parts/cc) of constituent species along lines of sight. Therefore, it is assumed that flow quantities have been extracted from flow fields computed using other tools, such as CFD codes like DPLR or LAURA, and that lines of sight have been constructed and written out in the format required by NEQAIR v14.0. There are two principal modes for running NEQAIR v14.0. In the first mode NEQAIR v14.0 is used as a tool for creating synthetic spectra of any desired resolution (including convolution with a specified instrument/slit function). The first mode is typically exercised in simulating/interpreting spectroscopic measurements of different sources (e.g. shock tube data, plasma torches, etc.). In the second mode, NEQAIR v14.0 is used as a radiative heat flux prediction tool for flight projects. Correspondingly, NEQAIR has also been used to simulate the radiance measured on previous flight missions. This report summarizes the database updates, corrections that have been made to the code, changes to input files, parallelization, the current usage recommendations, including test cases, and an indication of the performance enhancements achieved.

Radiation Solver↗

NEQAIR v15.0 Release Notes: Nonequilibrium and Equilibrium Radiative Transport and Spectra Program

NEQAIR v15.0 provides the first steps to improved coupling between NEQAIR and the DPLR CFD code, which will be fully realized in v15.1. The plan is to release NEQAIR v15.1 and DPLR 4.05 at the same time. The improvements implemented in NEQAIR v15.0 have focused on improving stability, solution robustness, usability and providing different options for running the code. It is also the first version of the code to have a new input file and line of sight format since 2009. Backward compatibility with previous formats of the input files (neqair.inp and LOS.dat) has also been provided. NEQAIR v15.0 supersedes the prerelease of this version, as well as NEQAIR v14.0, v13.2, v13.1 and the suite of NEQAIR2009 versions. These updates have predominantly been performed by Brett Cruden and Aaron Brandis from AMA Inc at NASA Ames Research Center between 2016 and 2018. NEQAIR v15.0 is a standalone software tool for line-by-line spectral computation of radiative intensities and/or radiative heat flux, with one-dimensional transport of radiation. In order to accomplish this, NEQAIR v15.0, as in previous versions, requires the specification of distances (in cm), temperatures (in K) and number densities (in parts/cc) of constituent species along lines of sight. Therefore, it is assumed that flow quantities have been extracted from flow fields computed using other tools, such as CFD codes like DPLR or LAURA, and that lines of sight have been constructed and written out in the format required by NEQAIR v15.0. There are two principal modes for running NEQAIR v15.0. In the first mode NEQAIR v15.0 is used as a tool for creating synthetic spectra of any desired resolution (including convolution with a specified instrument/slit function). The first mode is typically exercised in simulating/interpreting spectroscopic measurements of different sources (e.g. shock tube data, plasma torches, etc.). In the second mode, NEQAIR v15.0 is used as a radiative heat flux prediction tool for flight projects. Correspondingly, NEQAIR has also been used to simulate the radiance measured on previous flight missions. This report summarizes the database updates, corrections that have been made to the code, changes to input files, parallelization, the current usage recommendations, including test cases, and an indication of the performance enhancements achieved.

Brandis, Aaron M.↗

GenASiS: General Astrophysical Simulation System. II. Self-gravitating Baryonic Matter*

GenASiS (General Astrophysical Simulation System) is a code being developed initially and primarily, though not exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the second in a series---documents capabilities for Newtonian self-gravitating fluid dynamics, including tabulated microphysical equations of state treating nuclei and nuclear matter (`baryonic matter'). Computation of the gravitational potential of a spheroid, and simulation of the gravitational collapse of dust and of an ideal fluid, provide tests of self-gravitation against known solutions. In multidimensional computations of the adiabatic collapse, bounce, and explosion of spherically symmetric pre-supernova progenitors---which we propose become a standard benchmark for code comparisons---we find that the explosions are prompt and remain spherically symmetric (as expected), with an average shock expansion speed and total kinetic energy that are inversely correlated with the progenitor mass at the onset of collapse and the compactness parameter.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗

Numerical simulation of turbulent trailing edge flows

Numerical simulations of the time-dependent, Reynolds-averaged, Navier-Stokes equations, employing a two-equation turbulence model, are presented and compared with measurements from a series of trailing edge experiments at transonic Mach numbers. The test flows include an asymmetric flow with no separation, an asymmetric flow with a small region of separation and a symmetric flow with a large shock-wave induced separated zone. Comparisons are made for mean surface quantities as well as for mean and fluctuating flow-field quantities. For the trailing-edge flows with little or no separation, the solutions correctly predict all the major features of the flow field. Treatment of the viscous-inviscid interaction was found to be important for predicting these test cases. Two-equation eddy-viscosity turbulence models were found to be adequate for these flows. However, for the shock-wave induced separation case, these turbulence models were inadequate to predict this flow field. Modifications of the turbulence model to correct these deficiencies are discussed.

Horstman, C. C.↗

Evaluation of Industry Standard Turbulence Models on an Axisymmetric Supersonic Compression Corner

Reynolds-averaged Navier-Stokes computations of a shock-wave/boundary-layer interaction (SWBLI) created by a Mach 2.85 flow over an axisymmetric 30-degree compression corner were carried out. The objectives were to evaluate four turbulence models commonly used in industry, for SWBLIs, and to evaluate the suitability of this test case for use in further turbulence model benchmarking. The Spalart-Allmaras model, Menter's Baseline and Shear Stress Transport models, and a low-Reynolds number k- model were evaluated. Results indicate that the models do not accurately predict the separation location; with the SST model predicting the separation onset too early and the other models predicting the onset too late. Overall the Spalart-Allmaras model did the best job in matching the experimental data. However there is significant room for improvement, most notably in the prediction of the turbulent shear stress. Density data showed that the simulations did not accurately predict the thermal boundary layer upstream of the SWBLI. The effect of turbulent Prandtl number and wall temperature were studied in an attempt to improve this prediction and understand their effects on the interaction. The data showed that both parameters can significantly affect the separation size and location, but did not improve the agreement with the experiment. This case proved challenging to compute and should provide a good test for future turbulence modeling work.

Turbulent Boundary Layer↗

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗

Numerical simulation of transient hypervelocity flow in an expansion tube

Several numerical simulations of the transient flow of helium in an expansion tube are presented in an effort to identify some of the basic mechanisms which cause the noisy test flows seen in experiments. The calculations were performed with an axisymmetric Navier-Stokes code based on a finite volume formulation and upwinding techniques. Although laminar flow and ideal bursting of the diaphragms was assumed, the simulations showed some of the important features seen in experiments. In particular, the discontinuity in tube diameter of the primary diaphragm station introduced a transverse perturbation to the expanding driver gas and this perturbation was seen to propagate into the test gas under some flow conditions. The disturbances seen in the test flow can be characterized as either small amplitude, low frequency noise possibly introduced during shock compression or large amplitude, high frequency noise associated with the passage of the reflected head of the unsteady expansion.

Jacobs, P. A.↗

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

In the last two decades, there have been extensive developments in computational aerodynamics, which constitutes a major part of the general area of computational fluid dynamics. Such developments are essential to advance the understanding of the physics of complex flows, to complement expensive wind-tunnel tests, and to reduce the overall design cost of an aircraft, particularly in the area of aeroelasticity. Aeroelasticity plays an important role in the design and development of aircraft, particularly modern aircraft, which tend to be more flexible. Several phenomena that can be dangerous and limit the performance of an aircraft occur because of the interaction of the flow with flexible components. For example, an aircraft with highly swept wings may experience vortex-induced aeroelastic oscillations. Also, undesirable aeroelastic phenomena due to the presence and movement of shock waves occur in the transonic range. Aeroelastically critical phenomena, such as a low transonic flutter speed, have been known to occur through limited wind-tunnel tests and flight tests. Aeroelastic tests require extensive cost and risk. An aeroelastic wind-tunnel experiment is an order of magnitude more expensive than a parallel experiment involving only aerodynamics. By complementing the wind-tunnel experiments with numerical simulations the overall cost of the development of aircraft can be considerably reduced. In order to accurately compute aeroelastic phenomenon it is necessary to solve the unsteady Euler/Navier-Stokes equations simultaneously with the structural equations of motion. These equations accurately describe the flow phenomena for aeroelastic applications. At Ames a code, ENSAERO, is being developed for computing the unsteady aerodynamics and aeroelasticity of aircraft and it solves the Euler/Navier-Stokes equations. The purpose of this contract is to continue the algorithm enhancements of ENSAERO and to apply the code to complicated geometries. During the last year, the geometric capability of the code was extended to simulate transonic flows, a wing with oscillating control surface. Single-grid and zonal approaches were tested. For the zonal approach, a new interpolation technique was introduced. The key development of the algorithm was an interface treatment between moving zones for a control surface using the virtual-zone concept. The work performed during the period, 1 Apr. 1992 through 31 Mar. 1993 is summarized. Additional details on the various aspects of the study are given in the Appendices.

Obayashi, Shigeru↗

Challenging Pneumatic Requirements for Acoustic Testing of the Cryogenic Second Stage for the New Delta 3 Rocket

The paper describes the unique pneumatic test requirements for the acoustic and shock separation testing of the Second Stage for the new Delta III Rocket at the Goddard Space Flight Center in Greenbelt, Maryland. The testing was conducted in the 45,000 cu ft (25-feet wide by 30-feet deep by 50-foot high) Acoustic Facility. The acoustic testing required that the liquid oxygen (LOX) and liquid hydrogen (LH2) tanks be filled with enough liquid nitrogen (LN2) to simulate launch fuel masses during testing. The challenge for this test dealt with designing, procuring, and fabricating the pneumatic supply systems for quick assembly while maintaining the purity requirements and minimizing costs. The pneumatic systems were designed to fill and drain the both LOX and LH2 tanks as well as to operate the fill/drain and vent valves for each of the tanks. The test criteria for the pneumatic sub-systems consisted of function, cleanliness, availability, and cost. The first criteria, function, required the tanks to be filled and drained in an efficient manner while preventing them from seeing pressures greater than 9 psig which would add a pressure cycle to the tank. An LN2 tanker, borrowed from another NASA facility, served as the pre-cool and drain tanker. Pre-cooling the tanks allowed for more efficient and cost effective transfer from the LN2 delivery tankers. Helium gas, supplied from a high purity tube trailer, was used to pressurize the vapor space above the LN2 pushing it into the drain tanker. The tube trailer also supplied high pressure helium to the vehicle for valve control and component purges. Cleanliness was maintained by proper component selection, end-use particle filtration, and any on-site cleaning determined necessary by testing. In order to meet the availability/cost juggling act, products designed for LOX delivery systems were procured to ensure system compatibility while off the shelf valves and tubing designed for the semiconductor industry were procured for the gas systems.

Webb, Andrew T.↗

Assessment of Near-Field Sonic Boom Simulation Tools

A recent study for the Supersonics Project, within the National Aeronautics and Space Administration, has been conducted to assess current in-house capabilities for the prediction of near-field sonic boom. Such capabilities are required to simulate the highly nonlinear flow near an aircraft, wherein a sonic-boom signature is generated. There are many available computational fluid dynamics codes that could be used to provide the near-field flow for a sonic boom calculation. However, such codes have typically been developed for applications involving aerodynamic configuration, for which an efficiently generated computational mesh is usually not optimum for a sonic boom prediction. Preliminary guidelines are suggested to characterize a state-of-the-art sonic boom prediction methodology. The available simulation tools that are best suited to incorporate into that methodology are identified; preliminary test cases are presented in support of the selection. During this phase of process definition and tool selection, parallel research was conducted in an attempt to establish criteria that link the properties of a computational mesh to the accuracy of a sonic boom prediction. Such properties include sufficient grid density near shocks and within the zone of influence, which are achieved by adaptation and mesh refinement strategies. Prediction accuracy is validated by comparison with wind tunnel data.

Casper, J. H.↗

Integration Assessment of Visiting Vehicle Induced Electrical Charging of the International Space Station Structure

The International Space Station (ISS) undergoes electrical charging in low Earth orbit (LEO) due to positively biased, exposed conductors on solar arrays that collect electrical charges from the space plasma. Exposed solar array conductors predominately collect negatively charged electrons and thus drive the metal ISS structure electrical ground to a negative floating potential (FP) relative to plasma. This FP is variable in location and time as a result of local ionospheric conditions. ISS motion through Earth s magnetic field creates an addition inductive voltage up to 20 positive and negative volts across ISS structure depending on its attitude and location in orbit. ISS Visiting Vehicles (VVs), such as the planned Orion crew exploration vehicle, contribute to the ISS plasma charging processes. Upon physical contact with ISS, the current collection properties of VVs combine with ISS. This is an ISS integration concern as FP must be controlled to minimize arcing of ISS surfaces and ensure proper management of extra vehicular activity crewman shock hazards. This report is an assessment of ISS induced charging from docked Orion vehicles employing negatively grounded, 130 volt class, UltraFlex (ATK Space Systems) solar arrays. To assess plasma electron current collection characteristics, Orion solar cell test coupons were constructed and subjected to plasma chamber current collection measurements. During these tests, coupon solar cells were biased between 0 and 120 V while immersed in a simulated LEO plasma. Tests were performed using several different simulated LEO plasma densities and temperatures. These data and associated theoretical scaling of plasma properties, were combined in a numerical model which was integrated into the Boeing Plasma Interaction Model. It was found that the solar array design for Orion will not affect the ISS FP by more than about 2 V during worst case charging conditions. This assessment also motivated a trade study to determine acceptable plasma electron current levels that can be collected by a single or combined fleet of ISS-docked VVs.

Kramer, Leonard↗

Examination of Mars2020 Shock-Layer Conditions via Infrared Laser Absorption Spectroscopy of CO2 and CO

Amid-infrared laser absorption diagnostic was deployed to study a simulated Mars2020 shock layer in the Electric Arc Shock Tube (EAST) facility at NASA’s Ames Research Center. Rapid RF-diplexing techniques enabled quantitative temperature and number density measurements of CO2 and CO with μs-resolution over an incident shock velocity range of 1.39 – 3.75 km/s. Two interband cascade lasers were utilized at 4.17 and 4.19 μm to resolve eight CO2 asymmetric stretch fundamental band (𝜈3) transitions from two different vibrational levels: 0000 (ground state) and 0110 (first excited bending mode). The probed rotational levels span across J” = 58 to J” = 140. Results are compared to DPLR simulations of the shock layer using kinetic mechanisms of Johnston et al. and Cruden et al. At shock velocities below 3.1 km/s, the agreement between the measurements and the Johnston mechanism is typically within 5% for temperature and within 10% for number density. At shock velocities above 3.1 km/s, the CO2 measurement becomes sensitive to a thin boundary layer and corrections of this effect are presented. On test cases with enough energy to dissociate CO2, a quantum cascade laser scanned the P(2, 20), P(0, 31), and P(3, 14) transitions of the CO fundamental band at 4.98 μm. CO formation rate is measured to be close to the Johnston kinetic mechanism at low velocities, and then trending towards the Cruden kinetic mechanism at high velocities. On a few low velocity test cases, rovibrational relaxation of the Martian atmosphere is probed with μs resolution.

Laser Absorption↗

Development of three-dimensional code for the analysis of jet mixing problem. Part 1: Laminar solution

Future aircraft will eventually feature nonaxisymmetric or rectangular nozzles. Developing a three-dimensional code to stimulate the characteristics of the jet exhaust plume, issuing from nonaxisymmetric nozzles, in general, at different flight conditions, is very important. Two three-dimensional codes were developed to simulate the shock-cell structure of circular nozzles. These codes were developed to solve the parabolized and simplified Navier-Stokes equations respectively. Both codes are based on a method previously developed by Newsome et al. These codes are fully vectorized on the VPS 32 at NASA Langley Research Center. The axisymmetric underexpanded supersonic jet flow problem, exhausting into still air, was used as a test case for developing an efficient three-dimensional problems and preserving crossplane symmetry of the flow downstream of the jet exit.

Abdol-Hamid, Khaled S.↗

Numerical simulation of transient hypervelocity flow in an expansion tube

Several numerical simulations of the transient flow of helium in an expansion tube are presented. The aim of the exercise is to provide further information on the operational problems of the NASA Langley expansion tube. The calculations were performed with an axisymmetric Navier-Stokes code based on a finite-volume formulation and upwinding techniques. Although laminar flow and ideal bursting of the diaphragms was assumed, the simulations showed some of the important features seen in the experiments. In particular, the discontinuity in the tube diameter at the primary diaphragm station introduced a transverse perturbation to the expanding driver gas, and this perturbation was seen to propagate into the test gas under some flow conditions. The disturbances seen in the test flow can be characterized as either 'small-amplitude' noise possibly introduced during shock compression or 'large-amplitude' noise associated with the passage of the reflected head of the unsteady expansion.

Jacobs, P. A.↗