SEARCH · Engineering Papers

Results for “Conservation equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1995

The development of concentration gradients in a suspension of chemotactic bacteria

When a suspension of bacterial cells of the species Bacillus subtilis is placed in a chamber with its upper surface open to the atmosphere complex bioconvection patterns are observed. These arise because the cells: (1) are denser than water; and (2) usually swim upwards, so that the density of an initially uniform suspension becomes greater at the top than the bottom. When the vertical density gradient becomes large enough, an overturning instability occurs which ultimately evolves into the observed patterns. The reason that the cells swim upwards is that they are aerotactic, i.e., they swim up gradients of oxygen, and they consume oxygen. These properties are incorporated in conservation equations for the cell (N) and oxygen (C) concentrations, and these are solved in the pre-instability phase of development when N and C depend only on the vertical coordinate and time. Numerical results are obtained for both shallow- and deep-layer chambers, which are intrinsically different and require different mathematical and numerical treatments. It is found that, for both shallow and deep chambers, a thin boundary layer, densely packed with cells, forms near the surface. Beneath this layer the suspension becomes severely depleted of cells. Furthermore, in the deep chamber cases, a discontinuity in the cell concentration arises between this cell-depleted region and a cell-rich region further below, where no significant oxygen concentration gradients develop before the oxygen is fully consumed. The results obtained from the model are in good qualitative agreement with the experimental observations.

NASA Discipline Cell Biology↗

NASA NTRS2005

Multi-Spacecraft Observations of Interplanetary Shocks

Multi-spacecraft solar wind observations reveal that many interplanetary shocks deviate significantly from exact planarity on scale length of the magnetospheric cross section. A number of different IP shock observations with four spacecraft will be presented to demonstrate quantitatively the angular deviations between shock normals obtained from 4-spacecraft methods, using only the time and position information of shock observations but assuming a exactly planar geometry, and those obtained from a non-linear least squares fitting of the "Rankine-Hugoniot" conservation equations at each spacecraft. Moreover, the curvature of the shock fronts is strongly related to its driver, typically magnetic clouds. It will be demonstrated that small and slower moving magnetic clouds drive shocks with significantly more irregular surface geometries.

Szabo, Adam↗

NASA NTRS2002

Dynamics of Diffusion Flames in von Karman Swirling Flows Studied

Von Karman swirling flow is generated by the viscous pumping action of a solid disk spinning in a quiescent fluid media. When this spinning disk is ignited in an oxidizing environment, a flat diffusion flame is established adjacent to the disk, embedded in the boundary layer (see the preceding illustration). For this geometry, the conservation equations reduce to a system of ordinary differential equations, enabling researchers to carry out detailed theoretical models to study the effects of varying strain on the dynamics of diffusion flames. Experimentally, the spinning disk burner provides an ideal configuration to precisely control the strain rates over a wide range. Our original motivation at the NASA Glenn Research Center to study these flames arose from a need to understand the flammability characteristics of solid fuels in microgravity where slow, subbuoyant flows can exist, producing very small strain rates. In a recent work (ref. 1), we showed that the flammability boundaries are wider and the minimum oxygen index (below which flames cannot be sustained) is lower for the von Karman flow configuration in comparison to a stagnation-point flow. Adding a small forced convection to the swirling flow pushes the flame into regions of higher strain and, thereby, decreases the range of flammable strain rates. Experiments using downward facing, polymethylmethacrylate (PMMA) disks spinning in air revealed that, close to the extinction boundaries, the flat diffusion flame breaks up into rotating spiral flames (refs. 2 and 3). Remarkably, the dynamics of these spiral flame edges exhibit a number of similarities to spirals observed in biological systems, such as the electric pulses in cardiac muscles and the aggregation of slime-mold amoeba. The tail of the spiral rotates rigidly while the tip executes a compound, meandering motion sometimes observed in Belousov-Zhabotinskii reactions.

Nayagam, Vedha↗

NASA NTRS2005

Dynamic Deployment Simulations of Inflatable Space Structures

The feasibility of using Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method in LSDYNA to simulate the dynamic deployment of inflatable space structures is investigated. The CV and ALE methods were used to predict the inflation deployments of three folded tube configurations. The CV method was found to be a simple and computationally efficient method that may be adequate for modeling slow inflation deployment sine the inertia of the inflation gas can be neglected. The ALE method was found to be very computationally intensive since it involves the solving of three conservative equations of fluid as well as dealing with complex fluid structure interactions.

