Research Article  Open Access
M. N. Mahmud, R. Idris, I. Hashim, "Effects of Magnetic Field and Nonlinear Temperature Profile on Marangoni Convection in Micropolar Fluid", International Journal of Differential Equations, vol. 2009, Article ID 748794, 11 pages, 2009. https://doi.org/10.1155/2009/748794
Effects of Magnetic Field and Nonlinear Temperature Profile on Marangoni Convection in Micropolar Fluid
Abstract
The combined effects of a uniform vertical magnetic field and a nonuniform basic temperature profile on the onset of steady Marangoni convection in a horizontal layer of micropolar fluid are studied. The closedform expression for the Marangoni number for the onset of convection, valid for polynomialtype basic temperature profiles upto a third order, is obtained by the use of the singleterm Galerkin technique. The critical conditions for the onset of convection have been presented graphically.
1. Introduction
Convective flow in a thin layer of fluid, free at the upper surface and heated from below, is of fundamental importance and a prototype to a more complex configuration in experiments and industrial processes. The convective flows in a liquid layer can be driven by buoyancy forces due to temperature gradients and/or thermocapillary (Marangoni) forces caused by surface tension gradients. Thermal convective problems have long been studied extensively since the pioneering experimental and theoretical works of Bénard [1], Rayleigh [2], and Pearson [3]. The instability problems have been studied in several other directions (cf. [4–18]).
Most of the previous studies were concerned with convection in Newtonian fluids. However, much less work has been done on convection in nonNewtonian fluids such as the micropolar fluids. The theory of micropolar fluids, as developed by Eringen [19], has been a field of sprightly research for the last few decades especially in many industrially important fluids like paints, polymeric suspensions, colloidal fluids, and also in physiological fluids such as normal human blood and synovial fluids. Rama Rao [20] studied the effect of a magnetic field on convection in a micropolar fluid. The onset of convection as overstable motions in a micropolar fluid was examined in [21]. Sharma and Gupta [22] studied convection in micropolar fluids in a porous medium. Ramdath [23] considered buoyancyand thermocapillarydriven (BénardMarangoni) convection in a layer of micropolar fluid. The effect of throughflow on Marangoni convection in micropolar fluids was analyzed in [24]. Siddheshwar and Sri Krishna [25] presented both linear and nonlinear analyses of convection in a micropolar fluid occupying a porous medium. Sunil et al. [26] studied the effect of rotation on convection in a micropolar ferrofluid.
There has also been much less work focused on the effect of nonuniform temperature gradient on convection. Friedrich and Rudraiah [27] studied the combined effects of nonuniform temperature gradients and rotation on Marangoni convection. The combined effects of nonuniform temperature gradients and a magnetic field on Marangoni convection were investigated by Rudraiah et al. [28]. The work of Friedrich and Rudraiah [27] was further extended to include the effect of buoyancy by Rudraiah and Ramachandramurthy [29]. Dupont et al. [30] studied the effect of a cubic quasisteady temperature profile on Marangoni convection. The effects of nonuniform temperature gradients on the onset of oscillatory Marangoni and BénardMarangoni convection in a magnetic field were analyzed in [31, 32], respectively. Chiang [33] investigated the effect of Dupont et al. [30] temperature profile on the onset of stationary and oscillatory BénardMarangoni convection.
Thermal convection in micropolar fluids has also been studied. Rudraiah and Siddheshwar [34] analyzed the effects of nonuniform temperature gradients of parabolic and stepwisetypes on the onset of Marangoni convection in a micropolar fluid. This study was later extended by Siddheshwar and Pranesh [35] to include the effect of a magnetic field and buoyancy forces. Very recently, Idris et al. [36] studied the effect of Dupont et al. [30] cubic temperature profile on the onset of BénardMarangoni convection in a micropolar fluid.
In this paper, we shall investigate the combined effects of Dupont et al. [30] cubic temperature profile and a magnetic field on the onset of Marangoni convection in a micropolar fluid. The singleterm Galerkin technique [37] is employed to obtain a closedform expression of (Marangoni number) for the onset of convection. Comparisons with the other polynomialtype temperature profiles normally used by previous investigators shall be undertaken.
2. Mathematical Formulation
We wish to examine the stability of a horizontal layer of quiescent micropolar fluid of thickness in the presence of a magnetic field. We assume that the layer is bounded below by a rigid boundary, which is kept at a constant temperature, and above by a perfectly insulated, flat free surface. Moreover, the spinvanishing boundary condition is assumed at the boundaries.
The governing equations for the problem are the continuity equation, conservation of momentum, conservation of angular momentum, conservation of energy, and magnetic induction, compare [19, 34, 35]: where is the velocity, is the spin, is the temperature, is the magnetic field, is the hydromagnetic pressure, is the coupling viscosity coefficient, is the shear kinematic viscosity coefficient, is the moment of inertia, and are the bulk and shear spin viscosity coefficients, is the micropolar heat conduction coefficient, is the specific heat, is the thermal conductivity, and is the magnetic viscosity (where electrical conductivity and magnetic permeability). All the viscosity coefficients, heat conduction coefficient and thermal conductivity are thermodynamically restricted on the assumption of ClausiusDuhem inequality (see Eringen [19]) and are all positive quantities.
The surface tension at the free upper surface is
where is the unperturbed value of and . The perturbation (2.1) are nondimensionalised using the following definition:
Following the classical lines of linear stability theory, the linearised and dimensionless governing equations are
where , , , and are, respectively, the amplitudes of the infinitesimal perturbations of velocity, spin, temperature, and magnetic field, is the coupling parameter ( ), is the couple stress parameter (, : finite, real), is the micropolar heat conduction parameter (, : finite, real), is the Chandrasekhar number, is the Prandtl number, is the magnetic Prandtl number, and is a nondimensional basic temperature gradient satisfying the condition .
The infinitesimal perturbations , , , and are assumed to be periodic waves and hence these permit a normal mode solution in the following form:
where and are horizontal components of the wave number .
Substituting (2.5) into (2.4), we get where .
Eliminating between (2.6) and (2.9), we obtain
Equations (2.7), (2.8), and (2.10) are solved subject to the linearized and dimensionless boundary conditions: where is the Marangoni number (where is the temperature difference between the two boundaries).
Following [30], we consider the steady state temperature profile given by
which precisely represents an experimental data, where (^{}) denotes dimensional quantities, is the temperature at the upper free surface, and , are constants. In nondimensional form, the in this case is given by
The case , , and recovers the classical linear basic state temperature distribution. The different temperature gradients studied in this paper are listed in Table 1.

