Ph.D. Proposal Presentation by Dingkang Zhang
Thursday, May 6, 2004

(Dr. Weston M. Stacey, Chair)

"Neutral Particle Transport in Plasma Edge Using Transmission/Escape
Probability (TEP) Method "


Transmission/Escape probability (TEP) method has been demonstrated to be an efficient and accurate calculation for neutral transport in the tokamak edge and divertor configuration. The obvious advantages of the TEP method are ability to simulate neutral transport problems with complicated configurations, widely varying mean free path and strong particle density gradient. However, the comparisons with Monte Carlo for specially designed problems have also shown that a number of approximations to simplify the TEP methodology limit the applications of TEP method, especially for very short or very large mean free path cases.

In this proposed research, the original TEP methodology, which assumed isotropic inward and outward neutral distributions (DP-0) in each computation region for the purpose of calculating transmission probabilities across the region, is extended to calculate linearly (DP-1) and quadratically (DP-2) anisotropic distributions incident at each interface for the purpose of calculating transmission probabilities. To address the non-uniformity of change-exchange neutral sources when the mean free path is small compared to the dimension of the computation region, two different methods will be developed to take into account that the escape of scattered or charge-exchange neutrals is preferentially across the incident surface. One approach is to use more than one spatial expansion functions per region, so that local flux gradients can be taken into account. The alternative approach is to use diffusion theory to calculate the total and directional escape probabilities. The accuracy and effectiveness of these two methods will be investigated by comparison with Monte Carlo method. In addition, A linear approximation will be used to take into account of the gradient of angular flux at interfaces. The research is expected to provide a fast and accurate method for calculation of neutral particle transport in the complex tokamak edge and divertor configuration.