Search NASA⌕ Search

Engineering topics

Myint, Philip C.

Publications and source records attributed to Myint, Philip C..

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fluid-phase helium: Shock-compression experiments, quantum molecular dynamics simulations, and development of an equation of state

Helium (He) plays a critical role in numerous areas ranging from the study of celestial objects like brown dwarfs and gas giants to modern-day technologies like nuclear energy and rocket propulsion. For many of these applications, it is essential to have a reliable equation of state (EOS) for He that yields an accurate representation of its thermodynamic behavior. To help constrain and develop such EOS models, we have performed a series of shock-compression experiments on cryogenic liquid He to pressures exceeding 100 GPa using a magnetically accelerated flyer plate on Sandia National Laboratories' Z-machine. We have also performed quantum molecular dynamics simulations that are consistent with our shock measurements. None of the previously available EOSs agree with our experimental and simulation results, motivating the development of a fluid-phase He EOS that we present in this study. Here, we show that our EOS yields good agreement with published data that span temperatures and pressures encountered across a diverse array of applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Model for the Solid–Liquid Interfacial Free Energy at High Pressures

The free energy involved in the formation of an interface between two phases (e.g., a solid–liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid–liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or over-pressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid–liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Measuring the melting curve of iron at super-Earth core conditions

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth’s interior. Here, we used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth’s inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth’s mass have the longest dynamos, which provide important shielding against cosmic radiation.

58 GEOSCIENCES↗
DOE OSTI2020

Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes

In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.

42 ENGINEERING↗
DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Efficient and robust phase-split computations in the internal energy, volume, and moles ( UVN ) space

Phase-split computations in an isolated system, which in general may be defined as one where the total internal energy (U), volume (V), and the number of moles (N = N 1 , N 2 ,..., N n ) of the components are fixed at some specified set of values, involves the determination of the temperature (T), pressure (P), plus the amount and composition of the various phases that constitute the system. A simpler, but analogous problem is one where T, P, and N are specified instead. In TPN space, one may first perform a stability analysis to determine whether the system is stable, meaning whether at equilibrium it will split up into multiple phases or remain in single phase. If the single-phase state is unstable, the stability analysis reliably provides a good set of initial guesses in the subsequent phase-split computations. In UVN space, however, we demonstrate that the stability analysis (which is the main subject of our earlier study [1]) may not in general provide good enough initial guesses; we offer alternative strategies for setting up good initial guesses. Furthermore, we show that a combination of successive substitution iteration (SSI) and Newton's method---two iterative methods that are prevalent in the literature in TPN space---facilitates a robust and efficient algorithm for phase-split computations in isolated systems. Finally, this combination has so far not been applied in UVN space.

02 PETROLEUM↗
DOE OSTI2019

Equation of state for a chemically dissociative, polyatomic system: Carbon dioxide

A notorious challenge in high-pressure science is to develop an equation of state (EOS) that explicitly treats chemical reactions. For instance, many materials tend to dissociate at high pressures and temperatures where the chemical bonds that hold them together break down. We present an EOS for carbon dioxide (CO 2 ) that allows for dissociation and captures the key material behavior in a wide range of pressure– temperature conditions. Carbon dioxide is an ideal prototype for the development of a wide-ranging EOS that allows for chemical-dissociation equilibria since it is one of the simplest polyatomic systems and because it is of great interest in planetary science and in the study of detonations. Here, we show that taking dissociation into account significantly improves the accuracy of the resulting EOS compared to other EOSs that either neglect chemistry completely or treat CO 2 dissociation in a more rudimentary way.

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

Numerical modeling of solid-cluster evolution applied to the nanosecond solidification of water near the metastable limit

Classical nucleation theory (CNT) is a promising way to predictively model the submicrosecond kinetics of phase transitions that occur under dynamic compression, such as the suite of experiments performed over the past two decades on the solidification of liquid water to the high-pressure ice VII phase. Myint et al. [Phys. Rev. Lett. 121, 155701 (2018)] presented the first CNT-based model for these types of rapid phase transitions, but relied on an empirical scaling parameter in their transient induction model to simulate the lag time that occurs prior to the onset of significant formation of ice VII clusters in the system. To build on that study, we model the liquid water–ice VII phase transformation using a numerical discretization scheme to solve the Zel’dovich–Frenkel partial differential equation, which is a fundamental CNT-based kinetic equation that describes the statistical time-dependent behavior of solid cluster formation. The Zel’dovich–Frenkel equation inherently accounts for transience in the nucleation kinetics and eliminates the need for the empirical scaling factor used by Myint et al. One major result of this research is that transience is found to play a relatively small role in the nucleation process for the dynamic-compression time scales of the liquid water–ice VII experiments being simulated. Instead, we show that it is possible to properly model the lag time using steady-state CNT by making small refinements to the interfacial free energy value. In conclusion, we have also developed a new dimensionless parameter that may be applied a priori to predict whether or not transient nucleation will be important in a given dynamic-compression experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Two-phase equation of state for lithium fluoride

Here, we present an equation of state for the solid and liquid phases of lithium fluoride that covers a wide range of conditions from ambient pressure and temperature to the high pressures and temperatures exhibited in shock- and ramp-compression studies. The particular solid phase we have focused on in this work is the B1 phase. We have followed an approach where the pressure and heat-capacity functions of both phases are fit to experimental data and our own quantum molecular dynamics simulations and are then integrated in a thermodynamically consistent way to obtain the corresponding free-energy functions. This approach yields a two-phase equation of state that provides better overall agreement with experimental data than other equations of state for lithium fluoride, such as SESAME 7271v3, LEOS 2240, and the model presented by Smirnov. The last of these is a three-phase equation of state that predicts a B1–B2 transition along the shock Hugoniot at a pressure of about 140 GPa. This solid–solid transition has been a topic of speculation and debate in the literature for over 50 years, culminating in the work of Smirnov, who has developed the only potentially viable equation of state that allows for this transition. We explain why the proposed B1–B2 transition at 140 GPa is not consistent with recent velocimetry data.

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

Nanosecond Freezing of Water at High Pressures: Nucleation and Growth near the Metastability Limit

The fundamental study of phase transition kinetics has motivated experimental methods toward achieving the largest degree of undercooling possible, more recently culminating in the technique of rapid, quasi-isentropic compression. This approach has been demonstrated to freeze water into the high-pressure ice VII phase on nanosecond timescales, with some experiments undergoing heterogeneous nucleation while others, in apparent contradiction, suggest a homogeneous nucleation mode. We show through a combination of theory, simulation, and analysis of experiments that these seemingly contradictory results are in agreement when viewed from the perspective of classical nucleation theory. We find that, perhaps surprisingly, classical nucleation theory is capable of accurately predicting the solidification kinetics of ice VII formation under an extremely high driving force ( | Δ μ / k B T | ≈ 1 ) but only if amended by two important considerations: (i) transient nucleation and (ii) separate liquid and solid temperatures. Finally, this is the first demonstration of a model that is able to reproduce the experimentally observed rapid freezing kinetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Rapid freezing of water under dynamic compression

Understanding the behavior of materials at extreme pressures is a central issue in fields like aerodynamics, astronomy, and geology, as well as for advancing technological grand challenges such as inertial confinement fusion. Dynamic compression experiments to probe high-pressure states often encounter rapid phase transitions that may cause the materials to behave in unexpected ways, and understanding the kinetics of these phase transitions remains an area of great interest. We examine experimental and theoretical/computational efforts to study the freezing kinetics of water to a high-pressure solid phase known as ice VII. We first present a detailed analysis of dynamic compression experiments in which water has been observed to freeze on sub-microsecond time scales to ice VII. This is followed by a discussion of the limitations of currently available molecular and continuum simulation methods in modeling these experiments. We then describe how our phase transition kinetics models, which are based on classical nucleation theory, provide a more physics-based framework that overcomes some of these limitations. We give suggestions on future experimental and modeling work on the liquid–ice VII transition, including an outline of the development of a predictive multiscale model in which molecular and continuum simulations are intimately coupled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Free energy models for ice VII and liquid water derived from pressure, entropy, and heat capacity relations

Here, we present equations of state relevant to conditions encountered in ramp and multiple-shock compression experiments of water. These experiments compress water from ambient conditions to pressures as high as about 14 GPa and temperatures of up to several hundreds of Kelvin. Water may freeze into ice VII during this process. Although there are several studies on the thermodynamic properties of ice VII, an accurate and analytic free energy model from which all other properties may be derived in a thermodynamically consistent manner has not been previously determined. We have developed such a free energy model for ice VII that is calibrated with pressure-volume-temperature measurements and melt curve data. Furthermore, we show that liquid water in the pressure and temperature range of interest is well-represented by a simple Mie-Grüneisen equation of state. Our liquid water and ice VII equations of state are validated by comparing to sound speed and Hugoniot data. Although they are targeted towards ramp and multiple-shock compression experiments, we demonstrate that our equations of state also behave reasonably well at pressures and temperatures that lie somewhat beyond those found in the experiments.

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

Thermodynamics of HMX Polymorphs and HMX/RDX Mixtures

In this paper, we present thermodynamic models for the five most commonly studied phases of the energetic material octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX): liquid HMX and four solid polymorphs (α-, β-, γ-, and δ-HMX). We show results for the density, heat capacity, bulk modulus, and sound speed, as well as a phase diagram that illustrates the temperature and pressure regions over which the various HMX phases are most thermodynamically stable. The models are based on the same equation of state presented in our recently published paper [Myint et al., Ind. Eng. Chem. Res., 2016, 55, 2252] on another energetic material, hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). We combine our HMX and RDX models together so that the equation of state can also be applied to liquid and solid mixtures of HMX/RDX. This allows us to generate an HMX/RDX phase diagram and calculate the enthalpy change associated with a few different kinds of phase transitions that these mixtures may undergo. Our paper is the first to present a single equation of state that is capable of modeling both pure HMX and HMX/RDX mixtures. A distinct feature of HMX is the strongly metastable nature of its polymorphs. This has caused some ambiguity in the literature regarding the thermodynamic stability of α-HMX. Finally, by examining possible arrangements for the relative order of the six different solid-solid transition (α–β, α–γ, α–δ, β–γ, β–δ, and γ–δ) temperatures, we conclude that α-HMX must be thermodynamically stable so that the HMX phase diagram must have an α phase region.

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

