Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on Wave Propagation

*Ravinder Goyal and R.P. Sharma*

## **Abstract**

The damping phenomenon is studied due to the collisions of ions and neutral particles and Landau approach on the turbulent spectra of kinetic Alfvén wave (KAW) in magnetized plasma which is inhomogeneous as well. The localization of waves is largely affected by inhomogeneities in plasma which are taken in transverse as well as parallel directions to the ambient magnetic field. There is significant effect of damping on the wave localization and turbulent spectra. Numerical solutions of the equations governing kinetic Alfvén waves in the linear regime give the importance of wave damping phenomena while retaining the effects of Landau (collisionless) damping and ion-neutral collisional damping. A comparative study of the two damping effects reveals that the Landau damping effect is more profound under similar plasma conditions.

**Keywords:** kinetic Alfvén wave, inhomogeneous plasma, ion-neutral collisional damping, Landau damping, laboratory plasma, turbulence

## **1. Introduction**

As far as enormous space plasma phenomena like solar wind turbulence, acceleration of solar wind, heating of solar coronal loops, solar flares, etc. are concerned, kinetic Alfvén wave (KAW) has a vital role to play [1]. The dispersion characteristics of KAW make it distinguishable from its parent Alfvén wave [2]. When Alfvén wave propagation develops a large value for wave number perpendicular ð Þ *k*<sup>⊥</sup> to the ambient magnetic field, then it gives rise to kinetic Alfvén wave (KAW). This mode conversion may lead to transportation of large amount of electromagnetic energy across geomagnetic field lines [3]. In Tokamak/ITER plasma where KAW plays an important role, the anomalous diffusion coefficient which depends upon the turbulence level is inversely related to the confinement time. In the scenario where the electric fields are nonstationary, this wave plays a vital role in the phenomenon of particle acceleration [4]. The acceleration and heating [5] of plasma particles are dependent on perpendicular wave number and is attributed to KAW carrying finite perturbations parallel to electric field. Also, the plasma heating is caused by dissipation of turbulence and is essential to interpret the observations of most astronomical [6] and laboratory [7] systems. The turbulent behavior of KAW

can be thought of as the probable mechanism to describe the magnetic fluctuations which are observed near earth's magnetopause [8]. The theoretical and experimental studies [9, 10] support the importance of KAW in solar wind turbulence energy cascade and particle heating. Earth's magnetopause is also known to have a vital role played by KAW in the mechanism of stochastic ion heating [11]. Nykyri et al. [12] observed the turbulent spectra in high altitude cusp and revealed that break in observed spectra might be due to damping of obliquely propagating KAW. Scime [17] to experimentally study the collisional impact due to collisions between ion and neutral particles on the KAW propagation in plasma. Recently, Goyal and Sharma [35] have studied the effect of ion-neutral collisional damping on KAW

*Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on…*

The present research deals with a comparative study of the two damping phe-

(collisionless) damping on the propagation of KAW which have been studied separately by Sharma et al. [25] and Goyal and Sharma [35] in the light of inhomogeneous magnetized plasma. This comparison is very much important with respect to scrape-off regions of the fusion energy devices like Tokamak and ITER. The plasma inhomogeneity is taken into account by considering spatial (in x-z plane) inhomogeneities having scale lengths both in longitudinal and transverse directions to background magnetic field. Distinct model equations for KAW have been taken by considering damping factors corresponding to collisional and collisionless damping, and the results have been drawn, respectively, using two fluid simulation

The governing equations of KAW come into their present form by taking into account the propagation of the wave in two dimensions (x-z surface), which

the component of the wave vector in perpendicular (parallel) direction to ambient

implies the same in parallel direction. Initially, the wave field profile is considered to follow hollow Gaussian behavior [36]. As the current numerical study is inspired from the experimental technique [17] where the plasma environment was developed experimentally to study the wave propagation, we have used the same laboratory parameters in our model which are as follows: *B0* ≈ 560 G, plasma density

, neutral particle density *nn* <sup>≈</sup> <sup>5</sup> � <sup>10</sup><sup>12</sup> cm�<sup>3</sup>

Using the well-established methods [26–28], Maxwell's equations and Faraday's

mental values, the other characteristic plasma parameters may be obtained like:

law, the two dynamical equations for finite amplitude KAW incorporating the effects of ion-neutral collisional damping [35] and Landau damping [25] can,

*<sup>υ</sup>ti* <sup>≈</sup> 2.445 � <sup>10</sup><sup>5</sup> cm/s, <sup>ω</sup>ci = 2.68 � 106 rad/s, *<sup>λ</sup><sup>e</sup>* = 0.1 cm, *<sup>ρ</sup><sup>i</sup>* <sup>≈</sup> 0.09 cm,

*teλ*2 *<sup>e</sup>* <sup>þ</sup> *<sup>υ</sup>*<sup>2</sup> *Aρ*<sup>2</sup> *i <sup>∂</sup>*<sup>4</sup>*By*

2ð Þ *nnυti=niυtn* <sup>1</sup> <sup>þ</sup> <sup>ω</sup><sup>2</sup>*=ν*<sup>2</sup> *ni*

*∂*<sup>4</sup>*By*

<sup>0</sup> along z direction. Also, the plasma having background density *n*<sup>0</sup>

*<sup>o</sup>* ¼ *koxx*^ þ *kozz*^where *k0x(k0z)* symbolizes

!

, ion density

<sup>≈</sup> 0.002, *<sup>υ</sup>te* <sup>≈</sup> 1.1 � <sup>10</sup><sup>8</sup> cm/s,

<sup>A</sup> exp *x*<sup>2</sup>

*=L*<sup>2</sup> *<sup>x</sup>* <sup>þ</sup> *<sup>z</sup>*<sup>2</sup> *=L*<sup>2</sup> *z*

*<sup>∂</sup>*<sup>2</sup>

*By ∂z*<sup>2</sup>

(1)

<sup>0</sup> while *Lz*

*<sup>x</sup>* � *<sup>z</sup>*<sup>2</sup>*=L*<sup>2</sup> *z* , where *Lx* implies

,*T*<sup>e</sup> = 81,200 K,*Tn* ≈ *Ti* = 2900 K. Using these experi-

*<sup>s</sup> =υ*<sup>2</sup> *A*

, and *k0z* ≈ 0.16 cm�<sup>1</sup> for ω/ωci = *f/fci* = 0.9.

*<sup>∂</sup>x*<sup>2</sup>*∂z*<sup>2</sup> <sup>þ</sup> *<sup>υ</sup>*<sup>2</sup>

*∂z*<sup>2</sup>*∂t*<sup>2</sup>

!

is considered to be inhomogeneous in x as well as in z directions and hence is

the inhomogeneity scale length perpendicular to ambient magnetic field *B*

nomena discussed above viz. ion-neutral collisional damping and Landau

**2. Equations governing the propagation of KAW**

