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 73 records · Page 4

EMISApproximateEquilibrium.jl [SWR-19-56]

The Electricity Markets Investment Suite Approximate Equilibrium (EMIS-AE) package is developed at NREL and provides a solution capability to solve generation expansion equilibrium problems in the electricity markets. EMIS-AE is designed to capture the evolution of the electricity generation portfolio resulting from the interactions of heterogeneous investors under different policy and market designs. The investment problem for each generation company (GENCO) is a bi-level problem with the investment decision made in the upper level and market clearing condition in the lower level, which traditionally is represented as a Mathematical Program with Equilibrium Constraint (MPEC). EMIS-AE provides a predictive model to be trained for estimating the system-wide revenues for each technology type across energy, ancillary services and capacity markets given the amount of installed capacity on the grid. The profit maximization investment problem for each GENCO is solved using a global search algorithm, which uses the predictive model to evaluate the objective function. To solve for the strategic equilibrium, each GENCO’s problem is plugged into a diagonalization algorithm that is generally used in multi-leader, single-follower bi-level problems.

Dalvi, Sourabh↗

Resource Adequacy in Electricity Markets With Renewable Energy

This paper presents a market-based resource adequacy assessment framework to analyze the system generation portfolio that results in a competitive market environment. In particular, we apply the modeling framework to investigate the impact of high penetration levels of variable renewable energy resources and different market designs on achieving resource adequacy. The model captures the strategic capacity investment and retirement decision-making of profit-maximizing generation companies. In contrast to many previous studies, this paper considers markets for capacity, energy, and reserve products in detail; thus, the modeling framework allows improved analysis of generation expansion and revenue sufficiency in a competitive market environment. The modeling framework is based on Stackelberg leader-follower games and is formulated as a bi-level optimization problem, which is then transformed into a mathematical program with equilibrium constraints. Furthermore, we employ a Lagrangian decomposition algorithm to enhance computational performance. In a case study of the electric reliability council of texas (ERCOT) system, we compare the results from the proposed model with those obtained from a traditional centralized least-cost planning model. Finally, the results demonstrate the effectiveness of the proposed modeling framework and highlight the importance of considering strategic behavior in a competitive market framework when assessing resource adequacy.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

The static nonlinear analysis of shells of revolution (SNASOR II)

Utilizing stiffness matrices and supplying as input the loading and boundary conditions, program generates equilibrium equations for structure. Nonlinear strain energy terms result in pseudogeneralized forces which are combined with applied generalized forces. Resulting set of nonlinear algebraic equilibrium equations is solved by one of several methods.

Stricklin, J. A.↗

Polynomial coefficients of thermochemical data for the C-H-O-N system

Thermodynamic data are required input for the finite kinetics and equilibrium computer programs needed for modeling the combustion of hydrocarbons in the fields of energy and pollution research. Least squares determined coefficients of the curve-fitted thermodynamic data for 193 species in the C-H-O-N system are presented in card image form and are of suitable format for use by common computer programs.

Wakelyn, N. T.↗

NEQAIR96,Nonequilibrium and Equilibrium Radiative Transport and Spectra Program: User's Manual

This document is the User's Manual for a new version of the NEQAIR computer program, NEQAIR96. The program is a line-by-line and a line-of-sight code. It calculates the emission and absorption spectra for atomic and diatomic molecules and the transport of radiation through a nonuniform gas mixture to a surface. The program has been rewritten to make it easy to use, run faster, and include many run-time options that tailor a calculation to the user's requirements. The accuracy and capability have also been improved by including the rotational Hamiltonian matrix formalism for calculating rotational energy levels and Hoenl-London factors for dipole and spin-allowed singlet, doublet, triplet, and quartet transitions. Three sample cases are also included to help the user become familiar with the steps taken to produce a spectrum. A new user interface is included that uses check location, to select run-time options and to enter selected run data, making NEQAIR96 easier to use than the older versions of the code. The ease of its use and the speed of its algorithms make NEQAIR96 a valuable educational code as well as a practical spectroscopic prediction and diagnostic code.

Whiting, Ellis E.↗

Thermal/Pyrolysis Gas Flow Analysis of Carbon Phenolic Material

