Lighthill stress tensor ($T_{ij}=\rho u_i{u}_j+(p^{\prime }-c_0^2{\rho }^{\prime }){\delta }_{ij}-{\sigma }_{ij}$)

This tensor contains convective momentum flux, departure from the reference linear pressure-density relation, and viscous stress. Its double divergence is an acoustic quadrupole distribution in Lighthill acoustic analogy.

Write a piecewise smooth Lighthill stress tensor as $T_{ij}=T_{ij}^{-}+D_{ij}H(S)$, where $D_{ij}=T_{ij}^{+}-T_{ij}^{-}$ and the regular level surface $S=0$ is a shock wave. The distributional derivative of the Heaviside step function givesHere repeated indices are summed, $S_i={\partial }_{i}S$ and $S_{ij}={\partial }_i{\partial }_{j}S$. Keep the coefficients as smooth extensions before multiplying distributions: prematurely replacing the coefficient of ${\delta }^{\prime }(S)$ by its surface trace can lose a $\delta (S)$ contribution. The extra surface delta distributions and their derivatives are supported entirely on the discontinuity.