Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Neural-Network-Enhanced COTSIM: Advancing Predictive Capabilities for Fast DIII-D Simulations

Sustaining fusion reactions in tokamaks requires heating plasma to thermonuclear temperatures while maintaining confinement and stability. Neutral beam injection (NBI) provides heating, current drive, torque, and fueling, while electron cyclotron (EC) waves are widely used for heating and current drive; together, these actuators shape the plasma current, temperature, and density profiles. The control-oriented tokamak simulator (COTSIM), a predictive, control-oriented code, has been enhanced with neural-network surrogates for transport and sources. Turbulent transport is predicted by MMMnet—a neural-network version of the updated multimode model (MMM 9.0.10)—with significantly reduced computation time relative to MMM; neoclassical transport follows the Chang–Hinton model. NUBEAMnet, a surrogate of the Monte Carlo NUBEAM module, predicts beam-driven heating, current, and torque. EC heating and current drive use a control-oriented, empirically scaled source model; plasma resistivity follows the Spitzer formulation; bootstrap current uses the Sauter model. Equilibrium is computed using both prescribed and fixed-boundary solvers (FBSs), and the pedestal structure is modeled with an empirical pedestal model. For a representative DIII-D discharge, COTSIM predicts electron and ion temperature and safety-factor profiles in close agreement with TRANSP predictive and interpretive simulations while extending predictions through the pedestal region to the plasma edge (versus 80% of the minor radius in TRANSP). Furthermore, the equivalent COTSIM simulation runs in under 3 min compared to about 2 h for TRANSP, enabling rapid scenario planning, optimization of tokamak operation, and between-pulse control design.

Control-oriented tokamak simulator (COTSIM)↗
NASA NTRS1999

Model Atmospheres for Novae in Outburst: Summary of Research

This paper presents a final report and summary of research on Model Atmospheres for Novae in Outburst. Some of the topics include: 1) Detailed NLTE (non-local thermodynamic equilibrium) Model Atmospheres for Novae during Outburst: II. Modeling optical and ultraviolet observations of Nova LMC 1988 #1; 2) A Non-LTE Line-Blanketed Stellar Atmosphere Model of the Early B Giant epsilon CMa; 3) Spectroscopy of Low Metallicity Stellar atmospheres; 4) Infrared Colors at the Stellar/Substellar Boundary; 5) On the abundance of Lithium in T CrB; 6) Numerical Solution of the Expanding Stellar Atmosphere Problem; and 7) The NextGen Model Atmosphere grid for 3000 less than or equal to T (sub eff) less than or equal to 10000K.

Hauschildt, Peter H.↗
DOE OSTI2024

Reductive Dynamic and Static Excited State Quenching of a Homoleptic Ruthenium Complex Bearing Aldehyde Groups

