Article Versions
Export Article
Cite this article
  • Normal Style
  • MLA Style
  • APA Style
  • Chicago Style
Research Article
Open Access Peer-reviewed

Stability Analysis of a Magneto Micropolar Fluid Layer by Variational Method

Joginder Singh Dhiman, Gurdeep Singh
American Journal of Applied Mathematics and Statistics. 2019, 7(3), 86-92. DOI: 10.12691/ajams-7-3-1
Received January 11, 2019; Revised February 28, 2019; Accepted April 04, 2019

Abstract

This Paper deals with the Stability of a micropolar fluid layer heated from below in the presence of uniform magnetic field. The basic hydrodynamic equations of magneto-microploar fluid layer heated from below are modified by using Boussinesq approximation and then the linearized perturbation equations are converted into a characteristic value problem with the help of Normal mode analysis. The expressions for Rayleigh number are obtained by using the variational principle. The effect of magnetic field and micropolar parameter on the Rayleigh numbers are discussed andupper bounds of critical Rayleigh number for all type of boundaries are obtained by using variational method.

1. Introduction

Thermal convection occurs in nature in so many forms and over such a wide range of scales that it could be claimed with some justification that convection represents the most common fluid flow in the universe.

The experiments by Bénard 1 particular have attracted great attention and are today considered as classical in fluid mechanics. The Bénard stability problem was first formulated and mathematically solved by Loard Rayleigh 2 for the case of free boundaries with a linear temperature gradient. A comprehensive account of the linearized stability theory of Rayleigh-Bénard convection problem in the presence of uniform rotation has been given in chanderasekhar 3.

Some years ago, Eringen 4 developed the theory of micropolar fluids in which local effects arising from micro-Structure and intrinsic motions of the fluid elements are taken into account. Ahmadi 5 studies the stability of a layer of micropolar fluid heated from below using linear theory as well as energy method the derived lower bounds for critical Rayleigh number. The profiles on Benard convection in micropoler fluids is also discussed by Narasimha, MY 6. Dattav and Sastry 7 also discussed the theory of micropolar fluid layer heated from below and obtained exact solution of the eigen value problem. Joginder Singh Dhiman, Praveenkumar Sharma and Gurdeep Singh.

8 studied the stability of micropolar fluid layer heated from below by variational principle.

The stability of magneto-micropolar fluid motion was studied by G. Ahmadi and M. Shahinpoor 9. The stability method employed by G. Ahmadi and M. Shahinpoor was an energy technique due to James Serrin 10.

Our aim in this paper is to verify the stability of magneto-micropolar fluid motion by using the variational method and we have established this by finding the upper bounds for critical Rayleigh numbers for all combinations of different type of boundaries.

2. Mathematical Analysis

2.1. Physical Problem

A viscous finitely heat and electrically conducting, micropolar fluid is statically confined between two horizontal boundaries Z = 0 and Z = d of infinite extension and finite vertical depth which are maintained at uniform temperatures T0 and T1 (T0> T1) respectively in the presence of uniform magnetic field acting antiparallel to the force field of gravity.

2.2. Hydrodynamic Equations and Boundary Conditions

The basic thermodynamic equations of the problem of thermal stability of magnetic micropolar fluid layer heated from below are modified by using the usual steps of Boussinesq of approximation and Normal mode. Analysis the non dimensional linearized perturbation equations with boundary conditions are as follows:

(1)
(2)
(3)
(4)

The boundary conditions are

(5)

(When both the boundaries are rigid).

(6)

(When both boundaries are free) and

(7)

(When both boundaries are either rigid or free, but boundaries are conducting).

In the foregoing equations;

k is dynamic micro rotation viscosity, μ is the dynamic Newtonain viscosity, p is hydrostatic pressure, J is microinertia, γis constant stands for coefficient of angular viscosity.

