Abstract

The onset of thermal convection of a Boussinesq fluid located in an unbounded layer heated from below and subject simultaneously to rotation and magnetic field, whose vectors act in different directions, is presented. To the knowledge of the authors, the convective thermal instability analysis for this complex problem has not been previously reported. In this paper, we use the Tau Chebyshev spectral method to calculate the value of the critical parameters (wave number and Rayleigh number at the onset of convection) as a function of (i) different kinds of boundaries, (ii) angle between the three vectors, and (iii) different values of the Taylor number (rate of rotation) and magnetic parameter (strength of the magnetic force). For the classical problems previously reported in the literature, we compare our calculations with Chandrasekhar’s variational method results and show that the present method is applicable.

1. Introduction

The thermal instability of a horizontal layer of fluid heated from below (Rayleigh-Bénard problem) is an important and classical fluid dynamics problem in science and engineering. The linear stability analysis of this problem, in terms of normal modes, has been carried out and reported in several investigations.

The rotating Rayleigh-Bénard problem has also been widely studied and it has been found that the effect of the Coriolis force due to rotation of the infinite layer of fluid, when the rotation vector acts about the vertical direction, that is, parallel to the gravity force vector , inhibits the onset of convection and thus induces a stabilizing effect [13]. The amount of the stabilizing effect depends on the rate of rotation (the Taylor number ). Results for the problem in which the vectors and act in different directions, with an angle between them, were not explicitly reported by [1, 4]; instead an analogy with the closely related problem of thermal convection subject to the action of a magnetic field was presented by [1].

Moreover, it has also been reported that, in the magnetic Rayleigh-Bénard problem, a stabilizing effect is obtained when the layer of an electrically conducting fluid is immersed in a uniform magnetic field vector that is parallel to the gravity force vector [5]. The extent of the inhibition depends on the value of the nondimensional magnetic parameter . It has been found that when and are parallel, convection at marginal stability is characterized by a cellular pattern; hence, longitudinal and transverse rolls appear simultaneously. Results for the problem in which the vectors and act in different directions, with an angle between them, and when instability sets in as stationary convection have been previously obtained by considering two cases: (i) solution of the perturbation equations which are independent of and (ii) solutions of the perturbation equations that conduct to a more general patterns of motion. In the former case, the onset of instability is characterized by rolls in the direction, whereas in the latter case, by using a variational procedure, [1] shows that as the parameter , defined as (where is the wave number in the direction), increases, keeping the Chandrasekhar number (defined by , where , , , , , and are the magnetic permeability, the magnitude of the magnetic field, the distance between the two plates, the coefficient of electrical conductivity, the density of reference, and the kinematic viscosity, resp.) as a constant value, both the critical Rayleigh number and the wave number (where is the wave number in the direction) at the onset of stationary convection shift to higher values. It has been shown that when is very slightly inclined to the direction of gravity, the extent to which transverse rolls are suppressed at marginal stability is also very slight. According to [1], when longitudinal rolls and transverse rolls appear simultaneously, a cellular pattern of convection emerges. Therefore if and (as well as and ) are not parallel, convection at marginal stability occurs as longitudinal rolls. On the other hand, when is impressed with an inclination angle to the direction of the vertical and lies in the -plane, the vertical component of the magnetic field () is equal to zero. The authors in [6] found that a magnetic field, impressed in a horizontal direction (), did not inhibit convection, even though the magnetic field was five times stronger than that one needed to suppress convection when acting in the vertical direction (); hence, it was experimentally confirmed that the vertical component of the magnetic field () is the critical parameter to inhibit the onset of convection as it is included in the Chandrasekhar number . According to the experimental and theoretical evidence, it is well known that rotation and magnetic field have a stabilizing effect on the Rayleigh-Bénard problem when they act separately [2, 4, 5, 7].