Wang, John T.↗

NASA NTRS2006

A Theoretical Solid Oxide Fuel Cell Model for System Controls and Stability Design

As the aviation industry moves towards higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The Hybrid Solid Oxide Fuel Cell system combines the fuel cell with a microturbine to obtain up to 70 percent cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multi-discipline system, and design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and distribution system and the fuel cell and microturbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. A novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled, but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗

NASA NTRS2008

A Theoretical Solid Oxide Fuel Cell Model for Systems Controls and Stability Design

As the aviation industry moves toward higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The hybrid solid oxide fuel cell system combines the fuel cell with a micro-turbine to obtain up to 70% cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multidiscipline system and the design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and a distribution system, and the fuel cell and micro-turbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. The novelty in this paper is that, first, the case is made why a high fidelity fuel cell mode is needed for systems control and stability designs. Second, a novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗

NASA NTRS2010

Analytic Modeling of Pressurization and Cryogenic Propellant

An analytic model for pressurization and cryogenic propellant conditions during all mission phases of any liquid rocket based vehicle has been developed and validated. The model assumes the propellant tanks to be divided into five nodes and also implements an empirical correlation for liquid stratification if desired. The five nodes include a tank wall node exposed to ullage gas, an ullage gas node, a saturated propellant vapor node at the liquid-vapor interface, a liquid node, and a tank wall node exposed to liquid. The conservation equations of mass and energy are then applied across all the node boundaries and, with the use of perfect gas assumptions, explicit solutions for ullage and liquid conditions are derived. All fluid properties are updated real time using NIST Refprop.1 Further, mass transfer at the liquid-vapor interface is included in the form of evaporation, bulk boiling of liquid propellant, and condensation given the appropriate conditions for each. Model validation has proven highly successful against previous analytic models and various Saturn era test data and reasonably successful against more recent LH2 tank self pressurization ground test data. Finally, this model has been applied to numerous design iterations for the Altair Lunar Lander, Ares V Core Stage, and Ares V Earth Departure Stage in order to characterize Helium and autogenous pressurant requirements, propellant lost to evaporation and thermodynamic venting to maintain propellant conditions, and non-uniform tank draining in configurations utilizing multiple LH2 or LO2 propellant tanks. In conclusion, this model provides an accurate and efficient means of analyzing multiple design configurations for any cryogenic propellant tank in launch, low-acceleration coast, or in-space maneuvering and supplies the user with pressurization requirements, unusable propellants from evaporation and liquid stratification, and general ullage gas, liquid, and tank wall conditions as functions of time.

Corpening, Jeremy H.↗

NASA NTRS2010

Turbofan Volume Dynamics Model for Investigations of Aero-Propulso-Servo-Elastic Effects in a Supersonic Commercial Transport

A turbofan simulation has been developed for use in aero-propulso-servo-elastic coupling studies, on supersonic vehicles. A one-dimensional lumped volume approach is used whereby each component (fan, high-pressure compressor, combustor, etc.) is represented as a single volume using characteristic performance maps and conservation equations for continuity, momentum and energy. The simulation is developed in the MATLAB/SIMULINK (The MathWorks, Inc.) environment in order to facilitate controls development, and ease of integration with a future aero-servo-elastic vehicle model being developed at NASA Langley. The complete simulation demonstrated steady state results that closely match a proposed engine suitable for a supersonic business jet at the cruise condition. Preliminary investigation of the transient simulation revealed expected trends for fuel flow disturbances as well as upstream pressure disturbances. A framework for system identification enables development of linear models for controller design. Utilizing this framework, a transfer function modeling an upstream pressure disturbance s impacts on the engine speed is developed as an illustrative case of the system identification. This work will eventually enable an overall vehicle aero-propulso-servo-elastic model

Connolly, Joseph W.↗

NASA NTRS2011

A Simple, Analytical Model of Collisionless Magnetic Reconnection in a Pair Plasma

A set of conservation equations is utilized to derive balance equations in the reconnection diffusion region of a symmetric pair plasma. The reconnection electric field is assumed to have the function to maintain the current density in the diffusion region, and to impart thermal energy to the plasma by means of quasi-viscous dissipation. Using these assumptions it is possible to derive a simple set of equations for diffusion region parameters in dependence on inflow conditions and on plasma compressibility. These equations are solved by means of a simple, iterative, procedure. The solutions show expected features such as dominance of enthalpy flux in the reconnection outflow, as well as combination of adiabatic and quasi-viscous heating. Furthermore, the model predicts a maximum reconnection electric field of E(sup *)=0.4, normalized to the parameters at the inflow edge of the diffusion region.

