Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic simulation”

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 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Development and application of a uranium mononitride (UN) potential: Thermomechanical properties and Xe diffusion

We report that atomic-scale modeling of thermophysical and defect properties of uranium mononitride (UN) plays an important role in establishing a better understanding and improved models of UN fuel performance. Having an accurate interatomic potential is crucial for generating reliable data at finite temperatures using molecular dynamic simulations. We report a new interatomic potential for UN, based on a combination of many-body and pairwise interactions, a simple form that we later show could be easily adapted to include Xe-U and Xe-N interactions, i.e., generating a U-N-Xe interatomic potential. The potential was fitted to experimental thermal expansion and single crystal elastic constants, as well as Frenkel, Schottky, anti-Schottky, and antisite pair reaction energies from density functional theory (DFT) calculations. Using the potential, we successfully reproduced experimental lattice parameters, thermal expansion, single crystal elastic constants, and temperature dependent heat capacity. The potential also performs reasonably well in reproducing the energy of the aforementioned stoichiometric defect reactions and defect migration barriers calculated using DFT. However, the potential underestimates the energy difference between the tetrahedral and dumbbell uranium interstitials, and a more complex potential form might be needed to overcome this issue. The potential was also used to predict UN single crystal elastic constants and elastic properties at different temperatures, showing that UN becomes softer and more compressible with increasing temperature. We also compare our potential against literature data from previous empirical potentials, demonstrating similar or better behavior depending on the property of interest. To enable the simulation of Xe in UN a Buckingham potential has been fitted to DFT-derived Xe incorporation energies. The potential was then used to determine the activation energy for Xe diffusion due to various Xe-containing defects, with {XeU : VU} exhibiting the lowest activation energy.

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

Accurate Fe–He machine learning potential for studying He effects in BCC-Fe

Nanostructured Ferritic Alloys (NFAs) containing oxide nanoparticles are prospective advanced structural materials in future fusion reactors. Lacking fusion neutron sources with an adequate flux, modeling is critical to understand radiation damage in the NFAs, including helium effects. Machine learning interatomic potentials (MLPs) have been demonstrated to be capable of describing atomic interactions in chemically complex systems with comparable accuracy to density functional theory (DFT) but with much less computational costs. Hence, MLPs are chosen to study He bubble accumulation in chemically complex NFAs. As a preliminary step, we develop a Fe–He potential, based on ~10,000 atomic configurations. The developed MLP accurately predicts the bulk properties of BCC-Fe compared to DFT, including the lattice constant (2.832 Å), elastic constants (c 11 = 271 GPa, c 12 = 141, c 44 = 93), and phonon frequencies (maximum error < 3.5%). The MLP predicts He binding energies in He n V bubbles and He n clusters with a mean absolute error (MAE) of only ~70 meV, which is ~3 to ~5 times smaller than the MAE of binding energies estimated using empirical potentials. The MLP also accurately predicts the migration barrier of interstitial He (64 meV). Furthermore, the MLP correctly predicts the maximum number of He atoms (critical size) at which a He n V bubble (He n cluster) can kick-out a self-interstitial atom (SIA), forming a bigger He n V 2 bubble (He n V bubble). Subsequently, we use the MLP to perform molecular dynamics simulations to investigate the effect of temperature on the critical size of He n V bubbles and He n clusters. We find that the critical size decreases with increasing temperature — indicating that smaller bubbles and clusters can kick-out an SIA at higher temperatures, thereby escalating the He bubble growth. We report the current Fe–He MLP can be further developed to include all chemical interactions in NFAs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An ab initio molecular dynamics study of varied compositions of the $\text{LiF-NaF-KF}$ molten salt