When rotation and magnetic field act simultaneously and when both of them are parallel to the gravity vector, the stabilizing effects have conflicting tendencies; that is, the results reveal some very unexpected features showing the complex behaviour of the flow; thus, it has been found that the total effect is not always stabilizing with respect to both fields [1, 4, 8, 9]. The conflictual behaviour depends on the values of the parameters of the system: rotation rate (Taylor number ) and magnetic parameter (Chandrasekhar number ). Reference [4] only shows explicit results on the dependence of the critical Rayleigh number on and , for the problem in which the conducting liquid is located between two free boundaries, and when the medium adjoining the fluid is nonconducting, because in this situation the system of equations can be solved explicitly.

To the knowledge of the authors, the solution of the problem, in which the three implicated vectors (, , and ) are not parallel, has not been previously reported. In this paper we show, for the rotating-magnetic Rayleigh-Bénard problem with rigid and free surfaces and nonconducting boundaries, the influence on the critical Rayleigh number of the cooperative work (or competition), between the strength of the magnetic field and the Coriolis force, when they act simultaneously, but in different directions relatively with the gravity vector. In the literature, very few investigations have been carried out aimed to study the cooperative work of perturbations in the flow to inhibit the thermal convection of an infinite layer of fluid; for example, [7, 10, 11] have found that when the layer of fluid is simultaneously subject to rotation and to a salt concentration field, the stabilizing effects are cumulative; that is, rotation and salt concentration (salted from below) show a cooperative behaviour.

By using the Tau Chebyshev method, we find that depending on the value of the parameters and , a critical situation may appear, at which the competition between and conduct to an unstable system; however as the values of or are modified, a cooperative work between and may appear which leads to a system in which the thermal convection of the layer of fluid is inhibited.

The paper is organized as follows. In Section 2 we briefly describe the physical model of the problem under study. In Section 3 we formulate the general perturbation equations, from which we get the solution of the problems 1 to 6 described in Table 1. Note that the classical linear problems 1 to 5 of this table have been previously solved by [1], while the linear problem 6 of Table 1 is the subject of the present paper that is, thermal instability analysis of a fluid layer subject to gravity , rotation , and magnetic field vectors, with angles and different from zero. In Section 4 we describe the Tau Chebyshev spectral method that is used to solve the set of perturbation equations. Section 5 shows the solution of the six formulated linear problems by using the Tau Chebyshev method, together with a comparison of our calculations with the solution of the classical problems previously published in the literature; see problems 1 to 5 of Table 1. Finally, in Section 6 we present our concluding remarks.

2. Physical Situation

Figure 1 shows a schematic representation of the geometry under study. A layer of fluid of constant thickness parallel to the -plane with no lateral boundaries is considered. The upper and lower boundaries are located at . The temperatures and (where ) are applied at the lower and upper boundaries, respectively. The whole layer can rotate with an angular velocity vector , which can be inclined at an angle to the direction of the vertical -axis, in such a way that lies in the -plane. Gravity acts in the negative -direction . The fluid layer is assumed to be located in a magnetic field represented by the vector , which can also be inclined at an angle to the direction of the vertical -axis, in such a way that also lies in the -plane. The surfaces are electrically nonconducting; thus, no currents can cross the boundary, and hence we must require , where is the third component of the current density vector .

3. Thermal Instability Analysis in the Presence of a Rotation Vector and a Magnetic Field

In this section we analyse the linear thermal instability of a rotating infinite layer of an incompressible, electrically conducting fluid upon which a uniform magnetic field is impressed. The gravity vector , the rotation vector , and the magnetic field vector may act in different directions.

3.1. Governing Equations

Let us consider an infinite fluid layer of an electrically conducting fluid upon which a uniform magnetic field is impressed. The fluid is also subjected to rotation about a fixed axis. For an incompressible conducting fluid, the linear convective thermal stability of a fluid layer (initially at rest and with a steady adverse temperature gradient) is analyzed by solving the following perturbation equations in the Boussinesq approximation and in a rotating frame of [1, page 199]: where is the time, are the components of the velocity vector, denotes the kinematic viscosity, are the components of the rotation vector, is the coefficient of volume expansion, is the thermal diffusivity, is the perturbation in the temperature, is the acceleration of gravity, is the magnetic permeability, denotes the resistivity, is the impressed magnetic field, is the perturbation of the magnetic field, is the unit vector in the direction of the vertical, and is the imposed adverse temperature gradient defined as where and are the temperature of the lower hot surface and the temperature of the upper cold surface, respectively. The last two terms on the right hand side of (1) are the Coriolis force and the Lorentz force, respectively, and the term is defined as where is the centrifugal force, is the perturbation in the pressure, and is the density at a mean temperature .

