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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale CFD simulation of biomass fast pyrolysis with a machine learning derived intra-particle model and detailed pyrolysis kinetics

Coupling particle and reactor scale models is as essential as reactor fluid dynamics and particle motion for accurate Computational Fluid Dynamic (CFD) simulations of biomass fast pyrolysis reactors due to intraparticle heat transfer and chemical reactions controlling conversion time and product distributions. Direct online coupling of a particle model with a reactor model is computationally expensive, while offline coupling is case-dependent. In this research, solutions from a series of particle pyrolysis simulations were regressed with Artificial Neural Network (ANN). Furthermoer, this machine learning-derived model predicted the same temperature and conversion profiles compared with particle resolved simulation while the isothermal approach overpredicted the temperature by 130 K and underpredicted the conversion time by 30 s. The ANN model was then integrated into CFD simulations of fluidized bed biomass fast pyrolysis with varied feedstocks via coupling PyTorch and MFiX. The averaged error of simulation predicted bio-oil yields with four feedstocks is 6.4%. This multi-scale approach provides an efficient tool for the coupled particle and reactor scale simulations of biomass pyrolysis.

09 BIOMASS FUELS↗
DOE OSTI2024

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Physics↗
DOE OSTI2022

EXPERIMENTAL AND MODELLING ANALYSIS OF A LARGE-SCALE TWO-PHASE LOOP THERMOSYPHON

Liquid pumping requires large quantities of electrical energy, including about 7% of the energy of building heating, ventilation, and air conditioning (HVAC) systems. To replace pumped condenser-cooling tower water loops with a passive alternative system, we implemented a commercial-scale two-phase loop thermosyphon (TPLT). The unit consists of a 13 m riser integrated with a commercially available cooling tower and circulation water heater that simulates heat loads up to 25 kW. In addition to providing passive cooling capabilities, the cooling tower unit is also maintenance free, reliable, and can operate in both dry and wet modes. This study characterizes the performance (temperature difference between the evaporator and condenser) and the two-phase flow behavior of the loop under various refrigerant charges. Sight glasses installed throughout the loop are used to identify the operating flow regimes in the riser and downcomer. Over the range of operating conditions, we identified that there is an optimal refrigerant charge range for a specific heat load at which near-isothermal operation can be achieved. We further developed a model to predict the operating flow rate and gravitational height in the downcomer and compared it with the experimental data. The results show that the model agrees with the experimental data, in particular the threshold gravitational head height at which there will be subcooled liquid at the condenser exit, indicating that it can be used as a predictive tool for identifying the optimal loop charge for a given heat load.

Loop thermosyphon, passive two phase↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling the thermodynamics of the FeTi hydrogenation under para-equilibrium: An ab-initio and experimental study

We report FeTi-based hydrides have recently re-attracted attention as stationary hydrogen storage materials due to favorable reversibility, good sorption kinetics and relatively low costs compared to alternative intermetallic hydrides. Employing the OpenCalphad software, the thermodynamics of the (FeTi) 1-x H x (0 ≤ x ≤ 1) system were assessed as a key basis for modeling hydrogenation of FeTi-based alloys. New thermodynamic data were acquired from our experimental pressure-composition-isotherm (PCI) curves, as well as first-principles calculations utilizing density functional theory (DFT). The thermodynamic phase models were carefully selected based on critical analysis of literature information and ab-initio investigations. Key thermodynamic properties such as dissociation pressure, formation enthalpies and phase diagrams were calculated in good agreement to our performed experiments and literature-reported data. This work provides an initial perspective, which can be extended to account for higher-order thermodynamic assessments and subsequently enables the design of novel FeTi-based hydrides. In addition, the assessed thermodynamic data can serve as key inputs for kinetic models and hydride microstructure simulations.

08 HYDROGEN↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

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

Deep learning of accurate force field of ferroelectric HfO 2

The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
DOE OSTI2022

The formation and evolution of Ni 2 Cr precipitates in Ni–Cr model alloys as a function of stoichiometry characterized by synchrotron x-ray diffraction

