Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Bubble Point Measurements of Three Binary Mixtures of Refrigerants: R-32/1234yf, R-32/1234ze(E), and R-1132a/1234yf

In this article, bubble point pressures of three binary mixtures at two compositions each have been measured utilizing a static method. The mixtures studied were difluoromethane + 2,3,3,3-tetrafluoroprop-1-ene (R-32/1234yf), difluoromethane + (E)-1,3,3,3-tetrafluoroprop-1-ene (R-32/1234ze(E)), and 1,1-difluoroethene + 2,3,3,3-tetrafluoroprop-1-ene (R-1132a/1234yf). For all of the mixtures, the minimum temperature measured was 270 K. The maximum temperatures measured for the R-32/1234ze(E), R-32/1234yf, and R-1132a/1234yf mixtures were 360, 345, and 325 K, respectively. A total of 100 bubble point pressures are reported and compared to available literature data relative to the REFPROP database. Lemmon, E. W. [ NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, 10.0; National Institute of Standards and Technology, Standard Reference Data Program: Gaithersburg, 2018]. Additionally, the bubble point pressures obtained in this study are modeled using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state. The PC-SAFT model parameters, as well as deviations of this work and literature data from the model, are reported.

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

Vapor and Liquid ( p –ρ– T–x ) Measurements of Binary Refrigerant Blends Containing R-32, R-152a, R-227ea, R-1234yf, and R-1234ze(E)

In this article, The pressure–density–temperature–composition (p–ρ–T–x) data of binary refrigerant mixtures containing R-32 (difluoromethane), R-152a (1,1,-difluoroethane), R-227ea (1,1,1,2,3,3,3-heptafluoropropane), R-1234yf (2,3,3,3-tetrafluoropropene), and R-1234ze(E) (trans-1,3,3,3-tetrafluoropropene) were measured in both the vapor and liquid phases using a two-sinker, magnetic suspension densimeter. The specific samples in this study comprised two compositions of approximately (0.3/0.7) and (0.7/0.3) mole fraction for each of the following four binary refrigerant blends: R-32 + R-1234yf, R-32 + R-1234ze(E), R-1234yf + R-152a, and R-1234ze(E) + R-227ea. Single-phase vapor densities were measured over a temperature range of approximately 253 to 293 K and pressures from 0.05 to 0.98 MPa. Single-phase liquid and supercritical densities were measured over a temperature range of approximately 230 to 400 K and pressures up to 22 MPa; for refrigerant blends containing R-1234yf, the maximum pressure was limited to 14 MPa. Overall relative combined, expanded (k = 2) uncertainties in density ranged from 0.025 to 0.191%, with an average uncertainty of approximately 0.05%. Here, we present measurement results, along with comparisons to available literature data and to default equations of state and mixture models included in REFPROP.

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

Speed of Sound Measurements of Binary Mixtures of trans -1,2-Difluoroethylene (R-1132(E)) with Difluoromethane (R-32) or 2,3,3,3-Tetrafluoropropene (R-1234yf)

The speed of sound of two binary mixtures containing 0.336 mole fraction trans-1,2- difluoroethylene (R-1132(E)) with difluoromethane (R-32) and 0.435 mole fraction R-1132(E) with 2,3,3,3-tetrafluoropropene (R-1234yf) were measured with a dual-path pulse-echo instrument. The speed of sound was measured along several pseudo-isochores for each blend at temperatures ranging from 230 K to 325 K for blends with R-32 and to 342 K for blends with R- 1234yf. Measurements started at pressures just above each mixtures bubble point pressure and were limited to 8 MPa to avoid potential disproportionation reactions of the R-1132(E). The data were compared to multi-fluid models incorporating Helmholtz-energy-explicit equations of state (EOS) for each pure fluid.

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

Evaluation of Force Fields for Molecular Dynamics Simulations of Platinum in Bulk and Nanoparticle Forms

