用有限差分做磁流体模拟但是出现了异常的震荡

[nothingts]

最近在研究用有限差分解MHD（磁流体）方程，贴一下公式：
\begin{equation}
\frac{\partial\rho}{\partial t}+\nabla\cdot\left[\left(\rho_0+\rho\right){v}\right]-\nabla\cdot\left(d\nabla\rho\right)=0\
\frac{\partial{v}}{\partial t}+\nabla{{v}}\cdot{{v}}-\frac{1}{\rho+\rho_0}\left[{j}_0\times{B}+\left(\nabla\times{B}\right)\times\left({B}+{B}_0\right)\right]-\nabla\cdot\left(\nu\nabla{v}\right)=0\
\frac{\partial{B}}{\partial t}-\nabla\times\left({v}\times{B}\right)-\eta\nabla^2{{B}}=0
\end{equation}
，时间上是一阶前差全隐求解，空间上用了二阶中心差分，尝试算了一下，在计算前期看似正常，但是后期出现了一个异常的震荡，最后导致整个计算都发散了，想请问一下这有没有可能是选取了太简单的离散格式导致了数值色散呢？

[nothingts]

好像不能发动图，我试试截几个图吧（前两张是计算“正常”时，后面两张是出现了震荡和最后崩溃时的情况）