A coarse projective method for simulation of multi-dimensional macroscopic systems,Primal Equation Free Projective Integration(p-EFPI)method,was developed and applied to simulation of plasma ion acoustic wave.The method consists of two parts: Time stepper and time projection.For election velocities, probability distribution function or cumulative distrihution function are recorded in time to project,while,corresponding electron positions are reloaded(restricted)via electron charge density.Ion motion is assumed inherently coarse-grained as compared to electron dynamics,to describe possibly macro scale plasma evolution sufficiehtly well. To preserve full ion kinetic effects,individual ion orbits are tracked in time and simply extrapolate to pmject.This is one attempt to test method in p-EFPI framework on plasma physics phenomenon.Our calculations for ion acoustic wave paradigm indicate that P-EFPI simulation follows PIC results nice before ion trapping in phase space appears.However,as ion trapping takes place,main characteristic still keeps,except for ion trapping characteristic.

Thermonuclear reaction rate describing thermonuclear burning efficiency in thermonuclear fusion is studied.As one species atomic nucleus with Maxwell velocity distribution,another species atomic nucleus velocity is divided into several groups.Velocity group thermonuclear reaction rate and averaged thermonuclear reaction rate are deduced,and simulation codes are obtained.Taking deuterium-tritium thermonuclear fusion reaction as an example,velocity group thermonuclear reaction rates and corresponding averaged thermonuclear reaction rates are calculated and analyzed comparing with exact direct integration results.Finally,relation between averaged thermonuclear reaction rate and velocity group number is studied.

It reports the investigation on calculation methods of ground-state correlation energy of double-electron system atoms (or ions) by using B-splines basic sets. These methods belong to the application category of calculus of variations, and adopt B-splines techniques to account for the radial bases, and at the same time, deal with the angular bases by using Goldman's mixed L method. The calculation results are given for He, H^-,Li⁺,Be⁺⁺,B⁺⁺⁺ atoms (or ions) and compared with those of others' theoratical methods. This comparison demonstrates that our calculation methods have excellent computing accuracy and efficiency.

Based on the multi-configuration Hartree-Fock program code designed by Fischer, a program code is presented for treating the bound states as well as the continuum states, with transforming the radial mesh from ρ=lnZr to ρ=αr+βlnr. As an illustration of the completeness of the basis sets obtained from this program, the results of Thomas-Reiche-Kuhn sums under the Coulomb potential, are given as well.