The Effects of Rotation and Salt Concentration on Thermal Convection in a Linear Magneto-Fluid Layer Overlying a Porous Layer

A linear stability analysis is applied to a system consisting of a linear magneto-fluid layer overlying a porous layer affected by rotation and salt concentration on both layers. The flow in the fluid layer is governed by Navier-Stokes’s equations and while governed by Darcy-Brinkman’s law in the porous medium. Numerical solutions are obtained using Legendre polynomials. These solutions are studied through two modes of instability: stationary instability and overstability when the heat and the salt concentration are effected from above and below.


Motivations and Goals
Thermal instability theory has attracted considerable interest and has been recognized as a problem of fundamental importance in many fields of fluid dynamics.The earliest experiments to demonstrate the onset of thermal instability in fluids are those of Bernard's [1,2].Benard worked with very thin layers of an incompressible viscous fluid standing on a levelled metallic plate maintained at a constant temperature.The upper surface which was usually free and, being in contact with the air, was at a lower temperature.In his experiments, Benard deduced that a certain critical adverse temperature gradient must be exceeded before instability can set in.The instability of a layer of fluid heated from below and subjected to Coriolis forces has been studied by Chandrasekhar [3,4] for a stationary and overstability case.He showed that the presence of these forces usually has the effect of inhibiting the onset of thermal convection.Nield [5] considered the onset of salt-finger convection in a porous layer.Taunton et al. [6] considered the thermohaline instability and salt-finger in a porous medium and solved the boundary value problem.Sun [7] was the first to consider such a problem, and he used a shooting method to solve the linear stability equations.Sun [7] and Nield [8] used Darcy's law in formulating the equations of porous layer.In Darcy's law of motion in porous mediums, the Darcy resistance term took the place of the Navier-stokes viscosity term, while in the modified Darcy's law (Brinkman model), suggested by Brinkman [9], the Navier-stokes viscosity term still exists.Chen & Chen [10] considered the multi-layer problem when the above layer is heated and salted from above, and the solution of the problem is obtained using a shooting method.Lindsay & Ogden [11] worked in the implementation of spectral methods resistant to the generation of spurious eigenvalues.Lamb [12] used expansion of Chebyshev polynomials to investigate an eigenvalue problem arising from a model discussing a finitely conducting inner core of the earth on magnetically driven instability.Bukhari [13] studied the effects of surface-tension in a layer of conducting fluid with an imposed magnetic field and obtained results when the free surface is deformable and non-deformable.He solved that by using Chebyshev spectral method, and thus obtained some different results from that of Chen & Chen [10].Straughan [14] studied the thermal convection in fluid layer overlying a porous layer and considered the problem of lower layer heated from below and surface tension driven on the free top boundary of upper layer.In [15], he also dealt with the same problem considering the ratio depth of the relative layer and investigated the effect of the variation of relevant fluid and porous material properties.Allehiany [16] studied Benard convection in a horizontal porous layer

The Governing Equations
We consider a fluid layer overlying a porous layer so that the top of the porous layer touches the bottom of the fluid layer.The plane interface between the two layers is where the subscripts f and m denote the fluid layer and porous medium layer respectively.We suppose that the upper layer is filled with an incompressible thermally and electrically conducting viscous fluid consisting of melted salt which flows in it and governed by Navier-Stokes equations.However, the lower layer is occupied by a porous medium permeated by the fluid flowing in it and governed by Darcy-Brinkman's law.Both layers subjected to a constant vertical linear magnetic field and affected by a rotation around 3 x with a constant angular velocity .Gravity g acts in the negative direction of  3 x (see Figure 1).Convection is driven by the temperature depending on the fluid density and salting, and damped by viscosity.The Oberbeck-Boussineq approximation is used as the density of fluid is constant everywhere except in the body force term where the density is linearly proportional to temperature and salt concentration, i.e.

   
where T denotes the Kelvin temperature of the fluid, S is the salt concentration, 0  is the density of fluid at 0 and 0 , T S  (constant) is the thermal coefficient of volume expansion of the fluid and  (constant) is the salting coefficient of volume expansion of the fluid.Let be the solenoidal velocity of the fluid.

