Viscous effects in timedependent planar reconnection
Department of Mathematics, University of Waikato, P. B. 3105, Hamilton, New Zealand
email: cka3@waikato.ac.nz; i.craig@waikato.ac.nz; yuril@waikato.ac.nz
Received: 8 July 2011
Accepted: 23 August 2011
Context. Viscous dissipation is expected to play a significant part in energy release in solar flares, yet the role of viscosity in a weakly resistive plasma of the solar corona remains unclear.
Aims. We attempt to clarify the role of viscous effects in magnetic reconnection by performing simulations of reconnection in planar periodic geometry in an incompressible viscous resistive plasma.
Methods. We consider magnetic reconnection, driven by largescale vortical flows. We use both the classical shear viscosity and the Braginskii form for the ion parallel viscosity in a magnetised plasma. We determine the scalings of the current sheet parameters and the global rates of resistive and viscous dissipation. We use steadystate exact solutions and scaling arguments to interpret the numerical results.
Results. We show that, regardless of the form of viscosity, the resistive nonviscous analytical solutions for flux pileup merging provide a very good approximation of the numerical results in the reconnecting current sheet. We find no evidence for a viscoresistive scale. Numerical results for a highly sheared magnetic field, however, appear to deviate from the analytical predictions in the case of the Braginskii viscosity.
Key words: magnetic reconnection / magnetohydrodynamics (MHD) / Sun: flares
© ESO, 2011
1. Introduction
Viscous effects are expected to play a significant part in a wide variety of dynamic phenomena in the solar convection zone (e.g, Rüdiger & Kitchatinov 1997), in the solar atmosphere (e.g., Hollweg 1985), and in the solar wind (Montgomery 1992). Of particular interest is the role of viscosity in the process of magnetic reconnection in the solar corona (Hollweg 1986). The classical electric resistivity η of the hightemperature coronal plasma is so low that additional nonideal processes are needed, for instance, to explain the observed rapid energy release in solar flares (for a review, see Priest & Forbes 2000).
Despite a number of studies, the role of viscosity in magnetic reconnection remains unclear. On the one hand, SweetParker (Parker 1957) style arguments suggest that a sufficiently large classical shear viscosity ν leads to thicker reconnecting current sheets, slower reconnection inflows, and slower reconnection outflows (Park et al. 1983). Notably, the scaling ~ (ην)^{1/4} is predicted for the thickness of the current sheet, which should be contrasted with the usual resistive SweetParker scaling ~η^{1/2}. The viscoresistive scaling also emerges in an analytical model for slow magnetic reconnection at a twodimensional neutral point (Titov & Priest 1997). On the other hand, a number of solutions for steady magnetic merging in viscous resistive plasmas neither contain a viscoresistive scale nor predict the slowdown of the reconnection flow due to viscosity (e.g., Sonnerup & Priest 1975; Besser et al. 1990; Jardine et al. 1992).
The picture is further complicated by the fact that the classical isotropic viscosity, while adequately representing viscous effects in a weakly magnetised plasma (e.g., Rüdiger & Kitchatinov 1997), is inaccurate in a strongly magnetised plasma such as the solar corona. For a coronal plasma in which the mean time between momentumchanging collisions exceeds the inverse proton gyrofrequency, the Braginskii (1965) viscosity that accounts for the anisotropies introduced by the magnetic field is more accurate (Hollweg 1986). Exact analytical solutions and scaling arguments, based on the Braginskii (1965) viscosity, predict that the current sheet parameters are independent of the viscosity coefficient (Litvinenko 2005; Craig & Litvinenko 2009, 2010).
The contradictory results of the steady analytical models motivate us to investigate the structure of the reconnection region in an incompressible viscoresistive plasma, using numerical simulations of timedependent magnetic reconnection. Numerical results and scaling arguments suggest the effect of compressibility is relatively minor for the flow driven merging solutions discussed here (Litvinenko & Craig 2003). Accordingly, we investigate reconnection in a “closed” twodimensional doubly periodic reconnection region, thus simulating advectiondriven reconnection. We also use scaling arguments to interpret our numerical results.
The paper is organised as follows. We introduce the incompressible magnetohydrodynamic (MHD) equations and discuss resistive and viscous dissipation processes in Sect. 2. We summarise the available analytical solutions and scaling law arguments for steady magnetic merging in Sect. 3. In Sect. 4, we describe our simulations of timedependent magnetic reconnection. We use the numerical results to obtain scalings for both the local current sheet parameters and the global rates of resistive and viscous dissipation. We specifically compare the effects of the classical viscosity and the Braginskii viscosity. We also present an exploratory analysis of the case of strongly sheared reconnecting fields. In Sect. 5 we present our conclusions.
2. Formulation of the problem
2.1. MHD equations
The MHD system for an incompressible resistive viscous plasma consists of the momentum and induction equations, supplemented by the divergencefree constraints for the velocity v and magnetic field B. In dimensionless form we have Here p is the plasma pressure and η the normalised resistivity. The viscous force is defined by the stress tensor whose form is discussed below. The MHD equations have been scaled using the typical coronal parameters: magnetic field strength B_{c} = 10^{2} G, length scale l_{c} = 10^{9.5} cm, and number density n_{c} = 10^{9} cm^{3}. Velocities are expressed in units of the Alfvén speed v_{A} = B_{c}/(4πρ_{c})^{1/2} ≃ 10^{9} cm s^{1}, where ρ_{c} = m_{p}n_{c} is the mass density and m_{p} is the proton mass.
Below we consider two alternative forms for the viscosity tensor: the classical isotropic viscosity (4)and the Braginskii (1965) form for the ion parallel viscosity (5)As usual, summation over repeated indices is assumed. Although the form (4) is commonly employed, it becomes inaccurate in a magnetised plasma such as the solar corona. In such cases the Braginskii form that accounts for the anisotropies introduced by the magnetic field is more accurate (Hollweg 1986).
Two dimensionless parameters control the rate of energy dissipation, the resistivity η and the viscosity ν. Notably, although both are ≪ 1, the coronal temperature interval (2 − 10) × 10^{6} K implies that ν ≫ η, and so viscous dissipation is likely to dominate resistive dissipation in solar flares (e.g., Craig & Litvinenko 2009, and references therein). We use to scale the global energy dissipation rate. The dimensionless rate of resistive (Ohmic) dissipation is given by (6)where J = ∇ × B is the current density. The viscous dissipation rate is given by (7)
2.2. Planar momentum and induction equations
Planar solutions for the magnetic and velocity fields can be expressed using flux and stream functions: (8)Since these forms satisfy Eq. (3) identically, the MHD equations reduce to where we use the Poisson bracket notation typified by
Viscous effects are represented by (11)
3. Analytical models of magnetic merging
3.1. An exact viscoresistive model: classical viscosity
The simplest exact solution for magnetic merging is that of Sonnerup & Priest (1975). A stagnation point flow supports the steady annihilation of straight field lines advected towards a strong current layer centred on x = 0. Specifically the MHD equations are satisfied by where the Dawson function is defined by
and μ^{2} = α/2η. Here α determines the amplitude of the velocity field and E is a uniform electric field in the zdirection. The Dawson function identifies (14)as the thickness of the current layer. This result implies the inflow and outflow speeds associated with the current sheet and v_{out} ≃ α. This solution can be generalised to include vortical and shear flows (Besser et al. 1990; Jardine et al. 1992). Another generalisation, valid only if ν = 0, is to allow the merging of curved field lines, leading to Xpoint reconnection (Craig & Henton 1995). None of these refinements significantly changes the scaling laws derived below.
Using Eqs. (4) and (7), the classical viscous dissipation rate in volume V is given by (15)(In our numerical simulations, the volume is that of the cube 0 ≤ x, y, z ≤ 1.) The resistive dissipation is easy to estimate if the thinness of the sheet is exploited. Thus if B_{s} is the peak field strength in the current layer, then J ≃ B_{s}/x_{s} is the current density, and we have from Eq. (6) that (16)Similarly, the rate of flux annihilation is (17)The key aspect of these results is the presence of a single small scale x_{s} determined by the resistivity. Although viscous losses do occur, these arise from maintaining the global velocity field that drives the merging. By comparing Eqs. (15) and (16) and treating α and B_{s} as parameters of order unity, we see that viscous dissipation should dominate when ν > η^{1/2}. More physically, if we assume (as in SweetParker merging) that exhaust speeds are related to the strength of the current layer, then it follows that α ≃ B_{s} (see Litvinenko & Craig 2000, for details), and so (18)A similar result is obtained when the Braginskii viscosity is considered, as we now demonstrate.
3.2. An exact viscoresistive model: Braginskii viscosity
Exact solutions for magnetic merging, supported by vortical flows, have been obtained for the case of Braginskii viscosity. Both planar merging (Litvinenko 2005) and merging driven by threedimensional flows (Craig & Litvinenko 2009) have been analysed.
We consider only the simplest planar solution, namely Here the parameter 0 ≤ ω < π determines the velocity gradient on the inflow boundary. Note that, as ω → 0, the stagnationpoint flow of the previous section is recovered. The parameters of the current layer (namely its thickness x_{s} and the inflow v_{in} and outflow v_{out} speeds) are all found to scale as those in the merging solution based on the classical viscosity.
The dimensionless dissipation rate, associated with the Braginskii viscosity, follows from Eqs. (5) and (7): (21)which reduces to 3να^{2} V as ω → 0. Since the reconnection rate and Ohmic dissipation rate are essentially unchanged from the previous model (at least provided ω does not approach π), we again recover Eq. (18) within factors of order unity.
3.3. The viscoresistive current layer
The previous two cases are based on exact solutions for largescale velocity fields that drive the merging. By way of contrast, the following argument depends mainly on the local properties of the current layer. In particular, as in Park et al. (1983) and Craig & Litvinenko (2010), SweetParker style arguments can be used to infer the emergence of a viscoresistive scale in the case ν ≫ η.
The argument is based on mass continuity, v_{in} = x_{s}v_{out}, steady flux transfer, v_{in}B_{s} = ηB_{s}/x_{s}, and pressure balance along and across the sheet. For a sufficiently large ν, the classical shear viscosity slows down the outflow and thickens up the current layer with the result that (22)The scaling for viscous dissipation is now based on the assumption that viscous damping is effective only within the narrow current layer–what happens outside the sheet is effectively neglected. The resulting dissipation scalings (23)are quite different from the previous exact analytical models. These scalings (which correct a misprint in Eq. (20) of Craig & Litvinenko 2010) suggest that a large classical viscosity ν ≫ η should slow down the rate of magnetic merging. Note that when ν = η we obtain the traditional resistive SweetParker scalings.
It is interesting that a similar analysis based on the Braginskii viscosity does not lead to the emergence of a viscoresistive scale, as pointed out by Craig & Litvinenko (2010). However, the scale x_{s} ~ (ην)^{1/4} has also been obtained in an analytical model of slow magnetic reconnection at a twodimensional Xpoint (Titov & Priest 1997), as well as in numerical studies of Xpoint collapse for both classical and Braginskii viscosities (Craig et al. 2005; Craig 2010).
To summarise, exact reconnection solutions exist for both classical and Braginskii viscosities, which do not contain a viscoresistive scale. By contrast, some analytical models and Xpoint simulations predict the emergence of the viscoresistive scale x_{s} ~ (ην)^{1/4}. Moreover, SweetParker style scaling arguments imply that this scale should almost universally emerge in a viscous plasma with ν > η.
4. Viscoresistive simulations of magnetic reconnection
4.1. Numerical simulation
In this section we describe a series of sidebyside, flowdriven reconnection experiments using the Braginskii and classical viscosities. Typically the viscous coefficient ν is fixed at some representative value and scalings are derived by systematically reducing the resistivity η ≪ ν. To provide a further check on the results we also perform a “control” resistive simulation, based on the common numerical expedient of setting ν = η. That is, we would like to compare the viscoresistive reconnection scalings with those in the case ν = η, but, as is well known, setting ν = 0 is susceptible to numerical difficulties. Accordingly we set ν = η to approximate the purely resistive case (Biskamp 1994; Craig & Litvinenko 2010). Our aim is twofold: first, to explore the emergence, or otherwise, of a viscoresistive scale; second, to quantify the differences that develop when classical and Braginskii viscosities are employed.
The results are obtained by numerically solving Eqs. (9) and (10), using a doubly periodic code over the region 0 ≤ x, y ≤ 1. Initial conditions are derived by assuming a vortical flow velocity field (clockwise in the first quadrant) given by (24)where α_{0} > 0 sets the amplitude of the initial flow. We assume an initial magnetic disturbance of the form (25)where g_{0} > 0 sets the amplitude of the magnetic disturbance and β_{0} ≥ 0 controls the level of magnetic shear in the merging. In the simplest case of “head on” reconnection (β_{0} = 0), initially straight field lines, washed together by the inflow, rapidly evolve, forming a well defined current layer centered on the origin.
To obtain resistive scalings, the parameters α_{0}, β_{0} and ν are held fixed while η and g_{0} are systematically reduced. Note that simply reducing η while keeping g_{0} fixed would lead to high pressure, flux pileup current sheets that would feed back on the flow and stall it. Accordingly, g_{0} is adjusted to ensure that peak fields B_{s} in the reconnecting current layers have magnitudes of order unity. This allows us to obtain scalings at resistivities that are limited only by numerical resolution in the computation (the lowest values of resistivity considered here require mesh sizes ≲ 10^{3}). In practice we take α_{0} = 1 and ν = 0.004 and allow resistivities in the range 10^{4.5} ≤ η ≤ 10^{2}. Most of our simulations apply to headon reconnection, β_{0} = 0, but some preliminary results for sheared reconnection are given in Sect. 4.4.
As a final point we note that the Braginskii viscosity tensor (5) cannot be applied to fields that are very weak. Since our computation involves null points, it is convenient to adopt a form that effectively interpolates between the classical and Braginskii viscosities. In practice we adopt the “Liley” form of Hosking & Marinoff (1973).
Fig. 1 Contour plot of J taken at the time of maximum current density at the origin. The parameters are α_{0} = 1, β_{0} = 0 and η = ν = 0.004. 