A new homoleptic Ru polypyridyl complex bearing two aldehyde groups on each bipyridine ligand, [Ru(dab) 3 ](PF 6 ) 2 , where dab is 4,4′-dicarbaldehyde-2,2′-bipyridine, was synthesized, characterized, and utilized for iodide photo-oxidation studies. In acetonitrile (CH 3 CN) solution, the complex displayed an intense metal-to-ligand charge transfer (MLCT) absorbance maximum at 475 nm (ε = 22,000 M –1 cm –1 ) and an infrared (IR) band at 1712 cm –1 assigned to the pendent aldehyde groups. Visible light excitation in air-saturated solution resulted in room temperature photoluminescence (PL) with a maximum at 675 nm, a quantum yield, ϕ PL = 0.048, and an excited state lifetime, τ ο = 440 ns, from which radiative and nonradiative relaxation rate constants were extracted, k r = 9.1 × 10 4 s –1 and knr = 1.8 × 10 6 s –1 . Pulsed visible light excitation yielded transient UV–vis and IR absorption spectra consistent with an MLCT excited state; relaxation occurred with the maintenance of two isosbestic points in the visible region, and a lifetime that agreed with that measured by time-resolved PL. Cyclic voltammetry studies in a CH 3 CN solution with 0.1 M TBAPF 6 electrolyte revealed a quasi-reversible oxidation, E°(Ru III/II ) = +1.25 V vs. Fc +/0 , and three sequential one-electron reductions at −1.10, −1.25, and −1.54 V vs. Fc +/0 . Here, an excited state reduction potential of E°(Ru *2+/+ ) = +0.89 V vs. Fc +/0 was estimated with the Rehm–Weller expression. Titration of tetrabutylammonium iodide, TBAI, into a CD 3 CN solution of [Ru(dab) 3 ](PF 6 ) 2 resulted in significant shifts in the aldehyde H atom and 3,3′-biypridyl resonances that were analyzed with a 1:1 equilibrium model, from which K eq = 460 M –1 was extracted, increasing to 5800 M –1 when the solvent was changed to acetone-d 6 . Iodide titrations resulted in a significant quenching of the [Ru(dab) 3 ] *2+ lifetime and quantum yield in both CH 3 CN and acetone solvents. In CH 3 CN, the quenching was mainly dynamic and well described by the Stern–Volmer model, from which a quenching rate constant, k q , of 4.5 × 10 10 M –1 s –1 and an equilibrium constant, K eq , of 8.3 × 10 3 M –1 were obtained. In acetone, the static quenching pathway by iodide was greatly enhanced, with a K eq of 1.2 × 10 4 M –1 and a higher k q of 9.2 × 10 10 M –1 s –1 .

Aldehydes↗
DOE OSTI2020

Pervasive Warming Bias in CMIP6 Tropospheric Layers

Abstract The tendency of climate models to overstate warming in the tropical troposphere has long been noted. Here we examine individual runs from 38 newly released Coupled Model Intercomparison Project Version 6 (CMIP6) models and show that the warm bias is now observable globally as well. We compare CMIP6 runs against observational series drawn from satellites, weather balloons, and reanalysis products. We focus on the 1979–2014 interval, the maximum span for which all observational products are available and for which models were run using historically observed forcings. For lower‐troposphere and midtroposphere layers both globally and in the tropics, all 38 models overpredict warming in every target observational analog, in most cases significantly so, and the average differences between models and observations are statistically significant. We present evidence that consistency with observed warming would require lower model Equilibrium Climate Sensitivity (ECS) values.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2006

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters (or closwre models) for an initial thermai non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two CFD codes currently being used at Glenn Research Center (GRC) for Stirling engine modeling are Fluent and CFD-ACE. The porous-media models available in each of these codes are equilibrium models, which assmne that the solid matrix and the fluid are in thermal equilibrium at each spatial location within the porous medium. This is believed to be a poor assumption for the oscillating-flow environment within Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, we non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location end time during the cycle. A NASA regenerator research grant has been providing experimental and computational results to support definition of various empirical coefficients needed in defining a noa-equilibrium, macroscopic, porous-media model (i.e., to define "closure" relations). The grant effort is being led by Cleveland State University, with subcontractor assistance from the University of Minnesota, Gedeon Associates, and Sunpower, Inc. Friction-factor and heat-transfer correlations based on data taken with the NASAlSunpower oscillating-flow test rig also provide experimentally based correlations that are useful in defining parameters for the porous-media model; these correlations are documented in Gedeon Associates' Sage Stirling-Code Manuals. These sources of experimentally based information were used to define the following terms and parameters needed in the non-equilibrium porous-media model: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity (including themal dispersion and estimate of tortuosity effects}, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity (including the effect of tortuosity) was also estimated. Determination of the porous-media model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Convertor (TDC), which uses a random-fiber regenerator matrix. The non-equilibrium porous-media model presented is considered to be an initial, or "draft," model for possible incorporation in commercial CFD codes, with the expectation that the empirical parameters will likely need to be updated once resulting Stirling CFD model regenerator and engine results have been analyzed. The emphasis of the paper is on use of available data to define empirical parameters (and closure models) needed in a thermal non-equilibrium porous-media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates. However, it is anticipated that a thermal non-equilibrium model such as that presented here, when iacorporated in the CFD codes, will improve our ability to accurately model Stirling regenerators with CFD relative to current thermal-equilibrium porous-media models.

