A&A 494, 317-327 (2009)

DOI: 10.1051/0004-6361:200811106

**C. T. M. Clack - I. Ballai - M. S. Ruderman**

Solar Physics and Space Plasma Research Centre (*SP*^{2}*RC*), Department of Applied Mathematics, University of Sheffield, Hicks Building, Hounsfield Road, Sheffield, S3 7RH,
UK

Received 8 October 2008 / Accepted 21 November 2008

**Abstract**
*Aims.* In the approximation of linear dissipative magnetohydrodynamics (MHD), it can be shown that driven MHD waves in magnetic plasmas with high Reynolds number exhibit a near resonant behaviour if the frequency of the wave becomes equal to the local Alfvén (or slow) frequency of a magnetic surface. This behaviour is confined to a thin region, known as the dissipative layer, which embraces the resonant magnetic surface. Although driven MHD waves have small dimensionless amplitude far away from the resonant surface, this near-resonant behaviour in the dissipative layer may cause a breakdown of linear theory. Our aim is to study the nonlinear effects in Alfvén dissipative layer

*Methods.* In the present paper, the method of simplified matched asymptotic expansions developed for nonlinear slow resonant waves is used to describe nonlinear effects inside the Alfvén dissipative layer.

*Results.* The nonlinear corrections to resonant waves in the Alfvén dissipative layer are derived, and it is proved that at the Alfvén resonance (with isotropic/anisotropic dissipation) wave dynamics can be described by the linear theory with great accuracy.

**Key words: **magnetohydrodynamics (MHD) - methods: analytical - Sun: atmosphere - Sun: oscillations

Magnetic fields are ubiquitous in solar and space plasmas. For regions where plasma-beta (the ratio of the kinetic and magnetic pressures) is less than one, magnetism controls the dynamics, topology and thermal state of the plasma. The magnetic field in the solar atmosphere is not dispersed, but it tends to accumulate in thinner or thicker entities often approximated as magnetic flux tubes. These magnetic flux tubes serve as an ideal medium for guided wave propagation.

One particular aspect of the solar physics that has attracted much attention since the 1940s is the very high temperature of the solar corona compared with the much cooler lower regions of the solar atmosphere requesting the existence of some mechanism(s) that keeps the solar corona hot against the radiative cooling. One of the possible theories proposed is the transfer of omnipresent waves' energy into thermal energy by resonant absorption or resonant coupling of waves (see e.g. Goossens et al. 1995; Poedts et al. 1990; Sakurai et al. 1991).

Waves which were initially observed sporadically mainly in radio wavelengths (see e.g., Aschwanden et al. 1992; Kai & Takayanagi 1973) are now observed in abundance in all wavelengths, especially in (extreme) ultraviolet (see e.g., Aschwanden et al. 1999; King et al. 2003; Nakariakov et al. 1999; Erdélyi & Taroyan 2008; Mariska et al. 2008; DeForest & Gurman 1998; Robbrecht et al. 2001). Since the plasma is non-ideal, waves can lose their energy through transport processes, however, the time over which the waves dissipate their energy is far too long. In order to have an effective and localized energy conversion, the plasma must exhibit transversal inhomogeneities relative to the direction of the ambient magnetic field. It was recognised a long time ago that solar and space plasmas are inhomogeneous, with physical properties varying over length scales much smaller than the scales determined by the gravitational stratification. Homogenous plasmas have a spectrum of linear eigenmodes which can be divided into slow, fast and Alfvén subspectra. The slow and fast subspectra have discrete eigenmodes whereas the Alfvén subspectrum is infinitely degenerated. When an inhomogeneity is introduced the three subspectra are changed. The infinite degeneracy of the Alfvén point spectrum is lifted and replaced by the Alfvén continuum along with the possibility of discrete Alfvén modes occurring, the accumulation point of the slow magnetoacoustic eigenvalues is spread out into the slow continuum and a number of discrete slow modes may occur, and the fast magnetoacoustic point spectrum accumulates at infinity (see e.g., Goedbloed 1975,1984).

According to the accepted wave theories, effective energy transfer between an energy carrying wave and the plasma occurs if the frequency of the wave matches one of the frequencies in the slow or
Alfvén continua, i.e. at the slow or Alfvén resonances. The Alfvén resonance has been more frequently associated with heating of coronal structures given the low- regime of the solar corona.
Nevertheless, slow resonance cannot be ruled out as an additional source of energy transfer. From a mathematical point of view, a resonance is equivalent to regular singular points in the equations
describing the dynamics of waves, but these singularities can be removed by, e.g., dissipation. Recently resonant absorption has acquired a new applicability when the observed damping of waves and
oscillations in coronal loops has been attributed to resonant absorption. Hence, resonant absorption has become a fundamental constituent block of one of the newest branches of solar physics,
called *coronal seismology* (see e.g., Goossens et al. 2002; Nakariakov et al. 1999; Arregui et al. 2007; Goossens et al. 2008; Terradas et al. 2008; Ballai et al. 2008; Ruderman & Roberts 2002) when applied to corona and solar magneto-seismology
when applied to the entire coupled solar atmosphere (see e.g., Verth et al. 2007; Erdélyi et al. 2007; Verth 2007).