supposed to follow the density profile *<sup>n</sup>*<sup>0</sup> exp �*x*<sup>2</sup>*=L*<sup>2</sup>

*<sup>υ</sup><sup>A</sup>* <sup>≈</sup> 2.5 � <sup>10</sup><sup>7</sup> cm/s, *cs* <sup>¼</sup> <sup>1</sup>*:*<sup>3</sup> � <sup>10</sup><sup>6</sup> cm/s, *<sup>β</sup>* <sup>¼</sup> *<sup>c</sup>*<sup>2</sup>

*∂*<sup>4</sup>*By <sup>∂</sup>x*<sup>2</sup>*∂t*<sup>2</sup> � *<sup>υ</sup>*<sup>2</sup>

1 þ

þ *υ*2 A ω2 *ci*

imitates the propagation vector to be *k*

*DOI: http://dx.doi.org/10.5772/intechopen.89706*

!

turbulence.

technique.

magnetic field *B*

*<sup>n</sup>*<sup>0</sup> <sup>≈</sup> <sup>6</sup> � <sup>10</sup><sup>12</sup> cm�<sup>3</sup>

*∂*2 *By <sup>∂</sup>*t2 <sup>þ</sup> <sup>2</sup>*γ<sup>n</sup>*

**141**

*ni = nn*/2 <sup>≈</sup> 2.5 � 1012 cm�<sup>3</sup>

*ρ<sup>s</sup>* ≈ 0.49 cm, *k0x* ≈ 0.85 cm�<sup>1</sup>

*∂By* <sup>∂</sup><sup>t</sup> <sup>¼</sup> *<sup>λ</sup>*<sup>2</sup> *e*

respectively, be written as follows:

In general, both space as well as laboratory plasmas incorporate various kinds of inhomogeneities, few of these being magnetic field [13], temperature, and density perturbations. The coronal heating has been theoretically studied by Davila [14] by assuming Alfven absorption when Alfven velocity is dependent on background plasma density in transverse direction to ambient magnetic field. Also, the magnetic field fluctuations due to shear Alfvén waves have been observed experimentally [15]. The Alfven wave energy changes along radial direction have also been studied by Vincena et al. [16]. The present work is inspired from the experimental studies of formation of varying magnetic fields and electric fields by propagation of KAW in inhomogeneous plasmas [17].

A comparative study of experimental results with theoretical ones [18] fully support the launching of KAW in inhomogeneous plasma. There may be qualitative and quantitative deviation of ideal behavior [19] of KAW from that when inhomogeneities or nonlinearities are introduced due to kinetic theory. Therefore, it is always better to have understanding of magnetic fluctuations with an introduction of Landau damping, which is considered to be an important phenomenon as far as KAW propagation is concerned [20, 21]. Based on Landau fluid model, several theories [5, 22–24] have been proposed regarding the propagation of magnetohydrodynamic waves in collisionless plasma and Alfven filamentation process due to density channels. Taking into consideration the experimental observations by Houshmandyar and Scime [17], Sharma et al. [25] have studied the effect of Landau damping on KAW propagation in inhomogeneous magnetized plasma and have concluded that there is considerable damping effect on wave propagation due to phenomenologically incorporated Landau damping factor.

The damping refers to gradual decrease in wave intensity and wave amplitude. Not only collisionless damping (like Landau damping) but collisional damping also plays an important role in wave energy decay phenomenon and hence marks its role in Tokamak/ITER fusion plasma processing. The electron-neutral and ion-neutral damping is dominant in scrape-off region of these fusion plasma devices where energy transportation of wave takes place, and hence, neutral particles are the one among the important factors playing a vital role [26] in drawing the turbulence level. Also, the wave dispersion characteristics in space plasmas may get affected by the inclusion of neutral particles in the wave dynamics. Some of these characteristics include particle heating, spicules formation [27], solar wind turbulence and acceleration [28], and chromospheric energy balance [29]. As a consequence of the charged particles of plasma colliding with the neutral ones, the plasma waves may get absorbed by interstellar clouds, photosphere, and chromosphere [30–33]. As far as effect of neutral-dominated collisions on Alfvén waves is concerned, several research groups have worked in this direction and have given their fruitful theories in their respective domains. Based on the theory of chromospheric heating by Alfvén waves, the ion-neutral collisional damping has been studied analytically by De Pontieu et al. [29] in partially ionized chromosphere. Alfvénic turbulence has been considered as the base by Krishan and Gangadhara [34] to get intercoupling mechanism among three plasma species which are neutrals, ions, and electrons. The mass loading of ions by neutral particles has been used by Houshmandyar and

*Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on… DOI: http://dx.doi.org/10.5772/intechopen.89706*

Scime [17] to experimentally study the collisional impact due to collisions between ion and neutral particles on the KAW propagation in plasma. Recently, Goyal and Sharma [35] have studied the effect of ion-neutral collisional damping on KAW turbulence.

The present research deals with a comparative study of the two damping phenomena discussed above viz. ion-neutral collisional damping and Landau (collisionless) damping on the propagation of KAW which have been studied separately by Sharma et al. [25] and Goyal and Sharma [35] in the light of inhomogeneous magnetized plasma. This comparison is very much important with respect to scrape-off regions of the fusion energy devices like Tokamak and ITER. The plasma inhomogeneity is taken into account by considering spatial (in x-z plane) inhomogeneities having scale lengths both in longitudinal and transverse directions to background magnetic field. Distinct model equations for KAW have been taken by considering damping factors corresponding to collisional and collisionless damping, and the results have been drawn, respectively, using two fluid simulation technique.

## **2. Equations governing the propagation of KAW**

The governing equations of KAW come into their present form by taking into account the propagation of the wave in two dimensions (x-z surface), which imitates the propagation vector to be *k* ! *<sup>o</sup>* ¼ *koxx*^þ *kozz*^where *k0x(k0z)* symbolizes the component of the wave vector in perpendicular (parallel) direction to ambient magnetic field *B* ! <sup>0</sup> along z direction. Also, the plasma having background density *n*<sup>0</sup> is considered to be inhomogeneous in x as well as in z directions and hence is supposed to follow the density profile *<sup>n</sup>*<sup>0</sup> exp �*x*<sup>2</sup>*=L*<sup>2</sup> *<sup>x</sup>* � *<sup>z</sup>*<sup>2</sup>*=L*<sup>2</sup> *z* , where *Lx* implies the inhomogeneity scale length perpendicular to ambient magnetic field *B* ! <sup>0</sup> while *Lz* implies the same in parallel direction. Initially, the wave field profile is considered to follow hollow Gaussian behavior [36]. As the current numerical study is inspired from the experimental technique [17] where the plasma environment was developed experimentally to study the wave propagation, we have used the same labora-

