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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Monte Carlo Hauser-Feshbach computer code system to model nuclear reactions: YAHFC

A computer program framework, YAHFC, to model low-energy nuclear reactions is presented. The framework allows for reactions with incident particles ranging from protons/neutrons to alphas and is designed to address reactions that ultimately lead to the formation of compound nuclear systems that then decay statistically as outlined in concepts of Hauser and Feshbach. Additionally, instead of a reaction, it is also possible to model the decay of a nuclear system with an initial excitation and population. The code models nuclear decays with a Monte Carlo process that tracks the decay of each state. This allows for an exact representation of the spectra for all emitted particles in each of the final exit channels and the possibility of generating reaction data for simulation purposes. The program is interfaced with the optical model code system FRESCOX to calculate transmission coefficients as well as the effects of coupled channels and other direct excitations via the distorted wave Born approximation (DWBA). Modules are included to account for nuclear processes such as width corrections, pre-equilibrium emission, and fission. The program is controlled by a series of input commands and while a set of input parameters exists for each projectile and target, the input commands allow for complete control over each input parameter. Extensive data files are produced and a program is provided that converts YAHFC data files into nuclear data library entries in the generalized nuclear data structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Perturbative approaches in relativistic kinetic theory and the emergence of first-order hydrodynamics

Hydrodynamics can be formulated in terms of a perturbative series in derivatives of the temperature, chemical potential, and flow velocity around an equilibrium state. Different formulations for this series have been proposed over the years, which consequently led to the development of various hydrodynamic theories. In this work, we discuss the relativistic generalizations of the perturbative expansions put forward by Chapman and Enskog, and Hilbert, using general matching conditions in kinetic theory. This allows us to describe, in a comprehensive way, how different out-of-equilibrium definitions for the hydrodynamic fields affect the development of the hydrodynamic perturbative series. We provide a perturbative method for systematically deriving the hydrodynamic formulation recently proposed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK) from relativistic kinetic theory. The various transport coefficients that appear in BDNK (at first-order) are explicitly computed using a new formulation of the relaxation time approximation for the Boltzmann equation. Assuming Bjorken flow, we also determine the hydrodynamic attractors of BDNK theory and compare the overall hydrodynamic evolution obtained using this formulation with that generated by the Israel-Stewart equations of motion and also kinetic theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

The effects of finite rate chemical processes on high enthalpy nozzle performance - A comparison between SPARK and SEAGULL

The generalized chemistry version of the computer code SPARK is extended to include two higher-order numerical schemes, yielding fourth-order spatial accuracy for the inviscid terms. The new and old formulations are used to study the influences of finite rate chemical processes on nozzle performance. A determination is made of the computationally optimum reaction scheme for use in high-enthalpy nozzles. Finite rate calculations are compared with the frozen and equilibrium limits to assess the validity of each formulation. In addition, the finite rate SPARK results are compared with the constant ratio of specific heats (gamma) SEAGULL code, to determine its accuracy in variable gamma flow situations. Finally, the higher-order SPARK code is used to calculate nozzle flows having species stratification. Flame quenching occurs at low nozzle pressures, while for high pressures, significant burning continues in the nozzle.

Carpenter, M. H.↗
NASA NTRS1995

Program Helps To Determine Chemical-Reaction Mechanisms

General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code developed for use in solving complex, homogeneous, gas-phase, chemical-kinetics problems. Provides for efficient and accurate chemical-kinetics computations and provides for sensitivity analysis for variety of problems, including problems involving honisothermal conditions. Incorporates mathematical models for static system, steady one-dimensional inviscid flow, reaction behind incident shock wave (with boundary-layer correction), and perfectly stirred reactor. Computations of equilibrium properties performed for following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. Written in FORTRAN 77 with exception of NAMELIST extensions used for input.

Bittker, D. A.↗
NASA NTRS2000

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗
DOE OSTI2024

Probabilistic measures for biological adaptation and resilience

