Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equilibrium calculations”

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

Emission lines from X-ray-heated accretion disks in low-mass X-ray binaries

We investigate the structure of accretion disks illuminated by X-rays from a central compact object in a binary system. X-rays can photoionize the upper atmosphere of the disk and form an accretion disk corona (ADC) where emission lines can form. We construct a model to calculate the vertical structure and the emission spectrum of the ADC with parameters appropriate to low-mass X-ray binaries. These models are made by nonlocal thermodynamic equilibrium calculations of ion and level populations and include a large number of atomic processes for 10 cosmically abundant elements. Transfer of radiation is treated by using the escape probability formalism. The vertical temperature profile of the ADC consists of a Compton-heated region and a mid-T zone where the temperature is approximately 10(exp 6) K. A thermal instability occurs close to the disk photosphere and causes the temperature of the ADC to drop abruptly from 10(exp 6) K to several times 10(exp 4) K. The emission spectrum in the optical, ultraviolet, extreme ultraviolet, and X-ray range is discussed and compared with the observations.

Ko, Yuan-Kuen↗
NASA NTRS1992

Chemistry of atmosphere-surface interactions on Venus and Mars

Earth-based, earth-orbital, and spacecraft observational data are used in the present evaluation of Venus atmosphere-surface interactions to quantitatively characterize the reactions between C, H, S, Cl, F, and N gases and plausible surface minerals. Calculation results are used to predict stable minerals and mineral assemblages on the Venus surface, in order to ascertain which (if any) of the atmospheric gases are buffeted by mineral assemblages. Chemical equilibrium calculations using extant thermodynamic data on scapolite minerals predict that carbonate-bearing scapolite and sulfate meionite are unstable on the surface of Venus, while chloride-bearing scapolite is stable.

Fegley, Bruce, Jr.↗
DOE OSTI2024

Charge self-consistent density functional theory plus ghost rotationally invariant slave-boson theory for correlated materials

We present a charge self-consistent density functional theory combined with the ghost rotationally invariant slave-boson (DFT+gRISB) formalism for studying correlated materials. Here, this method is applied to SrVO 3 and NiO, representing prototypical correlated metals and charge-transfer insulators. For SrVO 3 , we demonstrate that DFT+gRISB yields an accurate equilibrium volume and effective mass close to experimentally observed values. Regarding NiO, DFT+gRISB enables the simultaneous description of charge-transfer and Mott-Hubbard bands, significantly enhancing the accuracy of the original DFT+RISB approach. Furthermore, the calculated equilibrium volume and spectral function reasonably agree with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1993

Evaporation and combustion of LOX under supercritical and subcritical conditions

The objective is to study the evaporation and combustion of LOX under supercritical and subcritical conditions both experimentally and theoretically. In the evaporation studies, evaporation rate and surface temperature were measured when LOX vaporizing in helium environments at pressures ranging from 5 to 68 atm. A Varian 3700 gas chromatograph was employed to measure the oxygen concentration above the LOX surface. For the combustion tests, high-magnification video photography was used to record direct images of the flame shape of a LOX/H2/He laminar diffusion flame. The gas composition in the post-flame region is also being measured with the gas sampling and chromatography analysis. These data are being used to validate the theoretical model. A comprehensive theoretical model with the consideration of the solubility of ambient gases as well as variable thermophysical properties was formulated and solved numerically to study the gasification and burning of LOX at elevated pressures. The calculated flame shape agreed reasonably well with the edge of the observed luminous flame surface. The effect of gravity on the flame structure of laminar diffusion flames was found to be significant. In addition, the predicted results using the flame-sheet model were compared with those based upon full equilibrium calculations (which considered the formation of intermediate species) at supercritical pressures. Except at the flame front where temperature exceeded 2,800 K, the flame-sheet and equilibrium solutions in terms of temperature distributions were in very close agreement. The temperature deviation in the neighborhood of the flame front is caused by the effect of high-temperature dissociation.

Yang, A. S.↗
DOE OSTI2025

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗
NASA NTRS1989

Collisional rate coefficients of C3H2 and the determination of physical conditions in molecular clouds

Collisional excitation rates for C3H2, calculated using the coupled states approximation at temperatures of 10-30 K, are presented. C3H2 produces a number of spectral line pairs whose members are close together in frequency but arise from levels with different excitation energies. The rates are used in statistical equilibrium calculations to illustrate the excitation properties and density-dependent behavior of various C3H2 line ratios.

Avery, L. W.↗
NASA NTRS1979

Ammonia photolysis and the greenhouse effect in the primordial atmosphere of the earth

Photochemical calculations indicate that in the prebiotic atmosphere of earth ammonia would have been irreversibly converted to N2 in less than 40 years if the ammonia surface mixing ratio were no more than 0.0001. However, if a continuous outgassing of ammonia were maintained, radiative-equilibrium calculations indicate that a surface mixing ratio of ammonia of 0.0001 or greater would provide a sufficient greenhouse effect to keep the surface temperature above freezing. With a 0.0001 mixing ratio of ammonia, 60% to 70% of the present-day solar luminosity would be adequate to maintain surface temperatures above freezing. A lower limit to the time constant for accumulation of an amount of nitrogen equivalent to the present day value is 10 my if the outgassing were such as to provide a continuous surface mixing ratio of ammonia of at least 0.00001.

