Rayleigh-Taylor Instability of Magnetized Plasma through Darcy Porous Medium ()
1. Introduction
The hydromagnetics stability of a magnetized plasma of varying density is of considerable importance in several astrophysical situations such as supernova explosions, in heating in solar corona, theories of sunspot magnetic fields, the formation and mixing of clouds and the stability of the stellar atmospheres in magnetic fields.
The classical study of the equilibrium of an incompressible, inviscid fluid of variable density was first undertaken by Rayleigh [1] , and later applied to all accelerated fluids by Taylor [2] . Rayleigh showed the equilibrium of a horizontal layer of incompressible, ideal fluid is stable or unstable according as the density increases or decreases anywhere in the vertically upward direction. Since then the problem is called Rayleigh-Taylor instability (RTI). RTI occurs naturally in many phenomenons of astrophysics, geophysics, and laboratories. It derives its character from the adverse density distribution of the matter, where the investigation of Rayleigh-Taylor instabilities in a magnetized plasma is a problem of considerable interest in space (ionospheric spread-F), fusion (curvature induced instabilities like interchange, ballooning, etc.) and the astrophysical plasmas.
Under various physical effects, the Rayleigh-Taylor instability problem of a finite layer of a fluid has been studied by several authors in hydrodynamics and in magnetohydrodynamics domain; the stabilizing effect of magnetic field on RTI problem for an incompressible plasma has been demonstrated by Kruskal and Schwarzschild [3] for a horizontal orientation of the magnetic field and by Hide [4] for a vertical orientation. The effects of viscosity and compressibility on the RTI of stratified plasma in the presence of magnetic field have been studied by Bhatia [5] . The RTI in a rotating plasma of variable density including simultaneously the effects of viscosity and the finiteness of the ion Larmor radius have been investigated by Bhatia and Steiner [6] . The effect of Hall current and finite electrical resistivity on the RTI of viscous, incompressible, finitely conducting plasma in a downward gravitational field under the influence of a uniform magnetic field normal to gravity has been studied by Kamla and Srivastav [7] . The RTI of an infinitely conducting stratified dusty plasma medium including the effects of FLR corrections in the presence of a horizontal magnetic field has been studied by Kamal and Chhajlani [8] . 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 Donald [9] . The effects of Hall currents and viscosity on the RTI of an incompressible infinitely conducting stratified plasma permeated by a two-dimensional horizontal magnetic field have been investigated by Ahsan and Bhatia [10] . The effects of Hall currents on the RTI of a finitely conducting stratified partially ionized plasma, where the plasma is permeated by a two dimensional horizontal magnetic field have been studied by Aiyub and Bhatia [11] . In the presence of magnitude of the gravitational acceleration, the RTI of stratified incompressible plasma has been studied by Goldston and Rutherford [12] . The RTI in the presence of horizontal magnetic field of incompressible plasma has been studied by Wu et al. [13] .
The RTI of magnetized plasma through porous medium problem has a great scientific interest, where this problem corresponds physically (in astrophysics) to the Rayleigh-Taylor instability of an equatorial section of a planetary magnetosphere or of a stellar atmosphere where the magnetic field is perpendicular or parallels to gravity. So, the RTI of a stratified plasma through porous medium in the presence or absence of magnetic field has been studied by a number of researchers (Chhajlani and Vaghela, Vyas and Chhajlani, Sharma and Bhardwaj, Sharma and Sharma, Sharma and Trilok, Sharma and Sunil, Shikha and Bhatia, Opara, Sharma and Sunil, Sunil and Sharma, Sharma and Thakur and Sharma and Rajput). In this case (Darcy’s model), the usual viscous term in the equation of motion is replaced by the resistive term
, where
is the fluid viscosity,
is the medium permeability and
is the Darcian (filter) velocity of the fluid.
In all the above-mentioned studies, the behaviour of growth rates is considered with respect to the porosity of porous medium and the medium permeability in the presence of an variable magnetic field in
direction only or in
direction only. Here, we will discuss the role of resistive term (Darcy’s term) besides the components of magnetic field in both
and
direction on growth rates of RTI of plasma layer
2. Formulation of the Problem
We consider the strata of incompressible and inviscous plasma as a fluid of electrons and immobile ions through Darcy porous medium in the presence of magnetic field
, where the relevant equations may be written, respectively (see references [12] -[25] ),
, (1)
, (2)
. (3)
For incompressible flow the fluid elements move without changing density is to say that the Lagrangian total derivative of density is zero, that is (see reference [1] )
(4)
where
is the velocity of the fluid,
is the density,
thermal pressure,
magnetic permeability,
coefficient of dynamic viscosity and
is the gravitational acceleration.
One can see that the set of Equations (1)-(4) is complete for describing the magnetic field effects on the R-T instability of incompressible plasma, since its number of equations exactly equals its number of unknown quantities: Two unknown vector quantities
and
and two unknown scalar quantities
and
. For the equilibrium profiles can be expressed in the form
and
.
Now, we assume a small perturbation in the system of Equations (1)-(4), where the perturbations in the velocity
, pressure
, magnetic field
, and density
, respectively, are
,
, and
. Then, the linearized equations can be easily derived from Equations (1)-(4) in the form
(5)
(6)
(7)
(8)
Now, let
,
,
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 (4)-(8) become
(9)
(10)
(11)
(12)
, (13)
(14)
If we assume that the perturbation in any physical quantity takes the form
, (15)
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 thee initial state. Using the expression (15) in the system of Equations (9)-(14), we have
, (16)
, (17)
, (18)
, (19)
, (20)
, (21)
Now, if we eliminate some of the variables from the system of Equations (16)-(21), we have a differential equation in ![](https://www.scirp.org/html/htmlimages\4-7501683x\42708472-191b-4085-897f-f5ff211ea2e2.png)
, (22)
![](https://www.scirp.org/html/htmlimages\4-7501683x\9f3ac9e9-36d8-4d9e-b7f7-e6b0e0a9f340.png)
(23)
3. A Continuously Stratified Plasma
In this section we consider the case of incompressible continuously stratified plasma layer of thickness
confined between two rigid boundaries, in which the density and magnetic field distribution are given, respectively, by
,
and ![](https://www.scirp.org/html/htmlimages\4-7501683x\82c3c9c1-48c1-4453-8de5-28c032104e0f.png)
where
,
and
and
(the density-scale length) are constants, then Equation (22) takes the form
(24)
where
and
are Alfvén velocity.
Now, if we choose
in the form
and by substituting in Equation (24), we will have an equation in both
and
. Then coefficients both
and
, respectively, are given by:
(25)
(26)
Now, we define the dimensionless quantities
(27)
Then Equations (25) and (26), respectively, take the form
(28)
(29)
Now, we put
and for
(stable oscillations), then Equations (28) and (29) may be given by:
(30)
(31)
Now, if we rearrange the above two equations (Equations (30) and (31)), we will have
(32)
(33)
From Equations (32) and (33) maybe we can specialize the next special cases:
(i) In the case of
and ![](https://www.scirp.org/html/htmlimages\4-7501683x\98eae95d-efa4-4a97-8482-06a502f540d6.png)
From Equation (33) we get
, and substituting in Equation (32) we find that the normalized growth rate given by
(34)
This case is considered by Goldston and Rutherford (see reference [12] ).
(ii) In the case of
and ![](https://www.scirp.org/html/htmlimages\4-7501683x\c9abcce1-5c58-467d-9e77-ad423d1e3886.png)
A second time, from Equation (33) we get
, and substituting in Equation (32), then the normalized growth rate given by
. (35)
This case studied in reference [12] . It is clarified that, the horizontal magnetic field has stabilizing effect on RTI problem. This influence is obvious from Equations (34) and (35), where
.
(iii) In the case of ![](https://www.scirp.org/html/htmlimages\4-7501683x\238d7761-0100-4dae-b8cb-8ae21183b383.png)
and ![](https://www.scirp.org/html/htmlimages\4-7501683x\fe78e258-2876-4406-9efe-3b62ef4849c1.png)
A third time, from Equation (33) we get
, and substituting in Equation (32), the normalized growth rate given by
. (36)
Now, comparing between Equations (34) and (36), someone can observe that, the stabilizing role for the vertical magnetic field on the considerable system, where
.
(iv) In the case of ![](https://www.scirp.org/html/htmlimages\4-7501683x\38e7a081-2d48-4bad-bf7d-730a2bdacef1.png)
and ![](https://www.scirp.org/html/htmlimages\4-7501683x\ecadb6a5-157d-479b-930d-cee029383ac5.png)
A fourth time, from Equation (33) we get
, and substituting in Equation (32), the dispersion rate given by
(37)
Then, the normalized growth rate becomes
. (38)
From Equations (34) and (38) it is very clear that,
.
The stabilizing effects of the horizontal, vertical magnetic and resistive term, unaccompanied, (above cases (i)-(iv)) on the RTI have been numerically presented in Figure 1.
(v) For the general case (
and
), if we eliminate the term
between Equations (32) and (33) the normalized dispersion relation takes in the form.
(39)
In this case the normalized growth rate given as.
(40)
where
(41)
(42)
(43)
In fact, the square of normalized growth rate
in Equation (40) is a function in the dimensionless quantities of horizontal
and vertical
components of the magnetic field, the dimensionless resistive term
, dimensionless Darcy term
, the wave number
and
(
, where
is constant and
is the density-scale length). The dimensionless quantities
,
and
are the parameters of problem that maybe take a different values. But the constant
is unknown in the general case (v), while in the absence of magnetic field or in the presence of either them (horizontal or vertical magnetic field components) we note that
(special cases (i)-(iv)). So, firstly we will discuss the role of constant
on the square of normalized growth rate
in the presence of horizontal, vertical magnetic field components and resistive term.
The role of constant
in the presence of both horizontal and vertical magnetic field components
and resistive term
is plotted in Figure 2, where the square normalized growth
rate
is plotted against
at
. For the values
that is less than
, one can see that, the magnitude of
decrease with decreasing of
. For the values ![](https://www.scirp.org/html/htmlimages\4-7501683x\ad192673-65ca-4487-96d6-baf67d229565.png)
that is greater than
the magnitude of
decrease while the magnitude of
increases. These implies that the maximum instability in the presence of both horizontal and vertical components of magnetic field happens at
(in the presence or absence of resistive term).
The general case (Equation (40)), that gives the effects of horizontal, vertical magnetic field and resistive term together on the instability of the considered system has been presented in Figure 3. where the square normalized growth rate
is plotted against the square normalized wave number
at
,
and for different values of
. Second time, one can see that, the maximum instability (maximum square normalized growth rate
) happens at
and at
the magnitudes of
are less than their counterpart at
. Also, the magnitudes of
at
are less than their counterpart at
. Moreover the magnitudes of
at
are less than their counterpart at
.
In the case
(maximum instability) the system of Equations (39)-(43) take the form:
(44)
Then the maximum normalized growth rate gives by
(45)
Finally, Figure 4 Shows the role of constant
, where the maximum happens at
. If we move
![](https://www.scirp.org/html/htmlimages\4-7501683x\6b35044c-8d8a-44e8-a85a-a7542b35cb35.png)
Figure 4. Sketch of the
through the range
.
toward the point
from the right or left hand the system moves toward instability and vice versa if we move from the point
to the left or right hand, the system tends to stability.
In closing this paper, the Rayleigh-Taylor instability in stratified plasma in the presence of combined effect of horizontal and vertical magnetic field through Darcy porous medium is considered. The solution of the system leads to a dispersion relation where the physical parameters are put in the dimensionless form. Some special cases are particularized to explain the roles that play the variables of the problem and numerical solutions are made. Some stability diagrams are plotted and discussed. The results show that, as the growth rate depends on the horizontal and vertical components of magnetic field and resistive term (Darcy’s term) also depends on the parameter
. Numerically the maximum instability (normalized growth rate) happens at
and then analytically the maximum instability gives in Equation (45). The system will be more stable if we select
such that to be different than −0.5. Finally, for the behavior of our selected system with respect to the parameter
we have not clear interpreting at this stage.