CORONAL HEATING BY THE RESONANT ABSORPTION OF ALFVEN WAVES - THE EFFECT OF VISCOUS STRESS TENSOR

The time-dependent linearized MHD equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tenser (Braginskii 1965) is included with the five parameters eta(i), i = 0-4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(i) are eta(0) = 0.65 g cm(-1) s(-1), eta(1) = 3.7 x 10(-12) g cm(-1) s(-1), eta(2) = 4 eta(1), eta(3) = 1.4 x 10(-6) g cm(-1) s(-1), eta(4) = 2 eta(3), with n = 10(10) cm(-3), T = 2 x 10(6) K, and B = 100 G. When the Lundquist number S = 10(4) and R(1) much greater than S(where R(1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(r) is 0.22a (where a is the density gradient length scale) and d(r) similar to S--1/3. When S much greater than R(1) the shear viscous dissipation layer width d(v) scales as R(1)(-1/3) The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(1)(1/3) It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(0) = 560) is negligible compared to the resistive and the shear viscous (R(1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximate to 10(14) and R(1) x 10(14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.