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

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗
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 OSTI2026

Cryogenic-Refined MOSFET Modeling for Oscillator, Frequency Divider, and Amplifier Designs Below 4 K

Capturing device characteristic changes at cryogenic temperatures is crucial for cryo-CMOS circuit designs. In this work, we present an isothermal cryogenic-refined modeling approach for CMOS transistors that is simple, low overhead, and easy to implement while offering the required accuracy for predicting circuit performance at the designated temperatures. Guided by die-level measurement data and circuit design principles, the model introduces corrections to only five critical parameters: threshold voltage, carrier mobility, elevated low-frequency flicker noise, dominant high-frequency shot noise, and subthreshold swing (SS). These refinements are implemented around the foundry-provided SPICE model, which is typically validated only down to about 200 K. With these adjustments, the proposed cryogenic-refined model achieves less than 5% error in both large-signal metrics (I–V characteristics) and small-signal parameters (e.g., transconductance) when compared with device measurements at deep-cryogenic temperatures. The methodology is validated in two advanced technologies: TSMC 40-nm CMOS and GlobalFoundries (GF) 45-nm RF-SOI. We further demonstrate its applicability in three representative RF circuits: a 30-GHz LC oscillator, a high-speed current-mode-logic (CML) frequency divider (FD), and a subthreshold Gb/s amplifier, all showing close agreement between simulated predictions and measurements performed at 4 and 2.5 K. Finally, we believe that the proposed approach is implementation-friendly and can significantly accelerate the development of cryo-CMOS integrated circuits.

circuit modeling↗
DOE OSTI2021

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗
DOE OSTI2021

Analytical modeling of the evolution of the nonlinearity parameter of sensitized stainless steel

Austenitic stainless steels are subject to the precipitation of chromium carbides (M 23 C 6 ) during exposure to high temperatures, causing these alloys to be susceptible to intergranular corrosion due to chromium depletion along grain boundaries. The acoustic nonlinearity parameter, β, shows sensitivity to the formation of carbides in these alloys. The Thermo-Calc TC-PRISMA module was used to model the nucleation and growth of grain boundary M 23 C 6 carbides. The model was verified with scanning transmission electron microscopy analysis that allowed measurements of the grain boundary precipitates. The paper introduces a reduced-order model of the acoustic nonlinearity based on the formation of misfit dislocations at the interface of the grain boundary precipitate and matrix to explain the change in β during isothermal aging. Furthermore, a direct relationship between the radius of the M 23 C 6 grain boundary carbides and β was observed and verified with nonlinear ultrasound measurements on 304L and 316L stainless steels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Towards a heat- and mass-balanced kinetic model of TATB decomposition

We report TATB (1,3,5-triamino-2,4,6-trinitrobenzene) was thermally degraded by two small-scale analytical methods – simultaneous differential scanning calorimetry and thermogravimetric analysis (SDT) and a hot-stage microscope with Fourier Transform Infrared (FTIR) analysis capabilities. SDT used ramped heating, isothermal soaking, and thermal pretreatment at various conditions. The heat flow and mass loss were monitored during various treatment conditions to derive chemical decomposition kinetics and Arrhenius parameters. FTIR experiments used isothermal heating, and changes were monitored spectroscopically. Solid samples generated at specific conditions were collected from both methods and were analyzed by DMSO extraction followed by chemical speciation by optical and mass spectrometric methods. Characterization provided the following reaction insights: TATB decreases in a sigmoidal pattern in isothermally heated samples. Other soluble products gradually increase in concentration and then abruptly decline in concentration during the second exotherm, such as diamino-dinitro-benzofurazan and amino-nitro-benzodifurazan. FTIR showed gradual changes in the amino and nitro functionality, shifting positions and decreasing intensity for the first 40 min. Then the solid gradually appeared more like an amorphous C with N incorporated, similar to previous studies on thermally degraded TATB-type materials. Extracted residues (DMSO-soluble components removed) examined by FTIR showed an abrupt change in chemical composition between 40- and 45-min isothermal treatment, indicating early forming solids are different than later forming residues. A reliable mass- and energy-balanced global reaction network must include at least two autocatalytic reactions, either in parallel or series, and at least one must have an explicit initiation reaction having a low activation energy.

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

