Proceedings of the World Congress on Engineering 2017 Vol I

WCE 2017, July 5-7, 2017, London, U.K.

Irreversibility Analysis of a Radiative MHD

Poiseuille Flow through Porous Medium with
Slip Condition
A. A. Opanuga* Member, IAENG, H.I. Okagbue, O.O. Agboola


Abstract—In this article, irreversibility analysis of thermal II. MATHEMATICAL FORMULATION
radiation with slip condition on MHD Poiseuille flow through
The assumptions made include:
porous medium is investigated. The upper and lower walls are
kept constant with the same temperature. The radiative heat The flow is steady, electrically conducting and
flux in the energy equation is assumed to follow Roseland incompressible; the fluid is viscous and flow through parallel
approximation. Semi-analytical solutions of the non-linear porous medium; both plates are fixed and maintained at
boundary value problems obtained from the governing uniform temperature; uniform transverse magnetic field B0
equations is constructed using Adomian decomposition method,
and the effects of some fluid parameters on fluid motion, is applied neglecting the induced magnetic field and the Hall
temperature, entropy generation and Bejan number are effect; Navier slip boundary condition is assumed at the
presented. fluid-solid interface; the fluid is optically thick following
Roseland approximation.
Index Terms— Irreversibility, radiation, MHD, Poiseuille The governing equations are given as [8, 9]
flow, slip condition
d 2u  B02u bu dp
    (1)
d 2  K d
I. INTRODUCTION
2
ecent researches reveal that more attention has been d 2T   du  B02u2 bu2
 
R devoted to the preservation of scarce resources. This
has led to the investigation of the causes of
k
d 2
  
 d 



K

(2)
irreversibility in various flow systems; some of these are dqr
found in Refs. [1-5]. In addition, Arikoglu [6] submitted 0
that, all energy producing, converting and consuming d
systems must be re-examined carefully and possible 2
k  dT     du   B02u2
available-work destruction mechanisms be removed. EG  2      
Available research works show that the effect of T0  d   T0  d   T0
velocity slip on entropy generation of plane Poiseuille flow (3)
has not been fully addressed. Few investigations on this  u 2
subject are [7-9]. Motivated by [8, 9], this article examines T0 K
the entropy generation due to thermal radiation and velocity
du(0) du(0)
slip on MHD Poiseuille flow through porous medium. u (0)   1 , u ( h)   2 ;
Numerous semi-analytical methods for solving d  d 
boundary value problems are found in literature, most of
these techniques have difficulties in relation to the size of
T (0)  T0 , T (h)  Th (4)
computational work and convergence. However the The Roseland approximation term for optimally thick fluid
technique of Adomian Decomposition Method (ADM) [10-, is written as
12] applied in this article is easy to apply with high accuracy 4 c dT 4
and rapid convergence. qr  (5)
3k c d 4
The temperature term (T  ) in equation (5) can be
4

expressed in term of its linearity function as given by Raptis

et al. [13], then the expansion in Taylor series about T0
Manuscript received February, 13, 2017; revised March 10, 2017. This
work was supported by the Centre for Research and Innovation, Covenant gives
T 4  T 4  4T03 T   T0  6T02 T   T0  
University, Ota, Nigeria. 2
A. A. Opanuga, H.I. Okagbue, and O.O. Agboola are with the Department

4T0 T   T0  T   T0

of Mathematics,Covenant University, Ota, Nigeria.(e-mail: 3 4 (6)
[email protected],
[email protected],
[email protected]).

ISBN: 978-988-14047-4-9 WCE 2017

ISSN: 2078-0958 (Print); ISSN: 2078-0966 (Online)
Proceedings of the World Congress on Engineering 2017 Vol I
WCE 2017, July 5-7, 2017, London, U.K.

Using equations (5) and (6) in equation (2) and neglecting The Bejan number assumes values between 0 and 1.
higher order terms, we obtain Be  0 for VDI  , Be  1 for  HTI  and Be  0.5 is
2 when both VDI and HTI contribute equally to entropy
d 2T   du   B02u2 u2
k        generation. Then setting
d 2  d   K N1 1 N
(7) Be   ,  2 (14)
16 cT03d 2T  Ns 1   N1
3k c d 2 where
 4 d 
2
The dimensionless expressions for the present problem are: N1  1  R  2 ,
 u T   T0 h 2 dp  3  d
 ,u  ,  ,A ,