Tew, Roy↗
DOE OSTI2020

Truck platooning in the U.S. national road network: A system-level modeling approach

Truck platooning enables a group of trucks to move close together, which helps reduce truck fuel use and increase effective road capacity. In this paper, a system-level equilibrium model is developed to characterize spontaneous truck platooning with coexistence of non-platooning vehicles in a network, by explicitly accounting for the interlocking relationship among platoon formation time, truck fuel saving, and increase in effective road capacity. To equilibrate the relationships, an algorithm is proposed which involves a diagonalization approach and a bush based algorithm to solve decomposed subproblems. The condition of proportionality is imposed to obtain unique traffic flows for each class of vehicles on road links. In addition, a spatially constrained multivariate clustering technique is employed to construct origin/destination zones that are smaller than the coarse Freight Analysis Framework (FAF) zones, while maintaining reasonable computational burden for network traffic assignment. Model implementation in the U.S. shows that platooning could lead to 7.9% fuel saving among platoonable trucks in 2025 and a comparable increase in effective capacity of platoonable road links, which would account for 60% of rural interstate roads. The fuel saving and road capacity improvement translate into an annual cost reduction of $\$$868 million for the U.S. intercity trucking sector and reduced road infrastructure investment needs worth $\$$4.8 billion. Extensive sensitivity analysis further reveals that fuel saving of platoonable trucks increases with platoon size but decreases with inter-truck distance in a platoon. Fuel saving potential suggests that priority should be given to rural rather than urban roads in deploying platooning technologies. As we expected, greater market penetration of platooning technologies means higher fuel saving and greater increase in effective road capacity.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1986

Diagnostic analysis of two-dimensional monthly average ozone balance with Chapman chemistry

Chapman chemistry has been used in a two-dimensional model to simulate ozone balance phenomenology. The similarity between regions of ozone production and loss calculated using Chapman chemistry and those computed using LIMS and SAMS data with a photochemical equilibrium model indicate that such simplified chemistry is useful in studying gross features in stratospheric ozone balance. Net ozone production or loss rates are brought about by departures from the photochemical equilibrium (PCE) condition. If transport drives ozone above its PCE condition, then photochemical loss dominates production. If transport drives ozone below its PCE condition, then photochemical production dominates loss. Gross features of ozone loss/production (L/P) inferred for the real atmosphere from data are also simulated using only eddy diffusion. This indicates that one must be careful in assigning a transport scheme for a two-dimensional model that mimics only behavior of the observed ozone L/P.

Stolarski, Richard S.↗
NASA NTRS2017

On the Development of a New Nonequilibrium Chemistry Model for Mars Entry

This paper represents a summary of results to date of an on-going effort at NASA Ames Research Center to develop a physics-based non-equilibrium model for hypersonic entry into the Martian atmosphere. Our approach is to first compute potential energy surfaces based on accurate solutions of the electronic Schroedinger equation and then use quasiclassical trajectory calculations to obtain reaction cross sections and rate coefficients based on these potentials. We have presented new rate coefficients for N2 dissociation and CO dissociation and exchange reactions. These results illustrate shortcomings with some of the rate coefficients in Parks original T-Tv model for Mars entries and with some of the 30-45 year old shock tube data. We observe that the shock tube experiments of CO + O dissociation did not adequately account for the exchange reaction that leads to formation of C + O2. This reaction is actually the primary channel for CO removal in the shock layer at temperatures below 10,000 K, because the reaction enthalpy for exchange is considerably lower than the comparable value for dissociation.

quasiclassical trajectory calculation↗
NASA NTRS1993

The embedded young stars in the Taurus-Auriga molecular cloud. II - Models for scattered light images

