Effect of Non-Uniform Basic Concentration Gradient on the Onset of Double-Diffusive Convection in Micropolar Fluid

The effect of non-uniform basic concentration gradient on the onset of double diffusive convection in a micropolar fluid layer heated and saluted from below and cooled from above has been studied. The linear stability analysis is performed. The eigen value of the problem is obtained using Galerkian method. The eigen values are obtained for 1) free-free; 2) rigid-free; 3) rigid-rigid velocity boundary combination with isothermal temperature condition on spin-vanishing permeable boundaries. The influence of various micropolar parameters on the onset of convection has been analyzed. One linear and five non linear concentration profiles are considered and their comparative influence on onset is discussed and results are depicted graphically. It is observed that fluid layer with suspended particles heated and soluted from below is more stable compare to the classical fluid without suspended particles.


Introduction
The convection driven by two different density gradients with differing rates of diffusion is widely known to as "double-diffusive convection" and is an important fluid dynamics phenomenon (see Mojtabi and Charrier-Mojtabi [1]).The study of double-diffusive convection has attracted attention of many researchers during the recent past due to its occurrence in nature and industry.Oceanography is the root of double-diffusive convection in natural settings.The existence of heat and salt concentrations at different gradients and the fact that they diffuse at different rates lead to spectacular double-diffusive instabilities known as "salt-fingers" (see Stern [2,3]).The formation of salt-fingers can also be observed in laboratory settings.Double-diffusive convection occurs in the sun where temperature and Helium diffusions take place at different rates.Convection in magma chambers and sea-wind formations are among other manifestations of double-diffusive convection in nature.The theory of double-diffusive convection both theoretically and experimentally is investigated by Turner [4], Chen et al. [5] and recently Malashetty and Bharati [6].
The theory of micro-fluid introduced by Eringen [7] deals with a class of fluids, which exhibits certain microscopic effects arising from the local structure and mi-cro-motions of the fluid elements.Consequently new principles must be added to the basic principle of continuous media which deals with 1) Conservation of micro inertia moments, and 2) Balance of first stress moments.The theory of micro fluids naturally give rise to the concept of inertial spin, body moments, micro-stress averages and stress moments which have no counter part in the classical fluid theories.A detailed survey of the theory of micropolar fluid and its applications are considered in the books of Erigen [8,9], Lukasazewicz [10] and power [11].The theory of thermomicropolar convection was studied by many authors Datta and Sastry [12], Ahmadi [13], Rama Rao [14], Lebon and Perez-Gracia [15], Bhattacharya and Jena [16], Siddheshwar and Pranesh [17][18][19][20], Y. N. Murthy [21], and more recently by Pranesh and Kiran [22] and Alloui et al. [23].
The main object of this paper is to study the onset of double diffusive convection in a micropolar fluid using one linear and five non-linear concentration profiles as a mechanism to control the onset of convection.ture and species concentration of the fluid between lower and upper plates.Appropriate single-phase heat and solute transport Equations are chosen with effective heat and capacity ratio and effective thermal diffusivity.A Cartesian coordinate system (x, y, z) is used with the origin at the bottom of the fluid layer and the z-axis vertically upwards.

Mathematical Formulation
The governing Equations for the Rayleigh-Benard situation in a Boussinesquian micropolar fluid are; Continuity Equation: Conservation linear momentum: Conservation of angular velocity: Conservation of energy: Conservation soluted concentration: Equation of state: Here, is the velocity, o q  is density of the fluid at temperature T = T 0 , p is the pressure,  is the density, is acceleration due to gravity, g  is coupling viscosity coefficient or vortex viscosity,  and  are the bulk and shear spin-viscosity coefficients,  is the an- gular velocity, I is moment of inertia,  and  are bulk and shear spin-viscosity coefficients, T is the temperature,  is the thermal conductivity,  is micropolar heat conduction coefficient, s  is solute thermal conductivity, C is concentration,  is coefficient of thermal expansion,  is electrical conductivity.

Basic State
The basic state of the fluid quiescent and is given by where the subscript "b" denotes the basic state.The non-dimensional basic concentration gradient g(z) which is non-negative satisfies the condition One linear and five non-linear basic concentration gradients are considered in this paper to study the onset of convection and are given below in Table 1.

