Next Article in Journal
Fault Location of Distribution Network Based on Back Propagation Neural Network Optimization Algorithm
Next Article in Special Issue
Direct Numerical Simulation of Bubble Cluster Collapse: Shape Evolution and Energy Transfer Mechanisms
Previous Article in Journal
Quantitative Characteristics of Micro Bedding Fractures in the Wufeng–Longmaxi Formation Based on High-Resolution Map Imaging Technology
Previous Article in Special Issue
Droplet Based Estimation of Viscosity of Water–PVP Solutions Using Convolutional Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Analysis of Cross Rolls of Electrically Conducting Fluid under an Applied Magnetic Field with Rotation

1
Department of Mathematics, University College of Science, Osmania University, Hyderabad 500007, India
2
Instituto de Alta Investigación, Universidad de Tarapacá, Casilla 7 D, Arica 1000000, Chile
3
Department of Mathematics, Shivagangotri, Davangere University, Davangere 577007, India
*
Author to whom correspondence should be addressed.
Processes 2023, 11(7), 1945; https://doi.org/10.3390/pr11071945
Submission received: 18 April 2023 / Revised: 13 June 2023 / Accepted: 22 June 2023 / Published: 27 June 2023
(This article belongs to the Special Issue Complex Fluid Dynamics Modeling and Simulation)

Abstract

:
The proposed planer layer dynamo physical model has real-world applications, especially in the Earth’s liquid core. Thus, in this paper, an attempt is made to understand the finite amplitude convection when there exists a coupling between the Lorentz force and the Coriolis force. In particular, the effect of a horizontally applied magnetic field is studied on the Rayleigh–Bénard convection (RBC) that contains the electrically conducting fluid and rotates about its vertical axis. Free–free boundary conditions are assumed on the geometry. Attention is focused on the nonlinear convective flow behavior during the occurrence of cross rolls which are perpendicular to the applied magnetic field and parallel to the rotation axis. The visualization of cross rolls is achieved using the Fourier analysis of perturbations up to the O( ε 8 ). The relationship of the Nusselt number ( N u ) with respect to the Rayleigh number (R), the Ekman number (E), and the Elsasser number ( Λ ) is investigated. It is observed that E generates a strong damping effect on the flow velocity and on the heat transfer at high rotation rates. Using the heatline concept, it is observed that the temperature within the central regime is enhanced as the Λ increases. The results show that either E decreases or Λ increases, then the heat transfer rate increases.

1. Introduction

In astrophysical and geophysical models related to the Sun, stars, and the outer core of Earth, the convection is affected by both the Coriolis and Lorenz forces. Such a rotating magnetoconvection model with Boussinesq approximation has been studied by many authors, for example, Chandrasekhar [1], Roberts [2], and Cox and Matthews [3], etc. The linear stability analysis of this model shows the system is unstable with respect to either stationary convection or oscillatory convection and depends on the governing physical parameters, namely, Rayleigh number (R), Chandrasekhar number (Q), Taylor number ( T a ), thermal Prandtl number ( P r ), and magnetic Prandtl number ( P m ) or Roberts number. Most of the experimental studies on Rayleigh–Bénard convection ( R B C ), rotating Rayleigh–Bénard convection, magnetoconvection, and rotating magnetoconvection ( R M C ) have focused on heat transfer laws. When the control parameter R exceeds its critical value, a cellular regime of steady convection starts to appear. In addition, the motion increases its intensity but remains laminar and steady for a large range of values of R, followed by unsteady turbulent convection. Finite amplitude cellular convection with better approximate solutions has been studied for R B C by Malkus and Veronis [4] and Kuo [5].
In R B C with T a 0 and Q 0 , the critical R ( = R c ) remains approximately a constant until Q reaches a certain value. When Q increases further, R c starts to decrease, reaches a minimum and again starts to increase (Chandrasekhar, [1]). Braginsky [6] also studied these rotating magnetic systems and stressed the importance of the Archimedean, magnetic, and Coriolis forces. He revealed that these forces, together with the pressure gradient, would determine the dynamic balance and inertial forces. Eltayeb [7] considered the linear stability analysis to analyze the convection in the hydromagnetic rotating layer. He observed when the principle of exchange of stabilities is valid, four different models can be classified based on the relative directions of the constant magnetic field, B 0 and Ω , ( B 0 , Ω ) namely, (vertical, vertical); (horizontal, vertical); (horizontal, horizontal); (vertical, horizontal) and also for different types of boundaries. The numerical results indicated that the asymptotic dependence of R c on T a and Q are equal and independent of the nature of the boundary conditions considered. Later, Eltayeb [8] extended his previous model [7] to study the convective motions near the onset of oscillatory convection. He classified three different motions near the onset, namely: (i) T a > Q , the results for the rotating non-magnetic case, which are retained to leading order; (ii) Q > T a , the results are similar to those for the magnetic non-rotating case to leading order; and (iii) T a and Q are of same order, the minimum temperature gradient required for the instability is greatly reduced. When these solutions were examined in detail, it was observed that the motions that follow the onset of instability depended primarily on the electrical conductivity rather than on the kinematic properties of boundaries. In addition, in the leading order, the boundary conditions to be applied to the mainstream solutions depended on the conductivity of the boundary but not on the no-slip conditions. When either of the magnetic field or rotation were dominant, there was a possibility of the occurrence of two-dimensional motion. In this case, the Taylor–Proudman theorem was satisfied, while when both magnetic field and rotation were influential, this theorem was no longer valid and the motions were essentially of three-dimensional in nature.
The important laboratory experiments on R B C , rotating R B C and magnetoconvection R B C using the liquid gallium ( P r = 0.025 ) as the working fluid have been carried out by Aurnou and Olson [9]. The properties of liquid gallium are similar to those of the liquid iron in the Earth core. The studies of magnetoconvection have vast industrial applications too [10,11]. The R c for the magnetoconvection is experimentally determined as a function of Q and T a . At low rotation rates, the R c increases linearly with magnetic field intensity. At moderate rotation rates, coherent thermal oscillations were detected by Aurnou and Olson [9] near the onset of convection. These experimental results at the onset were compared with the theoretical predictions of Chandrasekhar [1]. In nearly all of the experimental results, it was mentioned that no well-defined steady convective regime was found. Instead, unsteady or turbulent convection was detected just after onset. Later, these experimental predictions were reproduced by using direct numerical simulations near the onset (Rani et al. [12]). These simulations showed the occurrence of interesting cell patterns. The most relevant geodynamo models were given by Roberts and Jones [13], who extended the model of Braginsky [14] with two sets of boundary conditions.
The present nonlinear convection problem is studied based on the plane layer model proposed by Roberts and Jones [13]. In astrophysics and even in planetary physics, the model considered by Roberts and Jones [13] is yet sufficiently close to reality and is, really, heuristic. This model is very convenient for laboratory experiments, too, which have not yet been done. The linear planer layer dynamo model was considered by Roberts and Jones [13] with the limiting case of Prandtl number tending to infinity. This limiting case enabled the removal of the inertial terms and thus the resultant linearized equation of motion filtered the fast-inertial modes and Alfven waves. The main reason for considering this limiting case was that it simplified the analysis considerably and could able to exhibit an amazingly rich structure. The other motivation was that it is an important limit for geodynamo modelling, in which fast modes are believed to be relatively unimportant.
In the linear stability analysis of the planer layer problem [13], the effect of physical parameters such as E, Λ , R, and q have been thoroughly studied for the occurrence of parallel rolls, cross rolls, and oblique rolls at the onset of convection. In general, E and q 1 represent the geophysical models but with these values, it is very difficult to simulate these models and obtain the converged results. This difficulty can be overcome by using one of the approximate methods, such as weakly nonlinear analysis, which yields comparable results with the experimental observations. This method of weakly nonlinear analysis was applied for the Braginsky [14] model by Roberts and Stewartson [15] near the onset of oscillatory convection using the small finite amplitude equations. Further, they have analyzed the linear stability analysis of the cubic and cubic-quintic amplitude equations. These amplitude equations are valid only when R is close to the R c . In addition, the nonlinear analysis proposed by Malkus and Veronis [4] is valid when the R is near the threshold R c . Later Kuo [5] proposed a different nonlinear approach, which is valid even for large values of R. The advantages of this approach are that the solutions are found to be valid even for a large range of the imposed temperature differences across the fluid layer and also the rapid convergence of solutions. This solution provides a quantitative theory for the convective heat transport as a function of the temperature difference in the range of laminar flow. Using this approach by Kuo [5], Rameshwar et al., [16], have studied the finite amplitude cellular convection in R B C under the influence of a vertical magnetic field and analyzed the magnetohydrodynamics ( M H D ) of electrically conducting fluid. Linear and nonlinear properties of thermohaline convection at the onset with the stress-free boundary conditions were investigated using perturbation analysis relevant to oceanic water and groundwater by Rawoof Sayeed and Rameshwar [17]. The stationary and oscillatory finite amplitude convections of a binary mixture with a porous medium were thoroughly investigated by Rameshwar et al., [18,19]. Baklouti et. al. [20] studied the dynamics of incompressible homogeneous turbulence by numerical simulations. Gupta et. al. [21] have analytically examined the effect of rotational speed modulation on the onset of magneto-thermal convection.
The extension studies related to geodynamo models proposed by Eltayeb [7,8], Roberts and Jones [13], and Jones and Roberts [22] have been studied by Šoltis and Brestenský [23]. The authors Šoltis and Brestenský [23] studied the influence of anisotropic diffusive coefficients (thermal diffusion and viscosity) on marginal stability of the horizontal fluid planar layer rotating about the vertical axis and permeated by a horizontal homogeneous magnetic field. The linear stability analysis was thoroughly investigated by the authors for two different types of anisotropic diffusive coefficients. This model was further extended to study the linear stability of the model of rotating magnetoconvection in the horizontal planar layer dynamo by Filippi et al. [24], which is influenced by three anisotropic diffusivities, such as viscosity, thermal diffusivity, and magnetic diffusivity.
In the present study, the nonlinear analysis was employed as proposed by Kuo [5] to study the behavior of cross rolls of electrically conducting fluid in a rotating magnetic system, given by Roberts and Jones [13]. Thus, the objectives of the present problem were as follows:
  • Investigate theoretically the nonlinear behavior of cross rolls which occur in the vertically rotating Rayleigh–Bénard convective system of planar layer of electrically conducting fluid in the presence of horizontal magnetic field;
  • Solve the nonlinear partial differential equations using the perturbation method proposed by Kuo [5], until the O( ε 8 ) and obtain the approximate solutions;
  • Obtain the combined effect of Lorenz and Coriolis forces with stress-free boundaries;
  • Find the local ( N L ) and average ( N u ) Nusselt numbers on the hot wall to understand the development of heat flow and the rate of heat transfer, respectively;
  • Obtain the cellular pattern of the fluid flow (streamlines) and hot regions (isotherms) from the eigenfunctions related to stream function and temperature, respectively;
  • Study the heatline patterns of the flow by using the heat function.