Hesse, Michael↗

NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

NASA NTRS2003

Noncircular Cross Sections Could Enhance Mixing in Sprays

A computational study has shown that by injecting drops in jets of gas having square, elliptical, triangular, or other noncircular injection cross sections, it should be possible to increase (relative to comparable situations having circular cross section) the entrainment and dispersion of liquid drops. This finding has practical significance for a variety of applications in which it is desirable to increase dispersion of drops. For example, in chemical-process sprays, increased dispersion leads to increases in chemical- reaction rates; in diesel engines, increasing the dispersion of drops of sprayed fuel reduces the production of soot; and in household and paint sprays, increasing the dispersion of drops makes it possible to cover larger surfaces. It has been known for some years that single-phase fluid jets that enter flow fields through noncircular inlets entrain more fluid than do comparable jets entering through circular inlets. The computational study reported here was directed in part toward determining whether and how this superior mixing characteristic of noncircular single phase jets translates to a similar benefit in cases of two-phase jets (that is, sprays). The study involved direct numerical simulations of single- and two-phase free jets with circular, elliptical, rectangular, square, and triangular inlet cross sections. The two-phase jets consisted of gas laden with liquid drops randomly injected at the inlets. To address the more interesting case of evaporating drops, the carrier gas in the jets was specified to be initially unvitiated by the vapor of the liquid chemical species and the initial temperature of the drops was chosen to be smaller than that of the gas. The mathematical model used in the study was constructed from the conservation equations for the two-phase flow and included complete couplings of mass, momentum, and energy based on thermodynamically self-consistent specification of the enthalpy, internal energy, and latent heat of vaporization of the vapor.

Bellan, Josette↗

NASA NTRS2011

Large Eddy Simulation Study for Fluid Disintegration and Mixing

A new modeling approach is based on the concept of large eddy simulation (LES) within which the large scales are computed and the small scales are modeled. The new approach is expected to retain the fidelity of the physics while also being computationally efficient. Typically, only models for the small-scale fluxes of momentum, species, and enthalpy are used to reintroduce in the simulation the physics lost because the computation only resolves the large scales. These models are called subgrid (SGS) models because they operate at a scale smaller than the LES grid. In a previous study of thermodynamically supercritical fluid disintegration and mixing, additional small-scale terms, one in the momentum and one in the energy conservation equations, were identified as requiring modeling. These additional terms were due to the tight coupling between dynamics and real-gas thermodynamics. It was inferred that if these terms would not be modeled, the high density-gradient magnitude regions, experimentally identified as a characteristic feature of these flows, would not be accurately predicted without the additional term in the momentum equation; these high density-gradient magnitude regions were experimentally shown to redistribute turbulence in the flow. And it was also inferred that without the additional term in the energy equation, the heat flux magnitude could not be accurately predicted; the heat flux to the wall of combustion devices is a crucial quantity that determined necessary wall material properties. The present work involves situations where only the term in the momentum equation is important. Without this additional term in the momentum equation, neither the SGS-flux constant-coefficient Smagorinsky model nor the SGS-flux constant-coefficient Gradient model could reproduce in LES the pressure field or the high density-gradient magnitude regions; the SGS-flux constant- coefficient Scale-Similarity model was the most successful in this endeavor although not totally satisfactory. With a model for the additional term in the momentum equation, the predictions of the constant-coefficient Smagorinsky and constant-coefficient Scale-Similarity models were improved to a certain extent; however, most of the improvement was obtained for the Gradient model. The previously derived model and a newly developed model for the additional term in the momentum equation were both tested, with the new model proving even more successful than the previous model at reproducing the high density-gradient magnitude regions. Several dynamic SGS-flux models, in which the SGS-flux model coefficient is computed as part of the simulation, were tested in conjunction with the new model for this additional term in the momentum equation. The most successful dynamic model was a "mixed" model combining the Smagorinsky and Gradient models. This work is directly applicable to simulations of gas turbine engines (aeronautics) and rocket engines (astronautics).

Bellan, Josette↗

NASA NTRS2012

3D Drop Size Distribution Extrapolation Algorithm Using a Single Disdrometer

