World Journal of Mechanics
Vol.4 No.8(2014), Article
ID:49285,13
pages
DOI:10.4236/wjm.2014.48027
Rayleigh-Taylor Instability in Magnetized Plasma
G. A. Hoshoudy
Department of Applied Mathematics and Computer Science, Faculty of Science, South Valley University, Kena, Egypt
Email: g_hoshoudy@yahoo.com
Copyright © 2014 by author and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
Received 12 April 2014; revised 8 May 2014; accepted 5 June 2014
ABSTRACT
The Rayleigh-Taylor instability in stratified plasma has been investigated in the presence of combined effect of horizontal and vertical magnetic field. The linear growth rate has been derived for the case where plasma with exponential density distribution is confined between two rigid planes by solving the linear MHD equations into normal mode. Some special cases have been particularized to explain the roles the variables of the problem play; numerical solutions have been made and some stability diagrams are plotted and discussed. The results show that, the growth rate depends on the horizontal and vertical components of magnetic field and also depends on the parameter (is constant and is the density-scale length). The maximum instability happens at and to get more stability model we select such that it is different than. The vertical magnetic field component have a greater effect than the horizontal magnetic field component in the case of large wavelength, while in the case of short wavelength, the horizontal magnetic field components have greater effect than the vertical magnetic field component.
Keywords:Rayleigh-Taylor Instability, Plasma, Horizontal and Vertical Magnetic Field
1. Introduction
It is well-known that, plasma is a hot ionized gas consisting of approximately equal numbers of positively charged ions and negatively charged electrons. The characteristics of plasmas are significantly different from those of ordinary neutral gases so that plasmas are considered a distinct fourth state of matter. There are numerous everyday uses for plasmas, for example, fluorescent lights and neon signs work because of plasma. Circuit features on micro-processor chips in computers also contain plasma. Even the exhaust that is emitted during rocket launches is actually plasma, not gas. Plasma is also used to manufacture industrial diamonds and superconducting films. Plasma is also a key technology in the development of alternative energy sources. Nuclear fusion, which is plasma based, is one of the most promising candidates for the energy needs of the future when fossil fuels finally run out. One of the important models rises in hydrodynamic plasma called the RayleighTaylor instability (RTI) problem [1] [2] . RTI can occur when dense plasma is supported against the gravity. Studies on RTI show that this instability is troublesome because it obstructs the realization of ICF [3] [4] . In astrophysics, this instability is related to the stellar structure and evolution [5] [6] . Hence, it is important to understand the physical mechanisms that can affect such instability, especially its suppression as well as the details of this instability. Such knowledge will aid our understanding of the origin of white dwarfs and type-Ia supernovas.
Also, it is well-known that, the plasmas are strongly influenced by magnetic fields (where the plasma parameters may vary with the application of magnetic field), where the behavior of plasmas in the presence of a magnetic field is among the oldest problems in plasma physics. It has been central in plasma fusion research since the early experiments on plasma confinement by a magnetic field in the 1950s, and it has remained great interest in plasma fusion studies that use contemporary sophisticated devices [7] . It is also an important problem for many plasma discharges used in processing semiconductor materials where the application of a magnetic field results in enhancement of some desirable features of specific plasma sources [8] . The linear growth rate of RTI has been obtained by Goldston and Rutherford [9] under fixed boundary conditions. The hydromagnetic stability of a magnetized plasma of variable density is of considerable importance in several astrophysical situations, e.g. in theories of sunspot magnetic fields, heating of solar corona and the stability of stellar atmospheres in magnetic fields. Ariel [10] investigated the stability of an inviscid compressible fluid of variable density in the presence of a uniform vertical magnetic field and viscous plasma has been investigated by Bhatia [11] . The instability of stratified plasmas in the presence of horizontal magnetic field of compressible plasmas is studied by Bhimsen [12] . The RTI of a plasma layer in the presence of a horizontal magnetic field is investigated, taking into account the effects of Hall-currents and an arbitrarily large density gradient by Ariel [13] . The effects of horizontal magnetic field, Hall currents and viscosity have been studied on the RTI of an incompressible infinitely conducting stratified plasma by Ahsan and Bhatia [14] . The effect of horizontal magnetic field of Hall currents have been investigated on the RTI of a finitely conducting stratified partially ionized plasma by Aiyub and Bhatia [15] . The instability of stratified plasmas in the presence of horizontal magnetic field of incompressible plasmas with the effect of a transverse velocity shear was studied by Wu et al. [16] . The effect horizontal magnetic field on RTI of a plasma layer in the presence of quantum mechanism was studied by Jintao et al. [17] . The effects of vertical magnetic field in the presence of quantum mechanism on RTI inviscous and viscous plasma were studied by Hoshoudy [18] [19] . The effects of magnetic field gradient on the Rayleigh-Taylor instability (RTI) with continuous magnetic field and density profiles were analytically investigated by Yang et al. [20] . The analytically and numerically investigated stabilization of the linear growth of the RTI from density gradients, magnetic fields, and quantum effects, in an ideal incompressible magnetized plasma was studied by Wang et al. [21] . Because of the great scientific interest in magnetized plasma problem, it is attempted to discuss the RTI problem of a stratified plasma layer in the presence of an variable magnetic field, where this problem corresponds physically (in astrophysics) to the RTI of an equatorial section of a planetary magnetosphere or of a stellar atmosphere where the magnetic field are perpendicular or parallels to gravity. In all the above-mentioned studies the behaviour of growth rates is considered in the presence of a variable magnetic field in direction only or in direction only. Here, the effect of magnetic field in both and direction on RTI problem for a finite thickness layer of incompressible plasmas are studied. The dispersion relation is obtained analytically and numerically analyzed.
2. Linearized Equations
We consider the strata of incompressible and inviscous plasma as a fluid of electrons and immobile ions. The plasma is immersed in a magnetic field, where the relevant liner perturbation equations may be written as (Refs. [9] -[17] )
, (1)
(2)
(3)
(4)
Here, , and, are the perturbations in the velocity, pressure, magnetic field, and density, respectively. While, , and the fluid is arranged in horizontal strata, then is a function of the vertical coordinate only (i.e.) and . Then the system of Equations (1)-(4) become:
, (5)
, (6)
, (7)
, (8)
(9)
. (10)
If we assume that the perturbation in any physical quantity takes the form:
, (11)
where and are horizontal components of the wave-number vector such that and (may be complex) is the frequency of perturbations or the rate at which the system departs from equilibrium. Using the Expression (11) in the system of Equations (5) - (10), we have:
, (12)
, (13)
, (14)
(15)
(16)
(17)
Eliminating some variables from the system of Equations (12) - (17) we have:
(18)
where,
(19)
3. A Continuously Stratified Plasma Layer
In this section we consider the case of incompressible continuously stratified plasma layer of thickness units confined between two rigid boundaries, in which the density and magnetic field distribution are given, respectively, by, and where, and and (the density-scale length) are constants, then Equation (18) takes the form:
(20)
where and are Alfvén velocity. Now, if we choose in the form
and by substituting in Equation (26), we will have an equation in both and. Then coefficients both and, respectively, are given by:
(21)
(22)
Now, we define the dimensionless quantities:
(23)
Then Equations (21) and (22), respectively, take the form:
(24)
(25)
Now, we put and for (stable oscillations), then Equations (24) and (25) may be given by:
(26)
(27)
To discuss the role of parameter’s problem we consider the next special cases from Equations (26) and (27).
(i) In the case of.
From Equation (27) we get, and substituting in Equation (26) we find that the square normalized growth rate given by:
, (28)
This case is considered by Goldston and Rutherford (see Ref. [9] ), which represents an exponentially growing perturbation (instability case).
(ii) In the case of.
A second time, from Equation (27) we get, and substituting in Equation (26), then the square normalized growth rate given by:
. (29)
This case studied in Refs. ([16] [17] ). It is clarified that, the horizontal magnetic field has stabilizing effect on RTI problem. This influence is obvious from Equations (28) and (29), where,
. Also, one can see that the square normalized growth rate decreases as
increases and the system arrives to complete stability case at. Now, if we use the relation under the condition , then,. Then Equation (29) becomes:
. (30)
At, the critical point becomes.
This implies that the complete stability happens at.
Also, if the perturbation in the Equation (11) in the form. Then
and the square normalized growth rate becomes, and the critical point of stability becomes this results correspond with Ref. [17] .
(iii) In the case of. Again, from Equation (27) we get, and substituting in Equation (26), the square normalized growth rate given by:
. (31)
Now, comparing between Equations (28) and (31), someone can observe that, the stabilizing role for the vertical magnetic field on the considerable system, where and the system arrives to complete stability case at. This indicates that, in the presence of, the stability role stratifies under the condition.
(iv) For the general case, if we can eliminate the term between Equations (26) and (27). Then the square normalized growth rate given in the from
(32)
In the general case, if we wish to look into the effect of both horizontal and vertical magnetic field together on the instability of the considered system, Equation(32) is to be numerically solved, where is function in the dimensionless quantities of horizontal and vertical components of the magnetic field, the wave number and, where is constant and is the density-scale length. In these figures, , , and.
The role of as a function horizontal in the presence of both horizontal and vertical magnetic field components is plotted in Figure 1 and Figure 2, where the square normalized growth rate
(a)(b)
Figure 1. The square normalized growth rate against the square normalized wave number at (a) at (b) at.
is plotted against the square normalized wave number. Figure 1(a) shows the role of, for example. One can see that decrease with increasing of. While Figure 1(b) shows the role of or, where decrease with decreasing of. These implies that the maximum instability in the presence of both vertical and horizontal magnetic field happens at. The same phenomenon hold in Figure 2(a) and Figure 2(b) at. In this case, if we put in Equation (32), then the maximum square normalized growth rate is:
(33)
From Equation (33) it can see that, the maximum square normalized growth rate (in the presence of both horizontal and vertical magnetic field components) given as we add the second term in Equations (29) and (31) to Equation (28). Figure 3 shows the role of both horizontal and vertical components’ magnetic field on the considered system at, where the square normalized growth rate is plotted against the square normalized wave number. One can see that, the magnitudes of in the presence of these parameters or either them is less than their magnitudes in the Goldston and Rutherford case, which means that these factors have a stabilizing effect on the considered system. Also, it can seen that, no mode of maximum instability exists when as the square normalized growth rate usually increases by increase with the square normalized wave number values. While, in the presence of quantum term, there is a mode of maximum instability, where the square normalized growth rate increases with increases through the range (at the square normalized growth rate arrives to the maximum instability, and when the square normalized growth rate starts to decreases as increases and then goes to the complete stable at (is the critical value for stability, at this point the square normalized growth rate goes to zero). This means that the horizontal component of magnetic field has a crucial capability to suppress the instability that satisfies for the large wave number (short wavelength), while the vertical component of magnetic field has this strength just for the very short wave number (large wavelength). Also, the considerable model is more stability in the presence of both horizontal and vertical components of the magnetic field. This case shows in Figure 4, where the magnitudes of decrease with increasing of both horizontal and vertical components of magnetic field
Figure 3. The square normalized growth rate against the square normalized wave number at in the presence of, and.
Figure 4. The square normalized growth rate against the square normalized wave number at for different values of the magnetic field.
that are less than their counterpart at (Goldston and Rutherford results). Figure 5(a) and Figure 5(b) show the role of magnetic field when have
(a)(b)
Figure 5. The square normalized growth rate against the square normalized wave number for different values of the magnetic field. (a) at (b) at.
different values than (at, respectively). In Figure 5(a), one can see that, the maximum instability (maximum square of growth rate) at in Figure 5(a), respectively, are, they are less than their counterpart in Figure 4. Also, the critical point for stability in Figure 5(a), respectively, are, they are less than their counterpart in Figure 4. This implies that the magnitudes of in Figure 5(a) at are less than their counterpart in Figure 4 at. The same phenomenon holds in Figure 5(b) but at, where both the maximum instability and the critical point for stability in Figure 5(b) are less than their counterpart in Figure 4 at, which also indicates that, the magnitudes of at are less than their counterpart at.
4. Conclusions
Finally, we have investigated the effect of magnetic field in both horizontal and vertical direction on the Rayleigh-Taylor instability of stratified incompressible plasmas layer. We can summarize the results as follows:
(i) In the presence of horizontal magnetic field only or vertical magnetic field only the growth rate depends on
of plasma layer and takes one value. While, in the presence of both horizontal and vertical magnetic field, the square growth rate (Equation (32)) depends on the dimensionless thickness
and also depends on the parameter.
(ii) The constant plays an important stabilizing role in the growth rate, where the maximum instability happens at and to get more stability model, the values must be different than (i.e. or). In other words, in studying the combined effect of horizontal and vertical magnetic field on RTI of stratified plasma, the system will be more stable when we accept that is different than −0.5 (i.e.).
(iii) Our model is more stable than those considered in previous study (horizontal magnetic field only or vertical magnetic field only). This divergence ascribes to the stabilizing role that the magnetic field plays on RTI problem in the presence of both horizontal and vertical components, where the presence of both these parameters dissipates the energy of any disturbance and thereby the system becomes more stability.
References
- Rayleigh, L. (1882) Investigation of the Character of the Equilibrium of an Incompressible Heavy Fluid of Variable Density. Proceedings of the London Mathematical Society, 14, 170-177.
- Taylor, G.I. (1950) The Instability of Liquid Surfaces When Accelerated in a Direction Perpendicular to Their Planes. Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 201, 192-196.
- Sanz, J. (1994) Self-Consistent Analytical Model of the Rayleigh-Taylor Instability in Inertial Confinement Fusion. Physical Review Letters, 73, 2700-2703. http://dx.doi.org/10.1103/PhysRevLett.73.2700
- Lindl, J.D. (1995) Development of the Indirect-Drive Approach to Inertial Confinement Fusion and the Target Physics Basis for Ignition and Gain. Physics of Plasmas, 2, 3933. http://dx.doi.org/10.1063/1.871025
- Cabot, W.H. and Cook, A.W. (2006) Reynolds Number Effects on Rayleigh-Taylor Instability with Possible Implications for Type-Ia Supernovae. Nature Physics, 2, 562-568.
- Timmes, F.X. and Woosley, S.E. (1992) The Conductive Propagation of Nuclear Flames. I - Degenerate C + O and O + NE + MG White Dwarfs. Astrophysical Journal, 396, 649-667. http://dx.doi.org/10.1086/171746
- Blinnikov, S. and Sorokina, E. (2004) Type Ia Supernova Models: Latest Developments. Astrophysics and Space Science, 290, 13-28. http://dx.doi.org/10.1023/B:ASTR.0000022161.03559.42
- Chen, F.F. (1974) Introduction to Plasma Physics. Plenum, New York.
- Goldston, R.J. and Rutherford, P.H. (1997) Introduction to Plasma Physics. Institute of Physics, London.
- Ariel, P.D. (1971) Rayleigh-Taylor Instability of Compressible Fluids in the Presence of a Vertical Magnetic Field. Applied Scientific Research, 24, 294-304.
- Bhatia, P.K. (1974) Rayleigh-Taylor Instability of a Viscous Compressible Plasma of Variable Density. Astrophysics and Space Science, 26, 319-325. http://dx.doi.org/10.1007/BF00645614
- Shivamoggi, B.k. (1981) Rayleigh-Taylor Instability of Compressible Plasma in a Vertical Magnetic Field. Astrophysics and Space Science, 79, 3-9. http://dx.doi.org/10.1007/BF00655900
- Ariel, P.D. (1991) Rayleigh-Taylor Instability of a Hall Plasma with Arbitrary Density Gradient. Astrophysics and Space Science, 184, 205-219.
- Ali, A. and Bhatia P.K. (1993) Rayleigh-Taylor Instability of a Stratified Hall Plasma in Two-Dimensional Horizontal Magnetic Field. Physica Scripta, 47, 567-570. http://dx.doi.org/10.1088/0031-8949/47/4/016
- Khan, A. and Bhatia, P.K. (1993) Rayleigh-Taylor Instability of a Finitely Conducting Partially Ionized Hall Plasma. Physica Scripta, 48, 607-611. http://dx.doi.org/10.1088/0031-8949/48/5/017
- Wu, Z., Zhang, W., Li, D. and Yang. W. (2004) Effect of Magnetic Field and Equilibrium Flow on Rayleigh-Taylor Instability. Chinese Physics Letters, 21, 2001-2004.
- Cao, J.T., Ren, H.J., Wu, Z.W. and Chu, P.K. (2008) Quantum Effects on Rayleigh-Taylor Instability in Magnetized Plasma. Physics of Plasmas, 15, 012110. http://dx.doi.org/10.1063/1.2833588
- Hoshoudy, G.A. (2010) Quantum Effects on Rayleigh-Taylor Instability of Incompressible Plasma in a Vertical Magnetic Field. Chinese Physics Letters, 27, 125201. http://dx.doi.org/10.1088/0256-307X/27/12/125201
- Hoshoudy, G.H. (2011) Rayleigh-Taylor Instability in Quantum Magnetized Viscous Plasma. Plasma Physics Reports, 37, 775-784. http://dx.doi.org/10.1134/S1063780X11080046
- Yang, B.L., Wang, L.F., Ye, W.H. and Xue, C. (2011) Magnetic Field Gradient Effects on Rayleigh-Taylor Instability with Continuous Magnetic Field and Density Profiles. Physics of Plasmas, 18, 072111. http://dx.doi.org/10.1063/1.3609773
- Wang, L.F., Yang, B.L., Ye, W.H. and He, X.T. (2012) Stabilization of the Rayleigh-Taylor Instability in Quantum Magnetized Plasmas. Physics of Plasmas, 19, 072704. http://dx.doi.org/10.1063/1.4737162