tory parameters in our model which are as follows: *B0* ≈ 560 G, plasma density *<sup>n</sup>*<sup>0</sup> <sup>≈</sup> <sup>6</sup> � <sup>10</sup><sup>12</sup> cm�<sup>3</sup> , neutral particle density *nn* <sup>≈</sup> <sup>5</sup> � <sup>10</sup><sup>12</sup> cm�<sup>3</sup> , ion density *ni = nn*/2 <sup>≈</sup> 2.5 � 1012 cm�<sup>3</sup> ,*T*<sup>e</sup> = 81,200 K,*Tn* ≈ *Ti* = 2900 K. Using these experimental values, the other characteristic plasma parameters may be obtained like: *<sup>υ</sup><sup>A</sup>* <sup>≈</sup> 2.5 � <sup>10</sup><sup>7</sup> cm/s, *cs* <sup>¼</sup> <sup>1</sup>*:*<sup>3</sup> � <sup>10</sup><sup>6</sup> cm/s, *<sup>β</sup>* <sup>¼</sup> *<sup>c</sup>*<sup>2</sup> *<sup>s</sup> =υ*<sup>2</sup> *A* <sup>≈</sup> 0.002, *<sup>υ</sup>te* <sup>≈</sup> 1.1 � <sup>10</sup><sup>8</sup> cm/s, *<sup>υ</sup>ti* <sup>≈</sup> 2.445 � <sup>10</sup><sup>5</sup> cm/s, <sup>ω</sup>ci = 2.68 � 106 rad/s, *<sup>λ</sup><sup>e</sup>* = 0.1 cm, *<sup>ρ</sup><sup>i</sup>* <sup>≈</sup> 0.09 cm, *ρ<sup>s</sup>* ≈ 0.49 cm, *k0x* ≈ 0.85 cm�<sup>1</sup> , and *k0z* ≈ 0.16 cm�<sup>1</sup> for ω/ωci = *f/fci* = 0.9.

Using the well-established methods [26–28], Maxwell's equations and Faraday's law, the two dynamical equations for finite amplitude KAW incorporating the effects of ion-neutral collisional damping [35] and Landau damping [25] can, respectively, be written as follows:

$$\begin{split} \frac{\partial^2 B\_{\mathcal{V}}}{\partial t^2} + 2\gamma\_n \frac{\partial B\_{\mathcal{V}}}{\partial \mathbf{t}} &= \lambda\_\epsilon^2 \frac{\partial^4 B\_{\mathcal{V}}}{\partial \mathbf{x}^2 \partial t^2} - \left(\nu\_\text{in}^2 \lambda\_\epsilon^2 + \nu\_\text{A}^2 \rho\_i^2\right) \frac{\partial^4 B\_{\mathcal{V}}}{\partial \mathbf{x}^2 \partial \mathbf{z}^2} + \nu\_\text{A}^2 \exp\left(\mathbf{x}^2/\mathbf{L}\_x^2 + \mathbf{z}^2/\mathbf{L}\_x^2\right) \frac{\partial^2 B\_{\mathcal{V}}}{\partial \mathbf{z}^2} \\ &+ \frac{\nu\_\text{A}^2}{\text{o}\_{\text{c}i}^2} \left(\mathbf{1} + \frac{2(\eta\_\text{n} \eta\_{\text{i}i}/\mathbf{n}\_\text{i} \nu\_{\text{n}i})}{\mathbf{1} + \text{o}^2/\nu\_{\text{n}i}^2}\right) \frac{\partial^4 B\_{\mathcal{V}}}{\partial \mathbf{z}^2 \partial \mathbf{t}^2} \end{split} \tag{1}$$

can be thought of as the probable mechanism to describe the magnetic fluctuations which are observed near earth's magnetopause [8]. The theoretical and experimental studies [9, 10] support the importance of KAW in solar wind turbulence energy cascade and particle heating. Earth's magnetopause is also known to have a vital role played by KAW in the mechanism of stochastic ion heating [11]. Nykyri et al. [12] observed the turbulent spectra in high altitude cusp and revealed that break in observed spectra might be due to damping of obliquely propagating KAW. In general, both space as well as laboratory plasmas incorporate various kinds of inhomogeneities, few of these being magnetic field [13], temperature, and density perturbations. The coronal heating has been theoretically studied by Davila [14] by assuming Alfven absorption when Alfven velocity is dependent on background plasma density in transverse direction to ambient magnetic field. Also, the magnetic field fluctuations due to shear Alfvén waves have been observed experimentally [15]. The Alfven wave energy changes along radial direction have also been studied by Vincena et al. [16]. The present work is inspired from the experimental studies of formation of varying magnetic fields and electric fields by propagation of KAW

A comparative study of experimental results with theoretical ones [18] fully support the launching of KAW in inhomogeneous plasma. There may be qualitative and quantitative deviation of ideal behavior [19] of KAW from that when inhomogeneities or nonlinearities are introduced due to kinetic theory. Therefore, it is always better to have understanding of magnetic fluctuations with an introduction of Landau damping, which is considered to be an important phenomenon as far as KAW propagation is concerned [20, 21]. Based on Landau fluid model, several theories [5, 22–24] have been proposed regarding the propagation of magnetohydrodynamic waves in collisionless plasma and Alfven filamentation process due to density channels. Taking into consideration the experimental observations by Houshmandyar and Scime [17], Sharma et al. [25] have studied the effect of Landau damping on KAW propagation in inhomogeneous magnetized plasma and have concluded that there is considerable damping effect on wave propagation due to

The damping refers to gradual decrease in wave intensity and wave amplitude. Not only collisionless damping (like Landau damping) but collisional damping also plays an important role in wave energy decay phenomenon and hence marks its role in Tokamak/ITER fusion plasma processing. The electron-neutral and ion-neutral damping is dominant in scrape-off region of these fusion plasma devices where energy transportation of wave takes place, and hence, neutral particles are the one among the important factors playing a vital role [26] in drawing the turbulence level. Also, the wave dispersion characteristics in space plasmas may get affected by the inclusion of neutral particles in the wave dynamics. Some of these characteristics include particle heating, spicules formation [27], solar wind turbulence and acceleration [28], and chromospheric energy balance [29]. As a consequence of the charged particles of plasma colliding with the neutral ones, the plasma waves may get absorbed by interstellar clouds, photosphere, and chromosphere [30–33]. As far as effect of neutral-dominated collisions on Alfvén waves is concerned, several research groups have worked in this direction and have given their fruitful theories in their respective domains. Based on the theory of chromospheric heating by Alfvén waves, the ion-neutral collisional damping has been studied analytically by De Pontieu et al. [29] in partially ionized chromosphere. Alfvénic turbulence has been considered as the base by Krishan and Gangadhara [34] to get intercoupling mechanism among three plasma species which are neutrals, ions, and electrons. The mass loading of ions by neutral particles has been used by Houshmandyar and

phenomenologically incorporated Landau damping factor.

