中文

精确保持 $\nabla \cdot \boldsymbol{B} = 0$ 的 MHD 模型稳定有限元方法

数值分析 2014-10-07 v1

摘要

本文致力于设计和分析磁流体动力学 (MHD) 系统的若干保结构有限元格式。该方法的主要特征是自然地保持了重要的高斯定律,即 B=0\nabla\cdot\boldsymbol{B}=0。与大多数消除电场变量 E\boldsymbol{E} 并直接离散化磁场的现有方法不同,我们的新方法使用 H(curl)H(\mathrm{curl}) 的 N\'{e}d\'{e}lec 型边单元离散化电场 E\boldsymbol{E},同时使用 H(div)H(\mathrm{div}) 的 Raviart-Thomas 型面单元离散化磁场 B\boldsymbol{B}。因此,磁场的无散度条件在离散层面上精确成立。对于这种新的有限元方法,可以以与连续情形类似的方式自然地建立能量稳定性估计。此外,通过对连续和离散情形均使用 Brezzi 理论,严格建立了完全非线性系统的 Picard 和 Newton 线性化的适定性。这种适定性自然地导出了线性化系统的鲁棒(且最优)预条件子。

关键词

引用

@article{arxiv.1410.1095,
  title  = {Stable Finite Element Methods Preserving $\nabla \cdot \boldsymbol{B} = 0$ Exactly for MHD Models},
  author = {Kaibo Hu and Yicong Ma and Jinchao Xu},
  journal= {arXiv preprint arXiv:1410.1095},
  year   = {2014}
}