We describe NIR imaging observations of embedded young stars in the Taurus-Auriga molecular cloud. We find a large range in J-K and H-K colors for these class I sources. The bluest objects have colors similar to the reddest T Tauri stars in the cloud; redder objects lie slightly above the reddening line for standard ISM dust and have apparent K extinctions of up to 5 mag. Most of these sources also show extended NIR emission on scales of 10-20 arcsec which corresponds to linear sizes of 1500-3000 AU. The NIR colors and nebular morphologies for this sample and the magnitude of linear polarization in several sources suggest scattered light produces most of the NIR emission in these objects. We present modeling results that suggest mass infall rates that agree with predictions for cold clouds and are generally consistent with rates estimated from radiative equilibrium models. For reasonable dust grain parameters, the range of colors and extinctions require flattened density distributions with polar cavities evacuated by bipolar outflows. These results support the idea that infall and outflow occur simultaneously in deeply embedded bipolar outflow sources. The data also indicate fairly large centrifugal radii and large inclinations to the rotational axis for a typical source.

Kenyon, Scott J.↗
NASA NTRS1994

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗
DOE OSTI2020

Soft X-ray and EUV emission spectra of beryllium plasma produced by neodymium-glass laser radiation with broad frequency and angular spectra

We present the results of an experimental study of soft X-ray (SXR) and extreme ultraviolet emission spectra of the plasma produced by exposing a plane solid beryllium target to laser radiation with broad frequency and angular spectra. SXR lines up to 1s – 9p of Be IV as well as the plasma continuum are recorded for a laser focal-spot intensity of 5.3 × 10{sup 13} W cm{sup −2}. To model the SXR beryllium plasma spectra, simulations are carried out using the INDHAUS programme and the FLYCHK code in the framework of local thermodynamic equilibrium model, which agree nicely with experimentally obtained results. (laser plasma)

36 MATERIALS SCIENCE↗
NASA NTRS1990

Model of Jovian F region ionosphere (Jovian ionosphere model in offset dipole approximation)

The geomagnetic control of the Earth's atmosphere is well understood. In the F-region and the topside ionosphere, non-electrical forces transport plasma along the magnetic field lines only. In consequence, the worldwide distribution of ionization is strongly dependent on the dip angle. For example, the equatorial anomaly is roughly symmetrical about the dipole equator rather than the geographic. The same appears to be the case in the Jovian ionosphere (Mahajan, 1981). The influence of the magnetic field of Jupiter on its ionization pattern is one of several outstanding topics which need to be studied. Tan (1986) investigated the formation of the equatorial anomaly in the Jovian ionosphere under a centered dipole model. Tan (1988) further studied the effect of the tilt of the Jovian dipole. The results were in broad agreement with those of a diffusive equilibrium model (Tan and Wu, 1981). An off-centered dipole model is constructed and its effects on the ionization pattern are investigated.

Tan, A.↗
NASA NTRS2006

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters for an initial thermal non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two codes currently used at Glenn Research Center for Stirling modeling are Fluent and CFD-ACE. The codes porous-media models are equilibrium models, which assume solid matrix and fluid are in thermal equilibrium. This is believed to be a poor assumption for Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, use non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle. Experimentally based information was used to define: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity was also estimated. Determination of model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Converter (TDC), which uses a random-fiber regenerator matrix. Emphasis is on use of available data to define empirical parameters needed in a thermal non-equilibrium porous media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates.

Tew, Roy C.↗
DOE OSTI2020

Collinear dual-pulse laser optical breakdown and energy deposition

