A Linear Stability Analysis of Thermal Convection in a Fluid Layer with Simultaneous Rotation and Magnetic Field Acting in Different Directions

The onset of thermal convection of a Boussinesq fluid located in an unbounded layer heated from below and subject simultaneously to rotation andmagnetic 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 T (rate of rotation) and magnetic parameter Q (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.


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 g, inhibits the onset of convection and thus induces a stabilizing effect [1][2][3].The amount of the stabilizing effect depends on the rate of rotation (the Taylor number ).Results for the problem in which the vectors Ω and g 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 H that is parallel to the gravity force vector g [5].The extent of the inhibition depends on the value of the nondimensional magnetic parameter .It has been found that when H and g 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 g and H 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  1 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  1 direction, whereas in the latter case, by using a variational procedure, [1] shows that as the parameter   , defined as   =   tan   (where   is the wave number in the  1 direction), increases, keeping the Chandrasekhar number  Table 1: Linear problems under study.  is the angle between the rotation vector Ω and the gravity vector g.If   ̸ = 0, Ω lies in the  1  3 -plane.  is the angle between the magnetic field vector H and the gravity vector g.If   ̸ = 0, H lies in the  1  3 -plane.g acts along the  3 direction; see Figure 1 (defined by  =  2 ( cos   ) 2 D2 σ/  ], where , , D, σ,   , 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  = √ 2   +  2  (where   is the wave number in the  2 direction) at the onset of stationary convection shift to higher values.It has been shown that when H 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 H and g (as well as Ω and g) are not parallel, convection at marginal stability occurs as longitudinal rolls.On the other hand, when H is impressed with an inclination angle   = /2 to the direction of the vertical and H lies in the  1  3 -plane, the vertical component of the magnetic field ( cos   ) is equal to zero.The authors in [6] found that a magnetic field, impressed in a horizontal direction (  = /2), 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 (  = 0); hence, it was experimentally confirmed that the vertical component of the magnetic field ( cos   ) 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 (g, H, 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   =   tan   and   =   tan   , a critical situation may appear, at which the competition between H and Ω conduct to an unstable system; however as the values of   or   are modified, a cooperative work between H 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 g, rotation Ω, and magnetic field H 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.

Physical Situation
Figure 1 shows a schematic representation of the geometry under study.A layer of fluid of constant thickness D = 2 d parallel to the  1  2 -plane with no lateral boundaries is considered.The upper and lower boundaries are located at Figure 1: Thermal instability analysis in the presence of a rotation vector Ω and a magnetic field H.The gravity vector g acts in the negative  3 -direction.The temperatures  ℎ and   (where  ℎ >   ) are applied at the lower and upper boundaries, respectively.The distance between the two surfaces is D = 2 d.The angle between g and Ω is   , while the angle between g and H is   .Note that the vectors Ω and H lay in the  1  3 -plane.The surfaces are electrically nonconducting; thus, no currents can cross the boundary.
3 = ± d.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  3 -axis, in such a way that Ω lies in the  1  3plane.Gravity acts in the negative  3 -direction g = −i 3 .The fluid layer is assumed to be located in a magnetic field represented by the vector H, which can also be inclined at an angle   to the direction of the vertical  3 -axis, in such a way that H also lies in the  1  3 -plane.The surfaces are electrically nonconducting; thus, no currents can cross the boundary, and hence we must require  3 = 0, where  3 is the third component of the current density vector J.

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 g, the rotation vector Ω, and the magnetic field vector H may act in different directions.

Governing Equations.
Let us consider an infinite fluid layer of an electrically conducting fluid upon which a uniform magnetic field H 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, H is the impressed magnetic field, h is the perturbation of the magnetic field,  = (0, 0, 1) 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 k = V  i  = ∇ × h 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 ℎ 3 , V 3 ,  3 , and  3 are the  3 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  3 -axis.Additionally let the magnetic field vector H be also inclined at an angle   to the direction of the vertical  3 -axis.The  1 -axis of the Cartesian coordinate system is defined such as the vectors Ω and H lie in the  1  3 -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 H, respectively.Using (17) into ( 13)-( 16) we obtain )  3 . (18)

The Analysis into Normal Modes.
We analyse the perturbations , ℎ 3 , V 3 ,  3 , and  3 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  = ( 2  +  2  ) 1/2 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 D = 2 d.If we define the following nondimensional variables (20) is written as where   = D  .If instability sets in as stationary convection,  = 0, we can operate (25) by [4 2 −  2 ] and in the resulting equation we use (24) to get (let us assume that D = 2) 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 (4 2 −  2 ), 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 H = 0, 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  * 3 = −1 and at  * 3 = 1:  = 0,  = 0,  = 0,  = 0, (b) For free surfaces at  * 3 = −1 and at  * 3 = 1:  = 0,  = 0,  = 0,  2  = 0, (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 4 2 − 2 + 1/2 +(1/2) 1/2    = 0, 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.
(For this particular case in which Ω = 0, the  3 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).

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 ( * 3 ), ( * 3 ), and ( * 3 ) as well as their derivatives  = / * 3 ,  2 =  2 / * 2 3 ,  3 =  3 / * 3 3 , and  4 =  4 / * 4 3 , in the interval −1 ≤  * 3 ≤ 1 as follows: where  is an integer number defining the order of the derivative of the functions , , and , and   ( * 3 ) 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 [12][13][14]: where  0 = 2,   = 1 for  > 0. 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 Ω = 0; see (49), one equation for each  = 0, 1, 2, . . ., 2 + 2 − 6.The remaining six equations, those for 2 + 2 − 6 + 1 ≤  ≤ 2 + 2, 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 (Ω ̸ = 0) we generate, by substituting (52)-(57) into the three perturbation equations, one equation for each  = 0, 1, 2, . . ., 3 + 3 − 8.The remaining eight equations, those for 3 + 3 − 8 + 1 ≤  ≤ 3 + 3, 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 (62) 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 (Ax) = (Bx), where the vector x includes the coefficients   ,   (and   ),  is the Rayleigh number, and A and B are (2 + 2) × (2 + 2) matrices (or (3 + 3) × (3 + 3) matrices) whose first 2 + 2 − 6 (or 3 + 3 − 8) rows are defined by the two (or three) perturbation equations.The last six (or eight) rows of A are given by the boundary conditions, and the last six (or eight) rows of B 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 Ax = λBx, by using the  algorithm (see http://www.nag.co.uk/numeric/Fl/manual/ xhtml/F02/f02gjf.xml).

Thermal Instability Analysis in the Presence of a Rotation
Vector,   = 0, H = 0. 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.

Thermal Instability Analysis in the Presence of a Rotation
Vector,   ̸ = 0, H = 0. 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  = 100, as the value of the   parameter increases, the critical parameters (Rayleigh number   and wave number   ) become greater.Note in Table 4 that when   = 1,   = 3.15 and   = 1766.28;however, when   = 4,   = 3.35 and   = 1905.43.Also note in Table 5 that when   = 1,   = 2.66 and   = 855.0;however, when   = 4,   = 3.2 and   = 1145.92.Hence for a fixed Taylor number, the net effect of increasing   is to stabilize the system.

Thermal Instability Analysis in the Presence of a Magnetic
Field,   = 0, Ω = 0. 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 Table 10: Thermal instability analysis in the presence of a rotation and a magnetic field when both bounding surfaces are free: gravity vector g parallel to both rotation vector Ω and magnetic field vector H. Critical Rayleigh numbers and wave numbers at marginal stability for the onset of stationary convection.Taylor number  1 = / 4 = 10 5 for various values of the dimensionless magnetic parameter  1 = / 2 . and  are the Rayleigh number and the wave number, respectively, obtained by the Tau Chebyshev method. Ch and  Ch are the Rayleigh number and the wave number, respectively, obtained by Chandrasekhar [1].The superscripts 1 and 2 mean the two minima of the Rayleigh number.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.

Thermal Instability Analysis in the Presence of a Magnetic
Field,   ̸ = 0, Ω = 0. 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  = 100, 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   ̸ = 0, of the critical  number calculated by the Tau Chebyshev method occurs when   = 4.0 and  = 4.1, with an error of 3%, while the minimum percentage error occurs when   = 0.5 and  = 3.9, 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   = 1, the critical wave number  = 4; however, when   = 4, the critical wave number  = 4.1 Also note in Table 9 that when   = 1, the critical wave number  = 3.7; however, when   = 4, the critical wave number  = 4.1.Hence for a fixed Chandrasekhar number , the net effect of increasing   is to stabilize the system.

Thermal Instability Analysis in the Presence of a Rotation
Vector and a Magnetic Field,   = 0,   = 0. 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  1 = / 2 and for a fixed value of the modified Taylor number  1 = / 4 = 10 5 .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  1 = 10 5 ).We found (not shown in Table 10) that when  1 = 30,  = 3.93 and  = 1.125 × 10 6 , and when  1 = 200,  = 17.1 and  = 3.855 × 10 5 .Figure 2 shows that if  1 = 10 5 , when both bounding surfaces are free (see row (a)), and the magnetic parameter  1 is either  1 = 40 (curve labelled a),  1 = 100 (curve labelled b), or  1 = 200 (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  1 = 100 (curve labelled b) and  1 = 200 (curve labelled c); however, when  1 = 40 (curve labelled a), only one minimum value is present.Figure 2 shows on the left panel of both rows the case with pure rotation  1 = 10 5 , 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  1 = 200 (curve labelled c).This effect is clearly observed for low values of the modified wave number  =  2 / 2 .It is also observed that the system with  1 = 200 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  1 = 40 and the curve with circles), but for  1 = 100 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 ( 1 = 400 (see continuous lines),  1 = 1500 (see dashed lines),  1 = 2000 (see dotted lines), and  1 = 4000 (see dash-dotted lines)).Also shown on the right panel of both rows is the case with pure magnetic field  1 = 4000, that is, without rotation (see the curves with asterisks).It is pointed out that when the surfaces are free or rigid, when  1 = 400 (see the continuous lines), the system reaches its critical situation; that is, the minimum value of the Reynolds number is attained.However as  1 increases even more, the system becomes more stable.Now when the surfaces are free (row (a)), and when the with asterisks ( 1 = 4000, without rotation) is compared with the dashdotted line ( 1 = 4000,  1 = 10 5 ), 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 ( 1 = 4000, without rotation) is compared with the dash-dotted line ( 1 = 4000,  1 = 10 5 ), we observe that rotation has a significant stabilizing effect.According to the previous discussion we confirm the fact that when H 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  1 = 2000 and when the surfaces are free (left panel) and rigid (right panel).Notice that when the surfaces are free, the Coriolis force when  1 = 1000 (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  1 = 0).However when the surfaces are rigid (right panel), the Coriolis force when  1 = 1000 induces a stabilizing effect as compared with the case without rotation (see the asterisks curve for  1 = 0).It is clearly observed that as the modified Taylor number  1 acquires the values 5 × 10 5 and 10 6 , the system becomes more stable.

Thermal Instability Analysis in the Presence of a Rotation
Vector and a Magnetic Field,   ̸ = 0,   ̸ = 0. Figure 4 shows, for free surfaces, the dependence of the critical Rayleigh number  * = /100 and the critical wave number  =  2 / 2 , on the parameter   , when the Taylor number  1 = / 4 = 10 5 , the magnetic parameter  1 = / 2 = 40, and the gravity vector g is parallel to the magnetic field vector H (  = 0).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 0.5 ≤   ≤ 4 is observed.When   = 0 and   = 0, the critical Rayleigh number is  * = 8829.8;see left panel of Figure 2, while left panel of Figure 4 shows that when   = 0.5, the critical Rayleigh number has been lowered to  * = 3009.2.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 6 ≤   ≤ 10, 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   = 6 (see right panel).It is also possible to observe that the curves corresponding to   = 8 (curve labelled b) and   = 10 (curve labelled c) only show one minimum value.In general it is observed that as the   parameter increases, the wave number  =  2 / 2 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  1 = 10 5 , the magnetic parameter  1 = 40, and the gravity vector g is parallel to the rotation vector Ω (  = 0).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   = 0.5 resembles curve labelled a, of Figure 2 left panel, when   = 0 and   = 0. 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 0.5 ≤   ≤ 10; 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 15 ≤   ≤ 35, 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   = 15.We can also conclude that as the   parameter increases, the wave number  =  2 / 2 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  * = /100 on the wave number  =  2 / 2 when   is a fixed value equal to   = 6 and   has the values   = 2, 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   = 15 and   has the values   = 0.5, 2, 10, and 15.In both panels of Figure 6, the modified Taylor and magnetic numbers are  1 = 10 5 and  1 = 40, respectively.When   is fixed to   = 6 (see left panel), we observe that the curve   = 2 (continuous line) resembles the curve  (for   = 0) shown in right panel of Figure 4. However the case   = 2 is a bit more unstable than the case   = 0.When   increases to   = 10 (+ 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   = 20 ( * curve) and   = 40 (∘ curve).When   is fixed to   = 15 (see right panel), we observe that the curve   = 0.5 (continuous line) is very similar to the curve  (for   = 0) shown in Figure 5 right panel.When   increases to   = 2 (+ 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   = 10 ( * curve) and   = 15 (∘ 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.

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 H 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.

Table 2 :
[1]tical Rayleigh numbers and wave numbers at marginal stability for the onset of stationary convection when both bounding surfaces are rigid.Gravity vector g parallel to the rotation vector Ω.  is the Taylor number.is the Rayleigh number obtained by the Tau Chebyshev method.Ch is the Rayleigh number obtained by Chandrasekhar[1].

Table 3 :
[1]tical Rayleigh numbers and wave numbers at marginal stability for the onset of stationary convection when both bounding surfaces are free.Gravity vector g parallel to the rotation vector Ω.  is the Taylor number.is the Rayleigh number obtained by the Tau Chebyshev method.Ch is the Rayleigh number obtained by Chandrasekhar[1].

Table 4 :
Thermal instability analysis in the presence of a rotation vector when both bounding surfaces are rigid: gravity vector and rotation vector acting in different directions.Rayleigh numbers  for the Taylor number  = 100 and for various values of  and   .

Table 5 :
Thermal instability analysis in the presence of a rotation vector when both bounding surfaces are free: Gravity vector and rotation vector acting in different directions.Rayleigh numbers  for the Taylor number  = 100 and for various values of  and   .

Table 6 :
[1]tical Rayleigh numbers and wave numbers at marginal stability for the onset of stationary convection when both bounding surfaces are rigid.Gravity vector g parallel to the magnetic field vector H.  is the dimensionless magnetic parameter.isthe Rayleigh number obtained by the Tau Chebyshev method.Ch is the Rayleigh number obtained by Chandrasekhar[1].

Table 7 :
[1]tical Rayleigh numbers and wave numbers at marginal stability for the onset of stationary convection when both bounding surfaces are free.Gravity vector g parallel to the magnetic field vector H.  is the dimensionless magnetic parameter.isthe Rayleigh number obtained by the Tau Chebyshev method.Ch is the Rayleigh number obtained by Chandrasekhar[1].

Table 8 :
[1]rmal instability analysis in the presence of a magnetic field when both bounding surfaces are rigid: gravity vector and magnetic field vector acting in different directions.Rayleigh numbers for the dimensionless parameter  = 100 and for various values of  and   . is the Rayleigh number obtained by the Tau Chebyshev method.Ch is the Rayleigh number obtained by Chandrasekhar[1].

Table 9 :
[1]rmal instability analysis in the presence of a magnetic field when both bounding surfaces are free: gravity vector and magnetic field vector acting in different directions.Rayleigh numbers for the dimensionless parameter  = 100, and for various values of  and   . is the Rayleigh number obtained by the Tau Chebyshev method.Ch is the Rayleigh number obtained by Chandrasekhar[1].