The CPA Equation of State and an Activity Coefficient Model for Accurate Molar Enthalpy Calculations of Mixtures with Carbon Dioxide and Water/Brine

Thermodynamic property calculations of mixtures containing carbon dioxide (CO 2 ) and water, including brines, are essential in theoretical models of many natural and industrial processes. The properties of greatest practical interest are density, solubility, and enthalpy. Many models for density and solubility calculations have been presented in the literature, but there exists only one study, by Spycher and Pruess, that has compared theoretical molar enthalpy predictions with experimental data. In this report, we recommend two different models for enthalpy calculations: the CPA equation of state by Li and Firoozabadi, and the CO 2 activity coefficient model by Duan and Sun. We show that the CPA equation of state, which has been demonstrated to provide good agreement with density and solubility data, also accurately calculates molar enthalpies of pure CO 2 , pure water, and both CO 2 -rich and aqueous (H 2 O-rich) mixtures of the two species. It is applicable to a wider range of conditions than the Spycher and Pruess model. In aqueous sodium chloride (NaCl) mixtures, we show that Duan and Sun’s model yields accurate results for the partial molar enthalpy of CO 2 . It can be combined with another model for the brine enthalpy to calculate the molar enthalpy of H 2 O-CO 2 -NaCl mixtures. We conclude by explaining how the CPA equation of state may be modified to further improve agreement with experiments. This generalized CPA is the basis of our future work on this topic.

02 PETROLEUM↗
DOE OSTI2015

Thin liquid films in improved oil recovery from low-salinity brine

Low-salinity waterflooding is a relatively new method for improved oil recovery that has generated much interest. It is generally believed that low-salinity brine alters the wettability of oil reservoir rocks towards a wetting state that is optimal for recovery. The mechanism(s) by which the wettability alteration occurs is currently an unsettled issue. This study reviews recent studies on wettability alteration mechanisms that affect the interactions between the brine/oil and brine/rock interfaces of thin brine films that wet the surface of reservoir rocks. Of these mechanisms, we pay particular attention to double-layer expansion, which is closely tied to an increase in the thickness and stability of the thin brine films. Our review examines studies on both sandstones and carbonate rocks. We conclude that the thin-brine-film mechanisms provide a good qualitative, though incomplete, picture of this very complicated problem. Finally, we give suggestions for future studies that may help provide a more quantitative and complete understanding of low-salinity waterflooding.

02 PETROLEUM↗