This paper introduces an approach to quantifying ecological resilience in biological systems, particularly focusing on noisy systems responding to episodic disturbances with sudden adaptations. Incorporating concepts from nonequilibrium statistical mechanics, we propose a measure termed “ecological resilience through adaptation,” specifically tailored to noisy, forced systems that undergo physiological adaptation in the face of stressful environmental changes. Randomness plays a key role, accounting for model uncertainty and the inherent variability in the dynamical response among components of biological systems. Our measure of resilience is rooted in the probabilistic description of states within these systems and is defined in terms of the dynamics of the ensemble average of a model-specific observable quantifying success or well-being. Our approach utilizes stochastic linear response theory to compute how the expected success of a system, originally in statistical equilibrium, dynamically changes in response to a environmental perturbation and a subsequent adaptation. Importantly, the resulting mathematical derivations allow for the estimation of resilience in terms of ensemble averages of simulated or experimental data. Finally, through a simple but clear conceptual example, we illustrate how our resilience measure can be interpreted and compared to other existing frameworks in the literature. The methodology is general but inspired by applications in plant systems, with the potential for broader application to complex biological processes.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2005

Afterbody Heating Predictions for a Mars Science Laboratory Entry Vehicle

The Mars Science Laboratory mission intends to deliver a large rover to the Martian surface within 10 km of its target site. One candidate entry vehicle aeroshell consists of a 3.75-m diameter, 70-deg sphere-cone forebody and a biconic afterbody similar to that of Viking. This paper presents computational fluid dynamics predictions of laminar afterbody heating rates for this configuration and a 2010 arrival at Mars. Computational solutions at flight conditions used an 8-species Mars gas model in chemical and thermal non-equilibrium. A grid resolution study examined the effects of mesh spacing on afterbody heating rates and resulted in grids used for heating predictions on a reference entry trajectory. Afterbody heating rate reaches its maximum value near 0.6 W/sq cm on the first windward afterbody cone at the time of peak freestream dynamic pressure. Predicted afterbody heating rates generally are below 3% of the forebody laminar nose cap heating rate throughout the design trajectory. The heating rates integrated over time provide total heat load during entry, which drives thermal protection material thickness.

Edquist, Karl T.↗
NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗
NASA NTRS1988

Development of the general interpolants method for the CYBER 200 series of supercomputers

The General Interpolants Method (GIM) is a 3-D, time-dependent, hybrid procedure for generating numerical analogs of the conservation laws. This study is directed toward the development and application of the GIM computer code for fluid dynamic research applications as implemented for the Cyber 200 series of supercomputers. An elliptic and quasi-parabolic version of the GIM code are discussed. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and an implicit finite difference scheme are also included.

Stalnaker, J. F.↗
NASA NTRS2002

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗
NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗
DOE OSTI2023

Dynamics of skyrmion contraction and expansion in a magnetic film

Contraction and expansion of skyrmions in ferromagnetic films are investigated. In centrosymmetric systems, the dynamics of a collapsing skyrmion is driven by dissipation. The collapse time has a minimum on the damping constant. In systems with broken inversion symmetry, the evolution of skyrmions toward equilibrium size is driven by the Dzyaloshinskii–Moriya interaction. As a result, expressions describing the time dependence of the skyrmion size are derived and their implications for skyrmion-based information processing are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Three-dimensional dynamics of scientific balloon systems in response to sudden gust loadings

A mathematical model was developed of the three-dimensional dynamics of a high-altitude scientific research balloon system perturbed from its equilibrium configuration by an arbitrary gust loading. The platform is modelled as a system of four coupled pendula, and the equations of motion were developed in the Lagrangian formalism assuming a small-angle approximation. Three-dimensional pendulation, torsion, and precessional motion due to Coriolis forces are considered. Aerodynamic and viscous damping effects on the pendulatory and torsional motions are included. A general model of the gust field incident upon the balloon system was developed. The digital computer simulation program is described, and a guide to its use is given.

Dorsey, D. R., Jr.↗
NASA NTRS1959

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.↗
NASA NTRSDate not supplied

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↗
DOE OSTI2020

Structure and dynamics of hydrodynamically interacting finite-size Brownian particles in a spherical cavity: Spheres and cylinders

The structure and dynamics of confined suspensions of particles of arbitrary shape are of interest in multiple disciplines from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an immersed-boundary general geometry Ewald-like method to capture lubrication and long-range hydrodynamics and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction, and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than 10%, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion toward the walls, and a sub-diffusive regime-caused by crowding-in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.

36 MATERIALS SCIENCE↗
NASA NTRS2018

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↗