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 199 records · Page 11

Equilibrium gas flow computations. I - Accurate and efficient calculation of equilibrium gas properties

This paper treats the accurate and efficient calculation of thermodynamic properties of arbitrary gas mixtures for equilibrium flow computations. New improvements in the Stupochenko-Jaffe model for the calculation of thermodynamic properties of diatomic molecules are presented. A unified formulation of equilibrium calculations for gas mixtures in terms of irreversible entropy is given. Using a highly accurate thermo-chemical data base, a new, efficient and vectorizable search algorithm is used to construct piecewise interpolation procedures with generate accurate thermodynamic variable and their derivatives required by modern computational algorithms. Results are presented for equilibrium air, and compared with those given by the Srinivasan program.

Liu, Yen↗

Analysis of Aerothermal Heating Over Reusable Single-Stage to Orbit Vehicles

Using the General Aerodynamic Simulation Program (GASP) with finite-rate surface catalysis and radiative equilibrium surface temperature conditions, a computational study has been conducted to predict three-dimensional aero-convective heating over the proposed reusable single-stage to orbit vehicles. The solutions for three basic classes of usable single-stage-to-orbit vehicles, including wing body, lifting body and vertical lander, are presented and compared to provide information for vehicle design. The free stream conditions used in this study occur at an altitude of 57.5 km, ambient density of 4.129e-4 kg/cubic m, and velocity of 4,849 m/s. The convective heating over the vehicle surface in both laminar and turbulent flow regimes were computed. The Baldwin-Lomax algebraic model with a modified van Driest damping factor was applied in the turbulent flow calculations. A bifurcation model which guarantees the conservation of mass diffusive fluxes was adopted to account for the multi-species mass diffusion. The effects of turbulence and surface catalysis on the surface heating are discussed in detail The solutions indicate the degree of oxygen dissociation over the leeward side is sensitive to the shape of vehicle body. Predicted convective heating over the entire surface of lifting body is strongly affected by the presence of surface catalysis. For wing body and vertical lander, the surface catalysis effect does not appear to be important in the far downstream region, and turbulence is the dominant effect on the heating prediction in this region.

Chen, Y. K.↗

Equilibrium and nonequilibrium electron density measurements in transient shock-heated hydrogen-helium plasmas

In the NASA planetary program, atmospheric entry probe missions to Jupiter, Saturn, and Uranus are potential programs now being planned for the 1980's. The work reported in this paper is concerned with predicting the electron density, and resulting radiative emission from the shock layer which must be accommodated by the entry probe heatshield. Within the hydrogen-helium plasma surrounding the forward heatshield, there are regions of near thermochemical equilibrium whose properties are predictable - however, there are uncertainties in these electron density calculations necessitating further measurements by H-beta line broadening and holographic interferometer fringe shift as reported here. There is good agreement between the results obtained from the holographic measurements and the spectroscopic measurements which lends confidence to each of the two completely independent diagnostic techniques. The equilibrium electron density determined from the experimental measurements is somewhat higher than that predicted.

Livingston, F. R.↗

Quantitative Boltzmann–Gibbs Principles via Orthogonal Polynomial Duality

We study fluctuation fields of orthogonal polynomials in the context of particle systems with duality. We thereby obtain a systematic orthogonal decomposition of the fluctuation fields of local functions, where the order of every term can be quantified. This implies a quantitative generalization of the Boltzmann–Gibbs principle. In the context of independent random walkers, we complete this program, including also fluctuation fields in non-stationary context (local equilibrium). For other interacting particle systems with duality such as the symmetric exclusion process, similar results can be obtained, under precise conditions on the n particle dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Doppler-Free Saturation Spectroscopy (DFSS) for Magnetic and Electric Field Measurements in an FRC plasma