Establishing a Design Strategy for Corrosion Resistant Structural Materials in Molten Salt Technologies

To address the gap in material design strategies for structural components in molten salt-powered technologies, the current work aimed to identify key parameters governing corrosion induced degradation of materials during exposures in molten halide salts. The applicability of the data on dissolution rates of pure metals (e.g., Cr) in a given molten salt to enable predictions of corrosion behavior of multicomponent binary and ternary Fe-based model and Ni-based commercial alloys isothermally exposed in the purified KCl-MgCl 2 in quartz capsules at 700 and 800°C was evaluated. The influence of alloy chemistry on initial dissolution rates and the consequent impact on the time-dependent corrosion induced microstructural evolution in the alloy subsurface was predicted with a coupled thermodynamic-kinetic model. Finally, the model was well able to predict the Cr depletion in all investigated alloys demonstrating the ability of the proposed approach to enable a unique design strategy for structural materials in molten salts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

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

Unraveling the multi-step crystallization mechanism of polytetrafluoroethylene, modified polytetrafluoroethylene, and their nanocomposites with boron nitride nanobarbs: Experimental insights and theoretical analysis

The non-isothermal crystallization behavior and kinetics of polytetrafluoroethylene (PTFE) composites with boron nitride nanobarb (BNNB), a new generation nanostructure with unique surface morphology and mechanical “barbs” have been analyzed, understanding these properties is essential for their high-end applications as thermal interface materials (TIM) for microwave, 5G and microelectronic devices. The analysis of the crystallization parameters includes crystallization onset, peak and end temperatures, crystallization half-life and overall crystallinity of PTFE, modified PTFE and their BNNB composites. The results were further analyzed using theoretical models such as the combined Avrami-Ozawa model. It was found that BNNB supports crystallization in the modified PTFE but shows minimal effect on the crystallization of PTFE. Due to the limitation of the classical theoretical models used above in fully characterizing the multi-step crystallization process of PTFE, an in-depth analysis using the model-free advanced isoconversional computation was used to characterize the PTFE crystallization based on the evolution of activation energy with fractional crystallinity and for the first time with temperature. Three kinetic regions were identified in the crystallization mechanism. Here, this study investigated the molecular organization and microstructural evolution of PTFE, modified PTFE and their composites during non-isothermal crystallization using advanced X-ray scattering measurements. An insight into the changes undergone by the material's microstructural units including crystallite size and morphology, lamellar thickness and lamellar interfacial layer thickness, and crystallographic phase dynamics during non-isothermal cooling from the melt, was provided here in this work. The effect of copolymer modification of PTFE and the inclusion of pristine and functionalized BNNB (a thermally conductive and electrically insulating ceramic) are both new investigations that provide valuable knowledge for the development of materials with strong matrix-nanofiller interaction and guidance for optimizing sintering and cooling cycles, two key steps in PTFE processing that largely affect the material microstructural features. Overall, the result of the three-part investigation demonstrates that BNNB supports crystallization in the modified PTFE up to 20 wt% concentration and at low and high cooling rates typically used in the industrial processing of PTFE.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗
DOE OSTI2022

Thermoelastic properties of zircon: Implications for geothermobarometry

