Computational challenges for simulations related to the NASA electric arc shock tube (EAST) experiments

The goal of this study is to gain some physical insights and an understanding of the computational challenges for the simulations related to the hypersonic nonequilibrium multi-species and multi-reaction experiments on the NASA Electric Arc Shock Tube (EAST). While experimental measurement does not provide any information about the radial structure of this type of flow, accurate and reliable numerical simulations can provide more insight into the physical structure of the flow to aid the design of atmospheric entry spacecrafts. The paper focuses on the spurious numerics which take place in numerical simulations of the subject physics containing stiff source terms and discontinuities. This paper is based on the knowledge gained from Yee et al. on simple reacting test cases (Yee et al. 2013, [9]) as a guide to reveal the computational challenges involved for such an extreme flow type. The results of the 1D and 2D EAST viscous and inviscid simulations using a simplified physical model are presented. The computation reveals, for the first time, that the 2D viscous model which contains both shocks and shears exhibits Tollmien– Schlichting-like instability complex patterns at the boundary layer. In addition to exhibiting spurious numerical behavior of wrong propagation speed of discontinuities by typical shock-capturing methods, there is improved understanding on the cause of numerical difficulties by previous investigators. One example is that the relative distance between the shocks and shear/contact is different from one grid spacing to another for each considered high order shock-capturing scheme. The results presented can provide insight on the numerical instability observed by previous investigations and future algorithm development for this type of extreme flow

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

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

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

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures

Experimental Investigation of Nozzle/Plume Aerodynamics at Hypersonic Speeds

The work performed by D. W. Bogdanoff and J.-L. Cambier during the period of 1 Feb. - 31 Oct. 1992 is presented. The following topics are discussed: (1) improvement in the operation of the facility; (2) the wedge model; (3) calibration of the new test section; (4) combustor model; (5) hydrogen fuel system for combustor model; (6) three inch calibration/development tunnel; (7) shock tunnel unsteady flow; (8) pulse detonation wave engine; (9) DCAF flow simulation; (10) high temperature shock layer simulation; and (11) the one dimensional Godunov CFD code

Aeronautical engineering: A continuing bibliography with indexes (supplement 296)

This bibliography lists 592 reports, articles, and other documents introduced into the NASA scientific and technical information system in Oct. 1993. Subject coverage includes: design, construction and testing of aircraft and aircraft engines; aircraft components, equipment, and systems; ground support systems; and theoretical and applied aspects of aerodynamics and general fluid dynamics

Experimental and Computational Analysis of the Interaction of Carbon and Silicon Ablation Products in Expanding Hypersonic Flows

Thermal protection is required for vehicles entering planetary atmospheres to protect against the severe heating loads experienced. Characterization of candidate materials is often done utilizing plasma or arc-jet facilities, which provide steady-state testing of the thermal environments experienced during hypersonic flight, but do not correctly simulate hypersonic flowfields. Conversely, impulse facilities can reproduce flight velocities and enthalpies but have extremely short test times, prohibiting testing of thermal response. Modeling ablation and heating rates, particularly in the wake region, remains a significant challenge due to the complexity of the flowfield. To better understand this complex phenomenon and provide data to validate current computational models, experiments were conducted at the X2 expansion tunnel at the University of Queensland. Preheated strips of C-C and SiC-coated C-C were mounted in a two-dimensional compression wedge and tested in Earth entry ow. Calibrated spectral measurements were obtained in the near-stagnation and expansion regions targeting atomic Si, CN violet, C2 Swan, atomic N, atomic O, CO, and CO2 emissions for surface temperatures from approximately 1500 K to 2700 K. Emissions for C-C and SiC appeared similar in the near-stagnation region, increasing near the wall, while emissions for SiC-coated C-C displayed a distinct rise downstream of the shock, which suggests a higher concentration of ablative species. Comparisons were made to simulated results, which were conducted using the LAURA and HARA simulation codes following the process developed for this work. There was generally good agreement for CN emissions, which were most dominant, while the agreement was not as good for the other radiative phenomena investigated. It is believed that the underprediction of the ablation rate of the equilibrium-char model is a key factor

NASA SBIR abstracts of 1990 phase 1 projects

The research objectives of the 280 projects placed under contract in the National Aeronautics and Space Administration (NASA) 1990 Small Business Innovation Research (SBIR) Phase 1 program are described. The basic document consists of edited, non-proprietary abstracts of the winning proposals submitted by small businesses in response to NASA's 1990 SBIR Phase 1 Program Solicitation. The abstracts are presented under the 15 technical topics within which Phase 1 proposals were solicited. Each project was assigned a sequential identifying number from 001 to 280, in order of its appearance in the body of the report. The document also includes Appendixes to provide additional information about the SBIR program and permit cross-reference in the 1990 Phase 1 projects by company name, location by state, principal investigator, NASA field center responsible for management of each project, and NASA contract number

Monte Carlo sensitivity analyses of DSMC parameters for ionizing hypersonic flows

This work focuses on the development and sensitivity analyses of a direct simulation Monte Carlo (DSMC) code to understand the complex physical processes that occur during hypersonic entry into a rarefied atmosphere. Simulations are performed on 1-dimensional hypersonic shock scenarios that mimic the conditions of high altitude atmospheric entry to Earth and Saturn with the Computation of Hypersonic Ionizing Particles in Shocks (CHIPS) code. To model hypersonic entry problems accurately, the CHIPS code must resolve nonequilibrium flows and account for a number of complex gas dynamics processes at the molecular level. In this thesis, several high temperature models are added to the CHIPS code including charged particle models and electronic excitation. These models are refined using preliminary sensitivity analyses resulting in improved electronic excitation models and a new backward chemical reaction model. The CHIPS simulations completed in this work reproduce rarefied hypersonic shock tube experiments performed in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center. The CHIPS results are post-processed by the NEQAIR line-by-line radiative solver to compare directly to spectra measured experimentally in EAST. The DSMC techniques used to model hypersonic phenomena require numerous experimentally calibrated parameters. Many of these parameters are inferred from lower temperature experiments, resulting in an unknown amount of uncertainty in the simulated results at the extreme conditions of hypersonic flow. A global Monte Carlo sensitivity analysis is performed by simultaneously varying the CHIPS input parameter values to understand the sensitivity of experimentally measured quantities simulated by the CHIPS and NEQAIR codes. The sensitivity of several of these output quantities is used to rank the input parameters, identifying the most important parameters for the simulation of the hypersonic scenario. It was concluded that experimentally measured radiation intensity is most sensitive to the following key processes: N+e⁻⇌N⁺+e⁻+e⁻, NO+N⁺⇌N+NO⁺, N₂+N⇌N+N+N, N+O⇌NO⁺+e⁻, N+N⇌N₂⁺+e⁻, and Z [subscript elec] for N, O, and N₂⁺. In the future, this ranking can be used to identify which input parameters should be experimentally investigated, where model improvements could be beneficial, and aid in reducing the parameter space for DSMC calibrations to experimental data.Aerospace Engineerin

Bibliography of Lewis Research Center technical publications announced in 1985

This compilation of abstracts describes and indexes the technical reporting that resulted from the scientific and engineering work performed and managed by the Lewis Research Center in 1985. All the publications were announced in the 1985 issues of STAR (Scientific and Technical Aerospace Reports) and/or IAA (International Aerospace Abstracts). Included are research reports, journal articles, conference presentations, patents and patent applications, and theses