h U T f  T0  dx Br   du 
2

N2     H 2u 2   2u 2  (15)
v c p U 2 2
T h EG    d  
Pr  , Br  , Ns  0
, 
k k (T f  T0 ) k (T f  T0 ) 2
(8)
kT f  T0  B02 h 2 4 cT03
 ,H2  ,R  ,
T0  kk c IV. RESULTS AND DISCUSION
h 2
  In this article, the effect of Navier slip and thermal
 , 1  1 ,  2  2 radiation are investigated on the entropy generation of MHD
K h h Poiseuille flow through porous medium. The effects of some
Applying the above dimensionless variables in equations
parameters on fluid velocity, temperature, entropy
(1, 3, 4, 7) yields
generation and Bejan number are presented in this section.
d 2u Figs. 1 and 2 depict the effect of slip parameters on
 H 2u   2u  A  0 (9)
d 2 fluid velocity. It is observed from Fig. 1 that fluid velocity
increases with increase in lower wall slip parameter while
2
 4  d   du 
2
the situation is reversed with upper wall slip parameter in
 1  R  2    BrH u 
2 2
Fig. 2. In Fig. 3, we present the effect of radiation parameter
 3  d  d  (10)
on the temperature. It is obvious that fluid temperature is
Br u  0
2 2 lowered with increased values of radiation parameter. This is
caused by the absorption of heat emitted by the absorptivity
 du 2  parameter. Figs. 4 and 5 reveal that fluid temperature is
 4  d  Br 
2
  
Ns  1  R  2   d   (12) enhanced by increase in slip parameter.
 3  d 
2 2
Furthermore, Figs. 6 and 7 depict that entropy
H u   u 
2 2
generation is retarded at the lower wall while it is enhanced
at the upper wall. Also, Fig. 8 is the plot of thermal radiation
du (0) du (0)
u (0)  1 , u (1)   2 ; effect on entropy generation. The Figure shows that entropy
d d (13) generation is significantly increased with increase in
radiation parameter (R).
 (0)  0, (1)  1 Finally, Figs. 9 and 10 show similar results. In the plots
Solving equations (9-10) by ADM yields the solution of the the Bejan number increases at the lower wall while there a
boundary value problems. reduction in the middle and upper walls of the channel. This
is an indication that heat transfer irreversibility dominates
III. ENTROPY GENERATION entropy generation at the lower wall while viscous
The dimensionless entropy generation expression in dissipation irreversibility is the major contributor to
equation (11) provides four sources of irreversibility, that is irreversibility at the upper wall. In Fig. 11 a rise in thermal
equation (11) is of the form; radiation parameter leads to an increase in Bejan number
across the channel. This shows that heat transfer
 4 d 
2
HTI  TRI  1  R  2 heat transfer and thermal irreversibility is the dominant contributor to entropy
 3  d generation.
radiation irreversibility;
2
Br  du 
VDI    viscous dissipation irreversibility;
ACKNOWLEDGMENT
  d  The authors sincerely appreciate Covenant University for the
2 2 financial support, and the constructive comments of the
BrH u
MFI  magnetic field irreversibility and reviewers.

Br 2u 2
PI  porosity irreversibility.


ISBN: 978-988-14047-4-9 WCE 2017

ISSN: 2078-0958 (Print); ISSN: 2078-0966 (Online)
Proceedings of the World Congress on Engineering 2017 Vol I
WCE 2017, July 5-7, 2017, London, U.K.

Fig 1: Velocity versus lower wall slip parameter

Fig 5: Temperature versus upper wall slip parameter

Fig 2: Velocity versus upper wall slip parameter

Fig 6: Entropy generation versus lower wall slip

parameter

Fig 3: Temperature versus thermal radiation parameter

Fig 7: Entropy generation versus upper wall slip parameter

Fig 4: Temperature versus lower wall slip parameter

Fig 8: Entropy generation versus Thermal radiation
parameter

ISBN: 978-988-14047-4-9 WCE 2017

ISSN: 2078-0958 (Print); ISSN: 2078-0966 (Online)
Proceedings of the World Congress on Engineering 2017 Vol I
WCE 2017, July 5-7, 2017, London, U.K.

[4] S. O. Adesanya, S. O. Kareem, J. A. Falade and S. A. Arekete,