The novelty of the present work is the study of finite amplitude cellular convection when the stationary convection exists as a first instability. The dynamical behavior of the system depends on the type of instabilities that occur in that system. Interesting convective instabilities occur near the onset of convection which is analyzed from the linear stability analysis. At the onset of convection, the system is unstable to either stationary convection (at least one eigenvalue vanishes) or oscillatory convection (an eigenvalue with a purely imaginary part) as a first instability. When an eigenvalue vanishes the principle of exchange of stabilities occurs. In other words, a new steady state replaces the stable motionless state of the fluid. When stationary convection exists, the continuous release of potential energy is balanced by the viscous dissipation of mechanical energy and the convection always occurs in a fairly regular pattern. The results from the linear stability analysis of Roberts and Jones [13] show the occurrence of the modes such as parallel rolls, cross rolls, and oblique rolls based on the wave numbers. Only stationary convection exists as a first instability when the parallel rolls occur, but for the modes of cross rolls and oblique rolls, both stationary convection and oscillatory convection can occur as a first instability depending on the physical parameters. A detailed investigation of the linear stability analysis of the present physical model has been studied by Roberts and Jones [13]. Hence, the nonlinear dynamical behavior of the present considered physical model is investigated when stationary convection exists.
In Section 2, the basic governing equations that are considered by Roberts and Jones [13] are presented. The linear stability analysis is discussed in Section 3 to obtain critical Rayleigh numbers for steady cross roll modes. The nonlinear solutions for the field variables are presented in Section 4. In Section 5, the local Nusselt number ( N L ) and average Nusselt number ( N u ) are discussed. The distortion of streamlines and isotherms is discussed in Section 6. The heat flow visualization is discussed in Section 7. Finally, the conclusions are presented in Section 8.

2. Mathematical Model

In the present study the fluid with uniform density confined in an infinite horizontal layer was considered. It was assumed that the whole configuration rotates about the vertical axis ( O Z ) with angular velocity Ω ( = Ω 1 Z ) in the presence of a uniform gravitational field g ( = g 1 Z ) and the uniform magnetic field B ( = B 0 1 X ) applied in the horizontal direction where 1 X is the unit vector along the X-axis and 1 Z is the unit vector along Z-axis. The Prandtl number is assumed to be large, so as to ignore the inertial forces in the momentum equation in comparison to the Coriolis force [13]:
2 Ω ρ 0 × V = P + J × B μ m + g α ρ 0 T + μ 2 V ,
B t = × ( V × B ) + η 2 B ,
T t + V · ( T ) = κ 2 T ,
· V = 0 ,
· B = 0 ,
where P includes the centrifugal force, J = × B is the electric current density, and other notations are given the nomenclature.
We non-dimensionalized the Equations (1)–(5) using the corresponding length, time, velocity, temperature, magnetic field, and pressure scales as d , d 2 η , η d , β d , B 0 , and 2 Ω ρ 0 η , respectively.
Therefore, the non-dimensional governing equations are [13]:
1 Z × V = P + Λ J × B + q R T 1 Z + E 2 V ,
B t = × ( V × B ) + 2 B ,
T t + V · ( T ) = q 2 T ,
· V = 0 ,
· B = 0 ,
where R = β g α d 2 / 2 Ω k is the modified Rayleigh number, which measures the ratio of buoyancy force to Coriolis force, E = ν / 2 Ω d 2 is the ratio of viscous and Coriolis forces, Λ = B 0 2 / 2 Ω η μ m ρ 0 is the ratio of magnetic force and Coriolis force, and q = κ / η is the ratio between the thermal and magnetic diffusivities (the Roberts number). For the static solutions from Equations (6)–(10), we obtain
V s 0 , B s = 1 X , T s = Z .
After introducing the following perturbed quantities in the above static solutions, we obtain
V = V s + V * , B = B s + b * , T = T s + θ * q .
For convenience the asterisk symbols are omitted in the further analysis. The perturbed dimensionless governing equations are given by
E 4 W ω Z Z + R h 2 θ + Λ X 2 b Z = Λ 1 Z · × { × [ ( × b ) × b ] } ,
t + V · θ q W q 2 θ = 0 ,
t 2 b Z = W X + 1 Z · × ( V × b ) ,
t 2 J Z = ω Z X + 1 Z · × [ × ( V × b ) ] ,
E 2 ω Z + Λ J Z X + W Z + Λ 1 Z · × [ ( × b ) × b ] = 0 ,
where
h 2 = 2 X 2 + 2 Y 2 , 2 = 2 X 2 + 2 Y 2 + 2 Z 2 , V = ( U , V , W ) , b = ( b X , b Y , b Z ) , J Z = ( × b ) · 1 Z , ω Z = ( × V ) · 1 Z .
Eliminating θ , ω Z , b Z , and J Z from the linear part of Equations (13)–(17), we obtain
L W = N ,
L = L 1 + L 2 + L 3 + L 4 + L 5 + L 6 ,
and
N = N 1 + N 2 + N 3 + N 4 + N 5 ,
where
L 1 = t q 2 t 2 2 2 Z 2 ,
L 2 = E 2 t q 2 t 2 2 6 ,
L 3 = 2 E Λ t q 2 t 2 2 X 2 4 ,
L 4 = Λ 2 t q 2 4 X 4 2 ,
L 5 = R q E t 2 2 2 h 2 ,
L 6 = R q Λ t 2 2 X 2 h 2 ,
and
N 1 = E R t 2 2 2 h 2 + R Λ t 2 2 X 2 h 2 ( V · ) θ ,
N 2 = E Λ t q 2 t 2 4 X + Λ 2 t q 2 3 X 3 2 1 Z · × ( V × b ) ,
N 3 = Λ t q 2 t 2 2 X Z 1 Z · × [ × ( V × b ) ] ,
N 4 = Λ 2 t q 2 t 2 2 X 2 1 Z · × × [ ( × b ) × b ) ]
+ E Λ t q 2 t 2 2 2 1 Z · × × [ ( × b ) × b ) ] ,
N 5 = Λ t q 2 t 2 2 Z 1 Z · × [ ( × b ) × b ) ] .
Because the surfaces are maintained at a uniform temperature,
θ = 0 on Z = 0 and 1 , for all X , Y ,
and also normal component of the velocity should vanish on boundaries, i.e.,
W = 0 on Z = 0 and 1 , for all X , Y .
The conditions Equations (19) and (20) are independent of the nature of boundaries, such as free–free or rigid–rigid, etc. In the present work, we assumed stress-free boundary conditions [1], hence we obtain
2 W Z 2 = ω Z Z = 0 on Z = 0 and 1 , for all X , Y ,
J X = J Y = J Z Z = 0 on Z = 0 and 1 , for all X , Y .
Since the physical system is a triple diffusive system, it is unstable to either stationary convection or oscillatory convection near the onset.

3. Linear Stability Analysis

At the onset of convection, the existing perturbations in the system are very small. Hence, the nonlinear terms are smaller when compared to linear terms. The nonlinear contributions are neglected from Equation (18). We obtain a linear differential equation which is given as L W = 0 . This process is called linearization. We obtain
[ ( t q 2 ) ( t 2 ) 2 2 Z 2 + E 2 ( t q 2 ) ( t 2 ) 2 6 + 2 E Λ ( t q 2 ) ( t 2 ) ( 2 X 2 ) 4 + Λ 2 ( t q 2 ) ( 4 X 4 ) 2 + R ( q E ( t 2 ) 2 2 h 2 + q Λ ( t 2 ) ( 2 X 2 ) h 2 ) ] W = 0 .
The resulting Equation (23) is linear. The normal mode solution was considered as W ( X , Y , Z , t ) = W ( Z ) e i ( a X + l Y ) + p t , where a is the wavenumber along X direction and l is the wavenumber along Y direction. As such, a and l are real numbers and the growth rate (p) may be constant complex number [13]. The marginal state is obtained from R e ( p ) = 0 . The two types of modes are classified using the eigenvalue p, namely, if I m ( p ) = 0 , then the steady modes exist and if I m ( p ) = ω 0 , then the oscillatory convection exists. The preferred mode of convection depends on the physical parameters, which are relevant to the Earth’s outer core. In Earth’s outer core, the parameters E and q are considered as small and the Prandtl number is large. The orientation of rolls is classified based on the wavenumber. The modes are parallel rolls, if the wavenumber a = 0 (the axis of rolls are parallel to the applied magnetic field), if the wavenumber l = 0 , gives the cross rolls (the axis of the rolls are perpendicular to the applied magnetic field) and if both the wavenumbers a 0 and l 0 give the oblique rolls. The linear and nonlinear studies of the present physical model are based on l = 0 , i.e., cross rolls.
The physical parameters E, Λ , q, and R are used to study the linear and nonlinear behavior of the convective system. As the temperature gradient is increased, the unstable mode may be of stationary convection or oscillatory convection near the onset. We implemented the linear stability analysis using the normal mode analysis, i.e., by substituting W ( X , Y , Z , t ) = W ( z ) e ( i a X + p t ) in the linearized equation L W = 0 .

