Converging Shock Flows for a Mie-Gr\"uneisen Equation of State

TL;DR

This paper develops a scale-invariant approximation of the Mie-Gr"uneisen EOS to enable analytical shock solutions in 1D Euler flows, bridging a gap between idealized models and real solid material behavior.

Contribution

It introduces a scale-invariant surrogate for the Mie-Gr"uneisen EOS, allowing semi-analytical shock solutions for solid materials in Euler equations.

Findings

Constructed a scale-invariant adiabatic bulk modulus for Mie-Gr"uneisen EOS.

Enforced Rankine-Hugoniot conditions to determine model parameters.

Derived Guderley-like imploding shock solutions in a metal sphere.

Abstract

Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equations admit a wide variety of scale-invariant solutions (including the famous Noh, Sedov, and Guderley shock solutions) when the included equation of state (EOS) closure model assumes a certain scale-invariant form. However, this scale-invariant EOS class does not include even simple models used for shock compression of crystalline solids, including many broadly applicable representations of Mie-Gr\"uneisen EOS. Intuitively, this incompatibility naturally arises from the presence of multiple dimensional scales in the Mie-Gr\"uneisen EOS, which are otherwise absent from scale-invariant models that feature only dimensionless parameters (such as the adiabatic index in the ideal gas EOS). The current work extends previous efforts intended to rectify this inconsistency, by using a scale-invariant EOS model to approximate a Mie- Gr\"uneisen EOS form. To this end, the adiabatic bulk modulus for the Mie-Gr\"uneisen EOS is constructed, and its key features are used to motivate the selection of a scale-invariant approximation form. The remaining surrogate model parameters are selected through enforcement of the Rankine-Hugoniot jump conditions for an infinitely strong shock in a Mie-Gr\"uneisen material. Finally, the approximate EOS is used in conjunction with the 1D inviscid Euler equations to calculate a semi-analytical, Guderley-like imploding shock solution in a metal sphere, and to determine if and when the solution may be valid for the underlying Mie-Gr\"uneisen EOS.