Engineering topics

Perriot, Romain Thibault

Publications and source records attributed to Perriot, Romain Thibault.

4 records

DOE OSTI2022

Thermal conductivity tensor of γ and ε -hexanitrohexaazaisowurtzitane as a function of pressure and temperature

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). In this work, 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 heat 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, $\bar{κ}$, 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.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Pressure, temperature, and orientation dependent thermal conductivity of pentaerythritol tetranitrate (PETN)

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 1 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 κ 100 ~ κ 010 , with less than 4% average absolute error between the two quantities. We observe an anisotropic response for κ, with κ 100 > κ 001 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.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, 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. Our 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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Model for the electrical conductivity in dense plasma mixtures

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 for high temperature plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