Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com
|
customercare@cbspd.com
Department of Mathematics and Institute For Mathematical Research Universiti Putra Malaysia, 43400 UPM Serdang, Selangor, Malaysia
Haliza Rosali
Department of Mathematics and Institute For Mathematical Research Universiti Putra Malaysia, 43400 UPM Serdang, Selangor, Malaysia
Norfifah Bachok
Department of Mathematics and Institute For Mathematical Research Universiti Putra Malaysia, 43400 UPM Serdang, Selangor, Malaysia
*Address all correspondence to:
1. Introduction
Convection in a plane horizontal fluid layer heated from below, initially at rest and subject to an adverse temperature gradient, may be produced either by buoyancy forces or surface tension forces. These convective instability problems are known as the Rayleigh-Benard convection and Marangoni convection, respectively. The determination of the criterion for the onset of convection and the mechanism to control has been a subject of interest because of its applications in the heat and momentum transfer research. RayleighRayleigh 1916 was the first to solve the problem of the onset of thermal convection in a horizontal layer of fluid heated from below. His linear analysis showed that Benard convection occurs when the Rayleigh number exceeds a critical value. Theoretical analysis of Marangoni convection was started with the linear analysis by Pearson (1958) who assumed an infinite fluid layer, a nondeformable case and zero gravity in the case of no-slip boundary conditions at the bottom. He showed that thermocapillary forces can cause convection when the Marangoni number exceeds a critical value in the absence of buoyancy forces.
The determination of the criterion for the onset of convection and the mechanism to control convective flow patterns is important in both technology and fundamental Science. The problem of suppressing cellular convection in the Marangoni convection problem has attracted some interest in the literature. The effects of a body force due to an externally-imposed magnetic field on the onset of convection has been studied theoretically and numerically. The effect of magnetic field on the onset of steady buoyancy-driven convection was treated by ChandrasekharChandrasekhar 1961 who showed that the effect of magnetic field is to increase the critical value of Rayleigh number and hence to have a stabilising effect on the layer. The effect of a magnetic field on the onset of steady buoyancy and thermocapillary-driven (Benard-Marangoni) convection in a fluid layer with a nondeformable free surface was first analyzed by Nield Nield 1966. He found that the critical Marangoni number monotonically increased as the strength of vertical magnetic field increased. This indicates that Lorentz force suppressed Marangoni convection. Later, the effect of a magnetic field on the onset of steady Marangoni convection in a horizontal layer of fluid has been discuss in a series of paper by Wilson Wilson 1993, Wilson 1994. The influence of a uniform vertical magnetic field on the onset of oscillatory Marangoni convection was treated by Hashim & WilsonHashim & Wilson 1999 and Hashim & Arifin Hashim & Arifin 2003.
The present work attempts to delay the onset of convection by applying the control. The objective of the control is to delay the onset of convection while maintaining a state of no motion in the fluid layer. Tang and Bau Tang and Bau 1993, Tang and Bau 1994 and HowleHowle 1997)have shown that the critical Rayleigh number for the onset of Rayleigh-Bénard convection can be delayed. Or et al. Or et al. 1999 studied theoretically the use of control strategies to stabilize long wavelength instabilities in the Marangoni-Bénard convection. BauBau 1999 has shown independently how such a feedback control can delay the onset of Marangoni-Bénard convection on a linear basis with no-slip boundary conditions at the bottom. Recently, Arifin et. al.Arifin et. al. 2007 have shown that a control also can delay the onset of Marangoni-Bénard convection with free-slip boundary conditions at the bottom.
Therefore, in this paper, we use a linear controller to delay the onset of Marangoni convection in a fluid layer with magnetic field. The linear stability theory is applied and the resulting eigenvalue problem is solved numerically. The combined effect of the magnetic field and the feedback control on the onset of steady Marangoni convection are studied.
Consider a horizontal fluid layer of depth d with a free upper surface heated from below and subject to a uniform vertical temperature gradient. The fluid layer is bounded below by a horizontal solid boundary at a constant temperature T1 and above by a free surface at constant temperature T2 which is in contact with a passive gas at constant pressure P0 and constant temperature T∞ see Figure 1
Figure 1.
Problem set up.
We use Cartesian coordinates with two horizontal x- and y- axis located at the lower solid boundary and a positive z- axis is directed towards the free surface. The surface tension, τis assumed to be a linear function of the temperature
τ=τ0−γ(T−T0)E1
where τ0 is the value of τ at temperature T0 and the constant γ is positive for most fluids. The density of the fluid is given by
ρ=ρ0[1−α(T−T0)]E2
where αis the positive coefficient of the thermal liquid expansion and ρ0 is the value of ρ at the reference temperatureT0. Subject to the Boussinesq approximation, the governing equations for an incompressible, electrically conducting fluid in the presence of a magnetic field are expressed as follows:
Continuity equation:
∇⋅U=0E3
Momentum equation:
(∂∂t+U⋅∇)U=−1ρ∇Π+v∇2U+μ4πρ(H⋅∇)HE4
Energy equation :
(∂∂t+U⋅∇)T=κ∇2TE5
Magnetic field equations:
∇⋅H=0E6
(∂∂t+U.∇)H=(H.∇)U+η∇2HE7
where UTHρvκand ηdenote the velocity, temperature, magnetic field, pressure, density, kinematic viscosity, thermal diffusivity and electrical resistivity, respectively. Π=p+μ|H|2/8πis the magnetic pressure where p is the fluid pressure and μ is the magnetic permeability. When motion occurs, the upper free surface of the layer will be deformable with its position atz=d+f(xyt). At the free surface, we have the usual kinematic condition together with the conditions of continuity for the normal and tangential stresses. The temperature obeys the Newton's law of cooling, k∂T/∂n=h(T−T∞), where k and h are the thermal conductivity of the fluid and the heat transfer coefficient between the free surface and the air, respectively, and n is the outward unit normal to the free surface. The boundary conditions at the bottom wall, z = 0, are no-slip and conducting to the temperature perturbations.
To simplify the analysis, it is convenient to write the governing equations and the boundary conditions in a dimensionless form. In the dimensionless formulation, scales for length, velocity, time and temperature gradient are taken to be dκ/dd2/κandΔTrespectively. Furthermore, six dimensionless groups appearing in the problem are the Marangoni numberM=γΔTd/ρ0κv, the Biot number, Bi=hd/k, the Bond number, Bo=ρ0gd2/τ0, the Prandtl number,Pr=ν/κ , the Crispation number,Cr=ρ0νκ/τ0d and the internal heating,Q=qd2/2κΔT
Our control strategy basically applies a principle similar to that used by BauBau 1999, which is as follows:
Assumed that the sensors and actuators are continuously distributed and that each sensor directs an actuator installed directly beneath it at the same {x,y} location. The sensor detects the deviation of the free surface temperature from its conductive value. The actuator modifies the heated surface temperature according to the following rule BauBau 1999 :
T(xy0t)=1+BiBi−K(T(xy1t)−1Bi)E8
where K is the scalar controller gain. Equation 8 can be rewritten more conveniently as
T′(xy0t)=−K(T′(xy1t))E9
where T′ is the deviation of the fluid's temperature from its conductive value. The control strategy in equation9, in which K is a scalar will be used to demonstrate that our system can be controlled.
We study the linear stability of the basic state by seeking perturbed solutions for any quantity Φ(xyzt) in terms of normal modes in the form
Φ(xyzt)=Φ0(xyz)+ϕ(z)exp[i(αxx+αyy)+st]E10
where Φ0is the value of Φin the basic state, a=(αx2+αy2)1/2is the total horizontal wave number of the disturbance and s is a complex growth rate with a real part representing the growth rate of the instability and the imaginary part representing its frequency. At marginal stability, the growth rate s of perturbation is zero and the real part of s, ℜ(s)0represents unstable modes while ℜ(s)0 represents stable modes. Substituting equation 10 into equations 3 - 7 and neglecting terms of the second and higher orders in the perturbations we obtain the corresponding linearized equations involving only the z-dependent parts of the perturbations to the temperature and the z-components of the velocity denoted by T and w respectively,
[(D2−a2)2−H2D2]w=0E11
(D2−a2)T+w=0E12
subject to
sf−w(1)=0E13
P1Cr[(D2−3a2−H2)Dw(1)]−a2(a2+Bo)f=0E14
(D2+a2)w(1)+a2M(T(1)−(1+Q)f)=0E15
hz(1)=0E16
DT(1)+Bi(T(1)−(1+Q)f)=0E17
w(0)=0E18
Dw(0)=0E19
hz(0)=0E20
and
T(0)+KT(1)=0E21
on the lower rigid boundary z = 0. The operator D = d/dz denotes the differentiation with respect to the vertical coordinate z. The variables w, T and f denote respectively the vertical variation of the z-velocity, temperature and the magnitude of the free surface deflection of the linear perturbation to the basic state with total wave number a in the horizontal x-y plane and complex growth rates.
The effect of feedback control on the onset of Marangoni convection in a fluid layer with a magnetic field in the case of a deformable free surface (Cr≠0) is investigated numerically. The marginal stability curves in the (aM) plane are obtained numerically where M is a function of the parametersaBiBoCrand Q. For a given set of parameters, the critical Marangoni number for the onset of convection defined as the minimum of the global minima of marginal curve. We denote this critical value by Mc and the corresponding critical wave number,ac. The problem has been solved to obtain a detail description of the marginal stability curves for the onset of Marangoni convection when the free surface is perfectly insulated (Bi=0
Figure 2 shows the numerically calculated Marangoni number, M as a function of the wavenumber, a for different values of K in the caseCr=0. From Figure 4 it is seen that the critical Marangoni number increase with an increase of K. Thus, the magnetic always has a stabilizing effect on the flow. In the absence of controller gain, K = 0 and magnetic field, Q = 0, the present calculation reproduce closely the stability curve obtained by Pearson Pearson 1958. The present calculation are also reproduced the stability curve obtained by Wilson for K =0 and Q =100. It can been seen that the feedback control and magnetic field suppresses Marangoni convection. The critical Marangoni number, Mcincreases monotonically as the controller gain, K increases. In the case of non-deformable free surfaceCr=0, the controller can suppress the modes and maintain a no-motion state, but this situation is significantly different if the free surface is deformable,Cr≠0
When Crbecomes large the long-wavelength instability sets in as a primary one and the critical Marangoni numbers are at a = 0. Figure 3 shows the critical Marangoni number at the onset of convection as a function of the wave number, a, for a range of values of the controller gains, K when in the case of Cr=0.001andBo=0.1. At a = 0, the critical Marangoni number is zero and in this case, conductive state does not exist. Figure 3 shows that the controller is not effective at the wave number a = 0. Figure 4 shows the critical Marangoni number at the onset of convection as a function of the wave number, a for a range of values of the controller gains, K in the case Cr=0.0001 andBo=0.1. In this case, the global minimum occurs at a≠0 and as the controller gain K increases, the curve shifts upwards and most importantly, the controller increases the magnitude of the global minimum, thus it has a stabilizing effect.
Figure 2.
Numerically-calculated marginal stability curves for K = 0 and Q = 0 (solid line) and for various values of K (dashed line) in the case Q =100 andCr=0
Figure 3.
Numerically-calculated marginal stability curves for K = 0 and Q = 0 (solid line) and for various values of K (dashed line) in the case Q =100,Cr=0.001 and
Bo=0E22
Figure 4.
Numerically-calculated marginal stability curves for K = 0 and Q = 0 (solid line) and for various values of K (dashed line) in the case Q =100, Cr=0.0001andBo=0.1
The effect of the feedback control on the onset of steady Marangoni convection instabilities in a fluid layer with magnetic field has been studied. We have shown that the feedback control and magnetic field suppresses Marangoni convection We have also shown numerically that the effect of the controller gain and magnetic field is always to stabilize the layer in the case of a nondeforming surface. However, in the case of a deforming surface, the controller gain is effective depending on the parameter Cr and
Bo=0.1E24
References
1.ArifinN. M.NazarR.SenuN.2007Feedback control of the Marangoni-Bénard instability in a fluid layer with free-slip bottom.J. Phys. Soc. Japan, 76(1): 014401: 1- 4.
2.BauH. H.1999Contol of Marangoni-Bénard convection.Int. J.Heat Mass Transfer, 4213271341 .
3.ChandrasekharS.1961Hydrodynamic and Hydromagnetic Stability.
4.HashimI.ArifinN. M.2003Oscillatory Marangoni convection in a conducting fluid layer with a deformable free surface in the presence of a vertical magnetic field.Acta Mech, 164199-215.
5.HashimI.WilsonS. K.1999The effect of a uniform vertical magnetic field on the onset of oscillatory Marangoni convection in a horizontal layer of conducting fluid.Acta Mech. 132129- 146
6.HowleL. E.1997Linear stability analysis of controlled Rayleigh Bénard convection using shadowgraphic measurement.Physics of Fluids, 91131113113 .
7.NieldD. A.1966Surface tension and Buoyancy effects in cellular convection of an electrically conducting liquid in a magnetic field. Z. angew. Math. Mech, 17131-139.
8.OrA. C.KellyR. E.CortelezziL.SpeyerL.. (1999). Control of long wavelength Bénard-Marangoni convection.J. Fluid Mech., 387321341 .
9.PearsonJ. R. A.1958On convection cells induce by surface tension.J. Fluid Mech4489500 .
10.RayleighR.1916On convection currents in a horizontal layer of fluid, when the higher temperature is on the under side.Phil. Mag.32529546 .
11.Tang,J.BauH. H.1993Stabilization of the no-motion state in Rayleigh-Bénard convection through the use feedback control. Physical Review Letters, 7017951798 .
12.Tang,J.BauH. H.1994Stabilization of the no-motion state in Rayleigh-Bénard problem. Proceedings Royal Society, A447587607 .
13.WilsonS. K.1993The effect of a uniform magnetic field on the onset of steady Benard-Marangoni convection in a layer of conducting fluid.J. Engng Math, 27161-188.
14.WilsonS. K.1993The effect of a uniform magnetic field on the onset of Marangoni convection in a layer of conducting fluid.Q. Jl. Mech. Appl. Math, 46211-248.
15.WilsonS. K.1994 The effect of a uniform magnetic field on the onset of steady Marangoni convection in a layer of conducting fluid with a prescribed heat flux at its lower boundary.Phys. Fluid, 63591-3600.