Let
and be respectively the solenoiddal velocity of the fluid, the magnetic field, the magnetic induction, the current density and the electric field.Hence , , , , , H B J and connected by the relations where m  is the magnetic permeability and  is the electrical conductivity.And the Maxwell equations curl , curl t where the displacement current has been neglected in the second of these Maxwell equations as is customary in situation when free charge is instantaneously dispersed.By substituting from (3), (4) 2 in to (4) 1 obtain by where 1 The equation of motion is where P is hydrostatic pressure,  is the dynamic viscosity and 2  is three-dimensional Laplacian operator.And by subsituting from (3) 1 , (4) 2 in to the Lorentz force Hence the equation of motion becomes and so the governing equations of the fluid layer are and the governing equations of the porous medium layer are where f P , m are the modified pressure of the fluid and the porous medium layers respectively and are the centrifugal force of the fluid and porous medium layer respectively, , f m V V are the solenoidal and seepage velocity respectively,   is impenetrable and maintained at constant temperature l and constant salt concentration , then the boundary conditions can be written as on the upper boundary, and  The boundary conditions on the interface plane x 3 = 0 are based on the assumption that temperature, salt concentration, heat flux, salt flux, normal and tangential fluid velocity, normal stress and tangential stress are continuous so that Equations ( 10)-( 19) have an equilibrium solution satisfying the boundary conditions ( 20)-( 22) on the form and with the boundary conditions and the interface conditions the equilibrium temperature field, hydrostatic pressure and salt concentration in the fluid layer and porous medium layer are respectively , 0

Perturbed Equations
We apply the perturbation by following linear perturbation quantities to the governing equations in the fluid layer and porous medium layer respectively and to the boundary conditions.After perturbation, the non-dimensionlisation will be apply by using , for the fluid layer, and by using for the porous medium layer, here where is the unit vector in the , and the Equations ( 15)-( 19) becomes where 1 when heating from below 1 when heating from above, 1 when salt concentration from below, 1 when salt concentration from above.

The Linearized Equations
Linearization will be done by neglecting all products and powers (higher than the first) of the linear perturbation quantity, and by dropping the (•) superscript, then by taking the curl of the Equations ( 30) and ( 35) we obtain if return to the original Equations ( 30) and ( 35), but in this case we take the (curl curl).Thus and if we use the curl of Equations ( 34) and (39) with using (33) and (38), we obtain Now, the third components of Equations ( 31), (32), (34), (36), (37), ( 39) and ( 41)-( 46) are Copyright x       is tow-dimensional Laplacian operator and .We apply the normal modes solution in the form w ,θ ,s ,ξ , j ,b ,w ,θ ,s ξ , and .Thus the governing equations are (66) , , and The boundary conditions in the final form are

Numerical Solution
A Legendre polynomials (see Bukhari [5]) is applied to solve the Equations ( 59)-( 70) with the relevant boundary conditions (71)-( 73), and we map by the transformations

  
Then the Equations ( 59)-( 70) can be rewritten in a system of twenty ordinary differential equations of first order, since and if we put σ so the eigenvalue problem can be reformulated in the form where A and B are real 28 × 28 matrices.The final eigenvalue problem reduces to V σ V  E F where matrices E and F have the block forms.The boundary conditions replace the 1 Mth, 2 Mth, … 28 Mth rows of E and F.

Results and Remarks
Using Legendre polynomials, the eigenvalue problems (59)-(70) with the boundary conditions (71)-( 73) are transformed to a system of fourteen ordinary differential equations of first order in the fluid layer and a system of fourteen ordinary differential equations of first order in the porous layer with twenty eight boundary conditions.In this work, we will discuss the numerical results through two cases-when the heat and salt concentration affected from above and below.
Case (1): the heat and the salt concentration affected from above.
Here, we put      Rt decreases which means that the increase of the rotation causes the increase of the thermal convections, leading to an increase in the instability of the fluid, as shown in Figure 2 and its numerical results in Table 1. The presence of the linear magnetic field helps reduce the currents of the thermal convections, meaning that the stability will increase in the fluid, as shown in Figure 3 and its numerical results in Table 2.  The fluid becomes unstable when the permeability of the porous medium increases, as shown in Figure 4 and its numerical results in Table 3.  The increase of depth ratio between the two layers makes the fluid unstable meaning that as d increases Rt m decreases, as shown in Figure 5 and its numerical results in Table 4. Case (2): the heat and the salt concentration affected from below.
Here, we put and the value of the initial salt Rayleigh number of a porous medium Rs m = 10000, to find the thermal Rayleigh numbers of a porous medium, m corresponding to wave numbers, for different values of In this case, the eigenvalues are complex, and thus the overstability happens, as shown in the following Ta- bles 5-8 and Figures 6-9.Therefore, we got the following results: P  The increase of the rotation helps reduce the currents of the thermal convections meaning that the stability will increase in the fluid, as displayed in Figure 6 and its numerical results in Table 5.  The presence of the linear magnetic field makes the fluid more stable as displayed in Figure 7 and its numerical results in Table 6.