Open with DEXTER 
Fig. 2 Arrow plot of v taken at the time of maximum current density at the origin. The parameters are the same as in Fig. 1, namely α_{0} = 1, β_{0} = 0 and η = ν = 0.004. For visual clarity only selected velocity vectors are displayed. 

Open with DEXTER 
4.2. The reconnecting current sheet
Figures 1 and 2 show an example of a headon Braginskii viscosity simulation with parameters α_{0} = 1, β_{0} = 0 and η = ν = 0.004. The dominant feature, which reaches maximum development after approximately one Alfvén time in the present simulation, is a strong current sheet centred on the origin. As discussed by Heerikhuisen et al. (2000), there is a good deal of numerical evidence to suggest that the properties of the current layer, when taken at the time of maximum current, can be accurately described by the steady merging models such as those discussed in Sect. 3.
Fig. 3 The ycomponent of the magnetic field in the current layer in the case of headon merging. The parameters are the same as in Figs. 1 and 2. The solid line is the numerical magnetic field and the dashed line is the analytical solution from (13). 

Open with DEXTER 
More specifically, suppose we use the fields displayed in Figs. 1 and 2 to evaluate the rate of flux transfer E and the inflow velocity amplitude α_{0}. The analytical solution (13) based on these parameters can then be compared with the simulated current layer. Figure 3 displays a slice of the simulated current sheet alongside the analytical model. Despite the fact that the analytic model is purely resistive and the vorticity parameter ω = 0, we see that the model provides a good representation of the dynamically evolved field in the current layer. The most obvious variations occur in the outer regions dominated by the largescale vortical flowsthe regions not accurately represented by the analytic, linear inflow model of (13). The implication is that Braginskii viscosity simulations are unlikely to undermine the form of the purely resistive current layer.
4.3. Headon magnetic merging simulations
We now turn to quantifying the key parameters v_{out}, v_{in}, and x_{s} of the current layer and relating these to the global Ohmic and viscous energy release rates. Current sheet thickness is measured along the x axis at the level of half current maximum. Inflow and outflow velocities are measured on the x and y axes respectively at the current half maximum level. Figure 4 shows a comparison of the classical, Braginskii and ν = η control current sheet outflow speeds, calculated at the time of maximum current. Minor variations are apparent, and the outflow is generally slower in the case of the classical viscosity, yet the outflow speeds are essentially invariant with η. In particular, there is no evidence for a slowdown, v_{out} ~ (η/ν)^{1/2}, which would be caused by a viscoresistive scale in the case of classical viscosity (Sect. 3.3). The analytical scalings of Sect. 3.1 are clearly more accurate
Fig. 4 Current sheet outflow speed comparisons. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER 
Fig. 5 Current sheet inflow speed comparisons. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER 
This interpretation is borne out by an analysis of the inflow speeds v_{in} and the current sheet thickness x_{s}, which both scale as η^{1/2}, as confirmed by the plots in Figs. 5 and 6. In both cases it is very difficult to distinguish the Braginskii, classical and ν = η control plots. The classical shear viscosity does provide the thickest current layer but the effect is marginal. This should be contrasted with the prediction x_{s} ~ (ην)^{1/4} that follows from SweetParker style arguments (Park et al. 1983).
The relation x_{s} ~ η^{1/2} for the sheet thickness, coupled to the constraint B_{s} ≃ 1, leads to flux transfer and Ohmic dissipation rates that also follow the η^{1/2} trend. These scalings are confirmed numerically for all three regimes (the ν = η control, the classical and Braginskii viscosities).
More interesting is the behaviour of the global viscous losses W_{ν}. As indicated in Fig. 7, the Braginskii and classical viscous losses are effectively invariant, but the ν = η control strongly decreases with η. The decline of the ν = η control, however, is slower than the expected W_{ν} ~ η scaling. We have confirmed that this behaviour is due to strong vortical flows associated with the very weak levels of viscous damping. When we weaken the magnetic field, washed in by the flow, the distortion diminishes and the expected scaling is regained. Nevertheless, this behaviour does indicate the enhanced sensitivity of the computation in the case of weak damping of the global velocity field.
Figure 8 plots the ratio of the viscous and resistive energy dissipation rates W_{ν}/W_{η}. As anticipated by Eq. (18) in the case of fixed ν, W_{ν}/W_{η} scales as η^{1/2} as η is reduced. We see that, independent of the form of the viscosity, W_{ν} exceeds W_{η} by almost a factor of two at the lowest resistivity levels η ≤ 10^{4}.
In summary, all the resistive scalings presented in this section appear consistent with the simple analytical models of Sects. 3.1 and 3.2. That is, for the flowdriven reconnection simulations considered here, the resistive scaling laws are effectively identical for both the Braginskii and classical viscosities. The limited role of the traditional scaling arguments (Park et al. 1983) should be emphasised. Although they provide a reliable description of the reconnecting current sheet in a nonviscous plasma, they do not give a valid description of magnetic reconnection in a viscous plasma, at least in the present, flowdriven simulations.
Fig. 6 Current sheet thickness comparisons. The dotted line shows η^{1/2} scaling. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER 
Fig. 7 Viscous dissipation rate comparisons. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER 
Fig. 8 The ratio W_{ν}/W_{η} for classical (diamonds) and Braginskii (crosses) viscosity. The dotted line shows η^{ − 1/2} scaling. The parameters are α_{0} = 1, β_{0} = 0 and ν = 0.004. 