is the differentiation with respect to z, z in the real independent variable, a2 is the square of wave number, σis the thermal Prandtl number,is the heat diffusivity is the Rayleigh number, g is the gravitational Acceleration, α is the coefficient of tharmal gradient, d is the depth of layer, k1 is the thermometric conductivity, ν is the kinematic μ is viscosity, w, θ and G are the pertubations in vertical velocity, temperature and microrotation respectively and

The system of equations (1) — (4) together with boundary condition (5) — (7) constitutes an eigen value problem for R for given values of other parameters σ, c0, n1, a2K and Q.

2.3. Principle of Exchange of Stabilities (PES)

Firstly in the present problem, we shall see that whether the PES is valid or not for magneto micropolar fluid heated from below.

So multiplying equation (1) by w* (the complex conjugate of w) and integrating over the range 0 ≤z≤ 1, a suitable number of times by using relevant boundary conditions (5)-(7) we get

(8)

Taking the complex conjugate of equation (3) and multiplying the resulting equation by G on both sides and integrating it over the range of z, a suitable number of times by using boundary conditions (5) – (7), we get

(9)

By using equation (9) is equation (8) and integrating by parts the resulting equation, we get

(10)

Where

Now equating real and imaginary parts of both sides of equation (10) and canceling pi (≠ 0) (supposition) throughout from imaginary part, we get;

Now, multiplying the equation (4) by hz*(the complex conjugate of hz) and integrating the resulting equation over the vertical range of z, a suitable number of times and making use of either boundary conditions (5) – (7), we have.

(12)

Equating real parts from both sides of above equation, we have:

(13)

Since pr≥ 0 from above inequality, we have

(14)

Further, since w(0) = 0 = w(1), we have theReyleigh – Ritz inequality namely;

(15)

Using inequalities (14) and (15) is equation (11), we get

(16)

It clearly follows from above inequality that

Which is a contradiction to our suppositon that pi≠ 0. Hence, we must have pi = 0, In particular pr = 0 ⇒ pi = 0, hence PES is Valid.

Therefore, we have the result that PES is not valid or have over stability when Qσ1 > which can also be seen in classical magneto hydrodynamic Bénard convection.

2.4. Variational Principle

When instability sets is as stationary convection, the marginal state will be characterized by p =0 and basic equations (1) – (4) reduce to

(17)
(18)
(19)
(20)

From equation (17), letting

(21)

implies,

(22)

Operating the above equation lay (D2 – a2) and using equation (18), we get:

(23)

Multiplying the above equation by F and integrating the resulting equation over the range of z, we have

(24)

Since F(0) = 0 = F(1)because of equation (22), therefore, we have

(25)

Now integrating the equation (24) by using F from equation (22) and integrating a suitable number of times by using relevant boundary conditions, we get:

(26)

By using equation (19) in above equation and integrating the resulting equation by parts a suitable number of times and using boundary conditions on G, we get

(27)

Now, by using all these equations (25), (27),in the equation (24), we have :

Or

which may be written as

(28)

Where

And

2.5. Stationary Property

The stationary property is R given by equation (28) is checkedjust by giving a small variation δW to w, δG to G and δhz to hz, which are again compatible with boundary conditions i.e.

This can be easily checked that δR = 0 if (D2 – a2) F = and conversely if δR = 0 for any arbitrary Variation, then (D2 – a2) F =

2.6. Minimum Property

Now, we shall show that the lowest characteristic value of R is indeed, a true minimum.

Let Ri be a characteristic value and let corresponding characteristic function be distinguished by a subscript i. Then from equation (23), we have

(29)

and

(30)

and either

(31)

Multiplying equation (29) by Fj and integrating the resulting equation over vertical range of z, we have

By using equation (21), we have

Integrating the integrants on right hand side of above equation by parts, a suitable number of times with the help of boundary conditions(30) – (31), We have

(32)

Now, using this equation in (29) and interchanging i and j then subtracting the resulting equations, we have

Since, RiRjwe have;

The functions Wi form an orthogonal set.