3. Solution of the Linearized Problem
Equations (2.7), (2.8), and (2.10) subject to the boundary conditions (2.11) constitute an eigenvalue problem. To solve the resulting eigenvalue problem, a singleterm Galerkin expansion technique [37] is used to encompass a vast parameter space. Also, the technique employed yields sufficiently accurate and useful results for the purpose in hand with minimum of mathematics [37].
First we multiply (2.7), (2.8) and (2.10) by , , and , respectively. Then we integrate the resulting equations by parts with respect to from 0 to 1. By using the boundary conditions (2.11) and taking , , and , and in which , , and are constants and , , and are trial functions, yields the eigenvalue in the form
where Now with as given in (2.13), we rewrite the expression (3.1) in the closedform expression for :
where We remark that (3.3) is valid for all polynomialtype basic temperature profiles up to a third order. The critical Marangoni number, , for the onset of convection is the global minimum of over .
4. Discussion
The critical Marangoni number which attains its minimum at is computed from (3.3) for different volumes of , , , and and the results are depicted in Figures 1, 2, and 3. We recover the results of Rudraiah and Siddheshwar [34] for the linear and inverted parabolic temperature gradients when . We observe that as or increases, also increases. Obviously, the onset of convection will be delayed by increasing the concentration of the microelements or heat induced into the fluid by the microelements. But, an increase in leads to a decrease in microrotation, and hence the system becomes more unstable. Also it is observed that Model 4 (Cubic 2), with , , as used by Dupont et al. [30], is less stabilizing than Model 2 (Inverted parabolic), that is, . Based on our results, Model 3 (Cubic 1) with , , is shown to be the most stabilizing of all the considered types of temperature gradients, that is, .
Figures 4–6 illustrate the variations of the critical Marangoni number with the Chandrasekhar number for some assigned values of , , and , respectively. The results indicate that is generally an increasing function of . From Figure 4, we notice that the increase in the concentration of the microelements is to stabilize the system by superposing on the effect of the magnetic field. Figure 5 shows that the effect of on the system is very small compared to the effects of the other microelements. As before, Model 3 (Cubic 1) with , , is shown to be the most stabilizing of all the considered types of temperature gradients, that is, .
5. Conclusion
The problem of Marangoni convection in a micropolar fluid in the presence of a cubic basic state temperature profile and a vertical magnetic field has been studied theoretically. The results indicate that it is possible to delay the onset of convection by the application of a cubic basic state temperature profile. In addition, the presence of a magnetic field is to suppress Magnetomarangoni convection and hence leads to a more stable system. As expected, the presence of the micronsized suspended particles adds to the stabilizing effect of the magnetic field.
Acknowledgment
The authors acknowledge the financial support received under the Grant UKMGUPBTT0725173 and from Universiti Kuala Lumpur (UniKL MICET).
References
 H. Bénard, “Les tourbillons cellulaires dans une nappe liquide,” Revue Générale des Sciences Pures et Appliquées, vol. 11, pp. 1261–1271, 1900. View at: Google Scholar
 L. Rayleigh, “On convection currents in a horizontal layer of fluid when the higher temperature is on the other side,” Philosophical Magazine, vol. 32, pp. 529–546, 1916. View at: Google Scholar
 J. R. A. Pearson, “On convection cells induced by surface tension,” Journal of Fluid Mechanics, vol. 4, pp. 489–500, 1958. View at: Google Scholar
 S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, The International Series of Monographs on Physics, Clarendon Press, Oxford, UK, 1961. View at: Zentralblatt MATH  MathSciNet
 D. A. Nield, “Surface tension and buoyancy effects in the cellular convection of an electrically conducting liquid in a magnetic field,” Zeitschrift für Angewandte Mathematik und Physik, vol. 17, no. 1, pp. 131–139, 1966. View at: Publisher Site  Google Scholar
 M. Takashima, “Nature of the neutral state in convective instability induced by surface tension and buoyancy,” Journal of the Physical Society of Japan, vol. 28, p. 810, 1970. View at: Google Scholar
 S. H. Davis and G. M. Homsy, “Energy stability theory for freesurface problems: buoyancythermocapillary layers,” Journal of Fluid Mechanics, vol. 98, no. 3, pp. 527–553, 1980. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M.I. Char and K.T. Chiang, “Boundary effects on the BénardMarangoni instability under an electric field,” Applied Scientific Research, vol. 52, no. 4, pp. 331–354, 1994. View at: Publisher Site  Google Scholar
 I. Hashim and S. K. Wilson, “The effect of a uniform vertical magnetic field on the onset of oscillatory Marangoni convection in a horizontal layer of conducting fluid,” Acta Mechanica, vol. 132, no. 1–4, pp. 129–146, 1999. View at: Google Scholar
 I. Hashim, “On competition between modes at the onset of BénardMarangoni convection in a layer of fluid,” The Australian & New Zealand Industrial and Applied Mathematics Journal, vol. 43, no. 3, pp. 387–395, 2002. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 I. Hashim and N. Md Arifin, “The effect of a magnetic field on the linear growth rates of BénardMarangoni convection,” Microgravity Science and Technology, vol. 17, no. 2, pp. 5–8, 2005. View at: Google Scholar
 I. Hashim and Z. Siri, “Stabilization of steady and oscillatory marangoni instability in rotating fluid layer by feedback control strategy,” Numerical Heat Transfer A, vol. 54, no. 6, pp. 647–663, 2008. View at: Publisher Site  Google Scholar
 S. Awang Kechil and I. Hashim, “Control of Marangoni instability in a layer of variableviscosity fluid,” International Communications in Heat and Mass Transfer, vol. 35, no. 10, pp. 1368–1374, 2008. View at: Publisher Site  Google Scholar
 Z. Siri and I. Hashim, “Control of oscillatory Marangoni convection in a rotating fluid layer,” International Communications in Heat and Mass Transfer, vol. 35, no. 9, pp. 1130–1133, 2008. View at: Publisher Site  Google Scholar
 S. Awang Kechil and I. Hashim, “Oscillatory Marangoni convection in variableviscosity fluid layer: the effect of thermal feedback control,” International Journal of Thermal Sciences, vol. 48, no. 6, pp. 1102–1107, 2009. View at: Publisher Site  Google Scholar
 W.M. Yang, “Thermal instability of a fluid layer induced by radiation,” Numerical Heat Transfer A, vol. 17, no. 3, pp. 365–376, 1990. View at: Google Scholar
 A. Y. Gelfgat and I. Tanasawa, “Numerical analysis of oscillatory instability of buoyancy convection with the Galerkin spectral method,” Numerical Heat Transfer A, vol. 25, no. 6, pp. 627–648, 1994. View at: Google Scholar
 E. EvrenSelamet, V. S. Arpaci, and A. T. Chai, “Thermocapillarydriven flow past the Marangoni instability,” Numerical Heat Transfer A, vol. 26, no. 5, pp. 521–535, 1994. View at: Google Scholar
 A. C. Eringen, “Micropolar theory of liquid crystals,” in Liquid Crystals and Ordered Fluids, J. F. Johnson and R. S. Porter, Eds., vol. 3, Plenum, New York, NY, USA, 1978. View at: Google Scholar
 K. V. Rama Rao, “Thermal instability in a micropolar fluid layer subject to a magnetic field,” International Journal of Engineering Science, vol. 18, no. 5, pp. 741–750, 1980. View at: Google Scholar
 C. PérezGarcía and J. M. Rubí, “On the possibility of overstable motions of micropolar fluids heated from below,” International Journal of Engineering Science, vol. 20, no. 7, pp. 873–878, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. C. Sharma and U. Gupta, “Thermal convection in micropolar fluids in porous medium,” International Journal of Engineering Science, vol. 33, no. 13, pp. 1887–1892, 1995. View at: Google Scholar
 G. Ramdath, “BénardMarangoni instability in a layer of micropolar fluid,” Journal of NonEquilibrium Thermodynamics, vol. 22, no. 4, pp. 299–310, 1997. View at: Google Scholar
 Y. N. Murty and V. V. Ramana Rao, “Effect of throughflow on Marangoni convection in micropolar fluids,” Acta Mechanica, vol. 138, no. 3, pp. 211–217, 1999. View at: Google Scholar
 P. G. Siddheshwar and C. V. Sri Krishna, “Linear and nonlinear analyses of convection in a micropolar fluid occupying a porous medium,” International Journal of NonLinear Mechanics, vol. 38, no. 10, pp. 1561–1579, 2003. View at: Publisher Site  Google Scholar
 Sunil, P. Chand, P. K. Bharti, and A. Mahajan, “Thermal convection in micropolar ferrofluid in the presence of rotation,” Journal of Magnetism and Magnetic Materials, vol. 320, no. 34, pp. 316–324, 2008. View at: Publisher Site  Google Scholar
 R. Friedrich and N. Rudraiah, “Marangoni convection in a rotating fluid layer with nonuniform temperature gradient,” International Journal of Heat and Mass Transfer, vol. 27, no. 3, pp. 443–449, 1984. View at: Google Scholar
 N. Rudraiah, V. Ramachandramurthy, and O. P. Chandna, “Effects of magnetic field and nonuniform temperature gradient on Marangoni convection,” International Journal of Heat and Mass Transfer, vol. 28, no. 8, pp. 1621–1624, 1985. View at: Google Scholar
 N. Rudraiah and V. Ramachandramurthy, “Effects of nonuniform temperature gradient and Coriolis force on BénardMarangoni's instability,” Acta Mechanica, vol. 61, no. 1–4, pp. 37–50, 1986. View at: Publisher Site  Google Scholar
 O. Dupont, M. Hennenberg, and J. C. Legros, “MarangoniBénard instabilities under nonsteady conditions. Experimental and theoretical results,” International Journal of Heat and Mass Transfer, vol. 35, no. 12, pp. 3237–3244, 1992. View at: Google Scholar
 M.I. Char and C.C. Chen, “Effects of nonuniform temperature gradients on the onset of oscillatory Marangoni convection in a magnetic field,” Acta Mechanica, vol. 161, no. 12, pp. 17–30, 2003. View at: Publisher Site  Google Scholar
 M.I. Char and C.C. Chen, “Effect of a nonuniform temperature gradient on the onset of oscillatory BénardMarangoni convection of an electrically conducting liquid in a magnetic field,” International Journal of Engineering Science, vol. 41, no. 15, pp. 1711–1727, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 K.T. Chiang, “Effect of a nonuniform basic temperature gradient on the onset of BénardMarangoni convection: stationary and oscillatory analyses,” International Communications in Heat and Mass Transfer, vol. 32, no. 12, pp. 192–203, 2005. View at: Publisher Site  Google Scholar
 N. Rudraiah and P. G. Siddheshwar, “Effect of nonuniform basic temperature gradient on the onset of Marangoni convection in a fluid with suspended particles,” Aerospace Science and Technology, vol. 4, no. 8, pp. 517–523, 2000. View at: Google Scholar
 P. G. Siddheshwar and S. Pranesh, “Magnetoconvection in fluids with suspended particles under 1 g and $\mu $ g,” Aerospace Science and Technology, vol. 6, no. 2, pp. 105–114, 2002. View at: Publisher Site  Google Scholar
 R. Idris, H. Othman, and I. Hashim, “On effect of nonuniform basic temperature gradient on BénardMarangoni convection in micropolar fluid,” International Communications in Heat and Mass Transfer, vol. 36, no. 3, pp. 255–258, 2009. View at: Publisher Site  Google Scholar
 B. A. Finlayson, The Method of Weighted Residuals and Variational Principles, Academic Press, New York, NY, USA, 1972.
Copyright
Copyright © 2009 M. N. Mahmud et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.