# Shock Waves

Shocks are non-isentropic phenomena, i.e. entropy is not conserved, and necessitates a special formulation.

^{1}This value is added to the pressure and is computed by:

- l
- Is equal to $\sqrt[3]{\Omega}$ or to the characteristic length
- $\Omega $
- Volume
- $\frac{\partial {\epsilon}_{kk}}{\partial t}$
- Volumetric compression strain rate tensor
- $c$
- Speed of sound in the medium

- ${q}_{a}$ is a scalar factor on the quadratic viscosity to be adjusted so that the Hugoniot equations are verified. This value is defined by the user. The default value is 1.10.
- ${q}_{b}$ is a scalar factor on the linear viscosity that damps out the oscillations behind the shock. This is user specified. The default value is 0.05.

Default values are adapted for velocities lower than Mach 2. However, for
viscoelastic materials (LAW34, LAW35, LAW38) or honeycomb (LAW28), very small values are
recommended, that is, 10^{-20}.

^{1}Von Neumann J. and Richtmeyer R., “A method for the numerical calculation of hydrodynamical shocks”, Journal of applied physics, 1950.