SEARCH · Search NASA

Results for “Equilibria”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2025

MCCCS-MN

The MCCCS‒MN (Monte Carlo for Complex Chemical Systems‒Minnesota) software is developed by the Siepmann research group at the University of Minnesota. MCCCS‒MN allows for the simulation of multi-component molecular systems in the canonical, isobaric-isothermal (including constant stress for solids), grand-canonical, semi-grand, and Gibbs (NVT, NPT, and more than two simulation boxes) ensembles. It uses the configurational-bias Monte Carlo method to efficiently sample phase space for linear, branched and cyclic chain molecules, the adiabatic nuclear and electronic sampling Monte Carlo method to treat many-body polarization effects, and the aggregation-volume-bias Monte Carlo algorithm to efficiently sample the spatial distribution of associating molecules. MCCCS-MN employs a molecular representation of the system where force fields contain bonded and non-bonded terms. Funding for the development of MCCCS-MN through grants from the National Science Foundation (simulation of fluid phase equilibria and chromatography) and the Department of Energy (simulation of adsorption equilibria) is gratefully acknowledged.

Siepmann, J.Ilja [University of Minnesota - Twin C↗

DOE OSTI2024

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

IDAES-PSE 2.5.0 Release

AS↗

DOE OSTI2024

IDAES-PSE 2.6.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.

AS↗

DOE OSTI2020

Dissolution of Fluoride Salts in Hanford Tank Waste

The Direct Feed High-Level Waste (DFHLW) strategy seeks to bypass the Hanford Waste Treatment and Immobilization Plant Pretreatment Facility while retaining some processing functions to maximize waste feed loading and minimize high-level waste (HLW) waste volume. The DFHLW flowsheet needs leaching, washing, and solids concentration operations either in new or existing tanks. The effectiveness and efficiency of sludge washing has a substantial impact on DST space, mission duration, and the evaporation and low-activity waste (LAW) treatment operations required by these large wash-water additions. One target species requiring washing is fluoride. The HLW glass composition limits for fluorine drive operations to dissolve fluoride-bearing salts into the LAW fraction and thereby maximize waste loading in HLW glass. The fluoride in many high-level wastes at Hanford is predominantly in the form of fluoride-salt precipitates: villiaumite (NaF), kogarkoite (Na 3 FSO 4 ), and natrophosphate (Na 7 F(PO 4 ) 2 ·19H 2 O). Fluoride produces melter off-gas that creates corrosion risk in the off-gas system piping, while the sulfate and phosphate in the fluoride double salts kogarkoite and natrophosphate can be detrimental to glass waste loading. Fluoride salts are sparingly soluble, with solubilities ranging from approximately 40 to 130 kg per kL of pure water, and the dissolution kinetics of the three fluoride salts are not well known. Unexpected delays in a tank dissolution process could be encountered as a result of the lack of information about dissolution rates. In addition, if the double salts show transient non-stoichiometric dissolution of fluoride versus phosphate or sulfate, unexpectedly high concentrations of one of these other constituents could be produced. Washington River Protection Solutions authorized Pacific Northwest National Laboratory to collect the available data for fluoride salt dissolution rate, provide a scoping estimate of dissolution time if possible, and identify gaps in the understanding and predictive capability for estimating dissolution time. Open literature and Hanford reports were reviewed to document, understand, and (where possible) evaluate limitations on fluoride salt equilibria and dissolution kinetics, including both mass transport and surface reaction rate. Scoping estimates of dissolution time were made for mass-transfer-controlled dissolution of spherical particles of fluoride salts suspended in liquid. This is not the only potential governing mechanism; dissolution could be substantially slower if the surface reaction rate (the rate of release of ions from the surface) is the controlling mechanism. When the minimum amount of water for complete dissolution is used and the slip velocity between the liquid and suspended particles is less than or equal to the terminal settling velocity, the estimated mass-transfer rates allow 0.1 mm particles of fluoride salt to dissolve in minutes at 25 °C in water containing no other dissolved salts. Much larger solids, such as the 6-mm chunks that have been seen in heels, could take a few hours to more than a week to dissolve. The actual dissolution times will depend strongly on the actual slip velocity, the extent of particle suspension, constraint by surface reaction rates, the ratio of solvent to solid, and the presence of common ions that shift the solubility equilibria to restrict dissolution of fluoride salts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Physics-Based Real-time Analysis and Control to Achieve Transient-Free Operations for the ITER Era (Final Report)

This is the final report for the early career award for Prof. Egemen Kolemen that aimed to develop the physics and technical basis of real-time transient – i.e. disruptions and Edge Localized Modes (ELMs) – avoidance and feedback control at the DIII-D nuclear fusion facility for the ITER era. During this project, we reconstructed high quality real-time DIII-D equilibria using kinetic information from the motional Stark effect (MSE), Thomson Scattering (TS), and Charge Exchange Recombination (CER) spectroscopy. These equilibria were analyzed in real-time to determine the low-n ideal MHD stability and the pedestal stability thresholds, and to obtain plasma response models. These calculations were used to understand which physical process evolves to lead to the onset of disruption in ITER-relevant scenarios at DIII-D, and how disruption can be predicted using stability analysis as opposed to the current methods of correlation of various signals, which are descriptive but lack predictive power. In the disruption studies, we focused on tearing modes, since these constitute the vast majority of disruptions in this regime. With this knowledge, real-time disruption avoidance algorithms were developed and tested. The enhanced situation awareness and physics insight were used in advanced real-time feedback control of pedestal and core profiles, which allowed the pursuit and sustainment of transient-free high-confinement regimes that cannot be achieved via feed-forward design based on in-between shot analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1974

Ultrapure materials - Containerless evaporation and the roles of diffusion and Marangoni convection

Contamination from containers is a major problem in preparing ultrapure refractory materials. Space with its zero gravity and its high vacuum offers an opportunity for containerless purification of these materials. The evaporation of impurities from a melt will involve many complex chemical equilibria. Thermodynamic calculations have been modified to describe these equilibria when impurities in the melt evaporate into vacuum. The contributions of diffusion and Marangoni convection to mass transfer rates in the bulk liquid have been estimated. Calculations for the evaporative purification of molten alumina are given.

Dragoo, A. L.↗

NASA NTRS1975

Origin of titaniferous lunar basalts

Delineation of low pressure phase equilibria in the composition space relevant to titaniferous lunar basalts demonstrates a significant degree of control by those equilibria on the compositions of the basalts. The existence of two distinct chemical groups of basalts (high and low K) which cannot be related one to the other by fractional crystallization at any pressure, suggests that melting is responsible for the two groups. Consideration of the pressure shift required to produce the differences between groups constrains magma segregation to have occurred in the outer 150 km of the moon. It is difficult to relate low-Ti and high-Ti basalts to the same source region. The preferred source region of high-Ti basalts, based on phase equilibrium considerations, is a late ilmenite-rich cumulate produced from the residual liquid of the primordial differentiation of the outer portions of the moon. This ilmenite-rich layer is sandwiched between the lunar feldspathic crust and a complementary mafic cumulate.

Walker, D.↗

NASA NTRS1980

The heating of a thermally conducting stratified medium. I - Self-gravitating gas threaded by magnetic fields

Two theoretical problems are presented to illustrate the equilibrium structure of a gravitating gas threaded by magnetic fields. The gas is heated everywhere at a rate proportional to the local density and cooled by thermal conduction channeled along the magnetic field. Emphasis is placed on the influence of anisotropic thermal conduction on the distribution of the gas. In the first problem, the magnetic field has straight field lines that are nearly parallel. The simple field line geometry allows us to calculate, in closed form, the equilibria configurations of a gas slab which is cooled nonuniformly by thermal conduction along field lines which leave the gas slab at varying angles. The problem involves the gas having a free boundary, which is one of the unknowns to be determined. A discussion is given of the set of all possible equilibria so constructed. In the second problem there is no electric current flowing in the gas but there is a potential magnetic field having curved field lines. Exact solutions of equilibriums are presented to illustrate the heating of stellar and galactic coronas.

Lerche, I.↗

NASA NTRS1981

One-dimensional models of quasi-neutral parallel electric fields

Parallel electric fields can exist in the magnetic mirror geometry of auroral field lines if they conform to the quasineutral equilibrium solutions. Results on quasi-neutral equilibria and on double layer discontinuities were reviewed and the effects on such equilibria due to non-unique solutions, potential barriers and field aligned current flows using as inputs monoenergetic isotropic distribution functions were examined.

Stern, D. P.↗

NASA NTRS1982

Preliminary modeling of high pressure partial melting - Implications for early lunar differentiation

A quantitative model is presented for equilibrium crystallization of lunar magmas at high pressure, which, viewed in reverse, is equivalent to equilibrium partial melting. Data are extrapolated from simple systems, mare basalt systems, and from experiments on terrestrial compositions. Liquidus phase equilibria in a portion of the CaO-MgO-Al2O3-SiO2 system serve as a guide for predicting phase equilibria in the natural system. High pressure partial melting sequences are calculated for several proposed whole moon compositions at nominal pressure of 10 and 20 kb. Results show that compositions with chondritic Ca/Al have quite different melting paths from those with subchondritic Ca/Al, and as a result the types of magmas produced are quite different. If the moon possesses a chondritic Ca/Al ratio, models which rely on small degrees of partial melting to produce the early lunar magnetism are suspect.

Longhi, J.↗

NASA NTRS1982

Chemistry of uranium in aluminophosphate glasses

The U(VI)-U(V)-U(IV) redox equilibria are investigated in two sodium aluminophosphate base compositions at a variety of melt temperatures, imposed oxygen fugacities, and uranium contents. Results show that the higher redox states of uranium are quite soluble in the phosphate glasses, although U(IV) readily precipitates from the melts as UO2. In addition, comparisons of the uranium redox equilibria established in phosphate melts versus those in silicate melts shows that the coordination sites of the individual uranium species are generally the same in both solvent systems although they differ in detail.

Schreiber, H. D.↗

NASA NTRS1983

Relativistic thermal plasmas - Effects of magnetic fields

Processes and equilibria in finite, relativistic, thermal plasmas are investigated, taking into account electron-positron creation and annihilation, photon production by internal processes, and photon production by a magnetic field. Inclusion of the latter extends previous work on such plasmas. The basic relations for thermal, Comptonized synchrotron emission are analyzed, including emission and absorption without Comptonization, Comptonized thermal synchrotron emission, and the Comptonized synchrotron and bremsstrahlung luminosities. Pair equilibria are calculated, including approximations and dimensionless parameters, the pair balance equation, maximum temperatures and field strengths, and individual models and cooling curves.

Araki, S.↗

NASA NTRS1984

Long-term dynamical stability and transition in a spectral model

The relationship between blocking and multiple equilibria was investigated in a two layer model of the Northern Hemisphere January heating and boundary conditions are imposed. The truncation varies over global wavenumbers 4-8 and the coefficients of horizontal diffusion of heat and potential vorticity are varied about realistic values. It is shown that while severely truncated models can exhibit multiple quasi equilibria under realistic domain, forcing, and boundary conditions, higher resolution is required to allow realistic transitions to be mediated by strong baroclinic wave activity. It is indicated that for extended range forecasting blocking transitions are not likely to be forecast well unless the associated synoptic events can be forecasted statistically.

Wobus, R. L.↗

NASA NTRS1984

Global Biogeochemistry: an Overview

The dynamic biogeochemical equilibria among the major pools of carbon, nitrogen, sulfur, and phosphorus represented by terrestrial biomes, the world's oceans, and the troposphere are disturbed. Since even the most rapid processes of adjustments among the reservoirs take decades, new equilibria are far from established. These human-induced perturbations and the system's subsequent responses constitute an on-going biogeochemical experiment at the global level. Current and new information must be combined in a way that allows testing of various hypotheses about the workings of global biogeochemical systems. This enables assessment of current knowledge and evaluation of the gaps.

Moore, B., III↗

NASA NTRS1985

Was core formation violent enough to homogenize the early mantle?

The dynamics of iron, its thermal state and its phase in the accreting Earth probably played a major role in the Earth's early thermal evolution. Plausible impact thermal histories make it possible that pure iron was molten in the accreting Earth after it was about 10% grown. Hence, iron eutectic alloys (FeS, FeO) certainly were. Additionally, the initial temperature of the core is an important constraint on the secular cooling of the early Earth and on the strength of the early geodynamo. Whether iron is solid or molten would influence geochemical equilibria in the upper and lower mantle; the mode of core formation, by spherical or near-spherical blobs, stalk-like instabilities, or something more catastrophic would influence the partitioning of siderophiles between silicate and iron phases. Early descent of iron (during accretion) favors partitioning according to low-pressure phase equilibria, whereas late descent favors higher pressure. The later core formation occurs, the greater the heat pulse, due to the strong dependence of gravitational potential energy on planetary radius. The heat may homogenize the mantle if core formation is global; otherwise, heterogeneity of iron differentiation may leave some of the pre-archean mantle unaffected. The larger the chunks of proto-core (and hence smaller surface/volume ratios) the greater the heterogeneity.

Cooperman, S. A.↗

NASA NTRS1985

The utilization of satellite data and dynamics in understanding and predicting global weather phenomena

A two layer spectral quasi-geostrophic model is used to simulate the effects of topography on the equilibria, the stability, and the long term evaluation of incipient unstable waves. The flow is forced by latitudinally dependent radiational heating. The nature of the form drag instability of high index equilibria is investigated. The proximity of the equilibrium shear to a resonant value is essential for the instability, provided the equilibrium occurs at a slightly stronger shear than resonance. The properties of the steady Hadley and Rossby required for a thermally forced rotating fluid on a sphere are further explained. An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. Also documented is a steady, axisymmetric model of the general circulation developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. Included is further research on cloud street phenomena. Orientation angles and horizontal wavelengths of boundary layer rolls and cloud streets are determined from an analysis of a truncated spectral model of three dimensional shallow moist Boussinesq convection in a shearing environment is further explained. Relatively broadly spaced roll clouds have orientations for which the Fourier component of the roll perpendicular shear is nearly zero, but the second corresponds to narrowly spaced rolls having orientations for which the Fourier coefficients of both the perpendicular and the parallel components of the shear are nearly equal.

Shirer, H. N.↗

NASA NTRS1985

The effect of topography on the evolution of unstable disturbances in a baroclinic atmosphere

A two layer spectral quasi-geostrophic model is used to simulate the effects of topography on the equilibria, their stability, and the long term evolution of incipient unstable waves. The flow is forced by latitudinally dependent radiative heating. Dissipation is in the form of Rayleigh friction. An analytical solution is found for the propagating finite amplitude waves which result from baroclinic instability of the zonal winds when topography is absent. The appearance of this solution for wavelengths just longer than the Rossby radius of deformation and disappearance of ultra-long wavelengths is interpreted in terms of the Hopf bifurcation theory. Simple dynamic and thermodynamic criteria for the existence of periodic Rossby solutions are presented. A Floquet stability analysis shows that the waves are neutral. The nature of the form drag instability of high index equilibria is investigated. The proximity of the equilibrium shear to a resonant value is essential for the instability, provided the equilibrium occurs at a slightly stronger shear than resonance.

Clark, J. H. E.↗