Ni-based alloys containing significant amounts of Cr, may precipitate new phases during long-term service at elevated temperatures. In Ni–Cr model alloys, the formation of Ni 2 Cr has been found to impact the material properties, including the strength and ductility after isothermal aging below the critical temperature of 590°C. In this work, we quantify the formation and evolution of long-range ordering in four Ni–Cr binary model alloys (Ni/Cr = 1.8, 2.0, 2.2, 2.4) after isothermal aging up to 10,000h at temperatures between 373 and 475°C. The alloys were characterized by hardness testing and synchrotron-based x-ray diffraction to quantify the impact of Ni 2 Cr phase fraction and precipitate size on mechanical properties. After 500h of isothermal aging, the formation of Ni 2 Cr was detected in all four alloys at 475°C. In the stochiometric alloy samples (Ni/Cr = 2.0), the formation of Ni 2 Cr was found after 500 and 3000h aging at 418 and 373°C, respectively. We found that the matrix lattice contraction and Ni 2 Cr phase fraction both saturate after early aging times. This is in stark contrast to the hardness and Ni 2 Cr precipitate size that both continue to increase with increasing aging time. Our results highlight that changes in hardness correlate linearly with Ni 2 Cr precipitate size rather than phase fraction. This important structure-property relationship can potentially help define Ni–Cr-based component lifetimes directly through an understanding of how Ni 2 Cr formation impacts strength and ductility. Here, we find that a precipitation hardening model for critical resolved shear stress with weakly coupled dislocations shows good agreement with the material property changes quantified from experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A simple and practical process modeling methodology for pressure swing adsorption

Although many dynamic models exist for the design and simulation of pressure swing adsorption (PSA) processes, these models involve the solution of a complex system of coupled partial differential equations. Process engineers need a simple, practical, and yet robust short-cut model that helps decide whether to implement a PSA system in a process flowsheet. This work presents a “virtual” moving bed modeling methodology that considers only mass and energy balances and adsorption isotherms to describe the cyclic steady state behavior of PSA systems. Similar to tray efficiencies in distillation calculations, adsorption efficiencies are further introduced to account for system “non-ideality.” Finally, a lab-scale air separation system is used to illustrate the application of this modeling methodology.

42 ENGINEERING↗
DOE OSTI2023

Molecular Rotations, Multiscale Order, Hyperuniformity, and Signatures of Metastability during the Compression/Decompression Cycles of Amorphous Ices

We model, via large-scale molecular dynamics simulations, the isothermal compression of low-density amorphous ice (LDA) to generate high-density amorphous ice (HDA) and the corresponding decompression extending to negative pressures to recover the low-density amorphous phase (LDA HDA ). Both LDA and HDA are nearly hyperuniform and are characterized by a dynamical HBN, showing that amorphous ices are nonstatic materials and implying that nearly hyperuniformity can be accommodated in dynamical networks. In correspondence with both the LDA-to-HDA and the HDA-to-LDA HDA phase transitions, the (partial) activation of rotational degrees of freedom activates a cascade effect that induces a drastic change in the connectivity and a pervasive reorganization of the HBN topology which, ultimately, break the samples’ hyperuniform character. Key to this effect is the rapid rate at which changes occur, and not their magnitude. The inspection of structural properties from the short- to the long-range shows that signatures of metastability are present at all length-scales, hence providing further solid evidence in support of the liquid–liquid critical point scenario. LDA and LDA HDA differ in terms of HBN and structural properties, implying that they are distinct low-density glasses. Our work unveils the role of molecular rotations in the phase transitions between amorphous ices and shows how the unfreezing of rotational degrees of freedom generates a cascade effect that propagates over multiple length-scales. Our findings greatly improve our basic understanding of water and amorphous ices and can potentially impact the field of molecular network-forming materials at large.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Performance of a CO 2 sorbent for indoor air cleaning applications: Effects of environmental conditions, sorbent aging, and adsorption of co-occurring formaldehyde

Indoor air cleaning systems that incorporate CO 2 sorbent materials enable HVAC load shifting and efficiency improvements. This study developed a bench-scale experimental system to evaluate the performance of a sorbent under controlled operation conditions. A thermostatic holder containing 3.15 g sorbent was connected to a manifold that delivered CO 2 -enriched air at a known temperature and relative humidity (RH). The air stream was also enriched with 0.8-2.1 ppm formaldehyde. The CO 2 concentration was monitored in real-time upstream and downstream of the sorbent, and integrated formaldehyde samples were collected at different times using DNPH-coated silica cartridges. Sorbent regeneration was carried out by circulating clean air in countercurrent. Almost 200 loading/regeneration cycles were performed in the span of 17 months, from which 104 were carried out at reference test conditions defined by loading with air at 25°C, 38% RH, and 1000 ppm CO 2 , and regenerating with air at 80°C, 3% RH and 400 ppm CO 2 . The working capacity decreased slightly from 43-44 mg CO 2 per g sorbent to 39-40 mg per g over the 17 months. The capacity increased with lower loading temperature (in the range 15-35°C) and higher regeneration temperature, between 40 and 80°C. The CO 2 capacity was not sensitive to the moisture content in the range 6-9 g/m 3 , and decreased slightly when dry air was used. Loading isothermal breakthrough curves were fitted to three simple adsorption models, verifying that pseudo-first-order kinetics appropriately describes the adsorption process. The model predicted that equilibrium capacities decreased with increasing temperature from 15 to 35°C, while adsorption rate constants slightly increased. The formaldehyde adsorption efficiency was 80%-99% in different cycles, corresponding to an average capacity of 86 ± 36 µg/g. Formaldehyde was not quantitatively released during regeneration, but its accumulation on the sorbent did not affect CO 2 adsorption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Diffraction and thermal effect of a Bessel-Gaussian laser for Ag nanoparticle deposition

