skip to main content


Title: On coupling fluid plasma and kinetic neutral physics models

The coupled fluid plasma and kinetic neutral physics equations are analyzed through theory and simulation of benchmark cases. It is shown that coupling methods that do not treat the coupling rates implicitly are restricted to short time steps for stability. Fast charge exchange, ionization and recombination coupling rates exist, even after constraining the solution by requiring that the neutrals are at equilibrium. For explicit coupling, the present implementation of Monte Carlo correlated sampling techniques does not allow for complete convergence in slab geometry. For the benchmark case, residuals decay with particle number and increase with grid size, indicating that they scale in a manner that is similar to the theoretical prediction for nonlinear bias error. Progress is reported on implementation of a fully implicit Jacobian-free Newton–Krylov coupling scheme. The present block Jacobi preconditioning method is still sensitive to time step and methods that better precondition the coupled system are under investigation.
; ; ORCiD logo [1] ; ; ; ; ; ;
  1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Publication Date:
Report Number(s):
Journal ID: ISSN 2352-1791; PII: S2352179116302149
Grant/Contract Number:
AC52-07NA27344; 15-ERD-059; AC02-09CH11466
Published Article
Journal Name:
Nuclear Materials and Energy
Additional Journal Information:
Journal Volume: 12; Journal ID: ISSN 2352-1791
Research Org:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org:
Country of Publication:
United States
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Divertor modeling; Charge exchange; Ionization; Recombination; Implicit; Newton–Krylov; UEDGE; DEGAS2; EIRENE; SOLPS; 74 ATOMIC AND MOLECULAR PHYSICS; 70 PLASMA PHYSICS AND FUSION
OSTI Identifier:
Alternate Identifier(s):
OSTI ID: 1349230; OSTI ID: 1366930