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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

The tropical water and energy cycles in a cumulus ensemble model. Part 1: Equilibrium climate

A cumulus ensemble model is used to study the tropical water and energy cycles and their role in the climate system. The model includes cloud dynamics, radiative processes, and microphysics that incorporate all important production and conversion processes among water vapor and five species of hydrometeors. Radiative transfer in clouds is parameterized based on cloud contents and size distributions of each bulk hydrometeor. Several model integrations have been carried out under a variety of imposed boundary and large-scale conditions. In Part 1 of this paper, the primary focus is on the water and heat budgets of the control experiment, which is designed to simulate the convective - radiative equilibrium response of the model to an imposed vertical velocity and a fixed sea surface temperature at 28 C. The simulated atmosphere is conditionally unstable below the freezing level and close to neutral above the freezing level. The equilibrium water budget shows that the total moisture source, M(sub s), which is contributed by surface evaporation (0.24 M(sub s)) and the large-scale advection (0.76 M(sub s)), all converts to mean surface precipitation bar-P(sub s). Most of M(sub s) is transported verticaly in convective regions where much of the condensate is generated and falls to surface (0.68 bar-P(sub s)). The remaining condensate detrains at a rate of 0.48 bar-P(sub s) and constitutes 65% of the source for stratiform clouds above the melting level. The upper-level stratiform cloud dissipates into clear environment at a rate of 0.14 bar-P(sub s), which is a significant moisture source comparable to the detrained water vapor (0.15 bar-P(sub s)) to the upper troposphere from convective clouds. In the lower troposphere, stratiform clouds evaporate at a rate of 0.41 bar-P(sub s), which is a more dominant moisture source than surface evaporation (0.22 bar-P(sub s)). The precipitation falling to the surface in the stratiform region is about 0.32 bar-P(sub s). The associated latent heating in the water cycle is the dominant source in the heat budget that generates a net upward motion in convective regions, upper stratiform regions (above the freezing level), and a downward motion in the lower stratiform regions. The budgets reveal a cycle of water and energy resulted from radiation-dynamic-convection interactions that maintain equilibrium of the atmosphere.

Sui, C. H.↗
DOE OSTI2020

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Solar astrophysics - Ghettosis from, or symbiosis with, stellar and galactic astrophysics

The purpose of the paper is to show how the solar-stellar symbiotic approach has led to the modeling of a star as a concentration of matter and energy. By 'solar-stellar symbiosis' is meant the philosophy of investigation according to which one asks what change in our general understanding of stellar structure and of stellar spectroscopic diagnostics is required to satisfy both the sun and an unusual star when, for example, some feature of an unusual star is discovered. The evolution of stellar models is traced, from walled, thermodynamic-equilibrium models to de-isolated models featuring transition zones and nonlocal thermodynamic equilibrium.

Pecker, J.-C.↗
DOE OSTI2025

HD 143811 AB b: A Directly Imaged Planet Orbiting a Spectroscopic Binary in Sco-Cen

We present confirmation of HD 143811 AB b, a substellar companion to spectroscopic binary HD 143811 AB through direct imaging with the Gemini Planet Imager (GPI) and Keck NIRC2. HD 143811 AB was observed as a part of the GPI Exoplanet Survey in 2016 and 2019 and is a member of the Sco-Cen star formation region. The exoplanet is detected ∼430 mas from the host star by GPI. With two GPI epochs and one from Keck/NIRC2 in 2022, we confirm through common proper motion analysis that the object is bound to its host star. We derive an orbit with a semimajor axis of $64^{+32}_{-14}$ au and eccentricity $0.23^{+0.24}_{-0.16}$. Spectral analysis of the GPI H-band spectrum and NIRC2 L′ photometry provides additional proof that this object is a substellar companion. We compare the spectrum of HD 143811 AB b to PHOENIX stellar models and Exo-Radioactive-Convective Equilibrium Model (REM) exoplanet atmosphere models and find that Exo-REM models provide the best fits to the data. From the Exo-REM models, we derive an effective temperature of $1042^{+178}_{-132}$ K for the planet and translate the derived luminosity of the planet to a mass of 5.6 ± 1.1 M Jup assuming hot-start evolutionary models. HD 143811 AB b is the first directly imaged planet around a binary that is not on an ultrawide orbit. Future characterization of this object will shed light on the formation of planets around binary star systems.

