A simplified equilibrium hydrocarbon-air combustion gas model for convenient use in air-breathing engine cycle computer programs
Hydrocarbon-air combustion gas model for engine cycle computer programs
Mascitti, V. R.↗
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.
An energy-optimization method to study gel-swelling in confinement
Using our open-source programming environment Morpho, an energy-optimization method is developed to compute equilibrium shapes of an arbitrarily confined swollen hydrogel. The internal strains and contact forces are measured, probing its mechanics.
Joshi, Chaitanya↗
Crater saturation and equilibrium - A Monte Carlo simulation
Conditions of crater saturation and equilibrium have been simulated with a Monte Carlo computer program which utilizes a generating function of the form N = b x D to the -a. The following conclusions have been deduced from the results of the model: (1) saturation and equilibrium for large craters occur at crater densities substantially greater than those observed on the moon, Mars, and Mercury, an indication that the observed populations are essentially production populations, (2) generating functions with a greater than 2 lead to final crater populations with a = 2, (3) generating functions with a between 2 and 1.5 lead to crater populations with the parameter a approaching 1.0, and (4) changes in the slope of crater populations occur gradually and over a wide interval of crater densities and diameters as saturation is approached.
Woronow, A.↗
Dynamic nonlinear analysis of shells of revolution (DYNASOR II)
Equations of motion of shell are solved using Houbolt's numerical procedure with nonlinear terms being moved to right-hand side of equilibrium equations and treated as generalized loads. Program was written in FORTRAN IV for IBM 360 or CDC 6000 series computers.
Tillerson, J. R.↗
Recent update of the RPLUS2D/3D codes
The development of the RPLUS2D/3D codes is summarized. These codes utilize LU algorithms to solve chemical non-equilibrium flows in a body-fitted coordinate system. The motivation behind the development of these codes is the need to numerically predict chemical non-equilibrium flows for the National AeroSpace Plane Program. Recent improvements include vectorization method, blocking algorithms for geometric flexibility, out-of-core storage for large-size problems, and an LU-SW/UP combination for CPU-time efficiency and solution quality.
Tsai, Y.-L. Peter↗
Eleven-Species Thermochemical Model Of Air
SPECIES program computes thermodynamic and transport properties and equilibrium constants. Uses least-squares polynomial curve fits of most accurate data available to provide requested values more quickly than table-lookup methods. Subroutines computing transport coefficients and collision cross sections include additional code to correct for any electron pressure when working with ionic species. Species represented: diatomic nitrogen, diatomic oxygen, atomic nitrogen, atomic oxygen, nitric oxide, ionized nitric oxide, free electron, ionized atomic nitrogen, ionized atomic oxygen, ionized diatomic nitrogen, and ionized diatomic oxygen. Written in FORTRAN 77.
Thompson, Richard A.↗
Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools
The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.
pressurization↗
Computations of non-reacting and reacting viscous blunt body flows, volume 2
Computer programs for calculating the flow distribution in the nose region of a blunt body at arbitrary speed and altitude are discussed. The programs differ from each other in their ability to consider either thin shock or thick shock conditions and in the use of either ideal, equilibrium air, or nonequilibrium air chemistry. The application of the programs to analyzing the flow distribution around the nose of the shuttle orbiter during reentry is reported.
Li, C. P.↗
Extending the thermodynamic form factor bootstrap program: multiple particle-hole excitations, crossing symmetry, and reparameterization invariance
In this study, we further the thermodynamic bootstrap program which involves a set of recently developed ideas used to determine thermodynamic form factors of local operators in integrable quantum field theories. These form factors are essential building blocks for dynamic correlation functions at finite temperatures or non-equilibrium stationary states. In this work we extend this program in three ways. Firstly, we demonstrate that the conjectured annihilation pole axiom is valid in the low energy particle-hole excitations. Secondly, we introduce a crossing relation, which establishes a connection between form factors with different excitation content. Typically, the crossing relation is a consequence of Lorentz invariance, but due to the finite energy density of the considered states, Lorentz invariance is broken. Nonetheless a crossing relation involving excitations with both particles and holes can established using the finite volume representation of the thermodynamic form factors. Finally, we demonstrate that the thermodynamic form factors satisfy a reparameterization invariance, an invariance which encompasses crossing. Reparameterization invariance exploits the fact that the details of the representation of the thermodynamic state are unimportant. In the course of developing these results, we demonstrate the internal consistency of the thermodynamic form factor bootstrap program in a number of ways. Finally, we provide explicit computations of form factors of conserved charges and densities with crossed excitations and show our results can be used to infer information about thermodynamic form factors in the Lieb-Liniger model.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
High Power Flex-Propellant Arcjet Performance
A MW-class electrothermal arcjet based on a water-cooled, wall-stabilized, constricted arc discharge configuration was subjected to extensive performance testing using hydrogen and simulated ammonia propellants with the deliberate aim of advancing technology readiness level for potential space propulsion applications. The breadboard design incorporates alternating conductor/insulator wafers to form a discharge barrel enclosure with a 2.5-cm internal bore diameter and an overall length of approximately 1 meter. Swirling propellant flow is introduced into the barrel, and a DC arc discharge mode is established between a backplate tungsten cathode button and a downstream ringanode/ spin-coil assembly. The arc-heated propellant then enters a short mixing plenum and is accelerated through a converging-diverging graphite nozzle. This innovative design configuration differs substantially from conventional arcjet thrusters, in which the throat functions as constrictor and the expansion nozzle serves as the anode, and permits the attainment of an equilibrium sonic throat (EST) condition. During the test program, applied electrical input power was varied between 0.5-1 MW with hydrogen and simulated ammonia flow rates in the range of 4-12 g/s and 15-35 g/s, respectively. The ranges of investigated specific input energy therefore fell between 50-250 MJ/kg for hydrogen and 10-60 MJ/kg for ammonia. In both cases, observed arc efficiencies were between 40-60 percent as determined via a simple heat balance method based on electrical input power and coolant water calorimeter measurements. These experimental results were found to be in excellent agreement with theoretical chemical equilibrium predictions, thereby validating the EST assumption and enabling the utilization of standard TDK nozzle expansion analyses to reliably infer baseline thruster performance characteristics. Inferred specific impulse performance accounting for recombination kinetics during the expansion process implied nearly frozen flow in the nozzle and yielded performance ranges of 800-1100 sec for hydrogen and 400-600 sec for ammonia. Inferred thrust-to-power ratios were in the range of 30-10 lbf/MWe for hydrogen and 60-20 lbf/MWe for ammonia. Successful completion of this test series represents a fundamental milestone in the progression of high power arcjet technology, and it is hoped that the results may serve as a reliable touchstone for the future development of MW-class regeneratively-cooled flex-propellant plasma rockets.
Litchford, Ron J.↗
A computer program to generate equations of motion matrices, L217 (EOM). Volume 2: Supplemental system design and maintenance document
The equations of motion program L217 (EOM) is described. The program formulates the matrix coefficients for a second order linear differential equation which describes the motion of an airplane relative to its level equilibrium flight condition. Aerodynamic data from FLEXSTAB or Doublet Lattice (L216) programs are used to derive the equations for quasi-steady or complete unsteady aerodynamics.
Clemmons, R. E.↗
Solving Stochastic Dynamic Programming Problems: A Mixed Complementarity Approach
We detail a mixed complementarity problem (MCP) formulation of continuous state dynamic programming problems (DP-MCP). We write the solution to projection methods in value function iteration (VFI) as a joint set of optimality conditions that characterize maximization of the Bellman equation; and approximation of the value function. The MCP approach replaces the iterative component of projection based VFI with a one-shot solution to a square system of complementary conditions. We provide three numerical examples to illustrate our approach.
97 MATHEMATICS AND COMPUTING↗
Optimal Equilibrium Selection of Price-Maker Agents in Performance-Based Regulation Market
This paper analyzes the oligopolistic equilibrium of multiple price-maker agents in performance-based regulation (PBR) markets. In these markets, there are price-maker agents representing some frequency regulation (FR) providers and a number of independent price-taker FR providers. An equilibrium problem with equilibrium constraints (EPEC) model is employed in this paper to study the equilibria of a PBR market in the presence of price-maker agents and pricetaker FR providers. Due to the incorporation of the FR providers' dynamics, the proposed model is reformulated as a mixed-integer linear programming (MILP) problem over innovative mathematical techniques. An optimal equilibrium point is also selected for the market, where none of the agents is unique deviator and the dynamic performance of power system is improved simultaneously. The effectiveness of proposed optimal equilibrium point is evaluated in the numerical results section by comparing the outputs with conventional optimal dispatches of the FR providers.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
A computer program to generate equations of motion matrices, L217 (EOM). Volume 1: Engineering and usage
The equations of motion program L217 formulates the matrix coefficients for a set of second order linear differential equations that describe the motion of an airplane relative to its level equilibrium flight condition. Aerodynamic data from FLEXSTAB or Doublet Lattice (L216) programs can be used to derive the equations for quasi-steady or full unsteady aerodynamics. The data manipulation and the matrix coefficient formulation are described.
Kroll, R. I.↗
Off-of-equilibrium effects on Kurtosis Along Strangeness-Neutral Trajectories
The Beam Energy Scan program at the Relativistic Heavy Ion Collider (RHIC) is searching for the QCD critical point. The main signal for the critical point is the kurtosis of the distribution of proton yields obtained on an event by event basis where one expects a peak at the critical point. However, its exact behavior is still an open question due to out-of-equilibrium effects and uncertainty in the equation of state. Here we use a simplistic hydrodynamic model that enforces strangeness-neutrality, selecting trajectories that pass close to the critical point. We vary the initial conditions to estimate the effect of out-of-equilibrium hydrodynamics on the kurtosis signal.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Spin-selective photocatalysis and quantum transport using ab initio density-matrix dynamics
The overarching goal of this research program is to quantitatively predict out-of-equilibrium quantum dynamics and transport with electrons and spin degrees of freedom coupled to phonons, photons, defects and liquid environments, specifically targeting spin-selective chemistry. Spin polarized carriers, which can be generated at room temperature by circularly-polarized light, provide unique opportunities for creating new reaction pathways and increasing selectivity towards specific pathways. If realized, this can substantially increase the level of reaction control possible in the field of heterogeneous catalysis. Spin, a quantum mechanical entity without a classical counterpart, requires quantum dynamics. Capturing both coherent and incoherent dynamics and accurate correlation of spins, electrons and phonon environments are essential to quantitative description of generation, relaxation and transport, as well as chemical reaction of spin-polarized carriers in molecules and materials promising for spin-photocatalysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Thermal strain analysis of advanced manned spacecraft heat shields first quarterly status report, 1 sep. - 30 nov. 1963
Advanced manned spacecraft heat shield thermal strain analysis - equilibrium and stress equations, singular point, and computer program
HEAT SHIELD↗
Open system magnetic evolution of the taos plateau volcanic field, Northern New Mexico. I - The petrology and geochemistry of the servilleta basalt
MULTIFIT, an embodiment of the conceptual structure needed in modeling multisource and multiprocess magmatic systems, is described. This program, which uses familiar materials balance methodology and the equilibrium form of the Rayleigh equations, links evolutionary arrays, which is turn collectively relate the starting and final compositions of a given magmatic system. Moreover, MULTIFIT incorporates variations within major element data arrays; the linkage between them can be tested using an extension of the least squares algorithm, which selects the best branch point according to the minimum-sum-of-squared-residuals criterion. Advantages and disadvantages of the materials balance approach used in this program are discussed, an example is provided, and equations utilized by MULTIFIT are summarized. While MULTIFIT may not be the best approach for poorly constrained models involving partial melting for complex mixing, it may ultimately prove useful for ascertaining trace element partition coefficients in magnetic systems.
Dungan, M. A.↗