# Effects of Thermal Radiation, Heat Generation, and Induced Magnetic Field on Hydromagnetic Free Convection Flow of Couple Stress Fluid in an Isoflux-Isothermal Vertical Channel

Authors: Hasan Nihal Zaidi; Mohammed Yousif; S. Nazia Nasreen
Journal: Journal of Applied Mathematics
 (2020)

---

## Abstract

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.


---

## Body

## 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 distanceh. 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.
(1)
The flow is steady and fully developed(2)
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(3)
The flux of thermal radiation in they-direction is negligible compared to the x-directionFigure 1
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 fluxqr 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.Introducing dimensionless parameters(8)u=u′u0,y=y′h,θ=T−T0ΔT,M=B0′hvμeρ,B=Bx′u0μeρ,Rt=T2−T1ΔT,u0=gβΔTh2μ.Equations (1), (2), and (6) become
(9)σ2d4udy4−d2udy2−MdBdy−θ=0,(10)d2Bdy2+MPmdudy=0,(11)1+4R3d2θdy2+ϕθ=0,whereR=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,wherea1,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 numberNu 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.Figure 2
Velocity profiles for different values ofM and ϕ.Figure 3
Velocity profiles for different values ofR and ϕ.Figure 4
Velocity profile for different values ofσ and Pm.Figure 5
Induced magnetic field for different values ofϕ and R.Figure 6
Induced magnetic field for different values ofM and ϕ.Figure 7
Induced magnetic field for different values of the magnetic Prandtl numberPm.Figure 8
Temperature profiles for different values ofR at ϕ=2.0.Figure 9
Temperature profiles for different values ofϕ at R=0.5.Figure 10
Induced current density for different values of radiation parameterR.Figure 11
Induced current density for different values ofM and ϕ.Figure 12
Induced current density for different values of the magnetic Prandtl numberPm.Figure 13
Variation of Nusselt numberNu1 versus ϕ for different values of R.Figure 14
Variation of skin frictions0 versus ϕ for different values of M and R.Figure 15
Variation of skin frictions0 versus ϕ for different values of Pm and σ.Figure 16
Variation of skin frictions1 versus ϕ for different values of M and R.Figure 17
Variation of skins1 versus ϕ for different values of Pm and σ.Figure2 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.Figures5–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 parameterR 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.Figures10–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.Figure13 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:(1)
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(2)
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(3)
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(4)
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(5)
It is also found that at the isoflux wall skin friction varies inversely with the parametersR 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



---

*Source: 4539531-2020-09-07.xml*