Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computable general equilibrium”

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 145 records · Page 8

Estimated Performance of Radial-Flow Exit Nozzles for Air in Chemical Equilibrium

The thrust, boundary-layer, and heat-transfer characteristics were computed for nozzles having radial flow in the divergent part. The working medium was air in chemical equilibrium, and the boundary layer was assumed to be all turbulent. Stagnation pressure was varied from 1 to 32 atmospheres, stagnation temperature from 1000 to 6000 R, and wall temperature from 1000 to 3000 R. Design pressure ratio was varied from 5 to 320, and operating pressure ratio was varied from 0.25 to 8 times the design pressure ratio. Results were generalized independent of divergence angle and were also generalized independent of stagnation pressure in the temperature range of 1000 to 3000 R. A means of determining the aerodynamically optimum wall angle is provided.

Englert, Gerald W.↗

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↗

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

This work presents a general framework for model reduction of non-equilibrium energy transfer and dissociation processes. The multi-group maximum-entropy method is coupled with the quasi-classical trajectory method to directly construct a reduced order model for chemical non-equilibrium. Kinetic data is calculated by applying detailed balance at a microscopic level, overcoming the limitation of calculating recombination kinetic data. This approach enables the construction of a reduced order model for kinetics which bypasses the need to compute state-to-state kinetic data. This physics based reduced order model ensures that an equilibrium distribution is reached given infinite time, while allowing for non-equilibrium distributions during the relaxation and dissociation processes. A proof-of-concept test case demonstrates the applicability of this model by comparison with state-to-state kinetic data for the N2(Chi (sup 1) Sigma ((sub g), (sup plus sign)) minus N ((exp 4) S (sub u)) system with excellent agreement.

Maximum Entropy↗

Aerothermodynamic Analysis of Stardust Sample Return Capsule with Coupled Radiation and Ablation

An aerothermodynamic analysis of the forebody aeroshell of the Stardust Sample Return Capsule is carried out by using the axisymmetric viscous shock-layer equations with and without fully coupled radiation and ablation. Formulation of the viscous shock-layer equations with shoulder radius as the length scale and implementation of the Vigneron pressure condition allow resolution of the flowfield over the shoulder. With a predominantly supersonic outflow over the shoulder, a globally iterated solution or viscous shock-layer equations can be obtained. The stagnation-point results are obtained along a specified trajectory, whereas detailed calculations along the body are provided at the peak-heating point. The equilibrium calculations with ablation injection are the focus of the present study because of the lack of a general chemical nonequilibrium analysis that accounts for both surface and flowfield effect. The equilibrium calculations also provide a simple way to conserve surface (and flowfield) elemental composition for the current small ablation injection rates, where the surface elemental composition is a mixture of freestream and ablator elements. Therefore, the coupled laminar and turbulent flow solutions with radiation and ablation are obtained by using the equilibrium flow chemistry, whereas a nonequilibrium chemistry model is used for solutions without ablation and turbulence. Various computed results are compared with those obtained by the other researchers.

Gupta, Roop N.↗

A family of dynamic models for large-eddy simulation

Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Carati, D.↗

Plotting component maps in the Navy/NASA Engine Program (NNEP): A method and its usage

The Navy/NASA Engine Program (NNEP) and the new extended version which handles chemical equilibrium (NNEPEQ) are very general cycle analysis codes that have been used extensively to calculate design and off-design performance of a wide range of turbine engine cycles and configurations. Component maps are used to obtain the off-design engine performance and a matched engine cycle. This paper describes a method of plotting the scaled NNEP compressor and turbine maps as the user runs the NNEP code as well as plotting the operating line defined by all the cases that were computed in that particular NNEP run. Afterwards, an example demonstrates the use of this capability to help analyze an engine cycle model and then make improvements to that cycle.

Plencner, Robert M.↗

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗

Using Multiplicative Algebraic Reconstruction Techniques (MART) to Derive Instrument Requirements for Computed Tomography Imaging Spectrographs

Computed Tomography Imaging Spectrographs (CTISs), which are generally slitless or large-aperture spectrographs that observe a wide, dispersed field of view in multiple dispersion angles or diffraction orders, provide a unique opportunity to capture spectral information over an extended source, such as the Sun, but require software reconstruction techniques to be fully utilized. A useful data product that can be recovered from CTIS observations is intensity maps of solar features in single spectral lines. These intensity maps can then be used to determine temperature, density, abundance, and equilibrium properties of the emitting plasma. We apply the multiplicative algebraic reconstruction technique (MART) to example data from a variety of CTIS configurations to determine the capability and limitations of the method to return spectrally-pure maps. By completing this study, we aim to establish a path to derive requirements for CTIS instruments that rely on reconstruction techniques to meet their science objectives.

Amy Winebarger↗

Numerical techniques in radiative heat transfer for general, scattering, plane-parallel media

The study of radiative heat transfer with scattering usually leads to the solution of singular Fredholm integral equations. The present paper presents an accurate and efficient numerical method to solve certain integral equations that govern radiative equilibrium problems in plane-parallel geometry for both grey and nongrey, anisotropically scattering media. In particular, the nongrey problem is represented by a spectral integral of a system of nonlinear integral equations in space, which has not been solved previously. The numerical technique is constructed to handle this unique nongrey governing equation as well as the difficulties caused by singular kernels. Example problems are solved and the method's accuracy and computational speed are analyzed.

Sharma, A.↗

Turbulence modeling for non-equilibrium flow

The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Durbin, P. A.↗

Hypersonic blunt body computations including real gas effects

The recently developed second-order explicit and implicit total variation diminishing (TVD) shock-capturing methods of the Harten and Yee, Yee, and van Leer types in conjunction with a generalized Roe's approximate Riemann solver of Vinokur and the generalized flux-vector splittings of Vinokur and Montagne for two-dimensional hypersonic real gas flows are studied. A previous study on one-dimensional unsteady problems indicated that these schemes produce good shock-capturing capability and that the state equation does not have a large effect on the general behavior of these methods for a wide range of flow conditions for equilibrium air. The objective of this paper is to investigate the applicability and shock resolution of these schemes for two-dimensional steady-state hypersonic blunt body flows. The main contribution of this paper is to identify some of the elements and parameters which can affect the convergence rate for high Mach numbers or real gases but have negligible effect for low Mach number cases for steady-state inviscid blunt body flows.

Montagne, J.-L.↗

Nonlinear Dynamic Models in Advanced Life Support

To facilitate analysis, ALS systems are often assumed to be linear and time invariant, but they usually have important nonlinear and dynamic aspects. Nonlinear dynamic behavior can be caused by time varying inputs, changes in system parameters, nonlinear system functions, closed loop feedback delays, and limits on buffer storage or processing rates. Dynamic models are usually cataloged according to the number of state variables. The simplest dynamic models are linear, using only integration, multiplication, addition, and subtraction of the state variables. A general linear model with only two state variables can produce all the possible dynamic behavior of linear systems with many state variables, including stability, oscillation, or exponential growth and decay. Linear systems can be described using mathematical analysis. Nonlinear dynamics can be fully explored only by computer simulations of models. Unexpected behavior is produced by simple models having only two or three state variables with simple mathematical relations between them. Closed loop feedback delays are a major source of system instability. Exceeding limits on buffer storage or processing rates forces systems to change operating mode. Different equilibrium points may be reached from different initial conditions. Instead of one stable equilibrium point, the system may have several equilibrium points, oscillate at different frequencies, or even behave chaotically, depending on the system inputs and initial conditions. The frequency spectrum of an output oscillation may contain harmonics and the sums and differences of input frequencies, but it may also contain a stable limit cycle oscillation not related to input frequencies. We must investigate the nonlinear dynamic aspects of advanced life support systems to understand and counter undesirable behavior.

Jones, Harry↗

Systematic generation of nonlinear discretized dynamic equilibrium equations of spinning cantilevers

General nonlinear discretized governing equations of motion of spinning elastic solids and structures are adjusted for the case of a spinning cantilever with initial geometric imperfections. Consideration is given to second degree nonlinearities in the strain-displacement and velocity-displacement relationships. Parameters of the discretization are developed to include the type and number of the coordinate functions used in the admissible trial solution in order to unify the discretization approaches associated with stationarity principles. The coordinate functions comprise both sets of continuous and piecewise continuous functions employed in the Rayleigh-Ritz and the finite element methods, respectively. Coefficient matrices are provided which contain the energy density expressions and which are adaptable to computer programming.

El-Essawi, M.↗

Jacchia-Lineberry upper atmosphere density model

A series of increasingly accurate models which are a careful blend of empirical and theoretical formulae were developed. The exospheric temperature is assumed to be a function of: (1) the average and daily variations in the solar flux, (2) the average and three hourly variations in the geomagnetic index, (3) the angle between the position vector and the axis of the unsymmetric atmospheric bulge, and (4) the angle between the position vector and the geomagnetic pole. The exospheric temperature is related to the density by the solution of the diffusion equilibrium equations for the different constituents of the atmosphere as a function of altitude. Other variations are modeled directly as changes in the density. They are: (1) changes due to the semiannual effect, and (2) changes due to the seasonal latitudinal effect. The causes for these variations are not exactly known but may be modeled sufficiently by empirical formulae. The Jacchia model is assumed to be valid over the altitude range of 90 to 2500 km. The residuals between the observed density from satellite drag observations and the computed densities show the mean relative error to be generally less than 10 percent with occasional peak errors near 50 percent.

Mueller, A. C.↗

Viscous equilibrium computations using program LAURA

Modifications have been made to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) that enable it to compute viscous airflows under the assumption of thermal and chemical equilibrium. Equilibrium thermodynamic and transport property information are input to the code via curve fits. The periodic updating of this information enables the equilibrium algorithm to perform at a computational rate that is only a small percentage larger than the rate associated with the perfect-gas algorithm. Presented in this article are the results of the initial validation of the modified code. Solutions for surface pressure and heating are presented for the flow over slender and blunt cones at realistic reentry conditions. LAURA solutions are compared with those produced by a viscous shock-layer method, and, for one case considered, with heat transfer data from a flight experiment. For both pressure and heating, the agreement is good. In general, differences in pressure of a few percent were noted, while differences in heating rates were in the 5-10 percent range.

