Numerical Simulation Of Neutral Beam Heated Plasma

In order to realize the fusion ignition condition, auxiliary heating must be used toincrease the plasma temperature. It is necessary to simulate the neutral beam heatingwhich is one of the best effective ways to heat plasma.In this paper, the basic theory of heating with neutral beam injection issummarized. The five processes during the neutral beam injection heating: ionization,attenuation, drift motion and fast ion slowing-down. A group of neutral beam injectionheating basic equations has been solved numerically by the Gill method. The timeevolution of Electron and ion temperature in plasmas has been simulated. The resultsare consistent with the experimental results. The results indicate that the particletemperature and temperature climbing speed increase with the neutral beam energyand equivalent flow increasing.The finite volume method is the most successful way in the numericalcalculation of flow and heat transfer problems. Using the finite volume method fornumerical solution of the Fokker-Planck equation can study background particledistribution function evolution in the neutral beam heating conditions. This paperdescribes several methods for solving the Fokker-Planck equation, analysis of thecharacteristics of these methods. This article describes the common four methods forsolving partial differential equations: finite difference method, finite element method,boundary element method and finite volume method, and analysis the advantages anddisadvantages of these methods. A Fokker-Planck equation calculation format hasbeen obtained by the finite volume method, it can guarantee the conservation ofparticle number density. By one dimension numerical simulation of super-Gaussiandistribution of electron self-collision process tends to Maxwell distribution, finitevolume method can efficiently solving the Fokker-Planck equation, it can guaranteeparticle distribution function non-negative, the particle number density and energyconservation. The two-dimension Fokker-Planck equation procedure has beendeveloped, it need further debug.