Second law analysis for hydromagnetic couple stress fluid flow through a porous channel

Abstract

In this work, the combined effects of magnetic field and ohmic heating on the entropy generation rate in the flow of couple stress fluid through a porous channel are investigated. The equations governing the fluid flow are formulated, non-dimensionalised and solved using a rapidly convergent semi-analytical Adomian decomposition method (ADM). The result of the computation shows a significant dependence of fluid’s thermophysical parameters on Joule’s dissipation as well as decline in the rate of change of fluid momentum due to the interplay between Lorentz and viscous forces. Moreover, the rate of entropy generation in the flow system drops as the magnitude of the magnetic field increases.

Keywords

Magnetic field; Entropy generation; Slip flow; Irreversibility ratio; ADM; Ohmic heating

Nomenclature

${\textstyle {\beta }_{1}}$- dimensionless Navier slip parameter at the lower wall

${\textstyle {\beta }_{2}}$- dimensionless Navier slip parameter at the upper wall

${\textstyle \eta }$- coefficient of couple stress

${\textstyle {\gamma }_{1}}$- Navier slip coefficient at the lower plate

${\textstyle {\gamma }_{2}}$- Navier slip coefficient at the upper plate

${\textstyle \kappa }$- thermal conductivity

${\textstyle \mu }$- dynamic viscosity

${\textstyle \Omega }$- parameter that measures the temperature difference between the two heat reservoirs

${\textstyle \rho }$- fluid density

${\textstyle {\sigma }_{e}}$- electrical conductivity

${\textstyle \theta }$- dimensionless fluid temperature

${\textstyle a^{2}}$- couple stress parameter

${\textstyle B_{0}}$- magnetic field strength

Br- Brinkman number

${\textstyle c_{P}}$- specific heat at constant pressure

${\textstyle E_{G}}$- volumetric rate of entropy

G- dimensionless pressure gradient

h- width of the channel

${\textstyle H^{2}}$- Hartmann number

${\textstyle N_{1}}$- entropy generation due to heat transfer

${\textstyle N_{2}}$- entropy generation due to entropy generation due to fluid friction and ohmic heating

${\textstyle N_{s}}$- entropy generation number

Pr- Prandtl number

s- suction/injection parameter

T- fluid temperature

${\textstyle T_{0}}$- temperature at the lower plate

${\textstyle T_{f}}$- final fluid temperature

u- dimensionless fluid velocity