Stationary Convection ( ω = 0 )

By solving the linearized Equation (23), R = R s value is obtained for stationary convection and is given by
R = π d 2 2 + d 2 ( E d 2 2 + Λ a 2 ) 2 a 2 ( E d 2 2 + Λ a 2 ) ,
where d 2 = a 2 + π 2 . The critical value of R is obtained from R / a = 0 . The critical wavenumber is given by a 2 = a 2 c s = 2 π 2 , and the critical Rayleigh number for stationary convection is
R c s = 9 π 5 + 3 π 2 ( 9 E π 4 + 2 Λ π 2 ) 2 2 π 2 ( 9 E π 4 + 2 Λ π 2 ) .
From the above result the marginal Rayleigh number ( R c s ), E, and Λ values are obtained for high rotation rates and weak field [13] by the linear stability analysis. The critical values of control parameters were obtained to study the nonlinear behavior of cross rolls. For small values of q, stationary convection occurs and for large values of q oscillatory convection can occur. The value of R s is free from q while the finite amplitudes depend on q. For the nonlinear studies, a fixed q value as 0.01 was considered for which the stationary convection occurs as a first instability near the onset. For small values of the parameters E and Λ , the minimum value R c s decreases as E decreases and as Λ increases. Thus, the effect of E destabilizes the convective system when E decreases and the effect of Λ destabilizes the convective system when Λ increases (see Figure 1).
The linear stability theory adopts a less ambitious objective to ascertain when a flow is unstable to infinitesimal disturbances. It thus gives no prediction about transition promoted by sufficiently large disturbance. The ultimate consequence of the instability is never completely determined by linear theory. Thus, in the present study, an attempt was made to understand the nonlinear convection in the presence of the Coriolis force and magnetic field.

4. Method of Solution

The solutions of steady non-linear equations were obtained by following the method proposed by Kuo [5]. These solutions converge more rapidly and are valid for larger temperature differences. In this method, the dependent variables are first expressed as infinite series of a set of orthogonal space functions. This approach to the solution sheds light on the problem of transition to turbulent convection, which happens at a larger temperature difference. An expansion parameter ( ε ), is defined by [5]:
ε 2 = R R c s R .
Note that ε is less than one for all values of R . The solution of Equations (13)–(17) are written as
f = ε f 1 + ε 2 f 2 + ε 3 f 3 + ε 4 f 4 + ε 5 f 5 + ε 6 f 6 + ,
where
f = f ( U , W , θ , b X , b Z , ω Z , J Z ) .
According to Equation (25), R is given by
R = R c s 1 ε 2 ,
expanding Equation (27) in the power series of ε or by applying the finite formula
R = R c s + R o s ( ε 2 + ε 4 + ε 6 + + ε 2 s ) ,
where
R o s = R c s 1 ε 2 s , s = 1 , 2 , 3
By introducing Equations (26) and (28) in Equation (18), we obtain for the different ε orders, a sequence of linear non homogeneous differential equations as
O ( ε ) : ( L 1 + L 2 + L 3 + L 4 ) W 1 + R c s ( L 5 + L 6 ) W 1 = 0 ,
O ( ε 2 ) : ( L 1 + L 2 + L 3 + L 4 ) W 2 + R c s ( L 5 + L 6 ) W 2 = 0 ,
O ( ε 3 ) : ( L 1 + L 2 + L 3 + L 4 ) W i + R c s ( L 5 + L 6 ) W i + R o s ( L 5 + L 6 ) W i 2 + R o s ( L 5 + L 6 ) W i 4 = R c s N 1 + N 2 + N 3 + N 4 + N 5 , for i = 3 .
Similarly, at the orders O ( ε 4 ) , O ( ε 5 ) , we obtain
( L 1 + L 2 + L 3 + L 4 ) W i + R c s ( L 5 + L 6 ) W i + R o s ( L 5 + L 6 ) W i 2 + R o s ( L 5 + L 6 ) W i 4 = ( R c s + R o s ) N 1 + N 2 + N 3 + N 4 + N 5 , for i = 4 , 5 .
In general,
( L 1 + L 2 + L 3 + L 4 ) W i + R c s ( L 5 + L 6 ) W i + R o s ( L 5 + L 6 ) W i 2 + R o s ( L 5 + L 6 ) W i 4 = ( R c s + 2 R o s ) N 1 + N 2 + N 3 + N 4 + N 5 , for i 6 .
Here L i , i = 1 , 2 , 3 , 4 , 5 , 6 is the linear operator and N i , i = 1 , 2 , 3 , 4 , 5 represents the nonlinear terms and are functions of W i , θ i , b X i , b Z i , ω Z i , and J Z i . The auxiliary equations for temperature field are given by
t 2 θ 1 = q W 1 ,
t 2 θ 2 + ( V 1 · ) θ 1 = q W 2 .
In general,
t 2 θ i + l = 1 i 1 ( V l · ) θ i l = q W i , for i 3 .
The auxiliary equations for the magnetic field
t 2 b X 1 = U 1 X .
In general,
t 2 b X i = U i X + Y U 1 b Y i 1 + + U i 1 b Y 1 V 1 b X i 1 + V i 1 b X 1 + Z U 1 b Z i 1 + + U i 1 b Z 1 W 1 b X i 1 + W i 1 b X 1 , for i 2 ,
and
t 2 b Z 1 = W 1 X .
In general written as:
t 2 b Z i = W i X X U 1 b Z i 1 + + U i 1 b Z 1 W 1 b X i 1 + W i 1 b X 1 Y V 1 b Z i 1 + + V i 1 b Z 1 W 1 b Y i 1 + W i 1 b Y 1 for i 2 .
Likewise, the auxiliary equations for vorticity are
E 2 ω Z 1 = W 1 Z Λ J Z 1 X ,
in general,
E 2 ω Z i = W i Z Λ J Z i X Λ X b X 1 b Y i 1 X b X i 1 Y + Λ X b X i 1 b Y 1 X b X 1 Y + Λ X b Z 1 b Z i 1 Y b Y i 1 Z + + Λ X b Z i 1 b Z 1 Y b Y 1 Z + Λ Y b Z 1 b X i 1 Z b Z i 1 X + Λ Y b Z i 1 b X 1 Z b Z 1 X + Λ Y b Y 1 b Y i 1 X b X i 1 Y + + Λ Y b Y i 1 b Y 1 X b X 1 Y , for i 2 ,
and
t 2 J Z 1 = ω Z 1 X ,
in general written as:
t 2 J Z i = ω Z i X + 2 X Z V 1 b Z i 1 + + V i 1 b Z 1 W i 1 b Y 1 W 1 b Y i 1 + 2 Y Z W i 1 b X 1 + + W 1 b X 1 U 1 b Z i 1 U i 1 b Z 1 + h 2 V 1 b X i 1 + + V i 1 b X 1 U i 1 b Y 1 U 1 b Y i 1 , for i 2 .

4.1. Approximate Solutions

The approximate solutions U , W , θ , b X , and b Z are attained in terms of the amplitudes near the onset of stationary convection. The horizontal two boundaries are stress-free, all the space functions U , W , θ , b X , and b Z are sine and cosine functions. Thus, from Equations (30), (35), (42), and (44) we have, to the first order periodic solutions as,
W 1 = A 1 cos a X sin π Z , θ 1 = q A 1 π 2 + a 2 cos a X sin π Z ,
b X 1 = π A 1 π 2 + a 2 cos a X cos π Z , b Z 1 = a A 1 π 2 + a 2 sin a X sin π Z ,
where the nonlinear terms are used to calculate the amplitude A 1 . Normally, the terms in Equation (26) are written as
W i = A i cos a X sin π Z + p 1 , q 1 W p 1 q 1 ( i ) cos p 1 a X sin q 1 π Z ,
U i = 1 a 2 2 W i Z X , V i = 1 a 2 2 W i Z Y , ω Z i = V i X ,
θ i = q A i π 2 + a 2 cos a X sin π Z + p 1 , q 1 θ p 1 q 1 ( i ) cos p 1 a X sin q 1 π Z ,
b X i = π A i ( π 2 + a 2 ) cos a X cos π Z + p 1 , q 1 b X p 1 q 1 ( i ) cos p 1 a X cos q 1 π Z ,
b Z i = a A i ( π 2 + a 2 ) sin a X sin π Z + p 1 , q 1 b Z p 1 q 1 ( i ) sin p 1 a X sin q 1 π Z ,
where W p 1 q 1 ( i ) , θ p 1 q 1 ( i ) , b X p 1 q 1 ( i ) , and b Z p 1 q 1 ( i ) are nonlinear functions of A 1 , A 2 , A 3 ,… A i 1 . The unknown functions W p 1 q 1 ( i ) , θ p 1 q 1 ( i ) , b X p 1 q 1 ( i ) , and b Z p 1 q 1 ( i ) are calculated by substituting the Equations (47)–(51) in Equation (18).

4.2. Evaluation of Amplitude A 1