in inhomogeneous plasmas [17].

*Progress in Fine Particle Plasmas*

**140**

and

$$\begin{split} \frac{\partial^2 B\_\text{y}}{\partial t^2} + 2\chi\_L \frac{\partial B\_\text{y}}{\partial t} &= \left( \frac{\partial^2}{\partial t^2} + 2\chi\_L \frac{\partial}{\partial t} \right) \lambda\_\epsilon^2 \frac{\partial^2 B\_\text{y}}{\partial x^2} - \upsilon\_\text{A}^2 \rho\_s^2 \frac{\partial^4 B\_\text{y}}{\partial x^2 \partial x^2} \\ &+ \frac{\upsilon\_\text{A}^2}{\text{o}\_\text{ci}^2} \left( \text{o}\_{\text{ci}}^2 + \frac{\partial^2}{\partial t^2} \right) \exp\left[ \frac{x^2}{L\_x^2} + \frac{x^2}{L\_x^2} \right] \frac{\partial^2 B\_\text{y}}{\partial x^2} \end{split} \tag{2}$$

where *γ<sup>n</sup>* is damping factor representing ion-neutral collisions and *γ<sup>L</sup>* is phenomenologically incorporated Landau damping factor obtained using the dispersion relation used by Lysak and Lotko [20] and Hasegawa and Chen [21]. The perturbations being very small in comparison to other terms in the above equations; their spatial derivatives may not be considered. Now, in order to study the damping mechanism, wave localization phenomenon is very much important which is done by using an envelope solution [36] into Eqs. (1) and (2) to obtain a set of normalised equations [25, 35] which look mathematically similar to nonlinear Schrödinger (NLS) equation. So, on solving these equations individually by following the standard algorithm which has been developed in order to solve NLS equation, one can infer about the wave localization, energy dissipation, and turbulent cascading of wave energy. Post examining the accuracy for NLS invariants up to 10�<sup>5</sup> , these normalized equations are accommodated keeping in view the modifications being done in order to solve NLS equation. By considering hollow Gaussian profile (initial simulation condition) before the start of simulation process (at *t* = 0), these equations are solved numerically. The hollow Gaussian profile at *t* = 0 is given by [35]:

$$B(\mathbf{x}, \mathbf{z}, \mathbf{0}) = a\_0(\mathbf{x}^2/r\_{10}^2 + \mathbf{z}^2/r\_{20}^2) \exp\left(-\mathbf{x}^2/r\_{10}^2 - \mathbf{z}^2/r\_{20}^2\right) \tag{3}$$

magnetic field intensity peaks due to density inhomogeneity in plasma. The wave gets focused and defocused due to these inhomogeneities. At a fixed time, the normalized magnetic field intensity has spatial spreads (in x-z plane), in almost similar manner in undamped as well as damped (ion-neutral collisions) cases but its peak values differ for different cases viz. undamped, ion-neutral collisions, and Landau damping. As expected, the field intensity has lower peak values for damped cases (ion-neutral collisions and Landau damping). A comparative picture of ionneutral collisional damping and Landau damping reveals that the Landau damping effect seems to be more dominating due to smaller magnetic field intensity peaks but it can also be seen that the spatial spread is also more in Landau damping case. We can make an inference that due to Landau damping effect, the spatial energy spread in smaller peaks is more pronounced than that in ion-neutral damping case where the intensity peaks resemble almost the ones in undamped case. Also, the damping phenomena get more pronounced as we scale up the time domain.

*Landau damping is most profound on KAW propagation in inhomogeneous magnetized plasma.*

*KAW having variations in intensity profile of magnetic field (normalized) with the direction of ambient magnetic field (z) and transverse direction (x) at t = 0, t = 9, and t = 13 with and without damping. At initial time, the presumed magnetic field profile, i.e., hollow Gaussian is followed. But as there is progress in time, the respective damping effects come into play which is depicted from the decrease in amplitude of the wave and hence the intensity. At times greater than t = 0, there is an enhancement in the peak intensity of respective curves for undamped, ion-neutral collision-dominated damping, and Landau damping cases, but on comparison, one can easily state that the Landau damping peaks are the lowest for a particular time. There is formation of chaotic structures in case of Landau damping as is clear from smaller intensity peaks. Hence, one can infer that*

*Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on…*

*DOI: http://dx.doi.org/10.5772/intechopen.89706*

**Figure 1.**

**143**

The damping effect may also be observed on the temporal variation of magnetic field in **Figure 2** for a specific spatial location (x, z). It can be clearly observed from **Figure 2** that at the initial times, all the fluctuations are of high amplitude (magnetic field) but with the passage of time, the damping effects come into play as fluctuations die down for both ion-neutral collisional as well as Landau damping effect. This behavior is indicated from **Figure 1** as well. Comparatively, Landau damping effect is more pronounced than the collisional effect. The Fourier transform of undamped as well as damped magnetic field fluctuations as given in **Figure 2** results into the magnetic power spectra of the respective cases, which are indicated as the j j *<sup>B</sup>*<sup>ω</sup> <sup>2</sup> varying with frequency <sup>ω</sup> in **Figure 3**. These spectra have been obtained using a magnetic field fluctuations' time window �10 μs and indicate a clear difference due to damping effects for some specified values of inhomogeneity scale lengths. Once again, it is depicted from **Figure 3** that the Landau damping effect is more pronounced than ion-neutral collisional effect on the propagation of

KAW which is indicative from the initial spectral structures.

Here, *a*<sup>0</sup> and *r*10, respectively, dictate the amplitude of KAW and the scale size of wave field in perpendicular direction to the ambient magnetic field *B0* before start of the simulation and *r*2 0 is the same in parallel direction. *r*10ð Þ *r*<sup>20</sup> may be termed as the system length in 2-D geometry. These system lengths while spreading in 2-D geometry expand even beyond the stature of hollow Gaussian profile in transverse as well as longitudinal directions to *B0*. Also, as *kx*(*kz*) is the normalized propagation constant in transverse (longitudinal) direction to *B0*, one can infer *<sup>r</sup>*<sup>10</sup> <sup>¼</sup> *<sup>k</sup>*�<sup>1</sup> *<sup>x</sup>* and *<sup>r</sup>*<sup>20</sup> <sup>¼</sup> *<sup>k</sup>*�<sup>1</sup> *<sup>z</sup>* . The regulation of KAW evolution is done by initially taking the system constants to be *a*<sup>0</sup> = 0.02 and *kx = kz* = 0.2 where *a0* is computed by the use of experimental parameters [17] and the wave amplitude is normalized by ambient magnetic field *B0*. Sharma and Singh [37] have established a distinct algorithm in order to numerically solve Eqs. (1) and (2). This algorithm is framed in such a way that it relies upon a pseudo spectral mechanism for spatial integrals having recurring lengths *lx = r*10*, lz = r*20, and a grid having squared dimensions (grid points) given by 26 � 26 . And, the temporal evolution is studied by the use of the predictor corrector mechanism besides finite difference method where the evolution advances with a step size Δ*t*≈ 10�4.

