Mesh generation for confined fusion plasma simulation

• Published in: Engineering with computers 2016-04, Vol.32 (2), p.285-293
• Main Authors: Zhang, Fan, Hager, Robert, Ku, Seung-Hoe, Chang, Choong-Seock, Jardin, Stephen C., Ferraro, Nathaniel M., Seol, E. Seegyoung, Yoon, Eisung, Shephard, Mark S.
• Format: Article
• Language: English
• Subjects: 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Automatic mesh generation; CAE) and Design; Calculus of Variations and Optimal Control; Optimization; Classical Mechanics; Computational fluid dynamics; Computer Science; Computer simulation; Computer-Aided Engineering (CAD; Control; Finite element method; Fluid flow; Geometric model; Magnetohydrodynamic equations; Magnetohydrodynamic turbulence; Magnetohydrodynamics; Math. Applications in Chemistry; Mathematical analysis; Mathematical and Computational Engineering; MATHEMATICS AND COMPUTING; Mesh generation; Original Article; Particle in cell technique; Plasma; Plasma physics; Plasmas; Software development tools; Systems Theory; Tokamak fusion reactor
• ISSN: 0177-0667; 1435-5663
• DOI: 10.1007/s00366-015-0417-y

XGC1 and M3D- C 1 are two fusion plasma simulation codes being developed at Princeton Plasma Physics Laboratory. XGC1 uses the particle-in-cell method to simulate gyrokinetic neoclassical physics and turbulence (Chang et al. Phys Plasmas 16(5):056108, 2009 ; Ku et al. Nucl Fusion 49:115021, 2009 ; Admas et al. J Phys 180(1):012036, 2009 ). M3D- C 1 solves the two-fluid resistive magnetohydrodynamic equations with the C 1 finite elements (Jardin J comput phys 200(1):133–152, 2004 ; Jardin et al. J comput Phys 226(2):2146–2174, 2007 ; Ferraro and Jardin J comput Phys 228(20):7742–7770, 2009 ; Jardin J comput Phys 231(3):832–838, 2012 ; Jardin et al. Comput Sci Discov 5(1):014002, 2012 ; Ferraro et al. Sci Discov Adv Comput, 2012 ; Ferraro et al. International sherwood fusion theory conference, 2014 ). This paper presents the software tools and libraries that were combined to form the geometry and automatic meshing procedures for these codes. Specific consideration has been given to satisfy the mesh configuration and element shape quality constraints of XGC1 and M3D- C 1 .