A 3D immersed finite element method with non-homogeneous interface flux jump for applications in particle-in-cell simulations of plasma-lunar surface interactions

详细信息    查看全文

• 作者：Daoru Han^a ; ^{daoruhan@usc.edu" class="auth_mail" title="E-mail the corresponding author} ; Pu Wang^b ; ^{wangpu@vt.edu" class="auth_mail" title="E-mail the corresponding author} ; Xiaoming He^c ; ^{hex@mst.edu" class="auth_mail" title="E-mail the corresponding author} ; Tao Lin^d ; ^{tlin@math.vt.edu" class="auth_mail" title="E-mail the corresponding author} ; Joseph Wang^a ; ^{josephjw@usc.edu" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Interface problems ; Non-homogeneous flux jump conditions ; Immersed finite elements ; Particle-in-Cell ; Plasma environments ; Lunar surface
• 刊名：Journal of Computational Physics
• 出版年：2016
• 出版时间：15 September 2016
• 年：2016
• 卷：321
• 期：Complete
• 页码：965-980
• 全文大小：2062 K

文摘

Motivated by the need to handle complex boundary conditions efficiently and accurately in particle-in-cell (PIC) simulations, this paper presents a three-dimensional (3D) linear immersed finite element (IFE) method with non-homogeneous flux jump conditions for solving electrostatic field involving complex boundary conditions using structured meshes independent of the interface. This method treats an object boundary as part of the simulation domain and solves the electric field at the boundary as an interface problem. In order to resolve charging on a dielectric surface, a new 3D linear IFE basis function is designed for each interface element to capture the electric field jump on the interface. Numerical experiments are provided to demonstrate the optimal convergence rates in L²$L^2$ and H¹$H^1$ norms of the IFE solution. This new IFE method is integrated into a PIC method for simulations involving charging of a complex dielectric surface in a plasma. A numerical study of plasma–surface interactions at the lunar terminator is presented to demonstrate the applicability of the new method.