With increasing interest in molten salt reactors, there becomes a demand for investigations into thermophysical properties of salt systems. The LiF-NaF-KF (FLiNaK) salt system is a primary candidate for use in these reactors. However, the thermophysical properties of compositions outside the eutectic composition are still largely unknown. In this article, properties of ten unique compositions, including four ternary compositions, are investigated using ab initio molecular dynamics simulations across five temperatures between 900 K and 1300 K. The properties of interest are the density, thermal expansion, bulk modulus, compressibility, heat capacity, and enthalpy of mixing. In general, the results were found to be in good agreement with other literature and experimental results. The density and heat capacity had a tendency to be slightly underpredicted. No conclusions could be drawn about the bulk modulus and compressibility in terms of compositional dependence. The thermal expansion had a negative trend with respect to the LiF concentration and no trends were observed for the NaF or KF concentration. The enthalpy of mixing shows minima for the ternary compositions, with the near-equiatomic composition exhibiting the lowest values. Here this work shows the potential for compositional tailoring in the FLiNaK system to optimize thermophysical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

High-temperature structure, elasticity, and thermal expansion of ε-ZrH 1.8

Zirconium hydride is a promising candidate material for nuclear microreactor applications as a solid-state moderator component, owing to its favorable neutronics properties and good thermal stability over other metal hydrides. Here, in the present work, the crystal structure, thermal expansion, and elastic properties of the hydrogen-rich ε phase hydride were measured at elevated temperatures in the range 300–900 K. Samples were prepared by direct hydriding Zircaloy-4 metal – a nuclear-grade zirconium alloy. Room-temperature lattice parameters agree well with those reported from literature for unalloyed zirconium hydride and fall within an observed quadratic H-content dependence. The coefficients of thermal expansion, determined from lattice expansion and dilatometry, agree well within our work but were about 30 % lower than those reported by others for unalloyed hydrides. Density functional theory-based molecular dynamics simulations were used to compare with thermal expansion and elasticity measurements. Results showed lattice parameter temperature dependence and slope of thermal expansion align with those from measurements. Based on diffraction scans at select temperatures, ε phase remained stable in air up to at least 770 K. Likewise, dilatometry showed smooth thermal expansion up to the thermal decomposition temperature around 950 K. The precise decomposition temperature was not determined via diffraction due to sparse scanning. The complete elastic property measurements were gathered for ε-phase Ziracloy-4 hydride for the first time. Young's modulus was lower compared to the metal and δ hydride phases. High-temperature elasticity measurements were limited to <350 K due to acoustic dissipation effects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fast ionic conduction achieved through the design and synthesis of ceramic heterointerfaces

Lithium (Li) chloride and iron oxychloride (FeOCl), typically nonconductive, were combined to form a [Li 1+δ Cl] δ+ /[FeOCl] δ- heterointerface composite material (LFH), achieving ionic conductivities of >1 mS cm -1 . Analysis techniques (scanning transmission electron microscopy [STEM] and electron energy-loss spectroscopy [EELS]) indicated that the microstructure of LFH consisted of an amorphous LiCl-based shell surrounding a crystalline FeOCl-based core. Electrochemical measurements alongside solid-state 6,7 Li nuclear magnetic resonance (NMR) and molecular dynamic simulations revealed Li+ as the sole conductive species, with a diffusion barrier of ~0.25 eV. X-ray photoelectron spectroscopy (XPS) and X-ray absorption fine structure (XAFS) results further supported interstitial Li + diffusion at the heterointerface and within the LiCl phase, made possible by the heterointerface. Despite susceptibility to electronic conductivity, iron’s defects and multivalency (Fe³⁺, Fe²⁺) enable the Fe–O–Cl framework to accept Cl - , facilitating Li⁺ ionic conduction. Finally, a prototype solid-state cell (showing 97% Coulombic efficiency) demonstrated the viability of this heterointerface design for applications in energy storage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mapping of fracture and ionic conductivity changes in ion implanted solid electrolytes: Insights from molecular dynamics