## **3. Outcomes and results**

The damping effects on the wave propagation are plainly detectable by having a view of **Figure 1** where the time-dependent variations in normalized magnetic field intensity for specific values of parallel and transverse inhomogeneity scales are shown. The figure depicts the formation of localized structures or normalized

*Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on… DOI: http://dx.doi.org/10.5772/intechopen.89706*

#### **Figure 1.**

and

*<sup>r</sup>*<sup>10</sup> <sup>¼</sup> *<sup>k</sup>*�<sup>1</sup>

**142**

*<sup>x</sup>* and *<sup>r</sup>*<sup>20</sup> <sup>¼</sup> *<sup>k</sup>*�<sup>1</sup>

(grid points) given by 26 � 26

**3. Outcomes and results**

evolution advances with a step size Δ*t*≈ 10�4.

*∂*2 *By <sup>∂</sup>t*<sup>2</sup> <sup>þ</sup> <sup>2</sup>*γ<sup>L</sup>*

*Progress in Fine Particle Plasmas*

*∂By <sup>∂</sup><sup>t</sup>* <sup>¼</sup> *<sup>∂</sup>*<sup>2</sup>

*B x*ð Þ¼ , *<sup>z</sup>*, 0 *<sup>a</sup>*<sup>0</sup> *<sup>x</sup>*<sup>2</sup>

*=r* 2 <sup>10</sup> <sup>þ</sup> *<sup>z</sup>*<sup>2</sup> *=r* 2 20 � � exp �*x*<sup>2</sup>

*<sup>∂</sup>t*<sup>2</sup> <sup>þ</sup> <sup>2</sup>*γ<sup>L</sup>*

þ *υ*2 A ω2 ci

� �

ω2 ci þ *∂*2 *∂t*2 � �

where *γ<sup>n</sup>* is damping factor representing ion-neutral collisions and *γ<sup>L</sup>* is phenomenologically incorporated Landau damping factor obtained using the dispersion relation used by Lysak and Lotko [20] and Hasegawa and Chen [21]. The perturbations being very small in comparison to other terms in the above equations; their spatial derivatives may not be considered. Now, in order to study the damping mechanism, wave localization phenomenon is very much important which is done by using an envelope solution [36] into Eqs. (1) and (2) to obtain a set of normalised equations [25, 35] which look mathematically similar to nonlinear Schrödinger (NLS) equation. So, on solving these equations individually by following the standard algorithm which has been developed in order to solve NLS equation, one can infer about the wave localization, energy dissipation, and turbulent cascading of wave energy. Post examining the accuracy for NLS invariants up to 10�<sup>5</sup>

normalized equations are accommodated keeping in view the modifications being done in order to solve NLS equation. By considering hollow Gaussian profile (initial simulation condition) before the start of simulation process (at *t* = 0), these equations are solved numerically. The hollow Gaussian profile at *t* = 0 is given by [35]:

Here, *a*<sup>0</sup> and *r*10, respectively, dictate the amplitude of KAW and the scale size of wave field in perpendicular direction to the ambient magnetic field *B0* before start of the simulation and *r*2 0 is the same in parallel direction. *r*10ð Þ *r*<sup>20</sup> may be termed as the system length in 2-D geometry. These system lengths while spreading in 2-D geometry expand even beyond the stature of hollow Gaussian profile in transverse as well as longitudinal directions to *B0*. Also, as *kx*(*kz*) is the normalized propagation constant in transverse (longitudinal) direction to *B0*, one can infer

the system constants to be *a*<sup>0</sup> = 0.02 and *kx = kz* = 0.2 where *a0* is computed by the use of experimental parameters [17] and the wave amplitude is normalized by ambient magnetic field *B0*. Sharma and Singh [37] have established a distinct algorithm in order to numerically solve Eqs. (1) and (2). This algorithm is framed in such a way that it relies upon a pseudo spectral mechanism for spatial integrals having recurring lengths *lx = r*10*, lz = r*20, and a grid having squared dimensions

the predictor corrector mechanism besides finite difference method where the

The damping effects on the wave propagation are plainly detectable by having a view of **Figure 1** where the time-dependent variations in normalized magnetic field intensity for specific values of parallel and transverse inhomogeneity scales are shown. The figure depicts the formation of localized structures or normalized

*∂ ∂t* *λ*2 *e ∂*2 *By <sup>∂</sup>x*<sup>2</sup> � *<sup>υ</sup>*<sup>2</sup>

exp

A*ρ*2 *s*

" # *∂*<sup>2</sup>

*x*2 *L*2 *x* þ *z*2 *L*2 *z*

*=r* 2 <sup>10</sup> � *<sup>z</sup>*<sup>2</sup> *=r* 2 20 � � (3)

*<sup>z</sup>* . The regulation of KAW evolution is done by initially taking

. And, the temporal evolution is studied by the use of

*∂*<sup>4</sup>*By ∂z*<sup>2</sup>*∂x*<sup>2</sup>

> *By ∂z*<sup>2</sup>

(2)

, these

*KAW having variations in intensity profile of magnetic field (normalized) with the direction of ambient magnetic field (z) and transverse direction (x) at t = 0, t = 9, and t = 13 with and without damping. At initial time, the presumed magnetic field profile, i.e., hollow Gaussian is followed. But as there is progress in time, the respective damping effects come into play which is depicted from the decrease in amplitude of the wave and hence the intensity. At times greater than t = 0, there is an enhancement in the peak intensity of respective curves for undamped, ion-neutral collision-dominated damping, and Landau damping cases, but on comparison, one can easily state that the Landau damping peaks are the lowest for a particular time. There is formation of chaotic structures in case of Landau damping as is clear from smaller intensity peaks. Hence, one can infer that Landau damping is most profound on KAW propagation in inhomogeneous magnetized plasma.*

magnetic field intensity peaks due to density inhomogeneity in plasma. The wave gets focused and defocused due to these inhomogeneities. At a fixed time, the normalized magnetic field intensity has spatial spreads (in x-z plane), in almost similar manner in undamped as well as damped (ion-neutral collisions) cases but its peak values differ for different cases viz. undamped, ion-neutral collisions, and Landau damping. As expected, the field intensity has lower peak values for damped cases (ion-neutral collisions and Landau damping). A comparative picture of ionneutral collisional damping and Landau damping reveals that the Landau damping effect seems to be more dominating due to smaller magnetic field intensity peaks but it can also be seen that the spatial spread is also more in Landau damping case. We can make an inference that due to Landau damping effect, the spatial energy spread in smaller peaks is more pronounced than that in ion-neutral damping case where the intensity peaks resemble almost the ones in undamped case. Also, the damping phenomena get more pronounced as we scale up the time domain.

