Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “isotherm model parameters”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A quantitative comparison between velocity dependent SIDM cross-sections constrained by the gravothermal and isothermal models

One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross-section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest – the effects of these processes can be highly degenerate with those of dark matter self-interactions at small scales. In this work, we select 9 isolated galaxies and brightest cluster galaxies (BCGs) with baryonic components small enough such that the baryonic gravitational potentials do not significantly influence the halo gravothermal evolution processes. We then constrain the parameters of Rutherford and Møller scattering cross-section models with the measured rotation curves and stellar kinematics through the gravothermal fluid formalism and isothermal method. Cross-sections constrained by the two methods are consistent at $1\sigma$ confidence level, but the isothermal method prefers cross-sections greater than the gravothermal approach constraints by a factor of ~ 3.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Phase field modeling of discontinuous dynamic recrystallization in hot deformation of magnesium alloys

We have developed a quantitative phase field model of discontinuous dynamic recrystallization (PF-DDRX) for investigating grain structure evolution and its impact on mechanical response in magnesium alloy during hot deformation. A set of isothermal compression tests were conducted by Gleeble thermo-mechanical simulator and used to determine PF-DDRX model parameters such as critical stresses of recrystallization, activation energies of dislocation, grain boundaries and recrystallization grain nucleation. In the simulations, the experimentally measured average grain size and distribution were used to generate initial microstructures. The results show that the typical ‘necklace’ microstructures observed in DDRX forms at the earlier stage of recrystallization. The inhomogeneity of microstructures first increases rapidly and then decreases gradually and reaches a constant with the increase of dynamic recrystallization volume fraction. Both predicted grain structure and mechanical response are in good agreement with experimental results. This demonstrates that the model parameters determined by experimental data are reliable for the developed PF-DDRX model to predict the grain structure evolution and flow behavior of AZ80 magnesium alloy under hot deformation. This work presents a path to combine experiments and modeling to give quantitative prediction and minimize experimental cost.

magnesium alloy, Phase field method, dynamic recry↗
DOE OSTI2026

Chemical-Free Lithium Separation from High-Salinity Brines Using Model-Informed and Machine Learning-Optimized Multi-Column Zwitterionic Chromatography

Direct Lithium Extraction (DLE) technologies often struggle to produce high-purity lithium salts from high-salinity brines, as current approaches require chemical-based elution, regeneration, and precipitation steps, resulting in significant environmental footprints. A novel salt fractionation approach using carboxybetaine resin, known as zwitterionic chromatography (ZIC), has demonstrated that lithium ions can be separated from divalent cations under high-salinity conditions using only water as eluent, with no regeneration required. To enable continuous and scalable deployment of this approach, we developed a chemical-free Multi-column Zwitterionic Chromatography (MZC) process and its theoretical and process models. To predict and optimize this nontraditional separation system, we introduced a novel anti-Langmuir isotherm, and the isotherm parameters were estimated through a machine learning-driven optimization based on artificial neural network ensembles with numerical feasibility assessment. Using machine learning-driven optimization, the MZC process achieved 98.0% lithium recovery, 99.5 % Li/(Li + Mg + Ca) purity, a 31.3% productivity increase, and a 33% reduction in water use compared to batch operation. The proposed MZC process enables lithium separation at $0.6-1.2 kg-1 Li, with costs dominated by resin manufacturing, while offering lower separation costs and carbon footprint compared with conventional carbonation. Overall, these findings position the MZC process as an effective polishing step within scalable and sustainable lithium production pipelines.

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

Desolvation Processes in Channel Solvates of Niclosamide

The antiparasitic drug niclosamide (NCL) is notable for its ability to crystallize in multiple 1:1 channel solvate forms, none of which are isostructural. Here, using a combination of time-resolved synchrotron powder X-ray diffraction and thermogravimetry, the process-induced desolvation mechanisms of methanol and acetonitrile solvates are investigated. Structural changes in both solvates follow a complicated molecular-level trajectory characterized by a sudden shift in lattice parameters several degrees below the temperature where the desolvated phase first appears. Model fitting of kinetic data obtained under isothermal heating conditions suggests that the desolvation is rate-limited by the nucleation of the solvent-free product. The desolvation pathways identified in these systems stand in contrast to previous investigations of the NCL channel hydrate, where water loss by diffusion initially yields an anhydrous isomorph that converts to the thermodynamic polymorph at significantly higher temperatures. Taking the view that each solvate lattice is a unique “pre-organized” precursor, a comparison of the pathways from different starting topologies to the same final product provides the opportunity to reevaluate assumptions of how various factors (e.g., solvent binding strength, density) influence solid-state desolvation processes.

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

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Isothermal Versus Sensible Cold Storage: A Model-Based Performance Comparison for Pumped Thermal Electricity Storage