To obtain the second order solutions W 2 , θ 2 , b X 2 , and b Z 2 , the nonlinear term N is solved and Equation (31) is obtained. With N = 0, L W 2 = 0 and
W 02 ( 2 ) = 0 , θ 02 ( 2 ) = q A 1 2 8 π ( π 2 + a 2 ) a n d b X 02 ( 2 ) = A 1 2 4 ( π 2 + a 2 ) , b Z 02 ( 2 ) = 0 .
The unknown functions W 2 , θ 2 , b X 2 , and b Z 2 are obtained from Equations (31), (36), (39), and (41), respectively, and are given by
W 2 = A 2 cos a X sin π Z , θ 2 = q A 2 π 2 + a 2 cos a X sin π Z + θ 02 ( 2 ) A 1 2 sin 2 π Z , b X 2 = π A 2 π 2 + a 2 cos a X cos π Z + b X 02 ( 2 ) A 1 2 cos 2 π Z , b Z 2 = a A 2 π 2 + a 2 sin a X sin π Z .
From Equation (32) for i = 3 the amplitude A 1 is calculated. After using the first and second order solutions, we write Equation (32) as
( L 1 + L 2 + L 3 + L 4 ) W 3 + R c s ( L 5 + L 6 ) W 3 = R o s [ q E ( π 2 + a 2 ) 3 a 2 + q Λ ( π 2 + a 2 ) a 4 ] A 1 cos a X sin π Z R c s [ E d 2 3 a 2 + Λ d 2 a 4 ] π θ 02 ( 2 ) A 1 3 cos a X sin π Z
[ a 2 b X 02 ( 2 ) E Λ q d 2 4 + a 4 b X 02 ( 2 ) Λ 2 q d 2 2 ] A 1 3 2 cos a X sin π Z + R c s [ E ( a 2 + 9 π 2 ) 3 a 2 + Λ ( a 2 + 9 π 2 ) a 4 ] π θ 02 ( 2 ) A 1 3 cos a X sin 3 π Z + [ a 2 b X 02 ( 2 ) E Λ q ( a 2 + 9 π 2 ) 4 + a 4 b X 02 ( 2 ) Λ 2 q ( a 2 + 9 π 2 ) 2 ] A 1 3 2 cos a X sin 3 π Z [ 2 π 2 a 2 d 2 b X 02 ( 2 ) a 2 b X 02 ( 2 ) 2 ] A 1 3 cos a X sin π Z + [ 2 π 2 a 2 d 2 b X 02 ( 2 ) a 2 b X 02 ( 2 ) ( 2 ) ] A 1 3 cos a X sin 3 π Z ,
where d 2 = a 2 + π 2 . By solving the above equation, A 1 is given by
A 1 = 2 Λ b X 02 ( 2 ) π 2 Λ b X 02 ( 2 ) d 2 π R cs θ 02 ( 2 ) R os q 2 Λ b X 02 ( 2 ) π 2 Λ b X 02 ( 2 ) d 2 π R cs θ 02 ( 2 ) .
The unknown functions W 3 , θ 3 , b X 3 , and b Z 3 are obtained from Equations (32), (37), (39), and (41), respectively,
W 3 = A 3 cos a X sin π Z + W 13 ( 3 ) A 1 3 cos a X sin 3 π Z , θ 3 = q A 3 π 2 + a 2 + θ 11 ( 3 ) A 1 3 cos a X sin π Z + θ 13 ( 3 ) A 1 3 cos a X sin 3 π Z , b X 3 = π A 3 π 2 + a 2 + b X 11 ( 3 ) A 1 3 cos a X cos π Z + b X 13 ( 3 ) A 1 3 cos a X cos 3 π Z , b Z 3 = a A 3 π 2 + a 2 + b Z 11 ( 3 ) A 1 3 sin a X sin π Z + b Z 13 ( 3 ) A 1 3 sin a X sin 3 π Z .
where
W 13 ( 3 ) = 1 D 13 R cs q Λ d 13 a 4 + E q d 13 3 a 2 π θ 02 ( 2 )
+ 1 D 13 1 2 a 2 E Λ q d 13 4 + 1 2 a 4 Λ 2 q d 13 2 b X 02 ( 2 ) + 2 π 2 a 2 d 2 1 2 a 2 b X 02 ( 2 ) ,
θ 11 ( 3 ) = π θ 02 ( 2 ) d 2 , θ 13 ( 3 ) = π θ 02 ( 2 ) + W 13 ( 3 ) q d 13 ,
and
b X 11 ( 3 ) = π b X 02 ( 2 ) 2 d 2 , b X 13 ( 3 ) = 3 π W 13 ( 3 ) d 2 3 π b X 02 ( 2 ) 2 d 2 ,
b Z 11 ( 3 ) = a b X 02 ( 2 ) 2 d 2 , b Z 13 ( 3 ) = a W 13 ( 3 ) d 2 a b X 02 ( 2 ) 2 d 2 ,
D 13 = 9 q π 2 d 13 3 E 2 q d 13 6 2 E Λ q a 2 d 13 4
Λ 2 q 2 a 4 d 13 2 + R cs q a 2 E d 13 3 + q Λ a 4 d 13 ,
where d 13 = a 2 + 9 π 2 . This iterative procedure is continued to find A 2 to A 6 and the equivalent W i , θ i , b X i , b Z i , ω Z i , and J Z i , etc.

4.3. Calculation of A 2 and A 3

Initially, Equation (33) is simplified for i = 4 , i.e.,
( L 1 + L 2 + L 3 + L 4 + R c s ( L 5 + L 6 ) ) W 4 + R o s ( L 5 + L 6 ) W 2 = ( R c s + R o s ) N 1 + N 2 + N 3 + N 4 + N 5 .
The nonlinear terms N i , i = 1 , 2 , 3 , 4 , 5 are analyzed by using the Equations (46), (53) and (56). Solving the above Equation (57), we obtain A 2 = 0. This indicates that every second order approximate solutions vanish except for θ 02 ( 2 ) and b X 02 ( 2 ) as shown in Equation (52). Therefore, Equation (57) decreases to
( L 1 + L 2 + L 3 + L 4 + R c s ( L 5 + L 6 ) ) W 4 = K 1 A 1 4 cos 2 a X sin 4 π Z + K 2 A 1 4 cos 2 a X sin 2 π Z + K 3 A 1 A 3 cos 2 a X sin 2 π Z .
The approximate solutions W 4 , θ 4 , b X 4 , b Z 4 are evaluated by Equations (33), (37), (39) and (41) and those are given by
W 4 = A 4 cos a X sin π Z + K 4 A 1 4 cos 2 a X sin 4 π Z + K 5 A 1 4 cos 2 a X sin 2 π Z + K 6 A 1 A 3 cos 2 a X sin 2 π Z , θ 4 = q ( a 2 + π 2 ) A 4 cos a X sin π Z + K 7 A 1 A 3 sin 2 π Z + K 8 A 1 4 sin 2 π Z + K 9 A 1 4 cos 2 a X sin 2 π Z + K 10 A 1 4 cos 2 a X sin 4 π Z + K 11 A 1 4 sin 4 π Z + K 12 A 1 A 3 cos 2 a X sin 2 π Z , b X 4 = π ( a 2 + π 2 ) A 4 cos a X cos π Z + K 13 A 1 4 cos 2 a X cos 4 π Z + K 14 A 1 4 cos 2 a X cos 2 π Z + K 15 A 1 A 3 cos 2 a X cos 2 π Z , b Z 4 = a ( a 2 + π 2 ) A 4 sin a X sin π Z + K 16 A 1 4 sin 2 a X sin 4 π Z + K 17 A 1 4 sin 2 a X sin 2 π Z + K 18 A 1 A 3 sin 2 a X sin 2 π Z .
To determine the value of A 3 , Equation (33) is solved for i = 5, and is given by
( L 1 + L 2 + L 3 + L 4 + R c s ( L 5 + L 6 ) ) W 5 + R o s ( L 5 + L 6 ) ( W 3 + W 1 ) = ( R c s + R o s ) N 1 + N 2 + N 3 + N 4 + N 5 .
Evaluating A 3 from Equations (46), (53), (56), (59) and (60) we obtain,
A 3 = S 1 S 2 ,
where
S 1 = A 1 5 Λ 2 q d 2 2 a 3 + E Λ q d 2 4 a π K 17 4 a + K 14 4 b X 02 ( 2 ) W 13 3 2 3 π K 5 4 d 2 a
+ A 1 5 Λ 2 q d 2 2 a 2 + E Λ q d 2 4 a 2 π K 14 + 2 a K 17 4 d 2 π b X 02 ( 2 ) b Z 11 ( 3 ) π b X 02 ( 2 ) b Z 13 ( 3 ) + K 17 π a
+ A 1 5 a 2 Λ 2 q d 2 2 a 2 + E Λ q d 2 4 π 2 π K 14 + 2 a K 17 4 d 2 + b X 02 ( 2 ) 2 π b X 11 ( 3 ) + a b Z 11 ( 3 )
A 1 5 a 2 Λ 2 q d 2 2 a 2 + E Λ q d 2 4 b X 02 ( 2 ) 2 3 π b X 13 ( 3 ) + a b Z 13 ( 3 ) + π b X 02 ( 2 ) b X 11 ( 3 ) π b X 02 ( 2 ) b X 13 ( 3 )
+ q Λ d 2 a 4 + E q d 2 3 a 2 π θ 02 ( 2 ) W 13 ( 3 ) 3 π q K 5 4 d 2 π K 9 4 π K 8 A 1 5 R cs
+ q Λ d 2 a 4 + E q d 2 3 a 2 π θ 02 ( 2 ) W 13 ( 3 ) 3 π q K 5 4 d 2 π K 9 π K 8 A 1 5 R os
Λ R os a 4 d 2 q 2 E R os a 2 d 2 3 q 2 A 1 5 a 2 Λ 2 q d 2 2 a 2 + E Λ q d 2 4 K 14 ,
and
S 2 = q Λ d 2 a 4 + E q d 2 3 a 2 π K 8 π θ 02 ( 2 ) π q 3 K 5 4 A 1 2 R cs
q Λ d 2 a 4 + E q d 2 3 a 2 π K 8 π θ 02 ( 2 ) π q 3 K 5 4 A 1 2 R os
+ A 1 2 Λ 2 q d 2 2 a 3 + E Λ q d 2 4 a b X 02 ( 2 ) 2 π K 17 4 a + K 14 4 3 π K 5 4 d 2 a
A 1 2 Λ 2 q d 2 2 a 2 + E Λ q d 2 4 a 2 π K 14 2 a K 17 4 d 2 + π a b X 02 ( 2 ) d 2 + K 17 π a
A 1 2 Λ 2 q d 2 2 a 2 + E Λ q d 2 4 π 2 π K 14 + 2 a K 17 4 d 2 b X 02 ( 2 ) 2 a 2
A 1 2 Λ 2 q d 2 2 a 2 + E Λ q d 2 4 π 2 b X 02 ( 2 ) d 2 K 14 a 2 + Λ R os a 4 d 2 q 2 + E R os a 2 d 2 3 q 2 .
The amplitude A 3 is determined by using the first, second, third and fourth order approximate solutions. From Equations (60), (37), (39) and (41), the fifth order approximate solutions are obtained
W 5 = A 5 cos a X sin π Z + K 19 A 1 5 + K 20 A 1 2 A 3 cos 3 a X sin π Z + K 21 A 1 3 + K 22 A 1 5 + K 23 A 1 2 A 3 cos a X sin 3 π Z + K 24 A 1 5 + K 25 A 1 2 A 3 cos 3 a X sin 3 π Z + K 26 A 1 5 cos a X sin 5 π Z + K 27 A 1 5 cos 3 a X sin 5 π Z ,
θ 5 = q A 5 d 2 + K 28 A 1 2 A 3 cos a X sin π Z + K 29 A 1 3 + K 30 A 1 5 + K 31 A 1 2 A 3 cos a X sin 3 π Z + K 32 A 1 5 + K 33 A 1 2 A 3 cos 3 a X sin π Z + K 34 A 1 5 + K 35 A 1 2 A 3 cos 3 a X sin 3 π Z + K 36 A 1 5 cos a X sin 5 π Z + K 37 A 1 5 cos 3 a X sin 5 π Z ,
b X 5 = π A 5 d 2 + K 38 A 1 5 cos a X cos π Z + K 39 A 1 3 + K 40 A 1 5 + K 41 A 1 2 A 3 cos a X cos 3 π Z + K 42 A 1 5 + K 43 A 1 2 A 3 cos 3 a X cos π Z
+ K 44 A 1 5 + K 45 A 1 2 A 3 cos 3 a X cos 3 π Z + K 46 A 1 5 cos a X cos 5 π Z + K 47 A 1 5 cos 3 a X cos 5 π Z ,
b Z 5 = a A 5 d 2 + K 48 A 1 5 + K 49 A 1 2 A 3 sin a X sin π Z + K 50 A 1 3 + K 51 A 1 5 + K 52 A 1 2 A 3 sin a X sin 3 π Z + K 53 A 1 5 + K 54 A 1 2 A 3 sin 3 a X sin π Z + K 55 A 1 5 + K 56 A 1 2 A 3 sin 3 a X sin 3 π Z + K 57 A 1 5 sin a X sin 5 π Z + K 58 A 1 5 sin 3 a X sin 5 π Z .
Here the coefficients K i , i = 1 , 2 , 58 in Equations (58)–(65) are functions of a , E , Λ , and q . The simplifications become more critical as ε order increases. Similarly, the simplification was carried until the eighth order of Equation (34) to calculate for A 4 , A 5 , and A 6 . Proceeding as above, it can be observed that A 2 = A 4 = A 6 = 0 .

