Shocks are non-isentropic phenomena, that is, entropy is not conserved, and
necessitates a special formulation.
The missing energy is generated by an artificial
bulk viscosity as derived by von Neumann and Richtmeyer. 1 This
value is added to the pressure and is computed by:図 1.
Where,
l
Is equal to or to the characteristic length
Volume
Volumetric compression strain rate tensor
Speed of sound in the medium
The values of and are adimensional scalar factors defined as:
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.
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.
1Von Neumann J. and Richtmeyer R., A method for the numerical
calculation of hydrodynamical shocks, Journal of applied physics,
1950.