Pumped-thermal energy storage (PTES) systems consume and produce electrical energy using thermal storage media as an intermediate stage. PTES lends itself to long-duration energy storage to facilitate high penetration of intermittent electricity generation. This study presents a model-based comparison of two thermal storage types within a PTES system: a conventional, single-phase, stratified water-glycol sensible storage system (SGS), and an ideal isothermal, two-phase heat exchanger that freezes a water reservoir (isothermal heat exchanger (IHEX)). The SGS thermal storage capacity is based on the liquid’s sensible heat change with temperature, whereas the capacity of the IHEX is based on the latent heat of isothermally freezing and melting water. The idealized IHEX modeled here undergoes steady-state melting and freezing (in contrast to transient rates, as observed with ice-on-coil storage). A computational model of a complete PTES system is presented and used to evaluate the PTES system-level performance with each type of cold storage. Compared to SGS-based PTES, under nominal operating conditions, the IHEX-based PTES increased electrical round-trip efficiency from 61% to 82% and increased energy density from 1.13 to 8.09 kWh/m 3 In conclusion, the performance of the PTES configured with IHEX storage was also analyzed under varying operating parameters.

25 ENERGY STORAGE↗
DOE OSTI2025

Validation of Numerical Tools for Calculating Reactivity Feedback in Sodium Fast Reactors Using SEFOR Experimental Data

The Southwest Experimental Fast Oxide Reactor (SEFOR) was an experimental sodium-cooled fast breeder reactor operated from 1969 to 1972 with experiments designed to measure Doppler reactivity feedback in a wide temperature range from around 350 °F to temperatures approaching the melting point of mixed oxide fuel of around 5000 °F, providing valuable data for code validations. Co-supported by the Department of Energy (DOE) Fast Reactor Program (FRP) and the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the SEFOR benchmark project focused on using the experimental data to validate numerical tools that are used in industry and academia to design and license sodium-cooled fast reactors (SFRs). By the end of FY-25, substantial progress was achieved in the SEFOR benchmark study. A variety of numerical tools commonly used for modeling SFRs were applied to develop models for SEFOR core configurations I-D, I-E, I-I, and I-J. These included Monte Carlo codes such as MCNP, Serpent, and Shift; deterministic codes such as the legacy Argonne Reactor Computation (ARC) suite and the high-fidelity NEAMS code Griffin; and the system analysis code SAS4A/SASSYS-1 (SAS). Using these models, both SEFOR zero-power experiments and power-ascending tests were successfully simulated. Comparisons were performed against experimental measurements of core criticalities, reflector worth, kinetics parameters (Λ/βeff), isothermal reactivity feedback (from 350 °F to 760 °F at zero power), and power-ascending reactivity feedback (as power increased from 0.4 MW to 17 MW). In general, these comparisons demonstrated very good agreement between numerical results and experimental data. In Fiscal Year 26 (FY-26), the SEFOR benchmark project will continue to address the modeling issues identified in FY-25. Effort will focus on the simulation of reactivity insertion transients in SEFOR core II using the ARC/SAS model. Future work will also focus on incorporating BISON into the SEFOR core modeling process to enable the first Multiphysics simulations of the isothermal tests based on the MOOSE framework.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Yttrium solid phase equation of state with uncertainty quantification

We discuss the development of an equation of state (EOS) for solid yttrium up to roughly 80 GPa. The EOS makes use of different experimental datasets, including measurements of the isobaric density, isobaric specific heat, room temperature isotherm, and principal shock Hugoniot. The fitting procedure is carried out using Markov Chain Monte Carlo (MCMC), where we fit model parameters for both the cold curve and ion thermal models simultaneously. The results show close agreement with experimental data and provide meaningful uncertainty estimates on the model parameters. This work serves as a first step towards a multiphase EOS for yttrium, which will include higher pressure solid phases (> 80 GPa), as well as modeling of the liquid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling ethanol/water adsorption in all-silica zeolites using the real adsorbed solution theory