Linear Stability Analysis
We now superpose infinitesimal perturbations on the quiescent basic state and study the instability.Let the basic state be disturbed by an infinitesimal thermal perturbation, given by; The prime indicates that the quantities are infinitesimal perturbations.In the present problem we assume that principle of exchange of stability to be valid and hence deal with only stationary convection.Substituting Equation (9) into Equations ( 1) to ( 6) and using the basic state (8) we get, the linearized Equations governing the infinitesimal perturbations in the form: The perturbation Equations ( 10)-( 15) are non-dimensionalised using the following definitions: Using Equation (15) in Equation ( 11), operating curl twice on the resulting Equation, and operating curl on Equation (12), we get where, the asterisks have been dropped for simplicity and non-dimensional parameters are The infinitesimal perturbation , W , z  T and C are assumed to be periodic waves (see Chandrasekhar 1961) and hence these permit a normal mode solution in the form: where and are horizontal components of the wave number , substituting Equation ( 21) into Equations ( 17)-( 20), we get and a l m   .The sets of ordinary differential Equations ( 22)-( 25) are approximations based on physical considerations to the system of partial differential Equations ( 17)- (20).Although the relationship between the solutions of the governing partial differential Equations and the corresponding ordinary differential Equations has not been established, these linear models reproduce qualitatively the convective phenomena observable through the full system.
When the fluid layer is heated from below, the nonuniform concentration gradient g(z) is not only non-negative but also decreases monotonically.Thus our main interest here is to find out suitable non-uniform basic concentration profile that gives the maximum/minimum critical Rayleigh number, for this we have considered the various concentration profile as shown in the Table 1.
In the Galerkin procedure, we expand the velocity, microrotation, temperature, concentration by, are polynomials in z that generally have to satisfy the given boundary conditions.For the single term Galerkin expansion technique we take i = j = 1.Multiplying Equation ( 22) by W, Equation ( 23) by G, Equation (24) by T and, and Equation (25) by C, integrating the resulting integrals by parts with respect to z from 0 to 1 and taking 1 , where In the Equation ( 26),    denotes integration with respect to z between z = 0 and z = 1.We note here that R in Equation ( 30) is a functional and the Euler-Lagrange Equations for the extremisation of R are Equations ( 22) to (25).
The value of critical Rayleigh number depends on the boundaries.In this paper we consider following boundary combinations.
1) Free-Free Isothermal-Permeable No-Spin condition: 3) Rigid-Rigid Isothermal-Permeable No-Spin condition: The trial functions satisfying the above boundary conditions are presented in Table 2.

Result and Discussion
In the paper, we study the effect of basic non-uniform concentration gradient on the convective instability of a micropolar fluid.One uniform, five non-uniform basic concentration gradient are chosen for the study and this have been presented in Table 1.We find that,

B.Cs
Free-free Rigid-free Rigid-rigid for symmetric boundary combination(free-free and rigidrigid).In the case of non-symmetric boundary combination namely, rigid-free we find that, i.e., salting from above and step function profiles are most destabilizing and stabilizing respectively.Figures 2(a), 3(a) and 4(a) are the plots of critical Rayeigh number R c versus coupling parameter N 1 for different values of solutal Rayleigh number, R s and different non-uniform basic concentration profiles, for freefree, rigid-free and rigid-rigid isothermal permeable, nospin boundary condition.It is observed that as N 1 increases, R c also increases.Increase in N 1 indicates the increase in the concentration of the microelements.These microelements consume the greater part of the energy in developing gyrational velocity and as a result the onset of convection is delayed.From this we conclude that an increase in N 1 is to stabilize the system.
Figures 2(b), 3(b) and 4(b) are the plots of critical Rayeigh number R c versus couple stress parameter N 3 for different values of solutal Rayleigh number, R s and different non-uniform basic concentration profiles, for freefree, rigid-free and rigid-rigid isothermal permeable, nospin boundary condition.We note that the role played by the shear stress in the conservation of linear momentum is played by the couple stress in angular momentum Equations.It is observed that as N 3 increases the R c decreases, because, when N 3 increases the couple stress of the fluid increases, which causes the microrotation to decease.Therefore, increase in N 3 destabilize the system hence, we can conclude that couple stress are operative only for small values of N 3 .
Figures 2(c), 3(c) and 4(c) are the plots of critical Rayeigh number R c versus micropolar heat conduction parameter N 5 for different values of solutal Rayleigh number, R s and different non-uniform basic concentration profiles, for free-free, rigid-free and rigid-rigid isothermal permeable, no-spin boundary condition.When N 5 increases, the heat induced into the fluid due to these microelements also increases, thus reducing the heat transfer from bottom to top.The decrease in heat transfer is responsible for delaying the onset of instability.Thus Copyright © 2012 SciRes.AM N 5 stabilizes the system.Figure 5 is the plot of critical Rayeigh number R c versus diffusivity ratio г.The increase in the value of г decreases R c , thus advancing the onset of double diffusive convection.
From the above figures we observed that increase in solutal Rayleigh number R s increases R c , indicating that the additional diffusing component stabilizes the system.
From the figures we also observe that , where superscript denotes the velocity boundary combinations.

Conclusions
The following conclusions are drawn from the study: 1) By choosing appropriate non-uniform basic temperature profile it is possible to control the Rayleigh-Bénard convection.
2) The step function is the most stabilizing basic concentration distribution and salting from above is the most destabilizing basic concentration distribution.
3) The effect of increase in the value of solute Rayliegh number is to increase the value of R c , thus making the system more stable.
4) It is found that the increase in the value of diffusivity ratio, decrease the value of R c , thus advancing the onset of double diffusive convection.
5) The coupling parameter N 1 and micropolar heat conduction parameter N 5 stabilizes the system, whereas the couple stress parameter N 3 destabilizes the system.6) Rayleigh-Bénard double diffusive convection in Newtonian fluids may be delayed by adding micron sized suspended particles.

Figure 1 .
Figure 1.Schematic diagram for the problem.

Figure 2 .Figure 3 .Figure 4 .Figure 5 .
Figure 2. Plot of Critical Rayleigh number versus (a) Coupling parameter N 1 , (b) Couple stress parameter N 3 , (c) Micropolar heat conduction parameter N 5 for various Solutal Rayleigh number R S for the six models as in the Table 1 in the case of free-free isothermal-permeable no spin condition (---R S = 25 & -R S = 200).

Table 1 . Non-uniform basic concentration profiles.
AMare constants and W 1 , G 1 , T 1 , C 1 are trial functions, yield the following Equation for the Rayleigh number, which A, B, E and F Copyright © 2012 SciRes.