Kuhn, W. R.↗
NASA NTRS1976

The Langley thermal protection system test facility: A description including design operating boundaries

A description of the Langley thermal protection system test facility is presented. This facility was designed to provide realistic environments and times for testing thermal protection systems proposed for use on high speed vehicles such as the space shuttle. Products from the combustion of methane-air-oxygen mixtures, having a maximum total enthalpy of 10.3 MJ/kg, are used as a test medium. Test panels with maximum dimensions of 61 cm x 91.4 cm are mounted in the side wall of the test region. Static pressures in the test region can range from .005 to .1 atm and calculated equilibrium temperatures of test panels range from 700 K to 1700 K. Test times can be as long as 1800 sec. Some experimental data obtained while using combustion products of methane-air mixtures are compared with theory, and calibration of the facility is being continued to verify calculated values of parameters which are within the design operating boundaries.

Klich, G. F.↗
NASA NTRS1993

Molecular equilibria and condensation sequences in carbon rich gases

Chemical equilibria in stellar atmospheres have been investigated by many authors. Lattimer, Schramm, and Grossman presented calculations in both O rich and C rich environments and predicted possible presolar condensates. A recent paper by Cherchneff and Barker considered a C rich composition with PAH's included in the calculations. However, the condensation sequences of C bearing species have not been investigated in detail. In a carbon rich gas surrounding an AGB star, it is often assumed that graphite (or diamond) condenses out before TiC and SiC. However, Lattimer et al. found some conditions under which TiC condenses before graphite. We have performed molecular equilibrium calculations to establish the stability fields of C(s), TiC(s), and SiC(s) and other high temperature phases under conditions of different pressures and C/O. The preserved presolar interstellar dust grains so far discovered in meteorites are graphite, diamond, SiC, TiC, and possibly Al2O3.

Sharp, C. M.↗
NASA NTRS1979

Statistical equilibrium in cometary C2. II - Swan/Phillips band ratios

Statistical equilibrium calculations have been made for both the triplet and ground state singlets for C2 in comets, using the exchange rate as a free parameter. The predictions of the results are consistent with optical observations and may be tested definitively by accurate observations of the Phillips and Swan band ratios. Comparison with the one reported observation indicates compatibility with a low exchange rate and resonance fluorescence statistical equilibrium.

Swamy, K. S. K.↗
DOE OSTI2025

Self-consistent equilibrium and transport simulations for NSTX-U plasmas enhanced via machine learning surrogate models

The Control-Oriented Transport SIMulator (COTSIM) is an advanced equilibrium and transport code designed for simulating tokamak discharges at computational speeds suitable for control applications. COTSIM’s modular framework enables users to select models that balance accuracy with speed according to specific needs, allowing the code to operate from fast to faster-than-real-time performance levels. This work presents recent enhancements to COTSIM’s predictive accuracy for NSTX-U scenarios, achieved by integrating neural-network-based surrogate models and self-consistent equilibrium calculations. To improve source deposition predictions, a surrogate model for NUBEAM has been incorporated. Additionally, a surrogate model for the Multi-Mode Module (MMM) now supports predictions of anomalous thermal, momentum, and particle diffusivities—key factors for modeling the evolution of temperature and rotation. Each surrogate model was specifically trained for the NSTX-U operational regime to enhance COTSIM’s accuracy while maintaining computational efficiency. Moreover, COTSIM now couples fixed-boundary equilibrium solvers with its transport solvers, enabling self-consistent predictions of plasma profiles and equilibrium evolution over the discharge. Simulation results demonstrate strong agreement between COTSIM and TRANSP predictions for NSTX-U discharges. These substantial advancements expand COTSIM’s utility in model-based control applications for NSTX-U. Potential applications include simultaneous optimization of equilibrium and transport scenarios, integration into digital twins, real-time profile estimation (e.g., temperature and rotation) from limited or noisy measurements, and advanced feedback-based scenario control.

Equilibrium and transport modeling↗
NASA NTRS2012

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

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

Pohorille, Andrew↗
DOE OSTI2025

Catalytic hydrogenation of HMF to BHMF over copper catalysts