A comprehensive set of single-component and binary isotherms were collected for ethanol/water adsorption into the siliceous forms of 185 known zeolites using grand-canonical Monte Carlo simulations. Using these data, a systematic analysis of ideal/real adsorbed-solution theory (IAST/RAST) was conducted and activity coefficients were derived for ethanol/water mixtures adsorbed in different zeolites based on RAST. It was found that activity coefficients of ethanol are close to unity while activity coefficients of water are larger in most zeolites, indicating a positive excess free energy of the mixture. This observation can be attributed to water/ethanol interactions being less favorable than water/water interactions in the single-component adsorption of water at comparable loadings. The deviation from ideal behavior can be highly structure-dependent but no clear correlation with pore diameters was identified. Furthermore, our analysis also demonstrates the following: (1) accurate unary isotherms in the low-loading regime are critical for obtaining physically sensible activity coefficients; (2) the global regression scheme to solve for activity model parameters performs better than fitting activity models to activity coefficients calculated locally at each binary state point; and (3) including the dependence on adsorption potential offers only a minor benefit for describing binary adsorption at the lowest fugacities. Finally, the Margules activity model was found incapable of capturing the non-ideal adsorption behavior over the entire range of fugacities and compositions in all zeolites, but for conditions typical of solution-phase adsorption, RAST predictions using zeolite-specific or even bulk Margules parameters provide an improved description compared to IAST.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2021

Interface pinning causes the hysteresis of the hydride transformation in binary metal hydrides

Hydriding and dehydriding transitions in bulk and nanocrystalline binary metal hydrides were studied using the Pd-H model system by measuring pressure-composition isotherms with in situ x-ray diffractometry. Nanocrystalline Pd showed a smaller pressure hysteresis, solvus hysteresis, and hysteresis in lattice parameter, compared to bulk Pd. The time-dependence of pressure equilibration was measured after dosing with aliquots of hydrogen, giving equilibration times that were much faster in the single-phase regions than in the two-phase plateaus. In the broad two-phase plateaus, the pressure relaxations were exponential functions of time. An explanation of hysteresis is developed that is based on a dissipative potential barrier that impedes the motion of the interface due to interactions between lattice defects and the two-phase interface. The exponential pressure relaxations and hysteresis are consistent for this mechanism. For a simple model of the pinning potential, the potential barrier maximum is an order of magnitude less than typical grain boundary energies. These pinning effects are substantially different in the nanocrystalline Pd, suggesting differences in the hydriding mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Hydrogen Uptake Kinetics of 1,4-Bis(phenylethynyl)benzene (DEB) Rubberized Coating on Silicone Foam Substrate

The hydrogen uptake kinetics of 1,4-bis(phenylethynyl)benzene, or DEB, mixed with palladium (Pd) on activated carbon in a rubber matrix coating on top of a porous silicone foam substrate are investigated. First, isothermal isobaric hydrogenation experiments were performed under different temperatures and H 2 pressures to extract the uptake kinetics. The H 2 uptake models based on the measured kinetic parameters were then employed to investigate/simulate the performance of the getter under dynamic application environments. The actual hydrogenation characteristics in this type of getter is multifaceted and involve actual H 2 concentration in the getter matrix, micrometer scale diffusion of atomic hydrogen away from Pd sites, precipitation of hydrogenated DEB crystals at the coating surfaces, and mobility of fresh DEB molecules. The kinetic analysis/modeling methodology described in this report can serve as a template for other gas-solid reactions as well. Aside from possessing a good hydrogen capacity and excellent performance, this type of rubberized getter also offers some unique advantages over traditional solid getter: flexible structure and protection of Pd catalyst from exposure to the environment.

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

Engineered Polymer Architectures for Thermo–Responsive Desiccants in Dehumidification Applications

Humidity control is relevant to many processes in buildings and manufacturing and is responsible for over 600 million tons of CO 2 annually. As such, desiccants that can efficiently remove moisture from air can greatly reduce emissions. A new class of desiccant materials, thermo–responsive desiccants, has been of increased interest in the literature. These materials are generally based on thermo–responsive polymers and exhibit a temperature–dependent isotherm, an attribute that is small in traditional desiccants. In this work, seven key parameters needed to create an effective thermo–responsive polymer desiccant are identified using a combination of modeling and experimental data. Further, it is found that more thermo–responsiveness is not necessarily better for energy efficiency, and that an optimal, “medium” value is desired. Additionally, polymer composition and architecture play key roles in adjusting the seven key parameters. Finally, based on findings from both modeling and experimental work, guidelines on synthesizing future thermo–responsive polymer desiccants are presented.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal–organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

LLEPE

LLEPE (liquid-liquid equilibrium parameter estimator) is a tool for estimating thermodynamic parameters like standard enthalpies, standard molar volumes and Pitzer coefficients, by combining Gibbs minimization with regression to multiphase equilibrium isotherm data. It also includes visualization tools to access data quality and prediction model performance.

ILOEJE, CHUKWUNWIKE↗
DOE OSTI2020

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