# Damping of linear nonadiabatic MHD waves in a flowing prominence medium.

1. IntroductionA large number of observations of oscillations in large scale magnetic structures such as prominence have been observed from space-borne and ground based instruments [1-5]. Small amplitude oscillations observed in quiescent solar prominence can be interpreted in terms of linear magnetohydyrodynamic waves. To study these small amplitude oscillations, one of the damping mechanism nonadiabaticities has been extensively studied because linear nonadiabatic magnetoacoustic waves have interesting physical effects such as time and spatial damping of the disturbances. Carbonell et al. [6] studied the time damping of linear nonadiabatic magnetohydrodynamic waves in an unbounded plasma with coronal properties and found that behavior of the damping changes significantly when going from very small to very large wavenumbers and becomes almost constant for very large wavenumbers. Spatial damping of linear nonadiabatic magnetoacoustic waves in prominence medium has been studied and it is found that the spatial damping is dominated by thermal conduction at shorter periods while the spatial damping is dominated by radiation at the periods greater than 1s [7]. Soler et al. [8] showed that nonadiabatic effects are efficient damping mechanism for magnetoacoustic modes.

Ground and space-borne satellites observations have confirmed that the solar and space plasmas are always dynamic showing flow [9-11]. Therefore, all the theoretical models should include the presence of flow. Nakariakov and Roberts [12] theoretically investigated the effect of steady mass flow in coronal and photospheric slabs. Carbonell et al. [13] discussed the combined effects of nonadiabatic mechanism and steady flow on the time damping of slow modes and found that flow does not modify the damping time but increases the period of slow mode dramatically for the flow velocities very close to nonadiabatic sound speed. In all previous studies, shear flow has not been taken into account. Therefore, in the present study, we investigate the local behavior of time damping of linear nonadiabtic magnetoacoustic waves in the presence of shear flow.

2. Basic Equations and Dispersion Relation

The equilibrium configuration considered in the present work is made of a homogeneous slab that is infinite in all directions. The magnetic field is uniform and is oriented along the z-direction. The background flow in z-direction is assumed to vary linearly in the v-direction as [v.sub.0] = [v.sub.0] (x) [??] = ax[??], with a constant a [14]. This assumption of the linear profile for velocity shear is applicable for wide variety of astrophysical shear flows. Parallel thermal conduction to the magnetic field, heating, and radiative losses are considered as nonadiabatic effects. The gravitational effects have been neglected because wavelengths are assumed much smaller than gravitational scale height. The basic MHD equations describing the plasma motion are given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

Here, [rho], p, T, v, and B denote the density, gas pressure, temperature, velocity, and magnetic field, respectively. The term [nabla] x ([kappa] x [nabla]T) represents the thermal conduction. The symbol L is heat-loss function and is given by

L([rho], T) = [chi][rho][T.sup.[alpha]], (2)

where [chi] and a are piecewise continuous functions depending on temperature [6].

Considering small perturbations from the equilibrium, the field quantities can be written as

B = [B.sub.0] + [B.sub.1] (r,t), v = [v.sub.0] (x) + [v.sub.1] (r, t), P = [P.sub.0] + [P.sub.1] (r,t), [rho] = [[rho].sub.0] + [[rho].sub.1] (r, t), T = [T.sub.0] + [T.sub.1] (r,t), (3)

where the subscripts "0" and "1" refer to equilibrium and perturbed quantities, respectively. Here, the equilibrium magnetic field B0 is homogeneous and is taken along z-axis; that is, [B.sub.0] = [B.sub.0][??]. The assumption about a homogeneous magnetic field may be reasonable for the large regions of quite sun but does not imply appropriate results. However, the assumption of a homogeneous magnetic field provides us a relatively simple mathematical model to study MHD waves in solar prominence. Several authors [6, 15, 16] used homogeneous magnetic field to study MHD waves in different regions of solar atmosphere. In the future, we will consider inhomogeneous magnetic field to study such MHD waves in solar prominence model and the results will be reported elsewhere.

