Search NASA⌕ Search

SEARCH · Search NASA

Results for “isotherm 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Facet-Dependent Adsorption of Pb(II) on Hematite (001), (116), and (104) Surfaces

Hematite’s common (001) and (012) facets are frequently used in model studies of lead (Pb) adsorption behavior, but there is a lack of research on the high-energy facets, e.g. (104), present in nature. Also, few studies have attempted to connect the molecular details of facet-specific Pb adsorption to macroscopic uptake behavior. Here, to address these knowledge gaps, we investigated Pb(II) adsorption behaviors on facet-engineered hematite nanoparticles dominated by (001), (104), and (116). Adsorption experiments revealed significant variations in Pb(II) uptake among the three samples, with (001) demonstrating the highest capacity and (116) showing the best adsorption efficiency when normalized to specific surface area. Adsorption kinetics followed the pseudo-second-order model, indicating the adsorption process is governed mostly by chemisorption. Adsorption isotherms were well fitted by the Langmuir model, indicating uptake proceeds until roughly monolayer adsorption. Detailed characterization revealed Pb(II) was adsorbed as single atoms with complex inner-sphere binding modes that varied across different facets, indicating adsorption is both structurally and energetically facet-dependent. Co-adsorption experiments further demonstrated Cu 2+ , Zn 2+ , and humic acid significantly promoted Pb(II) adsorption. This study advances the understanding of hematite surface reactivity in controlling macroscopic wet adsorption behaviors, providing valuable insights into the environmental fate of Pb(II).

Contaminant removal↗
DOE OSTI2023

Effect of Mn on eutectic phase equilibria in Al-rich Al-Ce-Ni alloys

Microstructural analysis of additively manufactured (AM) Al-Ce-Ni-Mn alloys has identified phases not predicted from existing ternary liquidus projections in the Al-Ce-Ni system. Because the rapid cooling rate of AM is orders of magnitude above that of traditional casting, it is unclear if these additional phases arose from the non-equilibrium processing conditions of AM, a drastic shift in phase stability in the system due to the addition of 1 wt% Mn, or some combination of these two influences. The phases and microstructure of cast samples of Al-Ce-Ni and Al-Ce-Ni-Mn alloys were characterized for several annealing conditions which revealed the equilibrium phases at different temperatures. Phase analysis confirmed that minute levels of Mn substituted for Ni in the system drastically shifts the liquidus projection in the Al-rich corner of the ternary phase diagram such that the eutectic Al 3 Ni phase is suppressed in favor of the Al 23 Ni 6 (Ce,Mn) 4 phase. Further addition of Mn promotes the formation of Al 20 Mn 2 Ce and Al 10 Mn 2 Ce phases. The phase analysis data was then used to improve the CALPHAD modeling of the liquidus projection and isothermal sections for the Al-rich Al-Ce-Ni-Mn quaternary system. Thermodynamic modeling and experimental analysis on phases in the AM sample of Al-Ce-Ni with Mn confirmed that the phases present are consistent with Mn-containing Al-Ce-Ni cast samples. Here, this investigation demonstrates the potential for using secondary alloying elements to drastically alter phase stability and microstructure in alloy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗
DOE OSTI2022

Thermal Pressure in the Thermal Equation of State for Solid and a Proposed Substitute

Abstract The thermal equation of state (TEOS) for solids is a mathematic model among pressure, temperature and density, and is essential for geophysical, geochemical, and other high pressure–temperature (high P–T) researches. However, in the last few decades, there has been a growing concern about the accuracy of the pressure scales of the calibrants, and efforts have been made to improve it by either introducing a reference standard or building new thermal pressure models. The existing thermal equation of state, P ( V , T ) = P ( V , T 0 ) + P th ( V , T ), consists of an isothermal compression and an isochoric heating, while the thermal pressure is the pressure change in the isochoric heating. In this paper, we demonstrate that, for solids in a soft pressure medium in a diamond anvil cell, the thermal pressure can neither be determined from a single heating process, nor from the thermal pressure of its calibrant. To avoid the thermal pressure, we propose to replace the thermal pressure with a well-known thermal expansion model, and integrate it with the isothermal compression model to yields a Birch–Murnaghan-expansion TEOS model, called VPT TEOS. The predicted pressure of MgO and Au at ambient pressure from Birch–Murnaghan-expansion VPT TEOS model matches the experimental pressure of zero (0) GPa very well, while the pressure prediction from the approximated Anderson PVT TEOS exhibit a big deviation and a wrong trend.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A simple model for short-range ordering kinetics in multi-principal element alloys