The damping effect may also be observed on the temporal variation of magnetic field in **Figure 2** for a specific spatial location (x, z). It can be clearly observed from **Figure 2** that at the initial times, all the fluctuations are of high amplitude (magnetic field) but with the passage of time, the damping effects come into play as fluctuations die down for both ion-neutral collisional as well as Landau damping effect. This behavior is indicated from **Figure 1** as well. Comparatively, Landau damping effect is more pronounced than the collisional effect. The Fourier transform of undamped as well as damped magnetic field fluctuations as given in **Figure 2** results into the magnetic power spectra of the respective cases, which are indicated as the j j *<sup>B</sup>*<sup>ω</sup> <sup>2</sup> varying with frequency <sup>ω</sup> in **Figure 3**. These spectra have been obtained using a magnetic field fluctuations' time window �10 μs and indicate a clear difference due to damping effects for some specified values of inhomogeneity scale lengths. Once again, it is depicted from **Figure 3** that the Landau damping effect is more pronounced than ion-neutral collisional effect on the propagation of KAW which is indicative from the initial spectral structures.

size of the dust particles is much larger than the neutral atoms, their collisional cross section would be much larger and hence the collisional damping would have been enhanced. Also, the frequency of collisions would be increased by incorporating dust particles into the plasma, their negative charges may also put adverse effect on the wave propagation and hence, damping phenomenon. In this way, the dust particles' role is unavoidable when we talk about collisional effects on plasma, but in the present study, we have concentrated only on ion-neutral collisions. The effect of

*Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on…*

A comparative study about relative profoundness of damping effects has been done when KAW propagates in two-dimensional inhomogeneous magnetized plasma. The study is based on the numerical simulation technique that has been employed to solve the dynamical equations using two-fluid approach. These equations serve as the numerical representations of KAW propagation when ion-neutral collisional damping and Landau damping effects are separately incorporated into the KAW dynamics. In collisionless plasma, the prominent dissipation phenomenon, i.e., Landau damping is due to electrons unlike neutral particle-dominated plasma where the dissipation is primarily by ion-neutral collisions. Likewise, the neutral atoms colliding with other plasma species may cause dissipation of interplanetary magnetic field. Also, in scrape-off region of Tokamak/ITER, the major cause of particle energization is due to ion-neutral collisions because of presence of neutral atoms in plasma. But it is longitudinal component of electric field that causes particle heating in collisionless plasma because KAW follows this field due to its large transverse wavenumber and hence small wavelength. Due to dependency of Landau damping on wavenumber, it affects the localization amplitude, and hence, the localization of KAW happens due to wave packets which have varying wave numbers. The formation of these wave packets is attributed to the inclusion of inhomogeneities in KAW dynamics. It is clear from the figures that there is appreciable effect of damping phenomena on magnetic field fluctuations as well as on intensity profiles and the energy spread over a wide frequency range can be interpreted from the turbulence analysis. The collisional dynamics which considers ions and neutral particles to be at same temperatures has clear effect on wave dissipation but it is the Landau damping which shows relatively weighty effect as indicated by the wave amplitudes and the resultant frequency spectra. Furthermore, there is no nonlinearity included in KAW dynamics. Some interesting features of the dynamics are expected to be observed when temperature difference between ion and neutral particle dynamics are considered and/or when the

This work is partially supported by DST (India) and ISRO (India) under

dust particles will be studied in our future work.

*DOI: http://dx.doi.org/10.5772/intechopen.89706*

**4. Conclusion**

nonlinear dynamics are included.

**Acknowledgements**

RESPOND program.

**145**

#### **Figure 2.**

*The time scales of wave magnetic field fluctuations showing the effects of ion-neutral collisional damping and Landau damping at a particular location (x, z). As the time advances, the Landau damping is more dominant over ion-neutral collisional damping, which is indicated from the dying down of fluctuations in case of Landau damping.*

#### **Figure 3.**

*Taking the Fourier transform of the data obtained in Figure 2 depicts the power spectra of the respective data for undamped, ion-neutral collisional-dominated and Landau-dominated wave propagation. The power spectra are shown to be varying* ∣*Bw*∣ *<sup>2</sup> versus ω and can clearly indicate the dominance of Landau damping shown as least on energy scale.*

Also, it is clearly observed from **Figure 3** that on inclusion of damping effects, the spectral behavior (almost) retains its trend in the initial frequency regimes which indicates the linear propagation of KAW. But at the later stage (higher frequencies), the turbulent effects come into play and KAW propagation ceases to follow linear trend (undamped oscillations). Or this behavior might also be due to some unidentified nonlinearity or due to drawback of the numerical model used. The alteration in wave localization occurs due to the involvement of damping effects as is observed in **Figure 1**. The rates at which the spectral components are associated with spatially localized filamentary structures might emerge with time at enunciated rates. So, at a specific spatial location, the wave magnetic field is supposed to have fluctuation behavior with time, as is indicated in **Figure 2**.

In the present study, we deal only with the comparison of Landau damping and collisional damping on kinetic Alfvén waves in which plasma is considered to be inhomogeneous and magnetized. The purpose of the study is to look for the effect of these damping effects on wave propagation. In case of ion-neutral collisions, the inclusion of dust particles could have further affect the damping phenomena. As the *Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on… DOI: http://dx.doi.org/10.5772/intechopen.89706*

size of the dust particles is much larger than the neutral atoms, their collisional cross section would be much larger and hence the collisional damping would have been enhanced. Also, the frequency of collisions would be increased by incorporating dust particles into the plasma, their negative charges may also put adverse effect on the wave propagation and hence, damping phenomenon. In this way, the dust particles' role is unavoidable when we talk about collisional effects on plasma, but in the present study, we have concentrated only on ion-neutral collisions. The effect of dust particles will be studied in our future work.

## **4. Conclusion**

A comparative study about relative profoundness of damping effects has been done when KAW propagates in two-dimensional inhomogeneous magnetized plasma. The study is based on the numerical simulation technique that has been employed to solve the dynamical equations using two-fluid approach. These equations serve as the numerical representations of KAW propagation when ion-neutral collisional damping and Landau damping effects are separately incorporated into the KAW dynamics. In collisionless plasma, the prominent dissipation phenomenon, i.e., Landau damping is due to electrons unlike neutral particle-dominated plasma where the dissipation is primarily by ion-neutral collisions. Likewise, the neutral atoms colliding with other plasma species may cause dissipation of interplanetary magnetic field. Also, in scrape-off region of Tokamak/ITER, the major cause of particle energization is due to ion-neutral collisions because of presence of neutral atoms in plasma. But it is longitudinal component of electric field that causes particle heating in collisionless plasma because KAW follows this field due to its large transverse wavenumber and hence small wavelength. Due to dependency of Landau damping on wavenumber, it affects the localization amplitude, and hence, the localization of KAW happens due to wave packets which have varying wave numbers. The formation of these wave packets is attributed to the inclusion of inhomogeneities in KAW dynamics. It is clear from the figures that there is appreciable effect of damping phenomena on magnetic field fluctuations as well as on intensity profiles and the energy spread over a wide frequency range can be interpreted from the turbulence analysis. The collisional dynamics which considers ions and neutral particles to be at same temperatures has clear effect on wave dissipation but it is the Landau damping which shows relatively weighty effect as indicated by the wave amplitudes and the resultant frequency spectra. Furthermore, there is no nonlinearity included in KAW dynamics. Some interesting features of the dynamics are expected to be observed when temperature difference between ion and neutral particle dynamics are considered and/or when the nonlinear dynamics are included.