Provided in this study are predicted in-depth temperature and pyrolysis gas pressure distributions for carbon phenolic materials that are externally heated with a laser source. Governing equations, numerical techniques and comparisons to measured temperature data are also presented. Surface thermochemical conditions were determined using the Aerotherm Chemical Equilibrium (ACE) program. Surface heating simulation used facility calibrated radiative and convective flux levels. Temperatures and pyrolysis gas pressures are predicted using an upgraded form of the SINDA/CMA program that was developed by NASA during the Solid Propulsion Integrity Program (SPIP). Multispecie mass balance, tracking of condensable vapors, high heat rate kinetics, real gas compressibility and reduced mixture viscosity's have been added to the algorithm. In general, surface and in-depth temperature comparisons are very good. Specie partial pressures calculations show that a saturated water-vapor mixture is the main contributor to peak in-depth total pressure. Further, for most of the cases studied, the water-vapor mixture is driven near the critical point and is believed to significantly increase the local heat capacity of the composite material. This phenomenon if not accounted for in analysis models may lead to an over prediction in temperature response in charring regions of the material.

Clayton, J. Louie↗

NEQAIR v15.0 Release Notes: Nonequilibrium and Equilibrium Radiative Transport and Spectra Program

NEQAIR v15.0 provides the first steps to improved coupling between NEQAIR and the DPLR CFD code, which will be fully realized in v15.1. The plan is to release NEQAIR v15.1 and DPLR 4.05 at the same time. The improvements implemented in NEQAIR v15.0 have focused on improving stability, solution robustness, usability and providing different options for running the code. It is also the first version of the code to have a new input file and line of sight format since 2009. Backward compatibility with previous formats of the input files (neqair.inp and LOS.dat) has also been provided. NEQAIR v15.0 supersedes the prerelease of this version, as well as NEQAIR v14.0, v13.2, v13.1 and the suite of NEQAIR2009 versions. These updates have predominantly been performed by Brett Cruden and Aaron Brandis from AMA Inc at NASA Ames Research Center between 2016 and 2018. NEQAIR v15.0 is a standalone software tool for line-by-line spectral computation of radiative intensities and/or radiative heat flux, with one-dimensional transport of radiation. In order to accomplish this, NEQAIR v15.0, as in previous versions, requires the specification of distances (in cm), temperatures (in K) and number densities (in parts/cc) of constituent species along lines of sight. Therefore, it is assumed that flow quantities have been extracted from flow fields computed using other tools, such as CFD codes like DPLR or LAURA, and that lines of sight have been constructed and written out in the format required by NEQAIR v15.0. There are two principal modes for running NEQAIR v15.0. In the first mode NEQAIR v15.0 is used as a tool for creating synthetic spectra of any desired resolution (including convolution with a specified instrument/slit function). The first mode is typically exercised in simulating/interpreting spectroscopic measurements of different sources (e.g. shock tube data, plasma torches, etc.). In the second mode, NEQAIR v15.0 is used as a radiative heat flux prediction tool for flight projects. Correspondingly, NEQAIR has also been used to simulate the radiance measured on previous flight missions. This report summarizes the database updates, corrections that have been made to the code, changes to input files, parallelization, the current usage recommendations, including test cases, and an indication of the performance enhancements achieved.

Brandis, Aaron M.↗

A General Method for Automatic Computation of Equilibrium Compositions and Theoretical Rocket Performance of Propellants

A general computer program for chemical equilibrium and rocket performance calculations was written for the IBM 650 computer with 2000 words of drum storage, 60 words of high-speed core storage, indexing registers, and floating point attachments. The program is capable of carrying out combustion and isentropic expansion calculations on a chemical system that may include as many as 10 different chemical elements, 30 reaction products, and 25 pressure ratios. In addition to the equilibrium composition, temperature, and pressure, the program calculates specific impulse, specific impulse in vacuum, characteristic velocity, thrust coefficient, area ratio, molecular weight, Mach number, specific heat, isentropic exponent, enthalpy, entropy, and several thermodynamic first derivatives.

Gordon, Sanford↗

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗

Computer program to solve two-dimensional shock-wave interference problems with an equilibrium chemically reacting air model

