Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Electron behaviour in three-dimensional collisionless magnetic reconnection A. Perona1, D. Borgogno2 , D. Grasso2,3 1 CFSA, Department of Physics, University of Warwick, UK 2Burning Plasma Research Group, Politecnico di Torino, Italy 3Istituto dei Sistemi Complessi-CNR, Roma, Italy Abstract Aim of this project, still in progress, is to investigate the behaviour of an electron population during the evolution of a spontaneous collisionless magnetic reconnection event reproduced by a fluid formulation in a three-dimensional geometry. In this 3D setting the magnetic field lines become stochastic when islands with different helicities are present1. The reconstruction of the test electron momenta, in particular, can assess the small scale behaviour shown by the fluid vorticity and by the current density during the nonlinear phase of the reconnection process even in the presence of chaoticity. We present here preliminary results of the numerical tool developed on this purpose. 1D. Borgogno et al., “Aspects of three-dimensional magnetic reconnection”, Phys. Plasmas, Vol. 12, 032309, (2005). 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 The fluid model for collisionless magnetic reconnection We consider a fluid model that describes drift-Alfvén perturbations in a plasma immersed in a strong, uniform, externally imposed magnetic field. Effects related to the electron temperature, through the sonic Larmor radius, and to the electron density, through the electron inertia, are retained. Modes with different helicities can be present. They evolve independently during the linear phase, while a strong interaction occurs during the nonlinear stage of the process. The set of equations is closed by assuming constant temperatures for both the ions and the electrons. 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 The fluid model for collisionless magnetic reconnection B = B0 ez+ (x,y,z,t) ez magnetic field with B0=const. and |B0|>>|B| ( = 8p2/B2 <<1), magnetic flux function, stream function. de 2 t 2 t 2 , 2 d e 2 , 2 2 , s 2 2 2 , s z z 2 2 z T. J. Schep, F. Pegoraro, B. N. Kuvshinov, Phys. Plasmas, Vol. 1, 2843, (1994) Two scale lengths: de me / 0 nee2 c / pe electron skin depth Fluid velocity: v=ez+ S Te M i / (eB 2 ) sound Larmor radius 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 The magnetic reconnection model: simulation of the reconnection process nonlinear We adopt a static, linearly unstable magnetic equilibrium ln eq ( x) ln(cosh( x)) linear Single-helicity initial perturbation J x, y, z, t 2 x, y, z, t Jˆ x exp ik y y ik z z k y m / Ly k z n / Lz m, n mode wave numbers along y and z Island width vs time de 0.16 ne 1019 m 3 s 0.32 Te 1keV 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 Single-helicity case (kz0): the ‘fluid’ parallel electric field Linear phase: E|| monopole structure peaked at the X-point of the magnetic island. V /m E|| V /m 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 The electron model In order to verify whether the parallel electric field generated during the reconnection leads to suprathermal energetic generation, a relativistic Hamiltonian formulation* of the electron guiding-center dynamics has been chosen. The unperturbed relativistic guiding-centre phase-space Lagrangian can be written in terms of the guiding-centre coordinates as L p||b - eA R where b B / B m W e p|| m0 v b 1 p||2 p2 2 2 0 mc m m0 * A. J. Brizard, A. A. Chan, Phys. Plasmas 6, 4548 (1999) 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 W e The electron model The relativistic Hamiltonian is H magnetic moment: p||2 c 2 2 Bmc 2 m02 c 4 e p 2 2 m0 B The equations of motion follow from the Hamiltonian in the usual manner H y p y H z pz H py y H pz z where the conjugate momenta to the y and z coordinates are given by py L p||by eAy , y pz L p||bz eAz z 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 The electron equations Az Az bz bz Az bz by p z e y z p|| y z t y z t y z x by g z bz g y Ay by Ay Ay by by bz p y e y z p|| y z t y z t y z by g z bz g y b A B p||2 z p||e z bz m0 me x x x x y m g z by g y bz by Ay B 2 p|| p||e by m0 me x x x x z m g z by g y bz by by by Az Az Az g z p y g y pz g y e y z g z p|| y t y z t y z p|| by g z bz g y Ay by g y p|| e x x Az bz g z p|| e x x Bz Ax 0 Ay Bz x Az eq Ay Ay Ay z e y t y z 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 z The f approach A f approach has been adopted in order to reduce the number of particles required while resolving small fluctuations in the electron distribution function. The distribution function is decomposed in an analytically described background component and the remaining component f f 0 p f p , t analytic markers • In each cell labeled by i in the real space, we calculate • the electron density ne i • the current density ep ji J N j 1 w / (c ) j j (s) i (s) w / m || j j j i c • and the mean longitudinal kinetic energy of electrons Wi p|| p||0 w j (jc ) / 2m i( s ) N j 1 2 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 Kinetic simulations The electron equations have been implemented in a 3-D code, which reads at each time step the fluid fields provided by the reconnection code1. The kinetic code evolves the spatial coordinates, the parallel velocity and the change in the weight of the markers according to the fluid fields. We assume periodic boundary conditions along y and z for the flux and for the stream function. The results presented have been obtained loading 2.5x106 electrons in the 5-D phase space (x,y,z,p||, p). The code has been parallelized by distributing markers among processers using Message Passing Interface (MPI) librairies. 1 A. Perona, L-G Eriksson, D. Grasso, Phys. Plasmas 17, 042104 (2010) 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 Kinetic simulations: single particle trajectory The Poincaré plot of a single test electron crossing the X-point maps the magnetic island as expected. magnetic surface electron (kz=0 case) electron trajectory 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 Numerical simulations: kinetic results (kz0) During the linear phase the electron and fluid current density evolve in good agreement, while the velocity distribution becomes thinner and the amount of electrons in the tails increases. Linear phase: Kinetic current density J Fluid current density J 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 Numerical simulations: kinetic current density Linear growth rate: 2 10 4 A/ m 2 Linear phase 5 106 A/ m 2 Nonlinear phase I t 2 J 0 1 r 2 / a 2 rdr J 0 a 2 / 2 a 0 fluid current * kinetic current J 0 2 I t / a 2 6.4 106 A/m 2 Good agreement with the toroidal current of small Tokamaks such as T-3 (a ≈ 10-1 m, It=100 kA). 14th European Fusion Theory Conference, Frascati, September 26-29, 2011 Further steps The behaviour of the electron population will be followed during the nonlinear phase of the single-helicity case (kz0) in order to confirm the agreement with the fluid results, as already observed in the 2D case. The test electron distribution function will be reconstructed in the presence of multiple-helicity fluid fields. This will allow us to assess the structures of the electron and current density in the stochastic fields that develop during the nonlinear phase of the reconnection process. In order to analyse the influence of the magnetic chaoticity on the electron transport we plan to compare the test electron distribution and the magnetic barriers detected through the analysis of the finite-time Lyapunov exponent ridges1. 1D. Borgogno et al., “Barriers in the transition to global chaos in collisionless magnetic reconnection”, in press Phys. Plasmas (2011). 14th European Fusion Theory Conference, Frascati, September 26-29, 2011