我正在使用以下等式计算冲击波的熵产生： 和 我可以访问空间和时间网格

$$\begin{equation} \sigma = J'_q\frac{\partial}{\partial x}\left(\frac{1}{T}\right) +\frac{1}{T}\frac{4\eta}{3}\left(\frac{\partial v}{\partial x}\right)^2 \end{equation}$$ $$\begin{equation} \sigma = \frac{\partial \rho_s}{\partial t} +\frac{\partial J_s}{\partial x} \end{equation}$$ $[x_i,...x_N]$$[t_i,...t_K]$. 目前我在空间和时间网格上都使用了中心离散化。该方法适用于第一个方程，但对于第二个方程，冲击波前后似乎存在振荡。我可以使用更好的离散化方法吗？我应该在时间维度上使用不同的方法吗？