## **Acknowledgements**

This work is partially supported by DST (India) and ISRO (India) under RESPOND program.

Also, it is clearly observed from **Figure 3** that on inclusion of damping effects, the spectral behavior (almost) retains its trend in the initial frequency regimes which indicates the linear propagation of KAW. But at the later stage (higher frequencies), the turbulent effects come into play and KAW propagation ceases to follow linear trend (undamped oscillations). Or this behavior might also be due to some unidentified nonlinearity or due to drawback of the numerical model used. The alteration in wave localization occurs due to the involvement of damping effects as is observed in **Figure 1**. The rates at which the spectral components are associated with spatially localized filamentary structures might emerge with time at enunciated rates. So, at a specific spatial location, the wave magnetic field is sup-

*<sup>2</sup> versus ω and can clearly indicate the dominance of Landau damping*

*Taking the Fourier transform of the data obtained in Figure 2 depicts the power spectra of the respective data for undamped, ion-neutral collisional-dominated and Landau-dominated wave propagation. The power*

*The time scales of wave magnetic field fluctuations showing the effects of ion-neutral collisional damping and Landau damping at a particular location (x, z). As the time advances, the Landau damping is more dominant over ion-neutral collisional damping, which is indicated from the dying down of fluctuations in case of Landau*

In the present study, we deal only with the comparison of Landau damping and collisional damping on kinetic Alfvén waves in which plasma is considered to be inhomogeneous and magnetized. The purpose of the study is to look for the effect of these damping effects on wave propagation. In case of ion-neutral collisions, the inclusion of dust particles could have further affect the damping phenomena. As the

posed to have fluctuation behavior with time, as is indicated in **Figure 2**.

**Figure 3.**

**144**

**Figure 2.**

*Progress in Fine Particle Plasmas*

*damping.*

*spectra are shown to be varying* ∣*Bw*∣

*shown as least on energy scale.*

*Progress in Fine Particle Plasmas*

## **Author details**

Ravinder Goyal\* and R.P. Sharma Centre for Energy Studies, Indian Institute of Technology, Delhi, India

\*Address all correspondence to: ravig.iitd@gmail.com

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

**References**

2005;**95**:065002

1976;**37**:690

2007;**99**:175004

1975;**35**:370

2009;**16**:012304

2008;**35**:L17S08

**147**

L9

Science. 1984;**32**:1387

[1] Chaston CC, Phan TD, Bonnell JW, Mozer FS, Auna M, Goldstein ML, et al. Drift-Kinetic Alfvén waves observed near a reconnection X line in the Earth's magnetopause. Physical Review Letters.

*DOI: http://dx.doi.org/10.5772/intechopen.89706*

and dissipation of solar-wind turbulence at the electron gyroscale. Physical Review Letters. 2009;**102**:231102

[10] Howes GG, Dorland W, Cowley SC,

simulations of magnetized turbulence in astrophysical plasmas. Physical Review

[11] Johnson JR, Cheng CZ. Stochastic ion heating at the magnetopause due to kinetic Alfvén waves. Geophysical Research Letters. 2001;**28**:4421

[12] Nykyri K, Grison B, Cargill PJ, Lavraud B, Lucek E, Dandouras I, et al. Origin of the turbulent spectra in the high-altitude cusp: Cluster spacecraft observations. Annales de Geophysique.

[13] Chen HH, Liu CS. Solitons in nonuniform media. Physical Review

[14] Davila JM. Heating of the solar corona by the resonant absorption of Alfven waves. The Astrophysical

[16] Vincena S, Gekelman W, Maggs J. Shear Alfvén wave perpendicular propagation from the kinetic to the inertial regime. Physical Review Letters.

[17] Houshmandyar S, Scime E. Ducted kinetic Alfvén waves in plasma with steep density gradients. Physics of

[18] Scime EE, Keiter PA, Balkey MM, Kline JL, Sun X, Keesee AM, et al. The

Hammett GW, Quataert E, Schekochihin AA, et al. Kinetic

*Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on…*

Letters. 2008;**100**:065004

2006;**24**:1057

Letters. 1976;**37**:693

Journal. 1987;**317**:514

**102**:7225

2004;**93**:105003

Plasmas. 2011;**18**:112111

[15] Gekelman W, Vincena S, Leneman D, Maggs J. Laboratory experiments on shear Alfvén waves and their relationship to space plasmas. Journal of Geophysical Research. 1997;

[2] Hasegawa A, Mima K. Exact solitary Alfvén wave. Physical Review Letters.

[3] Chaston CC, Wilber M, Mozer FS, Fujimoto M, Goldstein ML, Acuna M, et al. Mode conversion and anomalous transport in Kelvin-Helmholtz vortices and kinetic Alfvén waves at the Earth's magnetopause. Physical Review Letters.

[4] Goertz CK. Kinetic Alfvén waves on auroral field lines. Planetary and Space

[5] Hasegawa A, Chen L. Kinetic process of plasma heating due to Alfvén wave excitation. Physical Review Letters.

[6] Salem CS, Howes GG, Sundkvist D, Bale SD, Chaston CC, Chen CHK, et al. Identification of kinetic Alfvén wave

Astrophysical Journal Letters. 2012;**745**:

[7] Carter TA, Maggs JE. Modifications of turbulence and turbulent transport associated with a bias-induced confinement transition in the large plasma device. Physics of Plasmas.

[8] Chaston C, Bonnell J, McFadden JP, Carlson CW, Cully C, Le CO, et al. Turbulent heating and cross-field transport near the magnetopause from THEMIS. Geophysical Research Letters.

[9] Sahraoui F, Goldstein ML, Robert P, Khotyaintsev YV. Evidence of a cascade

turbulence in the solar wind.

*Turbulence Generation in Inhomogeneous Magnetized Plasma Pertaining to Damping Effects on… DOI: http://dx.doi.org/10.5772/intechopen.89706*

## **References**

[1] Chaston CC, Phan TD, Bonnell JW, Mozer FS, Auna M, Goldstein ML, et al. Drift-Kinetic Alfvén waves observed near a reconnection X line in the Earth's magnetopause. Physical Review Letters. 2005;**95**:065002