Taking the curl of (1) we may write where are the components of the vorticity vector and the vector is the current density induced by the perturbation. Taking the curl of (9) once again, we obtain Taking the curl of (4) we have Multiplying (2), (4), (9), (10), and (11) by the components we get where , , , and are the components of the perturbation of the magnetic field, the current density, the vorticity, and the velocity, respectively. In this investigation the set of equations to be solved are (12)–(16). Let the rotation vector be inclined at an angle to the direction of the vertical -axis. Additionally let the magnetic field vector be also inclined at an angle to the direction of the vertical -axis. The -axis of the Cartesian coordinate system is defined such as the vectors and lie in the -plane; see Figure 1. Thus the unit vector in the direction of the vertical, the rotation vector, and the magnetic field vector are written as where and are the magnitude of the rotation vector and the magnitude of the magnetic field vector , respectively. Using (17) into (13)–(16) we obtain

3.1.1. The Analysis into Normal Modes

We analyse the perturbations , , , , and into two dimensional periodic waves characterized by a particular wave number . We assume the forms of the perturbations as If we define the wave number of the disturbance as and as a constant, (12) and (18) can be written as As it is shown in Figure 1, the origin of the coordinate system is located at the middle between the two surfaces; hence, the distance between the two surfaces is . If we define the following nondimensional variables (20) is written as where . If instability sets in as stationary convection, , we can operate (25) by and in the resulting equation we use (24) to get (let us assume that ) where is the Taylor number, If we define the resistivity as it is possible to write (31) as where the variable defined by Chandrasekhar [1] is written as

We can eliminate from (26) by using (23), then we obtain We operate (35) by , and in the resulting equation we use (22); hence, we obtain where the Rayleigh number is defined as From (36) we can define the variable as hence, (36) is written as By using the following relationships into (38) we can formulate the relevant equations to be solved as Notice that when the magnetic field is negligible , that is, pure rotation, and the relevant equations to be solved are the following: together with the following boundary conditions. (a)For rigid surfaces at   and at :  (b)For free surfaces at   and at : (It must be pointed out that due to the process we follow to derive the general relevant equations, the vorticity equation (46) has both different form and higher order than the one we obtained by following Chandrasekhar’s procedure; see Chandrasekhar [1], Chapter III, page 90. The equation for the vorticity that we obtained by considering only rotation without magnetic field is , which is the same as the one obtained by Chandrasekhar [1].)

On the other hand when the rotation vector is equal to zero, the relevant equations are together with the following boundary conditions. (a)For rigid surfaces at   and at :  (b)For free surfaces at   and at : (For this particular case in which , the component of vorticity vanishes, see Chandrasekhar [1], Chapter IV, page 166; hence the general equation (43) is not taken into account in the analysis.)

Therefore the boundary conditions needed to solve the general equations (41)–(43) are defined by (47) and (48).

4. The Tau Chebyshev Numerical Method

In order to find an approximate solution to the set of perturbations equations formulated in the previous section, together with the boundary conditions, we use expansions in Chebyshev polynomials. Therefore we expand the variables , , and as well as their derivatives , , , and , in the interval as follows: where is an integer number defining the order of the derivative of the functions , , and , and are the th-degree Chebyshev polynomials of the first kind (with as the maximum degree used for the polynomials). The Chebyshev polynomials satisfy the orthogonality relation [1214]: where , for . Equations for the expansion coefficients , , , (in terms of ), (in terms of ), and (in terms of ) are generated by substituting (52)–(57) into the set of perturbation equations shown in the previous section and applying the orthogonality property (58). In the literature it is possible to find explicit formulae for the expansion coefficients in terms of , in terms of , and in terms of [12, 15].