Astronomy and AstroPhysics↗
NASA NTRS1972

Role of gas surface interactions in the reduction of OGO-6 neutral gas particle mass spectrometer data

The gas-surface interaction effects observed by the quadrupole mass spectrometer are described, and the technique developed to account for them in determining ambient neutral densities is summarized. The total ion current and the ion currents for ions with molecular weights 2, 4, 16, 28, and 32 are sampled for 1.125 sec once every 9.216 sec, for 258 sec out of a 368 sec cycle. An equation is given for the number density of any constituent in the ion source region, and source density data are discussed. The mass 28 background gas is considered to be CO rather than N2, and a CO model is developed. A quasi-equilibrium model of the atomic oxygen interactions is constructed, and a set of surface parameters is determined which provides a reasonable fit to the mass 16 and 32 source densities consistent with the predicted ambient atomic oxygen.

Hedin, A. E.↗
NASA NTRS2005

Partial Model of Insulator/Insulator Contact Charging

Two papers present a two-phase equilibrium model that partly explains insulator/ insulator contact charging. In this model, a vapor of ions within a gas is in equilibrium with a submonolayer of ions of the same species that have been adsorbed on the surface of an insulator. The surface is modeled as having localized states, each with a certain energy of adsorption for an ion. In an earlier version of the model described in the first paper, the ions do not interact with each other. Using the grand canonical ensemble, the chemical potentials of both vapor and absorbed phases are derived and equated to determine the vapor pressure. If a charge is assigned to the vapor particles (in particular, if single ionization is assumed), then the surface charge density associated with adsorbed ions can be calculated as a function of pressure. In a later version of the model presented in the second paper, the submodel of the vapor phase is extended to include electrostatic interactions between vapor ions and adsorbed ones as well as the screening effect, at a given distance from the surface, of ions closer to the surface. Theoretical values of this model closely match preliminary experimental data on the discharge of insulators as a function of pressure.

Hogue, Michael↗
NASA NTRS2014

Study of Aerothermodynamic Modeling Issues Relevant to High-Speed Sample Return Vehicles

This paper examines the application of state-of-the-art coupled ablation and radiation simulations to highspeed sample return vehicles, such as those returning from Mars or an asteroid. A defining characteristic of these entries is that the surface recession rates and temperatures are driven by nonequilibrium convective and radiative heating through a boundary layer with significant surface blowing and ablation products. Measurements relevant to validating the simulation of these phenomena are reviewed and the Stardust entry is identified as providing the best relevant measurements. A coupled ablation and radiation flowfield analysis is presented that implements a finite-rate surface chemistry model. Comparisons between this finite-rate model and a equilibrium ablation model show that, while good agreement is seen for diffusion-limited oxidation cases, the finite-rate model predicts up to 50% lower char rates than the equilibrium model at sublimation conditions. Both the equilibrium and finite rate models predict significant negative mass flux at the surface due to sublimation of atomic carbon. A sensitivity analysis to flowfield and surface chemistry rates show that, for a sample return capsule at 10, 12, and 14 km/s, the sublimation rates for C and C3 provide the largest changes to the convective flux, radiative flux, and char rate. A parametric uncertainty analysis of the radiative heating due to radiation modeling parameters indicates uncertainties ranging from 27% at 10 km/s to 36% at 14 km/s. Applying the developed coupled analysis to the Stardust entry results in temperatures within 10% of those inferred from observations, and final recession values within 20% of measurements, which improves upon the 60% over-prediction at the stagnation point obtained through an uncoupled analysis. Emission from CN Violet is shown to be over-predicted by nearly and order-of-magnitude, which is consistent with the results of previous independent analyses. Finally, the coupled analysis is applied to a 14 km/s Earth entry representative of a Mars sample return. Although the radiative heating provides a larger fraction of the total heating, the influence of ablation and radiation on the flowfield are shown to be similar to Stardust.

Johnston, Christopher O.↗
DOE OSTI2021

Polarized imaging of visible synchrotron emission from runaway electron plateaus in DIII-D