5. Convective Heat Transport

The changes in the two-dimensional flow patterns are illustrated by the local Nusselt number, N L , distributions over the heated plate. The heat transport coefficient in terms of the N L is expressed as [16,25]
N L = T n ,
here n denotes the normal direction on a plane. The heat transport is measured by the average Nusselt number ( N u ), which is independent of Z and is given by
N u = W T ¯ T ¯ Z ,
here the bar indicates a horizontal mean. Using Equation (67), N u is obtained by integrating the expression over the boundary, Z = 0 [16].
N u = 1 L 0 L W T T Z | z = 0 d X ,
where L is the normalized horizontal cell width.

5.1. Local Nusselt Number ( N L )

Figure 2a illustrates the changes of N L for distinct values of E with respect to X. In selected regions, the number of peaks and the location of maximum and minimum of N L values depend on E. The maximum of N L value is constant, as X increases for defined E. However, as E decreases, the number of peaks is increased in a selected region. Figure 2b illustrates the variation of N L for different Λ concerning to X. The location of the maximum N L value is independent of the dimensionless plate length but depends on Λ . The existence of the number of peaks in the given region of X increases, as Λ increases. Figure 2b illustrates the heat that is transported from the fluid to the boundary is increased as Λ increase.

5.2. Average Nusselt Number ( N u )

The dependency of N u on the control parameters was studied near and away from the onset of stationary convection. Let N u ( 2 ) ( s = 1 ) , N u ( 4 ) ( s = 2 ) and N u ( 6 ) ( s = 3 ) , indicate the second-, fourth- and sixth-order approximations, for N u , respectively. The second order approximation is given by
N u ( 2 ) = 1 + 2 q R o s 2 Λ ( a 2 π 2 ) + R c s .
The approximations for N u ( 4 ) and N u ( 6 ) are lengthy, so it is not shown here to conserve space. The change of N u with respect to R is shown in Figure 3 for different values of E and Λ . Figure 3a shows the effect of E on N u in the R plane for fixed Λ = 0.2 and q = 0.01. It illustrates that the rate of heat transfer is enhanced for decreasing E. Figure 3b shows the effect of Λ on heat transfer rate for fixed E = 0.005 and q = 0.01. The enhancement of heat transfer is observed for increasing Λ values. The small values of q are relevant to Earth’s outer core. It is very difficult to perform numerical simulations for small values of the physical parameters. For q < 1 stationary convection is preferred [13]. Figure 3c shows that for q < 1 and as q increases, the N u increase. Thus the effect of q < 1 shows, the heat transfer rate are enhanced and accordingly the intensity of the flow rate also increase. The change of kinetic energy with respect to R / R c s is represented in Figure 4. Figure 4a demonstrates the change in N u for various E values as well as for a fixed value of Λ = 0.2 . The change in E produces small change in the potential energy, in comparison with the kinetic energy. Thus, the total energy decreases as E increase. Figure 4b shows the energy distribution for different Λ and for fixed E = 0.005 and q = 0.01 . The amplifying values of Λ show that the total energy is also increased.

6. Distortion of Streamlines and Isotherms