Ion implantation emerges as a promising technique to address the persistent challenge of lithium (Li) filament growth in solid-state electrolytes as it can induce compressive stresses inhibiting crack growth and deflect dendrites, de facto mitigating early electrolyte failure. In this study, we examine the potential paradox of ion implantation: while aiming to enhance electrolyte performance, the radiation damage associated with implantation might inadvertently compromise both the ionic conductivity and the intrinsic fracture toughness of the material, rendering the material unsuitable for battery applications. Specifically, we employed molecular dynamics simulations to examine the scope of the downsides of ion implantation, specifically: (i) reduced ionic conductivity (due to radiation-induced amorphization) and (ii) mechanical stability (due to radiation-induced embrittlement) in ion-implanted Li 7 La 3 Zr 2 O 12 (LLZO) solid-state electrolytes. We explore how radiation damage impacts LLZO’s crystalline structure, Li-ion diffusion, and fracture properties at various temperatures and radiation damage levels. The study aims to provide insights into the competing effects of ion implantation and suggest potential engineering strategies for developing more robust solid-state electrolytes with improved conductivity and dendrite resistance.

Conductivity↗
DOE OSTI2022

Structure of maize BZR1-type β-amylase BAM8 provides new insights into its noncatalytic adaptation

Plant β-amylase (BAM) proteins play an essential role in growth, development, stress response, and hormone regulation. Despite their typical (β/α) 8 barrel structure as active catalysts in starch breakdown, catalytically inactive BAMs are implicated in diverse yet elusive functions in plants. Further, the noncatalytic BAM7/8 contain N-terminal BZR1 domains and were shown to be involved in the regulation of brassinosteroid signaling and possibly serve as sensors of yet an uncharacterized metabolic signal. While the structures of several catalytically active BAMs have been reported, structural characterization of the catalytically inactive BZR1-type BAMs remain unknown. Here, we determine the crystal structure of β-amylase domain of Zea mays BAM8/BES1/BZR1-5 and provide comprehensive insights into its noncatalytic adaptation. Using structural-guided comparison combined with biochemical analysis and molecular dynamics simulations, we revealed conformational changes in multiple distinct highly conserved regions resulting in rearrangement of the binding pocket. Altogether, this study adds a new layer of understanding to starch breakdown mechanism and elucidates the acquired adjustments of noncatalytic BZR1-type BAMs as putative regulatory domains and/or metabolic sensors in plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Targeted inhibition of MASTL kinase activity induces apoptosis in breast cancer

Microtubule-associated serine/threonine kinase-like (MASTL) (or Greatwall kinase (GWL)) is an important cell cycle regulating kinase that regulates the G2-M transition. Uncontrolled MASTL activity is implicated in breast cancer progression. To date, very few inhibitors have been reported against this protein. Here, structure-based computational modeling indicates that the natural product flavopiridol (FLV) binds strongly to MASTL and these results are validated using molecular dynamics simulation studies. Further, an in vitro kinase assay reveals an EC 50 (effective concentration) value of FLV to be 82.1 nM and a better IC 50 compared to the positive reference compound, staurosporine. FLV is found to inhibit MASTL kinase activity, arresting the cell growth in the G1 phase and inducing apoptosis in breast cancer cells. Consistent with these results differential gene expression obtained using RNA sequencing studies, and validated by RT PCR and immunoblot analysis, indicate that MASTL inhibition induces cell cycle arrest and apoptotic-related genes. Furthermore, metastasis- and inflammation- related genes are downregulated. Thus, the deregulation of MASTL signaling pathways on targeted inhibition of its kinase activity is revealed. This study lays a strong foundation for investigating FLV as a lead compound in breast cancer therapeutics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Effects of Fe atoms on hardening of a nickel matrix: Nanoindentation experiments and atom-scale numerical modeling

Significant hardening effect due to Fe concentrations in Ni-based alloys with face-centered-cubic structure has been studied by using a combined experimental and atomistic-based computational approach via nanoindentation tests. The obtained experimental load–displacement data for the [0 0 1] crystal orientation reached a qualitative good agreement with molecular dynamics simulations results, leading to strong evidence that the main strengthening factors are associated to sluggish dislocation diffusion, reduced defect sizes and the nucleation of tetrahedral stacking faults. Here, interstitial type prismatic dislocation loops mainly formed by $\frac{1}{6}$$\langle$1 1 2$\rangle$ Shockley dislocations are nucleated during the loading process, where their interaction leads to the formation of pyramidal shaped stacking fault which are mainly created by $\frac{1}{3}$$\langle$1 0 0$\rangle$ Hirth dislocations lines. Observing both types of defects coexisting in the same plastic deformation zone by both approaches. Reported mechanical data, measured experimentally and interpreted numerically, are also in accordance with microstructural SEM and TEM investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nanostructures evolution assessment and spectroscopic properties modification induced by electronic energy loss in KTaO 3 crystal

