An algorithm allowing a fast solution of stationary and time dependent 1-d diffusion problems for the case of finite ionization and recombination rates is presented. Results of numerical computations and an analytic stability analysis show its unconditional stability. The computational effort involved rises only linearly with the number of ionization stages included. An application to the diffusion of oxygen and iron impurities in JET is described
No takes yet. Share an insight, caveat, or question.
Lackner et al. (1982) studied this question.