Here, understanding the size- and shape-dependent properties of platinum nanoparticles is critical for enabling the design of nanoparticle-based applications with optimal and potentially tunable functionality. Toward this goal, we evaluated nine different empirical potentials with the purpose of accurately modeling faceted platinum nanoparticles using molecular dynamics simulation. First, the potentials were evaluated by computing bulk and surface properties - surface energy, lattice constant, stiffness constants, and the equation of state - and comparing these to prior experimental measurements and quantum mechanics calculations. Then, the potentials were assessed in terms of the stability of cubic and icosahedral nanoparticles with faces in the {100} and {111} planes, respectively. Although none of the force fields predicts all the evaluated properties with perfect accuracy, one potential - the embedded atom method formalism with a specific parameter set - was identified as best able to model platinum in both bulk and nanoparticle forms.

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

Periodic GFN1-xTB Tight Binding: A Generalized Ewald Partitioning Scheme for the Klopman–Ohno Function

A novel formulation is presented for the treatment of electrostatics in the periodic GFN1-xTB tight-binding model. Periodic GFN1-xTB is hindered by the functional form of the second-order electrostatics, which only recovers Coulombic behavior at large interatomic distances and lacks a closed-form solution for its Fourier transform. We address this by introducing a binomial expansion of the Klopman–Ohno function to partition short- and long-range interactions, enabling the use of a generalized Ewald summation for the solution of the electrostatic energy. This approach is general and is applicable to any damped potential of the form |R n + c| –m . Benchmarks on the X23 molecular crystal dataset and a range of prototypical bulk semiconductors demonstrate that this systematic treatment of the electrostatics eliminates unphysical behavior in the equation of state curves. In the bulk systems studied, we observe a mean absolute error in total energy of 35 meV/atom, comparable to the machine-learned universal force field, M3GNet, and sufficiently precise for structure relaxation. These results highlight the promising potential of GFN1-xTB as a universal tight-binding parametrization.

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

High-Pressure Investigation of 2,4,6-Trinitro-3-bromoanisole (TNBA): Structural Determination and Piezochromism

Understanding phase transitions in energetic materials is crucial for developing predictive models of detonation. 2,4,6-Trinitro-3-bromoanisole (TNBA), an energetic material, was studied in its single-crystal form up to pressures of 45 GPa in a diamond anvil cell. The material was characterized by using X-ray, Raman, and optical transmission measurements. From single-crystal X-ray diffraction, the ambient structure of TNBA was determined which crystallizes in the P2 1 /c space group having four molecular units per unit cell. The X-ray data up to 9.2 GPa were fitted to a third-order Birch–Murnaghan equation of state by using the parameters K o = 13.2(2.4) GPa and K p = 5.1(1.4). Between 6.8 and 7.3 GPa, a phase transition was inferred in TNBA from concurrent fading of X-ray diffraction, disappearance of Raman peaks, increase in sample fluorescence, and discontinuous color change. The new phase was consistent with an amorphous state of at least partially intact molecules judging from the presence of higher-order Raman modes and irreversibility of the Raman spectra upon release. Piezochromism was observed with the translucent yellow TNBA gradually darkening and becoming opaque black at ~25 GPa. This correlated to the absorption edge gradually shifting to the red in the visible spectrum. Signs of two possible additional structural transitions were detected in the 32.4–41.0 GPa range as suggested by a jump in the absorption edge, the irreversible changes in the absorption spectrum upon release to ambient pressure, and by the lack of Raman modes in recovered samples. Here, the crystal and electronic structures of TNBA were also investigated up to 10 GPa by using DFT calculations and crystal structure prediction (CSP) simulations. In agreement with the experimentally observed transition at 7 GPa, the simulations at 10 GPa found a bevy of polymorphs lower in enthalpy and higher in density than P2 1 /c. The lowest calculated enthalpy structure was determined to be P2 1 2 1 2 1 , being in a different space group than the ambient experimental result.

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