The Cooperative Research and Development Agreement (CRADA) between Oak Ridge National Laboratory (ORNL) and TAE Technologies (TAE) focused on a feasibility study of Doppler-free saturation spectroscopy (DFSS) as a diagnostic strategy to measure the internal magnetic and electric fields of TAE’s C-2W field-reversed configuration (FRC) fusion device. Non-perturbative measurements of magnetic and electric fields in the core of FRCs can provide critical information needed to constrain equilibrium models and for the continued advancement of the fusion research program at TAE. DFSS is a well-established, laser-based absorption spectroscopy technique used to suppress inhomogeneous (Doppler) broadening linewidths and reveal Zeeman/Stark split hyperfine structures of target emission lines. The collaboration involved the use of ORNL’s suite of experimentally validated quantum mechanical modelling codes to simulate theoretical atomic spectra for a variety of machine parameters and laser beam configurations. TAE’s expertise was heavily used to provide measurements/models of the required machine parameters needed as inputs to the simulations. Results of an all-encompassing error analysis demonstrate the favourable capabilities of the diagnostic and provided invaluable insights as to how DFSS can be best used to accurately measure internal B/E-fields in C-2W.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of Magnetic Field using Doppler-Free Saturation Spectroscopy (DFSS) in C-2W FRC Plasma

The Cooperative Research and Development Agreement (CRADA) between Oak Ridge National Laboratory (ORNL) and TAE Technologies (TAE) focused on implementing an existing mobile Doppler-free saturation spectroscopy (DFSS) diagnostic to determine the direction of the equilibrium magnetic field within TAE’s C-2W field-reversed configuration (FRC) fusion device. Non-perturbative measurements of magnetic fields in the core of FRCs can provide critical information needed to constrain equilibrium models and for the continued advancement of the fusion research program at TAE. DFSS is a well-established, laser-based absorption spectroscopy technique used to suppress inhomogeneous (Doppler) broadening and reveal the Zeeman fine structure of the target emission lines. The collaboration involved the use of ORNL’s recently constructed mobile DFSS diagnostic. The DFSS diagnostic was shipped and installed on C- 2W in February 2023. TAE’s expertise was used to design and construct the mechanic support structures required to mount the optical hardware on C-2W. Additionally TAE provide critical input to minimize the noise imposed on the DFSS signal by vibrations of the mechanic support structure during operations. DFSS data was obtained for over 500 C-2W shots and will continue until the expiration of this CRADA on September 30 th 2023. Over the course of this work, significant progress was made to maximize the signal to noise (SNR) ratio. Analysis of recently acquired data is on-going, however, initial indications are highly suggestive that a DFSS signal above background has been successfully acquired in C-2W.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FORTRAN 4 computer program for calculation of thermodynamic and transport properties of complex chemical systems

A FORTRAN IV computer program for the calculation of the thermodynamic and transport properties of complex mixtures is described. The program has the capability of performing calculations such as:(1) chemical equilibrium for assigned thermodynamic states, (2) theoretical rocket performance for both equilibrium and frozen compositions during expansion, (3) incident and reflected shock properties, and (4) Chapman-Jouguet detonation properties. Condensed species, as well as gaseous species, are considered in the thermodynamic calculation; but only the gaseous species are considered in the transport calculations.

Svehla, R. A.↗

Combustion of hydrogen injected into a supersonic airstream (the SHIP computer program)

The mathematical and physical basis of the SHIP computer program which embodies a finite-difference, implicit numerical procedure for the computation of hydrogen injected into a supersonic airstream at an angle ranging from normal to parallel to the airstream main flow direction is described. The physical hypotheses built into the program include: a two-equation turbulence model, and a chemical equilibrium model for the hydrogen-oxygen reaction. Typical results for equilibrium combustion are presented and exhibit qualitatively plausible behavior. The computer time required for a given case is approximately 1 minute on a CDC 7600 machine. A discussion of the assumption of parabolic flow in the injection region is given which suggests that improvement in calculation in this region could be obtained by use of the partially parabolic procedure of Pratap and Spalding. It is concluded that the technique described herein provides the basis for an efficient and reliable means for predicting the effects of hydrogen injection into supersonic airstreams and of its subsequent combustion.

Markatos, N. C.↗

Nuclear Reaction Data Evaluations for Nuclear Applications Physics and modeling [Slides]

This presentation provides an evaluation workflow scheme for the project. It also discusses nuclear data libraries, including an overview of worldwide nuclear data libraries and current funding agencies. The presentation also discusses the theoretical framework commonly used in the evaluated work, including the R-matrix theory (≈ up to a few keVs for heavy nuclei) as well as pre-equilibrium (PE) and Hauser-Feshbach (HF) theory (⪆1 MeV for medium/heavy nuclei). Inclusion of experimental effects are also discussed. Examples of evaluated work, including theory, experiments, and applications are discussed, specifically n+ 28,29,30 Si evaluation (direct capture component), n+ 53 Cr evaluation (strong multiple scattering effects), n+ 239 Pu(n,γf) reaction (fluctuations in neutron multiplicities), and α+ 17,18 O evaluations (emerging neutron spectra).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Non-Equilibrium Thermodynamic Chemistry and the Composition of the Atmosphere of Mars

A high priority objective of the Mars Exploration Program is to Determine if life exists today (MEPAG Goal I, Objective A). The measurement of gases of biogenic origin may be an approach to detect the presence of microbial life on the surface or subsurface of Mars. Chemical thermodynamic calculations indicate that on both Earth and Mars, certain gases should exist in extremely low concentrations, if at all. Microbial metabolic activity is an important non-equilibrium chemistry process on Earth, and if microbial life exists on Mars, may be an important nonequilibrium chemistry process on Mars. The non-equilibrium chemistry of the atmosphere of Mars is discussed in this paper.

J S Levine↗

Calculation of radiative properties of nonequilibrium hydrogen plasma

A computer program called NEQRAP is described that calculates the radiative properties of nonequilibrium ionized hydrogen. From the given electron temperature, electron density, and atom density values (which do not necessarily satisfy the equilibrium relationship) and intensities of incident radiation, the non-Boltzmann populations of electronic states are computed by solving the equation of quasi-steady-state population distribution. Emission and absorption coefficients are determined as functions of wavelength by invoking the principle of detailed balance between the upper and lower states of each radiative transition. Radiative transport through the medium is computed assuming a one-dimensional uniform slab. The rate of ionic reaction is also computed. When used on a sample case, the program shows that there is a large difference between the calculated intensities of radiation emitted by a bulk of equilibrium and nonequilibrium hydrogen. The accuracy of the program is estimated to be better than 10%.

Park, C.↗

RAMSCRAM: A flexible ramjet/scramjet engine simulation program

With the resurgence of interest in high supersonic and hypersonic flight there is a need to simulate airbreathing engines which may be used in this flight regime. To meet this requirement the RAMSCRAM code was developed. The code calculates 1-D flow properties at each component interface and the overall performance of the engine. It uses equilibrium thermodynamics which accounts for dissociation and allows for any fuel or combination of fuels. The program can simulate ramjet, scramjet, rocket, and ducted rocket engines.

Burkardt, Leo A.↗

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗

Thermodynamic Investigation of the Effect of Target Compositions and Conditions of Their Magnetron Sputtering on Phase Compositions of the Resulting Films

A thermodynamic investigation of plasmochemical processes taking place in a magnetron during sputtering of the targets consisting of indium and dysprosium ions is performed, which allowed identifying the feasibility of forming under the given conditions certain amounts of the chemical compounds In{sub 2}O{sub 3}, Dy{sub 2}O{sub 3}, DyInO{sub 3}. In the thermodynamic analysis of the equilibrium compositions during magnetron sputtering of the films, a universal TERRA software program is used.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of neutron beamline irradiation on hydrogen equilibrium vapor pressure over hydrided Zircaloy-4

Here, the effect of neutron irradiation on the hydrogen equilibrium vapor pressure over hydrided samples of nickel-plated Zircaloy-4 was evaluated at the Washington State University Research Reactor’s epithermal beamline facility. To enable this study, a measurement capability was designed and assembled to obtain low pressure (<0.1 Pa) measurements over the temperature range of interest (i.e., ~300° - 400°C). Measurement results were in good agreement with tends extrapolated from previous data obtained at higher temperatures as well as recent thermodynamic data. No irradiation enhancement of vapor pressure was observed under the incident flux of epithermal neutrons over this temperature range.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration

In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ~6–8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. Finally, the combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polynomial approximations of thermodynamic properties of arbitrary gas mixtures over wide pressure and density ranges

Computer programs for flow fields around planetary entry vehicles require real-gas equilibrium thermodynamic properties in a simple form which can be evaluated quickly. To fill this need, polynomial approximations were found for thermodynamic properties of air and model planetary atmospheres. A coefficient-averaging technique was used for curve fitting in lieu of the usual least-squares method. The polynomials consist of terms up to the ninth degree in each of two variables (essentially pressure and density) including all cross terms. Four of these polynomials can be joined to cover, for example, a range of about 1000 to 11000 K and 0.00001 to 1 atmosphere (1 atm = 1.0133 x 100,000 N/m sq) for a given thermodynamic property. Relative errors of less than 1 percent are found over most of the applicable range.

Allison, D. O.↗