The study scrutinizes the effects of thermal radiation, heat generation, and induced magnetic field on steady, fully developed hydromagnetic free convection flow of an incompressible viscous and electrically conducting couple stress fluid in a vertical channel. The channel walls are maintained at an isoflux-isothermal condition, such that the left channel wall is maintained at a constant heat flux. In contrast, the right channel wall is maintained at a constant temperature. The governing simultaneous equations are solved analytically utilizing the method of undetermined coefficient, and closed form solutions in dimensionless form have been acquired for the velocity field, the induced magnetic field, and the temperature field. The expression for the induced current density has been also obtained. A parametric study for the velocity, temperature, and induced magnetic field profiles, as well as for the skin-friction coefficient, Nusselt number, and induced current density, is conducted and discussed graphically.
1. Introduction
The problem of hydromagnetic free convection and heat transfer flow through a vertical channel has drawn much attention due to its possible application in many engineering and technological activities. A magnetic field induced to the fluid flow generates an electromotive force and changes the velocity distribution. The behavior of flow firmly relies upon the orientation and strength of the applied magnetic field. The applied magnetic field manipulates the suspended particles and modifies their concentration in the fluid, which unequivocally changes heat transfer characteristics of the flow. Such flows are found in the fusion reactors, MHD generator, MHD pumps, MHD bearings, etc. The effect of a magnetic field on free convection heat transfer was discussed by Sparrow and Cess [1]. They investigated the simultaneous action of buoyancy and induced magnetic field on free convection heat transfer. In the presence of electrical or magnetic fields, Poots [2] considered the natural convection flow of an electrically conducting viscous fluid, such as mercury or liquid sodium. He examined natural convection flow between two long parallel plane surfaces with a uniform magnetic field applied normal to the surfaces, with and without heat sources, and natural convection flow through a horizontal circular tube. Ghosh et al. [3] have investigated hydromagnetic natural convection boundary layer flow in the presence of an induced magnetic field on an infinite vertical plate and obtained an exact solution. Singh et al. [4] numerically studied free convection flow in a vertical channel under the influence of an applied magnetic field, including the effect of magnetic induction. Magnetohydrodynamic free convection flow between two vertical parallel porous plates was analyzed by Sarvesanand and Singh [5] considering the effect of the induced magnetic field emanating from the motion of an electrically conducting fluid. MHD convection fluid flow and heat transfer in an inclined microchannel with heat generation were studied by Zaidi and Ahmad [6]. Kumar and Singh [7] have discussed unsteady MHD free convection flow past a semi-infinite vertical wall with the induced magnetic field and solved coupled nonlinear partial differential equations by the Crank-Nicolson, finite difference method. Jha and Aina [8, 9] analyzed magnetohydrodynamic natural convection flow in a vertical porous or nonporous microchannel in the presence of an induced magnetic field. Recently, Mohammad et al. [10] investigated magnetohydrodynamic- (MHD-) radiated nanomaterial viscous material flow through a curved surface with second-order slip and entropy generation considering the effects of viscous dissipation, thermal radiation, and Joule heating. In the modeling of flow problem, they utilized the Buongiorno model thermophoresis and Brownian diffusion.
The thermal radiation effect on the free convection flow is important in the context of space technology and many engineering applications, such as in advanced types of power plants for nuclear rockets, reentry vehicles, high-speed flights, and procedures involving high temperatures. The Rosseland diffusion approximation is valid for optically thick media. The radiative heat fluxes approximated for an optically dense medium by the Rosseland diffusion approximation have been used extensively in many studies related to radiation. Cess [11] examined the effect of radiation on boundary-layer development in free convection of an absorbing-emitting gas on a vertical plate, and he presented interactions between radiation and free convection. The combined radiative and convective heat transfer in a vertical channel has been studied by Bouali and Mezrhab [12] and Grosan and Pop [13]. Kaladhar et al. [14] have presented mixed convection flow of couple stress fluid in a vertical channel with radiation and Soret effect employing a quasilinearization method. The entropy generation in the flow of ferromagnetic nanofluid with nonlinear radiation and slip condition on a nonlinear stretching sheet is explored by Rashid et al. [15]. Recently, Gireesha et al. [16] have studied the impact of thermal radiation and natural convection on nanoparticle-embedded water-based hybrid nanofluid flow across a permeable longitudinal fin using Darcy’s model.
Studies related to the free convection channel flow and heat transfer characteristic of couple stress fluids find several applications in many industrial processes such as the extrusion of polymer fluids, solidification of liquid crystals, and cooling of metallic plates in bath colloidal solution. To describe the flow of microcontinuum fluids (complex fluids), various general theories as well as specific methods were developed from time to time. The key feature of couple stresses is the implementation of a size-dependent effect. The size effect of material particles within the continua is neglected by classical continuum mechanics. This is consistent with ignoring the rotational interaction between the fluid particles, which results in symmetry of the force-stress tensor. Nevertheless, this cannot be valid in some important cases such as fluid flow with suspended particles, and a size-dependent couple stress theory is required. Stokes [17, 18] introduced the theory of couple stress fluids. He presented the straightforward generalization of the classical theory that allows polar effects such as the presence of couple stresses and body couples. The couple stress fluids also have the mechanism to describe rheological complex fluid-like liquid crystal and human blood. The rheological properties of such fluids are important in the flow of nuclear fuel slurries, lubricants with heavy oil and greases, paper coating, and polymers. Soundalgekar and Aranake [19] analyzed the fully developed MHD couple stress fluid in a channel and observed that the couple stresses are more effective in the presence of a weak magnetic field. Free convection of a couple stress fluid in a vertical channel was considered by Umavathi et al. [20]. Kaladhar [21] has obtained a solution using a homotopy analysis method for the natural convection of a couple stress fluids between two vertical plates, taking into account the Hall and Joule heating effects. Numerically, Makinde and Eegunjobi [22] contemplated the combined effects of nonlinear thermal radiation, buoyancy forces, thermophoresis, and Brownian motion on the entropy generation rate in hydromagnetic couple stress nanofluid flow through a vertical channel utilizing a shooting technique combined with the Runge–Kutta–Fehlberg method.
To the best of the authors’ knowledge, the effects of thermal radiation, heat generation, and an induced magnetic field on hydromagnetic free convection flow of couple stress fluid in an isoflux-isothermal vertical channel have not been reported in the literature. The governing equations for the velocity field, induced magnetic field, and temperature field were solved analytically. A parametric study was conducted, and the results were portrayed and discussed with the help of graphs.
2. Mathematical Formulation
Consider a laminar free convection flow of an electrically conducting couple stress fluid in a vertical parallel plate channel separated by a distance h. The walls of the channel were maintained at the isoflux-isothermal condition. The physical configuration of the problem is shown in Figure 1. We choose a Cartesian coordinate system such that the x′-axis is taken vertically upward along the plates and the y′-axis is normal to the channel. A uniform magnetic field of strength B0′ is imposed in the direction normal to the flow direction. For a fluid with significant electrical conductivity σe, this in turn induces a magnetic field Bx′ along the x′-axis. Let u′ be the velocity of the fluid along the x′-axis, then the velocity u′y′,0,0 has only a vertical component and is a function of y′ only. In contrast, magnetic field Bx′,B0′,0 has a component in the x′- and y′-directions, respectively. The following assumptions are adopted to facilitate the solution of the governing equations.
The flow is steady and fully developed
The fluid physical properties are constant, except the density in the body force term in the momentum equation for which the Boussinesq approximation is invoked
The flux of thermal radiation in the y-direction is negligible compared to the x-direction
Physical configuration and coordinate system.
Under these assumptions, the governing equations of the system are [7–9, 22]
(1)νd2u′dy′2−ηρd4u′dy′4+μeB0′ρdBx′dy′+gβT−T0=0,(2)1σeμed2Bx′dy′2+B0′du′dy′=0,(3)kρcpd2Tdy′2+Q0ρcpT−T0−dqrdy′=0.
Using the Rosseland approximation, the radiative heat flux qr is given by
(4)qr=−4σ∗3k∗∂T4∂y,
where σ∗ and k∗ are the Stefan Boltzmann constant and Rosseland mean absorption coefficient, respectively. Consider that the temperature difference within the flow is sufficiently small such that T4 may be expressed as a linear function of temperature.
(5)T4≈4Th3T−3Th4.
Using Equations (4) and (5), Equation (3) becomes
(6)kρcp+16σ∗Th33k∗d2Tdy′2+Q0ρcpT−T0=0.
The associated boundary conditions are
(7a)u′=0,d2u′dy′2=0,Bx′=0,q1=−kdTdy′aty′=0,(7b)u′=0,d2u′dy′2=0,Bx′=0,T=T2aty′=h.
Equations (1), (2), and (6) become
(9)σ2d4udy4−d2udy2−MdBdy−θ=0,(10)d2Bdy2+MPmdudy=0,(11)1+4R3d2θdy2+ϕθ=0,
where R=16σ∗Th3/3k∗ is the radiation parameter, σ=1/hη/μ is the couple stress parameter, ϕ=Q0h2/k is heat generation coefficient, and Pm=νσμe is the magnetic Prandtl number.
The relevant dimensionless boundary conditions are
(12a)u=0,d2udy2=0,B=0,dθdy=−1aty=0,(12b)u=0,d2udy2=0,B=0,θ=Rtaty=1.
Solutions of Equations (9)–(11) are subjected to the boundary conditions (12a) and (12b) given by
(13)θ=c1sina1y+c2cosa1y,(14)u=c3+c4cosha2y+c5sinha2y+c6cosha3y+c7sinha3y+1a4c1sina1y+c2cosa1y,(15)B=c8−MPmc4a2sinha2y+c5a2cosha2y+c6a3sinha3y+c7a3cosha3y+MPma1a4c1sina1y+c2cosa1y,
The induced current density is given by
(16)J=−dBdy=MPmc4cosha2y+cs5inha2y+c6cosha3y+c7sinha3y+MPma4c1sina1y+c2cosa1y,
where a1,a2,a3,a4,c1,c2,c3,c4,c5,c6,c7,andc8 are constants and are given in the appendix.
The skin friction on both channel walls in dimensionless form is
(17a)s0=dudyy=0=c5a2+c7a3+c1a1a4,(17b)s1=dudyy=1=c4a2sinha2+c5a2cosha2+c6a3sinha3+c7a3cosha3+a1a4c1cosa1−c2sina1.
And the Nusselt number Nu in dimensionless form is
(18a)Nu0=−dθdyy=0=c1a1=−1,(18b)Nu1=−dθdyy=1=−a1c1cosa1−c2sina1.
3. Results and Discussion
In the present study, the effect of thermal radiation, heat generation, and an induced magnetic field on hydromagnetic free convection flow of couple stress fluid in an isoflux-isothermal vertical channel is analyzed. The governing equations solved analytically, the obtained closed form solution for Equations (9)–(11), and the computed results of the analytical solutions are presented in Figures 2–17.
Velocity profiles for different values of M and ϕ.
Velocity profiles for different values of R and ϕ.
Velocity profile for different values of σ and Pm.
Induced magnetic field for different values of ϕ and R.
Induced magnetic field for different values of M and ϕ.
Induced magnetic field for different values of the magnetic Prandtl number Pm.
Temperature profiles for different values of R at ϕ=2.0.
Temperature profiles for different values of ϕ at R=0.5.
Induced current density for different values of radiation parameter R.
Induced current density for different values of M and ϕ.
Induced current density for different values of the magnetic Prandtl number Pm.
Variation of Nusselt number Nu1 versus ϕ for different values of R.
Variation of skin friction s0 versus ϕ for different values of M and R.
Variation of skin friction s0 versus ϕ for different values of Pm and σ.
Variation of skin friction s1 versus ϕ for different values of M and R.
Variation of skin s1 versus ϕ for different values of Pm and σ.
Figure 2 shows the variation of velocity for different values of Hartmann number M and the heat generation coefficient ϕ. It is clear that the velocity of fluid particles decreases on the increase of the Hartmann number M because as the magnetic field intensity increases, the resisting Lorentz force will increase; consequently, the velocity profile decreases. Meanwhile, in the same figure, a rise in heat generation coefficient ϕ enhances the velocity profile. Figure 3 illustrates the effect of the radiation parameter R on the velocity profile. The velocity of fluid particles decreases with the increase in radiation parameter R. The maximum velocity occurs at y=0.4 for R=1, and it shifts towards the wall y=0.0. The effect of the magnetic Prandtl number Pm and the couple stress parameter σ on the velocity profile is displayed in Figure 4. It is evident from Figure 4 that an increase in the value of the couple stress parameter σ enhances the velocity of fluid particles while the velocity of the fluid particles decreases with increase of the magnetic Prandtl number Pm. It can also be noted that whenσ→0, the velocity becomes large; this shows that the velocity in the case of couple stress fluid is less than that in the Newtonian fluid case.
Figures 5–7 describe the behavior of the induced magnetic field with the various parameters. Figure 5 shows the effect of radiation parameter R and heat generation coefficient ϕ on the induced magnetic field profile. It is observed that in the first one-third of the channel, as the value of radiation parameterRincreases, the induced magnetic field decreases and the reverse is the case in the remaining part of the channel; conversely, as the value ofϕincreases, the induced magnetic field increases fromy=0.0toy=0.4and reverses its direction fromy=0.4toy=1.0. Figures 6 and 7 display the effect of the Hartmann number M and the magnetic Prandtl number Pm on the induced magnetic field, respectively. It is noted that the induced magnetic field is directly proportional to the Hartmann number as well as the magnetic Prandtl number near the channel wall at y=0.0 while it is inversely proportional near the channel wall at y=1.0.
The effect of radiation parameter R and heat generation coefficient ϕ on the temperature profile is illustrated in Figures 8 and 9, respectively. It was found that temperature decreases with the increase of the radiation parameter; meanwhile, a rise in the heat generation coefficient increases the temperature of the fluid.
Figures 10–12 depicts the behavior of induced current density with the various parameters occurring in the governing equations. Figure 10 exhibits the variation of induced current density concerning radiation parameter R. It is clear from the graph that the value of induced current density decreases with the increase in radiation parameter R in the central part of the channel, while the reverse is the case near the channel walls. The influence of the Hartmann number M and the heat generation coefficient ϕ on the induced current density is illustrated in Figure 11. The graph shows that the Hartmann number M reduces the induced current density and the heat generation coefficient ϕ enhances the induced current density in the center of the channel, while the converse is the case near the walls. The impact of the magnetic Prandtl number Pm on the induced current density is shown in Figure 12. It is observed that the induced current density increases on the increase of the magnetic Prandtl number in the central region. It is also clear from Figures 10–12 that there exist two points of intersection inside the vertical channel where the induced current density is independent of the radiation parameter, Hartmann number, and Prandtl number.
Figure 13 demonstrates the effects of the radiation parameter R and heat generation coefficient ϕ on the Nusselt number. It was noted that the heat flux at the right wall decreases with the increase of the radiation parameter R while the heat flux at the left wall is constant. The effect of different parameters on skin friction is illustrated in Figures 14–17. It was observed that the magnitude of the skin friction on the isoflux channel wall decreases with the increase of the Hartmann number M, magnetic Prandtl number Pm, and radiation parameter R, whereas an increase in the value of the couple stress parameter σ increases the skin friction (see Figures 14 and 15). The reverse is the case for the isothermal wall (see Figures 16 and 17).
4. Conclusion
In the current analysis, the effect of thermal radiation, heat generation, and induced magnetic field on hydromagnetic free convection flow of couple stress fluid in an isoflux- isothermal vertical channel has been investigated. The governing equations are solved analytically, and closed form solutions are obtained. The outcomes of the present study are summarized as follows:
The velocity of fluid particles decreases with the increase in the radiation parameter, Hartmann number, and magnetic Prandtl number. In contrast, the velocity profiles increase with an increase in the heat generation coefficient and couple stress parameter
The induced magnetic field increases with the increase of the magnetic Prandtl number, Hartmann number, and heat generation coefficient and decreases with an increase of the radiation parameter within almost the left one-third of the channel, while the trend reverses itself near the isothermally heated wall
It was observed that the increase in the magnetic Prandtl number and heat generation coefficient enhances the induced current density. In contrast, the radiation parameter and Hartmann number suppress the induced current density. The observations are for the middle part of the channel, whereas the opposite cases were observed near both walls
Temperature and rate of heat transfer increase with the increase in the heat generation parameter and have a decreasing tendency with the increase in radiation parameter
It is also found that at the isoflux wall skin friction varies inversely with the parameters R and M on the increase in heat generation coefficient while the reverse is the case at the isothermal wall. The effect of the magnetic field on the skin friction is useful in mechanical engineering or modeling a system. We can easily optimize the value of skin friction by obtaining a suitable value of the Hartmann number
The data used to support the finding of this study is included with in the article.
Conflicts of Interest
The authors declare that they have no competing interests.
SparrowE. M.CessR. D.The effect of a magnetic field on free convection heat transfer19613426727410.1016/0017-9310(61)90042-42-s2.0-50549188123PootsG.Laminar natural convection flow in magneto-hydrodynamics19613112510.1016/0017-9310(61)90002-32-s2.0-4644240813GhoshS. K.Anwar BégO.ZuecoJ.Hydromagnetic free convection flow with induced magnetic field effects201045217518510.1007/s11012-009-9235-x2-s2.0-77952240074SinghR. K.SinghA. K.SachetiN. C.ChandranP.On hydromagnetic free convection in the presence of induced magnetic field201046552352910.1007/s00231-010-0594-62-s2.0-77954564463SarvesanandSinghA. K.Magnetohydrodynamic free convection between vertical parallel porous plates in the presence of induced magnetic field20154133310.1186/s40064-015-1097-12-s2.0-84937389257ZaidiH. N.AhmadN.MHD convection fluid flow and heat transfer in an inclined microchannel with heat generation20175512413110.11648/j.ajam.20170505.11KumarA.SinghA. K.Unsteady MHD free convective flow past a semi-infinite vertical wall with induced magnetic field201322246247110.1016/j.amc.2013.07.0442-s2.0-84883077141JhaB. K.AinaB.Role of induced magnetic field on MHD natural convection flow in vertical microchannel formed by two electrically non-conducting infinite vertical parallel plates20165532087209710.1016/j.aej.2016.06.0302-s2.0-84978880629JhaB. K.AinaB.Magnetohydrodynamic natural convection flow in a vertical micro-porous-channel in the presence of induced magnetic field201864143410.1016/j.cnsns.2018.04.0042-s2.0-85045616894MohammadR.KhanM. I.KhanN. B.JameelM.Magnetohydrodynamics (MHD) radiated nanomaterial viscous material flow by a curved surface with second order slip and entropy generation2020189, article 10529410.1016/j.cmpb.2019.105294CessR. D.The interaction of thermal radiation with free convection heat transfer19669111269127710.1016/0017-9310(66)90119-02-s2.0-0001064252BoualiH.MezrhabA.Combined radiative and convective heat transfer in a divided channel20061618410610.1108/096155306106369732-s2.0-30344463271GrosanT.PopI.Thermal radiation effect on fully develop mixed convection flow in a vertical channel20072713747KaladharK.MotsaS. S.SrinivasacharyaD.Mixed convection flow of couple stress fluid in a vertical channel with radiation and soret effects2016914350RashidM.KhanM. I.HayatT.KhanM. I.AlsaediA.Entropy generation in flow of ferromagnetic liquid with nonlinear radiation and slip condition201927644145210.1016/j.molliq.2018.11.1482-s2.0-85058025059GireeshaB. J.SowmyaG.KhanM. I.ÖztopH. F.Flow of hybrid nanofluid across a permeable longitudinal moving fin along with thermal radiation and natural convection202018510516610516610.1016/j.cmpb.2019.105166StokesV. K.Couple stresses in fluids1966991709171510.1063/1.17619252-s2.0-36849106425StokesV. K.Effects of couple stresses in fluids on hydromagnetic channel flows19681151131113310.1063/1.16920562-s2.0-22244447615SoundalgekarV. M.AranakeR. N.Effects of couple stresses in MHD channel flow197744330130810.1016/0029-5493(77)90165-02-s2.0-0017719784UmavathiJ. C.MalashettyM. S.ChamkhaA. J.Al MudhafA.2004KuwaitInternational Mechanical Engineering conferenceKaladharK.Natural convection flow of couple stress fluid in a vertical channel with hall and joule heating effects20151271071107810.1016/j.proeng.2015.11.4652-s2.0-84971454035MakindeO. D.EegunjobiA. S.MHD couple stress nanofluid flow in a permeable wall channel with entropy generation and nonlinear radiative heat201712210.1299/jtst.2017jtst00332-s2.0-85040741287