Abstract A thermal-pressure equation of state has been determined for zircon (ZrSiO4) that characterizes its thermoelastic behavior at metamorphic conditions. New pressure-volume (P-V) data from a “Mud Tank” zircon have been collected from 1 bar to 8.47(1) GPa using X-ray diffraction, and elastic moduli were measured from room temperature up to 1172 K by resonance ultrasound spectroscopy. These data were fitted simultaneously with temperature-volume (T-V) data from the literature in EosFit7c using a new scaling technique. The parameters of a third-order Birch-Murnaghan EoS with a Mie-Grüneisen-Debye model for thermal pressure have compressional EoS parameters K0T = 224.5(1.2) GPa, K0T' = 4.90(31) with a fixed initial molar volume V0 = 39.26 cm3/mol and thermal parameters γ0 = 0.868(15), q = 2.37(80), and ΘD = 848(38) K. EoS parameters that describe the variation of unit-cell parameters with pressure and temperature were determined using an isothermal-type EoS. This new EoS confirms that zircons are stiffer than garnets and exhibit a much lower thermal expansion. This results in steep isomekes between zircon and garnets, which makes zircon trapped as inclusions in garnets at metamorphic conditions a good piezothermometer.

Geochemistry & Geophysics↗
DOE OSTI2024

Study on the magnetothermal properties of the D y 1 - x T b x A l 2 series of compounds

Here, in this work, we developed a theoretical model Hamiltonian, in the mean field approximation, to describe the magnetic and magnetocaloric behavior of the series of compounds Dy 1-x Tb x Al 2 (x = 0.00, 0.15, 0.25, 0.30, 0.35, 0.40 and 0.75). We adjusted the exchange parameters λ DyDy , λ TbTb , and λ DyTb to obtain the spin reorientation temperatures (T SR ) and the critical temperature (T C ) for each compound in the series. The results obtained by the Hamiltonian model agree satisfactorily with the experimental results. The heat capacity curves with and without an applied magnetic field, adiabatic temperature variation and isothermal entropy variation were modeled and compared with the experimental data. As the experimental results show, our model was also able to reproduce the change in the spin reorientation process: a first order spin reorientation transition appears for concentrations x = 0.15, 0.25, 0.30, 0.35, and no spin reorientation transitions after x = 0.40.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Automated Fitting of a Semi-empirical Multiphase Equation of State for Carbon

The equation of state (EOS) of carbon is important in high explosive, geophysical, and inertial confinement fusion applications. Within the semi-empirical Sesame framework, the EOS of each phase is represented by a sum of cold, vibrational, and thermal electronic Helmholtz free energy contributions. Each phase has ~5-10 independent parameters that are adjusted to reproduce single-phase data (e.g., thermal expansion, isothermal compression) as well as experimentally- and computationally- derived phase boundaries. Manual calibration of the full multiphase EOS is arduous. We present our progress in development of automated EOS parameterization based on minimization of an objective function. Here, this function encodes deviation of model EOS results from experimental/computational benchmarks. Optimization is implemented as a combination of global (particle swarm) and local optimization techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2022

Isotherm Model for Moisture-Controlled CO 2 Sorption

Moisture-controlled sorption of CO 2 , the basis for moisture-swing CO 2 capture from air, is a novel phenomenon observed in strong-base anion exchange materials. Prior research has shown that Langmuir isotherms provide an approximate fit to moisture-controlled CO 2 sorption isotherm data. However, this fit still lacks a governing equation derived from an analytic model. As such, in this paper we derive an analytic form for an isotherm equation from a bottom-up approach, starting with a fundamental theory for an alkali liquid. In the range of interest relevant to CO 2 capture from air, an isotherm equation for an alkali liquid reduces to a simple analytic form with a single parameter, $K_{\text{eq}}$. In the limit $K_{\text{eq}} \gg 1$, a 2nd order approximation simplifies to a Langmuir isotherm that, however, deviates from experimental data. The isotherm theory for an alkali liquid has been generalized to a strong-base anion exchange material. In a strong-base anion exchange material, water concentration inside a sorbent, [H 2 O], is not large enough to be regarded as constant, which allows us to extend $K_{\text{eq}}$ to $K_{\text{eq(AEM)eff}} = K_{\text{eq(AEM)}} \times$ [H 2 O]$^{–n}$ according to the law of mass action. The final isotherm formula has been validated by experimental data from the literature. For a moisture-controlled CO 2 sorbent, $K_{\text{eq(AEM)eff}}$ varies significantly with moisture content of the sorbent. Depending on moisture level, the observed $K_{\text{eq(AEM)eff}}$ in a specific sorbent ranges from a few times to a few thousand times the value of $K_{\text{eq}}$ of a 2 mol L –1 alkali liquid.

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