Given the complexity of the mathematical approach, most theories describing resonant waves are limited to the linear regime. Perturbations, in these theories, are considered to be just small deviations from an equilibrium despite the highly nonlinear character of MHD equations describing the dynamics of waves and the complicated interaction between waves and plasmas. Initial numerical investigations of resonant waves in a nonlinear limit (see e.g., Ofman & Davila 1995) unveiled that the account of nonlinearity introduces new physical effects which cannot be described in the linear framework.

The first attempts to describe the nonlinear resonant waves analytically appeared after the papers by Ruderman et al. (1997b,a) which were followed by further analysis by, e.g., Ballai et al. (1998); Ballai & Erdélyi (1999); Ruderman (2000); Clack & Ballai (2008); however, all these papers focused on the slow resonant waves only. These studies revealed that nonlinearity does affect the absorption of waves. In addition, the absorption of wave momentum generates a mean shear flow which can influence the stability of resonant systems.

The present paper is the first analytical study on the nonlinear resonant Alfvén wave, where we obtain governing equations using techniques made familiar from previous studies on nonlinear slow
resonant MHD waves. Before embarking on the actual derivation, let us carry out a qualitative discussion. First of all we should point out that in plasmas with high Reynolds numbers (as in the solar corona) efficient dissipation only operates in a thin layer embracing the resonant surface. This layer is called *the dissipation layer*. This restriction on the effect of dissipation
makes the problem more tractable from a mathematical point of view, as outside the dissipative layer the dynamics of waves is described by the ideal MHD. Dissipation is a key ingredient of the problem of resonance. As it was mentioned earlier dissipation removes singularities in mathematical solutions. From a physical point of view dissipation is important as it is the mechanism which relaxes the accumulation of energy at the resonant surface and eventually
contributes to the global process of heating.

It is important to stress that the choice of dissipation has to be related to the very physics which is described as different waves are sensitive to different dissipative mechanisms. Due to the dominant role of the magnetic field in the solar corona, transport processes are highly anisotropic. Possible dissipation mechanisms acting in coronal structures can be described within the framework of Braginskii's theory (Braginskii 1965) as it was shown in applications by, e.g. Erdélyi & Goossens (1995); Mocanu et al. (2008); Ofman & Davila (1995). Alfvén waves are incompressible and transversal (in polarization), therefore, it is sensible to adopt shear viscosity and magnetic resistivity. Despite both transport processes being described by rather small coefficients, the net effect of dissipation can be increased considerably when the dissipative coefficients are multiplied by large transversal gradients.

The paper is organized as follows. In the next section we introduce the fundamental equations and discuss the main assumptions. In Sect. 3, we derive the governing equation for wave dynamics inside the Alfvén dissipative layer. Section 4 is devoted to calculating the nonlinear corrections at Alfvén resonance. Finally, in Sect. 5 we summarise and draw our conclusions, pointing out a few applications and further studies to be carried out in the future.

For describing mathematically the nonlinear resonant Alfvén waves we use the visco-resistive MHD equations. In spite of the presence of dissipation we use the adiabatic equation as an approximation of the energy equation. Numerical studies by Poedts et al. (1994) in linear MHD have shown that dissipation due to viscosity and finite electrical conductivity in the energy equation does not alter significantly the behaviour of resonant MHD waves in the driven problem.

When the product of the ion (electron) gyrofrequency, , and the ion (electron) collision time, , is much greater than one (as in the solar corona) the viscosity and finite electrical conductivity become anisotropic and viscosity is given by the Braginskii viscosity tensor (see Appendix A). The components of the viscosity tensor that remove the Alfvén singularity are the shear components. The parallel and perpendicular components of anisotropic finite electrical conductivity only differ by a factor of 2, therefore, we will consider only one of them without loss of generality.

The dynamics of waves in our model is described by the visco-resistive MHD equations