A new approach to estimate the pitch angle of highly energetic postdisruption runaway electrons is presented by imaging the polarized runaway electron synchrotron emission. The vertical to horizontal polarization ratio is found to range between 3 and 14 in image-average studies looking at repeated shots as well as the polarization in a single discharge. Numerical studies with guiding center and full orbit simulations performed with the SOFT and KORC codes are given and compared to the experiment. The pitch angle evolution analyzed using a 1D impurity diffusion model in combination with kinetic simulations is also presented. Furthermore, these studies find agreement in the magnitude of the pitch angle, but not in the temporal rate of change of pitch angle. The measured pitch angle evolution timescale of order 50 ms falls between a modeled timescale of 20 ms in an equilibrium model and >>100 ms in a nonequilibrium model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Theoretical models of Beta Cephei stars constructed with new radiative opacities

Equilibrium models for the main-sequence and early post-main-sequence phases of evolution are constructed for stars of 10.9 and 15 solar masses with two different initial hydrogen and metals abundances. Carson's (1976) radiative opacities are used in the calculations. It is found that the models which lie off the zero-age main sequence have cooler effective temperatures than corresponding models based on other opacities and cross the observed Beta Cephei strip in the H-R diagram only during the main phase of core hydrogen burning. The stability of radial and nonradial pulsations is studied by applying linear nonadiabatic perturbations to the models. The results show that all the modes examined are stable and that the margin of stability decreases with advancing evolution. It is suggested that the closeness of the periods of the two lowest nonradial quadrupole p-modes may be responsible for the 'beat' phenomenon observed in many Beta Cephei stars, if these modes are somehow excited. Possible instability mechanisms are evaluated, particularly the ultimate ionization of CNO elements.

Stothers, R.↗
NASA NTRS1990

Plasma diffusion at the magnetopause? The case of lower hybrid drift waves

The diffusion expected from the quasilinear theory of the lower hybrid drift instability at the Earth's magnetopause is recalculated. The resulting diffusion coefficient is in principle just marginally large enough to explain the thickness of the boundary layer under quiet conditions, based on observational upper limits for the wave intensities. Thus, one possible model for the boundary layer could involve equilibrium between the diffusion arising from lower hybrid waves and various low processes. However, some recent data and simulations seems to indicate that the magnetopause is not consistent with such a soft diffusive equilibrium model. Furthermore, investigation of the nonlinear equations for the lower hybrid waves for magnetopause parameters indicates that the quasilinear state may never arise because coalescence to large wavelengths, followed by collapse once a critical wavelengths is reached, occur on a time scale faster than the quasilinear diffusion. In this case, an inhomogeneous boundary layer is to be expected. More simulations are required over longer time periods to explore whether this nonlinear evolution really takes place at the magnetopause.

Treumann, R. A.↗
NASA NTRS1992

Pulsatile instability in rapid directional solidification: Strongly-nonlinear analysis

In models of rapid directional solidification, non-equilibrium interfacial conditions are employed. As a result, there is an oscillatory mode of instability, as well as the steady cellular mode, found in the equilibrium model of Mullins and Sekerka. When the temperature field is decoupled from the problem, the preferred wave number for the oscillatory mode is zero, and the interface pulsates in time while remaining spatially uniform. Results from multiple scale analyses in the two limiting cases of the parameters are reported. In these limits, it is found that the instability is a bifurcation to relaxation oscillations; these nonlinear oscillations may be related to the observed microstructure that results from rapid solidification processes such as laser surface remelting.

Braun, Richard J.↗
NASA NTRS1995

Experience with turbulence interaction and turbulence-chemistry models at Fluent Inc.

This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.

Choudhury, D.↗
DOE OSTI2025

Aerosol Thermodynamics (Aerosol Liquid Water Concentration and Aerosol pH) at S2 during CoURAGE

These data report aerosol liquid water content (ALWC) and aerosol pH for the S2 (Mt. Airy) site during CoURAGE for April - June 2025. The ISORROPIA-II aerosol thermodynamic equilibrium model was used to compute ALWC and pH. The model was run in forward mode, with solids formation disabled (metastable mode), according to Pye et al. (Atmospheric Chemistry and Physics, 2020). Information about the model can be found at https://www.epfl.ch/labs/lapi/models-and-software/isorropia/, and in Fountoukis and Nenes (Atmospheric Chemistry and Physics, 2007). Inputs to the model were: (1) measured Temperature and relative humidity, both from the DoE ARM “metwxt” data products at S2, (2) aerosol chemical composition from the ARM ACSM, and (3) gas-phase ammonia concentrations. For ISORROPIA-II, ambient RH values above 0.995 were input to the model as 0.995; therefore, caution should be exercised interpreting ALWC when RH was above 0.995. Aerosol composition or NH3(g) concentrations below the LOD were input into the model as 0.5*LOD. Model runs were only conducted for data points at 30-min resolution with all three inputs (meteorology, ACSM, and NH3). If any of the three model inputs were missing, aerosol thermodynamic outputs are flagged as -999. Similarly, aerosol thermodynamic outputs are flagged as -999 if measured NH3(g) and aerosol ammonium (NH4+) were simultaneously below the respective LODs. The ALWC includes liquid water from inorganic aerosol components, calculated directly by ISORROPIA-II, as well as liquid water from organics, estimated from the ACSM organics assuming a kappa value of 0.12 according to Guo et al. (Atmospheric Chemistry and Physics, 2015).

