Numerical Method for Shock Front Hugoniot States

TL;DR

This paper introduces a novel Continuous Hugoniot Method for efficiently simulating shock wave fronts and mapping Hugoniot states as a continuous function of shock strength, demonstrated on Lennard-Jonesium.

Contribution

The paper presents the first continuous mapping of Hugoniot states using a new efficient simulation method for shock fronts.

Findings

Good agreement with prior simulations

Efficient generation of shock states

Applicable to Lennard-Jonesium along <100> direction

Abstract

We describe a Continuous Hugoniot Method for the efficient simulation of shock wave fronts. This approach achieves significantly improved efficiency when the generation of a tightly spaced collection of individual steady-state shock front states is desired, and allows for the study of shocks as a function of a continuous shock strength parameter, $v_p$. This is, to our knowledge, the first attempt to map the Hugoniot continuously. We apply the method to shock waves in Lennard-Jonesium along the $<100>$ direction. We obtain very good agreement with prior simulations, as well as our own benchmark comparison runs.