[2] Hasegawa A, Mima K. Exact solitary Alfvén wave. Physical Review Letters. 1976;**37**:690

[3] Chaston CC, Wilber M, Mozer FS, Fujimoto M, Goldstein ML, Acuna M, et al. Mode conversion and anomalous transport in Kelvin-Helmholtz vortices and kinetic Alfvén waves at the Earth's magnetopause. Physical Review Letters. 2007;**99**:175004

[4] Goertz CK. Kinetic Alfvén waves on auroral field lines. Planetary and Space Science. 1984;**32**:1387

[5] Hasegawa A, Chen L. Kinetic process of plasma heating due to Alfvén wave excitation. Physical Review Letters. 1975;**35**:370

[6] Salem CS, Howes GG, Sundkvist D, Bale SD, Chaston CC, Chen CHK, et al. Identification of kinetic Alfvén wave turbulence in the solar wind. Astrophysical Journal Letters. 2012;**745**: L9

[7] Carter TA, Maggs JE. Modifications of turbulence and turbulent transport associated with a bias-induced confinement transition in the large plasma device. Physics of Plasmas. 2009;**16**:012304

[8] Chaston C, Bonnell J, McFadden JP, Carlson CW, Cully C, Le CO, et al. Turbulent heating and cross-field transport near the magnetopause from THEMIS. Geophysical Research Letters. 2008;**35**:L17S08

[9] Sahraoui F, Goldstein ML, Robert P, Khotyaintsev YV. Evidence of a cascade and dissipation of solar-wind turbulence at the electron gyroscale. Physical Review Letters. 2009;**102**:231102

[10] Howes GG, Dorland W, Cowley SC, Hammett GW, Quataert E, Schekochihin AA, et al. Kinetic simulations of magnetized turbulence in astrophysical plasmas. Physical Review Letters. 2008;**100**:065004

[11] Johnson JR, Cheng CZ. Stochastic ion heating at the magnetopause due to kinetic Alfvén waves. Geophysical Research Letters. 2001;**28**:4421

[12] Nykyri K, Grison B, Cargill PJ, Lavraud B, Lucek E, Dandouras I, et al. Origin of the turbulent spectra in the high-altitude cusp: Cluster spacecraft observations. Annales de Geophysique. 2006;**24**:1057

[13] Chen HH, Liu CS. Solitons in nonuniform media. Physical Review Letters. 1976;**37**:693

[14] Davila JM. Heating of the solar corona by the resonant absorption of Alfven waves. The Astrophysical Journal. 1987;**317**:514

[15] Gekelman W, Vincena S, Leneman D, Maggs J. Laboratory experiments on shear Alfvén waves and their relationship to space plasmas. Journal of Geophysical Research. 1997; **102**:7225

[16] Vincena S, Gekelman W, Maggs J. Shear Alfvén wave perpendicular propagation from the kinetic to the inertial regime. Physical Review Letters. 2004;**93**:105003

[17] Houshmandyar S, Scime E. Ducted kinetic Alfvén waves in plasma with steep density gradients. Physics of Plasmas. 2011;**18**:112111

[18] Scime EE, Keiter PA, Balkey MM, Kline JL, Sun X, Keesee AM, et al. The

**Author details**

**146**

Ravinder Goyal\* and R.P. Sharma

*Progress in Fine Particle Plasmas*

provided the original work is properly cited.

Centre for Energy Studies, Indian Institute of Technology, Delhi, India

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

\*Address all correspondence to: ravig.iitd@gmail.com

hot helicon experiment: HELIX. Journal of Plasma Physics. 2015;**28**:125

[19] Chen L, Zonca F. Gyrokinetic theory of parametric decays of kinetic Alfvén waves. Europhysics Letters. 2011;**96**:35001

[20] Lysak RL, Lotko W. On the kinetic dispersion relation for shear Alfvén waves. Journal of Geophysical Research. 1996;**101**:5085

[21] Hasegawa A, Chen L. Kinetic processes in plasma heating by resonant mode conversion of Alfvén wave. Physics of Fluids. 1976;**19**:1924

[22] Borgogno D, Hellinger P, Passot T, Sulem PL, Trávníček PM. Alfvén wave filamentation and dispersive phase mixing in a high-density channel: Landau fluid and hybrid simulations. Nonlinear Processes in Geophysics. 2009;**16**:275

[23] Passot T, Sulem PL. Long-Alfvénwave trains in collisionless plasmas. II. A Landau-fluid approach. Physics of Plasmas. 2003;**10**:3906

[24] Passot T, Sulem PL. Filamentation instability of long Alfvén waves in warm collisionless plasmas. Physics of Plasmas. 2003;**10**:3914

[25] Sharma RP, Goyal R, Gaur N, Scime EE. Linear kinetic Alfvén waves in inhomogeneous plasma: Effects of Landau damping. Europhysics Letters. 2016;**113**:25001

[26] Houshmandyar S, Scime E. Enhanced neutral depletion in a static helium helicon discharge. Plasma Sources Science and Technology. 2012; **21**:035008

[27] Haerendel G. Weakly damped Alfvén waves as drivers of solar chromospheric spicules. Nature (London). 1992;**360**:241

[28] Allen LA, Habbal SR, Qiu HY. Thermal coupling of protons and neutral hydrogen in the fast solar wind. Journal of Geophysical Research. 1998;**103**:6551

[29] De Pontieu B, Martens PCH, Hudson S. Electromagnetic field equations for a moving medium with Hall conductivity. The Astrophysical Journal. 2001;**558**:859

[30] Piddington JH. Electromagnetic field equations for a moving medium with Hall conductivity. MNRAS. 1954; **114**:638

[31] Piddington JH. The motion of ionized gas in combined magnetic, electric and mechanical fields of force. MNRAS. 1954;**114**:651

[32] Lehnert B. Plasma physics on cosmical and laboratory scale. Nuovo Cimento. 1959;**13**(Suppl):59

[33] Osterbrock DE. The heating of the solar chromosphere, plages, and corona by magnetohydrodynamic waves. The Astrophysical Journal. 1961;**134**:347

Section 5

Granular Particles

**149**

[34] Krishan V, Gangadhara RT. Meanfield dynamo in partially ionized plasmas – I. Monthly Notices of the Royal Astronomical Society. 2008;**385**: 849

[35] Goyal R, Sharma RP. Effect of ionneutral collisions on the evolution of kinetic Alfvén waves in plasmas. Plasma Physics and Controlled Fusion. 2018;**60**: 035002

[36] Goyal R, Sharma RP, Scime EE. Time dependent evolution of linear kinetic Alfvén waves in inhomogeneous plasma. Physics of Plasmas. 2015;**22**: 022101

[37] Sharma RP, Singh HD. Density cavities associated with inertial Alfvén waves in the auroral plasma. Journal of Geophysical Research. 2009;**114**: A03109