Short-range ordering (SRO) in multi-principal element alloys influences material properties such as strength and corrosion. While some degree of SRO is expected at equilibrium, predicting the kinetics of its formation is challenging. We present a simplified isothermal concentration-wave (CW) model to estimate an effective relaxation time of SRO formation. Estimates from the CW model agree to within a factor of five with relaxation times obtained from kinetic Monte Carlo (kMC) simulations when above the highest ordering instability temperature. Further, the advantage of the CW model is that it only requires mobility and thermodynamic parameters, which are readily obtained from alloy mobility databases and Metropolis Monte Carlo simulations, respectively. The simple parameterization of the CW model and its analytical nature makes it an attractive tool for the design of processing conditions to promote or suppress SRO in multicomponent alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparison of Heterogeneous Langmuirian Models for Mixed-Gas Adsorption Equilibria

Mixed-gas adsorption remains a very active area of research, where newly proposed adsorption equilibrium models need to be compared to standards in the field. The generalized Langmuir (gL) isotherm, proposed by Hamid et al. to address adsorption of different size molecules on heterogeneous adsorbents without requiring evaluation of the Gibbs isotherm, has thus far been tested mainly against the adsorbed solution theory (AST). The gL isotherm was shown to be thermodynamically consistent and gave either better or similar representation than AST models. Here, the gL isotherm is contrasted with the Dual-Process Langmuir isotherm popularized by Ritter and his co-workers and the Loading Ratio Correlation popularized for simple mixtures, both of which are considered the “standard” models for process simulation. Furthermore, the comparison is based on their theoretical roots, interpretations, and representation of the equilibrium loading and the isosteric enthalpy of adsorption of unary, binary, and ternary mixtures of various gaseous adsorbates on different adsorbents over a range of experimental pressures and temperatures.

42 ENGINEERING↗
DOE OSTI2023

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Addressing spreading pressure dependence of real adsorbed solution theory with generalized Langmuir isotherm

This work addresses spreading pressure dependence of Real Adsorbed Solution Theory (RAST) for mixed-gas adsorption equilibria using generalized Langmuir (gL) isotherm. Considering vacant sites as an integral part of competitive multicomponent adsorption on a constant adsorbent surface area, the gL isotherm properly accounts for surface heterogeneity and loading, adsorbate composition, and temperature dependence. We show the spreading pressure dependence of adsorbate activity coefficient expression in the RAST framework can be generated from the gL isotherm. The procedure is illustrated with a spreading pressure dependent adsorption Nonrandom Furthermore, two-Liquid activity coefficient model, and the results are validated for ten binary mixed-gas adsorption equilibria systems including two highly nonideal azeotropic systems.

42 ENGINEERING↗
DOE OSTI2022

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Evaluating Pretreatment Strategies with Modeling for Reducing Scaling Potential of Reverse Osmosis Concentrate: Insights from Ion Exchange and Activated Alumina

Reverse osmosis concentrate (ROC) treatment is critical for enhancing water recovery and minimizing concentrate volume for disposal, especially in regions facing water scarcity. This study investigates the application of ion exchange (IX) resins and activated alumina (AA) as pretreatment strategies to mitigate scaling in ROC due to high concentrations of total dissolved solids, hardness (Ca 2+ and Mg 2+ ), and silica. Through a series of Langmuir isotherms, continuous column experiments, and model simulation, two types of strong acid cation IX resins and three types of strong base anion IX resins alongside three types of AA were evaluated. Results indicate that AA exhibits superior performance in silica removal, achieving up to a 65% reduction and maintaining performance for up to 800 bed volume without reaching saturation. Model simulation of a secondary reverse osmosis treating ROC after the IX and AA pretreatment indicated an additional water recovery of ~70% using antiscalants. This study demonstrates the potential for achieving higher water recovery while also identifying opportunities for pretreatment improvement. Challenges such as the limited IX capacity treating ROC, which requires frequent regeneration and increases operational costs, along with the restricted regeneration capacity of AA, underscore the importance of innovation. These findings emphasize the critical need for developing advanced materials and optimized strategies to further enhance the efficiency of ROC treatment processes.

activated alumina↗
DOE OSTI2026