The computer program EASI, an acronym for Equilibrium Air Shock Interference, was developed to calculate the inviscid flowfield, the maximum surface pressure, and the maximum heat flux produced by six shock wave interference patterns on a 2-D, cylindrical configuration. Thermodynamic properties of the inviscid flowfield are determined using either an 11-specie, 7-reaction equilibrium chemically reacting air model or a calorically perfect air model. The inviscid flowfield is solved using the integral form of the conservation equations. Surface heating calculations at the impingement point for the equilibrium chemically reacting air model use variable transport properties and specific heat. However, for the calorically perfect air model, heating rate calculations use a constant Prandtl number. Sample calculations of the six shock wave interference patterns, a listing of the computer program, and flowcharts of the programming logic are included.

Glass, Christopher E.↗

NNEPEQ: Chemical equilibrium version of the Navy/NASA Engine Program

The Navy NASA Engine Program, NNEP, currently is in use at a large number of government agencies, commercial companies and universities. This computer code has bee used extensively to calculate the design and off-design (matched) performance of a broad range of turbine engines, ranging from subsonic turboprops to variable cycle engines for supersonic transports. Recently, there has been increased interest in applications for which NNEP was not capable of simulating, namely, high Mach applications, alternate fuels including cryogenics, and cycles such as the gas generator air-turbo-rocker (ATR). In addition, there is interest in cycles employing ejectors such as for military fighters. New engine component models had to be created for incorporation into NNEP, and it was found necessary to include chemical dissociation effects of high temperature gases. The incorporation of these extended capabilities into NNEP is discussed and some of the effects of these changes are illustrated.

Fishbach, Laurence H.↗

NNEPEQ - Chemical equilibrium version of the Navy/NASA Engine Program

The Navy NASA Engine Program, NNEP, currently is in use at a large number of government agencies, commercial companies and universities. This computer code has been used extensively to calculate the design and off-design (matched) performance of a broad range of turbine engines, ranging from subsonic turboprops to variable cycle engines for supersonic transports. Recently, there has been increased interest in applications for which NNEP was not capable of simulating, namely, high Mach applications, alternate fuels including cryogenics, and cycles such as the gas generator air-turbo-rocker (ATR). In addition, there is interest in cycles employing ejectors such as for military fighters. New engine component models had to be created for incorporation into NNEP, and it was found necessary to include chemical dissociation effects of high temperature gases. The incorporation of these extended capabilities into NNEP is discussed and some of the effects of these changes are illustrated.

Fishbach, L. H.↗

Transpiration and film cooling boundary layer computer program. Volume 1: Numerical solutions of the turbulent boundary layer equations with equilibrium chemistry

A finite difference turbulent boundary layer computer program has been developed. The program is primarily oriented towards the calculation of boundary layer performance losses in rocket engines; however, the solution is general, and has much broader applicability. The effects of transpiration and film cooling as well as the effect of equilibrium chemical reactions (currently restricted to the H2-O2 system) can be calculated. The turbulent transport terms are evaluated using the phenomenological mixing length - eddy viscosity concept. The equations of motion are solved using the Crank-Nicolson implicit finite difference technique. The analysis and computer program have been checked out by solving a series of both laminar and turbulent test cases and comparing the results to data or other solutions. These comparisons have shown that the program is capable of producing very satisfactory results for a wide range of flows. Further refinements to the analysis and program, especially as applied to film cooling solutions, would be aided by the acquisition of a firm data base.

Levine, J. N.↗

Combustion of hydrogen-air jets in local chemical equilibrium: A guide to the CHARNAL computer program

A guide to a computer program, written in FORTRAN 4, for predicting the flow properties of turbulent mixing with combustion of a circular jet of hydrogen into a co-flowing stream of air is presented. The program, which is based upon the Imperial College group's PASSA series, solves differential equations for diffusion and dissipation of turbulent kinetic energy and also of the R.M.S. fluctuation of hydrogen concentration. The effective turbulent viscosity for use in the shear stress equation is computed. Chemical equilibrium is assumed throughout the flow.

Spalding, D. B.↗

Program Computes Performances Of Scramjets

SRGULL, a scramjet-cycle-analysis computer code, engineer's software tool capable of nose-to-tail simulation of hydrogen-fueled, scramjet engine integrated with air-frame. Simulated flow that of real gas with equilibrium thermodynamic properties. Program facilities initial estimates of performance of scramjet cycle by linking code for two-dimensional forebody, inlet, and nozzle with code for one-dimensional combustor. Five computer codes (SCRAM, SEAGUL, INLET, Program HUD, and GASH) integrated in program to provide capability for analysis of changing flow conditions. Written in FORTRAN 77.

Walton, J. T.↗