Next Article in Journal
Scaling up the Natural Mode of Action of Macrophyte Allelochemicals and Their Effect on Toxic Cyanobacteria Using a Nitrogen-Limited Chemostat
Next Article in Special Issue
The Influence of Pre-Lift Gate Opening on the Internal and External Flow Characteristics During the Startup Process of an Axial Flow Pump
Previous Article in Journal
A Method for Optimizing Production Layer Regrouping Based on a Genetic Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Rotation and Viscosity on Parallel Rolls of Electrically Conducting Fluid

1
Instituto de Alta Investigación, Universidad de Tarapacá, Casilla 7D, Arica 1000000, Chile
2
Department of Mathematics, University College of Science, Osmania University, Hyderabad 500007, India
*
Author to whom correspondence should be addressed.
Processes 2024, 12(9), 1882; https://doi.org/10.3390/pr12091882
Submission received: 16 July 2024 / Revised: 21 August 2024 / Accepted: 28 August 2024 / Published: 2 September 2024
(This article belongs to the Special Issue Complex Fluid Dynamics Modeling and Simulation, 2nd Edition)

Abstract

:
Rayleigh–Bénard convection is a fundamental fluid dynamics phenomenon that significantly influences heat transfer in various natural and industrial processes, such as geophysical dynamics in the Earth’s liquid core and the performance of heat exchangers. Understanding the behavior of conductive fluids under the influence of heating, rotation, and magnetic fields is critical for improving thermal management systems. Utilizing the Boussinesq approximation, this study theoretically examines the nonlinear convection of a planar layer of conductive liquid that is heated from below and subjected to rotation about a vertical axis in the presence of a magnetic field. We focus on the onset of stationary convection as the temperature difference applied across the planar layer increases. Our theoretical approach investigates the formation of parallel rolls aligned with the magnetic field under free–free boundary conditions. To analyze the system of nonlinear equations, we expand the dependent variables in a series of orthogonal functions and express the coefficients of these functions as power series in a parameter ϵ . A solution for this nonlinear problem is derived through Fourier analysis of perturbations, extending to O( ϵ 8 ), which allows for a detailed visualization of the parallel rolls. Graphical results are presented to explore the dependence of the Nusselt number on the Rayleigh number (R) and Ekman number (E). We observe that both the local Nusselt number and average Nusselt number increase as the Ekman number decreases. Furthermore, the flow appears to become more deformed as E decreases, suggesting an increased influence of external factors such as rotation. This deformation may enhance mixing within the fluid, thereby improving heat transfer between different regions.

1. Introduction

The magneto-convection issue explored in this paper is motivated by research into the magnetohydrodynamics of the Earth’s outer core, which involves the behavior of liquid metals. The flow of electric current across a magnetic field is associated with a body force called the Lorentz force, which influences the fluid flow. As the strength of the magnetic field increases, the Lorentz forces become more dominant compared to the Coriolis force. Given that convection in a rotating and magnetized environment requires a delicate balance of forces, interpreting the results of these models can be quite complex. The geodynamo functions at a low Ekman number (E), as this condition is pertinent to geophysical scenarios. Likewise, the growth rate of the parallel rolls remains unaffected by the magnetic field strength ( Λ ). When Λ , the absence of a magnetic field means there is no preferential horizontal direction, resulting in a horizontally isotropic layer that corresponds to the Λ -independent parallel mode. Instabilities that can develop in a rapidly rotating magnetohydrodynamic system, such as the Earth’s core, are often regarded as essential to dynamo theory. This area of research examines how the magnetic field is preserved within the system even in the presence of Ohmic dissipation.
In recent years, notable progress has been made in the creation of dynamos driven by convection. Chendrasekhar [1] examined the linear stability of the layer under the simultaneous action of rotation and magnetic field when both the angular velocity ( Ω ) and applied magnetic field ( B 0 ) are vertically aligned and the boundaries are stress-free. The numerical results indicated that if the Taylor number (T) is fixed and Chandrasekhar number (Q) increases from zero, the critical Rayleigh number ( R c ) remains constant until Q reaches a certain value before decreasing. When Q is increased further, R c decreases, reaching a minimum before increasing again. Roberts and Jones [2] examined how convection begins when the temperature difference across the layer is increased, focusing on the planform and time-dependent behavior of small amplitude convection. Their analysis was guided by a similar study that used specific boundary conditions for illustration. Eltayeb [3] studied the linear stability of a rotating and electrically conducting viscous layer in the Rayleigh–Bénard system with a uniform magnetic field, applying the Boussinesq approximation. He considered several orientations of the magnetic field and rotation axes under various surface conditions. For example, when the magnetic field is horizontal and the rotation is vertical, the rolls align parallel to the magnetic field, provided that the field is weak enough. Childress and Soward [4] presented a hydromagnetic dynamo by examining Bénard convection between rotating parallel planes. Their approach was based on an asymptotic expansion involving two spatial scales, which is suitable for scenarios with large Taylor numbers. Soward [5] has discussed various models of thermal convection in rapidly rotating fluids permeated by strong magnetic fields. Braginsky [6] highlighted the crucial influence of nonuniform rotation in the core. He analyzed small oscillations in electrically conducting fluid affected by magnetic, buoyancy (Archimedean), and Coriolis forces, incorporating Ohmic diffusion to reflect conditions similar to those in the Earth’s fluid core. Glatzmaier and Roberts [7] developed and analyzed the first three-dimensional self-consistent numerical solution of the MHD equations, which model thermal convection and magnetic field generation in a rotating spherical fluid shell with a conducting solid inner core. Their results demonstrate a self-sustaining dynamo that can sustain a magnetic field for three different magnetic diffusion times. Jones and Roberts [8] investigated three-dimensional convection-driven dynamos through numerical simulations. Their study uncovered a rich variety of dynamical behaviors that necessitate systematic examination. They focused on both the kinematic and magnetically saturated regimes, proposing a simplified model of the dynamo mechanism. Furthermore, they analyzed the influence of the parameters R and E and the diffusivity ratio q on the dynamo, while also considering the significance of the Taylor constraint in low E convection-driven dynamos. Their results included a snapshot of the flow in the (y,z) plane, revealing a strong correlation between vertical velocity and temperature, as anticipated. Taylor [9] has discussed the influence of magnetic force on the behavior of a rotating incompressible conducting fluid. Soward [10] studied weakly nonlinear convection in Braginsky’s model. Geiger and Busse [11] studied the removal of the degeneracy of the patterns of convective motion in a spherically symmetric fluid shell by the effects of rotation. Matthews [12] investigated simple steady flows generated by thermal convection to determine which of them can generate a magnetic field through dynamo action. Kono and Robert [13] described ways of avoiding ambiguity and presented a preferred definition of the R that is particularly well suited for geodynamo simulations. Cattaneo and Hughes [14] studied dynamo processes in a convective layer of Boussinesq fluid rotating about the vertical, and measured the average electromotive force using an externally imposed uniform field. Roberts and Kono [15] presented a new method for evaluating the total energy release. They showed that the gravitational source that stirs the fluid core is less than 30 percent of the total gravitational energy released through the contraction of the Earth as it cools. Tilgner [16] presented numerical simulations of dynamos in rotating Rayleigh–Bénard convection in plane layers.
Earlier research suggests that further investigation is needed into the relationship between the heat transfer rates and control parameters associated with magnetohydrodynamic stability. There is also potential for studying flow visualization in electrically conducting fluids rotating about a vertical axis under the influence of a magnetic field.
The current study is particularly relevant to the research conducted by Roberts and Jones [2], which examined the onset of convection. This model involves a vertically rotating Rayleigh–Bénard layer of electrically conducting fluid that generates a magnetic field (( B 0 )) in the horizontal direction. We analytically solve the nonlinear governing equations using the perturbation method proposed by Kuo [17], computing approximate solutions for the dependent variables up to O( ϵ 8 ) with free–free boundary conditions. The flow behavior in our dynamo models with stress-free boundary conditions closely resembles that observed with rigid boundary conditions. Our investigation centers on the onset of convection and assesses how various parameters influence fluid flow, temperature distribution, and heat transfer rates. To visualize the fluid flow, temperature distribution, and heat transfer, we utilize the streamline, isotherm, and heatline concepts initially introduced by Kimura and Bejan [18]. Recent studies, such as that by Rameshwar et al. [19], have performed nonlinear analysis of cross-rolls in electrically conducting fluids subjected to an applied magnetic field with rotation. Khedher et al. [20] explored the effects of heating and magnetohydrodynamics (MHD) on the periodic behavior of current density and heat transfer amplitude in electrically conducting fluids. Shuguang et al. [21] analyzed heat transfer in the MHD stagnation point flow of a Maxwell nanofluid over a rotating porous disk, noting that the rotating disk phenomenon leads to the formation of a circular arc at every point on the disk during uniform rotation about an axis. Additionally, Shamshuddin et al. [22] investigated entropy generation in the flow of a tangent hyperbolic magnetized hybrid nanofluid over a stretchable rotating disk, incorporating physical effects such as Ohmic dissipation, nonlinear thermal radiation, heat dissipation, and entropy analysis. Taghreed et al. [23] numerically assessed the combined influence of the wedge angle and melting energy transfer on the flow of tangent hyperbolic magnetohydrodynamic (MHD) nanofluid across a permeable wedge, and also examined the effects of thermal radiation and heat source/sink on the flow of tangent hyperbolic nanofluid (THNF). In their research, Ramanuja et al. [24] examined the steady-state behavior of an electrically conducting and compressible fully developed viscous Jeffrey fluid flowing through infinitely parallel porous vertical microchannels and subjected to a transverse magnetic field. Their study modeled the fluid dynamics using the Navier–Stokes equations combined with energy conservation principles, enabling a detailed analysis of the flow characteristics and thermal behavior under these specific conditions. Qianqian Tang et al. [25] addressed the challenge of enhancing the thermal storage stability of nanosuspension concentrates (SC) derived from low-melting-point pesticides. Their study focused on the preparation of a 25 wt% pyraclostrobin nanosuspension concentrate using a water-based grinding technique, with pyraclostrobin serving as the primary raw material. Jun-Xia Li et al. [26] reported the synthesis of two cobalt(II) complexes via a solvothermal reaction using cobalt chloride hexahydrate and triclopyr acetic acid. Single-crystal X-ray diffraction analysis indicated that both complexes crystallized in the triclinic system, specifically within a defined space group. Yun-Long Wu et al. [27] explored metal–organic framework (MOF)-supported metal nanoparticle (MNP) composites (MNPs@MOFs), which leverage the synergistic effects of both materials. They outlined four main synthesis strategies for these composites: ship-in-the-bottle, bottle-around-the-ship, one-pot, and sandwich assembly. Their review summarizes recent developments in these methods for designing and fabricating various MNPs@MOFs composites.
This version provides a clearer structure and emphasizes the significance of the research questions in relation to the introduction:
  • How does the rotation speed of the electrically conducting fluid affect the onset of convection in a vertically oriented Rayleigh–Bénard layer?
  • What role do control parameters such as fluid viscosity and thermal conductivity play in the stability of magnetohydrodynamic convection?
  • How can the streamline, isotherm, and heatline visualization techniques enhance our understanding of flow dynamics and heat transfer in MHD systems?
  • What are the implications of entropy generation in the flow of magnetized hybrid nanofluids over rotating disks with regard to energy efficiency in engineering applications?
This research work is organized as follows: Section 2 provides the magnetohydrodynamic equations in the Boussinesq approximation for the description of the flow; Section 3 presents the linear stability analysis; the method of solution is discussed in Section 4; Section 5 presents the derivation of approximate solutions; analytical approximate expressions for N u are obtained in Section 6; the patterns of streamlines and isotherms are presented in Section 7; Section 8 deals with the heatfunction and heatline patterns; and Section 9 deals with entropy generation. Finally, the conclusions are provided in Section 10.

2. Mathematical Model

Electrically conducting fluids of uniform density are positioned in infinite horizontal layers that rotate about the vertical ( O Z ) axis with an angular velocity expressed as Ω ( = Ω 1 z ) . A uniform magnetic field B 0 ( = B 0 1 y ) is assumed to be applied in the horizontal direction in the presence of uniform gravitational field g (= g 1 z ) (Figure 1). The momentum equation neglects inertial forces as compared to Coriolis forces. This consideration arises from the assumption of a large Prandtl number ( P r ). The mathematical equations for the model, based on the Boussinesq approximation, are presented as follows [2]:
2 Ω ρ 0 × V = P + J × B μ m + g α ρ 0 T + μ 2 V ,
B t = × ( V × B ) + η 2 B ,
T t + V · ( T ) = κ 2 T ,
· B = 0 ,
· V = 0 ,
where J = × B is the dimensional electric current density. To make Equations (1)–(4) non-dimensional, the following scale parameters are taken into consideration [2]:
X = 1 d X , Y = 1 d Y , Z = 1 d Z , t = η d 2 t , u = d η u ,
v = d η v , w = d η w , T = T β d , B = 1 B 0 B , p = 1 2 Ω ρ 0 η p .
Equations (1)–(4) then take the following dimensionless form:
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 J = × B . Equations (6)–(10) are employed to examine the linear stability of a static state:
V s = 0 , B s = 1 Y , T s = Z .
After perturbation, these variables are expressed as dependent variables, provided by
V = V s + V * , B = B s + b * , θ = T s + θ * q .
In the following, the asterisk symbol will be removed for simplicity in the analysis below, and the perturbed dimensionless equations become
E 4 w w Z Z + q 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 B Z = w Z X + 1 Z · × [ × ( V × b ) ] ,
E 2 w Z + Λ B Z X + w Z + Λ 1 Z · × [ ( × b ) × b ] = 0 ,
where
= 1 X X + 1 Y Y + 1 Z Z , h 2 = 2 X 2 + 2 Y 2 ,
2 = 2 X 2 + 2 Y 2 + 2 Z 2 , B z = ( b y X b x Y ) , w z = ( v X u Y ) .
Here, V = V ( u , v , w ) is the perturbed velocity, b = b ( b X , b Y , b Z ) is the perturbed magnetic field and θ is the perturbed temperature field, w = w ( w x , w y , w z ) is the vorticity vector, and the other notations are as explained in the Nomenclature. In Equations (13)–(17), the control physical parameters are as follows: R = β g α d 2 / 2 Ω κ is the Rayleigh number, accounting for buoyancy; E = ν / 2 Ω d 2 is the Ekman number, which measures the viscous force against the Coriolis force; and Λ = B 0 2 / 2 Ω η μ m ρ 0 is the Elsasser number, which measures the magnetic forces against Coriolis force. The fourth parameter, q = κ / η , denotes the ratio of the thermal diffusivity to the magnetic diffusivity.
By removing θ , w Z , b Z , and B Z from the linear section 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 .
The stress-free boundary conditions are specified on the horizontal planes at Z = 0 and Z = 1 . Because the surface temperature is kept constant, we have
θ = 0 on Z = 0 and Z = 1 for   every X , Y .
At the boundaries, the standard component of velocity should also be equal to zero, i.e.,
w = 0 on Z = 0 and Z = 1 for   every X , Y .
The conditions in (19) and (20) are independent of any type of boundaries. The stress-free boundary conditions are considered; hence, we have the following:
2 w Z 2 = w z Z = 0 on Z = 0 and Z = 1 for   every X , Y ,
B X = B Y = B Z Z = 0 on Z = 0 and Z = 1 for   every X , Y .

3. Linear Stability Analysis

At the onset of convection, the system has only slight disturbances, causing the nonlinear terms to be considerably smaller than the linear terms. As a result, nonlinear contributions from Equation (18) are neglected. We obtain the linear differential equation, with the solution being L w = 0 ; thus,
t q 2 t 2 2 2 Z 2 + E 2 t q 2 t 2 2 6 + R q E t 2 2 2 h 2 w = 0 .
The resulting equation is linear. We define w ( X , Y , Z , t ) = W ( Z ) e i k X + i a Y + λ t , where k and a are constant wavenumbers. The modes with k = 0 are called parallel rolls, which are independent of the magnetic field strength ( Λ ). When Λ = 0 , the magnetic force does not impart a preferred horizontal direction. The growth rate ( λ ) is a constant that remains real for parallel rolls, indicating that they are never overstable [2]. The linear stability analysis is explored by replacing w ( Y , Z , t ) = W ( Z ) e i a Y + λ t in the linearized equation L w = 0 , using the normal mode analysis.

Stationary Convection

Equation (23) provides the Rayleigh number for stationary mode by taking λ = 0 , and is provided by
R = π 2 + d 2 3 E 2 a 2 E ,
where d 2 = a 2 + π 2 . Figure 2 displays R as a function of a for various Ekman number (E) values. It can be seen that as E decreases, the critical value R c s also diminishes. This suggests that a lower Ekman number (E) contributes to a destabilizing effect on the flow. The critical value of R ( R c s ) is found by setting R / a 2 = 0 . The critical wavenumber is provided by a 2 = a 2 c s = 2 π 2 , and the critical value of ( R c s ) of R for stationary convection is defined as
R c s = 27 E 2 π 4 + 1 2 E .

4. Method of Solution

According to Kuo [17], a solution can be achieved for the nonlinear system of equations by means of two expansions. First, the dependent variables are expanded in a series of orthogonal functions. Second, the coefficients associated with these functions are expanded in a power series of a parameter ϵ < 1 as
ϵ 2 = R R c s R .
The solutions for Equations (13)–(17) are represented as follows:
f = ϵ f 1 + ϵ 2 f 2 + ϵ 3 f 3 + ϵ 4 f 4 + ϵ 5 f 5 + ϵ 6 f 6 + ,
where
f = f ( v , w , θ , b Y , b Z , w Z , B Z ) .
The value of R according to Equation (24) is expressed as
R = R c s 1 ϵ 2 ,
and is either expanded in a power series in ϵ or the finite formula is applied as
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 (25) and (27) into Equation (13), the following different-order equations are obtained:
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 c s N 1 + N 2 + N 3 + N 4 + N 5 , for i = 3 .
O ( ϵ 4 ) : ( 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, the linear operators in L i , i = 1 , 2 , 3 , 4 , 5 , 6 and the nonlinear terms in N i , i = 1 , 2 , 3 , 4 , 5 are the functions of v i , w i , θ i , b Y i and b Z i . Each v i , w i and θ i must satisfy the boundary conditions in Equations (19)–(21). In addition, θ i , b Y i and b Z i are evaluated from v i and w i using an auxiliary equation derived from the Equations (31)–(33). The spectral equation for the temperature is provided 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 .
Similarly, the equations for the magnetic field are formulated as follows:
t 2 b Y 1 = v 1 X ,
t 2 b Y i = v i X + X v 1 b X i 1 + + v i 1 b X 1 u 1 b Y i 1 + u i 1 b Y 1 + Z 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 ,
t 2 b Z 1 = w 1 X ,
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 .

5. Discussion of the Solution

The approximate solutions for w , θ , b Y , and b Z are calculated in terms of amplitude. For the case with two free boundaries, the spatial functions associated with the different modes w , θ , b Y , and b Z are represented by sine and cosine functions. Hence, from Equations (29), (34), (37) and (39), the first order of the approximate solution is provided by
w 1 = A 1 cos a Y sin π Z , v 1 = 1 a 2 2 w 1 z y , θ 1 = q π 2 + a 2 A 1 cos a Y sin π Z , b Y 1 = 0 , b Z 1 = 0 ,
where the amplitude A 1 is determined from the nonlinear terms. Generally, the terms in Equation (25) may be written as
w i = A i cos a Y sin π Z + p 1 , q 1 w p 1 q 1 ( i ) cos p 1 a Y sin q 1 π Z ,
v i = 1 a 2 2 w i z y ,
θ i = q A i π 2 + a 2 cos a Y sin π Z + p 1 , q 1 θ p 1 q 1 ( i ) cos p 1 a Y sin q 1 π Z ,
b Y i = 0 ,
b Z i = 0 ,
where w p 1 q 1 ( i ) and θ p 1 q 1 ( i ) are nonlinear functions of the amplitudes A 1 , A 2 , A 3 …… A i 1 . The unknown functions w p 1 q 1 ( i ) , θ p 1 q 1 ( i ) (and the amplitudes A 1 , A 2 , A 3 …… A i 1 ) are obtained by substituting Equations (25) and (27) in Equations (13)–(17). We utilize the principle that the coefficients of each power of ϵ must equal zero.

5.1. Evaluation of Amplitude A 1

To solve the second-order solutions w 2 , θ 2 , b X 2 , and b Z 2 , the nonlinear term N needs to be calculated. From Equation (30), it is observed that N = 0; this leads to L w 2 = 0 and
w 02 ( 2 ) = 0 and θ 02 ( 2 ) = q π A 1 2 8 π 2 ( π 2 + a 2 ) .
The unknown functions w 2 , θ 2 , b Y 2 , and b Z 2 are obtained from Equations (30), (35), (38) and (40), respectively, and are expressed as follows:
w 2 = A 2 cos a Y sin π Z , v 2 = 1 a 2 2 w 2 z y , θ 2 = q A 2 π 2 + a 2 cos a Y sin π Z + θ 02 ( 2 ) A 1 2 sin 2 π Z , b Y 2 = 0 , b Z 2 = 0 .
From Equation (31), for i = 3 , the amplitude A 1 is calculated as
( L 1 + L 2 + R c s L 5 ) w 3 + R o s L 5 w 1 = R c s [ E q ( a 2 + 9 π 2 ) 3 a 2 ] π θ 02 ( 2 ) A 1 3 cos a Y sin 3 π Z R c s [ E q ( a 2 + π 2 ) 3 a 2 ] π θ 02 ( 2 ) A 1 3 cos a Y sin π Z .
Solving the above equation provides A 1 , and can be expressed by
A 1 = R cs π θ 02 ( 2 ) R os q R cs π θ 02 ( 2 ) .
The unknown functions w 3 and θ 3 are obtained from Equations (31) and (36), respectively, which are provided by
w 3 = A 3 cos a Y sin π Z + w 13 ( 3 ) cos a Y sin 3 π Z , v 3 = 1 a 2 2 w 3 z y , θ 3 = q A 3 π 2 + a 2 cos a Y sin π Z + θ 11 ( 3 ) cos a Y sin π Z + θ 13 ( 3 ) cos a Y sin 3 π Z + θ 02 ( 3 ) sin 2 π Z ,
where
w 13 ( 3 ) = R cs E q 9 π 2 + a 2 2 a 2 π θ 02 ( 2 ) D 13
and
θ 11 ( 3 ) = π θ 02 ( 2 ) d 2 ,
θ 13 ( 3 ) = q 2 9 π 2 + a 2 R cs E a 2 π θ 02 ( 2 ) 9 q 9 π 2 + a 2 3 π 2 E 2 q 9 π 2 + a 2 6 + R cs q 2 E 9 π 2 + a 2 3 a 2 π θ 02 ( 2 ) 9 π 2 + a 2 .
This iterative process is continued to find A 2 , A 3 , A 4 , A 5 , A 6 and the corresponding w, θ , b X , b Z .

5.2. Evaluation of Amplitudes A 2 and A 3

Initially, Equation (32) is solved for i = 4 , i.e.,
( L 1 + L 2 + R c s L 5 ) w 4 + R o s L 5 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 determined by utilizing Equations (41), (48) and (51). From the solvability condition of Equation (52), we have A 2 = 0 . This means that all second-order solution vanish except for θ 02 ( 2 ) , which is shown in Equation (48). Therefore, Equation (52) reduces to
( L 1 + L 2 + R c s L 5 ) w 4 + R o s L 5 w 2 = K 1 A 1 4 cos 2 a Y sin 4 π Z + K 2 A 1 4 cos 2 a Y sin 2 π Z .
The coefficients K 1 and K 2 in Equation (53) are functions of a , E , and q , while w 4 , θ 4 are evaluated by Equations (32) and (36) and are provided by
w 4 = A 4 cos a Y sin π Z + w 24 ( 4 ) A 1 4 cos 2 a Y sin 4 π Z + w 22 ( 4 ) A 1 4 cos 2 a Y sin 2 π Z , v 4 = 1 a 2 2 w 4 z y ,
θ 4 = q ( a 2 + π 2 ) A 4 cos a Y sin π Z + θ 02 ( 4 ) A 1 A 3 sin 2 π Z + θ 02 ( 4 ) A 1 4 sin 2 π Z + θ 22 ( 4 ) A 1 4 cos 2 a Y sin 2 π Z + θ 24 ( 4 ) A 1 4 cos 2 a Y sin 4 π Z + θ 04 ( 4 ) A 1 4 sin 4 π Z
where
w 24 ( 4 ) = 4 E q d 24 3 a 2 R cs + R os D 24 1 2 π θ 13 ( 3 ) + 1 4 q a v 13 ( 3 ) d 2 + 1 4 q π w 13 ( 3 ) d 2 ,
w 22 ( 4 ) = 4 E q d 22 3 a 2 R cs + R os D 22 π θ 13 ( 3 ) 1 4 q a v 13 ( 3 ) d 2 + 1 4 q π w 13 ( 3 ) d 2 ,
θ 02 ( 4 ) = 1 4 1 π 2 1 2 π θ 13 ( 3 ) + 1 4 q a v 13 ( 3 ) d 2 + 1 4 q π w 13 ( 3 ) d 2 + 1 2 π θ 11 ( 3 ) ,
θ 04 ( 4 ) = 1 16 1 π 2 π θ 13 ( 3 ) 1 4 q a v 13 ( 3 ) d 2 + 1 4 q π w 13 ( 3 ) d 2 ,
θ 22 ( 4 ) = 1 d 22 q w 22 ( 4 ) + π θ 13 ( 3 ) + 1 4 q a v 13 ( 3 ) d 2 1 4 q π w 13 ( 3 ) d 2 ,
θ 24 ( 4 ) = 1 d 24 q w 24 ( 4 ) 1 2 π θ 13 ( 3 ) 1 4 q a v 13 ( 3 ) d 2 1 4 q π w 13 ( 3 ) d 2 ,
and
D 24 = 16 q d 24 3 π 2 E 2 q d 24 6 + 4 R cs q 2 E d 24 3 a 2 ,
D 22 = 4 q d 22 3 π 2 E 2 q d 22 6 + 4 R cs q 2 E d 22 3 a 2 ,
d 22 = 4 a 2 + 4 π 2 , d 24 = 4 a 2 + 16 π 2 .
To determine the value of A 3 , the fifth-order equation for i = 5 of Equation (32) is solved, which can be written as follows:
( L 1 + L 2 + R c s L 5 ) w 5 + R o s L 5 ( w 3 + w 1 ) = ( R c s + R o s ) N 1 .
Evaluating N 1 from Equations (41), (48), (51), (54) and (55) we obtain
N 1 = D 1 A 1 5 cos a Y sin 5 π Z + D 2 A 1 5 cos a Y sin 3 π Z + D 3 A 1 5 cos 3 a Y sin 5 π Z + D 4 A 1 5 cos 3 a Y sin 3 π Z + D 5 A 1 5 cos a Y sin π Z Z + D 6 A 1 5 cos 3 a Y sin π Z + D 7 A 1 2 A 3 cos a Y sin 3 π Z + D 8 A 1 2 A 3 cos a Y sin π Z .
In Equation (57), the coefficients D i , i = 1 , 2 , 8 are the functions of a , E , and q ( D i = D i ( a , E , q ) ). Thus, A 3 is provided by
A 3 = S 1 S 2 ,
where
S 1 = R os q 2 E d 2 3 a 2 A 1
+ E q d 2 3 a 2 π θ 02 ( 2 ) w 13 ( 3 ) + a q v 22 ( 4 ) 4 d 2 + π q w 22 ( 4 ) 4 d 2 + 1 4 π θ 22 ( 4 ) θ 02 ( 4 ) π
× R cs + R os A 1 5 ,
and
S 2 = R os q 2 E d 2 3 a 2 E q d 2 3 a 2 π θ 02 ( 2 ) θ 02 ( 4 ) π R cs + R os A 1 3 .
The amplitude A 3 is determined from the first-, second, third-, fourth-, and fifth-order approximate solutions; w 5 and θ 5 are evaluated from Equations (36) and (56), and are provided by
w 5 = A 5 cos a Y sin π Z + D 9 A 1 3 cos a Y sin 3 π Z + D 10 A 1 5 cos a Y sin 3 π Z + D 11 A 1 2 A 3 cos a Y sin 3 π Z + D 12 A 1 5 cos a Y sin 5 π Z + D 13 A 1 5 cos 3 a Y sin 5 π Z + D 14 A 1 5 cos 3 a Y sin 3 π Z + D 15 A 1 5 cos 3 a Y sin π Z , v 5 = 1 a 2 2 w 5 z y , θ 5 = A 5 q cos a Y sin π Z d 2 + D 16 A 1 3 cos a Y sin 3 π Z + D 17 A 1 5 cos a Y sin 3 π Z + D 18 A 1 2 A 3 cos a Y sin 3 π Z + D 19 A 1 5 cos a Y sin 5 π Z + D 20 A 1 5 cos 3 a Y sin 5 π Z + D 21 A 1 5 cos 3 a Y sin 3 π Z + D 22 A 1 5 cos 3 a Y sin π Z + D 23 A 1 2 A 3 cos a Y sin π Z + D 24 A 1 A 4 sin 2 π Z ,
where the coefficients D i , i = 9 , 10 , 24 are the functions of a , E , a n d q ( D i = D i ( a , E , q ) ).

5.3. Evaluation of Amplitudes A 4 , A 5 , and A 6

Proceeding as above, the sixth-order approximate solutions are obtained from Equations (33) and (36) and are provided by
w 6 = A 6 cos a Y sin π Z + cos 2 a Y ( E 1 sin 2 π Z + E 2 sin 4 π Z + E 3 sin 6 π Z ) + cos 4 a Y ( E 4 sin 2 π Z + E 5 sin 4 π Z + E 6 sin 6 π Z ) , v 6 = 1 a 2 2 w 6 z y , θ 6 = q A 6 π 2 + a 2 cos a Y sin π Z + cos 2 a Y ( E 7 sin 2 π Z + E 8 sin 4 π Z + E 9 sin 6 π Z ) + cos 4 a X ( E 10 sin 2 π Z + E 11 sin 4 π Z + E 12 sin 6 π Z ) + E 13 sin 2 π Z + E 14 sin 4 π Z + E 15 sin 6 π Z .
In Equation (60), the coefficients E i ( i = 1 , 2 , 15 ) are functions of a , E , and q ( E i = E i ( a , E , q ) ). The solution of the sixth-order equation provides A 4 = 0 . The above solution plays an important role in the calculation of A 5 , which can be found by solving the seventh order of Equation (33). The computation becomes more tedious as the order of ϵ increases. Proceeding as above, A 6 = 0 can be obtained; thus, it can be observed that A 2 = A 4 = A 6 = 0 .

6. Natural Convective Heat Transport

The local Nusselt number ( N L ), which represents the heat transfer coefficient, is defined as follows [18,19,20,21,22,23,24,25,26,27,28,29,30,31]:
N L = T n ,
where n stands for dimensionless derivation of the normal direction on the plane. The average Nusselt number N u for heat transfer, which is independent of the axial coordinate Z, can be expressed as follows:
N u = w T ¯ T ¯ Z .
In this context, the overbar indicates that we are using the mean value of the convective heat transfer coefficient across the horizontal plane. Using Equation (62), we obtain N u by integrating the expression over the range of interest, such as on the boundary, with Z = 0 [30].
N u = 1 L 0 L W T T Z | Z = 0 d Y ,
where L denotes the normalized horizontal cell width.

6.1. Local Nusselt Number ( N L )

The influence of the Rayleigh number and Ekman number on heat transfer across a boundary is analyzed through the local Nusselt number. Figure 3a,b displays the local Nusselt number distributions along the heated bottom wall for different Rayleigh and Ekman numbers. In Figure 3a, it is evident that an increase in the Rayleigh number leads to a higher heat transfer rate, which is attributed to improved fluid flow. At lower Rayleigh numbers the convection is minimal and heat transfer is primarily governed by conduction, resulting in the presence of small spikes in the Nusselt number. Figure 3b shows how the local Nusselt number varies along the bottom wall for various Ekman numbers while the Rayleigh number is kept constant at R = 8 R c s . This demonstrates how changes in the Ekman number affect the local heat transfer performance, complementing the findings from varying Rayleigh numbers. This phenomenon can be attributed to the reduced viscous damping effects at lower Ekman numbers, which allow for more vigorous fluid motion and enhanced convective heat transfer. As the Ekman number decreases, the flow becomes less influenced by viscous forces, facilitating stronger thermal gradients and increased heat flux across the fluid layer.

6.2. Average Nusselt Number ( N u )

The average Nusselt number ( N u ) can be calculated to characterize the overall heat transfer behavior in a system. For this analysis, we consider N u ( 2 ) ( s = 1 ) , N u ( 4 ) ( s = 2 ) , and N u ( 6 ) ( s = 3 ) , to represent the second-, fourth-, and sixth-order approximations of the Nusselt number, respectively; N u ( 2 ) is calculated, and is provided by
N u ( 2 ) = 1 + 0.787 1 π R c s · 27 E 2 π 4 2 R cs q E + 1 R cs q E .
Due to the lengthy expression of A 3 and A 5 , the approximations for N u ( 4 ) and N u ( 6 ) are lengthy, and are not included here to conserve space.
Figure 4 illustrates the influence of the Rayleigh number R on the Nusselt number N u for various Ekman numbers (E). From the figure, it can be seen that at low Rayleigh numbers the average Nusselt number is roughly equal to the Rayleigh number; however, when R surpasses 10 4 , the mean Nusselt number shows an increasing trend with rising Rayleigh numbers. Additionally, a transition from conduction to convection heat transfer occurs beyond R = 10 4 . It is important to note that the heat transfer coefficient is influenced not only by the Rayleigh number but also by the Ekman number, with the average Nusselt number rising as the Ekman number decreases. If the Ekman number is low, then the viscous force is weak, while if the Ekman number is large, then the viscous force is strong. Therefore, the combination of the largest Rayleigh number and smallest Ekman number in the selected region provides the maximum value of the average Nusselt number. Figure 5 shows the plot of energy versus R / R c s for several values of E. We find that the potential energy is almost constant for all combinations of the Rayleigh number and Ekman number, while the kinetic energy increases as the value of the Ekman number decreases. This is because a lower Ekman number reduces viscous forces and increases kinetic energy.

7. Distortion of Streamlines and Isotherms

Fluid motion is characterized by a stream function ( Ψ ) derived from the velocity components v and w. The relationship between the stream function Ψ and the velocity components v and w in a two-dimensional flow can be expressed as follows:
v = Ψ Z and w = Ψ Y ,
resulting in a single equation. The arbitrary constants of integration are eliminated by applying the following boundary conditions for Ψ [17]:
Ψ = 2 Ψ = 0 for Y = 0 , 2 π / a and Z = 0 , 1 .
Snapshots of the flow field and heat transfer near the onset of steady-state convection are plotted in the form of streamlines and isotherms, respectively, for various values of R and E. To find the absolute minimum or absolute maximum values of Ψ ( Y , Z ) , all of the critical points are initially calculated; then, each critical point is examined using sufficient conditions to decide whether Ψ has an absolute minimum or absolute maximum value at this point. It is assumed that Ψ is continuous and possesses first- and second-order partial derivatives at a point P ( a , b ) . If P is a critical point, then P ( a , b ) is a point of
absolute   minimum if r τ s 2 > 0 and r > 0 ,
absolute   maximum if r τ s 2 > 0 and r < 0 ,
where
r = Ψ Y Y ( a , b ) , s = Ψ Y Z ( a , b ) and τ = Ψ Z Z ( a , b ) .
No conclusion about an extremum can be drawn if r τ s 2 = 0 , in which case further investigation is needed. If r τ s 2 < 0 , then Ψ has no absolute minimum or absolute maximum at this point, in which case P ( a , b ) is called a saddle point.
A parametric investigation was carried out to examine how the Rayleigh number (R) and Ekman number (E) affect the fluid flow and temperature distribution, as shown in Figure 6, Figure 7 and Figure 8. Figure 6 illustrates the impact of E on the streamlines and isotherms near the critical Rayleigh number ( R c s ). In Figure 6a, the streamlines for E = 0.03 and q = 0.01 reveal that the absolute maximum and minimum values of the circulation strength are 0.4155 and 0.41549 , respectively. Figure 6c presents the streamlines for E = 0.07 and q = 0.01 , showing a reduction in the absolute maximum and minimum circulation intensities to − 0.00998 and − 0.01002 , respectively. Both subfigures in Figure 6a,c indicate that the streamlines exhibit concentric rolling patterns at R R c s . The temperature distribution is represented in Figure 6b,d with isotherms at R R c s . The viscous forces near the onset are comparable to the buoyant forces, and give rise to conduction modes of heat transfer. As E increases from 0.03 to 0.07 (Figure 6b,d), the isotherm remains smooth because the viscous force begins to dominate the buoyant force.
Figure 7 illustrates the impact of the Rayleigh number on the fluid flow and temperature distribution for a fixed Ekman number (=0.03). Figure 7a,b presents the streamlines and isotherms for R = 2 R c s . It is evident from the figure that the primary parallel rolls are slightly inclined to the right, with maximum and minimum circulation strengths of 1.37693 and 1.37692 , respectively. The isotherms exhibit wavy patterns, indicating that buoyant forces are beginning to dominate over viscous forces, marking the onset of convective heat transfer. Figure 7c,d depicts the streamlines and isotherms for R = 5 R c s . In this case, the tilted primary cell has been altered by the emergence of two additional eddies with a circulation strength of 0.823 near the upper left and lower right boundaries. Moreover, two new vortices, referred to as A and A , form within the primary cell, displaying absolute maximum and minimum circulation strengths of 1.920 and 1.920 , respectively. The corresponding isotherm distribution is shown in Figure 7d, which illustrates the convective mode of temperature distribution in the liquid layer. Notably, in Figure 7e,f the primary streamline cell, with maximum and minimum values of 6.20542 and 6.20541 , is observed to have split into two counter-clockwise rotating cells, labeled B and B . This change occurs as the buoyant force becomes significantly stronger relative to the viscous force. Under these conditions, convection-driven fluid motion is present between the two horizontal plates, with most of the temperature being transferred from the bottom plate to the top plate. As the Rayleigh number increases from 2 R c s to 8 R c s , the strength of fluid circulation increases significantly and the temperature gradient between two plates is steep, as can be observed from the streamlines and temperature contours.
Figure 8 highlights the effects of viscous forces on streamlines and isotherms for fixed parameters of R = 8 R c s and q = 0.01 . The influence of the Ekman number (E) is illustrated by a series of flow patterns, shown in Figure 7e,f and Figure 8a–d for E values of 0.03 , 0.05 , and 0.07 . In Figure 8a, the streamlines for E = 0.05 are depicted within the range of 0 Y 1 . The maximum and minimum circulation strengths are measured at 6.07826 and 6.07805 , respectively. This figure reveals two vortices, labeled B and B , located near the upper left and lower right boundaries, each exhibiting a circulation strength of 3.879 . The primary cell is also subdivided into two additional vortices, named A and A . Figure 8c, representing E = 0.07 , shows absolute maximum and minimum circulation strengths of 5.95964 and 5.95926 , respectively. In the range of 0 Y 1 , the flow pattern appears distorted due to the presence of vortices B and B , which are situated near the top and bottom boundaries with a circulation strength of 3.778 . Temperature distributions are shown in Figure 8b,d through isotherms for E = 0.05 and 0.07 , respectively. Overall, Figure 7e,f and Figure 8 indicate a decline in both the circulation strength of the primary cell and the intensity of the temperature distribution as the Ekman number increases. It is noteworthy that the dominance of the left counterclockwise rotating cell over the right clockwise rotating cell becomes more pronounced as the Ekman number decreases. This phenomenon is attributed to the increase in circulation strength of the right cell, which rises from 3.778 at E = 0.07 to 3.879 at E = 0.05 . The enhanced circulation strength in the right cell suggests that reduced viscous effects allow for more vigorous flow patterns, reinforcing the dominance of the left cell as the Ekman number decreases.

Topology of Flow

Topological fluid mechanics focuses on finding the structural characteristics of a fluid flow. This often begins with identifying crucial places in the flow where the vorticity may be zero, then studying the structure of the flow in the vicinity of these spots. The topology constraint of the flow relies on the Euler number, denoted by ζ . According to Jana et al. [32], ζ on the surface is calculated by summing the Poincare indices at ‘stagnation points’ where the fluid is instantaneously at rest. The number of critical points are related by Euler’s identity:
N E ( N H + 1 2 N P ) = ζ ,
where N E refers to the count of elliptic points, N H signifies the count of hyperbolic points, and N P denotes the count of parabolic points [33,34]. This is the most straightforward outcome of a topological character for a flow, and it remains true throughout the fluid flow. Figure 9 illustrates the vorticity contours for different values of R and fixed E = 0.03 and q = 0.01 . The vorticity contours for R R c s are illustrated in Figure 9a with N E = 2 , N H = 2 and N P = 0 . Figure 8b represents the vorticity contours for R = 8 R c s with N E = 8, N H = 8 and N P = 0. From Figure 9a,b, we can conclude that the topological relation specified in Equation (67) is satisfied.

8. Heatfunction

The heatline concept is similar to a streamline, except that it visualizes net energy flow in a convection or conduction heat transfer scenario. The visualization of heat transmission through fluid flow was revolutionized by Kimura and Bejan [18] with the introduction of the idea of heatlines. Building upon this breakthrough, Morega and Bejan [35] effectively applied the concept of heatlines. Several researchers [36,37,38,39,40,41,42,43] have applied this concept to various natural convection systems. The heat function ( H * ) is defined as
H * Y = w T T Z ,
H * Z = v T T Y ,
where T = T s + θ and T s = T o Z , with T 0 as the reference temperature. Differentiating Equations (68) and (69) with respect to Y and Z, respectively, then subtracting the two results yields the following equation:
2 H * Y 2 + 2 H * Z 2 = ( w T ) Y ( v T ) Z .
The boundary conditions on H * can be deduced from the Equations (68) and (69), which define the heat function [30]:
H * ( Y , 0 ) = H * ( 0 , 0 ) + 0 Y w T T Z d Y , at Z = 0 and 0 Y 2 π / a ,
H * ( Y , 1 ) = H * ( 0 , 1 ) + 0 Y w T T Z d Y , at Z = 1 and 0 Y 2 π / a ,
H * ( 0 , Z ) = H * ( 0 , 0 ) 0 Z v T T Z d Z , at Y = 0 and 0 Z 1 ,
H * ( 2 π / a , Z ) = H * ( 2 π / a , 0 ) 0 Z v T T Z d Z , at Y = 2 π / a and 0 Z 1 .

Results and Discussion for Heatlines

Figure 10a–d depicts the net energy flow through convection and conduction at various Rayleigh numbers (R) while keeping the Ekman number constant at ( E = 0.03 ). For R < R c s , conduction primarily governs the heat transfer across the layer. In Figure 10a, we observe the heat transfer when R R c s , with absolute maximum and minimum values of the heatlines recorded at 6.32101 and 0.07152 , respectively, within the range of 0 Y 5 . This indicates that nonlinear heat transfer begins to take place at R R c s , marking a transition from conduction to convection. At this point, the heat transfer across the horizontal layer remains low due to weak buoyancy forces. Figure 10b shows the results for R = 1.05 R c s , where the absolute maximum and minimum values of the heatlines are 6.34066 and 0.2862 , respectively. The influence of increased convection is more evident in the heatlines. As R increases from R c s to 1.05 R c s , the previously curved heatlines transform into closed loops, as illustrated in Figure 10b. With a further increase in the Rayleigh number to R = 1.15 R c s , the absolute maximum and minimum heatline values are 6.36755 and 0.50782 , respectively, as shown in Figure 10c. The maximum value of 6.36755 signifies a further enhancement in heat transfer at R = 1.15 R c s . Finally, Figure 10d presents data for R = 1.15 R c s , displaying absolute maximum and minimum values of 6.38648 and 0.60663 , respectively. Here, convection intensifies and becomes the dominant mechanism for heat transfer. A significant number of heatlines appear near the hot wall, signifying a high heat flux toward the cold wall as R increases. This clearly illustrates the critical role that convection plays in the heat transfer process as the Rayleigh number continues to rise.
Figure 10d and Figure 11a,b demonstrate the impact of the Ekman number E on the heat transfer for a fixed ratio R = 1.25 R c s . In Figure 11a, with E = 0.05 , the heatlines exhibit absolute maximum and minimum values of 5.92779 and 0.66694 , respectively. Conversely, Figure 11b for E = 0.07 shows maximum and minimum values of 4.3253 and 0.65158 , respectively, indicating that both the magnitude and number of closed loops in the heatlines decrease within the range of 0 Y 5 as E increases. Increasing E correlates with a reduction in both the magnitude and number of closed loops in the heatlines. This suggests that higher Ekman numbers, which are associated with more significant viscous forces, diminish the effectiveness of convective heat transfer. Moreover, analysis of Figure 10d and Figure 11a,b indicates that the highest heatline values observed in the corner region increase as E decreases. This further reinforces the idea that lower Ekman numbers are linked to improved heat flow, highlighting the critical role of viscous effects in influencing the efficiency of convective heat transfer in the studied system. This can be attributed to the lower viscosity of the fluid at smaller Ekman numbers, which enhances the convective flow and facilitates more effective heat transfer.

9. Total Entropy Generation

It is essential to grasp how to use energy resources efficiently while minimizing energy degradation, particularly in terms of reducing entropy generation during heat transfer. This study connects entropy generation to the irreversible aspects of heat transfer and the viscous effects that occur within the fluid and at the interfaces between the fluid and solid surfaces. The total nondimensional entropy generation rate, denoted as ( S g e n ), is determined by summing the contributions from heat transfer irreversibility ( S H T I ) and fluid friction irreversibility ( S F F I ). This methodology draws upon the approaches outlined in references [44,45,46,47,48]. The total entropy generation can be expressed mathematically as follows:
S g e n = S H T I + S F F I ,
S H T I = T Y 2 + T Z 2 ,
S F F I = χ 2 v Y 2 + 2 w Z 2 + v Z + w Y 2 .
In this study, the nondimensional velocity and temperature functions, referred to as v , w, and T, are discussed in Section 2, with the magnetic permeability parameter χ (magnetic permeability) fixed as 10 4 . The viscous dissipation model proposed by Baytas [49] serves as the foundation for deriving Equation (75). A parametric investigation has been carried out to assess how the Rayleigh number (R) and Ekmann number (E) influence local entropy generation. Figure 12a–d illustrates the local entropy generation related to heat transfer and fluid friction for various values of R. It can be observed that at lower R values, the maximum entropy generation remains relatively low due to limited heat transfer and fluid motion. This phenomenon occurs because the weak buoyancy forces at low R lead to a regime in which heat transfer is primarily conductive. In contrast, Figure 12b–d reveals that as R increases, the total entropy generation associated with heat transfer also rises, particularly in the vicinity of the hot plate. This observation indicates that enhanced buoyancy forces improve the convective heat transfer mechanism. Additionally, Figure 13a,b demonstrates that as the Ekman number E increases there is a corresponding increase in total entropy generation due to fluid friction, highlighting the relationship between viscous effects and entropy generation within the system. This increase can be attributed to the fact that entropy generation due to fluid friction is dependent on the velocity gradient, which in turn is influenced by the magnetic permeability of the fluid.

10. Conclusions

In this research, we have analytically examined the Rayleigh–Bénard system involving an electrically conducting fluid that rotates around a vertical axis while exposed to a uniform horizontal magnetic field. Our detailed analysis focused on fluid flow patterns, temperature distribution, and heat transfer mechanisms, utilizing streamlines, isotherms, and heatlines for visualization. This approach enhances our understanding of heat flow dynamics in the Earth’s liquid core, particularly when the flow rolls align with the magnetic field. We found that when the Rayleigh number (R) surpasses the critical threshold R > R c s , convection becomes the dominant heat transfer mode, driven by substantial heat flow from the bottom of the layer. Below, we highlight some significant outcomes of our study:
  • Linear stability analysis revealed that as the Ekman number (E) decreases, the critical Rayleigh number ( R c s ) also decreases, indicating that lower E values stabilize the system.
  • The nonlinear partial differential equations were solved using the perturbation method up to O ( ϵ 8 ) yielding approximate solutions for the system.
  • Analysis of the local Nusselt number indicated that its maximum value increases with rising Rayleigh number (R).
  • As E decreases, the number of peaks in the local Nusselt number ( N L ) increases.
  • From the results of the Nusselt number, it is observed that heat flux increases as E decreases.
  • A decrease in E corresponds to an increase in the total energy of the system.
  • Based on the trajectories of heatlines, streamlines, and isotherms, it was determined that decreasing E leads to greater deformation of the flow field and enhanced heat transfer with increased rotation.
  • At higher Rayleigh numbers, the primary source of entropy generation is attributed to irreversible heat transfer; conversely, at high Ekman numbers, irreversibility arises from fluid friction.
To further expand on these findings, we recommend the following avenues for future research:
  • Three-Dimensional Nonlinear Analysis: Future studies should explore three-dimensional nonlinear analyses of the system. While our current analysis is based on two-dimensional models, three-dimensional effects can significantly influence flow patterns, stability, and heat transfer. Investigating the interplay between axial and radial components of the flow may yield richer dynamics and insights into the complexities of the thermal system.
  • Oscillatory Convection: Incorporating oscillatory convection into the study will provide a deeper understanding of how time-dependent flow can affect heat transfer mechanisms. Oscillatory convection can occur due to various external influences, such as periodic heating or fluid motion, and examining its impact could reveal important behaviors in fluid dynamics and thermal transport, particularly in systems subjected to fluctuating thermal gradients.

Author Contributions

Conceptualization, G.S. and Y.R.; Methodology, Y.R.; Software, D.L.; Validation, Y.R., G.S. and D.L.; Visualization G.S.; Writing—Original Draft Preparation, G.S.; Writing—Review and Editing, Y.R., G.S. and D.L.; Funding Acquisition, D.L.; Project Administration, 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 ABF220001.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

The following abbreviations are used in this manuscript:
AAmplitude
B Magnetic field
B s Static magnetic field
B 0 Characteristic field strength
aWavenumber
pPressure
EEkman number
1 Z Unit vector along Z-axis
1 Y Unit vector along Y-axis
qRatio of thermal and magnetic diffusivities
dConvective zone depth
g Gravitational field
H*Heatfunction
N u Average Nusselt number
RRayleigh number
R c s Critical Rayleigh number for stationary convection
TTemperature
T s Static temperature
V s Static velocity
T 0 Reference temperature
Δ T Temperature difference between top and bottom layers
tTime
V Velocity vector
u , v , w Velocity components
X , Y , Z Cartesian coordinates
R B C Rayleigh-Bénard Convection
Greek symbol
Λ Elsasser number
β Adverse temperature gradient
θ Perturbed temperature
η Magnetic diffusivity
ρ Density
ρ 0 Reference density
κ Thermal diffusivity
ν Kinematic viscosity
αThermal expansion coefficient
μDynamic viscosity
μ m Magnetic permeability
Ω Angular velocity
λ Growth rate
Superscript
  Dimensional form
* Perturbed quantities

References

  1. Chandrasekar, S. Hydrodynamic and Hydromagnetic Stability; Oxford Clarendon Press: Oxford, UK, 1961. [Google Scholar]
  2. 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]
  3. Eltayeb, I.A. Hydromagnetic convection in a rapidly rotating fluid layer. Proc. R. Soc. Lond. A 1972, 326, 229–254. [Google Scholar]
  4. Childress, S.; Soward, A.M. Convection-Driven Hydromagnetic Dynamo. Phys. Rev. Lett. 1972, 29, 837. [Google Scholar] [CrossRef]
  5. Soward, A.M. Convection Driven Dynamos. Phys. Earth Planet. Inter. 1979, 20, 134–151. [Google Scholar] [CrossRef]
  6. Braginsky, S.I. Magnetic Waves in the Core of the Earth. Geophys. Asirophys. 1980, 14, 189–208. [Google Scholar] [CrossRef]
  7. Gary, A.; Paul, G.; Roberts, H. A three-dimensional convective dynamo solution with rotating and finitely conducting inner core and mantle. Phys. Earth Planet. Inter. 1995, 91, 63–75. [Google Scholar]
  8. Jones, C.A.; Roberts, P.H. Convection-driven dynamos in a rotating plane layer. J. Fluid Mech. 2000, 404, 311–343. [Google Scholar] [CrossRef]
  9. Taylor, J.B. The Magneto-Hydrodynamics of a Rotating Fluid and the Earth’s Dynamo Problem. Proc. R. Soc. Lond. A 1963, 274, 274–283. [Google Scholar]
  10. Soward, A.M. Finite-amplitude thermal convection and geostrophic flow in a rotating magnetic system. J. Fluid Mech. 1980, 98, 449–471. [Google Scholar] [CrossRef]
  11. Geiger, G.; Busse, F. On the Onset of Thermal Convection in Slowly Rotating Fluid Shells. Geophys. Asrrophys. Fluid Dyn. 1981, 18, 147–156. [Google Scholar] [CrossRef]
  12. Matthews, P.C. Dynamo action in simple convective flows. R. Soc. Lond. A 1999, 455, 1829–1840. [Google Scholar] [CrossRef]
  13. Kono, M.; Roberts, P.H. Definition of the Rayleigh number for geodynamo simulation. Phys. Earth Planet. Inter. 2001, 128, 13–24. [Google Scholar] [CrossRef]
  14. Cattaneo, F.; Hughes, D.W. Dynamo action in a rotating convective layer. J. Fluid Mech. 2006, 553, 401–418. [Google Scholar] [CrossRef]
  15. Roberts, P.H.; Kono, M. Gravitational energy release in an evolving Earth model. Earth Planets Space 2007, 59, 651–659. [Google Scholar] [CrossRef]
  16. Tilgner, A. Transitions in Rapidly Rotating Convection Driven Dynamos. Phys. Rev. Lett. 2012, 109, 248–501. [Google Scholar] [CrossRef]
  17. 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]
  18. Kimura, S.; Bejan, A. The heatline visualization of convective heat transfer. J. Heat Transf. 1983, 105, 916–919. [Google Scholar] [CrossRef]
  19. Rameshwar, Y.; Srinivas, G.; Laroze, D. Nonlinear analysis of cross rolls of electrically conductiong fluid under an applied magnetic field with rotation. Processes 2023, 11, 1945. [Google Scholar] [CrossRef]
  20. Khedher, N.B.; Ullah, Z.; Alturki, M.; Mirza, C.R.; Eldin, S.M. Effect of joul heating and MHD on periodical analysis of current density and amplitude of heat transfer of electrically conducting fluid along thermal magnetized cylinder. Ain Shams Eng. J. 2024, 15, 102–374. [Google Scholar]
  21. Li, S.; Faizan, M.; Ali, F.; Ramasekhar, G.; Muhammad, T.; El-Wahed Khalifa, H.A.; Ahmad, Z. Modelling and analysis of heat transfer in MHD stagnation point flow of Maxwell nanofluid over a porous rotating disk. Alex. Eng. J. 2024, 91, 237–248. [Google Scholar] [CrossRef]
  22. Shamshuddin, M.D.; Anwar, S.; Kanayo Kenneth, A.; Usman; Wasim, J. A semi-analytical approach to investigate the entropy generation in a tangent hyperbolic magnetized hybrid nanofluid flow upon a stretchable rotating disk. J. Magn. Magn. Mater. 2023, 574, 170–664. [Google Scholar] [CrossRef]
  23. Assiri, T.A.; Bilal, M.; Mahmoud, E.E.; Ali, A.; Asamoah, J.K.K.; Adnan. Numerical investigation of forced convective MHD tangent hyperbolic nanofluid flow with heat source/sink across a permeable wedge. AIP Adv. 2024, 14, 65–302. [Google Scholar] [CrossRef]
  24. Ramanuja, M.; Muni Sarala, G.; Kavitha, J.; Srinivasulu, A.; Gopi Krishna, G. A fully developed viscous electrically conducting fluid through infinitely parallel porous plates. Heat Transf. 2024; early view. [Google Scholar] [CrossRef]
  25. Tang, Q.; Sun, Y.; Li, J.; Zhou, M.; Yang, D.; Pang, Y. Preparations of 25 wt% of Pyraclostrobin Nanosuspension Concentrate (SC) Using Lignosulfonate-Based Colloidal Spheres to Improve Its Thermal Storage Stability. Molecules 2024, 29, 1419. [Google Scholar] [CrossRef] [PubMed]
  26. Li, J.-X.; Lu, Y.-J.; Quan, K.-Y.; Wu, L.-B.; Feng, X.; Wang, W.-Z. One-pot cocrystallization of mononuclear and 1D cobalt(II) complexes based on flexible triclopyr and 2, 2 bipyridine coligands: Structural analyses, conformation comparison, non-covalent interactions and magnetic properties. J. Mol. Struct. 2024, 1297, 136–830. [Google Scholar] [CrossRef]
  27. Wu, Y.-L.; Tang, P.-F.; Zhang, Q.; Yan, Y.-T.; Zhang, S.; Yang, G.-P.; Wang, Y.-Y. Nanoarchitectonics and catalytic performances of metal–organic frameworks supported metal nanoparticles. Appl. Organomet. Chem. 2023, 38, e7308. [Google Scholar] [CrossRef]
  28. 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]
  29. 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]
  30. 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]
  31. Batchelor, G.K. An Introduction to Fluid Dynamics; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  32. 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]
  33. 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]
  34. 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]
  35. 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]
  36. Bejan, A. Convection Heat Transfer; Wiley: New York, NY, USA, 1984; pp. 21–23. [Google Scholar]
  37. Morega, A.I.M. Magnetic field influence on the convective heat transfer in the solidification processes. Rev. Roum. Sci. Tech. Electrotechnol. Energy 1988, 33, 155–166. [Google Scholar]
  38. Nield, D.A.; Bejan, A. Convection in Porous Media, 2nd ed.; Springer: New York, NY, USA, 1999. [Google Scholar]
  39. Costa, V.A.F. Bejan’s heatlines and masslines for convection visualization and analysis. Appl. Mech. Rev. 2006, 59, 26–145. [Google Scholar] [CrossRef]
  40. Basak, T.; Roy, S. Role of Bejan’s heatlines in heat flow visualization and optimal thermal mixing for differentially heated square enclosures. Int. J. Heat Mass Transfer. 2008, 51, 3486–3503. [Google Scholar] [CrossRef]
  41. Basak, T.; Roy, S.; Matta, A.; Pop, I. Analysis of heatlines for natural convection within porous trapozoidal cnclosure: Effect of uniform and non-uniform heating of bottom wall. Int. J. Heat Mass Transfer. 2010, 53, 5947–5961. [Google Scholar] [CrossRef]
  42. Deng, Q.H.; Tang, G.F.; Li, Y. A combined temperature scale for analyzing natural convection in rectangular enclosures with discrete wall heat sources. Int. J. Heat Mass Transfer. 2002, 45, 3437–3446. [Google Scholar] [CrossRef]
  43. Varol, Y.; Oztop, H.F.; Mobedi, M.; Pop, I. Visualization of heat flow using Bejan’s heatline due to natural convection of water near 4 °C in thick walled porous cavity. Int. J. Heat Mass Transfer. 2010, 53, 1691–1698. [Google Scholar] [CrossRef]
  44. Bejan, A. Entropy Generation through Heat and Fluid Flow; Wiley: New York, NY, USA, 1982. [Google Scholar]
  45. Basak, T.; Anandalakshmi, R.; Roy, S.; Pop, I. Role of entropy generation on thermal management due to thermal convection in porous trapezoidal enclosures with isothermal and non-isothermal heating of wall. Int. J. Heat Mass Trans. 2013, 67, 810–828. [Google Scholar] [CrossRef]
  46. Basak, T.; Das, D.; Biswal, P. Heatlines: Modeling, visualization, mixing and thermal management. Prog. Energy Combust. Sci. 2017, 64, 157–218. [Google Scholar] [CrossRef]
  47. Al-Hadhrami, A.K.; Elliott, L.; Ingham, D.B. A New Model for Viscous Dissipation in Porous Media across a Range of Permeability Values. Trans. Porous Med. 2003, 53, 117–122. [Google Scholar] [CrossRef]
  48. Homan, K.; Gurgenci, H. Effects of viscous dissipation and boundary conditions on forced convection in a channel occupied by a saturated porous medium. Trans. Porous Med. 2007, 68, 301–319. [Google Scholar] [CrossRef]
  49. Baytas, A.C. Entropy generation for natural convection in an inclined porous cavity. Int. J. Heat Mass Trans. 2000, 43, 2089–2099. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the physical system.
Figure 1. Schematic diagram of the physical system.
Processes 12 01882 g001
Figure 2. Dependence of R on a for E = 0.03 , 0.05 and 0.07 .
Figure 2. Dependence of R on a for E = 0.03 , 0.05 and 0.07 .
Processes 12 01882 g002
Figure 3. Variation of N L with respect to Y: (a) E = 0.03, R c s = 5262.664915 , and q = 0.01 for different R; (b) R = 8 R c s and q = 0.01 for different E.
Figure 3. Variation of N L with respect to Y: (a) E = 0.03, R c s = 5262.664915 , and q = 0.01 for different R; (b) R = 8 R c s and q = 0.01 for different E.
Processes 12 01882 g003
Figure 4. Dependence of N u on R for E = 0.03 , 0.05 , 0.07 , and q = 0.01 .
Figure 4. Dependence of N u on R for E = 0.03 , 0.05 , 0.07 , and q = 0.01 .
Processes 12 01882 g004
Figure 5. Kinetic energy and potential energy against R / R c s for q = 0.01 .
Figure 5. Kinetic energy and potential energy against R / R c s for q = 0.01 .
Processes 12 01882 g005
Figure 6. The effect of E near R c s and q = 0.01 : streamlines are plotted for (a) E = 0.03 and (c) E = 0.07 and isotherms for (b) E = 0.03, R c s = 5262.664915 and (d) E = 0.07, R c s = 6755.126806 .
Figure 6. The effect of E near R c s and q = 0.01 : streamlines are plotted for (a) E = 0.03 and (c) E = 0.07 and isotherms for (b) E = 0.03, R c s = 5262.664915 and (d) E = 0.07, R c s = 6755.126806 .
Processes 12 01882 g006
Figure 7. The effect of R = 2 R c s , 5 R c s , 8 R c s : streamlines (a,c,e) and isotherms (b,d,f) are plotted for E = 0.03, R c s = 5262.664915 , and q = 0.01.
Figure 7. The effect of R = 2 R c s , 5 R c s , 8 R c s : streamlines (a,c,e) and isotherms (b,d,f) are plotted for E = 0.03, R c s = 5262.664915 , and q = 0.01.
Processes 12 01882 g007aProcesses 12 01882 g007b
Figure 8. The effect of E = 0.05, R c s = 5896.136279 , and E = 0.07, R c s = 6755.126806 : streamlines (a,c) and isotherms (b,d) are plotted for R = 8 R c s and q = 0.01 .
Figure 8. The effect of E = 0.05, R c s = 5896.136279 , and E = 0.07, R c s = 6755.126806 : streamlines (a,c) and isotherms (b,d) are plotted for R = 8 R c s and q = 0.01 .
Processes 12 01882 g008
Figure 9. Vorticity lines for E = 0.03, R c s = 5262.664915 , and q = 0.01 : (a) R R c s and (b) R = 8 R c s .
Figure 9. Vorticity lines for E = 0.03, R c s = 5262.664915 , and q = 0.01 : (a) R R c s and (b) R = 8 R c s .
Processes 12 01882 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.03 and R c s = 5262.664915 , 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.03 and R c s = 5262.664915 , T 0 = 1, and q = 0.01, respectively.
Processes 12 01882 g010
Figure 11. The effect of E = 0.05, R c s = 5896.136279 , E = 0.07, and R c s = 6755.126806 : heatlines (a,b) are plotted for R = 1.25 R c s , T 0 = 1, and q = 0.01 .
Figure 11. The effect of E = 0.05, R c s = 5896.136279 , E = 0.07, and R c s = 6755.126806 : heatlines (a,b) are plotted for R = 1.25 R c s , T 0 = 1, and q = 0.01 .
Processes 12 01882 g011
Figure 12. Total entropy generation due to heat transfer and fluid flow for E = 0.03. R c s = 5262.664915 , and q = 0.01 : ( a ) R R c s , (b) R = 2 R c s , (c) 5 R c s , and (d) 8 R c s .
Figure 12. Total entropy generation due to heat transfer and fluid flow for E = 0.03. R c s = 5262.664915 , and q = 0.01 : ( a ) R R c s , (b) R = 2 R c s , (c) 5 R c s , and (d) 8 R c s .
Processes 12 01882 g012
Figure 13. Total entropy generation due to heat transfer and fluid flow for R = 8 R c s and q = 0.01 : (a) E = 0.05, R c s = 5896.136279 and (b) E = 0.07, R c s = 6755.126806 .
Figure 13. Total entropy generation due to heat transfer and fluid flow for R = 8 R c s and q = 0.01 : (a) E = 0.05, R c s = 5896.136279 and (b) E = 0.07, R c s = 6755.126806 .
Processes 12 01882 g013
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

Srinivas, G.; Rameshwar, Y.; Laroze, D. Influence of Rotation and Viscosity on Parallel Rolls of Electrically Conducting Fluid. Processes 2024, 12, 1882. https://doi.org/10.3390/pr12091882

AMA Style

Srinivas G, Rameshwar Y, Laroze D. Influence of Rotation and Viscosity on Parallel Rolls of Electrically Conducting Fluid. Processes. 2024; 12(9):1882. https://doi.org/10.3390/pr12091882

Chicago/Turabian Style

Srinivas, G., Y. Rameshwar, and D. Laroze. 2024. "Influence of Rotation and Viscosity on Parallel Rolls of Electrically Conducting Fluid" Processes 12, no. 9: 1882. https://doi.org/10.3390/pr12091882

APA Style

Srinivas, G., Rameshwar, Y., & Laroze, D. (2024). Influence of Rotation and Viscosity on Parallel Rolls of Electrically Conducting Fluid. Processes, 12(9), 1882. https://doi.org/10.3390/pr12091882

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