A non-equilibrium model for laser-generated plasmas is used to represent collinear dual nano-second-pulse interactions. In the specific case considered, the breakdown is initiated with an ultraviolet (UV) laser pulse at 266 nm, which pre-ionizes the gas, and it is followed by a near-infrared (NIR) pulse at 1064 nm, which deposits significant energy into the ionized mixture. The model is validated against corresponding experiments, and simulation results are interrogated to understand key features of the plasma-kernel dynamics and the post-discharge hydrodynamics. The hydrodynamics of the non-equilibrium plasma is governed by the two-temperature Navier–Stokes equations accounting for both multiphoton ionization and inverse bremsstrahlung. The interaction between the laser beam and the plasma is modeled based on the radiative transfer equation. The temporal and spatial offsets of the two pulses can generate ionization kernels with different topology and dynamics. It is shown that the UV pre-ionization pulse can tailor the plasma region, leading to a larger ionized volume than would occur for a single-pulse breakdown, and it increases the efficiency of the energy deposition for the following NIR discharge. Vorticity in the early post-discharge phase is generated via baroclinic torque from the misalignment of the radial gradient of density (sudden gas expansion) and the strong pressure gradient (initiated by energy deposition from the laser beam). It is predicted that the collinear dual pulse can be designed to specify the plasma kernel size, shape, and maximum temperature in the breakdown phase, and the initial post-breakdown vorticity and plasma-core decay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Linear analysis of poststall gyrations

The poststall gyrations of a high-performance aircraft are investigated by approximating the equations of motion with a linear mathematical model. Equilibrium poststall spin conditions are determined, which are not limited to small angles or spin rates and are used as reference conditions about which the linear model is developed. As a result, complete coupling between lateral and longitudinal motions is retained. Three aerodynamic models are evaluated using this linear analysis. These models include pure rotary derivatives obtained from VPI&SU stability tunnel, forced-oscillation data obtained from NASA tests, and a combination of the two which separates the unsteady and pure rotary derivatives. Results indicate that this last model gives better predictions of aircraft motion during poststall gyrations at high angles of attack and large sideslip angles. Comparisons are made with full-scale flight test results.

Hreha, M. A.↗
NASA NTRS1994

A model for studying the composition and chemical effects of stratospheric aerosols

We developed polynomial expressions for the temperature dependence of the mean binary and water activity coefficients for H2SO4 and HNO3 solutions. These activities were used in an equilibrium model to predict the composition of stratospheric aerosols under a wide range of environmental conditions. For typical concentrations of H2O, H2SO4, HNO3, HCl, HBr, HF, and HOCl in the lower stratosphere, the aerosol composition is estimated as a function of the local temperature and the ambient relative humidity. For temperatures below 200 K, our results indicate that (1) HNO3 contributes a significant mass fraction to stratospheric aerosols, and (2) HCl solubility is considerably affected by HNO3 dissolution into sulfate aerosols. We also show that, in volcanically disturbed periods, changes in stratospheric aerosol composition can significantly alter the microphysics that leads to the formation of polar stratospheric clouds. The effects caused by HNO3 dissolution on the physical and chemical properties of stratospheric aerosols are discussed.

Tabazadeh, Azadeh↗
DOE OSTI2021

plutus: An R package to calculate electricity investments and stranded assets from the Global Change Analysis Model (GCAM)

plutus is an R package that post-processes outputs from the Global Change Analysis Model (GCAM) to calculate the power sector capital investment costs and stranded asset costs associated with GCAM projections of electricity generation by technology. GCAM is a market equilibrium model used to examine the dynamics of the coupled human-Earth system and the response of this system to global change, including socioeconomics, technology, climate, and policy. GCAM tracks electricity generation by technology and vintage over 32 geopolitical regions throughout the lifetime of each generating technology. plutus extends GCAM functionality by (1) estimating the foregone economic value of prematurely retired power plants as a result of economic-induced retirements compared to scheduled lifetimes; and (2) estimating the new installations and capital investments driven by future changes in economic, energy, agriculture, and land-use systems in GCAM.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Figuring Out Gas & Galaxies in Enzo (FOGGIE). VI. The Circumgalactic Medium of L * Galaxies Is Supported in an Emergent, Nonhydrostatic Equilibrium

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium, with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly resolved L* galaxies from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM (Δx ≲ 1 kpc h –1 and M ≃ 200M ⊙ ). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of ≲5 kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond ~0.25R 200 , or ~50 kpc at z = 0. Within ~0.25R 200 , turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the nonequilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies, even when averaging out small-scale variations.

79 ASTRONOMY AND ASTROPHYSICS↗