Open with DEXTER 
Fig. 9 Current contour plot for sheared reconnection in the Braginskii case at time of maximum current density at the origin. The parameters are α_{0} = 1, β_{0} = 0.7, η = 0.004, and ν = 0.004. 

Open with DEXTER 
4.4. Sheared magnetic reconnection simulations
As pointed out previously, it is possible to construct an analytical model for sheared magnetic reconnection only in the nonviscous case (Craig & Henton 1995). Guided by those exact resistive solutions, we expect that, compared with the headon case, viscoresistive reconnection of sheared magnetic field lines should occur in thicker current sheets, which should lead to reduced energy dissipation and flux transfer rates. Figure 9 shows the current contour plot for a typical Braginskii run with a high value of shear (β_{0} = 0.7). As in the headon case, there is a well defined current sheet in the vicinity of the magnetic null, but there is now significant warping of the ends of the sheet.
Our numerical results show that, even with relatively high levels of shearing, the classical viscosity results essentially follow the same scalings as in the headon case (x_{s} ~ η^{1/2}, W_{η} ~ η^{1/2}, W_{ν} ~ η^{0}). Systematic deviations, however, appear in the case of the Braginskii viscosity for high levels of shearing. These deviations from the analytical predictions imply that the structure of the current sheet is significantly modified by the anisotropic Braginskii viscous forces when the reconnecting magnetic field lines are strongly sheared. Unfortunately, as yet, viscoresistive analytical solutions are only available for headon reconnection.
Figure 10 illustrates the main new effect of strong magnetic shear (cf. Fig. 8). It shows W_{ν}/W_{η} for the classical and Braginskii viscosities in the case β_{0} = 0.7. For the classical viscosity, both the resistive dissipation rate W_{η} and the viscous dissipation rate W_{ν} agree with the analytical arguments, leading to Eq. (18). The ratio W_{ν}/W_{η} for the Braginskii viscosity, however, appears to follow a different scaling (closer to η^{1}). We confirmed that this deviation is due to a slow down in the resistive dissipation rate (this scales closer to η rather than η^{1/2}) in the strongly sheared Braginskii simulation.
Fig. 10 The ratio W_{ν}/W_{η} for classical (diamonds) and Braginskii (crosses) viscosity. The dotted line shows η^{ − 1/2} scaling. The parameters are α_{0} = 1, β_{0} = 0.7 and ν = 0.004. 