In Eqs. (1)-(5) and are the velocity and magnetic induction vectors, the plasma pressure, the density, the coefficient of magnetic diffusivity, the adiabatic exponent, and the magnetic permeability of free space. In addition, is the Braginskii viscosity (see Appendix A for full details). Note that even though anisotropy has been considered the Hall term has been neglected. We neglect the Hall term from the induction equation, which can be of the order of diffusion term in the solar corona, because the largest Hall terms in the perpendicular direction relative to the ambient magnetic field identically cancel. The components of the Hall term in the normal and parallel directions relative to the ambient magnetic field have no effect on the dynamics of Alfvén waves in dissipative layers, hence these too are neglected. For full details on the Hall term and the reasoning behind neglecting it, we refer to Appendix B.

We adopt Cartesian coordinates *x*, *y*, *z* and limit our analysis to a static background equilibrium (
). We assume that all equilibrium quantities depend on *x* only. The equilibrium magnetic field, ,
is unidirectional and lies in the *yz*-plane. The equilibrium quantities must satisfy the condition of total pressure balance,

For simplicity we assume that the perturbations of all quantities are independent of

The definition of the resonant position can now be written mathematically as

where is the Alfvén speed defined as

In addition we introduce the sound and cusp speeds as

In what follows we can take without loss of generality. The perturbations of the physical quantities are defined by

where

The dominant dynamics of resonant Alfvén waves, in linear MHD, resides in the components of the perturbed magnetic field and velocity that are perpendicular to the equilibrium magnetic field and to the *x*-direction. This dominant behaviour is created by an *x*^{-1} singularity in the spatial solution of these quantities at the Alfvén resonance (Goossens & Ruderman 1995; Sakurai et al. 1991); these variables are known as *large variables*. The *x*-component of velocity, the components of magnetic field normal and parallel to the equilibrium magnetic field, plasma pressure and density are also singular, however, their singularity is proportional to .
In addition, the quantities *P* and the components of
and
that are parallel to the equilibrium magnetic field are regular; all these variables are called *small variables*.

To make the mathematical analysis more concise and the physics more transparent we define the components of velocity and magnetic field that are in the *yz*-plane and are either parallel or perpendicular to the equilibrium magnetic field:

(12) |

where

Let us introduce the characteristic scale of inhomogeneity,
.
The classical viscous Reynolds number, ,
and the magnetic Reynolds number, ,
are defined as

where is a characteristic value of , and is the shear viscosity coefficient (see Appendix A). These two numbers determine the importance of viscosity and finite electrical conductivity. We introduce the total Reynolds number as

The aim of this paper is to study the nonlinear behaviour of driven Alfvén resonant waves in the dissipative layer. We are not interested in MHD waves that have large amplitude everywhere and require a nonlinear description in the whole space. We focus on waves that have small dimensionless amplitude far away from the ideal Alfvén resonant point

In nonlinear theory, when studying resonant behaviour in the dissipative layer we must scale the dissipative coefficients (see e.g., Ruderman et al. 1997b; Ballai et al. 1998; Clack & Ballai 2008). The
general scaling to be applied is

Linear theory predicts that the characteristic thickness of the dissipative layer, , is of the order of and we assume that this is true in the nonlinear regime, too. Hence, we must introduce a stretching transversal coordinate, , in the dissipative layer defined as

We can rewrite Eqs. (1)-(5) in the scalar form as

In the above equations and .

Equations (17)-(26) will be used in the following sections to derive the governing equation for the resonant Alfvén waves inside the dissipative layer and to find the nonlinear corrections.

In order to derive the governing equation for wave motions in the Alfvén dissipative layer we employ the method of matched asymptotic expansions (Nayfeh 1981; Bender & Országh 1991). This method requires to
find the so-called *outer* and *inner* expansions and then match them in the overlap regions. This nomenclature is ideal for our situation. The outer expansion corresponds to the
solution outside the dissipative layer and the inner expansion corresponds to the solution inside the dissipative layer. A simplified version of the method of matched asymptotic expansions,
developed by Ballai et al. (1998), is adopted here.

The typical largest *quadratic* nonlinear term in the system of MHD equations is of the form
while the typical dissipative term is of the form
,
where *g* is any ``large'' variable. Linear theory predicts that ``large'' variables have an ideal singularity *x*^{-1} in the vicinity of *x*=0. This implies that the ``large'' variables have dimensionless amplitudes in the dissipative layer of the order of
.
It is now straightforward to estimate the ratio of a typical quadratic nonlinear and dissipative term,

where the quantity can be considered as the

