Thermal conductivity tensor of β‐HMX as a function of pressure and temperature from equilibrium molecular dynamics simulations

Abstract

We apply the Green–Kubo (G–K) approach to obtain the thermal conductivity tensor of β‐1,3,5,7‐tetranitro‐1,3,5,7‐tetrazoctane (β‐HMX) as a function of pressure and temperature from equilibrium molecular dynamics (MD) simulations. Direct application of the G–K formula exhibits slow convergence of the integrated thermal conductivity values even for long (120 ns) accumulated simulation times. To partially mitigate this slow convergence, we apply a recently implemented numerical procedure that involves physically justified filtering of the MD‐calculated G–K heat current and fitting the integrated time‐dependent thermal conductivity to a physically motivated double‐exponential function. The thermal conductivity tensor was determined for pressures 0.1 MPa 30 GPa and temperatures 300 K 900 K. The thermal conductivity increases with increasing pressure, by approximately an order of magnitude over the interval considered, and decreases with increasing temperature. The MD predictions are compared to experimental and other theoretically determined values for the thermal conductivity of β‐HMX. A simple, semi‐empirical fitting form is proposed that captures the behavior of over the pressure and temperature intervals studied.

Document Details

Document Type
Pub Defense Publication
Publication Date
May 23, 2023
Source ID
10.1002/prep.202200345

Entities

People

  • Andrey Pereverzev
  • Tommy Sewell

Organizations

  • Air Force Office of Scientific Research
  • University of Missouri

Tags

Fields of Study

  • Physics

Readers

  • Calculus or Mathematical Analysis
  • Computational Modeling and Simulation
  • Thermal Physics or Thermal Science.