A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions
详细信息 查看全文
- 作者：L. Chacó ; n ; chacon@lanl.gov" class="auth_mail" title="E-mail the corresponding author ; G. Chen
- 关键词：Conservative discretization ; Darwin model ; Fully implicit algorithms ; Curvilinear meshes
- 刊名：Journal of Computational Physics
- 出版年：2016
- 出版时间：1 July 2016
- 年：2016
- 卷：316
- 期：Complete
- 页码：578-597
- 全文大小：644 K
文摘
We extend a recently proposed fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions (Chen and Chacón (2015) [1]) to curvilinear geometry. As in the Cartesian case, the approach is based on a potential formulation (ϕ, A), and overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. Conservation theorems for local charge and global energy are derived in curvilinear representation, and then enforced discretely by a careful choice of the discretization of field and particle equations. Additionally, the algorithm conserves canonical-momentum in any ignorable direction, and preserves the Coulomb gauge ∇⋅A=0 exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. We demonstrate the accuracy and efficiency properties of the algorithm with numerical experiments in mapped meshes in 1D-3V and 2D-3V.