Far away from the dissipative layer the amplitudes of perturbations are small, so we use linear ideal MHD equations in order to describe the wave motion. The full set of nonlinear dissipative MHD equations are used for describing wave motion *inside* the dissipative layer where the amplitudes can be large. We, therefore, look for solutions in the form of asymptotic expansions. The equilibrium quantities change only slightly across the dissipative layer so it is possible to approximate them by the first non-vanishing term in their Taylor series expansion with respect to *x*. Similar to linear theory, we assume that the expansions of equilibrium quantities are
valid in a region embracing the ideal resonant position, which is assumed to be much wider than the dissipative layer. This implies that there are two overlap regions, one to the left and one to the
right of the dissipative layer, where both the outer (the solution to the linear ideal MHD equations) and inner (the solution to the nonlinear dissipative MHD equations) solutions are valid. Hence, both solutions must coincide in the overlap regions which provides
the matching conditions.

Before deriving the nonlinear governing equation we ought to make a note. In linear theory, perturbations of physical quantities are harmonic functions of and their mean values over a period are zero. In nonlinear theory, however, the perturbations of variables can have non-zero mean values as a result of nonlinear interaction of different harmonics. Due to the absorption of wave momentum, a mean shear flow is generated outside the dissipative layer (Ofman & Davila 1995). This result is true for our analysis also, however, due to the length of this study we prefer to deal with this problem in a forthcoming paper.

We suppose that nonlinearity and dissipation are of the same order so we have
,
i.e.
.
We can, therefore, substitute
for *R* in Eq. (15) to rescale viscosity and
finite electrical resistivity as

We do not rewrite the MHD equations as they are easily obtained from Eqs. (17)-(26) by substitution of Eq. (28).

The first step in our description is the derivation of governing equations outside the dissipative layer where the dynamics is described by ideal (
)
and linear MHD. The linear
form of Eqs. (17)-(26) can be obtained by assuming a regular expansion of variables of the form

and collect only terms proportional to the small parameter . This leads to a system of linear equations for the variables with superscript ``1''. All variables can be eliminated in favour of

where

The quantities

Since Eq. (30) has regular singular points, the solutions can be obtained in terms of Fröbenius series with respect to

The coefficient functions depending on in the above expansions are, generally, different for

As the characteristic scale of dissipation is of the order of
and we have assumed that
we obtain that the thickness of the dissipative layer is
,
implying the introduction of a new stretched variable to replace the transversal coordinate in the dissipative layer, which is defined as
.
Again, for brevity, Eqs. (17)-(26) are not rewritten as they can be obtained by the substitution of

for all derivatives. The equilibrium quantities still depend on

where

We seek the solution to the set of equations obtained from Eqs. (17)-(26) by the substitution of
into variables in the form of power series of .
These equations contain powers of
,
so we use this quantity as an expansion parameter. To derive the
form of the inner expansions of different quantities we have to analyze the outer solutions. First, since
and *P* are regular at *x*=0 we can write their inner expansions in the
form of their outer expansions Eq. (29). The amplitudes of large variables in the dissipative layer are of the order of
,
so the inner expansion of the variables
and is

The quantities

We now substitute the expansion (29) for *P*, *u*, *b*_{x},
,
,
*p* and
and the expansion given by Eq. (42) for
and
into the set of equations obtained from Eqs. (17)-(26) after substitution of
.
The first order approximation (terms proportional to ), yields a linear homogeneous system of equations for the terms with superscript ``1''. The important result that follows from this set of equations is that

that is

(44) |

(47) |

(48) |

All equilibrium quantities are calculated at

In the second order approximation we only use the expressions obtained from Eqs. (19) and (22). Employing Eqs. (43)-(49), we replace the variables in the second order approximation which have superscript ``1''. The equations obtained in the second order are

(50) |

Once the variables with superscript ``2'' have been eliminated from the above two equations, the governing equation for resonant Alfvén waves inside the dissipative layer is derived as

where

It is clear that Eq. (52) does not contain nonlinear terms despite considering the full MHD system of equations. This result is in stark contrast with the results obtained for nonlinear slow resonance where the governing equation was found to be always nonlinear (see e.g., Ruderman et al. 1997b; Ballai et al. 1998; Clack & Ballai 2008). The governing Eq. (52) suggests that resonant Alfvén waves can be described by the linear theory unless their amplitudes inside the dissipative layer is of the order of unity.

As the quadratic nonlinear terms cancel each other out, it is natural to take into account cubic nonlinearity (the system of MHD equations contain cubic nonlinear terms), where the nonlinearity
parameter is defined as

Despite the higher order nonlinearity the governing equation is similar to the equation derived for quadratic nonlinearity (52). These results require finding an explanation to the linear behaviour of waves inside the dissipative layer. The following section will be devoted to the study of nonlinear corrections in the Alfvén dissipative layer.

Since we have assumed that waves have small dimensionless amplitude outside the dissipative layer, we will concentrate only on the solutions *inside* the dissipative layer.