The fluid flow behavior is visualized by the stream function ( Ψ ) which is obtained from the velocity components U and W. The relation between the velocity components and stream function ( Ψ ) is [26]
U = Ψ Z and W = Ψ X ,
which produce a single equation
Ψ = 2 Ψ = 0 for X = 0 , 2 π / a and Z = 0 , 1 .
The points with equal temperature connected with lines are called isotherms. The snapshots of the heat transport and flow field near the onset of stationary convection are expressed in terms of streamlines and isotherms.
The general attributes of the streamlines and isotherms with respect to the variation in R, E, and Λ are shown in Figure 5, Figure 6, Figure 7 and Figure 8. Figure 5 illustrates the pattern of streamlines and isotherms near the onset of convection ( R R c s ). Figure 5a, shows the pattern of streamlines for E = 0.005 , Λ = 0.2 , q = 0.01 . The absolute maximum and the absolute minimum values of circular strengths are 0.11215 and 0.11216 , respectively. Figure 5c shows the pattern of streamlines for E = 0.01 , Λ = 0.2 . This figure shows the absolute maximum and the absolute minimum values of circular strengths as 0.11243 and 0.11243 , respectively. From Figure 5a,c the maximum strength of rolls at R R c s are decreased as E increases. Thus, as E increases, R c s also increases, accordingly N u decreases and hence the absolute maximum of circulation strength decreases. Figure 5e illustrate the pattern of streamlines with E = 0.005 and Λ = 0.4 . These streamlines have the absolute minimum and maximum values with the circular strengths as 0.12116 and 0.121117 , respectively. By comparing Figure 5a,e, the periodic rectangular rolls are observed near the R R c s , but as Λ increases the maximum strength of rolls is increased and the minimum strength of roll decreases. Thus, as Λ increases, R c s decreases, accordingly N u increases, and hence the absolute maximum of circulation strength increases. Figure 5a–f are plotted for the values of R R c s and the flow pattern are rectangular rolls and follows the symmetric nature over the range of 0 X 1 . Since the stream function equations show the symmetric property. Similarly, for the same values of E and Λ , the isotherms formed as horizontal lines near R R c s , as shown in Figure 5b,d,f. At R R c s , the strength of isotherms is of small magnitude, representing the conduction dominant heat transport inside the considered region. These isotherms are smooth lines that span over the whole region.
The snapshot of streamlines and isotherms for different values of R and for fixed values of E = 0.005 , Λ = 0.2 , and q = 0.01 are displayed in Figure 6a–h. It is observed that for R = 10 R c s and for the cell lying between 0 X 1 , the absolute minimum and maximum values are with the circulation strengths 43.4965 and 43.4965 , respectively, as shown in Figure 6a. As R increases from R c s to 10 R c s , the basic cells become more deformed due to the growth of two vortices B and B located at the top right and bottom left boundaries with the circulation strength 16.315 . The basic cell with two vortices A and A has circulation strength 38.059 . The temperature profiles in terms of isotherms are illustrated in Figure 6b for same values of physical parameters that are considered in Figure 6a. The isotherms are of nearly in wavy shape with the absolute maximum and minimum values of 0.00212 and 1.00271 , respectively. It indicates the maximum of heat transfer process is occurred by convection. Figure 6c,d illustrate the streamlines and isotherms for R = 20 R c s . The temperature gradient and the gravitational buoyancy force act together and changes the flow structure as shown in Figure 6c. The bicellular patterns of streamlines turn out to be multicellular models and these cells divide the field of motion at the core for a cell lies between 0 X 1 with the absolute maximum and minimum values of circulation strengths 220.829 and 220.845 , respectively. The vortices B and B showed their presence with 82.822 as the circulation strength in the opposite direction of an original cell. For these considered values of physical parameters, the behavior of isotherms is shown in Figure 6d, which exhibit the mode of convective heat transport inside the fluid layer. In the fluid layer, the absolute maximum and minimum values of isotherms are respectively, 0.00275 and 1.0052 . When R is increased from 20 R c s to 30 R c s , the small vortices B and B shown in Figure 6c are increased with circulation strength 398.740 . Thus, the basic cell encountered more deformation (Figure 6e) and has the absolute maximum and minimum values at 638.172 and 638.027 , respectively. Accordingly, the isotherm curves develop more deformation. The absolute maximum and minimum values of isotherms in the layer are, respectively, 0.00449 and 1.01101 . As R is increased from 30 R c s to 80 R c s (Figure 6g), the two vortices B and B grow in size and split the basic cell into two vortices located on either side of the secondary cell with the absolute maximum ( 9976.78 ) and minimum ( 9977.74 ) strengths. The heat flow pattern becomes chaotic, which is shown in Figure 6h when R increases to 80 R c s . The absolute maximum and minimum values of isotherms in the layer are, respectively, 8.75285 and 9.67347 . From Figure 6, it is observed that as R increases from R c s to 80 R c s , the onset of turbulent flows are producible.
Figure 7a–d, illustrate the streamlines and isotherms for different values of E and for a fixed set of other parameters R = 20 R c s , Λ = 0.2 , and q = 0.01 . The behavior of the flow field was investigated by considering the flow pattern in the region 0 X 1 , as shown in Figure 6a–d (E = 0.005 ) and Figure 7a–d (E = 0.01 and 0.015 ). Figure 7a shows streamlines for E = 0.01 in the considered range of X. The absolute maximum and minimum values of circulation strengths are 162.15 and 162.149 , respectively. There exist two vortices B and B outside the basic cell with the circulation strength 82.806 as shown in Figure 7a. The basic cell also contains two vortices namely A and A with a circulation strength of 141.881 . Figure 7c is plotted for E = 0.015 , which has the absolute maximum and absolute minimum values of circulation strength as 140.453 and 140.471 , respectively. The flow pattern in the region 0 X 1 contains a deformed basic cell due to the growth of two vortices B and B that exist at either side of the basic cell and are located at the top and bottom boundaries with circulation strength 60.81 . The basic cell also has two vortices A and A with a circulation strength value of 122.90 . Finally from Figure 6c and Figure 7a,c it is observed that the strength of the basic cell and pattern deformation decrease as E increases. This implies that the effect of E stabilizes the convective system. The flow of heat transfer is shown in Figure 7b,d for E = 0.01 and 0.015 , respectively.
Figure 8a–d show the streamlines and isotherms for different Λ values and for E = 0.005 , R = 20 R c s , and q = 0.01 . The effect of Λ was studied from Figure 6c,d and Figure 8a–d. In Figure 8a the streamlines are plotted for Λ = 0.4 . By considering the flow pattern in the range of 0 X 1 , the absolute maximum and minimum values of circulation strengths are 436.721 and 436.747 , respectively. In this range, the basic cell is deformed by two vortices B and B , which are located at the top right and bottom left of the layer and on either side of the basic cell with circulation strength 272.97 . The basic cell also encloses two vortices A and A with strength 382.13 . Figure 8c shows the streamlines for Λ = 0.6 in the considered range of 0 X 1 , with the absolute maximum and minimum values of circulation strength 663.595 and 663.555 , respectively. In addition, there exist two vortices B and B with circulation strength 414.71 . The basic cell also enclosed two vortices A and A with circulation strength 580.65 . As Λ increases from 0.4 to 0.6 the deformation and circulation strength of cells ( A , A ) increase. This implies that the effect of Λ destabilizes the convective system. The isotherms are plotted in Figure 6d and Figure 8b,d for distinct values of Λ = 0.2 , 0.4 and 0.6 . The lines of isotherms change to a more circular form as Λ increases. Thus, the incremental values of Λ destabilize the convective system.

Topology of Flow

The topology constraint is based on the Euler number ( ζ ) of the flow. As described by Jana et al. [27], ζ on the surface is defined as the sum of the Poincare indices of the critical points on the surface and is given by
N E ( N H + 1 2 N P ) = ζ ,
here the N E represents the number of elliptic points, N H is the number of hyperbolic points, and N P is the number of parabolic points [28,29]. In Figure 9a, the vorticity contours are exhibited for R R c s , E= 0.005, Λ = 0.2 and q = 0.01 . The present simulated flow fulfils the topological rule given in Equation (70) with N P = 0 , N H = 2 , and N E = 2 . For R = 20 R c s , an equivalent investigation has been done for vorticity contours in Figure 9b and Equation (70) is satisfied with N P = 0, N E = 8, and N H = 8.

7. Heat Function

Heatlines depict the convective heat transport phenomenon, whereas the isotherms are mainly useful for visualizing heat transfer in the domain of conduction. The heat function and heatline analyzes were developed by Kimura and Bejan [30] to visualize heat transmission through the fluid flow, later Morega and Bejan [31] successfully used the concept of heatlines. Different researchers [32,33,34,35] used this concept for dissimilar applications of natural convective systems. The heat function ( H * ) is defined as
H * X = W T T Z ,
H * Z = U T T X ,
where T = T s + θ and T s = T 0 Z . The Equations (71) and (72) do not exhibit the symmetric property. Differentiating Equations (71) and (72) with respect to X and Z, respectively, and subtracting the resulting equations yields
2 H * X 2 + 2 H * Z 2 = ( W T ) X ( U T ) Z .
From the definition of heat function, Equations (71) and (72), the boundary conditions on H * follow [16]:
H * ( X , 0 ) = H * ( 0 , 0 ) + 0 X ( W T T Z ) d X , at Z = 0 and 0 X 2 π / a ,
H * ( X , 1 ) = H * ( 0 , 1 ) + 0 X ( W T T Z ) d X , at Z = 1 and 0 X 2 π / a ,
H * ( 0 , Z ) = H * ( 0 , 0 ) 0 Z ( U T T Z ) d Z , at X = 0 and 0 Z 1 ,
H * ( 2 π / a , Z ) = H * ( 2 π / a , 0 ) 0 Z ( U T T Z ) d Z , at X = 2 π / a and 0 Z 1 .

Results and Discussion for Heatlines

Figure 10a–d, illustrate the pattern of heatlines for different Rayleigh number values, R R c s , R = 1.05 R c s , 1.15 R c s and 1.25 R c s , respectively, for fixed values of E = 0.005 , Λ = 0.2 , and q = 0.01 . When the system is at a conduction state ( R < R c s ) heatlines are always parallel to Z-axis and perpendicular to isotherms. Figure 10a illustrates the heatlines for R R c s . It is observed that the heatline contours within the domain are normal to the Z = 0 and Z = 1 lines due to conduction dominant heat transfer. For R R c s , the absolute maximum and minimum values of heatlines are 10.2224 and 0.0241 , respectively, in the considered range 0 X 5 . In the neighborhood of X = 0 , the heatlines at the center of the system depict the structure which is similar to the parabolic structure. The curvature at the central part of the system increases as X increases. This shows that the nonlinear propagation of heat transfer occurs when R R c s . Hence, the transition takes place from the conduction state to the convection state at R R c s . Figure 10b is plotted for R = 1.05 R c s . The absolute maximum and minimum values of heatlines are 10.3081 and 0.2774 , respectively, in the considered range 0 X 5 . The heatline with a strength of 0.18 exist near the line X = 0 and the heatline with strength 9.69 exist at X = 5 . The strength of heatlines increases as X increases. Some heatlines occurred in the form of a closed path. As R increases from R c s the heatlines with same strength are changed to a closed path as shown in Figure 10a,b. For higher values, this indicates that the convective heat flow is more intense at the center. Figure 10c shows heatlines at R = 1.15 R c s and having the absolute maximum and minimum values of heatlines 10.3666 and 0.62655 , respectively. In this figure, the number of closed paths of heatlines at the center is increased for R = 1.15 R c s . The size of closed path of heatlines for R = 1.15 R c s is more than that for R = 1.05 R c s . Figure 10d shows the heatlines for R = 1.25 R c s with the absolute maximum and minimum values of 10.3972 and 0.88096 , respectively. The number of closed paths of heatlines increases at the center for R = 1.25 R c s . The size of closed path of heatlines for R = 1.25 R c s is increased in comparison with that of the heatlines for R = 1.15 R c s . For large R, the convective heat transmission is more intense. It is observed that the heatlines become denser with the increase in R. Figure 10a–d indicate that the heat transfer across the layer is increased as R increases. Heatlines will not exhibit periodic patterns due to the non-symmetry nature of Equations (71) and (72).
Figure 10b and Figure 11a,b are plotted with the same strength of heatlines so as to analyze the influence of E on heat flow for the fixed values of Λ , R, and q. In the considered range of 0 X 5 , for E = 0.01 , (Figure 11a) the absolute maximum and minimum values of heatlines are noted to be 8.4264 and 0.18554 , respectively, and for E = 0.015 (Figure 11b) these values are 7.51142 and 0.11574 , respectively. In both of these Figure 11a,b, the heatlines intensity decays with E. The size of the closed path and the number of closed paths with the same strength decreased as E increases. From Figure 10b and Figure 11a,b it is fascinating to observe the inhibition of temperature in the central regime as E increases.
Figure 10b and Figure 12a,b are plotted with the same strength of heatlines to analyze the effect of Λ on heat flow for the fixed values of E, R, and q. In the considered range of 0 X 5 for Λ = 0.4 (Figure 12a) with the absolute maximum and minimum values of heatlines as 8.55642 and 0.23723 , respectively. These values for Λ = 0.6 (Figure 12b) are 7.74043 and 0.18296 , respectively. In both Figure 12a,b, the heatlines are dominated by convection and form closed loops. The size and the number of closed loops with the same strength increase as Λ increases. From Figure 10b and Figure 12a,b it is observed that the temperature within the central regime is enhanced as Λ increase by observing heatlines.