Open with DEXTER 
5. Conclusions
Given the smallness of the electric resistivity in most astrophysical application, the search for a robust fast mechanism of magnetic reconnection has provided the motivation for a great deal of work on reconnection theory. Both analytical and numerical studies suggest that viscous dissipation may be capable of dominating resistive dissipation under a wide variety of coronal conditions (Hollweg 1986; Litvinenko 2005; Craig & Litvinenko 2009). However, as we discussed in Sect. 3, analytical studies of magnetic reconnection in a resistive viscous plasma lead to contradictory results. Some analytical solutions predict the formation of a current layer on the viscoresistive scale ~(ην)^{1/4}, while other solutions do not contain this scale. This key difference leads to different predictions for the parameters of the reconnecting current sheet and the dependence of the energy release rate on the viscous and resistive coefficients.
In this paper, we attempted to clarify the issue by performing numerical simulations of planar magnetic reconnection in periodic geometry. We considered headon reconnection, driven by largescale vortical flows in an incompressible plasma. We used both the classical shear viscosity and the Braginskii form for the ion parallel viscosity in a magnetised plasma. Somewhat surprisingly, our numerical results show that the parameters of the reconnecting current sheet (its thickness, the inflow and outflow speeds) are accurately described by purely resistive analytical scalings of flux pileup magnetic merging, regardless of the form of the viscous stress tensor. The computed global resistive and viscous energy dissipation rates also closely follow Eq. (18), confirming that the current sheet thickness is effectively independent of viscosity for both the classical and Braginskii forms.
We have also made a preliminary study of the properties of the reconnection region in the case of a strong magnetic shear. It appears that in the Braginskii case magnetic shearing can modify the dissipation scalings, implying that the properties of the reconnecting current sheet depend on the form of the viscosity. This potentially important result requires further investigation. Interpreting numerical results in strongly sheared reconnection is not straightforward, since the available viscoresistive analytical solutions are only able to describe headon merging.
To summarise, we find no evidence of the viscoresistive scale ~(ην)^{1/4} and the corresponding energy release rate scalings in our simulations of planar magnetic reconnection in periodic geometry. This finding contrasts sharply with predictions of some steady analytical models, SweetParker style scaling arguments, and timedependent simulations in a closed Xpoint geometry. There appears to be no simple criterion for the emergence of the viscoresistive scale.
Our results may be partly due to there not being enough time for the viscoresistive scale to develop: we took the current sheet parameters at the time of the first current maximum (first implosion), and it may be that the scale would develop after several implosions. Either way, the precise conditions required to ensure the emergence of a viscoresistive scale remain unclear. We believe that further analytic and numerical work is warranted, given the potential importance of viscous dissipation
arising from either classical or Braginskii viscosity in solar flares and other astrophysical applications.
Acknowledgments
This work was supported by NSF grant ATM0734032 and NASA grant NNX08AG44G.
References
 Besser, B. P., Biernat, H. K., & Rijnbeek, R.P. 1990, Planet. Space Sci., 38, 411 [NASA ADS] [CrossRef] (In the text)
 Biskamp, D. 1994, Phys. Rep., 237, 179 [NASA ADS] [CrossRef] (In the text)
 Braginskii, S. I. 1965, Rev. Plasma Phys., 1, 205 [NASA ADS] (In the text)
 Craig, I. J. D. 2010, A&A, 515, A96 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Craig, I. J. D., & Henton, S. M. 1995, ApJ, 450, 280 [NASA ADS] [CrossRef] (In the text)
 Craig, I. J. D., & Litvinenko, Y. E. 2009, A&A, 501, 755 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Craig, I. J. D., & Litvinenko, Y. E. 2010, ApJ, 725, 886 [NASA ADS] [CrossRef] (In the text)
 Craig, I. J. D., Litvinenko, Y. E., & Senanayake, T. 2005, A&A, 433, 1139 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Heerikhuisen, J., Craig, I. J. D., & Watson, P. G. 2000, Geophys. Astrophys. Fluid Dyn., 93, 115 [NASA ADS] [CrossRef] (In the text)
 Hollweg, J. V. 1985, J. Geophys. Res., 90, 7620 [NASA ADS] [CrossRef] (In the text)
 Hollweg, J. V. 1986, ApJ, 306, 730 [NASA ADS] [CrossRef] (In the text)
 Hosking, R. J., & Marinoff, G. M. 1973, Plasma Phys., 15, 327 [NASA ADS] [CrossRef] (In the text)
 Jardine, M., Allen, H. R., Grundy, R. E., & Priest, E. R. 1992, J. Geophys. Res., 97, 4199 [NASA ADS] [CrossRef] (In the text)
 Litvinenko, Y. E. 2005, Sol. Phys., 229, 203 [NASA ADS] [CrossRef] (In the text)
 Litvinenko, Y. E., & Craig, I. J. D. 2000, ApJ, 544, 1101 [NASA ADS] [CrossRef] (In the text)
 Litvinenko, Y. E., & Craig, I. J. D. 2003, Sol. Phys., 218, 173 [NASA ADS] [CrossRef] (In the text)
 Montgomery, D. 1992, J. Geophys. Res., 97, 4309 [NASA ADS] [CrossRef] (In the text)
 Park, W., Monticello, D. A., & White, R. B. 1983, Phys. Fluids, 27, 137 [NASA ADS] [CrossRef] (In the text)
 Parker, E. N. 1957, J. Geophys. Res., 62, 509 [NASA ADS] [CrossRef] (In the text)
 Priest, E. R., & Forbes, T. 2000, Magnetic reconnection: MHD theory and applications (Cambridge Univ. Press) (In the text)
 Rüdiger, G., & Kitchatinov, L. L. 1997, Astron. Nach., 318, 273 [NASA ADS] [CrossRef] (In the text)
 Sonnerup, B. U. Ö., & Priest, E. R. 1975, Plasma Phys., 14, 283 [NASA ADS] [CrossRef] (In the text)
 Titov, V. S., & Priest, E. R. 1997, J. Fluid Mech., 348, 327 [NASA ADS] [CrossRef] (In the text)
