An eigenmode equation of ion temperature gradient driven instability in a tokamak plasma with negative magnetic shear is solved numerically by the shooting method.Numerical calculations show that the toroidal rotation sheared flow can modify the sheared slab η_i instability and has a stabilizing η_i role in the plasma core near minimum q magnetic surface.

It proposes three kinds of new difference schemes(exponential type,oscillatory type and exponential-oscillatory type) based on the coefficient structural characteristic and the behaviour of solution in the ODEs.The schemes can describe excellently dependent variables with steep gradient,shear layer and/orhighly oscillations.

The molecular dynamics method is employed to study the sticking and desorption of Xenon atoms from a platinum surface.After constructing correct gas-surface interaction potential,the stochastic classic trajectory method is adopted to solve the motion equation of both gas and solid atoms,then the trajectory and behavior of gas atoms can be obtained.The results show that the sticking probability depends on such factors as wall temperature,gas atom incident enenrgy and interaction potential parameters etc.,the thermal desorption rate of gas atoms from platinum surface satisfies Arrhenius correlation.

The eigenvalue problem of the secular trapped ion temperature gradient mode is solved by employing invariant imbedding method. The numerical methods removing the singularity of the eigen equation at zero point are presented for even mode and odd mode, respectively. The calculation for a concrete physical problem is outlined.

The paper dercribes a numerical method for the resistivity gradient driven turbulence in the Tokamak edge. The mode selection for the mode-mode coupling or nonlinear interaction is discussed. The corresponding computer code for three-dimension magnetic hydrodynamics equations has been constructed, which can save the computation time effectively.

A set of nonlinear second order parabolic partial differential equations with initial-boundary values is approxinated by Crank-Nicholson difference method, and successfully solved by means of Loop-Iteration method. The computation of energy transport in tokamak plasma by ohmic heating provides numerical results in agreement with the theoretial analysis.