二维电磁模型 基本方程与无量纲化 基本方程
无量纲化 方程化为
二维时的方程
网格划分 时间上利用蛙跳格式
计算步骤
稳定性条件
例子：磁场重联的二维粒子模拟 Magnetic reconnection rapidly converts magnetic energy into plasma energy, which leads to heating and acceleration of ions and electrons.
Topology changes of magnetic field lines during magnetic reconnection
Particle simulations[Fu and Lu, 2006] With 2D particle-in-cell simulations we investigate the influence of the guide field on the electron acceleration near X-point and O-point.
Initial conditions and boundary conditions: Initial conditions 1D Harris current sheet in the (x,y) plane Initial flux perturbation is introduced Boundary conditions X direction: periodic Y direction: ideal conducting boundary condition for EM fields, Reflection condition for particles
Parameters
Time evolution of the reconnection flux for
Contours of (a) ,(b) , (c) ,(d) at for
Contours of (a) ,(b) , (c) ,(d) at for
Contours of (a) ,(b) , (c) ,(d) at for
Typical trajectories in (x,y) plane, one passes through X-point (from to , the other is trapped near O-point (from to for
The time evolution of (a) the kinetic energy, (b) , (c) ,(d) for The electron passes through X-point .
Comparison with others[Hoshino,2005; Pritchett, 2006]
The time evolution of (a) the kinetic energy, (b) , (c) ,(d) for The electron is trapped near O-point .
Typical trajectories in (x,y) plane, one passes through X-point (from to , the other is trapped near O-point (from to for
The time evolution of (a) the kinetic energy, (b) , (c) ,(d) for The electron passes through X-point .
The time evolution of (a) the kinetic energy, (b) , (c) ,(d) for The electron is trapped near O-point .
Typical trajectories in (x,y) plane, one passes through X-point (from to , the other is trapped near O-point (from to for
The time evolution of (a) the kinetic energy, (b) , (c) ,(d) for The electron passes through X-point .
The time evolution of (a) the kinetic energy, (b) , (c) ,(d) for The electron is trapped near O-point .
The positions of the energetic electrons
Comparison with others [Pritchett, 2006] no guide field
With guide field 1.0B0
Discussion 1.Observations of energetic electron in ion diffusion region in magnetotail [Oieroset, 2002]
Figure 5 shows the plasma temperature, magnetic field vectors, high-speed flows and energetic electron differential fluxes. The bottom four panels denote electron differential fluxes obtained from the RAPID on the four satellites from 35.1 to 244.1keV. A depletion in the energetic electron fluxes in the diffusion region was detected by all the four satellites. The duration is about 162s. Please note that the first two low energy channels of c3 is not well calibrated (Private communication from Q. G. Zong), while the other energy channels have the same depletion as other satellites. Similarly, a local minimum of the plasma temperature can also be found near the center of diffusion region
From Øieroset et al NATURE 2001 Fig 2. By WIND satellite, the fluxes of energetic electrons up to ~300keV peak near the center of the diffusion region and decrease monotonically away from this region. No secondary acceleration was found in the reconnection. Note: the initial guild field is about 50% of the total magnetic field magnitude during the magnetic reconnection. From Øieroset et al PRL 2002 Fig 1.
Potential applications in Solar atmosphere [Drake et al., 2006]?
