20 Search Results

Using reverse non-equilibrium molecular dynamics simulations, we have determined the dependences on temperature and pressure of the thermal conductivity tensors for the monoclinic γ and ɛ polymorphs of hexanitrohexaazaisowurtzitane (HNIW or CL20). A recently developed non-reactive force field [X. Bidault and S. Chaudhuri, RSC Adv. 9, 39649–39661 (2019)], designed to study polymorphism and phase transitions in CL20, is employed. The effects of temperature and pressure are investigated between 200 and 500 K and up to 0.5 GPa for γ-CL20 and 2 GPa for ɛ-CL20. In order to obtain the full thermal conductivity tensor, κ _ij , for the monoclinic crystals, four distinct heatmore » propagation directions are used. We find that κ _ij for both polymorphs is more isotropic than for other energetic molecular crystals, including α- and γ-RDX, β-HMX, and PETN, with a maximum difference of 9.8% between orientations observed at 300 K and 0 GPa for γ-CL20 and a maximum difference of 4.8% for ɛ-CL20. The average thermal conductivity, [Formula: see text], of ɛ-CL20 is 6.4% larger than that of γ-CL20 at 300 K and 0 GPa. Analytic linear functions of the inverse temperature and the pressure are provided, which fit the data well and can be used to predict the thermal conductivity of both polymorphs for any orientation, pressure, and temperature in and around the fitting range. Our predictions agree reasonably well with the limited available experimental data, for which the polymorph type is unknown.« less

Energetic materials undergo hundreds of chemical reactions during exothermic runaway, generally beginning with the breaking of the weakest chemical bond, the “trigger linkage.” Herein we report the syntheses of a series of pentaerythritol tetranitrate (PETN) derivatives in which the energetic nitrate ester groups are systematically substituted by hydroxyl groups. Because all the PETN derivatives have the same nitrate ester-based trigger linkages, quantum molecular dynamics (QMD) simulations show very similar Arrhenius kinetics for the first reactions. However, handling sensitivity testing conducted using drop weight impact indicates that sensitivity decreases precipitously as nitrate esters are replaced by hydroxyl groups. These experimental resultsmore » are supported by QMD simulations that show systematic decreases in the final temperatures of the products and the energy release as the nitrate ester functional groups are removed. To better interpret these results, we derive a simple model based only on the specific enthalpy of explosion and the kinetics of trigger linkage rupture that accounts qualitatively for the decrease in sensitivity as nitrate ester groups are removed.« less

We use reverse non-equilibrium molecular dynamics to determine the thermal conductivity tensor, κ, of tetragonal pentaerythritol tetranitrate (PETN). The most stable form under standard ambient conditions (PETN I, with space group P4¯2₁c) is considered in the temperature and pressure intervals 200–500 K and 0–8 GPa, which covers approximately the stability range for this polymorph. We compute κ along the [100] and [001] directions, which is enough to construct the full thermal conductivity tensor for a system with tetragonal symmetry. In addition, we also determine κ along the [010] direction and confirm that κ₁₀₀ ~ κ₀₁₀, with less than 4% averagemore » absolute error between the two quantities. We observe an anisotropic response for κ, with κ₁₀₀ > κ₀₀₁ across the whole (T,P) interval, and 37% difference at 300 K, 0 GPa. Furthermore, we provide analytical functions to interpolate κ(P,T) within the fitting interval and serve as input for continuum-scale simulations.« less

We have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Ourmore » results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.« less

Atomistic simulations reveal increased cation inversion at grain boundaries in spinel. As the grain size is reduced, the apparent level of inversion in the material will increase as the grain boundaries become an increasing fraction of the material.