8. Conclusions

The nonlinear natural convection was studied in a planer layer of electrically conducting fluid that rotates about the vertical axis in the presence of a uniform horizontal magnetic field and vertical temperature gradient. This problem has applications in Earth’s liquid core. The present results help to enhance understanding of the finite amplitude convection when the coupling between the Lorentz force and the Coriolis force present in nonlinear planar layer convection-driven dynamos.
  • Linear stability analysis showed that as the small values of E keep decreasing or Λ increasing, the R c s decreases, i.e., the effect of E stabilizes and Λ destabilizes the system.
  • Theoretically investigated the nonlinear behavior of cross rolls that occur in the Rayleigh–Bénard convective system of a planar layer dynamo of electrically conducting fluid rotating about the vertical axis in the presence of a horizontal magnetic field.
  • The nonlinear partial differential equations was solved using the perturbation method, until O ( ϵ 8 ) and obtained the approximate solutions.
  • Computed the local Nusselt number ( N L ) and averaged Nusselt number ( N u ) on the hot wall to understand the development of heat flow and the rate of heat transfer, respectively.
  • The number of peaks is fixed for a given E while the value of the peak is independent of X for a given E. The absolute peak values of N L increase as E increases. The number of peaks is fixed for a given Λ . The value of the peak is independent of X for a given Λ . The absolute peak value of N L increases as Λ decrease. From the N L results, it is noted that the heat flux is high for decreasing E or increasing Λ .
  • It is observed that the Ekman number (E) generates a strong damping effect on heat transfer at high rotation rates but the heat transport enhances as Λ increases. The Roberts number (q) < 1, enhances the heat transfer rate and accordingly the intensity of the flow rate also increases. Similarly, the total energy decays as E increases. The increment in the values of Λ show, the increase in the total energy.
  • Obtained the cellular pattern of fluid flow (streamlines) and the hot regions (isotherms) from the eigenfunctions related to the stream function and temperature, respectively. From the streamlines and isotherms trajectories, it is observed that, for the lower values of E the deformation of the fluid pattern is enhanced and more transfer of heat in the flow occurs due to the presence of lesser viscous force in comparison with the Coriolis force. Similarly, for the amplifying values of Λ , there is more deformation in the streamlines and isotherms. This result shows, in the presence of Coriolis force, the magnetic field destabilizes the system.
  • Studied the heatline patterns of the flow by using the heat function. The results show that the deformation in the trajectories of heatlines are enhanced as E decreases. A similar trend of deforming heatlines is observed with increasing Λ .

Author Contributions

Conceptualization, Y.R. and G.S.; Methodology, Y.R.; Software, U.S.M. and D.L.; Validation, Y.R., G.S, U.S.M., A.K.R. and D.L.; Visualization A.K.R. and U.S.M.; Writing—Original Draft Preparation, G.S.; Writing—Review and Editing, Y.R., A.K.R. and D.L.; Funding Acquisition, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

D.L. acknowledges the partial financial support from Centers of excellence with BASAL/ANID financing, Grant ANID ABF220001, CEDENNA.

Data Availability Statement

All data underlying the results are available as part of the article and no additional source data are required.

Acknowledgments

The authors thank H. P. Rani (National Institute of Technology, Warangal) for her critical reading, editing, and improving the manuscript’s English grammar.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AAmplitude
B Magnetic field
B s Static magnetic field
B 0 Characteristic field strength
aWavenumber
PModified pressure
EEkman number
1 Z Unit vector along Z-axis
1 X Unit vector along X-axis
qRatio of thermal and magnetic diffusivities
dThe convective zone depth
g Gravitational field
N u Nusselt number
H*Heat function
RModified Rayleigh number
R c Critical Rayleigh number
R c s Critical Rayleigh number for stationary convection
T 0 Reference temperature
TTemperature field
T s Static temperature
Δ T Temperature difference between top and bottom layers
V Velocity field
V s Static velocity
U , V , W Velocity components
X , Y , Z Cartesian coordinates
tTime
R B C Rayleigh–Bénard Convection
Greek symbol
Λ Elsasser number
β Adverse temperature gradient
θ Perturbed temperature
η Magnetic diffusivity
ρ Density
ρ 0 Reference density
κ Coefficient of thermal diffusivity
ν Kinematic viscosity
αThermal expansion coefficient
μDynamic viscosity
μ m Magnetic permeability
ω Vorticity field
Ω Angular velocity
ω Frequency of oscillations
Superscript
Dimensional form
*Perturbed quantities

References

  1. Chandrasekar, S. Hydrodynamic and Hydromagnetic Stability; Oxford Clarendon Press: Oxford, UK, 1961. [Google Scholar]
  2. Robert, P.H. An Introduction to Magnetohydrodynamics; American Elsevier: New York, NY, USA, 1967. [Google Scholar]
  3. Cox, S.M.; Mathews, P.C. New instabilities in two-dimensional rotating convection and magnetoconvection. Phys. D 2001, 149, 210. [Google Scholar] [CrossRef]
  4. Malkus, W.V.R.; Veronis, G. Finite Amplitude Cellular Convection. J. Fluid Mech. 1958, 4, 225–260. [Google Scholar] [CrossRef]
  5. Kuo, H.L. Solution of the non-linear equations of the cellular convection and heat transport. J. Fluid Mech. 1960, 10, 611–630. [Google Scholar] [CrossRef]
  6. Braginsky, S.I. Magnetohydrodynamics of the Earth’s core. Geomagn. Aeron. 1964, 4, 698–712. [Google Scholar]
  7. Eltayeb, I.A. Hydromagnetic convection in a rapidly rotating fluid layer. Proc. R. Soc. Lond. A 1972, 326, 229–254. [Google Scholar]
  8. Eltayeb, I.A. Overstable hydromagnetic convection in a rotating fluid layer. J. Fluid Mech. 1975, 71, 161–179. [Google Scholar] [CrossRef]
  9. Aurnou, J.M.; Olson, P.L. Experiments on Rayleigh–Bénard convection, magnetoconvection and rotating magnetoconvection in liquid gallium. J. Fluid Mech. 2000, 430, 283–307. [Google Scholar] [CrossRef] [Green Version]
  10. Raju, C.S.K.; Ameer Ahammad, N.; Sajjan, K.; Shah, N.A.; Yook, S.; Dinesh Kumar, M. Nonlinear movements of axisymmetric ternary hybrid nanofluids in a thermally radiated expanding or contracting permeable Darcy Walls with different shapes and densities: Simple linear regression. Int. Commun. Heat Mass Trans. 2022, 135, 106110. [Google Scholar] [CrossRef]
  11. Kumar, M.D.; Raju, C.S.K.; Sajjan, K.; El-Zahar, E.R.; Shah, N.A. Linear and quadratic convection on 3D flow with transpiration and hybrid nanoparticles. Int. Commun. Heat Mass Trans. 2022, 134, 105995. [Google Scholar] [CrossRef]
  12. Rani, H.P.; Rameshwar, Y.; Brestensky, J. Topology of Rayleigh-Bénard convection and magnetoconvection in plane layer. Geophys. Astrophys. Fluid Dyn. 2019, 113, 208–221. [Google Scholar] [CrossRef]
  13. Roberts, P.H.; Jones, C.A. The onset of magnetoconvection at large Prandtl number in a rotating layer 1. Finite magnetic diffusion. Geophys. Astrophys. Fluid Dyn. 2000, 92, 289–325. [Google Scholar] [CrossRef]
  14. Braginsky, S.I. Torsional magnetohydrodynamic vibrations in the Earth’s core and variations in the day length. Geomagn. Aeron. 1970, 10, 3–12. [Google Scholar]
  15. Robert, P.H.; Stewartson, K. On Finite Amplitude Convection in a Rotaiting Magnetic System. Philos. Trans. R. Soc. Lond. Ser. Math. Phys. Sci. 1974, 277, 287–315. [Google Scholar]
  16. Rameshwar, Y.; Rawoof Sayeed, M.A.; Rani, H.P.; Laroze, D. Finite amplitude cellular convection under the influence of a vertical magnetic field. Int. J. Heat Mass Transf. 2017, 114, 559–577. [Google Scholar] [CrossRef]
  17. Rawoof Sayeed, M.A.; Rameshwar, Y. Finite Amplitude Cellular Thermohaline Convection. J. Heat Transf. 2022, 114, 112602. [Google Scholar] [CrossRef]
  18. Rameshwar, Y.; Srinivas, G.; Laroze, D.; Rawoof Sayeed, M.A.; Rani, H.P. Convective instabilities in binary mixture 3He-4He in porous media. Chin. J. Phys. 2022, 77, 773–803. [Google Scholar] [CrossRef]
  19. Rameshwar, Y.; Srinivas, G.; Laroze, D. Finite amplitude oscillatory convection of binary mixture kept in a porous medium. Processes 2023, 11, 664. [Google Scholar] [CrossRef]
  20. Baklouti, F.S.; Khlifi, A.; Salhi, A.; Godeferd, F.; Cambon, C.; Lehner, T. Kinetic-magnetic energy exchanges in rotating magnetohydrodynamic turbulence. J. Turbul. 2019, 20, 263–284. [Google Scholar] [CrossRef]
  21. Gupta, V.K.; Keshri, O.P.; Kumar, A. Effect of rotational speed modulation on weakly nonlinear magneto convective heat transfer with temperature-dependent viscosity. Chin. J. Phys. 2021, 72, 487–498. [Google Scholar] [CrossRef]
  22. Jones, C.A.; Roberts, P.H. The onset of magnetoconvection at large Prandtl number in a rotating layer. II. Small magnetic diffusion. Geophys. Astrophys. Fluid Dyn. 2000, 93, 173–226. [Google Scholar] [CrossRef]
  23. Šoltis, T.; Brestenský, J. Rotating magnetoconvection with anisotropic diffusivities in the Earth’s core. Phys. Earth Planet. Int. 2010, 178, 27–38. [Google Scholar]
  24. Filippi, E.; Brestenský, J.; Šoltis, T. Effects of anisotropic diffusion on onset of rotating magnetoconvection in plane layer; stationary modes. Geophys. Astrophys. Fluid Dyn. 2019, 113, 80–106. [Google Scholar] [CrossRef]
  25. Sparrow, E.M.; Carlson, C.K. Local and average natural convection Nusselt numbers for a uniformly heated, shrouded or unshrouded horizontal plate. Int. J. Heat Mass Transf. 1986, 29, 369–379. [Google Scholar] [CrossRef]
  26. Batchelor, G.K. An Introduction to Fluid Dynamics; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  27. Jana, S.C.; Metcalfe, G.; Ottino, J.M. Experimental and computational studies of mixing in complex Stokes flows: The vortex mixing flow and multicellular cavity flows. J. Fluid Mech. 1994, 269, 199–249. [Google Scholar] [CrossRef]
  28. Tony Sheu, W.H.; Rani, H.P. Exploration of vortex dynamics for transitional flows in a three-dimensional backward facing step channel. J. Fluid Mech. 2006, 550, 61–83. [Google Scholar] [CrossRef]
  29. Sheu, T.; Rani, H.P.; Ten-China, T.; Tsai, S.F. Multilple states, topology and bifurcations of natural convection in a cubical cavity. Comput. Fluids 2008, 37, 1011–1028. [Google Scholar] [CrossRef]
  30. Kimura, S.; Bejan, A. The heatline visualization of convective heat transfer. J. Heat Transf. 1983, 105, 916–919. [Google Scholar] [CrossRef]
  31. Morega, A.I.; Bejan, A. Heatline visualization of forced convection laminar boundry layers. Int. J. Heat Mass Transf. 1993, 36, 3957–3966. [Google Scholar] [CrossRef]
  32. Bejan, A. Convection Heat Transfer; Wiley: New York, NY, USA, 1984; pp. 21–23. [Google Scholar]
  33. Komori, K.; Kito, S.; Naumara, T.; Inaguma, Y.; Inagai, T. Fluid flow and heat transfer in the transition process of natural convection over an inclined plate. Heat Trans. Asian Res. 2001, 30, 648–659. [Google Scholar] [CrossRef]
  34. Kimura, F.; Kitamura, K.; Yamaguchi, M.; Asami, T. Fluid flow and heat transfer of natural convection adjacent to upward facing inclined heated plates. Heat Trans. Asian Res. 2003, 32, 278–291. [Google Scholar] [CrossRef]
  35. Hooman, K.; Gurgechi, H.; Dincer, I. Heatline and energy-flux-vector visualization of natural convection in a porous cavity occupied by a fluid with temperature-dependent viscosity. J. Porous Media 2009, 12, 265–275. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The effect of Λ and E on R c s .