Let

(33)

Similarly using the expanded forms Witerms of basic set functions for G, F, and hz, the values of integrals I1 and I2and equation (28) having the following values,

also, I2 =

Using these values of I1 and I2 is equations (28), we have

(34)

When .

In follows from equation (34) that

Since

Thus, we have

(35)

This shows that the quantity on R.H.S has a true minimum.

2.7 Upper Bounds for Critical Raylight Number for All Combination of Boundary Conditions

Let

(36)

Which obviously satisfies the boundary conditions

where the origin has been shifted to midway for convenience in computations.

Using the value of F from equation (36) in equation (21), we obtain

(37)

Operating the above equation by and using the equations (19) – (20), We get

(38)

The general solution of above differential equation is given by

(39)

Where are roots of auxiliary equation of (38) and

Case I: when both the boundaries are dynamically free.

By using equation (39), we have

(40)
(41)

And

(42)

Since, both the boundaries are dynamically free, so boundary conditions are W = 0 D2 w at z = and at z =.

By using these boundary conditions is above equations, we have

B1 = B2 = B3 = 0.

From equation (39), we get when

Now, evaluating I1 and I2 is equation (28), we have

and

Using these values of I1 and I2 in equations (35), we get

Now, using the value of A, we have

letting k →∞ in the above inequality, we have :

(43)

Which is the upper bound for Rcfor free-free boundaries,


Case II One Rigid – One Free Boundary:

Let us take lower boundary as rigid and upper dynamically free, so using these boundary conditions in the equations (39) – (41) we get:

(44)
(45)
(46)

On solvingthe above equations, we have :

Where

Now using the values of W and F is integrals I1 & 12, we have

(47)

By using these values of I1 & I2 is inequality (35), We have

Now, by using the values of B1, B2, B3, Δ A and letting K →∞ in above inequality, we have

(48)

Which is upper bound for R, for one rigid-one free boundaries.


Case III. For both rigid-rigid boundaries:

Now by using equations (39), (40) and (42), with boundary conditions in w, Dw and θ(proceeding same as in the previous two cases), We have:

(49)

Which is upper bound for Rc for both rigid-rigid boundaries.

3. Conclusion

In this paper it is verified that PES is valid for magneto-micropolar fluid motion for all type of boundary conditions when heated from below. There is over stability only when Qσ1 >, which can also be seen in classical magneto hydrodynamic Bénard convection. The upper bounds for Critical Rayleigh number, that isare obtained by using variational principle having same value for all type of boundary conditions namely, both free-free, one rigid-one free and both rigid-rigid. Further, these upper bounds depend upon the value of micro rotation coefficient K in magneto-micropolar fluid heated from below. The inequality for upper bounds clearly shows that the increasing value of K increases the Rayleigh number R. Thus K plays a stabilizing role in the system.

References

[1]  Benard, H.; Les tourbillons cellulairs dons une nappe liquids; Revue, generade des sciences pures et Applications, 1900: 11, 1261 & 1309.
In article      
 
[2]  Lord Rayleigh: On the connective currents in a horiyontal layer of fluid when higher temperature is one the upper side; phil.plag., 1990, 32, 529.
In article      
 
[3]  Chandrasekher S.. Hydrodynamic and Hydromagnetic stability, Oxford University Press, Aman House, London, 1961.
In article      
 
[4]  Eringen, A.C.: ’Theory of Micropolar Bluids; J. Math. Mech. 1906, 16, 1, 1-18.
In article      View Article
 
[5]  Ahmadi, G.: Stability of a micropolar fluid Layer Heated from Below’, Int. J. Engng. Sci., 1976, 14, 81-89.
In article      View Article
 
[6]  Narasimha, M.Y.: ‘Analysis of non uniform temp/ profiles on Benard convection in micropolar fluids; App. Math. comput, 2003: 134, 473-486.
In article      View Article
 
