Reconnection of orthogonal (NN6) Baxial=1+cos(r*pi/R) flux tubes
Lundquist number S=R*vA/(pi*eta)=3750
Initial peak velocity magnitude: 4vA/30, as opposed to
vA/30
Field in vertical tube runs top to bottom, in horizontal tube
runs left to right (so mutual helicity is negative, reconnection
should generate left handed twist).
Isosurfaces of B=Bmax/3 (not normalized)
tvA*pi/R=[0,120,176],
tvA*pi/R=[200,251,415]
Click on image for movie
Below: fieldlines running through lagrangian trace
elements placed in the plasma.
Only fieldlines running through
elements initially in the vertical flux tube are shown.
The horizontal flux tube is in front of the vertical flux tube
from this viewpoint.
Fieldline area at any point along the fieldline
is proportional to |B| there. Any part of a fieldline at |B| less
than 1/10th of the global maximum is not plotted.
tvA*pi/R=[0,120,176],
tvA*pi/R=[200,251,415]
Click on image for movie
Energy vs time for the reconnection simulation and the reference
simulations.
The reference simulation is a single tube and the same S with no
initial velocity.
For this figure/movie, the fieldline radius at any point along the fieldline
is proportional to Eparallel=eta*J.B/|B| there, normalized by the peak Eparallel
in the low velocity (nn6f2) simulation.
As Eparallel is a measure of reconnection, this emphasizes those fieldlines
which are reconnecting at any point in time.
Diagonal slices of B, J, and v : the diagional direction
is along y=z (bottom left to top right in fieldline figure above)
so that the component of magnetic field
common to both flux tubes initially points straight into
the page.
The times for the panels in these three figures are:
tvA/R=[0.0,670.,870.,980.,1300.,2000.]
Current:
Velocity:
The peak velocities at each timestep, normalized by the peak Alfven speed at each
step, are: [0.13, 0.13, 0.17, 0.25, 0.35, 0.37]. The peak velocity (at crude
time resolution) for the whole simulation was 0.375, at a time halfway between
the last two frames shown.