Modeling glass degradation and release of radionuclides from vitrified waste for performance assessment simulations

The release of radionuclides initially encapsulated in a slowly degrading solid waste form and contained in an eventually corroding canister defines the source term for numerical simulations for the assessment of a geologic repository for high-level radioactive waste. While the details of waste degradation, canister corrosion, and dissolution and mobilization of the radionuclides in pore water include complex chemical reaction and transport processes that are coupled to the thermal, hydrological, microbiological, and mechanical conditions in the repository, the source-term model suitable for use in a numerical performance assessment model should be a defensible abstraction of these mechanisms. We developed a radiological source-term model and implemented it into a non-isothermal flow and transport simulator. While the proposed source-term model is applicable to various waste forms, canister systems, and disposal concepts, we specifically considered radionuclide releases from vitrified high-level waste placed in a cylindrical canister disposed in a deep vertical borehole repository. In this model, waste degradation is a function of temperature, and it can be adjusted to evaluate the influence of and propagate uncertainties in pH, passivation reactions, and chemical conditions as well as geometrical factors. The time-dependent, congruent release of safety-relevant radionuclides present in the decaying inventory is then calculated. Finally, the radionuclides are mobilized by diffusive and advective transport according to the thermo-hydraulic conditions prevailing in the near field of the repository, from where they migrate through the geosphere to the accessible environment. We examine the influence of the source-term model’s parameters on performance assessment calculations through sensitivity and uncertainty propagation analyses, identifying influential factors and confirming the upper bound of their impact. These considerations align with the overarching goal of repository design, which is to demonstrate that engineered and natural barriers can collectively delay radionuclide migration for timescales far exceeding human planning, thereby providing multiple, redundant barriers against environmental contamination.

iTOUGH2↗
DOE OSTI2023

Insights into co-pyrolysis of polyethylene terephthalate and polyamide 6 mixture through experiments, kinetic modeling and machine learning

The non-isothermal pyrolysis of polyethylene terephthalate (PET), polyamide 6 (PA6), and their mixtures was studied in a thermogravimetric analyzer at different heating rates. Temperature of maximum decomposition (T max ) decreased by 25–45 °C and 35–55 °C for the PET:PA6 mixtures (3:1, 1:1, 1:3) compared to PET and PA6, respectively. The kinetic analysis was initially carried out using isoconversional method. However, the dependency of activation energy on conversion was observed for the co-pyrolysis of PET and PA6 that suggested the occurrence of multi-step reactions in the mixtures. Distributed activation energy model (DAEM) was used in this study to describe the multistep reactions occurring during pyrolysis of PET:PA6 mixtures. Here, in this work, a four-parallel reaction DAEM was developed to describe the pyrolysis kinetics of PET:PA6 mixtures. The apparent mean activation energies (E o ) for PET, PA6, and mixtures varied in the range of 244–255, 140–215, and 138–255 kJ mol –1 , respectively. The mass loss profiles of PET and PA6 mixtures were also modeled using artificial neural network (ANN). Out of 155 ANN models, the best prediction was made by ANN511 with R 2 greater than 0.997 for both test and unseen data. The interaction effects observed through TGA experiments and subsequent kinetic analysis were further assessed in terms of product composition using analytical pyrolysis coupled with gas chromatograph/mass spectrometer (Py-GC/MS). Co-pyrolysis of PET and PA6 resulted in the formation of new aromatic compounds with nitrogen-containing functional groups, which were not detected when PET or PA6 were pyrolyzed individually.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Role of Water Volume Fraction on Water Adsorption in Anion Exchange Membranes

Water absorption plays a key role in the performance of polymeric anion exchange membranes. It influences important properties such as ionic conductivity and mechanical strength and alters their performance as solid electrolytes in hydrogen electrochemical devices for energy conversion. However, computational approaches that address the relationship between the polymer design and the absorption process are scarce. In this work, we introduce a simple thermodynamic model to predict the water absorption isotherms of polyelectrolyte membranes in contact with a water vapor reservoir that incorporates the specific chemical design of the polymers. The model accurately predicts the water content and macrostructural properties of polynorbornene membranes as a function of the water activity and successfully captures the effect of various polymer design parameters. The energy of pairwise attractive interactions predicted by our model provides a means to interpret the absorption process at the molecular level. The model also reveals the most significant favorable and unfavorable contributions to the free energy and indicates that their balance is solely governed by the water volume fraction, regardless of the polymer design. This universal behavior leads to important implications in the search for better ion exchange membranes.