${\textstyle u^{'}}$- fluid velocity

${\textstyle v_{0}}$- uniform suction/injection velocity

${\textstyle x^{'}{\mbox{,}}y^{'}}$- cartesian coordinates

${\textstyle x{\mbox{,}}y}$- dimensionless cartesian coordinates

1. Introduction

In recent times, there has been a renewed interest in the thermodynamic analysis involving channel fluid flows. This is due to the usefulness of the study in several renewable energy applications. For instance, it could be used to predict the efficiency of many thermal systems exchanging heat between two heat reservoirs and other Carnot systems. Moreover, in the energy generation, excessive energy is wasted or dissipated in the form of heat. Hence there is need to minimise the wastage by improving the energy of the system. Based on this fact, few research works have been reported in the literature. For example, Adesanya and Makinde [1] reported the entropy generation in couple stress fluid flowing steadily through a porous channel with slip at the isothermal walls. Similarly, Adesanya and Makinde [2] studied the entropy generation rate in the couple stress fluid flowing through a porous channel with convective heating at the walls. Also in the class of couple stress fluid, Makinde and Eegunjobi [3] investigated the inherent irreversibility of heat in steady flow of a couple stress fluid through a vertical channel filled with porous materials. Other important work on the entropy generation in a moving fluid includes [4], [5], [6] and [7].

In all the studies above on the thermodynamics analysis, the effect of magnetic field, placed in the transverse direction to the flow channel has been neglected. In reality, magnetic field plays a vital role in many industrial and thermal engineering applications. For instance, it is useful in controlling extremely hot moving fluid like molten steel and many more. Several researchers have investigated the entropy analysis in hydromagnetic fluid flow in recent times. For example, Das and Jana [8] presented the second law analysis for magnetohydrodynamic incompressible fluid flow through a porous channel by imposing Navier slip conditions at the walls. Adesanya and Falade [9] analysed the inherent irreversibility in the flow of hydrodynamic third grade fluid through a channel saturated with porous materials. Similarly, the effect of hall current was presented in a study by Das and Jana [10]. Interested readers can see more interesting work on the influence of magnetic field on entropy generation rate in Refs. [11], [12], [13], [14], [15], [16], [17], [18], [19] and [20].

Motivated by studies in [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19] and [20], the objective of the present study was to examine the influence of magnetic field and ohmic heating of the couple stress fluid on the entropy production within the flow channel, which has not been accounted for in the literature. The outcoming results are expected to enhance many industrial and thermal engineering processes whose working medium is a non-Newtonian fluid, with a view to minimise entropy generation which tends to deplete the amount of available energy for work.

To achieve this objective, flow governing equations are formulated, and non-dimensionalised and approximate solution of the dimensionless coupled non-linear boundary-value problem are obtained by using a semi-analytical Adomian decomposition method [21] and [22]. The choice of the method is due to the fact that the method does not require any linearisation, discretisation, use of initial guess or perturbation. These approximation solutions are used to compute the entropy generation rate and irreversibility ratio.

In the following section, the problem is formulated and non-dimensional analysis is also presented. Section 3 of the work contains the method of solution, results are presented and discussed in Section 4, while Section 5 concludes the paper.

2. Mathematical formulation

A fully developed hydromagnetic non-Newtonian fluid flow between two parallel plates is considered. The parallel plates are infinite, permeable and stationary relative to the fluid motion as shown in Fig. 1. We choose a 2-dimensional cartesian coordinate system with ${\textstyle x}$-axis along the flow direction and ${\textstyle y}$-axis orthogonal to the planes of the parallel plates, separated by width ${\textstyle y=h}$. Fluid injection occurs at the lower plate at a uniform rate ${\textstyle v_{0}}$, matched with a corresponding fluid suction at the upper plate. A constant magnetic field of strength ${\textstyle B_{0}}$ is applied perpendicular to the direction of fluid flow. For most industrial applications, a non-chaotic fluid flow is desired such that the magnetic Reynolds number is very small, and since no external voltage is applied to the fully developed flow system, the induced magnetic field and Hall effect are negligible. We further assumed that all the fluid properties are constant and the Stokes constitutive model for the couple stresses is used. Consequently, the momentum and energy balance equations, with the local volumetric entropy generation rate (${\textstyle E_{G}}$) for the fluid flow can be written as follows [1] and [10]:

$v_{0}\rho {\frac {{du}^{'}}{{dy}^{'}}}=-{\frac {dP}{{dx}^{'}}}+$$\mu {\frac {d^{2}u^{'}}{{dy}^{{'}2}}}-\eta {\frac {d^{4}u^{'}}{{dy}^{{'}4}}}-$${\sigma }_{e}B_{0}^{2}u^{'}{\mbox{,}}$

(1)

$\rho c_{P}v_{0}{\frac {dT}{{dy}^{'}}}=\kappa {\frac {d^{2}T}{{dy}^{{'}2}}}+$$\mu {\left({\frac {{du}^{'}}{{dy}^{'}}}\right)}^{2}+$$\eta {\left({\frac {d^{2}u^{'}}{{dy}^{{'}2}}}\right)}^{2}+$${\sigma }_{e}B_{0}^{2}u^{{'}2}{\mbox{,}}$

(2)

$E_{G}={\frac {\kappa }{T_{0}^{2}}}{\left({\frac {dT}{{dy}^{'}}}\right)}^{2}+$${\frac {\mu }{T_{0}}}{\left({\frac {{du}^{'}}{{dy}^{'}}}\right)}^{2}+$${\frac {\eta }{T_{0}}}{\left({\frac {d^{2}u^{'}}{{dy}^{{'}2}}}\right)}^{2}+$${\frac {{\sigma }_{e}}{T_{0}}}B_{0}^{2}u^{{'}2}{\mbox{,}}$

(3)

with the boundary conditions

$u^{'}(0)={\gamma }_{1}{\frac {{du}^{'}(0)}{{dy}^{'}}}{\mbox{,}}\quad {\frac {d^{2}u^{'}(0)}{{dy}^{{'}2}}}=$$0{\mbox{,}}\quad T(0)=T_{0}{\mbox{,}}$

(4)

$u^{'}(1)={\gamma }_{2}{\frac {{du}^{'}(1)}{{dy}^{'}}}{\mbox{,}}\quad {\frac {d^{2}u^{'}(1)}{{dy}^{{'}2}}}=$$0{\mbox{,}}\quad T(1)=T_{0}{\mbox{,}}$

(5)

where ${\textstyle u^{'}}$ and P   are the fluid velocity and pressure respectively, ${\textstyle v_{0}}$ is the injection/suction velocity at the channel walls, ${\textstyle \eta }$ is the coefficient of couple stress, ${\textstyle \mu }$ is the dynamic viscosity, ${\textstyle {\gamma }_{i}}$ corresponds to the Navier slip coefficients at the lower plate and upper plate, respectively for ${\textstyle i=1{\mbox{,}}2{\mbox{,}}\rho }$ is the fluid density, ${\textstyle {\sigma }_{e}}$ is the fluid electrical conductivity, ${\textstyle \kappa }$ is the coefficient of thermal conductivity, ${\textstyle c_{P}}$ is the isobaric specific heat, ${\textstyle T_{f}}$ and ${\textstyle T_{0}}$ are referenced fluid temperatures, and T is the fluid temperature. The present model has potential application in automated heat control mechanisms, especially in micro-channels. It is often desired in manufacturing processes that the rate of randomness (or entropy generation) be tamed to avoid system complexity.

Figure 1. The geometry of the model.

Introducing the following dimensionless parameters and variables,

$y={\frac {y^{'}}{h}}{\mbox{,}}\quad u={\frac {u^{'}}{v_{0}}}{\mbox{,}}\quad \theta =$${\frac {T-T_{0}}{T_{f}-T_{0}}}{\mbox{,}}\quad s={\frac {v_{0}h}{\nu }}{\mbox{,}}\quad G=$$-{\frac {h^{2}}{\mu v_{0}}}{\frac {dP}{{dx}^{'}}}{\mbox{,}}$ $a^{2}={\frac {\mu h^{2}}{\eta }}{\mbox{,}}\quad H^{2}={\frac {{\sigma }_{e}h^{2}B_{0}^{2}}{\mu }}{\mbox{,}}\quad Pr=$${\frac {\nu \rho c_{P}}{\kappa }}{\mbox{,}}\quad Br={\frac {v_{0}^{2}\mu }{\kappa (T_{f}-T_{0})}}{\mbox{,}}$

(6)

$\Omega ={\frac {T_{f}-T_{0}}{T_{0}}}{\mbox{,}}\quad N_{s}={\frac {T_{0}^{2}h^{2}E_{G}}{\kappa {\left(T_{f}-T_{0}\right)}^{2}}}{\mbox{,}}\quad {\beta }_{1}=$${\frac {{\gamma }_{1}}{h}}{\mbox{,}}\quad {\beta }_{2}={\frac {{\gamma }_{2}}{h}}{\mbox{,}}$

we obtain the following boundary value problems:

$s{\frac {du}{dy}}=G+{\frac {d^{2}u}{{dy}^{2}}}-{\frac {1}{a^{2}}}{\frac {d^{4}u}{{dy}^{4}}}-$$H^{2}u{\mbox{,}}$

(7)

$sPr{\frac {d\theta }{dy}}={\frac {d^{2}\theta }{{dy}^{2}}}+Br{\left({\frac {du}{dy}}\right)}^{2}+$${\frac {Br}{a^{2}}}{\left({\frac {d^{2}u}{{dy}^{2}}}\right)}^{2}+$${BrH}^{2}u^{2}{\mbox{,}}$

(8)

together with the boundary conditions

$u(0)={\beta }_{1}{\frac {du(0)}{dy}}{\mbox{,}}\quad {\frac {d^{2}u(0)}{{dy}^{2}}}=$$0{\mbox{,}}\quad \theta (0)=0$ $u(1)={\beta }_{2}{\frac {du(1)}{dy}}{\mbox{,}}\quad {\frac {d^{2}u(1)}{{dy}^{2}}}=$$0{\mbox{,}}\quad \theta (1)=0{\mbox{,}}$

(9)

while the entropy generation can be computed using

$N_{s}={\left({\frac {d\theta }{dy}}\right)}^{2}+{\frac {Br}{\Omega }}{\left({\frac {du}{dy}}\right)}^{2}+$${\frac {Br}{\Omega a^{2}}}{\left({\frac {d^{2}u}{{dy}^{2}}}\right)}^{2}+$${\frac {Br}{\Omega }}H^{2}u^{2}{\mbox{.}}$

(10)

Our model given by Eqs. (7), (8) and (10) reduces to the asymptotic case as ${\textstyle H^{2}\rightarrow 0}$ as obtained in Ref. [1].

3. Adomian decomposition method of solution

In this section, we seek the solution of the model Eqs. (7), (8) and (10) subject to the boundary condition (9), using the semi-analytical Adomian decomposition scheme. To achieve this, the boundary value problems are converted to the equivalent integral equations noting that ${\textstyle u^{''}(0)=u^{''}(1)=0}$. Thus,

$u(y)=f_{0}+f_{1}y+{\frac {y^{3}}{6}}f_{2}+a^{2}{\int }_{0}^{y}{\int }_{0}^{y}{\int }_{0}^{y}{\int }_{0}^{y}F_{4}(Y)\quad dY\quad dY\quad dY\quad dY{\mbox{,}}$

(11)

$\theta (y)=c_{0}+c_{1}y+{\int }_{0}^{y}{\int }_{0}^{y}\left(F_{2}(Y)-\right.$$\left.{BrH}^{2}u^{2}\right)\quad dY\quad dY$

(12)

where

$F_{4}(Y)=G+{\frac {d^{2}u}{{dY}^{2}}}-s{\frac {du}{dY}}-H^{2}u{\mbox{,}}$ $F_{Y}(Y)=sPr{\frac {d\theta }{dY}}-Br{\left({\frac {du}{dY}}\right)}^{2}-$${\frac {Br}{a^{2}}}{\left({\frac {d^{2}u}{{dY}^{2}}}\right)}^{2}{\mbox{.}}$

The unknown coefficients ${\textstyle f_{0}{\mbox{,}}f_{1}{\mbox{,}}f_{2}{\mbox{,}}c_{0}}$ and ${\textstyle c_{1}}$ are to be determined by using the method of undetermined coefficients. Following ADM we seek solutions in the form of an infinite series:

$u(y)=\sum _{0}^{\infty }u_{n}(y){\mbox{,}}\quad \theta (y)=$$\sum _{0}^{\infty }{\theta }_{n}(y){\mbox{.}}$

(13)

Substituting (13) into (11) and (12), we obtain

$\sum _{n=0}^{n=\infty }u_{n}(y)=f_{0}+f_{1}y+{\frac {y^{3}}{6}}f_{2}+$$a^{2}L_{4}(Y)\quad dY\quad dY\quad dY\quad dY{\mbox{,}}$

(14)

and

$\sum _{n=0}^{n=\infty }{\theta }_{n}(y)=c_{0}+c_{1}y+{\int }_{0}^{y}{\int }_{0}^{y}\left(Z_{1}-\right.$$\left.Z_{2}\right)\quad dY\quad dY{\mbox{;}}$

(15)

where

$L_{4}(Y)={\int }_{0}^{y}{\int }_{0}^{y}{\int }_{0}^{y}{\int }_{0}^{y}\left(G+\right.$$\left.\sum _{n=0}^{n=\infty }{\frac {d^{2}u_{n}}{{dY}^{2}}}-\right.$$\left.s\sum _{n=0}^{n=\infty }{\frac {{du}_{n}}{dY}}-H^{2}\sum _{n=0}^{n=\infty }u_{n}\right){\mbox{,}}$

$Z_{1}=sPr\sum _{n=0}^{n=\infty }{\frac {d{\theta }_{n}}{dY}}-$$Br\sum _{n=0}^{n=\infty }{\left({\frac {{du}_{n}}{dY}}\right)}^{2}{\mbox{,}}$

and

$Z_{2}={\frac {Br}{a^{2}}}\sum _{n=0}^{n=\infty }{\left({\frac {d^{2}u_{n}}{{dY}^{2}}}\right)}^{2}-$${BrH}^{2}\sum _{n=0}^{n=\infty }u_{n}^{2}{\mbox{.}}$

Evidently (13) may be written in the recursive form:

$u_{0}(y)=f_{0}+f_{1}y+{\frac {y^{3}}{6}}f_{2}+a^{2}{\int }_{0}^{y}{\int }_{0}^{y}{\int }_{0}^{y}{\int }_{0}^{y}G\quad dY\quad dY\quad dY\quad dY{\mbox{,}}$

(16)

$u_{n+1}(y)=a^{2}{\int }_{0}^{y}{\int }_{0}^{y}\left({\frac {d^{2}u_{n}}{{dY}^{2}}}-\right.$$\left.s{\frac {{du}_{n}}{dY}}-H^{2}u_{n}\right)\quad dY\quad dY{\mbox{.}}$

(17)

In order to reduce the computational load in (15), the modified recurrent relation

${\theta }_{0}(y)=c_{1}y{\mbox{,}}$

(18)

is used, so that

${\theta }_{1}(y)={\int }_{0}^{y}{\int }_{0}^{y}\left(sPr\sum _{n=0}^{n=\infty }{\frac {d{\theta }_{0}}{dY}}-\right.$$\left.{\frac {Br}{a^{2}}}{\left({\frac {d^{2}u}{{dY}^{2}}}\right)}^{2}\right)\quad dY\quad dY{\mbox{.}}$

(19)

${\theta }_{2}(y)={\int }_{0}^{y}{\int }_{0}^{y}\left(sPr\sum _{n=0}^{n=\infty }{\frac {d{\theta }_{1}}{dY}}-\right.$$\left.Br{\left({\frac {du}{dY}}\right)}^{2}\right)\quad dY\quad dY{\mbox{.}}$

(20)

${\theta }_{3}(y)={\int }_{0}^{y}{\int }_{0}^{y}sPr\sum _{n=0}^{n=\infty }\left({\frac {d{\theta }_{2}}{dY}}-\right.$$\left.{BrH}^{2}u^{2}\right)\quad dY\quad dY{\mbox{.}}$

(21)

${\theta }_{n+1}(y)={\int }_{0}^{y}{\int }_{0}^{y}\left(sPr\sum _{n=0}^{n=\infty }{\frac {d{\theta }_{n}}{dY}}\right)\quad dY\quad dY{\mbox{;}}\quad (n\geqslant 2){\mbox{,}}$

(22)

Due to convergence of the ADM series solution conducted in Ref. [6] only few terms of the series will be required to obtain the approximate solutions of the problem. We set the number of iteration to m so that the approximate solutions (13) may be written as the finite series as follows:

$u(y)=\sum _{0}^{m}u_{n}(y){\mbox{,}}\quad \theta (y)=\sum _{0}^{m}{\theta }_{n}(y){\mbox{.}}$

(23)

The irreversibility in the heat flow to the viscous fluid is analysed using (10) by expressing the entropy generation number ${\textstyle N_{s}}$ as partial sum of the entropy generation due to heat transfer, and the irreversibility resulting from fluid friction and ohmic heating of the fluid. Therefore we set ${\textstyle N_{1}}$ and ${\textstyle N_{2}}$ as follows:

$N_{1}={\left({\frac {d\theta }{dy}}\right)}^{2}{\mbox{,}}\quad N_{2}=$${\frac {Br}{\Omega }}{\left({\frac {du}{dy}}\right)}^{2}+{\frac {Br}{\Omega a^{2}}}{\left({\frac {d^{2}u}{{dy}^{2}}}\right)}^{2}+$${\frac {Br}{\Omega }}H^{2}u^{2}{\mbox{.}}$

(24)

The Bejan number Be, which is a parameter that measures the irreversibility ratio in the heat flow, may be defined as

$Be={\frac {N_{1}}{N_{s}}}={\frac {1}{1+\Phi }}{\mbox{,}}\quad \Phi =$${\frac {N_{2}}{N_{1}}}{\mbox{.}}$

(25)

where ${\textstyle \Phi }$ is the irreversibility ratio, a parameter that measures the rate of destruction of available work in the flow system. The Bejan number in Eq. (25) is bounded in the interval ${\textstyle 0\leqslant Be\leqslant 1}$. The irreversibility due to heat transfer is dominant when ${\textstyle Be=1}$, while the irreversibility due to viscosity and magnetic field is dominant when ${\textstyle Be=0}$. The dimensionless Eqs. (7), (8), (9) and (10), with the boundary conditions (9) were solved using the algorithm in (16), (17), (18), (19), (20), (21), (22) and (23), coded in MATHEMATICA symbolic package. Using the numerical procedure discussed above, we computed the dimensionless velocity, temperature, entropy generation and irreversibility ratio. In what follows, we discussed some interesting results from our findings.

4. Results and discussion

Let us begin by examining the effect of the applied magnetic field and couple stresses on the flow velocity, the variations of the dimensionless velocity with the width of the flow channel as presented in Fig. 2, for varying values of magnetic parameter (Hatman number, ${\textstyle H^{2}}$), couple stress inverse parameter ${\textstyle a^{2}}$, suction/injection parameter s  , upper wall Navier-slip parameter ${\textstyle {\beta }_{2}}$, and lower wall Navier-slip parameter ${\textstyle {\beta }_{1}}$. We showed in Table 1 that the ADM could guarantee a result comparable with the exact solution of the problem, with absolute error of order ${\textstyle {10}^{-7}}$ for ${\textstyle m=4}$. It is clear from Fig. 2a that by increasing the magnetic field the flow velocity would be inhibited. This could be explained by the action of Lorentz force, whose presence in the flow system constitutes resistance to the momentum of fluid parcels in the adjacent fluid layers. Fig. 2b depicts the implication of increasing couple stress inverse parameter on the flow velocity. Here, we find that an increase in couple stress inverse parameter corresponds to increase in fluid friction in the moving fluid. This friction arises from the effect of particle additives, constituting size dependent effect on couple stress fluids. Moreover, rotational field of the velocity is generated in couple stress fluid. Hence, an increase in couple stress results in decrease in the velocity profile of the fluid. Fig. 2c shows the effect of suction/injection on the velocity profile of the fluid flow. Clearly, the increase in suction/injection parameter breaks the symmetry of the fluid flow. This is because the continual injection of fluid into the channel, from a direction different from that of channel fluid flow, changes the momentum of (or slows down) the fluid parcels at the lower wall-fluid interface. This effect is propagated beyond the neighbourhood of the lower plate, which breaks the symmetry of the channel flow system. However, increasing the suction/injection parameter does not show any significant effect at the upper wall of the channel. Fig. 2d and e depicts the effects of Navier-slip parameters on the velocity profile of the fluid flow. Evidently, the fluid slippage has significant effect on the fluid motion, in that any change in any of the parameters impacts on the viscous drag on the flow, which scales the maximum velocity of the flow up or down. In clear terms, increasing the upper wall Navier-slip parameter decreases the fluid velocity significantly at the upper wall as well as the maximum velocity of the fluid flow. In contrast, increasing the slip parameter at the lower wall increases the velocity at the wall as well as the maximum velocity of the fluid flow.

Figure 2. Parameterised velocity profiles of the fluid flow: (a) at varying Hatman number H  , with ${\textstyle a=s=1{\mbox{,}}\quad {\beta }_{1}={\beta }_{2}=0.1}$; (b) at varying couple stress inverse a  , with ${\textstyle H=s=1{\mbox{,}}\quad {\beta }_{1}={\beta }_{2}=0.1}$; (c) at varying suction/injection parameter s  , with ${\textstyle H=1{\mbox{,}}\quad a=5{\mbox{,}}\quad {\beta }_{1}={\beta }_{2}=}$$0.1$; (d) at varying lower wall Navier-slip parameter ${\textstyle {\beta }_{1}}$, with ${\textstyle H=a=s=1{\mbox{,}}\quad {\beta }_{2}=0.1}$; (e) at varying upper wall Navier-slip parameter ${\textstyle {\beta }_{2}}$, with ${\textstyle H=a=s=1{\mbox{,}}\quad {\beta }_{1}=0.1}$.

Table 1. Comparison between analytical and exact solution for the velocity profile for the parameters ${\textstyle a=1{\mbox{,}}\quad H=0.2{\mbox{,}}\quad s={\beta }_{1}={\beta }_{2}=}$$0.1$.

u	${\textstyle u_{Exact}}$	${\textstyle u_{ADM}}$	Absolute error
0	0.00303054	0.00303070	${\textstyle 1.56912\times {10}^{-7}}$
0.1	0.00598827	0.00598858	${\textstyle 3.11550\times {10}^{-7}}$
0.2	0.00854151	0.00854196	${\textstyle 4.52517\times {10}^{-7}}$
0.3	0.01037320	0.01037370	${\textstyle 5.66015\times {10}^{-7}}$
0.4	0.01125880	0.01125950	${\textstyle 6.38135\times {10}^{-7}}$
0.5	0.01106440	0.01106510	${\textstyle 6.55205\times {10}^{-7}}$
0.6	0.00974517	0.00974577	${\textstyle 6.04755\times {10}^{-7}}$
0.7	0.00734506	0.00734554	${\textstyle 4.77423\times {10}^{-7}}$
0.8	0.00399756	0.00399783	${\textstyle 2.70091\times {10}^{-7}}$
0.9	−0.00007304	−0.00007305	${\textstyle 9.47647\times {10}^{-9}}$
1.0	−0.00454921	−0.00454955	${\textstyle 3.36326\times {10}^{-7}}$

Next we describe the temperature response of the fluid with respect to flow to changes in the values of the magnetic parameter, couple stress parameter and the Brinkman number Fig. 3. Fig. 3a shows that increasing the magnetic field intensity (or parameter) translates into a decrease in the channel temperature. Fig. 3b reveals that increasing couple stresses decrease the temperature profile in the channel. Expectedly, couple stresses enhance the fluid’s intermolecular cohesion. This increases the fluid’s resistance to shear stress, resulting in temperature rise within the flow channel. Fig. 3c shows that an increase in Brinkman number results in the increase in the fluid temperature.

Figure 3. Parameterised temperature profiles of the fluid flow: (a) at varying values of Hatman number H  , with ${\textstyle a=s=Br=1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$; (b) at varying values of couple stress inverse a  , with ${\textstyle H=s=Br=1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$; (c) at varying values of the Brinkman number Br  , with ${\textstyle H=a=s=1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$.

To complete the analysis, the effects of the magnetic parameter, ${\textstyle H^{2}}$, couple stress parameter, ${\textstyle a^{2}}$, and the Brinkman number, Br on the entropy generation are shown in Fig. 4. Interestingly, it is observed from Fig. 4a that the entropy generation, ${\textstyle N_{s}}$ decreases with increase in magnetic parameter ${\textstyle H^{2}}$ in the middle of the channel. The interplay between the two thermophysical phenomena, namely the intermolecular cohesion and entropy generation rate in viscous fluid becomes imperative. Viscosity increases as the intermolecular cohesion within the fluid rises, or vice versa. The applied Lorentz force reduces the intermolecular cohesion within the fluid resulting in decrease in entropy generation ${\textstyle N_{s}}$ within the channel flow. Fig. 4b and c shows that the rate of entropy generation increases with increase in both couple stress parameter and the Brinkman number. Increase in Brinkman number raises both dissipations due to Joule heating and viscosity, which in turn increases the entropy generation rate. Fig. 5 shows the dominance of fluid viscosity on the irreversibility ratio, over heat transfer at the centreline of the channel. In Fig. 5a, it is observed that as the magnetic parameter ${\textstyle H^{2}}$ increases, heat transfer dominates the irreversibility ratio at the walls, while Fig. 5a and b revealed that as the couple stress inverse and the Brinkman number increase, fluid viscosity dominates the irreversibility ratio.

Figure 4. Entropy generation rate: (a) at varying values of Hatman number H  , with ${\textstyle a=s=Br=\Omega =1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$; (b) at varying values of couple stress inverse a  , with ${\textstyle H=s=Br=\Omega =1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$; (c) at varying values of the Brinkman number Br  , with ${\textstyle H=a=s=\Omega =1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$.

Figure 5. Irreversibility ratio: (a) at varying values of Hatman number H  , with ${\textstyle a=s=Br=\Omega =1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$; (b) at varying values of couple stress inverse a  , with ${\textstyle H=s=Br=\Omega =1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$; (c) at varying values of the Brinkman number Br  , with ${\textstyle H=a=s=\Omega =1{\mbox{,}}\quad Pr=0.71{\mbox{,}}\quad {\beta }_{1}=}$${\beta }_{2}=0.1$.

Following from the discussion above, at specific ${\textstyle a^{2}{\mbox{,}}s{\mbox{,}}{\beta }_{1{\mbox{,}}2}}$ and Br, varying the magnetic field automatically changes the fluid’s rate of change of momentum, which in turn changes the fluid velocity. Reducing velocity in micro-channel fluid flow effectively reduces fluid temperature, rate of entropy production in the fluid, and the irreversibility ratio as shown in Figure 2, Figure 3, Figure 4 and Figure 5.

5. Conclusion

The effects of magnetic field and couple stresses on entropy generation in an MHD flow through a porous channel, with suction/injection and Navier slip have been examined. The Adomian decomposition method was employed to obtain semi-analytical solution that approximates the velocity and temperature profiles, which are used to obtain the entropy generation production as well as the Bejan number. It was found that the entropy generation rate reduces with increasing magnetic parameter. The magnetic field acts in such a way as to tame the degree of randomness in the fluid’s particles, resulting in the lowering of the fluid’s velocity and entropy generation rate in the system. In effect, the fluid velocity, the fluid temperature and entropy generation showed significant dependence on Joule dissipation. This result suggests a balanced mix between heat transfer irreversibility and fluid friction irreversibility, which would enhance entropy minimisation and hence system design and manufacturing applications. Our result would motivate applications in micro-channel automatic heat control systems.

Acknowledgement

UEV is supported by the Royal Society of London, through their Newton International Fellowship Alumni scheme.

References

1. [1] S.O. Adesanya, O.D. Makinde; Entropy generation in couple stress fluid flow through porous channel with fluid slippage; Exergy, 15 (3) (2014), pp. 344–362
2. [2] S.O. Adesanya, O.D. Makinde; Effects of couple stresses on entropy generation rate in a porous channel with convective heating; Comput. Appl. Math., 34 (2015), pp. 293–307
3. [3] O.D. Makinde, A.S. Eegunjobi; Entropy generation in a couple stress fluid flow through a vertical channel filled with saturated porous media; Entropy, 15 (11) (2013), pp. 4589–4606
4. [4] J.A. Esfahani, M. Modirkhazeni; Entropy generation of forced convection film condensation on a horizontal elliptical tube; C. R. Mcanique, 340 (7) (2012), pp. 543–551
5. [5] Y. Wang, Z. Chen, X. Ling; Entropy generation analysis of particle suspension induced by Couette flow; Int. J. Heat Mass Transfer, 90 (2015), pp. 499–504
6. [6] S.O. Adesanya, O.D. Makinde; Thermodynamic analysis for a third grade fluid through a vertical channel with internal heat generation; J. Hydrodyn., 27 (2) (2015), pp. 264–272
7. [7] S.O. Adesanya, O.D. Makinde; Irreversibility analysis in a couple stress film flow along an inclined heated plate with adiabatic free surface; Physica A, 432 (2015), pp. 222–229
8. [8] S. Das, R.N. Jana; Entopy generation due to MHD flow in a porous channel with Navier slip; Ain Shams Eng. J., 5 (2014), pp. 575–584
9. [9] S.O. Adesanya, J.A. Falade; Thermodynamics analysis of hydromagnetic third grade fluid flow through a channel filled with porous medium; Alexandria Eng. J., 54 (3) (2015), pp. 615–622
10. [10] S. Das, R.N. Jana; Effect of hall current on entropy generation in porous channel with suction/injection; Int. J. Energy Technol., 5 (25) (2013), pp. 1–11
11. [11] C.K. Chen, B.S. Chen, C.C. Liu; Entropy generation in mixed convection magnetohydrodynamic nanofluid flow in vertical channel; Int. J. Heat Mass Transfer, 91 (2015), pp. 1026–1033
12. [12] N. Saeid; Magnetic field effects on entropy generation in heat and mass transfer in porous cavity; Acad. J., 8 (17) (2013), pp. 728–739
13. [13] M. Magherbi, A. El-Jery, N. Hidouri, A.B. Brahim; Evanescent magnetic field effects on entropy generation at the on set of natural convection; Indian Acad. J., 35 (2) (2010), pp. 163–176
14. [14] A. El-Jery, N. Hidouri, M. Magherbi, A.B. Brahim; Effect of external oriented magnetic field on entropy generation in natural convection; Entropy, 12 (2010), pp. 1391–1417
15. [15] S. Salas, S. Cuevas, M.L. Haro; Entropy analysis of magnetohydrodynamic induction devices; J. Phys. D: Appl. Phys., 32 (1999), pp. 2605–2608
16. [16] S.H. Tasnim, M.A.H. Mahmud; Entropy generation in a porous channel with hydromagnetic effects; Exergy, 2 (2002), pp. 300–308
17. [17] S. Mahmud, S.H. Tasnim, H.A.A. Mamun; Thermodynamics analysis of mixed convection in a channel with transverse hydromagnetic effect; Int. J. Therm. Sci., 42 (731–740) (2003)
18. [18] D. Theuri, O.D. Makinde; Thermodynamic analysis of variable viscosity MHD unsteady generalized Coutte flow with permeable walls; Appl. Comput. Math., 3 (1) (2014), pp. 1–8
19. [19] O. Mahian, H.F. Oztop, I. Pop, S. Mahmud, S. Wongwises; Design of vertical annulus with MHD flow using entropy generation analysis; Therm. Sci., 17 (4) (2013), pp. 1013–1022
20. [20] S.O. Adesanya, S.O. Kareem, J.A. Falade, S.A. Arekete; Entropy generation analysis for a reactive couple stress fluid flow through a channel saturated with porous material; Energy, 93 (1) (2015), pp. 1239–1245
21. [21] G. Adomian; Solving Frontier Problems of Physics: The Decomposition Method; Kluwer Academic Publishers, Dordrecht (1994)
22. [22] J.S. Duan, R. Rach; A new modification of the Adomian decomposition method for solving boundary value problems for higher order nonlinear differential equations; Appl. Math. Comput., 218 (2011), pp. 4090–4118