High-Pressure Structural and Thermodynamic Properties of Cerium Orthosilicates (CeSiO 4 )

Pressure-induced phase transitions from the zircon structure-type (I4 1 /amd) to the scheelite structure type (I4 1 /a) are known for many ternary oxides systems (ABO 4 ). In this work, we present the first high-pressure study on synthetic stetindite (CeSiO 4 ) by a combination of in situ high-pressure synchrotron powder X-ray diffraction up to 36 GPa, implemented with and without dual sided laser heating, and in situ high-pressure Raman spectroscopy up to 43 GPa. Two phase transitions were identified: zircon to a high-pressure low-symmetry (HPLS) phase at 15 GPa and then to a scheelite at 18 GPa. The latter from HPLS scheelite phase was found irreversible; i.e., scheelite is fully quenchable at ambient conditions, as in other zircon-type phases. The bulk moduli (K 0 ) of stetindite, HPLS, and high-pressure scheelite phases were determined, respectively, as 171(5), 105(4), and 221(40) GPa by fitting to a second-order Birch-Murnaghan equation of state. The pressure derivatives of vibrational modes and Gru''neisen parameters of the zircon-structured polymorph are similar to those of other orthosilicate minerals. In conclusion, due to the larger ionic radii of Ce 4+ , with respect to Zr 4+ , stetindite was found to possess a softer bulk modulus and undergo the phase transitions at a lower pressure than zircon (ZrSiO 4 ), such observations are consistent with what were found in coffinite (USiO 4 ).

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

Phase Transition of Heptane under Pressure

A solid-to-solid phase transition was discovered in n-heptane with an onset at a pressure of 1.7 ± 0.1 GPa, evidenced by a discontinuity in the pressure–volume equation of state obtained from synchrotron X-ray diffraction experiments conducted at room temperature. The phase transition also resulted in discontinuous shifts of vibrational peaks in the infrared spectra. No additional phase transitions were detected up to 10 GPa, contradicting an earlier claim of a phase transition at 3 GPa. Furthermore, the hydrostatic limit of heptane in a diamond anvil cell was estimated to be 3.5 GPa. Above that pressure, the stress state is not well-defined, which may explain prior reports of a phase transition at 7.5 GPa.

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

Structure, Bonding, and Vibrational Dynamics of a Triamine High Energy Density Material under Pressure

High energy density materials have complex intermolecular interactions which influence their stability and performance. We used a combination of synchrotron X-ray diffraction, synchrotron infrared spectroscopy, and Raman vibrational spectroscopy, supplemented by density functional theory calculations, to probe pressure-induced changes in structure and intermolecular interactions of 1H,4'H-[3,3'-bis(1,2,4-triazole)]-4',5,5'-triamine as a model high energy density material up to 40 GPa. We find that compression of the triamine is accompanied by increased intermolecular interactions that give rise to an interesting evolution of the structure and bonding with pressure. Analysis of the equation of state determined from the X-ray diffraction indicates a change in compression mechanism near 19 GPa consistent with changes in vibrational spectra that provide evidence for a structural rearrangement associated with changes in hydrogen bonding near that pressure. As a result, the overall compressional behavior calculated theoretically agrees with that observed experimentally though differences are found that indicate the need for improved treatment of the intermolecular interactions including hydrogen bonding under pressure.

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

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2021

X-ray Free Electron Laser-Induced Synthesis of ε -Iron Nitride at High Pressures

The ultrafast synthesis of ε-Fe 3 N 1+x in a diamond-anvil cell (DAC) from Fe and N 2 under pressure was observed using serial exposures of an X-ray free electron laser (XFEL). When the sample at 5 GPa was irradiated by a pulse train separated by 443 ns, the estimated sample temperature at the delay time was above 1400 K, confirmed by in situ transformation of α- to γ-iron. Ultimately, the Fe and N 2 reacted uniformly throughout the beam path to form Fe 3 N 1.33 , as deduced from its established equation of state (EOS). We thus demonstrate that the activation energy provided by intense X-ray exposures in an XFEL can be coupled with the source time structure to enable exploration of the time-dependence of reactions under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multi-reward Reinforcement Learning Based Bond-Order Potential to Study Strain-Assisted Phase Transitions in Phosphorene