Modifications of micro/nanostructure and photoelectric properties under swift ion irradiation (645 MeV Xe 35+ ) of KTaO 3 crystal with varying electronic energy losses (7.2–31.4 keV/nm) and ion velocities (0.18–5.00 MeV/u) have been studied by combining experimental and calculated approaches. The i-TS calculations combined with molecular dynamics simulations are compared with the experimental observations, revealing the inner track fine structures from individual spherical defects to continuous ion tracks with core–shell morphologies, and a quantitative relationship, including the melting (0.42 eV/atom), damage (0.75 eV/atom), and amorphous (1.71 eV/atom) thresholds, is established to successfully predict the inner disorder/amorphous proportions and track damage morphologies. The surface nanostructures of nanohillocks (isolated or partial overlaps) and nanopits (serious overlaps) directly depend on the combined action of the deposited potential and kinetic energies of incident ions, which induce local melting and sublimation in the near-surface region. Owing to the decreased recombination and the increased separation efficiency between electrons and holes, the higher photogenerated charge carrier mobility enhanced the photocurrent effect, further optimizing the photoconductivity performance in Xe 35+ -irradiated KTaO 3 . Therefore, controlled defect engineering using the ion irradiation technique, as an effective strategy, could design tailored nanostructure systems and regulate photoelectric properties, further promoting the development of novel technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Defect annihilation in heavy ion irradiated polycrystalline gold

In this study, we explore the interaction of electron wind force (EWF) with defects originating from ion irradiation in-situ inside a transmission electron microscope. Nanocrystalline gold specimens were self-ion irradiated to a dose of 5 × 10 15 ions/cm 2 (45 displacement per atom) to generate a high density of displacement damage. We also developed a molecular dynamics simulation model to understand the associated atomic scale mechanisms. Both experiments and simulations show that the EWF can impart significant defect mobility even at low temperatures, resulting in the migration and elimination of defects in a few minutes. Overall, we propose that the EWF interacts with defects to create highly glissile Shockley partial dislocations, which makes the fast and low temperature defect annihilation possible.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhanced ionic conductivity and lack of paddle-wheel effect in pseudohalogen-substituted Li argyrodites

Superionic conductors are key to the development of safe and high-energy-density, all-solid-state batteries. Using a combined theoretical and experimental approach, we explore the feasibility of increasing the ionic conductivity through pseudohalogen substitution in the Li argyrodite structure. Under the guidance of calculated thermodynamic stability, BH 4 -substituted Li argyrodite, Li 5.91 PS 4.91 (BH 4 ) 1.09 , was successfully synthesized via a mechanochemical method. As-synthesized BH 4- substituted Li argyrodite displays an ionic conductivity of 4.8 mS/cm at 25°C. Ab initio molecular dynamics simulation trajectory analysis was used to investigate how BH 4 facilitates Li-ion diffusion and indicates only a weak correlation with the B–H bond motion. Here, we find that the enhanced conductivity mainly originates from the weak interaction between Li and BH 4 and find no evidence of a paddle-wheel effect from the polyanion. This work provides insight into how cluster ions enhance Li diffusion and systematically describes how to explore superionic conductors with pseudohalogen substitution.

25 ENERGY STORAGE↗
DOE OSTI2020

Electrolyte design for Li metal-free Li batteries