[7]  Datta, A. B. and Sastry, V.U.K. ‘Thermal Notability of a Horigontal Layer of, Micropolar fluid Heated from Bolows, Int. J. Engng. Sci. 1976:14, 631-637.
In article      View Article
 
[8]  Joginder S. Dhiman, Praveen K. Sharma and Gurdeep Singh: ‘Convective stability analysis of a micropolar fluid layer by variational method’ The Chinese Society of Theoretical and Applied Machanics Letters; 2011.
In article      View Article
 
[9]  Ahmadi, G. and Shahinpoor, M.; ‘Universal stability of magneto-micropolar fluid motions; Int. J. Engn. Sci, 1947, 12, 657-663.
In article      View Article
 
[10]  J.S. ERRIN, Arch, ration. Mech. Analysis 1950, 3, 1.
In article      
 

Published with license by Science and Education Publishing, Copyright © 2019 Joginder Singh Dhiman and Gurdeep Singh

Creative CommonsThis work is licensed under a Creative Commons Attribution 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Cite this article:

Normal Style
Joginder Singh Dhiman, Gurdeep Singh. Stability Analysis of a Magneto Micropolar Fluid Layer by Variational Method. American Journal of Applied Mathematics and Statistics. Vol. 7, No. 3, 2019, pp 86-92. http://pubs.sciepub.com/ajams/7/3/1
MLA Style
Dhiman, Joginder Singh, and Gurdeep Singh. "Stability Analysis of a Magneto Micropolar Fluid Layer by Variational Method." American Journal of Applied Mathematics and Statistics 7.3 (2019): 86-92.
APA Style
Dhiman, J. S. , & Singh, G. (2019). Stability Analysis of a Magneto Micropolar Fluid Layer by Variational Method. American Journal of Applied Mathematics and Statistics, 7(3), 86-92.
Chicago Style
Dhiman, Joginder Singh, and Gurdeep Singh. "Stability Analysis of a Magneto Micropolar Fluid Layer by Variational Method." American Journal of Applied Mathematics and Statistics 7, no. 3 (2019): 86-92.
Share
[1]  Benard, H.; Les tourbillons cellulairs dons une nappe liquids; Revue, generade des sciences pures et Applications, 1900: 11, 1261 & 1309.
In article      
 
[2]  Lord Rayleigh: On the connective currents in a horiyontal layer of fluid when higher temperature is one the upper side; phil.plag., 1990, 32, 529.
In article      
 
[3]  Chandrasekher S.. Hydrodynamic and Hydromagnetic stability, Oxford University Press, Aman House, London, 1961.
In article      
 
[4]  Eringen, A.C.: ’Theory of Micropolar Bluids; J. Math. Mech. 1906, 16, 1, 1-18.
In article      View Article
 
[5]  Ahmadi, G.: Stability of a micropolar fluid Layer Heated from Below’, Int. J. Engng. Sci., 1976, 14, 81-89.
In article      View Article
 
[6]  Narasimha, M.Y.: ‘Analysis of non uniform temp/ profiles on Benard convection in micropolar fluids; App. Math. comput, 2003: 134, 473-486.
In article      View Article
 
[7]  Datta, A. B. and Sastry, V.U.K. ‘Thermal Notability of a Horigontal Layer of, Micropolar fluid Heated from Bolows, Int. J. Engng. Sci. 1976:14, 631-637.
In article      View Article
 
[8]  Joginder S. Dhiman, Praveen K. Sharma and Gurdeep Singh: ‘Convective stability analysis of a micropolar fluid layer by variational method’ The Chinese Society of Theoretical and Applied Machanics Letters; 2011.
In article      View Article
 
[9]  Ahmadi, G. and Shahinpoor, M.; ‘Universal stability of magneto-micropolar fluid motions; Int. J. Engn. Sci, 1947, 12, 657-663.
In article      View Article
 
[10]  J.S. ERRIN, Arch, ration. Mech. Analysis 1950, 3, 1.
In article