A three dimensional object oriented parallel particle simulation code PLASIM3D (using the first three letters of the words "plasma" and "simulator") is developed with the purpose of studying the interaction of laser and plasma. The arithmetic of the code is introduced and its performance is tested and analyzed. The linear speedup is acquired in the present test.

The Lindman's absorption boundary condition, which has been widely used in the FDTD method, is generalized to a special form suitable to the electromagnetic field calculation from the starting point of the absorption principle. Our 2D numerical simulations indicate that the P-polarized and S-polarized electromagnetic waves with a wide varity of incident angles can be well absorbed by left and right boundaries. The reflectivity of the electromagnetic energy at either boundary is less than 0.5%.

A 2(1/2) dimensional particle simulation code PLASIM (Plasma and Simulator) is developed with the purpose of studying the interaction of laser and plasma.The algorithm of the code is introduced.A main characteristics of the code is that it can simulate collision effects by adopting β-particle cloud.Another characteristics is that electrons and ions are both relativistic. So the computation results are more believable when incident laser is ultra-intense.

Two kinds of parallel 2(1/2)-dimensional particle simulation algorithm is designed.One of the algorithm codes is constructed under the distributed-memory process system.One 3-Dimensional Maxwellian velocity algorithm is presented and the Lindman absorbent boundary condition is generalized.The parallel efficiency is tested and analyzed. As an example,the result of the channel effect in plasma by ultrahigh laser is given.

The pseudo-current method is outlined proposed by Marder B for eliminating charge conservation errors in electromagnetic particle-in-cell codes. The new method is shown to be effective in removing charge errors by the choice of current deposition algorithms. In addition, a new model for field emission from conducting boundary surfaces is described. This model is proved to be practicable by the results of the particle simulation of plasma erosion opening switch involving complex field emission boundaries.

The 3-D Fast Multipole Method, a kind of gridless particle simulation method, is described and improved to obtain higher accuracy. A spherical particle with reduced dimensional motion model is presented. As examples, the spatial distribution of the potential in dust plasma and the instability produced by electrons streaming through ions are simulated. The results of simulations are in agreement with the theory.

A stretched one and two halves dimensional relativistic electromagnetic particle code is used to simulate the Cerenkov emission process from a relativistic electron beam when it is passing through the hole of a dielectric-lined waveguide at a speed more than the speed of light in the dielectric. Simulation results show that Cerenkov process is a potential high-power source of short-wavelength radiation.