Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-equilibrium properties”

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 37 records · Page 2

A time-accurate algorithm for chemical non-equilibrium viscous flows at all speeds

A time-accurate, coupled solution procedure is described for the chemical nonequilibrium Navier-Stokes equations over a wide range of Mach numbers. This method employs the strong conservation form of the governing equations, but uses primitive variables as unknowns. Real gas properties and equilibrium chemistry are considered. Numerical tests include steady convergent-divergent nozzle flows with air dissociation/recombination chemistry, dump combustor flows with n-pentane-air chemistry, nonreacting flow in a model double annular combustor, and nonreacting unsteady driven cavity flows. Numerical results for both the steady and unsteady flows demonstrate the efficiency and robustness of the present algorithm for Mach numbers ranging from the incompressible limit to supersonic speeds.

Shuen, J.-S.↗

Effect of non-equilibrium flow chemistry and surface catalysis on surface heating to AFE

The effect of nonequilibrium flow chemistry on the surface temperature distribution over the forebody heat shield on the Aeroassisted Flight Experiment (AFE) vehicle was investigated using a reacting boundary-layer code. Computations were performed by using boundary-layer-edge properties determined from global iterations between the boundary-layer code and flow field solutions from a viscous shock layer (VSL) and a full Navier-Stokes solution. Surface temperature distribution over the AFE heat shield was calculated for two flight conditions during a nominal AFE trajectory. This study indicates that the surface temperature distribution is sensitive to the nonequilibrium chemistry in the shock layer. Heating distributions over the AFE forebody calculated using nonequilibrium edge properties were similar to values calculated using the VSL program.

Stewart, David A.↗

On the Boundary Condition for Water at a Hydrophobic, Dense Surface

We study the no-slip boundary conditions for water at a hydrophobic (graphite) surface using non-equilibrium molecular-dynamics simulations. For the planar Couette flow, we find a slip length of 64 nm at 1 bar and 300 K, decreasing with increasing system pressure to a value of 31 nm at 1000 bar. Changing the properties of the interface to from hydrophobic to strongly hydrophilic reduces the slip to 14 nm. Finally, we study the flow of water past an array of carbon nanotubes mounted in an inline configuration with a spacing of 16.4 x 16.4 nm. For tube diameters of 1.25 and 2.50 nm we find drag coefficients in good agreement with the macroscopic, Navier-Stokes values. For carbon nanotubes, the no-slip condition is valid to within the definition of the position of the interface.

Walther, J. H.↗

Boundary Layer Transition Correlations and Aeroheating Predictions for Mars Smart Lander

Laminar and turbulent perfect-gas air, Navier-Stokes computations have been performed for a proposed Mars Smart Lander entry vehicle at Mach 6 over a free stream Reynolds number range of 6.9 x 10(exp 6)/m to 2.4 x 10(exp 7)/m (2.1 x 10(exp 6)/ft to 7.3 x 10(exp 6)/ft) for angles-of-attack of 0-deg, 11-deg, 16-deg, and 20-deg, and comparisons were made to wind tunnel heating data obtained a t the same conditions. Boundary layer edge properties were extracted from the solutions and used to correlate experimental data on the effects of heat-shield penetrations (bolt-holes where the entry vehicle would be attached to the propulsion module during transit to Mars) on boundary-layer transition. A non-equilibrium Martian-atmosphere computation was performed for the peak heating point on the entry trajectory in order to determine if the penetrations would produce boundary-layer transition by using this correlation.

Hollis, Brian R.↗

Boundary Layer Transition Correlations and Aeroheating Predictions for Mars Smart Lander

Laminar and turbulent perfect-gas air, Navier-Stokes computations have been performed for a proposed Mars Smart Lander entry vehicle at Mach 6 over a free stream Reynolds number range of 6.9 x 10(exp 6/m to 2.4 x 10(exp 7)m(2.1 x 10(exp 6)/ft to 7.3 x 10(exp 6)ft) for angles-of-attack of 0-deg, 11-deg, 16-deg, and 20-deg, and comparisons were made to wind tunnel heating data obtained at the same conditions. Boundary layer edge properties were extracted from the solutions and used to correlate experimental data on the effects of heat-shield penetrations (bolt-holes where the entry vehicle would be attached to the propulsion module during transit to Mars) on boundary-layer transition. A non-equilibrium Martian-atmosphere computation was performed for the peak heating point on the entry trajectory in order to determine if the penetrations would produce boundary-layer transition by using this correlation.

Hollis, Brian R.↗

Constraints on Solar Wind Plasma Properties Derived from Coordinated Coronal Observations

The goal of the proposed research is to increase the understanding of coronal plasma phenomena by making use of different observational approaches and combine the observations with the necessary theoretical considerations. We continued to study the formation of spectral lines in the corona/transition region under different non-equilibrium conditions. In addition to Mg and Ne we have also studied some cases involving Si ions and spectral lines. Due to the fact that the sun was at the maximum phase of the solar cycle, we spent some time on observing coronal mass ejections. Observations of the H I Lyman-alpha spectral line and the line pair 0 VI 1031.91 AA and 1037.61 AA were carried out with the UVCS instrument in the northern polar region of the sun at position angle 270 deg. The region was monitored at about 2 RS for about 5 hours on March 04, and for about 8 hours on March 05. During that time interval a major Coronal Mass Ejection developed in the northern hemisphere. Density, velocity and temperature maps of the ejected plasma have been obtained from the UVCS data. This event was also seen in the white light Large Angle and Spectrometric Coronagraph Experiment (LASCO) images, and its evolution at lower heights can be followed in Extreme Ultraviolet Imaging Telescope (EIT). The LASCO images are essential in providing the larger scale context for this event which is unique in the sense that it developed almost due North and had very little interactions with adjacent regions. The combination of UVCS velocity maps and LASCO images which were reduced using advanced image processing techniques, show very clearly how the mass ejection evolved from the solar surface to several solar radii, the twisting of the flux ropes, which are seen in UVCS as blue and red shifted velocities. First results were presented at the AGU Meeting in Boston. To study the quieter side of the coronal plasma, we carried out an experiment during the past eclipse, in June 01. We measured the corona in several iron spectral lines, and polarized white light. The results from that experiment look promising so far.

Esser, Ruth↗

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗

Metastable phase formation in undercooled liquid lead alloys

During non-equilibrium processes metastable phases are formed instead of stable phases due to the operation of various kinetic or structural constraints. By removing the most effective nucleation sites for the stable phase using emulsified droplet samples, stable phase formation can be prohibited in a broad range of undercooling and the phase space available in the metastable regime can be mapped out. With this method the thermodynamic properties of the undercooled liquid and the metastable phase boundaries corresponding to reversible metastable (solid + liquid) phase equilibria are examined from experimental data. The analysis reveals important implications for the nucleation kinetics and the choice of the kinetically most favored solidification path.

Fecht, Hans J.↗

High speed turbulent reacting flows: DNS and LES

Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.

Givi, Peyman↗

Analysis of Flow From Arc-Jet Spectra

Materials testing is carried out at the JSC Atmospheric Reentry Materials and Structures Facility. A flow diagnostics program is under development to characterize the energy distribution in arc-jet flows used to simulate atmospheric entry. Heat transfer to materials depends on flow properties which includes the composition of and energy distribution among the atoms, ions, molecules, and molecular ions which make up the flow. This project involves analysis of shock layer characteristics from the radiation emitted and experimentally collected from near the front of the shock to near the surface of the material. The analysis has yielded information on relative populations of neutral molecules and molecular ions within the layer. In determining non-equilibrium temperatures within the layer, some insight into the spectral constants used to compute radiative emission has been gained.

Blackwell, H. E.↗

Molecular Dynamics Simulations of the Thermal Conductivity of Single-Wall Carbon Nanotubes

Carbon nanotubes (CNT) have very attractive electronic, mechanical. and thermal properties. Recently, measurements of thermal conductivity in single wall CNT mats showed estimated thermal conductivity magnitudes ranging from 17.5 to 58 W/cm-K at room temperature. which are better than bulk graphite. The cylinderical symmetry of CNT leads to large thermal conductivity along the tube axis, additionally, unlike graphite. CNTs can be made into ropes that can be used as heat conducting pipes for nanoscale applications. The thermal conductivity of several single wall carbon nanotubes has been calculated over temperature range from l00 K to 600 K using non-equilibrium molecular dynamics using Tersoff-Brenner potential for C-C interactions. Thermal conductivity of single wall CNTs shows a peaking behavior as a function of temperature. Dependence of the peak position on the chirality and radius of the tube will be discussed and explained in this presentation.

Osman, M.↗

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗

What if chondritic porous interplanetary dust particles are not the real McCoy

To select a target comet for a Comet Nucleus Sample Return Mission (CNSRM) it is necessary to have an experimental data base to evaluate the extent of diversity and similarity of comets. For example, the physical properties (e.g., low density) of chondritic porous (CP) interplanetary dust particles (IDPs) are believed to resemble these properties of cometary dust although it is yet to be demonstrated that the porous structure of CP IDPs is inherent to presolar dust particles stored in comet nuclei. Porous structures of IDPs could conceivably form during sublimation at the surface of active comet nuclei. Porous structures are also obtained during annealing of amorphous Mg-SiO smokes which initially forms porous aggregates of olivine + platey tridymite and which, upon continued annealing, react to fluffy enstatite aggregates. It is therefore uncertain that CP IDPs are entirely composed of unmetamorphosed presolar dust. Conceivably, new minerals and textures may form in situ in nuclei of active comets as a function of their individual thermal history. Unmetamorphosed comet dust is probably structurally amorphous. Thermal annealing of this dust can produce ultra fine-grained minerals and this ultrafine grain size of CP IDPs should be considered in assessments of aqueous alterations that could affect presolar dust in comet nuclei between 200 and 400 K. Devitrification and hydration may occur in situ in ice-dust mixtures and the mantle of active comet nuclei. Devitrification, or uncontrolled crystallization, of amorphous precursor dust can produce a range of chemical compositions of ultrafine-grained minerals and (non-equilibrium) mineral assemblages and textures in dust contained in comet nuclei as a function of period and trajectory of orbit and number of perihelion passages (not considering internal heating). Thus, experimental data on relevant processes and reaction rates between 200 and 400 K are needed in order to evaluate comet selection, penetration depth for sampling device and curation of samples for CNSRM.

Rietmeijer, Frans J. M.↗

Metamorphism of Cosmic Dust: Processing from Circumstellar Outflows to the Cometary Regolith

Nucleation is a non-equilibrium process: the products of this process are seldom the most thermodynamically stable condensates but are instead those which form fastest. It should therefore not be surprising that grains formed in a circumstellar outflow will undergo some degree of metamorphism if they are annealed or are exposed to a chemically active reagent. Metamorphism of refractory particles continues in the interstellar medium (ISM) where the driving forces are sputtering by cosmic ray particles, annealing by high energy photons and grain destruction in supernova generated shocks. Studies of the depletion of the elements from the gas phase of the interstellar medium tell us that if grain destruction occurs with high efficiency in the ISM, then there must be some mechanism by which grains can be formed in the ISM. Various workers have shown that refractory mantles could form on refractory cores by radiation processing of organic ices. A similar process may operate to produce refractory inorganic mantles on grain cores which survived the supernova shocks. Most grains in a cloud which collapses to form a star will be destroyed; many of the surviving grains will be severely processed. Grains in the outermost regions of the nebula may survive relatively unchanged by thermal processing or hydration. It is these grains which we hope to find in comets. However, only those grains encased in ice at low temperature can be considered pristine since a considerable degree of hydrous alteration might occur in a cometary regolith if the comet enters the inner solar system. Some discussion of the physical, chemical and isotopic properties of a refractory grain at each stage of its life cycle will be attempted based on the limited laboratory data available to date. Suggestions will be made concerning the types of experimental data which are needed in order to better understand the processing history of cosmic dust.

Nuth, Joseph A., III↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Upwind schemes and bifurcating solutions in real gas computations

The area of high speed flow is seeing a renewed interest due to advanced propulsion concepts such as the National Aerospace Plane (NASP), Space Shuttle, and future civil transport concepts. Upwind schemes to solve such flows have become increasingly popular in the last decade due to their excellent shock capturing properties. In the first part of this paper the authors present the extension of the Osher scheme to equilibrium and non-equilibrium gases. For simplicity, the source terms are treated explicitly. Computations based on the above scheme are presented to demonstrate the feasibility, accuracy and efficiency of the proposed scheme. One of the test problems is a Chapman-Jouguet detonation problem for which numerical solutions have been known to bifurcate into spurious weak detonation solutions on coarse grids. Results indicate that the numerical solution obtained depends both on the upwinding scheme used and the limiter employed to obtain second order accuracy. For example, the Osher scheme gives the correct CJ solution when the super-bee limiter is used, but gives the spurious solution when the Van Leer limiter is used. With the Roe scheme the spurious solution is obtained for all limiters.

Suresh, Ambady↗

GRO: Black hole models for gamma-ray bursts

The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) has established that the distribution of gamma-ray bursts (GRB's) is isotropic but is bound radially. This finding suggests that the bursts are either cosmological or they originate from an extended Galactic halo. The implied luminosities and the observed variability of the GRB's on time scales as short as one millisecond suggest that they originate from compact objects. We are presently studying black hole models for GRB's. Any such model must produce a non-thermal photon spectrum to agree with the observed properties. For a wide range of burst parameters the assumed bursting source consists of a non-thermal electron-positron-photon plasma of very high density. It seems possible to produce such a plasma in accretion onto black holes. In our on-going work, we are developing the kinetic theory for a non-equilibrium pair plasma. The main new features of our work are as follows: (1) We do not assume the presence of a thermal electron bath. (2) Non-thermal, high-energy pairs are allowed to have an arbitrary concentration and energy distribution. (3) There is no soft photon source in our model; initially all the photons in the plasma are either energetic X-rays or gamma-rays. (4) The initial energy distribution of the pairs as well as photons is arbitrary. (5) We collect the analytical expressions for the kinetic kernels for all relevant processes. And (6) we present a different approach to finding the time-evolution of pair and photon spectra, which is a combination of the kinetic-theory and the non-linear Monte-Carlo schemes. We have developed many Monte-Carlo programs to model various process, to take into account the time evolution, and to incorporate various physical effects which are unique to non-thermal plasmas. The hydrodynamics of fireballs in GRB's was studied before. Applying results from kinetic theory will improve our understanding of these systems.

Ruderman, Malvin↗

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