Nanoparticles are known to sinter at much lower temperatures than the corresponding bulk or micro size particles. A laser-assisted sintering process is considered in this study to sinter Ag nanoparticles by dispensing Ag paste onto an indium tin oxide-coated Si substrate. The Gaussian beam of a CO 2 laser source is propagated through axicon and biconvex lenses, and the resulting hollow beam is focused on the Ag paste with a hollow parabolic mirror. A Bessel-Gaussian irradiance distribution is obtained at the focal plane of the parabolic mirror due to the interference of the hollow laser cone. The Fresnel diffraction approximation is considered to determine the phasor of the laser and an analytical approach is implemented to calculate the irradiance distribution of the Bessel-Gaussian beam. This irradiance distribution is utilized as a heat source in a heat conduction model and the temperature distribution is analyzed for thin Ag films formed during the laser sintering of Ag nanoparticles. An analytical expression is obtained for the temperature distribution by solving the heat conduction equation using Fourier transform for finite media. The widths of the deposited Ag lines are predicted from the temperature profiles and the model predictions compare well with the experimental results. The isotherms are found to be geometrically noncongruent with convex and concave tips depending on the locally maximum and minimum irradiances of the Bessel-Gaussian beam, respectively. The convex and concave tips, however, appear in the same isotherm for sufficiently high substrate speed relative to the laser beam.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

PRODUCTION OF GLASS-TO-METAL SEALS USING ADDITIVE MANUFACTURING TECHNIQUES

Glass-to-metal (GtM) seals are hermetic barriers between glass and metal components, often used in the electronic and vacuum industries. Creating the seals traditionally requires heating all the components to the glass melting temperature, where the glass will flow and bond to the metal. Due to the manufacturing constraints of the individual components and sealing conditions there are geometric and material restrictions. Additive manufacturing techniques were used to make GtM seals to reduce these restrictions. Hermetic, single pin seals were produced using printed metal shells, by selective laser melting (SLM), and digital light processing (DLP) printed glass preforms. Glass preforms were printed with photosensitive resin mixed with a silicate sealing glass powder in a 1:2 weight ratio. Thermal analysis, TGA/DSC, and screening experiments were used to determine an organic removal schedule. Video techniques determined a sintering schedule by producing a densification curve at different isothermal holds which were fit with a viscous sintering model. Single pin seals were also produced use a laser assisted manufacturing (LAM) process as an alternative to the conventional furnace sealing process. A Nd:YAG laser was focused on a solid glass preform, heating the glass to form the glass-metal bond. The sealing glass was doped with a transition metal oxide, Fe2O3 or CuO, to increase absorption of the laser wavelength. The interface of hermetic LAM seals were analyzed using electron microscopy, compared to conventional prepared seals, finding a reduction of interfacial morphologies.

Read, Aaron Reece↗
DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

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

Characterization of Precipitation-Strengthening Heat-Resistant Austenitic Stainless Steels for Life-Prediction Modeling

In this study, the role of minor alloying additions in 347H stainless steels (UNS34709, ASTM A240/240M) on creep-rupture properties at 650-750°C and microstructure evolution during isothermal exposure at 750°C has been investigated, aiming to provide the experimental dataset as boundary conditions of physics-based modeling for material/component life prediction. Four different 347H heats containing various amounts of boron and nitrogen additions were prepared and evaluated. The combined additions of B and N are found to stabilize the strengthening secondary M23C6 carbides and retarding the transition from M23C6 to sigma phase precipitates during thermal exposure. The observed kinetics of microstructure evolution reasonably explains the improvement of creep-rupture properties of 347H stainless steels with the B and N additions.

Yamamoto, Yukinori↗
DOE OSTI2022

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