Li metal, with the lowest thermodynamically achievable negative electrochemical potential and the highest specific capacity (3860 mAh g -1 ), is the ultimate anode choice for Li batteries. However, the highest reported Li plating/stripping Coulombic efficiency (CE) of 99.5% after extensive efforts is still too low for the Li metal-free (all the Li metal in cycling comes from cathode, without anode pre-lithiation) Li metal batteries. The low CE is attributed to both non-uniform Li plating/stripping on the lithiophobic Cu current collector and Li dendrite growth through lithiophilic organic–inorganic solid electrolyte interphase (SEI) formed in carbonate electrolytes. Here, we use a lithiophilic Bismuth graphite blend (Bi–Gr) substrate to replace lithiophobic Cu current collector to seed a uniform Li nucleation, and form a lithiophobic LiF-rich SEI rather than lithiophilic organic-rich SEI to suppress Li dendrite growth. Molecular dynamics simulations reveal the preferential reduction of anions in 2.0 M LiPF 6 in tetrahydrofuran/2-methyl tetrahydrofuran (2.0 M LiPF 6 –mixTHF) electrolyte to generate LiF-rich SEI on plated Li. Bi–Gr substrate and 2.0 M LiPF 6 –mixTHF electrolyte enable the Li anodes to achieve a record high CE of 99.83% at a high capacity of 1.0 mAh cm -2 and current of 0.5 mA cm -2 . The Bi particles serve as dispersed nucleation centers that promote uniform Li deposition with strong adhesion to the substrate to avoid dead Li, while the lithiophobic LiF-rich SEI promotes lateral Li growth and suppresses the vertical Li dendrite growth even at a high current density of 3.0 mA cm -2 and high areal capacities of 3.0 mAh cm -2 . The regulation of Li nucleation and growth enables the Li metal-free LiFePO 4 full cells to achieve 100 cycles at a practical areal capacity of >2.0 mAh cm -2 . This article highlights the benefits of simultaneous substrate design to improve Li nucleation and electrolyte design to promote lithiophobic SEI growth, enabling a promising and practical route Li metal-free Li metal batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of static local distortions vs. dynamic motions on the stability and band gaps of cubic oxide and halide perovskites

Ternary ABX 3 perovskites made of corner-sharing BX 6 octahedra have long featured prominently in solid-state chemistry and condensed matter physics. Still, the joint understanding of their two main subgroups—halides and oxides—has not been fully developed. Indeed, unlike the case in simpler compounds having a single, robust repeated motif (“monomorphous”), certain cubic perovskites can manifest a non-thermal (= intrinsic) distribution of local motifs (“polymorphous networks”). Such static deformations can include positional degrees of freedom (e.g., atomic displacements and octahedral tilting) or magnetic moment degrees of freedom in paramagnets. Unlike thermal motion, such static distortions do not time-average to zero, being an expression of the intrinsic symmetry breaking preference of the chemical bonding. Here, the present study compares electronic structure features of oxide and halide perovskites starting from the static polymorphous distribution of motifs described by Density Functional Theory (DFT) minimization of the internal energy, continuing to finite temperature thermal disorder modeled via finite temperature DFT molecular dynamics. We find that (i) different oxide vs. halide ABX 3 compounds adopt different energy-lowering symmetry-breaking modes. The calculated pair distribution function (PDF) of SrTiO 3 from the first-principles agrees with recently measured PDF. (ii) In both oxides and halides, such static distortions lead to band gap blueshifts with respect to undistorted cubic Pm-3m structure. (iii) For oxide perovskites, high-temperature molecular dynamics simulations initiated from the statically distorted polymorphous structures reveal that the thermally-induced distortions can lead to a band gap redshift. (iv) In contrast, for cubic halide perovskite CsPbI 3 , both the intrinsic distortions and the thermal distortions contribute in tandem to band gap blueshift, the former, intrinsic effect being dominant. (v) In the oxide SrTiO 3 and CaTiO 3 (but not in halide) perovskites, octahedral tilting leads to the emergence of a distinct Γ–Γ direct band gap component as a secondary valley minimum to the well-known indirect R–Γ gap. Understanding such intrinsic vs. thermal effects on oxide vs. halide perovskites holds the potential for designing target electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗

Refine your search to explore more results.