Here, we introduce a multi-reward reinforcement learning (RL) approach to train a flexible bond-order potential (BOP) for 2D phosphorene based on ab initio training data sets. Our approach is based on a continuous action space Monte Carlo tree search algorithm that is general and scalable and presents an efficient multiobjective optimization scheme for high-dimensional materials design problems. As a proof-of-concept, we deploy this scheme to parametrize multiple structural and dynamical properties of 2D phosphorene polymorphs. Our RL-trained BOP model adequately captures the structure, energetics, transformation barriers, equation of state, elastic constants, and phonon dispersions of various 2D P polymorphs. We use this model to probe the impact of temperature and strain rate on the phase transition from black (α-P) to blue phosphorene (β-P) through molecular dynamics simulations. A decrease in critical strain for this phase transition with increase in temperature is observed, and the underlying atomistic mechanisms are discussed.

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

Pressure-Dependent Phase Behavior of Chemically Synthesized Ytterbium Hydride Nanoparticles and Confinement in CMK-3

In this work, we developed a chemical synthesis protocol to prepare both spherical and rod shaped YbH x nanoparticles and infiltrate a mesoporous CMK-3 structure with YbH x nanoparticles. The phase evolution of the YbH x nanoparticles was monitored as a function of the pressure in a diamond anvil cell. The cell volume of the Yb atom showed an unexpected positive deviation from the calculated equation of state (EoS) for the α-phase (space group Pnma) of YbH x . Additionally, the β-phase (space group Fm3̅m) YbH x nanoparticles appeared to become amorphous at approximately 27 GPa, although the same behavior was not observed in the Pnma (α-phase). The Pmna phase remained stable (measured up to 59.3 GPa) and showed a significantly higher than expected volume per Yb atom relative to the calculated EoS. This deviation could be attributed to the nanoscale effect on the thermodynamics of ytterbium hydride formation, which could provide a path forward to accessing higher order ytterbium hydrides (YbH x , x > 3) at pressures much lower than in the bulk metal.

diamond anvil cell↗
DOE OSTI2020

Anisotropic Electron Fluid Closure Validated by in Situ Spacecraft Observations in the far Exhaust of Guide-field Reconnection

Anisotropic electron heating predictions of the equations of state of Le et al. (2009); https://doi.org/10.1103/PhysRevLett.102.085001 (“the EoS”) are applied to reconnection in the Earth's magnetosheath. The model is applicable to open systems where electrons are streaming along magnetic field lines into the reconnection region, sourced by fixed external reservoirs of plasma ambient to the reconnection region. While spacecraft observations have previously shown the EoS to hold in the region near the X-line, we find that for an event observed far downstream (~100 d i ) from the X-line both inflows and the exhaust follow the predictions of the EoS. Furthermore, the model underlying the EoS is extended to include additional perpendicular heating terms relevant to the considered event undergoing active magnetic compression with local electron trapping. Finally, the agreement between the spacecraft observations and the EoS emphasizes the roles of the ambient plasma sources and the dynamics of trapped and passing electrons in setting and controlling the anisotropic electron heating at large scale in naturally occurring plasma configurations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Thermal Pressure in the Laser‐Heated Diamond Anvil Cell: A Quantitative Study and Implications for the Density Versus Mineralogy Correlation of the Mantle