We linearise the basic equations (1) to obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

where = [[kappa].sub.[parallel]] = [10.sup.-11] [T.sup.5/2] [Wm.sup.-1] [K.sup.-1] is thermal conduction parallel to magnetic field. If background flow is spatially inhomogeneous in the vertical direction, we can not apply Fourier transformation and a nonlocal approach would be more appropriate in this context. Though the local approximation has limited applicability, in the present study we are interested in investigating the behavior of nonadiabatic waves in the initial stage of the evolution of wave perturbations. Since local approximation is restricted to short wavelength fluctuations, (4) can be analyzed using Fourier transform over space and time for the study of local behavior of waves [17-20]. So we perform Fourier analysis in plane waves and assume that all perturbations are of the form [e.sup.i(k,r-[omega]t)] where k = [k.sub.x][??] + [k.sub.z][??] is the wavenumber and [omega] to is the frequency. From linear equations (4), we obtain the following dispersion relation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (6)

Here, [v.sub.A] is the Alfven velocity, [c.sub.s] is the sound velocity, and [v'.sub.0] is the derivative of [v.sub.0] with respect to x. If we simplify the dispersion relation (5), we can obtain a sixth order polynomial in [omega] due to flow. Since we are solving this dispersion relation numerically, there is no need to simplify it in [omega].

In the absence of flow, the dispersion relation (5) reduces to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (7)

Equation (7) yields the following fifth order dispersion relation in [omega]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (8)

which is similar to that obtained by Carbonell et al. [6] and it has been discussed by them in detail. So we will here solve and discuss the dispersion relation 5) to study the damping time and damping per period of fast and slow waves.

3. Results and Discussion

We solve the dispersion relation (5) numerically taking density, temperature, and magnetic field for the prominence region [21] corresponding to quiescent prominence conditions of solar atmosphere. The values of equilibrium parameters for the above mentioned prominence region are taken as [T.sub.0] - 8000 K, [[rho].sub.0] - 5 x [10.sup.-11] kg [m.sup.-3], [B.sub.0] - 5 G, [alpha] = 7.4, and [chi] - 1.76 x [10.sup.-13] (Hildner [21]). Here we consider the real wavenumber and find the complex solutions of the frequency as [[omega].sub.r] + i[[omega].sub.i],. We obtain two pairs of roots representing the fast and slow modes whereas purely imaginary root corresponds to thermal mode. The roots thus obtained do not represent instability because of small values of flow. For each pair ([[omega].sub.1], [[omega].sub.2]), [[omega].sub.1] = [[omega].sub.r] - i[[omega].sub.i] and [[omega].sub.2] = - [[omega].sub.r] - i[[omega].sub.i], [[omega].sub.1] corresponds to the parallel wave propagating towards the positive z-direction, whereas [[omega].sub.2] corresponds to antiparallel wave propagating towards negative z-direction. The damping time is [[tau].sub.d] - 1/[[omega].sub.i] and damping per period is [D.sub.p] - [[omega].sub.i]/ [[omega].sub.r]. We plot Figures 1-6 to study the effect of flow on damping time and damping per period of these waves with wavenumber. Figure 1 shows the variation of damping time of fast wave with wavenumber whereas Figure 2 shows the variation of damping time of slow wave with wavenumber for different values of flow. Damping time of both slow and fast waves decreases with the increase in the value of flow. For smaller values of wavenumber k damping time is higher than the larger values of wavenumber and becomes constant after a certain value of wavenumber. Damping time of slow wave is very much higher than the fast wave corresponding to the same set of values of flow as considered in fast wave. In the absence of flow, Figure 3 has been drawn to see the variation of damping time for slow and fast wave by changing X- and Y-scales from linear to log according to N. Kumar and P. Kumar [22] and it is found that results shown in Figure 3 are in agreement with results shown in Figure 1(b) of [22].

