Application of hybrid particle in cell & lattice

0 downloads 0 Views 56MB Size Report
are effectively accelerated via shock acceleration process trapped in a plasmid ..... Diffusive shock acceleration (DSA, 1st Fermi Acceleration). 3.4 Fermi and ...
MR 2017 The US-Japan Workshop and School on Magnetic Reconnection March 19~23, 2017 Dogo Matsuyama and Miyajima Hiroshima, Japan

Application of hybrid particle in cell & lattice Boltzmann method (HPIC-LBM) in large scale turbulent reconnecting current sheet Bojing Zhu (朱伯靖)1,2 , Jun Lin (林隽)1,2 1 Yunnan 2

Observatories, Chinese Academy of Sciences

Center for Astronomical Mega-Science, Chinese Academy of Sciences

OUT LINE 1. Background and Motivation 2. HPIC-LBM Method 3. Numerical Simulation 1

Background and Motivation

1.1 Large scale turbulent reconnecting current sheet Energy=1032-1033erg=1025-1026 J

Direct observation-LASCO

Time= 1-101 minute

Eparticle=KeV-GeV

HXR ( MR threshold criteria II: CS thinning in nonlinear due to tearing or global resistive instability. III: KineticCS thinning stops —> Sweet-Parker Model. IV: Secondary tearing, and CS becomes fractal.

1.Bubble collapse 2.Decreasing scale 3.Energy density decreasing 4.Vortex shrink

V: Smallest bubble—> Bigger bubble VI: Bigger bubble—> Largest bubble MHDVII: Bubble stop expanding and carrying energy as constant velocity

1.Fragmenting coalescence 2.Increasing scale 3.Energy density increasing 4.Vortex stretching

VIII: Greatest energy release occurs when largest bubble is ejected. IX: HydroGlobal instability appears

1.Energy density constant

Numerical Simulation

3.3 Turbulence evolution

The distance between two adjacent magnetic field line

The bubble’s collapse —> smallest bubble I: Initial CS —> MR threshold criteria

Kinetic-

Ideal MHD left-handed Alfven wave and right-handed Alfven wave on adjacent parallel magnetic field line.------Ion Whistler Beyond ideal MHD II:CS thinning in nonlinear due to tearing or global resistive instability.

Mixed ion-electron inertial length

Numerical Simulation

3.3 Turbulence evolution

The distance between two adjacent magnetic field line

The bubble’s collapse —> smallest bubble III:CS thinning stops —> Sweet-Parker Model.

IV: Secondary tearing, and CS becomes fractal.

Kineticleft-handed Alfven wave and right-handed Alfven wave on adjacent parallel magnetic field line.------Electron electrostatic wave (Low fre)

Note: 1.The CS thickness of LSTMR arrange from (~100m) to (~ 0.2m) 2.If there are several CS, Magnetic energy dissipation efficiency will be greatly increased

Numerical Simulation

3.3 Turbulence evolution The bubble’s collapse —> smallest bubble IV: Secondary tearing, and CS becomes fractal.

KineticFine structure-

Numerical Simulation

3.3 Turbulence evolution

Three bubble coalesce O – point MR happen around adjacent bubble coalescence!

The bubble’s coalescence —> biggest bubble

Kinetic- & MHD-

V: Smallest bubble—> Bigger bubble Evolution time

Domain size

Alfven speed

33

Numerical Simulation

3.3 Turbulence evolution The bubble’s coalescence —> biggest bubble

Kinetic- & MHD-

V: Smallest bubble—> Bigger bubble

The fluctuating component of velocity and magnetic field parts-----pseudo scaler-----ω

The fluctuating component of velocity and magnetic field parts-----pseudo scaler-----J

The bubble size reach to ~500km

34

Numerical Simulation

3.3 Turbulence evolution

The distance between two adjacent magnetic field line

The bubble’s coalescence —> biggest bubble VI: Bigger bubble—> Largest bubble

Kinetic- & MHD-

Numerical Simulation

3.3 Turbulence evolution The bubble’s coalescence —> biggest bubble

Kinetic- & MHD-

VI: Bigger bubble—> Largest bubble

The fluctuating component of velocity and magnetic field parts-----pseudo scaler-----J

The mean component of velocity and magnetic field parts-----pseudo scaler-----B

The bubble size reach to 0.3 wdith of the simulaiton domain

36

Numerical Simulation

3.3 Turbulence evolution VII-IX: biggest bubble -> ejected

MHD-

It is clearly seen that the electrons from the energized group come in to the center of the electron diffusion region and stay longer in the area. In contract, the electrons from non energized group only skim the outer part of the diffusion regions

Numerical Simulation

3.4 Correlation between instabilities and energy release-translation

Numerical Simulation

3.4 Fermi and Betatron acceleration

Numerical Simulation

3.4 Fermi and Betatron acceleration Diffusive shock acceleration (DSA, 1st Fermi Acceleration)

40

Numerical Simulation

3.4 Fermi and Betatron acceleration Diffusive shock acceleration (DSA, 1st Fermi Acceleration)

Numerical Simulation

3.4 Fermi and Betatron acceleration Hybrid Acceleration Diffusive shock acceleration (DSA, 1st Fermi) & Stochastic acceleration( 2nd Fermi)

Numerical Simulation

3.4 Fermi and Betatron acceleration Hybrid Acceleration Diffusive shock acceleration (DSA, 1st Fermi) & Stochastic acceleration( 2nd Fermi)

Simulation Domain

Simulation Domain

Simulation Domain

Simulation Domain

Numerical Simulation

3.4 Fermi and Betatron acceleration Point (14555,16000)

D

E

B C A

Numerical Simulation

3.4 Fermi and Betatron acceleration 8 points View field 800x4200

http://ddl.escience.cn/f/BOwZ

Numerical Simulation

3.4 Fermi and Betatron acceleration

Numerical Simulation

3.4 Acceleration-turbulence evolution

Thank you for your attention!

48

Suggest Documents