SEARCH · Engineering Papers
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.
CO 2 BOLs are promising CCS solvents, whose mass transfer rates improve upon cooling. We show this unintuitive behavior is due to slow transport of reaction products from the gas–liquid interface, caused by shifts in reaction equilibria.
It has been recently suggested that hydroxide ions can be formed in the electrolyte of molten carbonate fuel cells when water vapor is present. The hydroxide can replace carbonate in transporting electrons across the electrolyte, thereby reducing the CO 2 separation efficiency of the fuel cell although still producing electricity. In this work, we obtain the equilibrium concentration of hydroxide in five molten alkali carbonate salts from molecular simulations. The results reveal that there can be a substantial amount of hydroxide in the electrolyte at low partial pressures of CO 2 . In addition, we find that the equilibrium concentration of molecular water dissolved in the electrolyte is over two orders of magnitude higher than that of CO 2 . Increasing the size and polarizability (or in other words reducing the “hardness”) of the cations present in the electrolyte can reduce the hydroxide fraction, but at the cost of lowering ionic conductivity.
Partial molar volumes, energies, and enthalpies can be computed from NpT-Gibbs ensemble simulations through a post-processing procedure that leverages fluctuations in composition, total volume, and total energy of a simulation box. By recording the instantaneous box volumes V and instantaneous number of molecules N i of each of n species for M frames, a large M ×n matrix N is constructed, as well as the M × 1 vector V. The 1 × n vector of partial molar volumes V¯ may then be solved using N • V¯=V. A similar construction permits calculation of partial molar energies using M instantaneous measurements of the total energy of the simulation box, and N • U¯=U. Partial molar enthalpies may be derived from U¯, V¯, and pressure p. These properties may be used to calculate enthalpy and entropy of transfer (absorption, extraction, and adsorption) for species in complex mixtures. The method is demonstrated on three systems in the NpT-Gibbs ensemble: a highly compressible natural gas condensate of methane, n-butane, and n-decane, the liquid-phase adsorption of 1,5-pentanediol and ethanol onto the MFI zeolite, and a relatively incompressible mixture of ethanol, n-dodecane, and water at liquid-liquid equilibrium. Property predictions are compared to those from numerical differentiation of simulation data sequentially changing the composition and from equations of state. Here, the method can also be extended to reaction equilibria in a closed system and is applied to a reactive first-principles Monte Carlo simulation of compressed nitrogen/oxygen.
Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.
Macromolecular crowding is known to modulate chemical equilibria, reaction rates, and molecular binding events, both in aqueous solutions and at lipid bilayer membranes, natural barriers that enclose the crowded environments of cells and their subcellular compartments. Previous studies on the effects of macromolecular crowding in aqueous compartments on conduction through membranes have focused on single-channel ionic conduction through previously formed pores at thermodynamic equilibrium. Furthermore, the effects of macromolecular crowding on the mechanism of pore formation itself were studied using the droplet interface bilayer (DIB) technique with the voltage-dependent pore-forming peptide alamethicin (alm). Macromolecular crowding was varied using 8 kDa molecular weight polyethylene glycol (PEG8k) or 500 kDa dextran (DEX500k) in two aqueous droplets on both sides of the bilayer membrane. In general, voltage thresholds for pore formation in the presence of crowders in the droplets decreased compared to their values in the absence of crowders, due to excluded volume effects, water binding by PEG, and changes in the ordering of water molecules and hydrogen-bonding interactions involving the polar lipid headgroups. In addition, asymmetric crowder loading (e.g., PEG8k–DEX500k on either side of the membrane) resulted in transmembrane osmotic pressure gradients that either enhanced or degraded the ionic conduction through the pores.
In the dynamic environment of multi-component reactive molten salts, speciation unfolds as a complex process, involving multiple competing reaction pathways that are likely to face free energy barriers before reaching the reaction equilibria.
Water chemistry plays a critical role in the design and operation of water treatment processes. Detailed chemistry modeling tools use a combination of advanced thermodynamic models and extensive databases to predict phase equilibria and reaction phenomena. The complexity and formulation of these models preclude their direct integration in equation-oriented modeling platforms, making it difficult to use their capabilities for rigorous water treatment process optimization. Neural networks (NN) can provide a pathway for integrating the predictive capability of chemistry software into equation-oriented models and enable optimization of complex water treatment processes across a broad range of conditions and process designs. Herein, we assess how NN architecture and training data impact their accuracy and use in equation-oriented water treatment models. We generate training data using PhreeqC software and determine how data generation and sample size impact the accuracy of trained NNs. The effect of NN architecture on optimization is evaluated by optimizing hypothetical black-box desalination processes using a range of feed compositions from USGS brackish water data set, tracking the number of successful optimizations, and testing the impact of initial guess on the final solution. Our results clearly demonstrate that data generation and architecture impact NN accuracy and viability for use in equation-oriented optimization problems.
Computer methods were developed and demonstrated that make it possible to predict the combustion and pyrolysis chemistry of proposed biofuels, even before the biofuels are ever manufactured or experimentally tested. The methods are based on first-principles calculations of reaction rates and equilibria. A large number of complicated calculations are required, so we developed a computer workflow that automates many of the steps. With the chemical structure of the biofuel and the system conditions (temperature, pressure, compositions, etc.) as the input, the computer methods can construct a detailed kinetic model with refined parameters without intense human intervention. The automation of model construction makes it practical for a single engineer to make the predictions in a reasonable time period. The accuracy of the computer predictions was tested by measuring the time-profiles of various reaction intermediates (carbon monoxide, carbon dioxide, water, formaldehyde, methane, and ethylene) formed as the biofuels pyrolyze or ignite using sophisticated laser shock tube methods. Several biofuels, including methyl propyl ether, cyclopentanone, butyl acetate isomers, etc., were carefully studied to demonstrate the validity of the proposed computational and experimental methodology. This research project opens up the possibility for assessing some proposed biofuels on the computer without expensive experimental campaigns. The new computer methods are open source and well-documented, allowing a knowledgeable modeler to easily access the developed methods and implement them in their work routine. Moreover, they are also convenient and powerful, promising for developing the fuel chemistry models needed for use in engine simulators, facilitating co-optimization of new engine designs with new fuels.
Alkali metal doping and grain boundaries (GB) have been at the center of attention within the Cu(In,Ga)(S,Se) 2 photovoltaics community for years. This study provides the first experimental evidence that the GB of sodium-doped CuInSe 2 thin films may undertake reversible oxidation even at room temperature, whereas undoped films may not. The findings are corroborated by cathodoluminescence imaging, secondary ion mass spectrometry, and Kelvin probe force microscopy on air-exposed films subsequently subject to vacuum. A thermochemical assessment identifies the likely solid-gas equilibria involved. These reactions open new research questions with respect to the beneficial role played by alkali metal dopants in chalcopyrite solar cells and may steer the community toward new breakthroughs.
Rare earth elements (REE) and many other critical minerals are necessary for the manufacturing of modern everyday technologies, including microchips, batteries and electric motors. Recovery of these materials typically involves aqueous systems which can be modeled as chemical equilibrium problems. One common method for solving these problems is the law of mass action approach (LMA), where a system of non-linear equations involving the equilibrium constants is solved. However, despite being theoretically simple, these problems are in practice very difficult to solve. Currently, the use of iterative heuristics based on saturation indices to decide on which species and reactions to include in the calculations is the state of the art to arrive at a solution. Here, we present an optimization-based alternative to solve chemical equilibria problems involving precipitation/dissolution reactions without the need for such heuristics. Our approach is first validated against the LMA software MINTEQ and PHREEQC for a number of case studies, and then applied to a novel REE recovery process reported in the literature. Overall, our approach was found to have close agreement with MINTEQ and PHREEQC, and we were able to successfully replicate the reported yield and purity for the published REE process.
BioSTEAM is a fast and flexible package for the design, simulation, and techno-economic analysis of biorefineries under uncertainty. BioSTEAM’s framework is built to streamline and automate early-stage technology evaluations and to enable rigorous sensitivity and uncertainty analyses. Complete biorefinery configurations are available at the Bioindustrial-Park GitHub repository, BioSTEAM’s premier repository for biorefinery models and results. The long-term growth and maintenance of BioSTEAM is supported through both community-led development and the research institutions invested in BioSTEAM. Through the open-source and community-lead platform, BioSTEAM aims to foster communication and transparency within the biorefinery research community for an integrated effort to expedite the evaluation of candidate biofuels and bioproducts. Additionally, an agile life cycle assessment (LCA) platform has been designed to interface with BioSTEAM, BioSTEAM-LCA. This open-source, installable package allows users to perform streamlined LCAs of biorefineries. The focus of BioSTEAM-LCA is to streamline and automate early-stage environmental impact analyses of processes and technologies, and to enable rigorous sensitivity and uncertainty analyses linking process design, performance, economics, and environmental impacts. ThermoSTEAM is a standalone thermodynamic engine capable of estimating mixture properties, solving thermodynamic phase equilibria, and modeling stoichiometric reactions. ThermoSTEAM builds upon chemicals, the chemical properties component of the Chemical Engineering Design Library, with a robust and flexible framework that facilitates the creation of property packages. The Biorefinery Simulation and Techno-Economic Analysis Modules (BioSTEAM) is dependent on ThermoSTEAM for the simulation of unit operations.
Thermosteam is a thermodynamic engine capable of solving mass and energy balances, estimating mixture properties, solving thermodynamic phase equilibria, and modeling stoichiometric reactions. All chemical data in Thermosteam is imported from the chemicals library, an open-source compilation of data and functions for the estimation of pure component chemical and mixture properties. Thermosteam’s fast and flexible platform has enabled the evaluation of conceptual and emerging biochemical production processes. The Biorefinery Simulation and Techno-Economic Analysis Modules (BioSTEAM) — capable of modeling reactors, distillation columns, heat exchangers, and other unit operations — has adopted Thermosteam as its premier thermodynamic engine. Published biorefinery designs modeled in BioSTEAM implement thermodynamic property packages created with Thermosteam, including a cornstover biorefinery for the production of cellulosic ethanol, a lipid-cane biorefinery for the co-production of ethanol and biodiesel, and a wheatstraw biorefinery for the production of cellulosic ethanol.
In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction—the dissociation of molecular nitrogen on ruthenium—at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an “on-top” adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. Finally, these findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.
Effects related to alkali metal doping and grain boundaries have puzzled the chalcopyrite photovoltaics community for a long time. This study is the first to report reversible oxidation of grain boundaries in CuInSe2 thin films. The phenomenon is observed in sodium-doped films, but not in undoped ones. Cathodoluminescence imaging, secondary ion mass spectrometry and Kelvin probe force microscopy analyses are performed on CuInSe2 thin films before and after exposure to vacuum. The findings suggest the existence of yet unidentified solid-gas equilibria. Resolving the nature of such reactions will provide new insights into the mechanism of alkali metal doping and passivation in chalcopyrite solar cells.
Topics presented include the adsorption and reaction of trichlorofluoromethane on various particles, equilibria of the marine multiphase ammonia system, a novel ozone sensor for direct eddy flux measurements, and characterization of the carbonate content of atmospheric aerosols. Also presented are variations in heavy metals concentrations in Antarctic snows, sources of continental dust over Antarctica during the last glacial cycle, an inventory of anthropogenic emissions and air pollution in the USSR, and atmospheric nuclei in the remote free-troposphere.
The global methanol budget is currently unbalanced, with source terms significantly larger than the sinks terms. To evaluate possible losses of gaseous methanol to sulfate aerosols, the solubility and reactivity of methanol in aqueous sulfuric acid solutions representative of upper tropospheric and lower stratospheric aerosols is under investigation. Methanol will partition into sulfate aerosols according to its Henry's law solubility. Using standard uptake techniques in a Knudsen cell reactor, we have measured the effective Henry's law coefficient, H*, for cold (196 - 220 K) solutions ranging between 45 and 70 wt % H2SO4. We have found that methanol solubility ranges from approx. 10(exp 5) - 10(exp 7) M/atm for UT/LS conditions. Solubility increases with decreasing temperature and with increasing sulfuric acid content. Although methanol is slightly more soluble than are acetone and formaldehyde, current data indicate that uptake by clean aqueous sulfuric acid particles will not be a significant sink for methanol in the UT/LS. These solubility measurements include uptake due to physical solvation and any rapid equilibria which are established in solution. Reaction between primary alcohols and sulfuric acid does occur, leading to the production of alkyl sulfates. Literature values for the rate of this reaction suggest that formation of CH3OSO3H is not significant over our experimental time scale for solutions below 80 wt % H2SO4. To confirm this directly, results obtained using a complementary equilibrium measurement technique will also be presented.
Analysis of the Skylab II M553 Experiment samples resulted in the hypothesis that the reduced gravity environment was altering the melting and solidification reactions. A theoretical study was conducted to define the conditions under which such alteration of phase relations is feasible, determine whether it is restricted to space processing, and, if so, ascertain which alloy systems or phase reactions are most likely to demonstrate such effects. Phase equilibria of unary and binary systems with a single solid phase (unary and isomorphous) were considered.
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.