An asymptotic preserving method for strongly anisotropic diffusion equations based on field line integration
详细信息 查看全文
- 作者：Min Tang ; tangmin@sjtu.edu.cn ; Yihong Wang
- 关键词：Anisotropic diffusion ; Asymptotic preserving ; Uniform convergence ; Field line integration
- 刊名：Journal of Computational Physics
- 出版年：2017
- 出版时间：1 February 2017
- 年：2017
- 卷：330
- 期：Complete
- 页码：735-748
- 全文大小：2978 K
- 卷排序：330
文摘
In magnetized plasma, the magnetic field confines the particles around the field lines. The anisotropy intensity in the viscosity and heat conduction may reach the order of 1012. When the boundary conditions are periodic or Neumann, the strong diffusion leads to an ill-posed limiting problem. To remove the ill-conditionedness in the highly anisotropic diffusion equations, we introduce a simple but very efficient asymptotic preserving reformulation in this paper. The key idea is that, instead of discretizing the Neumann boundary conditions locally, we replace one of the Neumann boundary condition by the integration of the original problem along the field line, the singular 1/ϵ1/ϵ terms can be replaced by O(1)O(1) terms after the integration, which yields a well-posed problem. Small modifications to the original code are required and no change of coordinates nor mesh adaptation are needed. Uniform convergence with respect to the anisotropy strength 1/ϵ1/ϵ can be observed numerically and the condition number does not scale with the anisotropy.