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

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] ),

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] )

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

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

If we assume that the perturbation in any physical quantity takes the form

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

Now, if we eliminate some of the variables from the system of Equations (16)-(21), we have a differential equation in

(23)

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

where, and and (the density-scale length) are constants, then Equation (22) takes the form

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:

Now, we define the dimensionless quantities

Then Equations (25) and (26), respectively, take the form

Now, we put and for (stable oscillations), then Equations (28) and (29) may be given by:

Now, if we rearrange the above two equations (Equations (30) and (31)), we will have

From Equations (32) and (33) maybe we can specialize the next special cases:

(i) In the case of and

From Equation (33) we get, and substituting in Equation (32) we find that the normalized growth rate given by

This case is considered by Goldston and Rutherford (see reference [12] ).

(ii) In the case of and

A second time, from Equation (33) we get, and substituting in Equation (32), then the normalized growth rate given by

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 and

A third time, from Equation (33) we get, and substituting in Equation (32), the normalized growth rate given by

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 and

A fourth time, from Equation (33) we get, and substituting in Equation (32), the dispersion rate given by

Then, the normalized growth rate becomes

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.

In this case the normalized growth rate given as.

where

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

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:

Then the maximum normalized growth rate gives by

Finally, Figure 4 Shows the role of constant, where the maximum happens at. If we move

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.