CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    Algorithm
    1
    2
    926
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • N
      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}
      ,时间上是一阶前差全隐求解,空间上用了二阶中心差分,尝试算了一下,在计算前期看似正常,但是后期出现了一个异常的震荡,最后导致整个计算都发散了,想请问一下这有没有可能是选取了太简单的离散格式导致了数值色散呢?
      diver.gif

      1 条回复 最后回复 回复 引用
      • N
        nothingts 最后由 编辑

        好像不能发动图,我试试截几个图吧(前两张是计算“正常”时,后面两张是出现了震荡和最后崩溃时的情况)
        diver1.jpg diver2.jpg diver3.jpg diver4.jpg

        1 条回复 最后回复 回复 引用
        • First post
          Last post