用有限差分做磁流体模拟但是出现了异常的震荡
-
最近在研究用有限差分解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}
，时间上是一阶前差全隐求解，空间上用了二阶中心差分，尝试算了一下，在计算前期看似正常，但是后期出现了一个异常的震荡，最后导致整个计算都发散了，想请问一下这有没有可能是选取了太简单的离散格式导致了数值色散呢？