Greene, Francis A.↗

A general stagnation-point convective heating equation for arbitrary gas mixtures

The stagnation-point convective heat transfer to an axisymmetric blunt body for arbitrary gases in chemical equilibrium was investigated. The gases considered were base gases of nitrogen, oxygen, hydrogen, helium, neon, argon, carbon dioxide, ammonia, and methane and 22 gas mixtures composed of the base gases. Enthalpies ranged from 2.3 to 116.2 MJ/kg, pressures ranged from 0.001 to 100 atmospheres, and the wall temperatures were 300 and 1111 K. A general equation for the stagnation-point convective heat transfer in base gases and gas mixtures was derived and is a function of the mass fraction, the molecular weight, and a transport parameter of the base gases. The relation compares well with present boundary-layer computer results and with other analytical and experimental results. In addition, the analysis verified that the convective heat transfer in gas mixtures can be determined from a summation relation involving the heat transfer coefficients of the base gases. The basic technique developed for the prediction of stagnation-point convective heating to an axisymmetric blunt body could be applied to other heat transfer problems.

Sutton, K.↗

Nonlinear analysis of compressively loaded linked-plate structures

A method is presented for the geometrically nonlinear analysis of the equilibrium behavior of a general class of compressively loaded, prismatic, linked-plate structures. The semi-analytical method accounts for geometric shape imperfections, and allows the use of laminated composite plate properties. Buckling eigenfunctions are used as displacement shape functions, and the associated second-order displacement fields are incorporated. The nonlinear plate equilibrium equations are expressed in terms of the assumed displacement form, and the Galerkin method is used to obtain a set of algebraic equations which are nonlinear in the modal amplitudes and the load parameter. Buckling eigensolutions for the class of structures under consideration can be obtained from the VIPASA computer code; thus, subject to some additional restrictions in loading and plate anisotropy, structural configurations which can be analyzed with regard to buckling response by VIPASA can be analyzed with regard to geometrically nonlinear response using the new method.

Stoll, Frederick↗