Abstract Thermal pressure is an inevitable thermodynamic consequence of heating a volumetrically constrained sample in the diamond anvil cell. Its possible influences on experimentally determined density‐mineralogy correlations are widely appreciated, yet the effect itself has never been experimentally measured. We present here the first quantitative measurements of the spatial distribution of thermal pressure in a laser‐heated diamond anvil cell (LHDAC) in both olivine and AgI. The observed thermal pressure is strongly localized and closely follows the distribution of the laser hotspot. The magnitude of the thermal pressure is of the order of the thermodynamic thermal pressure ( αK T Δ T ) with gradients between 0.5 and 1.0 GPa/10 μm. Remarkably, we measure a steep gradient in thermal pressure even in a sample that is heated close to its melting line. This generates consequences for pressure determinations in pressure‐volume‐temperature (PVT) equation of state measurements when using an LHDAC. We show that an incomplete account of thermal pressure in PVT experiments can lead to biases in the coveted depth versus mineralogy correlation. However, the ability to spatially resolve thermal pressure in an LHDAC opens avenues to measure difficult‐to‐constrain thermodynamic derivative properties, which are important for comprehensive thermodynamic descriptions of the interior of planets.

Yen, Connor Ethan↗
DOE OSTI2020

Spin Transitions and Compressibility of ε-Fe 7 N 3 and γ'-Fe 4 N: Implications for Iron Alloys in Terrestrial Planet Cores

Iron nitrides are possible constituents of the cores of Earth and other terrestrial planets. Pressure-induced magnetic changes in iron nitrides and effects on compressibility remain poorly understood. In this work, we report synchrotron X-ray emission spectroscopy (XES) and X-ray diffraction (XRD) results for ε-Fe 7 N 3 and γ'-Fe 4 N up to 60 GPa at 300 K. The XES spectra reveal completion of high- to low-spin transition in ε-Fe 7 N 3 and γ'-Fe 4 N at 43 and 34 GPa, respectively. The completion of the spin transition induces stiffening in bulk modulus of ε-Fe 7 N 3 by 22% at ~40 GPa, but has no resolvable effect on the compression behavior of γ'-Fe 4 N. Fitting pressure-volume data to the Birch-Murnaghan equation of state yields V 0 = 83.29 ± 0.03 (Å 3 ), K 0 = 232 ± 9 GPa, K 0 ' = 4.1 ± 0.5 for nonmagnetic ε-Fe 7 N 3 above the spin transition completion pressure, and V 0 = 54.82 ± 0.02 (Å 3 ), K 0 = 152 ± 2 GPa, K 0 ' = 4.0 ± 0.1 for γ'-Fe 4 N over the studied pressure range. By reexamining evidence for spin transition and effects on compressibility of other candidate components of terrestrial planet cores, Fe 3 S, Fe 3 P, Fe 7 C 3 , and Fe 3 C based on previous XES and XRD measurements, we located the completion of high- to low-spin transition at ~67, 38, 50, and 30 GPa at 300 K, respectively. The completion of spin transitions of Fe 3 S, Fe 3 P, and Fe 3 C induces elastic stiffening, whereas that of Fe 7 C 3 induces elastic softening. Changes in compressibility at completion of spin transitions in iron-light element alloys may influence the properties of Earth's and planetary cores.

58 GEOSCIENCES↗
DOE OSTI2023

Time‐Dependent Weakening of Granite at Hydrothermal Conditions

Abstract The evolution of a fault's frictional strength during the interseismic period is a critical component of the earthquake cycle, yet there have been relatively few studies that examine the time‐dependent evolution of strength at conditions representative of seismogenic depths. Using a simulated fault in Westerly granite, we examined how frictional strength evolves under hydrothermal conditions up to 250°C during slide‐hold‐slide experiments. At temperatures ≤100°C, frictional strength generally increases with hold duration but, at 200 and 250°C, an initial increase in strength transitions to rapid time‐dependent weakening for holds longer than 14 hr. Forward modeling of long hold periods at 250°C using the rate and state friction constitutive equations requires a second, strongly negative, state variable with a long evolution distance. This implies that significant hydrothermal alteration is occurring at 250°C, consistent with microstructural observations of dissolution and secondary mineral precipitation.

Geology↗