Pseudo-transient Continuation Based Variable Relaxation Solve in Nonlinear Magnetohydrodynamic Simulations

Chen, Jin

Princeton University. Plasma Physics Laboratory.

关键词: Efficiency;    Thermal Conductivity;    70 Plasma Physics And Fusion Technology;    Plasma Simulation;    Simulation;   

DOI  :  10.2172/969310
RP-ID  :  PPPL-4476
RP-ID  :  DE-ACO2-09CH11466
RP-ID  :  969310

美国|英语

来源: UNT Digital Library

【 摘 要 】

Efficient and robust Variable Relaxation Solver, based on pseudo-transient continuation, is developed to solve nonlinear anisotropic thermal conduction arising from fusion plasma simulations. By adding first and/or second order artificial time derivatives to the system, this type of method advances the resulting time-dependent nonlinear PDEs to steady state, which is the solution to be sought. In this process, only the stiffness matrix itself is involved so that the numerical complexity and errors can be greatly reduced. In fact, this work is an extension of integrating efficient linear elliptic solvers for fusion simulation on Cray XIE. Two schemes are derived in this work, first and second order Variable Relaxations. Four factors are observed to be critical for efficiency and preservation of solution's symmetric structure arising from periodic boundary condition: refining meshes in different coordinate directions, initializing nonlinear process, varying time steps in both temporal and spatial directions, and accurately generating nonlinear stiffness matrix. First finer mesh scale should be taken in strong transport direction; Next the system is carefully initialized by the solution with linear conductivity; Third, time step and relaxation factor are vertex-based varied and optimized at each time step; Finally, the nonlinear stiffness matrix is updated by just scaling corresponding linear one with the vector generated from nonlinear thermal conductivity.