2,5-Bis(hydroxymethyl)furan (BHMF) is a bio-derived building block for polyester production, obtained via the hydrogenation of 5-hydroxymethylfurfural (HMF). First-principles thermodynamic equilibrium calculations indicate that this reaction is not thermodynamically limited under relevant conditions (e.g., 100 °C and high H 2 partial pressure). In this work, crude HMF was employed as the feedstock for BHMF synthesis. Initially, acidic impurities and humins were removed from unrefined HMF through filtration using a packed bed of γ-alumina. A comprehensive study of the filtration process is presented, including filtration kinetics, breakthrough curve analysis, and mathematical modeling. The purified HMF was subsequently hydrogenated over a 10 wt% CuZrO 2 catalyst, using ethanol as the reaction solvent. Batch reactions were first performed for collection of kinetic data to guide the transition to continuous flow operation. Kinetic data was collected in a fixed bed reactor at varying contact time, time on stream, temperature, and HMF concentration. This data was used to develop a kinetic model for HMF hydrogenation. Maximum BHMF production rates were achieved at 130 °C, accompanied by minor formation of byproducts from BHMF ring-opening reactions. The BHMF selectivity was 100 % at 100 °C although with lower reaction rates. Furthermore, catalyst stability tests revealed a loss of up to 50 % in catalytic activity within the first 24 h, likely due to the adsorption of HMF-derived oligomers that are not easily removed by filtration.

Crude HMF filtration↗
DOE OSTI2026

Assessing the Accuracy of Property Model Predictions for Cost Optimization of Desalination Technologies

Accurate modeling of seawater thermophysical and thermodynamic properties is critical for optimizing desalination processes. This study compares three seawater property models, a Reaktoro multicomponent model, the thermophysical seawater properties library from the Massachusetts Institute of Technology, and a simplified sodium chloride model, in the context of levelized cost of water (LCOW) minimization for reverse osmosis (RO) and mechanical vapor compression systems. Process simulations and cost optimizations reveal that although all three models yield comparable LCOW and specific energy consumption (SEC) estimates under baseline conditions, deviations among their predictions increase with salinity. Relative differences in LCOW and SEC reach up to 6% and 8%, respectively. RO results show greater variability due to differences in osmotic pressure predictions, which affect pressure constraints at high recoveries. Computational performance varies substantially; specifically, Reaktoro simulations are up to 28 times slower than empirical models due to their detailed equilibrium calculations. These results suggest that empirical models offer acceptable accuracy for routine desalination process design, while Reaktoro provides advantages in scenarios requiring detailed speciation, such as scaling or pH adjustment studies. These findings underscore the importance of selecting appropriate property models based on the modeling objective of desalination applications and motivate future work integrating thermodynamic rigor with empirical efficiency.

Physical and chemical properties↗
DOE OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2024

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Emissions Characterization for Ammonia Fuel Blends in an Enclosed Swirl-Stabilized Diffusion Burner

This study investigates ammonia flames using the enclosed Sydney swirl burner (ESSB), focusing on detailed global emissions measurements. Emissions were measured via Fourier-transform-infrared (FTIR) spectrometer, utilizing a heated, long-path absorption cell. Hot, wet measurements of pertinent species were collected, and concentrations were determined via lineshape fitting in conjunction with the HITRAN database. Results show that partially cracked ammonia compositions yield lower emissions compared to pure NH3/H2, and small amounts of NH3 addition to CH4-containing fuel blends exhibit high levels of NOx and CO, with measurable HCN and unburnt CH4. Equilibrium calculations suggest trade-offs between chemical timescale, heat loss, and mixing. Future work will explore emissions sampling procedures and expand analysis using chemical reactor network modelling.

ammonia combustion↗
DOE OSTI2024

An assessment of ocean alkalinity enhancement using aqueous hydroxides: kinetics, efficiency, and precipitation thresholds

Abstract. Ocean alkalinity enhancement (OAE) is a promising approach to marine carbon dioxide removal (mCDR) that leverages the large surface area and carbon storage capacity of the oceans to sequester atmospheric CO2 as dissolved bicarbonate (HCO3-). One OAE method involves the conversion of salt in seawater into aqueous alkalinity (NaOH), which is returned to the ocean. The resulting increase in seawater pH and alkalinity causes a shift in dissolved inorganic carbon (DIC) speciation toward carbonate and a decrease in the surface ocean pCO2. The shift in the pCO2 results in enhanced uptake of atmospheric CO2 by the seawater due to gas exchange. In this study, we systematically test the efficiency of CO2 uptake in seawater treated with NaOH at aquarium (15 L) and tank (6000 L) scales to establish operational boundaries for safety and efficiency in advance of scaling up to field experiments. CO2 equilibration occurred on the order of weeks to months, depending on circulation, air forcing, and air bubbling conditions within the test tanks. An increase of ∼0.7–0.9 mol DIC per mol added alkalinity (in the form of NaOH) was observed through analysis of seawater bottle samples and pH sensor data, consistent with the value expected given the values of the carbonate system equilibrium calculations for the range of salinities and temperatures tested. Mineral precipitation occurred when the bulk seawater pH exceeded 10.0 and Ωaragonite exceeded 30.0. This precipitation was dominated by Mg(OH)2 over hours to 1 d before shifting to CaCO3,aragonite precipitation. These data, combined with models of the dilution and advection of alkaline plumes, will allow the estimation of the amount of carbon dioxide removal expected from OAE pilot studies. Future experiments should better approximate field conditions including sediment interactions, biological activity, ocean circulation, air–sea gas exchange rates, and mixing zone dynamics.

Ringham, Mallory C. (ORCID:0000000348020185)↗