Aerosol Liquid Water Content↗
NASA NTRS1978

Five micron limb-darkening and the structure of the Jovian atmosphere

A combined radiative, dynamical, and chemical model of the Jovian atmosphere is presented which reproduces the observed IR limb-darkening curves for several regions on the disk of Jupiter. Two approaches are used to perform the inversion required for determination of the source function, and a power-series-expansion model is examined, along with a convective-equilibrium model, cloud and intermediate-zone models, and a thin-shell model. The statistical methods used in the analysis are outlined, and model fits to the observations are discussed. The chemistry of phase transitions in the hot spots of Jupiter's equatorial belt is briefly investigated.

Newman, W. I.↗
NASA NTRS1987

A comparison of turbulence closure models for transonic flows about airfoils

Navier-Stokes transonic airfoil calculations using three different turbulence closure models are presented and compared with available experimental data. Two of the models, Cebeci-Smith and Baldwin-Lomax, are equilibrium models in the sense that the turbulent shear stress is assumed to depend only on the local properties of the mean flow. The third model, Johnson and King, is a nonequilibrium model in which an ordinary differential equation is employed to account for the streamwise development of the maximum Reynolds shear stress. An eddy viscosity distribution across the boundary layer is assumed which is functionally dependent on this stress. For attached flows with little inviscid-viscous interaction, little difference in the results with the different models is noted. With stronger interactions and with separation occurring on the airfoils, the nonequilibrium model is shown to perform better.

King, Lyndell S.↗
NASA NTRS1991

Axisymmetric model of the ionized gas in the Orion Nebula

New ionization and thermal equilibrium models for the ionized gas in the Orion Nebula with an axisymmetric two-dimensional 'blister' geometry/density distribution are presented. The HII region is represented more realistically than in previous models, while the physical detail of the microphysics and radiative transfer of the earlier spherical modeling is maintained. The predicted surface brightnesses are compared with observations for a large set of lines at different positions to determine the best-fitting physical parameters. The model explains the strong singly ionized line emission along the lines of sight near the Trapezium.

Rubin, R. H.↗
NASA NTRS2022

Economic incentives modify agricultural impacts of nuclear war

A nuclear war using less than 1% of the current global nuclear arsenal, which would inject 5 Tg of soot into the stratosphere, could produce climate change unprecedented in recorded human history and significant impacts on agricultural productivity and the economy. These effects would be most severe for the first five years after the nuclear war and may last for more than a decade. This paper calculates how food availability would change by employing the Environmental Impact and Sustainability Applied General Equilibrium model. Under a robust world trading system, global food availability would drop by a few percentage points. If the war would destabilize trade, it would magnify by several times the negative ramifications of land productivity shocks on food availability. If exporting countries redirect production to domestic consumption at the expense of importing countries, it would lead to the destabilization of international trade. The analysis suggests that economic models aiming to inform policymakers require both economic behavior analysis and biophysical drivers. Policy lessons derived from a crop model can be significantly nuanced when coupled with economic feedback derived from economic models. Through the impact on yield, farmers could shift production among crops and reallocate land use to maximize profits, showing the importance of general equilibrium effects such as product and input substitution and international trade. Although the global impact on corn and soybean production would be significant when just considering crop production, it could be considerably smaller under the economic model. However, this would be at the expense of other sectors, including livestock. In addition, the costs borne from disruptions to climate would vary significantly across regions, with significant adverse effects in high latitude regions. The severity of the shocks in the high-latitude areas would marginalize the farmers' product and input substitution ability.

Nuclear war↗
DOE OSTI2025

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium E h values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9.3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-1 and Table ES-2, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