An implicit electrostatic particle-in-cell/Monte Carlo (PIC/MC) algorithm is developed for the magnetized discharging device simulation. The inductive driving force can be considered. The direct implicit PIC algorithm (DIPIC) and energy conservation scheme are applied together and the grid heating can be eliminated in most cases. A tensor-susceptibility Poisson equation is constructed. Its discrete form is made up by a hybrid scheme in one-dimensional (1D) and two- dimensional (2D) cylindrical systems. A semi-coarsening multigrid method is used to solve the discrete system. The algorithm is applied to simulate the cylindrical magnetized target fusion (MTF) pre-ionization process and get qualitatively correct results. The potential application of the algorithm is discussed briefly.
The theory of slow backward-wave amplifications is developed based on electron cyclotron maser (ECM) mechanism employing an initially rectilinear beam, A nonlinear evolution equation is derived to describe the electron energy. Numerical calculations show that the saturated interaction efficiency in this system may exceed 20~, and the saturated interaction length spans 3-6 centimeters. The distinctive interaction mechanism is promising for the design of compact backward microwave amplification devices, Numerical studies are also presented for the slow-wave ECM efficiency with inclusion of Gaussian beam electron velocity spread. It is shown that the velocity spread reduces the interaction e^ciency.
With the saddle point analysis method for the Bessel function structure and property, the convergence problem and the scaling laws of Thomson backscattering spectra are solved and studied in both cases that are for the plane wave laser field without and with applied external constant magnetic field. Some unclear points appeared in previous work are clarified. The extension of the method to a general situation for the laser field with an arbitrary polarization is discussed. We also make a simple analysis and discussion about the optimal spectra dependence of field parameters and its implication to practical applications.
Enhanced electron–positron pair production by frequency chirping in one- and two-color laser pulse fields is investigated by solving the quantum Vlasov equation. A small frequency chirp shifts the momentum spectrum along the momentum axis. The positive and negative frequency chirp parameters play the same role in increasing the pair number density. The sign change of the frequency chirp parameter at the moment t = 0 leads the pulse shape and momentum spectrum to be symmetric, and the number density to be increased. The number density of produced pairs in the two-color pulse field is much higher than that in the one-color pulse field and the larger frequency chirp pulse field dominates more strongly. In the two-color pulse fields, the relation between the frequency ratio of two colors and the number density is not sensitive to the parameters of small frequency chirp added in either a low frequency strong field or a high frequency weak field but sensitive to the parameters of large frequency chirp added in a high frequency weak field.
Based on the quantum Vlasov equation, the effect of frequency chirp on electron-positron pair production is investigated. The cycle parameter, which characterizes the laser field cycle degree within the pulse, is also considered. In both supercycle and subcycle laser pulses the frequency chirp can greatly enhance the momentum distribution function of created pairs and the pair number density. The pair number density created by a supercycle laser pulse is larger than that by a subcycle pulse under the same laser frequency and chirping. There exists an optimal cycle parameter corresponding to the maximum value of the created pair number density for different chirp rates. It is found that the pair number density is sensitive/insensitive to chirping rate when the cycle parameter lies below/above the optimal one.
We benchmark and analyze the error of energy conservation (EC) scheme in particle-in-cell/Monte Carlo (PIC/MC) algorithms by simulating the radio frequency discharge. The plasma heating behaviors and electron distributing functions obtained by one-dimensional (1D) simulation are analyzed. Both explicit and implicit algorithms are checked. The results showed that the EC scheme can eliminated the self-heating with wide grid spacing in both cases with a small reduction of the accuracies. In typical parameters, the EC implicit scheme has higher precision than EC explicit scheme. Some "numerical cooling" behaviors are observed and analyzed. Some other errors are also analyzed. The analysis showed that the EC implicit scheme can be used to qualitative estimation of some discharge problems with much less computational resource cost without much loss of accuracies.