Zaldivar, Gervasio [New York University, Brooklyn,↗
DOE OSTI2022

Oxidation kinetics of SPS-densified U 3 Si 2 fuels—Microstructure impact

Here, U 3 Si 2 is a potential candidate for accident tolerant fuels because of its high uranium density and excellent thermal conductivity in comparison to UO 2 . However, U 3 Si 2 suffers from oxidation, steam corrosion, and subsequent disintegration/pulverization. The detailed investigation of kinetics that incorporates fundamental treatment of oxidation of U 3 Si 2 is scarcely reported, and the oxidation mechanisms have not been fully elucidated. In this paper, the oxidation behavior of microcrystalline (mc - ) and nanocrystalline (nc - ) U 3 Si 2 have been systematically investigated using a thermogravimetric analysis (TGA) apparatus through a series of isothermal and non-isothermal kinetic studies. The isothermal kinetic study with a model-fitting approach indicates oxidation activation energy of 85 kJ/mol for dense mc - U 3 Si 2 and 96.4 kJ/mol for nc - U 3 Si 2 pellets, while the isoconversional approach leads to an activation energy in the range of 70–85 kJ/mol for mc - U 3 Si 2 and 75–86 kJ/mol for nc - U 3 Si 2 with three most common model-free methods, including Kissinger–Akahira–Sunose, Flynn–Wall–Ozawa, and Friedman methods. The derivation of oxidation activation energies using both isothermal and isoconversional methods highlights the approach to evaluate the oxidation resistance of nuclear materials using TGA quantitatively and makes it possible to compare among various nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Phase change material integrated core–shell catalyst for in situ thermal control in methanol synthesis from syngas

A model-guided core–shell catalyst design is presented for methanol synthesis, featuring a phase change material (PCM) core encapsulated by a Cu–Zn–Al 2 O 3 (CZA) catalytic shell. The PCM enables in situ thermal management by absorbing reaction heat at its melting point, mitigates the kinetic decline at high temperatures and therefore avoids low conversion, prevents hot spots, and stabilizes the reaction temperature. A two-dimensional axisymmetric, non-isothermal packed-bed reactor model (COMSOL 6.3) was developed for a 10 g system. Simulations evaluate three PCM candidates, that is, LiNO 3 , 9 wt% LiCl + 91 wt% LiNO 3 , and commercial H250, with melting points near 244–250°C. Results indicate that CO conversion can increase from 34.4% to 52.4%, and methanol production can improve by 69% compared to a conventional packed-bed reactor. Beyond methanol synthesis, the PCM-integrated core–shell concept provides a scalable approach for thermal control in exothermic reactions, improving reactor efficiency and safety.

core–shell catalyst↗
DOE OSTI2021

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗
DOE OSTI2025

Aggressively Dissipative Dark Dwarfs: The Effects of Atomic Dark Matter on the Inner Densities of Isolated Dwarf Galaxies

We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of cold dark matter and a strongly dissipative subcomponent. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as atomic dark matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of ∼10 10 M ⊙ systems that vary over the aDM microphysics and the dwarf’s evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts, and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark “clumps,” which are clusters of subresolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction (∼5%) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.

cold dark matter↗
DOE OSTI2024

Development of a near-isothermal transcritical CO 2 compression system with a liquid piston compressor

Compressors are critical components in vapor compression cycle systems, significantly contributing to energy consumption. As global demand for HVAC&R systems rises, enhancing compressor efficiency becomes increasingly vital. This paper introduces a novel liquid piston compressor integrated with a gas cooler for the transcritical CO 2 refrigeration cycle. The liquid piston enables CO 2 refrigerant compression within various types of heat exchangers, facilitating the transfer of compression heat to the heat transfer fluid. By releasing significant heat, this compressor allows for removing or downsizing the traditional gas cooler in the refrigeration system. This paper presents the experimental performance of the first near-isothermal compressor utilizing a liquid piston in a vapor compression cycle. The critical parameters affecting heat transfer are analyzed by using a 1-D simulation model to achieve a near-isothermal compression process. The results show that the developed prototype successfully reduced the compression temperature increase from 95 K to 10 K, achieving 90 % isothermal efficiency. Furthermore, the 1-D simulation results suggest the smaller internal diameter tubes benefit the isothermal efficiency the most.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