Figures 4 and 5 are drawn for damping per period versus wavenumber for fast and slow waves for the same set of values of flow, respectively. Damping per period of slow waves is very much higher than the damping per period of fast waves corresponding to similar values of flow. The patterns of increasing damping per period are almost similar for both slow and fast waves. Damping per period for both waves has linear dependence on wavenumber except for smaller values of k in case of slow waves. Figure 6 drawn in log scale depicts that the pattern of variation of damping per period is similar to that shown in Figure 1(c) by N. Kumar and P. Kumar [22] in the absence of flow.

4. Conclusions

Quiescent solar prominences and filaments are those features in solar atmosphere which show material flows and oscillations. So, in this paper, we have analyzed the time damping of MHD waves in nonadiabatic solar prominence plasma taking into account the effects of thermal conduction and shear flow. The nonadiabaticity is modeled by a general heat-loss function of thermodynamic parameters of plasma. We have obtained the dispersion relation using local approximation and solved it numerically to see the variations of damping time and damping per period of slow and fast MHD waves with wavenumber for different values of flow parameter. It is found that damping time of slow and fast waves decreases with the increase in flow and damping time of the slow wave is very much higher than that of the fast wave. It is observed that damping per period also decreases with the increase in flow.

http://dx.doi.org/10.1155/2014/541376

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

N. Kumar would like to acknowledge the support from IUCAA, Pune, India.

References

[1] C. E. DeForest and J. B. Gurman, "Observation of quasiperiodic compressive waves in solar polar plumes," Astrophysical Journal Letters, vol. 501, no. 2, pp. L217-L220,1998.

[2] L. Ofman, M. Romoli, G. Poletto, G. Noci, and J. L. Kohl, "Ultraviolet Coronagraph Spectrometer observations of density fluctuations in the solar wind," Astrophysical Journal Letters, vol. 491, no. 2, pp. Lm-L114,1997

[3] I. De Moortel, J. Ireland, A. W. Hood, and R. W. Walsh, "The detection of 3 & 5 min period oscillations in coronal loops," Astronomy and Astrophysics, vol. 387, no. 1, pp. L13-L16, 2002.

[4] I. De Moortel, J. Ireland, R. W. Walsh, and A. W. Hood, "Longitudinal intensity oscillations in coronal loops observed with trace: I. Overview of measured parameters," Solar Physics, vol. 209, no. 1, pp. 61-88, 2002.

[5] V. M. Nakariakov, L. Ofman, E. E. Deluca, B. Roberts, and J. M. Davila, "Trace observation of damped coronal loop oscillations: implications for coronal heating," Science, vol. 285, no. 5429, pp. 862-864, 1999.

[6] M. Carbonell, R. Oliver, and J. L. Ballester, "Time damping of linear non-adiabatic magnetohydrodynamic waves in an unbounded plasma with solar coronal properties," Astronomy and Astrophysics, vol. 415, no. 2, pp. 739-750, 2004.

[7] M. Carbonell, J. Terradas, R. Oliver, and J. L. Ballester, "Spatial damping of linear non-adiabatic magnetoacoustic waves in a prominence medium," Astronomy and Astrophysics, vol. 460, no. 2, pp. 573-581, 2006.

[8] R. Soler, R. Oliver, and J. L. Ballester, "The effect of the solar corona on the attenuation of small-amplitude prominence oscillations I. Longitudinal magnetic field," Astronomy and Astrophysics, vol. 471, no. 3, pp. 1023-1033, 2007

[9] A. H. Gabriel, F. Bely-Dubau, and P Lemaire, "The contribution of polar plumes to the fast solar wind," Astrophysical Journal Letters, vol. 589, no. 1, pp. 623-634, 2003.

[10] J. M. Krijger, T. Roudier, and M. Rieutord, "Photospheric flows measured with TRACE," Astronomy and Astrophysics, vol. 387, no. 2, pp. 672-677, 2002.