Determining the Z-R relationship (where Z is the radar reflectivity factor and R is rainfall rate) from disdrometer data has been and is a common goal of cloud physicists and radar meteorology researchers. The usefulness of this quantity has traditionally been limited since radar represents a volume measurement, while a disdrometer corresponds to a point measurement. To solve that problem, a 3D-DSD (drop-size distribution) method of determining an equivalent 3D Z-R was developed at the University of Central Florida and tested at the Kennedy Space Center, FL. Unfortunately, that method required a minimum of three disdrometers clustered together within a microscale network (.1-km separation). Since most commercial disdrometers used by the radar meteorology/cloud physics community are high-cost instruments, three disdrometers located within a microscale area is generally not a practical strategy due to the limitations of these kinds of research budgets. A relatively simple modification to the 3D-DSD algorithm provides an estimate of the 3D-DSD and therefore, a 3D Z-R measurement using a single disdrometer. The basis of the horizontal extrapolation is mass conservation of a drop size increment, employing the mass conservation equation. For vertical extrapolation, convolution of a drop size increment using raindrop terminal velocity is used. Together, these two independent extrapolation techniques provide a complete 3DDSD estimate in a volume around and above a single disdrometer. The estimation error is lowest along a vertical plane intersecting the disdrometer position in the direction of wind advection. This work demonstrates that multiple sensors are not required for successful implementation of the 3D interpolation/extrapolation algorithm. This is a great benefit since it is seldom that multiple sensors in the required spatial arrangement are available for this type of analysis. The original software (developed at the University of Central Florida, 1998.- 2000) has also been modified to read standardized disdrometer data format (Joss-Waldvogel format). Other modifications to the software involve accounting for vertical ambient wind motion, as well as evaporation of the raindrop during its flight time.

Lane, John↗

NASA NTRS2012

Ice Accretion Modeling using an Eulerian Approach for Droplet Impingement

A three-dimensional Eulerian analysis has been developed for modeling droplet impingement on lifting bodes. The Eulerian model solves the conservation equations of mass and momentum to obtain the droplet flow field properties on the same mesh used in CFD simulations. For complex configurations such as a full rotorcraft, the Eulerian approach is more efficient because the Lagrangian approach would require a significant amount of seeding for accurate estimates of collection efficiency. Simulations are done for various benchmark cases such as NACA0012 airfoil, MS317 airfoil and oscillating SC2110 airfoil to illustrate its use. The present results are compared with results from the Lagrangian approach used in an industry standard analysis called LEWICE.

Kim, Joe Woong↗

NASA NTRS2013

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

NASA NTRS2013

Global Propagation of Gravity Waves Generated with the Whole Atmosphere Transfer Function Model

A brief review is presented of the Transfer Function Model (TFM) [e.g., Mayr et al., Space Science Reviews, 1990], which describes acoustic gravity waves (AGW) that propagate across the globe in a dissipative and static (no winds) background atmosphere with globally uniform temperature and density variations extending from the ground to 700 km. Unique among existing models, the TFM can be placed between the analytical approach on one end, and the rigorous numerical approach of general circulation models (GCM). The time consuming numerical integration of the conservation equations is restricted to compute the transfer function (TF) for a broad range of frequencies and spherical harmonics. Given TF, the atmospheric response for a chosen source configuration is then obtained in short order. Computationally efficient, the model is well suited to serve as experimental and educational tool for simulating propagating wave patterns across the globe. By design, the TFM is also semi-analytical and therefore well suited to explore the different wave modes that can be generated under different dynamical conditions.

Transfer Function Model↗

NASA NTRS2017

Modeling of Liquefaction of Cryogenic Propellant in a Tank

Over the past decades NASA has been focusing to develop technology that would to allow for production of cryogenic propellants on the surface of Mars. The in-situ propellant production reduces the amount of propellants needed to be taken to Mars and ultimately to reduce mission cost. Utilizing Martian resources, the produced gaseous propellants (i.e., oxygen and methane) are liquefied and stored prior to use on the Mars ascent vehicle. In this paper, a model for the liquefaction process of gaseous propellants in a cryogenically refrigerated tank is presented. The tank is considered to be cylindrical with elliptical top and bottom domes. A multi-node transient model is developed based on the mass and energy conservation equations and wall-gas and liquid-gas interface mass and heat transfer correlations. Description of the model and predicted results will be presented in the final paper.

Hedayat, A.↗

NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