Figure 1 .
Figure 1.Schematic diagram of the problem.
are the coriolis acceleration of the fluid and porous medium layer respectively, , f m T T are the Kelvin temperature of the fluid and porous medium layer respectively, , f m S S are salt concentration of the fluid and porous medium layer respectively, f m are the mass diffusivity of the fluid and porous medium layer respectively, and porous layer respectively, f m are the thermal and overall thermal conductivity of fluid and porous layers respectively, the heat and overall heat capacity per unit volume of the fluid and porous medium layers at constant pressure.In fact ) on the lower boundary, where f w and m are the normal axial velocity components of the fluid in fluid layer and porous medium layer respectively, 3 w f ζ and 3 m are the normal axial vortecity components of the fluid in.fluid layer and porous medium layer respectively. ζ Da, P rm , Rt m , Rs m , Ta m , Q m , Le m and m m P are non-dimensional numbers denote the Darcy number, viscous Prandtl number, thermal Rayleigh number, salt Rayleigh number, Taylor number, Chandarsekhar number, Lewis number and magnetic Prandtl number of the porous medium layer and given by 1and the value of the initial salt Rayleigh number of the porous medium Rs m = 5000 to find the thermal Rayleigh number of the porous medium m corresponding to the wave numbers m for the different values of m , m , m , this case, the eigenvalues are real, and thus the stationary instability happens, as shown in the following Tables 1-4 and Figures2-5.Therefore, we concluded that: P

Figure 2 .Figure 3 .
Figure 2. The relation between a m and Rt m for different values of Ta m when , Da m = 0.0003, Q m = 100, d . 0 01 f r f P Le   1 and

Figure 4 .
Figure 4.The relation between a m and Rt m for different values of Da m when Q m = 100, Ta m = 5, , d . 0 01 f r f P Le   1 and f m P  3 .

Figure 5 .
Figure 5.The relation between a m and Rt m for different values of when Q m = 100, Da m = 0.0003, Ta m = 5, and .d f 1 means that the increase of the permeability causes the increase of the thermal convections, leading to an increase in the instability of the fluid, as shown in Fig- ure8and its numerical results in Table7. As d increase m Rt increases meaning that the in-

Table 6 .
The relation between a m and Rt m for different values of Q m when Ta m = 0

Figure 6 .Figure 7 .Figure 8 .
Figure 6.The relation between a m and Rt m for different values of Ta m when , Q m = 100, Da m = 0.0001, and .d  0.2 f m P  3

Figure 9 .
Figure 9.The relation between a m and Rt m for different values of when Ta m = 0.5, Da m = 0.001, Q m = 100, d e the thermal diffusivity of the fluid sional numbers denote the viscous Prandtl number, thermal Rayleigh number, salt Rayleigh number, Taylor number, Chandraskhar number, Lewis number and magnetic Prandtl number of the fluid layer and given by f m are non-dimen- © 2012 SciRes.JEMAA 2012 SciRes.JEMAA

Table 1 . The relation between a m and Rt m for different va- lues of Ta m when , Da m = 0.0003, Q m = 100,
f m P  3 .m Rt 15

Table 2 . The relation between a m and Rt m for different va- lues of Q m when Ta m = 5, , Da m = 0.0003,
f m P  3 .m Rt 500

Table 5 . The relation between a m and Rt m for different va- lues of Ta m when , Q m = 100, Da m = 0.0003,
f m P  3 .m Rt 1

Table 7 . The relation between a m and Rt m for different va- lues of Da m when Ta m = 0.5, Q m = 50
f Le  1 f m P  3