Diffusion of fission gas in UO₂ nuclear fuel impacts several important performance metrics, such as fission gas release, swelling, and thermal conductivity. Current empirical models of fission gas release have significant uncertainty, some of which derives from the bulk diffusion rate and its dependence on, for example, fuel chemistry and irradiation. In this work, we have applied the previously-developed Free Energy Cluster Dynamics (FECD) methodology in the code Centipede to calculate xenon cluster concentrations in UO₂ under intrinsic (high temperature) and irradiation-enhanced (intermediate temperature) conditions in order to develop a model of the xenon diffusion coefficient based on the atomicmore » scale mechanisms responsible for transport. While the diffusion mechanism for xenon in UO₂ is adequately described by the Xe + U₂O vacancy cluster for intrinsic conditions, a similar process is not capable of capturing measured in-pile fission gas diffusivity at intermediate temperatures. Therefore, a different diffusion mechanism must dominate under this regime. Using calculated atomistic data, we have shown that irradiation-enhanced diffusion at intermediate temperatures occurs via the larger Xe + U₄O_y vacancy clusters, which have lower migration barriers and increase in concentration by several orders of magnitude compared to intrinsic conditions. This mechanism is enabled by the increased uranium vacancy concentration under irradiation due to Frenkel pair production. In addition, the fast migration of uranium interstitials with two attached oxygen interstitials lowers the total uranium interstitial concentration through reactions with sinks. This allows the extended defects, such as Xe + U₄O_y vacancy clusters, to maintain high concentrations by limiting annihilation with attached vacancies. Furthermore, predictions using the Xe + U₄O_y diffusion mechanism are in good agreement with experiment, albeit with some differences in the Arrhenius slope, which we believe may be related to either experimental or model parameter uncertainty. Lastly, an analytical expression suitable for application in fuel performance simulations was derived to capture the predictions of the Centipede simulations.« less

A new density functional theory, average atom based model for the electrical conductivity of dense plasmas with a mixture of ion species, containing no adjustable parameters, is presented herein. The model takes the temperature, mass density and relative abundances of the species as input. It takes into account partial ionization, ionic structure, and core-valence orthogonality, and uses quantum mechanical calculations of cross sections. Comparison to an existing high fidelity but computationally expensive method reveals good agreement. The new model is computationally efficient and can reach high temperatures. A new mixing rule is also presented that gives reasonably accurate conductivities formore » high temperature plasma mixtures.« less

A new Sesame-type table for the electrical conductivity of aluminium is dicussed. The table is built upon density functional theory calculations and ranges from 10^–3 to 1 times solid density (2.7 g/cm3), and from 10^–2 to 10³ eV in temperature. The table is compared with other those of simulations and to experiments and is generally in good agreement. The high-temperature, classical limit of the conductivity is recovered for the highest temperatures and lowest densities. The table is critically evaluated, and directions for improvements are discussed.

While fission fragments pass through UO₂ nuclear fuel, a considerable concentration of Frenkel pair defects (i.e. vacancies and interstitials) are created. The steady-state concentration of these defects leads to enhanced uranium self-diffusion, one of several fundamental kinetic parameters that control key engineering properties such as creep and fission gas swelling in UO₂ nuclear fuel. A cluster dynamics method to track point defects and defect clusters has been implemented in the MARMOT phase-field code in order to predict as-measured out-of-pile and irradiation enhanced thermal diffusivity. Here, the calculated uranium self-diffusion coefficient compares well with non-irradiated fuel measurements, and shows similar trendsmore » to those observed in irradiated fuel, which is a good result given the complexities introduced by non-stoichiometric compositions.« less

Using temperature accelerated dynamics, an accelerated molecular dynamics method, we examine the relationship between composition and cation ordering and defect transport in the mixed pyrochlore Gd₂(Ti_1-xZr_x)₂O₇, using the oxygen vacancy as a representative defect structure. In this study, we find that the nature of transport is very sensitive to the cation structure, transitioning, as a function of composition, from three-dimensional migration to two-dimensional to pseudo-one-dimensional to becoming essentially immobile before reverting back to three-dimensional as the Zr content is increased. The rates of migration are also affected by the cation structure in the various compositions. This behavior is driven bymore » the connectivity of Ti polyhedra in the material, with more extensive networks of Ti ions leading to a greater ability of the vacancy to traverse the material. In conclusion, our results indicate that the nature of transport is dictated by the cation structure of the material and that, conversely, the cation structure could be used to control transport and potentially other functionalities in mixed pyrochlores.« less