In our analysis we use the assumptions and equations presented in Sect. 2, however, we will not impose any relation between
and *R*. Equations (15) and (16) will be used to define the scaled dissipative coefficients and stretching transversal coordinate in the dissipative layer. For simplicity we denote
.
This means that our scaled dissipative coefficients and stretched transversal coordinate become

The first step to accomplish our task is to rewrite Eqs. (17)-(26) by substituting

All equilibrium quantities (which are still dependent on

The substitution of Eqs. (41), (55) and (56) will transform Eqs. (17)-(26) into

The only condition we need to impose when deriving the nonlinear corrections to resonant Alfvén waves in the dissipative layer is imported from the linear theory which predicts that in the dissipative layer ``large'' variables have dimensionless amplitude of the order of (see Sect. 3). We assume that the dimensionless amplitudes of the linear approximation of ``large'' variables ( and ) in the dissipative layer are small, so that

This condition ensures that the oscillation amplitude remains small inside the dissipative layer. From a naive point of view the linear theory is applicable as soon as the oscillation amplitude is small. The example of slow resonant waves clearly shows that this is not the case. The nonlinear effects become important in the slow dissipative layer as soon as , i.e. as soon as the oscillation amplitude in the dissipative layer, which is of the order of , is of the order of . For example, in the corona perturbations with dimensionless amplitudes less than 10

We now assume that all perturbations can be written as a regular asymptotic expansion of the form

where represents the equilibrium value. Substitution of expansion (67) into the system (57)-(65) leads to a system of equations which contains the small parameter . This observation inspires us to look for the solution in the form of expansions with respect to . In order to cast large and small variables in this description we are going to use the following expansion for small variables (

while large variables ( and ) will be expanded according to

The bar notation is used here to distinguish between these expansions and the expansions used in the previous section. From this point on we drop the bar notation.

Substituting Eqs. (68), (69) into the system (57)-(65), taking terms proportional to and then only retaining terms with the lowest power of , results in the set of linear equations

In these equations all equilibrium quantities are calculated at

Using these equations we can express all dependent variables in terms of *u*_{1}^{(1)},
and *P*_{1}^{(1)},

It follows from Eq. (71) that

Finally, we obtain the relation between

Note that Eqs. (79)-(85) are formally identical to Eqs. (43)-(49) for the linear approximation in Sect. 3. This is not surprising as both methods are designed to replicate linear theory in the first order approximation.

Once the first order terms are known we can proceed to derive the second and third order approximations with respect to
(i.e. terms from the expansion of Eqs. (57)-(65) that are proportional to
and
,
respectively). First, we write out the second order approximations and substitute for all first order terms (i.e. terms of the form *f*_{1}^{(1)}) using Eqs. (79)-(85). Secondly, we find (by solving the inhomogeneous system) the expansions of second order terms (terms with *subscript* ``2''). Thirdly, we derive the second order relations between all variables, similar to the ones obtained
in the first order approximation.

The equations representing the second order approximation with respect to (with variables in the first order substituted) are

It is clear that nonlinear terms appear from this order of approximation and they are expressed in terms of variables obtained in the first order.

The analysis of the system of Eqs. (86)-(94) reveals that the expansions with respect to
has to be written in the form

for

for , ,

Here we need to make a note. It follows from Eqs. (95) and (96) that the ratio of
to
is of the order of
,
and the same is true for
,
and *p*_{2}. It seems to be inconsistent with the regular perturbation method where it is assumed that the next
order approximation is always smaller than the previous one. However, this problem is only apparent. To show this we need to clarify the exact mathematical meaning of the statement ``in the
asymptotic expansion each subsequent term is much smaller than the previous one''. To do this we introduce the nine-dimensional vector
and consider it as an element of a Banach space. The norm in this space can be introduced in different ways. One possibility is

where

Once the expansions (95) and (96) are substituted into Eqs. (86)-(94), we can express the variables in this order of approximation as

For the total pressure we obtain that

In addition, we obtain that the equation which determines

Since the large variables in this order of approximation are and , we can deduce that the linear order of approximation of resonant Alfvén waves in the dissipative layer excite

We now calculate the third order approximation with respect to . On analyzing the third order system of equations we deduce that the large variables in this order of approximation are Alfvénic ( , ), so we only need the perpendicular components of momentum and induction equations given by Eqs. (59) and (62), respectively.

First, since some of the first order approximation terms contribute to the third order approximation in integral form we must introduce a new notation

The third order approximation of the perpendicular component of momentum is

while the perpendicular component of magnetic induction equation is

Note that in obtaining the third order approximations we have employed all the relations we have for variables in the first and second order of approximation.

Equations (104) and (105) clearly show that the nonlinear terms on the right-hand sides do not cancel. This implies that the expansion of
and
should be of the form

We should state, for completeness, that if we derive the third order approximation for all the Eqs. (57)-(65) we obtain the expansions for

The expansions calculated for all the variables can now be collected together and we can write the expansions for ``large'' and ``small'' variables in the dissipative layer when studying resonant Alfvén waves. Large variables ( and ) have the expansion

and the expansion of small variables (

Since Eq. (66) is the only condition enforced in the dissipative layer, we can state that

Therefore, since both Eqs. (108) and (109) contain successive higher powers of the parameter we can deduce that, considering Eq. (110), higher orders of approximation of large and small variables become increasingly insignificant in comparison to the linear order of approximation, so resonant Alfvén waves in the dissipative layer can be described accurately by linear theory if condition (66) is satisfied.

In the present paper we have investigated the nonlinear behaviour of resonant Alfvén waves in the dissipative layer in one-dimensional planar geometry in plasmas with anisotropic dissipative coefficients, a situation applicable to solar coronal conditions. The plasma motion outside the dissipative layer is described by the set of linear, ideal MHD equations. The wave motion inside the dissipative layer is governed by Eq. (52). This equation is linear, despite taking into consideration (quadratic and cubic) nonlinearity. The Hall terms of the induction equation in the perpendicular direction relative to the ambient magnetic field cancel each other out.

The nonlinear corrections were calculated to explain why Eq. (52), describing the nonlinear behaviour of wave dynamics, is always linear. We found that, in the second order of approximation, magnetoacoustic modes are excited by the perturbations of the linear order of approximation. These secondary waves act to counteract the small pressure and density variations created by the first order terms. In addition, these waves are not resonant in the Alfvén dissipative layer. In the third order approximation the perturbations become Alfvénic, however, these perturbations are much smaller than those in the linear order of approximation. Equations (108) and (109) describe the expansion of large and small variables, respectively, and demonstrate that all higher order approximations of both large and small variables at the Alfvén resonance are smaller than the linear order approximation, provided condition (66) is satisfied. This condition ensures that the oscillation amplitude remains small inside the dissipative layer. From a naive point of view the linear theory is applicable as soon as the oscillation amplitude is small. The example of slow resonant waves clearly shows that this is not the case. The nonlinear effects become important in the slow dissipative layer as soon as , i.e. as soon as the oscillation amplitude in the dissipative layer, which is of the order of , is of the order of . We also found that any dispersive effect due to the consideration of ions' inertial length (Hall effect) is absent from the governing equation.

This calculation of nonlinear corrections to resonant Alfvén waves in dissipative layers allows us to apply the already well-known linear theory for studying resonant Alfvén waves in the solar corona with great accuracy, where the governing equation, jump conditions and the absorption of wave energy are already derived (see e.g., Goossens et al. 1995; Erdélyi 1998; Sakurai et al. 1991).

It is interesting to note that this work can be transferred to isotropic plasma rather easily. Shear viscosity, supplied by Braginskii's viscosity tensor (see Appendix A), acts exactly as isotropic viscosity. Therefore, replacing by in Eq. (52) provides the required governing equation for resonant Alfvén waves in isotropic plasmas. Moreover, the work on the nonlinear corrections presented in this paper is also unaltered by anisotropy. This implies that we can consider resonant Alfvén waves in dissipative layers throughout the solar atmosphere and still use linear theory if condition (66) is satisfied.

C.T.M.C. would like to thank STFC (Science and Technology Facilities Council) for the financial support provided. I.B. was financially supported by NFS Hungary (OTKA, K67746) and The National University Research Council Romania (CNCSIS-PN-II/531/2007).

In this Appendix, we shall derive the largest terms of Braginskii's viscosity tensor inside the Alfvén dissipative layer to be used to study the nonlinearity effects. Braginskii's viscosity tensor
comprises of five terms. Its divergence can be written as (Braginskii 1965)

Note that the terms proportional to , and in Eq. (A.1) describe viscous dissipation, while terms proportional to and are non-dissipative and describe the wave dispersion related to the finite ion gyroradius. They will be ignored in what follows. For simplicity in the body of the paper, we have taken .

The quantities ,
and
are given by

Here is the velocity, , is the unit tensor and indicates the dyadic product of two vectors. The superscript ``T'' denotes a transposed tensor.

The first viscosity coefficient,
,
(*compressional viscosity*) has the following approximate expression (see e.g., Ruderman et al. 2000)

where and

The viscosity described by the sum of the second and third terms in Eq. (A.1) is

First, we shall calculate the components of the compressional viscosity. We will use the notation of parallel and perpendicular components as defined in the paper. It is straightforward to obtain that

The shear viscosity, as stated above, is the sum of the second and third terms of Eq. (A.1). To evaluate these terms we use the approximate expression for and given by Eq. (A.7). As a result we obtain

Equations (A.8)-(A.13) are complicated, but we can simplify them further by taking the largest term only in each equation. For the viscosity in the parallel direction it would, at first, seem obvious that the largest term will be proportional to rather than . However, some of the variables proportional to have derivatives with respect to

Equations (A.14)-(A.16) give an appropriate approximation to Braginskii's viscosity tensor when studying nonlinear resonant Alfvén waves in dissipative layers. It is interesting to note that the terms in Eqs. (A.14)-(A.16) are identical to the largest terms when considering isotropic viscosity. Obviously, compressional viscosity cannot remove the Alfvén singularity since Eq. (A.9) is identically zero.

In this Appendix we will derive the components of the Hall term in the induction equation and show that neglecting the Hall effect at the Alfvén resonance is a good approximation for typical conditions throughout the solar atmosphere. The main reasons qualitatively are as follows. When we are in the lower solar atmosphere (e.g. solar photosphere) the Hall conduction is much smaller than the direct conduction since the product of the electron gyrofrequency, , and collision time, , is less than unity (see e.g., Priest 1984). For the upper atmosphere (e.g. chromosphere, corona), where the product is greater than unity, the Hall conduction has to be considered. However, when the Hall terms are derived, the largest terms in the perpendicular direction relative to the ambient magnetic field cancel leaving only higher order approximation terms which are far smaller than the direct conduction. As the dominant dynamics of resonant Alfvén waves in dissipative layer resides in the components of velocity and magnetic field perturbation in the perpendicular direction relative to the background magnetic field we can neglect the Hall conduction completely from the analysis without affecting the governing equation.

In order to estimate the relative importance of the Hall term and resistive term in the dissipative layer we follow the sophisticated analysis which was presented by Ruderman et al. (1997b) and Clack & Ballai (2008). We do not write down all the steps of the analysis, but rather give the salient points specific to the Hall effect at the Alfvén resonance.

Equations (29) and (42) provide the following estimations in the dissipative layer:

where still denotes the dimensionless amplitude of oscillations far away from the dissipative layer. The thickness of the dissipative layer divided by the characteristic scale of inhomogeneity is . This gives rise to the estimations

where

We now need to calculate the components of the vectors of the resistive term and the Hall term from the induction equation normal to the magnetic surfaces (the

where the dots indicate terms much smaller than those shown explicitly. With the aid of Eqs. (B.1)-(B.3) we obtain the ratios

For the Hall conduction to be significant in the direction of the dominant dynamics of resonant Alfvén waves (i.e. in the perpendicular direction) we must have . This is plausible for the solar upper atmosphere. If this condition holds, then we must consider the Hall term in the induction equation. However, if Eq. (B.8) is expanded using Eq. (29) from Sect. 3 we obtain

It should be noted that in deriving Eq. (B.13) we have used the assumption that . The terms inside the braces are of the same order as the direct conduction. Hence, they would be expected to appear in the governing equation. When substituting for

Equation (B.14) shows that the Hall term in the perpendicular component of induction is always smaller than the direct conduction in the solar atmosphere. The normal and parallel components of the Hall conduction are, in fact, larger than the perpendicular component. Nevertheless they play no role in derivation of the governing equation of resonant Alfvén waves in dissipative layer. The parallel component is the largest of the three components and this is to be expected as the Hall effect is strongest at right angles to the dominant wave motion. This is in complete agreement with the study on resonant slow waves by Clack & Ballai (2008) which found the largest Hall effect was in the perpendicular component of Hall conduction, which is at right angles to the dominant dynamics of resonant slow waves.

In summary, it is a good approximation to neglect the Hall term in the induction equation when studying resonant Alfvén waves in dissipative layer. This approximation holds throughout the entire solar atmosphere.

- Arregui, I., Andries, J., van Doorsselaere, T., Goossens, M., & Poedts, S. 2007, A&A, 463, 333 [NASA ADS] [CrossRef] [EDP Sciences]
- Aschwanden, M. J., Bastian, T. S., & Gary, D. E. 1992, in Am. Astron. Soc., 180th AAS Meeting, #45.05, BAAS, 24, 802
- Aschwanden, M. J., Fletcher, L., Schrijver, C. J., & Alexander, D. 1999, ApJ, 520, 880 [CrossRef]
- Ballai, I., & Erdélyi, R. 1999, 8th SOHO Workshop: Plasma Dynamics and Diagnostics in the Solar Transition Region and Corona, 155
- Ballai, I., Douglas, M., & Marcu, A. 2008, A&A, 488, 1125 [NASA ADS] [CrossRef] [EDP Sciences]
- Ballai, I., Ruderman, M. S., & Erdélyi, R. 1998, Phys. Plasmas, 5, 252 [CrossRef]
- Bender, C. M., & Országh, S. A. 1991, Advanced Methods for Scientists and Engineers: Asymptotics (New York: Springer-Verlag)
- Braginskii, S. I. 1965, Rev. Plasmas Phys., 1, 205 [NASA ADS] (In the text)
- Clack, C. T. M., & Ballai, I. 2008, Phys. Plasmas, 15, 082310 [NASA ADS] [CrossRef] (In the text)
- DeForest, C. E., & Gurman, J. B. 1998, ApJ, 501, L217 [NASA ADS] [CrossRef]
- Erdélyi, R. 1998, Sol. Phys., 180, 213 [CrossRef]
- Erdélyi, R., & Taroyan, Y. 2008, ApJ, 489
- Erdélyi, R., & Goossens, M. 1995, A&A, 294, 575 [NASA ADS]
- Erdélyi, R., Pintér, B., & Malins, C. 2007, Astron. Nachr., 328, 305 [NASA ADS] [CrossRef]
- Goedbloed, J. P. 1975, Phys. Fluids, 18, 1258 [CrossRef]
- Goedbloed, J. P. 1984, Physica D, 12, 107 [NASA ADS] [CrossRef]
- Goossens, M., & Ruderman, M. S. 1995, Phys. Scr., 60, 171
- Goossens, M., Ruderman, M. S., & Hollweg, J. V. 1995, Sol. Phys., 157, 75 [NASA ADS] [CrossRef]
- Goossens, M., Andries, J., & Aschwanden, M. J. 2002, A&A, 394, L39 [NASA ADS] [CrossRef] [EDP Sciences]
- Goossens, M., Arregui, I., Ballester, J., & Wang, T. 2008, A&A, 484, 851 [NASA ADS] [CrossRef] [EDP Sciences]
- Kai, K., & Takayanagi, A. 1973, Sol. Phys., 29, 461 [NASA ADS] [CrossRef]
- King, D. B., Nakariakov, V. M., Deluca, E. E., Golub, L., & McClements, K. G. 2003, A&A, 404
- Mariska, J., Warren, H., Williams, D., & Watanabe, T. 2008, ApJ, 681
- Mocanu, G., Marcu, A., Ballai, I., & Orza, B. 2008, Astron. Nachr, 88, 789
- Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., & Davila, J. M. 1999, Science, 285, 862 [NASA ADS] [CrossRef]
- Nayfeh, A. H. 1981, Introduction to Perturbation Techniques (New York: Wiley-Interscience)
- Ofman, L., & Davila, J. M. 1995, J. Geophys. Res., 100, 23427 [NASA ADS] [CrossRef] (In the text)
- Poedts, S., Belien, A. J. C., & Goedbloed, J. P. 1994, Sol. Phys., 151, 271 [NASA ADS] [CrossRef] (In the text)
- Poedts, S., Goossens, M., & Kerner, W. 1990, ApJ, 360, 279 [NASA ADS] [CrossRef]
- Priest, E. R. 1984, Solar Magnetohydrodynamics (Berlin: Springer) (In the text)
- Robbrecht, E., Verwichte, E., Berghmans, D., et al. 2001, A&A, 370, 591 [NASA ADS] [CrossRef] [EDP Sciences]
- Ruderman, M. S. 2000, J. Plasma Phys., 63, 43 [NASA ADS] [CrossRef] (In the text)
- Ruderman, M. S., Goossens, M., & Hollweg, J. V. 1997a, Phys. Plasmas, 4, 91 [NASA ADS] [CrossRef]
- Ruderman, M. S., Hollweg, J. V., & Goossens, M. 1997b, Phys. Plasmas, 4, 75 [NASA ADS] [CrossRef]
- Ruderman, M. S., Oliver, R., Erdélyi, R., Ballester, J. L., & Goossens, M. 2000, A&A, 354, 261 [NASA ADS] (In the text)
- Ruderman, M. S., & Roberts, B. 2002, ApJ, 577, 475 [NASA ADS] [CrossRef]
- Sakurai, T., Goossens, M., & Hollweg, J. V. 1991, Sol. Phys., 133, 227 [NASA ADS] [CrossRef]
- Terradas, J., Arregui, I., Oliver, R., et al. 2008, ApJ, 679, 1611 [NASA ADS] [CrossRef]
- Verth, G. 2007, Astron. Nachr., 328, 764 [NASA ADS] [CrossRef]
- Verth, G., van Doorsselaere, T., Erdélyi, R., & Goossens, M. 2007, A&A, 475, 341 [NASA ADS] [CrossRef] [EDP Sciences]

Copyright ESO 2009