The method that we have selected to obtain the equations for the coefficients , , and is the Tau Chebyshev method. The Tau Chebyshev method is an efficient numerical procedure for finding eigenvalues and eigenfunctions of systems in fluid and solid mechanics theories. One advantage of the algorithm is that it is able to calculate as many eigenvalues as required, not just the dominant one [13, 14]. Using this numerical methodology we generate, by substituting (52), (53), (55), and (56) into the couple of the general perturbation equations for the case when ; see (49), one equation for each . The remaining six equations, those for , are generated by using the six boundary conditions. Hence for rigid surfaces, the set of boundary conditions are written as whereas for free surfaces the set of boundary conditions are written as Using the Tau Chebyshev method in the cases with rotation () we generate, by substituting (52)–(57) into the three perturbation equations, one equation for each . The remaining eight equations, those for , are generated by using the eight boundary conditions. Hence for rigid surfaces, the set of boundary conditions are written as whereas for free surfaces the set of boundary conditions are written as Notice that the high frequency behaviour (i.e., high ) of the solution is not governed by the perturbation equations but by the boundary conditions [12]. The Tau equations can be written as a generalized eigenvalue problem , where the vector includes the coefficients , (and ), is the Rayleigh number, and and are matrices (or matrices) whose first (or ) rows are defined by the two (or three) perturbation equations. The last six (or eight) rows of are given by the boundary conditions, and the last six (or eight) rows of vanish [16]. The system of equations has been solved numerically by using the routine F02GJF from the NAG library. The routine F02GJF calculates all the eigenvalues of the eigenproblem , by using the algorithm (see http://www.nag.co.uk/numeric/Fl/manual/xhtml/F02/f02gjf.xml).

5. Results

5.1. Thermal Instability Analysis in the Presence of a Rotation Vector, ,

Tables 2 and 3 show the critical Rayleigh number and wave number of the unstable modes at marginal stability for the onset of stationary convection when both bounding surfaces are rigid and when both bounding surfaces are free, respectively. It is possible to observe that the results obtained by the Tau Chebyshev method are in full agreement with the variational formulation; see [1]. It is also observed the inhibition effect of the Coriolis force as the Taylor number increases. As becomes greater, the value of the critical wave number also increases.

5.2. Thermal Instability Analysis in the Presence of a Rotation Vector, ,

Tables 4 and 5 show the results obtained by the Tau Chebyshev method when both bounding surfaces are rigid and when both bounding surfaces are free, respectively. It is observed that when the Taylor number is , as the value of the parameter increases, the critical parameters (Rayleigh number and wave number ) become greater. Note in Table 4 that when , and ; however, when , and . Also note in Table 5 that when , and ; however, when , and . Hence for a fixed Taylor number, the net effect of increasing is to stabilize the system.

5.3. Thermal Instability Analysis in the Presence of a Magnetic Field, ,

Tables 6 and 7 show the critical Rayleigh number and wave number of the unstable modes at marginal stability for the onset of stationary convection when both bounding surfaces are rigid and when both bounding surfaces are free, respectively. It is possible to observe that the results obtained by the Tau Chebyshev method are in full agreement with the variational formulation; see [1]. It is also observed the inhibition effect of the magnetic field, as the magnetic number increases. As becomes greater, the value of the critical wave number also increases.

5.4. Thermal Instability Analysis in the Presence of a Magnetic Field, ,

Tables 8 and 9 show the results obtained by the Tau Chebyshev method when both bounding surfaces are rigid and when both bounding surfaces are free, respectively. It is observed that when the Chandrasekhar number is fixed to , as the value of the parameter increases, the critical Rayleigh number becomes greater. In Table 8 we present a comparison between the results obtained by [1] and the results of the Tau Chebyshev method. We observe that if the data obtained by Chandrasekhar are the exact values, the maximum percentage error (% error), when , of the critical number calculated by the Tau Chebyshev method occurs when and , with an error of 3%, while the minimum percentage error occurs when and , with an error of 0.3%. Thus we may conclude that our calculations are in good agreement with Chandrasekhar’s results. Note in Table 8 that when , the critical wave number ; however, when , the critical wave number Also note in Table 9 that when , the critical wave number ; however, when , the critical wave number . Hence for a fixed Chandrasekhar number , the net effect of increasing is to stabilize the system.

5.5. Thermal Instability Analysis in the Presence of a Rotation Vector and a Magnetic Field, ,

Table 10 shows a comparison between the results obtained by the Tau Chebyshev method and those of Chandrasekhar’s variational formulation when both bounding surfaces are free. The results have been obtained for several values of the modified magnetic number and for a fixed value of the modified Taylor number . Notice that the agreement is encouraging. Also note that the two minima of the Rayleigh number are successfully reproduced by the Tau Chebyshev method. It is interesting to mention that by using the Tau Chebyshev method we found two additional critical Rayleigh numbers that were not reported by Chandrasekhar (see [1], Table XVIII, page 205, for ). We found (not shown in Table 10) that when , and , and when , and . Figure 2 shows that if , when both bounding surfaces are free (see row (a)), and the magnetic parameter is either (curve labelled a), (curve labelled b), or (curve labelled c), see left panel, the two minima of the Rayleigh number curves reported by Chandrasekhar are recovered (see [1], Figure 49, page 208). Row (b), left panel of Figure 2, shows that when the bounding surfaces are rigid, the two minima only slightly appear when (curve labelled b) and (curve labelled c); however, when (curve labelled a), only one minimum value is present. Figure 2 shows on the left panel of both rows the case with pure rotation , that is, without magnetic field (see the curves with circles). It is pointed out that when the surfaces are free (row (a), left panel), the effect of the magnetic field is to destabilize the system, particularly when (curve labelled c). This effect is clearly observed for low values of the modified wave number . It is also observed that the system with is more unstable than the pure rotation case. However when the surfaces are rigid (see row (b), left panel), it is possible to observe that as the magnetic field becomes greater, it firstly induces a stabilizing effect, with respect to the pure rotation case (see curve labelled a for and the curve with circles), but for and 200, it initiates a destabilizing effect; however, the pure rotation case (see the curve with circles) is more unstable. On the other hand Figure 2 shows, on the right panel of both rows, the influence of greater magnetic numbers ( (see continuous lines), (see dashed lines), (see dotted lines), and (see dash-dotted lines)). Also shown on the right panel of both rows is the case with pure magnetic field , that is, without rotation (see the curves with asterisks). It is pointed out that when the surfaces are free or rigid, when (see the continuous lines), the system reaches its critical situation; that is, the minimum value of the Reynolds number is attained. However as increases even more, the system becomes more stable. Now when the surfaces are free (row (a)), and when the curve with asterisks (, without rotation) is compared with the dash-dotted line (, ), it is possible to observe that the Coriolis force induces a marginal stabilizing effect. However when the surfaces are rigid (row (b)) and when the curve with asterisks (, without rotation) is compared with the dash-dotted line (, ), we observe that rotation has a significant stabilizing effect. According to the previous discussion we confirm the fact that when and act simultaneously and when both are parallel to the gravity vector, the stabilizing effects have conflicting tendencies. Figure 3 shows the influence of the Coriolis force when the magnetic number is fixed to and when the surfaces are free (left panel) and rigid (right panel). Notice that when the surfaces are free, the Coriolis force when (curve labelled a) does not produce any effect on the system as compared to that when there is no rotation (see the asterisks curve for ). However when the surfaces are rigid (right panel), the Coriolis force when induces a stabilizing effect as compared with the case without rotation (see the asterisks curve for ). It is clearly observed that as the modified Taylor number acquires the values and , the system becomes more stable.

5.6. Thermal Instability Analysis in the Presence of a Rotation Vector and a Magnetic Field, ,

Figure 4 shows, for free surfaces, the dependence of the critical Rayleigh number and the critical wave number , on the parameter , when the Taylor number , the magnetic parameter , and the gravity vector is parallel to the magnetic field vector (). On the left panel the competition between the stabilizing effect of the Coriolis force and the strong destabilizing properties of the parameter when it is in the range is observed. When and , the critical Rayleigh number is ; see left panel of Figure 2, while left panel of Figure 4 shows that when , the critical Rayleigh number has been lowered to . Note, on the left panel, the presence of the two minima in the curves of the Rayleigh number . On the right panel of Figure 4 it is shown that when the parameter is in the range , the fluid layer becomes more stable as increases. Therefore, we have found that there is a critical value of the parameter at which the destabilizing behaviour is changed to a stabilizing performance, and the critical value occurs around (see right panel). It is also possible to observe that the curves corresponding to (curve labelled b) and (curve labelled c) only show one minimum value. In general it is observed that as the parameter increases, the wave number also increases and a tendency to get only one minimum is apparent.

Figure 5 shows, for free surfaces, the dependence of the critical Rayleigh number and the critical wave number on the parameter , when the Taylor number , the magnetic parameter , and the gravity vector is parallel to the rotation vector (). On the left panel it is shown that the influence of the parameter , when its value is 0.5 (curve labelled a), is almost negligible; hence, curve for resembles curve labelled a, of Figure 2 left panel, when and . On the left panel of Figure 5, it is possible to observe the competition between the stabilizing effect of the Coriolis force and the destabilizing influence of the parameter when its values are in the range ; that is, as the parameter increases, the fluid layer gets more unstable. Right panel of Figure 5 shows the cooperative behaviour between the Coriolis force and the parameter; that is, as the parameter increases even more in the range , the fluid layer becomes more stable. Hence we have found again that a critical value of the parameter exists at which the destabilizing performance of is transformed to a stabilizing behaviour. The critical is around . We can also conclude that as the parameter increases, the wave number also increases and the tendency to get only one minimum is clear.

Finally left panel of Figure 6 shows for free surfaces, the dependence of the Rayleigh number on the wave number when is a fixed value equal to and has the values , 10, 20, and 40. Right panel of Figure 6 also shows for free surfaces the dependence of the Rayleigh number on the wave number when is a fixed value equal to and has the values , 2, 10, and 15. In both panels of Figure 6, the modified Taylor and magnetic numbers are and , respectively. When is fixed to (see left panel), we observe that the curve (continuous line) resembles the curve (for ) shown in right panel of Figure 4. However the case is a bit more unstable than the case . When increases to (+ curve), the system becomes even more unstable; however, it is clearly observed that there exists a critical value of (around 10), at which the system shifts its unstable behaviour to a more stable performance, see the curves for ( curve) and ( curve). When is fixed to (see right panel), we observe that the curve (continuous line) is very similar to the curve (for ) shown in Figure 5 right panel. When increases to (+ curve), the system becomes more unstable, however it is clearly observed that there also exists a critical value of (around 2), at which the system shifts its unstable behaviour to a more stable performance; see the curves for ( curve) and ( curve). In both panels of Figure 6 it is observed that before the critical situation there exists a competition between and parameters, and the result is that the fluid layer is transformed into a more unstable system. However after the critical situation, there is a cooperative work of and to inhibit the thermal convection of the infinite layer of fluid.

6. Conclusions

The Tau Chebyshev method has been used to solve the classical problem of thermal convection on an infinite fluid layer with rotation and an impressed magnetic field. Our methodology has been verified with results previously published in the literature, through the solution of the first five cases shown in Table 1. By using the Tau Chebyshev method, we find that, for the sixth case shown in Table 1, a critical value of the and parameters exists at which the competition between and , which leads to an unstable system, can be transformed to a cooperative work that promotes a system in which the thermal convection of the layer of fluid is inhibited.

Acknowledgments

Financial support for this investigation was provided by the Center for Aerospace Research & Education, University of California, Irvine, CA, USA. R. Avila thanks Professor. Dr. Satya Atluri for his encouragement.