Entropy generation analysis for a reactive couple stress fluid
flow through a channel saturated with porous material,
Energy, 93(2015): 1239-1245
[5] S. O. Adesanya, Second law analysis for third-grade fluid
with variable properties, Journal of Thermodynamics Volume
2014 (2014), 8pages http://dx.doi.org/10.1155/2014/452168
[6] A. Arikoglu, I. Ozkol and G. Komurgoz, effect of slip on
entropy generation in a single rotating disk in MHD flow,
Applied Energy 85 (2008): 1225–1236
[7] T. Hayat, S. Hina and N. Ali, Simultaneous effects of slip and
heat transfer on the peristaltic flow, Commun Nonlinear Sci
Numer Simulat 15 (2010):1526–1537
[8] S.O. Adesanya and O. D. Makinde, Entropy generation in
couple stress fluid flow through porous channel with fluid
slippage, International Journal of Exergy, 15(3)(2014):344 –
362.
Fig 9: Bejan number versus lower wall parameter [9] A.S. Eegunjobi and O.D. Makinde , Effects of Navier slip on
entropy generation in a porous channel with suction/injection,
Journal of Thermal Science and Technology, 7(4)(2012): 522-
535
[10] A.A. Opanuga, O.O. Agboola and H.I. Okagbue,
“Approximate solution of multipoint boundary value
problems”, Journal of Engineering and Applied Sciences, vol.
10, no 4, pp. 85-89, 2015.
[11] A.A. Opanuga, O.O. Agboola , H.I. Okagbue and J.G.
Oghonyon (2015a), “Solution of differential equations by
three semi-analytical techniques”, International Journal of
Applied Engineering Research, vol. 10, no 18, pp. 39168-
39174, 2015.
[12] A.A. Opanuga, H.I. Okagbue, E.A. Owoloko, and O.O.
Agboola,“Modified Adomian decomposition method for
thirteenth order boundary value problems”, Gazi University
Journal of Science, (in press)
[13] A. Raptis, Perdikis, H.S. Takhar, “Effect of thermal radiation
Fig 10: Bejan number versus upper wall parameter on MHD flow”, Applied Mathematics and Computation, vol.
153: pp. 645–649

u axial velocity
 dynamic viscosity
p fluid pressure
h channel width
 fluid density
T fluid temperature
T0 initial fluid temperature
Tf final fluid temperature
k thermal conductivity of the fluid
Fig 11: Bejan number versus thermal radiation parameter cp specific heat at constant pressure
 electrical conductivity of the fluid Navier slip
coefficients
B0 uniform transverse magnetic field
REFERENCES
[1] S.O. Adesanya and O. D. Makinde, Effects of couple stress on qr radiative heat flux
entropy generation rate in a porous channel with convective u dimensionless velocity
heating, Comp. Appl. Math.,34(2015):293–307  dimensionless temperature
[2] A.O. Ajibade, B.K. Jha and A. Omame, Entropy generation
under the effect of suction/injection, Applied Mathematical Pr Prandtl number
Modelling, 35(2011):4630-4646 Br Brinkman number
[3] A.A. Opanuga, J.A. Gbadeyan, S.A. Iyase and H.I. Okagbue,  parameter that measures the temperature
Effect of Thermal Radiation on the Entropy Generation of difference between the two heat reservoirs
Hydromagnetic Flow Through Porous Channel, The Pacific
Journal of Science and Technology, 17( 2)(2016):59-68 b empirical constant in the second order (porous
inertia resistance)

ISBN: 978-988-14047-4-9 WCE 2017

ISSN: 2078-0958 (Print); ISSN: 2078-0966 (Online)
Proceedings of the World Congress on Engineering 2017 Vol I
WCE 2017, July 5-7, 2017, London, U.K.

H magnetic field parameter

Be Bejan number
A axial pressure gradient
R thermal radiation parameter
K porous media permeability
1,2 Navier slip parameters respectively
 porous media shape parameter
EG local volumetric entropy generation rate
Ns dimensionless entropy generation rate
 is the kinematic viscosity
c Stefan-Boltzman constant
kc mean absorption coefficient for thermal radiation.

ISBN: 978-988-14047-4-9 WCE 2017

ISSN: 2078-0958 (Print); ISSN: 2078-0966 (Online)