Flexocaloric effect in superelastic materials

We present a combined theoretical–experimental study of flexocaloric effects in superelastic materials exhibiting structural transitions. We study a Ginzburg–Landau model combined with a vibrational model for a beam near a ferroelastic transition loaded with a three-point bending setup. We also perform experiments on a Cu–Al–Ni single crystal undergoing a martensitic transition. We measure bent beam profiles, vertical force vs vertical deflection during a slow isothermal process, time evolution of the bending and unbending amplitudes, and the evolution of temperature profiles. We also compute the evolution of heat source and heat sink profiles. Finally, we study the location of acoustic emission events during the bending/unbending experiment. Our observations are consistent with the model predictions and allow us to identify the main physical parameters relevant for flexocaloric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine Learning in Environmental Chemistry: Application to Surface Complexation Modeling

Environmental chemistry – or biogeochemistry – is the scientific discipline typically invoked when examining and quantifying groundwater or surface water contamination, and nutrient cycling in the environment. Over the last three decades, there have been significant advances in mechanistic model development to describe and predict these complex biogeochemical processes. In particular, surface complexation models (SCMs) have been developed to describe the rock/soil surface reactions of metals and radionuclides, and their partitioning between various mobile species in the aqueous phase or immobile species sorbed on solid surfaces. Often represented by a simplified linear isotherm constant – Kd – in reactive transport models, these reactions play a critical role in many environmental science applications; particularly in contamination risk assessments and nuclear waste disposal performance assessments. In the past several decades, efforts by various institutions across the world have focused on developing SCMs based on datasets from laboratory measurements, including the identification of key parameters such as equilibrium constants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

HOLISMOKES IX. Neural network inference of strong-lens parameters and uncertainties from ground-based images

Modeling of strong gravitational lenses is a necessity for further applications in astrophysics and cosmology. With the large number of detections in current and upcoming surveys, such as the Rubin Legacy Survey of Space and Time (LSST), it is pertinent to investigate automated and fast analysis techniques beyond the traditional and time-consuming Markov chain Monte Carlo sampling methods. Building upon our (simple) convolutional neural network (CNN), we present here another CNN, specifically a residual neural network (ResNet), that predicts the five mass parameters of a singular isothermal ellipsoid (SIE) profile (lens center x and y, ellipticity e x and e y , Einstein radius θ E ) and the external shear (γ ext, 1 , γ ext, 2 ) from ground-based imaging data. In contrast to our previous CNN, this ResNet further predicts the 1σ uncertainty for each parameter. To train our network, we use our improved pipeline to simulate lens images using real images of galaxies from the Hyper Suprime-Cam Survey (HSC) and from the Hubble Ultra Deep Field as lens galaxies and background sources, respectively. We find very good recoveries overall for the SIE parameters, especially for the lens center in comparison to our previous CNN, while significant differences remain in predicting the external shear. From our multiple tests, it appears that most likely the low ground-based image resolution is the limiting factor in predicting the external shear. Given the run time of milli-seconds per system, our network is perfectly suited to quickly predict the next appearing image and time delays of lensed transients. Therefore, we use the network-predicted mass model to estimate these quantities and compare to those values obtained from our simulations. Unfortunately, the achieved precision allows only a first-order estimate of time delays on real lens systems and requires further refinement through follow-up modeling. Nonetheless, our ResNet is able to predict the SIE and shear parameter values in fractions of a second on a single CPU, meaning that we are able to efficiently process the huge amount of galaxy-scale lenses expected in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