[11] R. Soler, R. Oliver, and J. L. Ballester, "Nonadiabatic magnetohydrodynamic waves in a cylindrical prominence thread with mass flow," Astrophysical Journal Letters, vol. 684, no. 1, pp. 725-735, 2008.

[12] V. M. Nakariakov and B. Roberts, "Magnetosonic waves in structured atmospheres with steady flows--I. Magnetic slabs," Solar Physics, vol. 159, no. 2, pp. 213-228,1995.

[13] M. Carbonell, R. Oliver, and J. L. Ballester, "Time damping of non-adiabatic MHD slow and thermal waves in a prominence medium: effect of a background flow," New Astronomy, vol. 14, no. 3, pp. 277-284, 2009.

[14] S. Poedts, A. D. Rogava, and S. M. Mahajan, "Shear-flow-induced wave couplings in the solar wind," Astrophysical Journal Letters, vol. 505, no. 1, pp. 369-375,1998.

[15] V. M. Nakariakov, C. A. Mendoza-Briceno, and M. H. Ibanez S., "Magnetoacoustic waves of small amplitude in optically thin quasi-isentropic plasmas," Astrophysical Journal Letters, vol. 528, no. 2, pp. 767-775, 2000.

[16] T. V. Zaqarashvili, R. Oliver, and J. L. Ballester, in Proceedings of the SOHO 13, Waves Oscillations and Small-Scale Transient Events in the Solar Atmosphere: A Joint View from SOHO and TRACE, no. ESASP-547, Palma de Mallorca, Spain, 2004.

[17] B. P Pandey and M. Wardle, "Hall instability of solar flux tubes in the presence of shear flow," Monthly Notices of the Royal Astronomical Society, vol. 426, pp. 1436-1443, 2012.

[18] B. P Pandey and M. Wardle, "Magnetic-diffusion-driven shear instability of solar flux tubes," Monthly Notices of the Royal Astronomical Society, vol. 431, pp. 570-581, 2013.

[19] H. Saleem, J. Vranjes, and S. Poedts, "Unstable drift mode driven by shear plasma flow in solar spicules," Astronomy and Astrophysics, vol. 471, no. 1, pp. 289-293, 2007

[20] G. D. Aburjania, G. Zimbardo, and O. A. Kharshiladze, "Effect of the shear flow in the generation and self-organization of internal gravity wave structures in the dissipative ionosphere," Plasma Physics Reports, vol. 38, no. 12, pp. 1055-1075, 2012.

[21] E. Hildner, "The formation of solar quiescent prominences by condensation," Solar Physics, vol. 35, no. 1, pp. 123-136,1974.

[22] N. Kumar and P Kumar, "Viscous damping of non-adiabatic MHD waves in an unbounded solar coronal plasma," Solar Physics, vol. 236, no. 1, pp. 137-153, 2006.

Nagendra Kumar, (1) Anil Kumar, (2) Himanshu Sikka, (2) and Pradeep Kumar (3)

(1) Department of Mathematics, M.M.H. College, Ghaziabad 201001, Uttar Pradesh, India

(2) Department of Mathematics, Maharaja Agrasen Institute of Technology, Rohini 110086, Delhi, India

(3) Department of Physics, K.G.K (P.G.) College, Moradabad 244001, Uttar Pradesh, India

Correspondence should be addressed to Anil Kumar; anilk_kang@yahoo.com

Received 23 December 2013; Revised 25 March 2014; Accepted 19 April 2014; Published 15 May 2014

Academic Editor: Wing Huen Ip

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Kumar, Nagendra; Kumar, Anil; Sikka, Himanshu; Kumar, Pradeep |

Publication: | Advances in Astronomy |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 2512 |

Previous Article: | Possibility of detection of exomoons with inclined orbits orbiting pulsar planets using the time-of-arrival analysis. |

Next Article: | Properties of the HII regions derived using integral field spectroscopy. |

Topics: |