Figure 1. The effect of Λ and E on R c s .
Processes 11 01945 g001
Figure 2. Variation of N L with respect to X. (a) Λ = 0.2 and q = 0.01 for different E, (b) E = 0.005 and q = 0.01 for different Λ .
Figure 2. Variation of N L with respect to X. (a) Λ = 0.2 and q = 0.01 for different E, (b) E = 0.005 and q = 0.01 for different Λ .
Processes 11 01945 g002
Figure 3. Dependence of Nussult number (Nu) on Rayleigh number ( R / R c s ). (a) Λ = 0.2 and q = 0.01 for different E, (b) E = 0.005 and q = 0.01 for different Λ , (c) E = 0.005 and Λ = 0.2 for different q.
Figure 3. Dependence of Nussult number (Nu) on Rayleigh number ( R / R c s ). (a) Λ = 0.2 and q = 0.01 for different E, (b) E = 0.005 and q = 0.01 for different Λ , (c) E = 0.005 and Λ = 0.2 for different q.
Processes 11 01945 g003
Figure 4. Dependence of kinetic energy and potential energy on Rayleigh number ( R / R c s ) are plotted. (a) Λ = 0.2 and q = 0.01 for different E, (b) E = 0.005 and q = 0.01 for different Λ .
Figure 4. Dependence of kinetic energy and potential energy on Rayleigh number ( R / R c s ) are plotted. (a) Λ = 0.2 and q = 0.01 for different E, (b) E = 0.005 and q = 0.01 for different Λ .
Processes 11 01945 g004
Figure 5. The Effect of E and Λ near R c s , streamlines (a) for E = 0.005, Λ = 0.2, and q = 0.01 , (c) E = 0.01, Λ = 0.2, and q = 0.01 , (e) E = 0.005, Λ = 0.4, and q = 0.01 and isotherms (b) E = 0.005, Λ = 0.2, and q = 0.01 , (d) E = 0.01, Λ = 0.2, and q = 0.01 , (f) E = 0.005, Λ = 0.4, and q = 0.01 are plotted.
Figure 5. The Effect of E and Λ near R c s , streamlines (a) for E = 0.005, Λ = 0.2, and q = 0.01 , (c) E = 0.01, Λ = 0.2, and q = 0.01 , (e) E = 0.005, Λ = 0.4, and q = 0.01 and isotherms (b) E = 0.005, Λ = 0.2, and q = 0.01 , (d) E = 0.01, Λ = 0.2, and q = 0.01 , (f) E = 0.005, Λ = 0.4, and q = 0.01 are plotted.
Processes 11 01945 g005aProcesses 11 01945 g005b
Figure 6. The Effect of R = 10 R c s , 20 R c s , 30 R c s and 80 R c s , streamlines (a,c,e,g) and isotherms (b,d, f, h) are plotted for E = 0.005, Λ = 0.2, and q = 0.01.
Figure 6. The Effect of R = 10 R c s , 20 R c s , 30 R c s and 80 R c s , streamlines (a,c,e,g) and isotherms (b,d, f, h) are plotted for E = 0.005, Λ = 0.2, and q = 0.01.
Processes 11 01945 g006aProcesses 11 01945 g006b
Figure 7. The Effect of E = 0.01 and 0.015, streamlines (a,c) and isotherms (b,d) are plotted for Λ = 0.2, R = 20 R c s , and q = 0.01 .
Figure 7. The Effect of E = 0.01 and 0.015, streamlines (a,c) and isotherms (b,d) are plotted for Λ = 0.2, R = 20 R c s , and q = 0.01 .
Processes 11 01945 g007
Figure 8. The Effect of Λ = 0.4 and 0.6, streamlines (a,c) and isotherms (b,d) are plotted for E = 0.005, R = 20 R c s and q = 0.01.
Figure 8. The Effect of Λ = 0.4 and 0.6, streamlines (a,c) and isotherms (b,d) are plotted for E = 0.005, R = 20 R c s and q = 0.01.
Processes 11 01945 g008
Figure 9. Vorticity lines for (a) R R c s , E = 0.005, Λ = 0.2, and q = 0.01 , (b) R = 20 R c s , E = 0.005, Λ = 0.2, and q = 0.01 .
Figure 9. Vorticity lines for (a) R R c s , E = 0.005, Λ = 0.2, and q = 0.01 , (b) R = 20 R c s , E = 0.005, Λ = 0.2, and q = 0.01 .
Processes 11 01945 g009
Figure 10. The effect of R R c s , R = 1.05 R c s , 1.15 R c s and 1.25 R c s , heatlines (ad) are plotted for E = 0.005, Λ = 0.2, T 0 = 1, and q = 0.01, respectively.
Figure 10. The effect of R R c s , R = 1.05 R c s , 1.15 R c s and 1.25 R c s , heatlines (ad) are plotted for E = 0.005, Λ = 0.2, T 0 = 1, and q = 0.01, respectively.
Processes 11 01945 g010
Figure 11. The effect of E = 0.01 and 0.015, heatlines (a,b) are plotted for Λ = 0.2, R = 1.05 R c s , T 0 = 1, and q = 0.01 .
Figure 11. The effect of E = 0.01 and 0.015, heatlines (a,b) are plotted for Λ = 0.2, R = 1.05 R c s , T 0 = 1, and q = 0.01 .
Processes 11 01945 g011
Figure 12. The effect of Λ = 0.4 and 0.6, heatlines (a,b) are plotted for E = 0.005, R = 1.05 R c s , T 0 = 1, and q = 0.01 .
Figure 12. The effect of Λ = 0.4 and 0.6, heatlines (a,b) are plotted for E = 0.005, R = 1.05 R c s , T 0 = 1, and q = 0.01 .
Processes 11 01945 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rameshwar, Y.; Srinivas, G.; Rao, A.K.; Mahabaleshwar, U.S.; Laroze, D. Nonlinear Analysis of Cross Rolls of Electrically Conducting Fluid under an Applied Magnetic Field with Rotation. Processes 2023, 11, 1945. https://doi.org/10.3390/pr11071945

AMA Style

Rameshwar Y, Srinivas G, Rao AK, Mahabaleshwar US, Laroze D. Nonlinear Analysis of Cross Rolls of Electrically Conducting Fluid under an Applied Magnetic Field with Rotation. Processes. 2023; 11(7):1945. https://doi.org/10.3390/pr11071945

Chicago/Turabian Style

Rameshwar, Y., G. Srinivas, A. Krishna Rao, U. S. Mahabaleshwar, and D. Laroze. 2023. "Nonlinear Analysis of Cross Rolls of Electrically Conducting Fluid under an Applied Magnetic Field with Rotation" Processes 11, no. 7: 1945. https://doi.org/10.3390/pr11071945

APA Style

Rameshwar, Y., Srinivas, G., Rao, A. K., Mahabaleshwar, U. S., & Laroze, D. (2023). Nonlinear Analysis of Cross Rolls of Electrically Conducting Fluid under an Applied Magnetic Field with Rotation. Processes, 11(7), 1945. https://doi.org/10.3390/pr11071945

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop