Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibrium programming”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Measurements of the H2(13)CO ortho/para ratio in cold dark molecular clouds

H2(13)CO has been detected for the first time toward cold dark molecular clouds using the NRAO 12 m telescope. The H2(13)CO ortho/para abundance ratio R for B335, which we report as R approximately 1.7, suggests equilibrium at the local kinetic temperature and appears to be distinctly different from that for both TMC-1 and L134N, where R is close to or higher than the statistical value 3. Since only B335 among the observed positions includes an imbedded IR source, this difference may result from heating of the grain surfaces, providing the energy necessary for desorption of formaldehyde formed on the grains.

NASA Discipline Exobiology↗
NASA NTRSDate not supplied

2020 International Astronautical Congress Conference Presentation - ANALYTICAL TECHNIQUES FOR ASSESSING GATEWAY AND OTHER SPACECRAFT ANTENNA LINE-OF-SIGHT FOR THE ARTEMIS PROGRAM

NASA’s Artemis program is committed to landing the next man and the first woman on the Moon by 2024. The Gateway, a critical piece of infrastructure for the long-term Artemis mission profiles, will inaugurate a sustainable crewed presence beyond low-Earth orbit in cis-lunar space and serve as a staging point for a lunar landing system. The Near Rectilinear Halo Orbit of the Gateway offers numerous operational benefits in terms of its Earth access via the Orion Multipurpose Crew Vehicle, its degree of Earth and lunar visibility, the versatile surface access it provides, and favorable conditions for station-keeping. To minimize burns during station-keeping and conserve propellant, the Gateway will continually align itself to a solar pressure equilibrium attitude with the Sun; this is convenient for solar array orientation and power generation, but causes difficulties in terms of communication line-of-sight with the Earth, Moon, and visiting vehicles, which appear to be continually rotating about the Gateway’s local frame. A fixed antenna on the Gateway is continually subject to a moving target and can find itself obstructed by the Gateway’s own geometry frequently if not placed carefully. To understand the performance of either Earth-pointing or Moon-pointing antenna placements for a given Gateway geometry, a transient solar system simulation of the Gateway is used to step through orbital ephemeris data and determine the periods in which line-of-sight to a chosen target is achieved and lost. Antenna placements can be tested in the simulation and optimized to increase the average line-of-sight exposure to the target, while also minimizing the maximum duration drop-out in line-of-sight communication experienced; long durations without a successful communication link could pose threats to operations, safety, and mission success. It may be necessary to deviate from the Gateway’s solar pressure equilibrium attitude for short periods of time in order to continue a successful line-of-sight link to a chosen target. Communication coverage spheres can offer insight into what these necessary attitude adjustments may be and are generated by using a 3D mesh of the Gateway and an implementation of the Möller-Trumbore intersection algorithm. This paper outlines the analytical techniques used to perform such antenna placements on the Gateway, offering examples of how line-of-sight strength can be enhanced with an understanding of the Gateway geometry and attitude constraints. These techniques have had direct impact on the writing of communication requirements for NASA and international partner-owned Gateway elements.

Gateway↗
NASA NTRS1980

Thermal elastohydrodynamic lubrication of spur gears

An analysis and computer program called TELSGE were developed to predict the variations of dynamic load, surface temperature, and lubricant film thickness along the contacting path during the engagement of a pair of involute spur gears. The analysis of dynamic load includes the effect of gear inertia, the effect of load sharing of adjacent teeth, and the effect of variable tooth stiffness which are obtained by a finite-element method. Results obtained from TELSGE for the dynamic load distributions along the contacting path for various speeds of a pair of test gears show patterns similar to that observed experimentally. Effects of damping ratio, contact ratio, tip relief, and tooth error on the dynamic load were examined. In addition, two dimensionless charts are included for predicting the maximum equilibrium surface temperature, which can be used to estimate directly the lubricant film thickness based on well established EHD analysis.

Wang, K. L.↗
DOE OSTI2018

Non-inductively driven tokamak plasmas at near-unity β t in the Pegasus toroidal experiment

A major goal of the spherical tokamak (ST) research program is accessing a state of low internal inductance ℓ i , high elongation κ , and high toroidal and normalized beta ( β t and β N ) without solenoidal current drive. Local helicity injection (LHI) in the Pegasus ST [Garstka et al., Nucl. Fusion 46, S603 (2006)] provides non-solenoidally driven plasmas that exhibit these characteristics. LHI utilizes compact, edge-localized current sources for plasma startup and sustainment. It results in hollow current density profiles with low ℓ i . The low aspect ratio ( R 0 / a ~ 1.2 ) of Pegasus allows access to high κ and high normalized plasma currents I N = I p / a B T > 14 ). Magnetic reconnection during LHI provides auxiliary ion heating. Together, these features provide access to very high β t plasmas. Equilibrium analyses indicate that β t up to ~100% is achieved. Finally, these high β t discharges disrupt at the ideal no-wall β limit at β N ~ 7.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Simultaneous concentration and purification through gradient deformation chromatography

Mobile-phase additives, commonly used to modulate absorbate retention in gradient elution chromatography, are usually assumed to be either linearly retained or unretained. Previous theoretical work from our laboratory has shown that these modulators, such as salts in ion-exchange and hydrophobic interaction chromatography and organic modifiers in reversed-phase chromatography, can absorb nonlinearly, giving rise to gradient deformation. Consequently, adsorbate peaks that elute in the vicinity of the head of the deformed gradient may exhibit unusual shapes, form shoulders, and/or be concentrated. These effects for a reversed-phase sorbent with aqueous acetonitrile (ACN) as the modulator are verified experimentally. Gradient deformation is demonstrated experimentally and agrees with simulations based on ACN isotherm parameters that are independently determined from batch equilibrium studies using the layer model. Unusual absorbate peak shapes were found experimentally for single-component injections of phenylalanine, similar to those calculated by the simulations. A binary mixture of tryptophan and phenylalanine is used to demonstrate simultaneous concentration and separation, again in agreement with simulations. The possibility of gradient deformation in ion-exchange and hydrophobic interaction chromatography is discussed.

NASA Discipline Number 93-10↗
NASA NTRS2019

An Analytic Model for Estimating the First Contact Resistance Needed to Avoid Damaging ESD During Spacecraft Docking in GEO

NASA’s Gateway program will involve spacecraft (s/c) docking in the outer radiation belt in order to transfer Gateway elements between s/c for transport to lunar orbit. The charging of these s/c to different potentials prior to docking raises the possibility of a damaging electrostatic discharge (ESD) at the time of first contact between the s/c. A proposed mitigation strategy is for first contact to occur prior to docking through a resistor with resistance R that would lower the potential difference at an optimal rate to a sufficiently low value to prevent ESD damage. The coupling of s/c by a resistor can be modeled by SPIS, but for realistic two s/c models SPIS can take hours to simulate the evolution of the s/c surface charges and potentials to an equilibrium state. Our objective is to develop a simpler model of s/c resistive coupling that runs orders of magnitude faster while providing useful first design estimates of the time variation of the s/c potentials, current through the resistor, and how these vary with R and s/c configuration This configuration is defined by the relative separation and orientation of the s/c, and their solar illumination. The configuration and geometry of the s/c determine their capacitive coupling. The s/c capacitances are computed using Nascap-2K. This poster presents the first version of such a model, and initial tests.

Goodman, Michael L.↗
DOE OSTI2018

Validation and Initial Characterization of the Long-period Planet Kepler-1654 b

Fewer than 20 transiting Kepler planets have periods longer than one year. Our early search of the Kepler light curves revealed one such system, Kepler-1654b (originally KIC 8410697b), which shows exactly two transit events and whose second transit occurred only five days before the failure of the second of two reaction wheels brought the primary Kepler mission to an end. A number of authors have also examined light curves from the Kepler mission searching for long-period planets and identified this candidate. Starting in 2014 September, we began an observational program of imaging, reconnaissance spectroscopy, and precision radial velocity (RV) measurements that confirm with a high degree of confidence that Kepler-1654b is a bona fide transiting planet orbiting a mature G5V star (T {sub eff} = 5580 K, [Fe/H] = −0.08) with a semimajor axis of 2.03 au, a period of 1047.84 days, and a radius of 0.82 ± 0.02 R {sub Jup}. RV measurements using Keck’s HIRES spectrometer obtained over 2.5 years set a limit to the planet’s mass of <0.5 (3σ) M {sub Jup}. The bulk density of the planet is similar to that of Saturn or possibly lower. We assess the suitability of temperate gas giants like Kepler-1654b for transit spectroscopy with the James Webb Space Telescope, as their relatively cold equilibrium temperatures (T {sub pl} ∼ 200 K) make them interesting from the standpoint of exoplanet atmospheric physics. Unfortunately, these low temperatures also make the atmospheric scale heights small and thus transmission spectroscopy challenging. Finally, the long time between transits can make scheduling JWST observations difficult—as is the case with Kepler-1654b.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Sea-level response to ice sheet evolution: An ocean perspective

The ocean's influence upon and response to Antarctic ice sheet changes is considered in relation to sea level rise over recent and future decades. Assuming present day ice fronts are in approximate equilibrium, a preliminary budget for the ice sheet is estimated from accumulation vs. iceberg calving and the basal melting that occurs beneath floating ice shelves. Iceberg calving is derived from the volume of large bergs identified and tracked by the Navy/NOAA Joint Ice Center and from shipboard observations. Basal melting exceeds 600 cu km/yr and is concentrated near the ice fronts and ice shelf grounding lines. An apparent negative mass balance for the Antarctic ice sheet may result from an anomalous calving rate during the past decade, but there are large uncertainties associated with all components of the ice budget. The results from general circulation models are noted in the context of projected precipitation increases and ocean temperature changes on and near the continent. An ocean research program that could help refine budget estimates is consistent with goals of the West Antarctic Ice Sheet Initiative.

Jacobs, Stanley S.↗
DOE OSTI2019

Orbital angular momentum constraints in the variational optimization of the two-electron reduced-density matrix

The direct variational determination of the two-electron reduced-density matrix (2-RDM) corresponding to an atomic or molecular system is usually carried out in a basis of real-valued atom-centered Gaussian basis functions, under the assumption that the 2-RDM is a real-valued quantity. Furthermore, for systems that possess orbital angular momentum symmetry, the description of states with a well-defined, nonzero z projection of the orbital angular momentum requires a computational framework generalized to include either complex basis functions or a complex-valued 2-RDM. We consider a semidefinite program suitable for the direct optimization of a complex-valued 2-RDM and explore the role of orbital angular momentum constraints in systems that possess the relevant symmetries. For atomic systems, constraints on the expectation values of the square and z projection of the orbital angular momentum operator allow one to optimize 2-RDMs for multiple orbital angular momentum states. Similarly, in linear molecules, orbital angular momentum projection constraints enable the description of multiple electronic states and, moreover, the application of such constraints is essential for a qualitatively correct description of the electronic structure. For example, in the case of molecular oxygen, we demonstrate that orbital angular momentum constraints are necessary to recover the correct energy ordering of the lowest-energy singlet and triplet states near the equilibrium geometry. However, care must still be taken in the description of the dissociation limit, because the 2-RDM-based approach is not size consistent, and the size-consistency error varies dramatically, depending on the z projections of the spin and orbital angular momenta.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1984

Production of large-particle-size monodisperse latexes

The research program achieved two objectives: (1) it has refined and extended the experimental techniques for preparing monodisperse latexes in quantity on the ground up to a particle diameter of 10 microns; and (2) it has demonstrated that a microgravity environment can be used to grow monodisperse latexes to larger sizes, where the limitations in size have yet to be defined. The experimental development of the monodisperse latex reactor (MLR) and the seeded emulsion polymerizations carried out in the laboratory prototype of the flight hardware, as a function of the operational parameters is discussed. The emphasis is directed towards the measurement, interpretation, and modeling of the kinetics of seeded emulsion polymerization and successive seeded emulsion polymerization. The recipe development of seeded emulsion polymerization as a function of particle size is discussed. The equilibrium swelling of latex particles with monomers was investigated both theoretically and experimentally. Extensive studies are reported on both the type and concentration of initiators, surfactants, and inhibitors, which eventually led to the development of the flight recipes. The experimental results of the flight experiments are discussed, as well as the experimental development of inhibition of seeded emulsion polymerization in terms of time of inhibition and the effect of inhibitors on the kinetics of polymerization.

Vanderhoff, J. W.↗
NASA NTRS1991

The density of the cell sap and endoplasm of Nitellopsis and Chara

We measured the densities of the cell sap, endoplasm and cell wall of Nitellopsis obtusa and Chara corallina using interference microscopy, refractometry, immersion refractometry, equilibrium sedimentation and chemical microanalysis techniques. These values are important for the determination of many rheological properties of the cytoplasm as well as for understanding buoyancy regulation, dispersal mechanisms and how cells respond to gravity. The average densities of the cell sap, endoplasm and cell wall are 1,006.9, 1,016.7 and 1,371 kg m-3 for Nitellopsis and 1,005.0, 1,013.9, and 1,355.3 kg m-3 for Chara.

NASA Discipline Number 40-50↗
DOE OSTI2020

Atomic kinetics of laboratory photoionized plasmas relevant to astrophysics. Final report

The major goal of this project was to study the atomic and radiation physics of photoionized plasmas driven by a broadband x-ray flux and its impact on plasma x-ray heating and ionization, emissivity and opacity, and radiation-hydrodynamics. This problem is relevant to the physics of accretion disks surrounding black holes, x-ray binaries, and active galactic nuclei in astrophysics. In addition, it is also important for understanding non-equilibrium laboratory plasmas driven by a distribution of photons. To this end, and motivated by previous and ongoing gas-cell experiments performed by the PI at the Z facility of Sandia National Laboratories, we developed a new experimental platform at the 1MA Zebra pulsed power accelerator of the University of Nevada, Reno based on using a gas jet produced by a supersonic nozzle. The gas jet is turned into a photoionized plasma driven by the x- ray flux from a wire-array z-pinch. Establishing this new experimental platform at Zebra was a major goal of this project. In addition, we have also dedicated effort to continue complementary experiments at the Z facility using the gas-cell photoionized plasma set up with support from the Fundamental Science Program at Z as well as to perform data processing and analysis, and theory and modeling simulations, including electron kinetics and radiation-hydrodynamics. Since we seek to test and establish what physics models are needed to describe photoionized plasma guided by experimental observation, these complementary experimental efforts were central to achieving the goals of this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Design Optimization for a Maglev System Employing Flux Eliminating Coils

Flux eliminating coils have received no little attention over the past thirty years as an alternative for realizing lift in a MAGLEV system. When the magnets on board the vehicle are displaced from the equilibrium or null flux point of these coils, they induce current in those coils which act to restore the coil to its null flux or centerline position. The question being addressed in this paper is that of how to choose the best coil for a given system. What appears at first glance to be an innocent question is in fact one that is actually quite involved, encompassing both the global economics and physics of the system. The real key in analyzing that question is to derive an optimization index or functional which represents the cost of the system subject to constraints, the primary constraint being that the vehicle lift itself at a certain threshold speed. Outlined in this paper is one scenario for realizing a total system design which uses sequential quadratic programming techniques.

Davey, Kent R.↗
NASA NTRS2019

An Analytic Model for Estimating the First Contact Resistance Needed to Avoid Damaging ESD During Spacecraft Docking in GEO

NASA’s Gateway program is to involve spacecraft (s/c) docking in the outer radiation belt in order to transfer Gateway elements between s/c for transport to lunar orbit. The charging of these s/c to different potentials prior to docking raises the possibility of a damaging electrostatic discharge (ESD) at the time of first contact between the s/c. A proposed mitigation strategy is for first contact to occur prior to docking through a resistor with resistance R that would lower the potential difference at an optimal rate to a sufficiently low value to prevent a damaging ESD. The coupling of s/c by a resistor can be modeled by SPIS (Spacecraft Plasma Interaction System), but for realistic two s/c models SPIS can take hours to simulate the evolution of the s/c surface charges and potentials to an equilibrium state. Our objective is to develop a simpler model of s/c resistive coupling that runs orders of magnitude faster while providing useful first design estimates of the time variation of the s/c potentials, current through the resistor, and how these vary with R and s/c configuration. This configuration is defined by the relative separation and orientation of the s/c, and their solar illumination. The configuration and geometry of the s/c determine their capacitive coupling. The s/c capacitances are computed using Nascap-2K. This abstract and the associated poster describe the first version of such a model, and initial tests.

Goodman, Michael L.↗
DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Radionuclide Migration at a Plowshare Program Site: the Gasbuggy Underground Nuclear Test Site, San Juan County, New Mexico - 17019

Project Gasbuggy in northwest New Mexico was the first of three subsurface nuclear tests the U.S. Atomic Energy Commission conducted to evaluate the potential to enhance the production of natural gas from low-permeability geologic formations. In 1967, a 29-kiloton nuclear device was detonated in the Lewis Shale (marine shale and siltstone), 1,292 m below the land surface. The detonation vaporized a cavity into which overlying rock collapsed, forming a 101-m high rubble-filled chimney with a radius of 25 m. Explosion-generated fractures extended between 120 m and 150 m from the working point (detonation location). Most of the radionuclides generated from the explosion remain in the cavity melt glass or have decayed to negligible amounts; tritium radioactivity, however, is still high enough and potentially mobile enough to be of concern. Tritium is able to exchange with stable hydrogen and form liquid water and water vapor; under two-phase (liquid and gas) reservoir conditions the majority of tritium is stored as liquid water and volatilises only to restore equilibrium conditions as gas is removed. A numerical model was developed to examine the migration of tritium in both phases through the fractured shale and sandstone under present conditions. The results suggest that tritium has migrated only 110 m from the working point between 1967 and 2014. Although tritium is able to exchange between the gas and aqueous phases, the faster-diffusing gas phase is responsible for most of the tritium transport from the chimney. A hypothetical producing gas well located 340 m from the working point was added to the model to determine if tritiated water vapor (HTO) could migrate toward the well; results show that, because the low gas velocities could not compete with the high rate of radioactive decay, tritium would not reach the well. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
NASA NTRS1991

Application of program LAURA to thermochemical nonequilibrium flow through a nozzle

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.

Gnoffo, Peter A.↗