All Figures
Fig. 1 Contour plot of J taken at the time of maximum current density at the origin. The parameters are α_{0} = 1, β_{0} = 0 and η = ν = 0.004. 

Open with DEXTER  
In the text 
Fig. 2 Arrow plot of v taken at the time of maximum current density at the origin. The parameters are the same as in Fig. 1, namely α_{0} = 1, β_{0} = 0 and η = ν = 0.004. For visual clarity only selected velocity vectors are displayed. 

Open with DEXTER  
In the text 
Fig. 3 The ycomponent of the magnetic field in the current layer in the case of headon merging. The parameters are the same as in Figs. 1 and 2. The solid line is the numerical magnetic field and the dashed line is the analytical solution from (13). 

Open with DEXTER  
In the text 
Fig. 4 Current sheet outflow speed comparisons. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER  
In the text 
Fig. 5 Current sheet inflow speed comparisons. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER  
In the text 
Fig. 6 Current sheet thickness comparisons. The dotted line shows η^{1/2} scaling. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER  
In the text 
Fig. 7 Viscous dissipation rate comparisons. Crosses, diamonds and circles refer to Braginskii, classical and ν = η control results respectively. The parameters are α_{0} = 1 and β_{0} = 0; ν = 0.004 for classical and Braginskii viscosity. 

Open with DEXTER  
In the text 
Fig. 8 The ratio W_{ν}/W_{η} for classical (diamonds) and Braginskii (crosses) viscosity. The dotted line shows η^{ − 1/2} scaling. The parameters are α_{0} = 1, β_{0} = 0 and ν = 0.004. 

Open with DEXTER  
In the text 
Fig. 9 Current contour plot for sheared reconnection in the Braginskii case at time of maximum current density at the origin. The parameters are α_{0} = 1, β_{0} = 0.7, η = 0.004, and ν = 0.004. 

Open with DEXTER  
In the text 
Fig. 10 The ratio W_{ν}/W_{η} for classical (diamonds) and Braginskii (crosses) viscosity. The dotted line shows η^{ − 1/2} scaling. The parameters are α_{0} = 1, β_{0} = 0.7 and ν = 0.004. 

Open with DEXTER  
In the text 