InTech uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Technology » "Advanced Technologies", book edited by Kankesu Jayanthakumaran, ISBN 978-953-307-009-4, Published: October 1, 2009 under CC BY-NC-SA 3.0 license. © The Author(s).

Feedback Control of Marangoni Convection with Magnetic Field

By Norihan Arifin, Haliza Rosali and Norfifah Bachok
DOI: 10.5772/8232

Article top

Overview

Figure 1. Problem set up.

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 and C r = 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, C r = 0.001 and

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, C r = 0.0001 and B o = 0.1

Feedback Control of Marangoni Convection with Magnetic Field

Norihan Md. Arifin1, Haliza Rosali1 and Norfifah Bachok1

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.

2. Problem Formulation

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) (1)

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)] (2)

where α is the positive coefficient of the thermal liquid expansion and ρ0 is the value of ρ at the reference temperature T0 . 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=0 (3)

Momentum equation:

 (∂∂t+U⋅∇​)U=−1ρ∇Π+v∇2U+μ4πρ(H⋅∇)H (4)

Energy equation :

 (∂∂t+U⋅∇)T=κ∇2T (5)

Magnetic field equations:

 ∇⋅H=0 (6)
 (∂∂t+U.∇)H=(H.∇)U+η∇2H (7)

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 at z=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, kT/n=h(TT) , 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ΔT respectively. Furthermore, six dimensionless groups appearing in the problem are the Marangoni number M=γΔ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) (8)

where K is the scalar controller gain. Equation 8 can be rewritten more conveniently as

 T′(xy0t)=−K(T′(xy1t)) (9)

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.

3. Linearized Problem

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] (10)

where Φ0 is the value of Φ in the basic state, a=(αx2+αy2)1/2 is 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)0 represents 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=0 (11)
 (D2−a2)T+w=0 (12)

subject to

 sf−w(1)=0 (13)
 P1Cr[(D2−3a2−H2)Dw(1)]−a2(a2+Bo)f=0 (14)
 (D2+a2)w(1)+a2M(T(1)−(1+Q)f)=0 (15)
 hz(1)=0 (16)
 DT(1)+Bi(T(1)−(1+Q)f)=0 (17)
 w(0)=0 (18)
 Dw(0)=0 (19)
 hz(0)=0 (20)
and
 T(0)+KT(1)=0 (21)

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.

4. Results and disussion

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 (Cr0) is investigated numerically. The marginal stability curves in the (aM) plane are obtained numerically where M is a function of the parameters a Bi Bo Cr and 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 case Cr=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, Mc increases monotonically as the controller gain, K increases. In the case of non-deformable free surface Cr=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, Cr0

When Cr becomes 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.001 and Bo=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 and Bo=0.1 . In this case, the global minimum occurs at a0 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 and Cr=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=0 (22)

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.0001 and

 Bo=0.1 (23)

5. Conclusion

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.1 (24)

References

1 - N. M. Arifin, R. Nazar, N. Senu, 2007 Feedback 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 - H. H. Bau, 1999 Contol of Marangoni-Bénard convection. Int. J.Heat Mass Transfer, 42 1327 1341 .
3 - S. Chandrasekhar, 1961 Hydrodynamic and Hydromagnetic Stability.
4 - I. Hashim, N. M. Arifin, 2003 Oscillatory Marangoni convection in a conducting fluid layer with a deformable free surface in the presence of a vertical magnetic field. Acta Mech, 164 199-215.
5 - I. Hashim, S. K. Wilson, 1999 The effect of a uniform vertical magnetic field on the onset of oscillatory Marangoni convection in a horizontal layer of conducting fluid. Acta Mech. 132 129- 146
6 - L. E. Howle, 1997 Linear stability analysis of controlled Rayleigh Bénard convection using shadowgraphic measurement. Physics of Fluids, 9 11 3111 3113 .
7 - D. A. Nield, 1966 Surface tension and Buoyancy effects in cellular convection of an electrically conducting liquid in a magnetic field. Z. angew. Math. Mech, 17 131-139.
8 - A. C. Or, R. E. Kelly, L. Cortelezzi, L. Speyer, . (1999). Control of long wavelength Bénard-Marangoni convection. J. Fluid Mech., 387 321 341 .
9 - J. R. A. Pearson, 1958 On convection cells induce by surface tension. J. Fluid Mech 4 489 500 .
10 - R. Rayleigh, 1916 On convection currents in a horizontal layer of fluid, when the higher temperature is on the under side. Phil. Mag. 32 529 546 .
11 - J. Tang,, H. H. Bau, 1993 Stabilization of the no-motion state in Rayleigh-Bénard convection through the use feedback control. Physical Review Letters, 70 1795 1798 .
12 - J. Tang,, H. H. Bau, 1994 Stabilization of the no-motion state in Rayleigh-Bénard problem. Proceedings Royal Society, A 447 587 607 .
13 - S. K. Wilson, 1993 The effect of a uniform magnetic field on the onset of steady Benard-Marangoni convection in a layer of conducting fluid. J. Engng Math, 27 161-188.
14 - Wilson S. K. 1993 The effect of a uniform magnetic field on the onset of Marangoni convection in a layer of conducting fluid. Q. Jl. Mech. Appl. Math, 46 211-248.
15 - S. K. Wilson, 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, 6 3591-3600.