Next Article in Journal
Applications of Information Theory in Solar and Space Physics
Next Article in Special Issue
Effects of Radiative Electro-Magnetohydrodynamics Diminishing Internal Energy of Pressure-Driven Flow of Titanium Dioxide-Water Nanofluid due to Entropy Generation
Previous Article in Journal
A Neighborhood Rough Sets-Based Attribute Reduction Method Using Lebesgue and Entropy Measures
Previous Article in Special Issue
Entropy Generation of Carbon Nanotubes Flow in a Rotating Channel with Hall and Ion-Slip Effect Using Effective Thermal Conductivity Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Generation in MHD Mixed Convection Non-Newtonian Second-Grade Nanoliquid Thin Film Flow through a Porous Medium with Chemical Reaction and Stratification

1
Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Khyber Pakhtunkhwa, Pakistan
2
Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 736464, Vietnam
3
Department of Mechatronics and System Engineering, College of Engineering, Majmaah University, Al-Majmaah 11952, Saudi Arabia
4
Department of Mechanical and Industrial Engineering, College of Engineering, Majmaah University, Al-Majmaah 11952, Saudi Arabia
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(2), 139; https://doi.org/10.3390/e21020139
Submission received: 8 December 2018 / Revised: 23 January 2019 / Accepted: 30 January 2019 / Published: 1 February 2019
(This article belongs to the Special Issue Entropy Generation in Nanofluid Flows II)

Abstract

:
Chemical reaction in mixed convection magnetohydrodynamic second grade nanoliquid thin film flow through a porous medium containing nanoparticles and gyrotactic microorganisms is considered with entropy generation. The stratification phenomena, heat and mass transfer simultaneously take place within system. Microorganisms are utilized to stabilize the suspended nanoparticles through bioconvection. For the chemical reaction of species, the mass transfer increases. The governing equations of the problem are transformed to nonlinear differential equations through similarity variables, which are solved through a well known scheme called homotopy analysis method. The solution is expressed through graphs and illustrations which show the influences of all the parameters. The residual error graphs elucidate the authentication of the present work.

1. Introduction

Entropy is a thermodynamic quantity representing the unavailability of a system’s thermal energy for conversion into mechanical work, often interpreted as the degree of disorder or randomness in the system. Entropy generation is an important part of the mechanical system. In the context of the present situation, irreversibility (entropy generation) is reliant to heat transmission, viscous dissipation and effects of magnetic field as well as diffusion. Irreversibility is studied with a sharp rate. The present study also discusses bioconvection and chemical reaction. Bio-convection is associated with a macroscopic flow of convecting fluid under the effect of density gradients due to hydrodynamic propulsion such as swimming of motile microorganisms. Mixed convection has several applications including heat dismissal into the atmosphere like lakes, rivers and seas. It is also utilized in the storage of thermal energy like solar ponds and heat transfer thermal sources such as the condensers of power plants. The applications of biosynthesized nanoparticles in a wide spectrum of potential areas are targeted drug delivery, cancer treatment gene therapy and DNA analysis, antibacterial agents, biosensors, enhancing reaction rates, separation science and magnetic resonance imaging (MRI). Entropy generation and bioconvection is a hot subject for scientists and researchers. Some published works and their few findings are given in Table 1.
The study of applied magnetic field plays a significant role in heat and mass transfer, controlling momentum in boundary layer nanofluid flow through a shrinking/stretching surface. Moreover, magnetohydrodynamics is utilized in diverse engineering applications like cooling of nuclear reactor, electromagnetic casting and plasma confinement. In microfluidics complexities, magnetohydrodynamics is helpful to create non-pulsating and continuous movement in the structure of complex micro-channels. In biomedical engineering, magnetohydrodynamics is applicable for magnetic drug targeting in cancer diseases. On behalf of these applications, Khan et al. [17] explored that due to the strong applied magnetic field, flow becomes slow and temperature becomes high. They also found the residual errors for the authentication of the solution. Palwasha et al. [18] discussed the intensified magnetic field which causes the Hall current effect in two dimensions case. They found that magnetic field effect is a highly dominant factor to the flow. Khan et al. [19] presented the study to show the effects of applied magnetic field on graphene nanoparticles in a thin film Eyring Powell nanofluid involving velocity slip condition. Their results reveal that the thermal conductivity increases with increasing the concentration of graphene nanoparticles. Zuhra et al. [20] also provided the study on graphene nanoparticles in two dimensional unsteady flow and heat transfer in thin liquid film non-Newtonian (Casson and Williamson) nanofluids suspension past a stretching sheet in the presence of transverse magnetic field and non-uniform heat source/sink. Their study shows that the temperature increases with temperature dependent heat source sink parameter. Khan et al. [21] analyzed the magnetohydrodynamic study in cylindrical coordinates for nanoliquid thin film using spraying phenomenon in thermal system. Apart from various effects of different parameters on different profiles, they also discussed the pressure role in the system. Focusing on the MHD flow, Khan et al. [22] presented Brownian motion and thermophoresis effects on MHD mixed convective thin film second-grade nanofluid flow with Hall effect and heat transfer past a stretching sheet, manifesting that the flow becomes three dimensional due to Hall effect.
The study of motion through a porous space is also very important in engineering and industrial processes. A porous media is a material having pores, i.e., oxides and metals, where these pores are covered with a fluid. Such type of motion has two main benefits, as its area of dissipation is higher as compared to conventional fins that improve the heat convection. Further, the regular flow of a fluid along the whole individual beads blends the fluid effectively. The nanofluid flow through a porous media reveals enhanced thermal properties, i.e., convective heat transfer coefficients as compared with a base material and higher thermal conductivity. Consequently, a porous media having nanoparticles plays a vital role in the enhancement of heat transfer and this mechanism is very important from an application and theoretical point of view. Such types of flow can be seen in filtration, oil flow, heat exchangers, drying process and geothermal systems. Moreover, nanofluid flow through a porous media has many applications in thermal engineering which are favorable to improve the thermal properties of lubricants and oil, electronic cooling equipment and vehicle cooling. Applying Navier-Stokes equations to the physical problems, Khan et al. [23] examined the important role of flow and transfer in a second-grade fluid through a porous medium past a stretching sheet in which the porous medium carries a leading role. Their results show a close agreement with the published work. It is very important to explain how nanofluid flow is possible in a porous medium. Thus, without special precautions, nanoparticles will be simply absorbed by the porous matrix. Basically, the porous matrix will work as a filter for nanoparticles. This situation has been described, explained and modeled in the paper by Wu et al. [24]. The physical situation described in this paper reveals that the work on porous media filled by nanofluids is not just a mathematical challenge, but is based on a deep physical understanding of nanofluid flows. This demonstrates that we are simulating here a real physics problem of mixed convection flow and heat transfer in a porous medium filled by a nanofluid. Kuznetsove [25] has compiled the studies on the developments of bioconvection in porous media caused by either gyrotactic or oxytactic microorganisms, with the limitation of average pore size is larger than the size of a microorganism, which allows the smooth movement of the microorganisms. Tham et al. [26] numerically studied the steady mixed convection boundary layer flow on a horizontal circular cylinder with a constant surface temperature embedded in a porous medium saturated by a nanofluid containing both nanoparticles and gyrotactic microorganisms in a stream flowing vertically upwards for both cases of a heated and cooled cylinder. The objective of that paper is to understand the behavior of nanofluid suspensions containing nanoparticles, gyrotactic microorganisms and entropy generation at the fundamental level.
Motivated from the above discussion, the present study discusses the entropy generation and mixed convection in gravity-driven MHD second grade nanofluid flow containing both nanoparticles and gyrotactic microorganisms through a porous medium with chemical reaction along a convectively heated vertical surface. Employing suitable transformations, the basic governing equations of the problem are converted to dimensionless form which have been solved using a powerful analytic tool called Homotopy Analysis Method (HAM) [27]. The influences of all the pertinents parameters on velocity, temperature, nanoparticle concentration, density motile microorganisms and entropy generation profiles have been shown graphically and illustrated.

2. Methods

Basic Equations

A motion of a magnetohydrodynamic two dimensional, time independent, laminar and an incompressible second grade nanofluid falling downwards through a porous medium along a vertical solid surface due to gravity is considered. The uniform incoming flow on the right side of the plate has a constant temperature T at x = 0 and the nanofluid on the left side of the plate has another constant temperature T f . To avoid the bioconvection instability, it is assumed that the nanofluid is dilute. In addition, the assumption is taken for the stability of the nanoparticles suspended in the base fluid so that the nanoparticles do not agglomerate in the fluid. The passive control of nanoparticles volume fraction at the boundary (on the solid wall) is taken. It is assumed that the microorganisms have constant distributions on the wall. It is most important to note that the base fluid is water so that the microorganisms can survive. The assumption is also maintained that the existence of nanoparticles has few effects on the motion of the microorganisms. A transverse magnetic field of strength B 0 is applied in the positive y-direction, normal to the plate surface. There is no applied voltage and the magnetic Reynolds number is small, so the induced magnetic field and Hall effects are negligible. The physical model and coordinate system are shown in Figure 1.
Applying the above-mentioned assumptions, the equations for the nano-bioconvection model are written in the following form as in [4,13,14,15]
u x + v y = 0 ,
ρ f u u x + v u y = ρ f U d U d x + μ f 2 u y 2 + α 1 x u 2 u y 2 u y 2 u x y + v 3 u y 3 + ( 1 C ) ρ f g β ( T T ) ( ρ p ρ f ) g ( C C ) + g ( N N ) γ a v Δ ρ σ B 0 2 u ν f k ( u U ) k F k ( u 2 U 2 ) ,
u T x + v T y = λ 2 T y 2 + τ D B C y T y + D T T T y 2 ,
u C x + v C y = D B 2 C y 2 + D T T 2 T y 2 k r ( C C ) ,
u N x + v N y + ( N v s ) y = D n 2 N y 2 .
Here
d p d x = ρ f U d U d x , U = ( 2 a x ) 1 2 ,
where U is the free stream velocity, a is constant, u and v are the velocity components in the x-and y-directions. α 1 is the normal stress moduli, v s = b W c Δ C C y is the average swimming velocity vector of the oxytactic microorganisms in which b is the chemotaxis constant and W c is the maximum cell swimming speed. The subscripts p, f and f denote respectively the solid particles, the nanofluid and the base fluid at far field. Δ ρ = ρ c e l l ρ f is the density difference between a cell and base fluid density at far field, μ f is the dynamic viscosity, γ a v is the average volume of microorganisms, σ is the electrical conductivity and ρ f is the density of the nanoliquid. ν f = μ f ρ f is the kinematic viscosity, β is the coefficient of volumetric volume expansion of a second grade nanofluid, g is the acceleration due to gravity, C is the nanopartical volume fraction, k = k 0 x is the permeability of porous medium and k F = k 0 F x is the Forchheimer resistance factor where k 0 is the initial permeability. N is the number density of motile microorganisms, C is the ambient nanofluid volume fraction, λ = k 1 ρ f is the thermal diffusivity of the nanofluid in which k 1 is the thermal conductivity, k r is the rate (constant) of chemical reaction, τ = ( ρ c ) p ( ρ c ) f is the ratio of nanoparticle heat capacity and the base fluid heat capacity, D B is the Brownian diffusion coefficient, D n is the diffusivity of microorganisms, T is the ambient temperature, D T is the thermophoretic diffusion coefficient and T is the temperature inside the boundary layer.
The boundary conditions for the nano-bioconvection flow are
u = 0 , v = 0 , k 1 T y = h f ( T f T ) , D B d C d y + D T T d T d y = 0 , N = N w a t y = 0 ,
u U ( x ) , u y 0 , T T , C C , N N a s y ,
where h f (x) denotes the heat transfer coefficient due to T f , N w is the wall concentration of microorganisms and N is the ambient concentration of microorganisms. Note that, in order to satisfy the boundary conditions at infinity, we have to set N = 0.
Introducing the transformations for nondimensional variables f, θ , ϕ , Ω and similarity variable ζ as
ψ ( x , y ) = 4 U ν f x 3 1 2 f ( ζ ) , θ ( ζ ) = T T T f T , ϕ = C C C w C , Ω = N N N w N , ζ = 3 U 4 ν f x 1 2 y ,
where ψ is the stream function such that u = ψ y and v = ψ x , x and y are the Cartesian coordinates along surface and normal to it. Equation (9) automatically satisfies mass conservation Equation (1). With the help of Equation (9), the Equations (2)–(5), (7) and (8) yield the following six Equations (10)–(15).
f + 2 3 1 f 2 + f f + γ 1 ( 2 f f f 2 3 f f i v ) + G r θ N r ϕ + R b Ω M f γ 3 f γ 4 f 2 = 0 ,
θ + P r f θ + N t ( θ ) 2 + N b θ ϕ = 0 ,
ϕ + L e f ϕ + N t N b θ γ 5 ϕ = 0 ,
Ω + S c f Ω Pe ( ϕ Ω + ϕ Ω ) = 0 ,
f = f = 0 , θ = γ 2 ( 1 θ ) , Nb ϕ + Nt θ = 0 , Ω = 1 a t ζ = 0 ,
f = 1 , f = 0 , θ = ϕ = Ω = 0 a t ζ = ,
where prime ( ) represents the derivative with respect to ζ and the various parameters in Equations (10)–(14) are explained in Table 2. For γ 1 = 0, the present study corresponds to viscous nanofluid case and if γ 1 = γ 3 = γ 4 = γ 5 = M = 0, then the study [14] is recovered.

3. Entropy Generation

The subject of the present analysis has a close relation with engineering and industrial processes so it is interesting to investigate the rate of entropy generation. The volumetric entropy generation of the fluid can be written as
E g e n = λ T 2 T y 2 + μ T u y 2 + μ α 1 T 2 u u y 2 u x y + v u y 2 u y 2 + σ B 0 2 T u 2 + R D C C y 2 + R D T T x C x + T y C y + R D N N y 2 + R D T T x N x + T y N y ,
where R denotes the ideal gas constant and D represents the diffusivity. In Equation (16), on the right hand side, the first term represents the irreversibility due to heat transfer, the second term is entropy generation due to viscous dissipation, the third term is second grade fluid friction irreversibility, the fourth term is the local entropy generation due to the effect of magnetic field and the fifth to eighth terms are the irreversibilities due to diffusion effect. The characteristic entropy generation rate is
E 0 = λ ( T f T ) 2 L 2 T 2 ,
where L is the characteristic length. The entropy generation E G is the ratio between the entropy generation rate E g e n and the characteristic entropy generation rate E 0 i.e.,
E G = E g e n E 0
Applying Equation (9), the non-dimensional form of Equation (18) for entropy generation number N G ( ζ ) is
N G ( ζ ) = R e ( θ ) 2 + R e B r ( θ w ) 2 ( f ) 2 + R e B r γ 1 ( θ w ) 2 f ( f ) 2 3 f f f + R e B r M ( θ w ) 2 ( f ) 2 + R e γ 6 ϕ w θ w 2 ( ϕ ) 2 + R e γ 6 ( γ 7 + 1 ) ϕ w θ w θ ϕ + R e γ 8 Ω w θ w 2 ( Ω ) 2 + R e γ 8 ( γ 7 + 1 ) Ω w θ w θ Ω ,
where Re = L 2 3 U 4 ν f x 1 4 is the Reynolds number, Br = μ U 2 λ T is the Brinkman number, γ 6 = R D C λ is the diffusive constant parameter due to nanoparticles concentration, γ 7 = ( y 4 x ) 2 is nondimensional positive number, γ 8 = R D N λ is the diffusive constant parameter due to gyrotactic microorganisms concentration, θ w = T f T w T , ϕ w = C w C C and Ω w = N w N N are respectively the temperature, nanoparticles concentration and microorganisms concentration differences. Due to entropy generation rate, the second-grade nanofluid parameter and magnetic field parameter, here become respectively γ 1 = α 1 U 4 T x and M = 4 σ ν f B 0 2 x 2 3 μ f U 2 .

4. Analytical Solution of the Problem by Homotopy Analysis Method

The initial approximations and the auxiliary linear operators for velocity, temperature, nanoparticles concentration and motile microorganism concentration are presented as
f 0 ( ζ ) = ζ exp ( ζ ) exp ( 2 ζ ) , θ 0 ( ζ ) = exp ( ζ ) exp ( 2 ζ ) 2 + γ 2 , ϕ 0 ( ζ ) = N t γ 2 N b ( γ 2 + 2 ) exp ( ζ ) , Ω 0 ( ζ ) = exp ( ζ ) ,
L f = f f , L θ = θ θ , L ϕ = ϕ ϕ , L Ω = Ω Ω
The following properties are related to the linear operators
L f C 1 + C 2 exp ( ζ ) + C 3 exp ( ζ ) = 0 , L θ C 4 exp ( ζ ) + C 5 exp ( ζ ) = 0 , L ϕ C 6 exp ( ζ ) + C 7 exp ( ζ ) = 0 , L Ω C 8 exp ( ζ ) + C 9 exp ( ζ ) = 0 ,
where C i (i = 1–9) are the arbitrary constants.

4.1. Zeroth-Order Deformation Problems

Introducing the nonlinear operator as
f [ f ( ζ , p ) , θ ( ζ , p ) , ϕ ( ζ , p ) , Ω ( ζ , p ) ] = 3 f ( ζ , p ) ζ 3 + 2 3 1 f ( ζ , p ) ζ 2 + f ( ζ , p ) 2 f ( ζ , p ) ζ 2 + γ 1 2 f ( ζ , p ) ζ 3 f ( ζ , p ) ζ 3 2 f ( ζ , p ) ζ 2 2 3 f ( ζ , p ) 4 f ( ζ , p ) ζ 4 + G r θ ( ζ , p ) N r ϕ ( ζ , p ) + R b Ω ( ζ , p ) M f ( ζ , p ) ζ γ 3 f ( ζ , p ) ζ γ 4 f ( ζ , p ) ζ 2 ,
θ [ f ( ζ , p ) , θ ( ζ , p ) , ϕ ( ζ , p ) ] = 2 θ ( ζ , p ) ζ 2 + P r f ( ζ , p ) θ ( ζ , p ) ζ + N t θ ( ζ , p ) ζ 2 + N b θ ( ζ , p ) ζ ϕ ( ζ , p ) ζ ,
ϕ [ f ( ζ , p ) , θ ( ζ , p ) , ϕ ( ζ , p ) ] = 2 ϕ ( ζ , p ) ζ 2 + L e f ( ζ , p ) ϕ ( ζ , p ) ζ + N t N b 2 θ ( ζ , p ) ζ 2 γ 5 ϕ ( ζ , p ) ,
Ω [ f ( ζ , p ) , ϕ ( ζ , p ) , Ω ( ζ , p ) ] = 2 Ω ( ζ , p ) ζ 2 + S c f ( ζ , p ) Ω ( ζ , p ) ζ P e ϕ ( ζ , p ) ζ Ω ( ζ , p ) ζ + 2 ϕ ( ζ , p ) ζ 2 Ω ( ζ , p ) ,
where p is an embedding parameter such that p ∈ [0, 1].
The zeroth-order deformation problems are expressed as
( 1 p ) L f [ f ( ζ , p ) f 0 ( ζ ) ] = p h f [ f ( ζ , p ) , θ ( ζ , p ) , ϕ ( ζ , p ) , Ω ( ζ , p ) ] ,
( 1 p ) L θ [ θ ( ζ , p ) θ 0 ( ζ ) ] = p h θ [ f ( ζ , p ) , θ ( ζ , p ) , ϕ ( ζ , p ) ] ,
( 1 p ) L ϕ [ ϕ ( ζ , p ) ϕ 0 ( ζ ) ] = p h ϕ [ f ( ζ , p ) , θ ( ζ , p ) , ϕ ( ζ , p ) ] ,
( 1 p ) L Ω [ Ω ( ζ , p ) Ω 0 ( ζ ) ] = p h Ω [ f ( ζ , p ) , ϕ ( ζ , p ) , Ω ( ζ , p ) ] ,
where h represents the auxiliary non-zero parameter.
Equations (27)–(30) have the boundary conditions in Equations (31)–(34) respectively
f ( 0 , p ) = 0 , f ( 0 , p ) = 0 , f ( , p ) = 1 .
θ ( 0 , p ) = γ 2 ( 1 θ ( 0 , p ) ) , θ ( , p ) = 0 .
N b ϕ ( 0 , p ) + N t θ ( 0 , p ) = 0 , ϕ ( , p ) = 0 .
Ω ( 0 , p ) = 1 , Ω ( , p ) = 0 .
Assuming p = 0 and p = 1, the following results hold
p = 0 f ( ζ , 0 ) = f 0 ( ζ ) a n d p = 1 f ( ζ , 1 ) = f ( ζ ) ,
p = 0 θ ( ζ , 0 ) = θ 0 ( ζ ) a n d p = 1 θ ( ζ , 1 ) = θ ( ζ ) ,
p = 0 ϕ ( ζ , 0 ) = ϕ 0 ( ζ ) a n d p = 1 ϕ ( ζ , 1 ) = ϕ ( ζ ) ,
Similarly
p = 0 Ω ( ζ , 0 ) = Ω 0 ( ζ ) a n d p = 1 Ω ( ζ , 1 ) = Ω ( ζ ) .
Using Taylor series expansion and Equations (27)–(30), gives the following results
f ( ζ , p ) = f 0 ( ζ ) + m = 1 f m ( ζ ) p m , f m ( ζ ) = 1 m ! m f ( ζ , p ) p m p = 0 ,
θ ( ζ , p ) = θ 0 ( ζ ) + m = 1 θ m ( ζ ) p m , θ m ( ζ ) = 1 m ! m θ ( ζ , p ) p m p = 0 ,
ϕ ( ζ , p ) = ϕ 0 ( ζ ) + m = 1 ϕ m ( ζ ) p m , ϕ m ( ζ ) = 1 m ! m ϕ ( ζ , p ) p m p = 0 ,
Ω ( ζ , p ) = Ω 0 ( ζ ) + m = 1 Ω m ( ζ ) p m , Ω m ( ζ ) = 1 m ! m Ω ( ζ , p ) p m p = 0 .
The convergence of the series is closely related to h. The value of h is taken in such a way that the series in Equations (39)–(42) converge at p = 1, then Equations (39)–(42) result in
f ( ζ ) = f 0 ( ζ ) + m = 1 f m ( ζ ) ,
θ ( ζ ) = θ 0 ( ζ ) + m = 1 θ m ( ζ ) ,
ϕ ( ζ ) = ϕ 0 ( ζ ) + m = 1 ϕ m ( ζ ) ,
Ω ( ζ ) = Ω 0 ( ζ ) + m = 1 Ω m ( ζ ) .

4.2. m-th Order Deformation Problems

By taking m times derivative with respect to p of Equations ((27), (31)), ((28), (32)), ((29), (33)) and ((30), (34)), then dividing by m! and substituting p = 0 in each computation, yield the following m-th order deformation problems
L f [ f m ( ζ ) χ m f m 1 ( ζ ) ] = h R m f ( ζ ) ,
f m ( 0 ) = f m ( 0 ) = f m ( ) = 0 ,
R m f ( ζ ) = f m 1 + 2 3 k = o m 1 1 f m 1 k f k + k = o m 1 f m 1 k f k + γ 1 k = 0 m 1 2 f m 1 k f k f m 1 k f k 3 f m 1 k f k i v + G r θ m 1 N r ϕ m 1 + R b Ω m 1 M f m 1 γ 3 f m 1 γ 4 k = 0 m 1 f m 1 k f k .
L θ [ θ m ( ζ ) χ m θ m 1 ( ζ ) ] = h R m θ ( ζ ) ,
θ m ( 0 ) = θ m ( ) = 0 ,
R m θ ( ζ ) = θ m 1 + P r k = o m 1 f m 1 k θ k + N t k = o m 1 θ m 1 k θ k + N b k = o m 1 θ m 1 k ϕ k .
L ϕ [ ϕ m ( ζ ) χ m ϕ m 1 ( ζ ) ] = h R m ϕ ( ζ ) ,
N b ϕ m ( 0 ) + N t θ m ( 0 ) = ϕ m ( ) = 0 ,
R m ϕ ( ζ ) = ϕ m 1 + L e k = o m 1 f m 1 k ϕ k + N t N b θ m 1 γ 5 ϕ m 1 .
L Ω [ Ω m ( ζ ) χ m Ω m 1 ( ζ ) ] = h R m Ω ( ζ ) ,
Ω m ( 0 ) = Ω m ( ) = 0 ,
R m Ω ( ζ ) = Ω m 1 + S c k = o m 1 f m 1 k Ω k P e k = o m 1 ϕ m 1 k Ω k + ϕ m 1 k Ω k ,
χ m = 0 , m 1 1 , m > 1 .
If f m ( ζ ), θ m ( ζ ), ϕ m ( ζ ) and Ω m ( ζ ) are the particular solutions, then the general solutions of Equations (47), (50), (54) and (56) are
f m ( ζ ) = f m ( ζ ) + C 1 + C 2 exp ( ζ ) + C 3 exp ( ζ ) ,
θ m ( ζ ) = θ m ( ζ ) + C 4 exp ( ζ ) + C 5 exp ( ζ ) ,
ϕ m ( ζ ) = ϕ m ( ζ ) + C 6 exp ( ζ ) + C 7 exp ( ζ ) ,
Ω m ( ζ ) = Ω m ( ζ ) + C 8 exp ( ζ ) + C 9 exp ( ζ ) .

5. Results

The non-linear differential Equations (10)–(13) with boundary conditions in Equations (14) and (15) have been solved with the help of symbolic computation software MATHEMATICA employing HAM program. The effects of embedded parameters on the velocity f( ζ ), temperature θ ( ζ ), concentration ϕ ( ζ ), motile gyrotactic microorganisms Ω ( ζ ) fields and entropy generation have been plotted in Figures 6–17, 18–25, 26–35 and 36–46 respectively. Similarly the residual error graphs have been shown in Figures 54–57. The schematic illustration of the problem is seen in Figure 1. Liao [27] introduced h curves for the convergence of the series solutions of the problems. Therefore, the admissible h-curves for f( ζ ), θ ( ζ ), ϕ ( ζ ) and Ω ( ζ ) are drawn in the ranges −1.3 ≤ h ≤ 0.0, −1.8 ≤ h ≤ 0.2, −2.1 ≤ h ≤ 0.1 and −2.0 ≤ h ≤ 0.0 in Figure 2, Figure 3, Figure 4 and Figure 5.

6. Discussion

6.1. Velocity Profile

This section depicts the influence of numerous dimensionless parameters on the velocity f( ζ ). Figure 6 shows that the velocity profile decreases significantly with an enhancement in the values of second grade nanofluid parameter γ 1 . Velocity is minimum for non-Newtonian nanofluids in general. Figure 7 demonstrates that the non-dimensional velocity profile f( ζ ) increases for the greater values of reduced heat transfer parameter γ 2 due to slip condition. Figure 8 reveals that the non-dimensional velocity profile f( ζ ) depreciates for the high values of porosity parameter γ 3 . The holes of the medium creates resistance to the flow thereby decreasing the motion. Figure 9 is prepared for the inertial parameter γ 4 . An increment in this parameter decreases the velocity f( ζ ). In Figure 10, it is seen that the non-dimensional velocity f( ζ ) amplifies with increasing magnitude of chemical reaction parameter γ 5 which shows that velocity enhances with the strengthening of chemical reaction. Figure 11 shows that the non-dimensional velocity f( ζ ) increases with the increasing magnitude of buoyancy parameter Gr. Gravitational force favors the second grade nanofluid flow along the vertical solid surface. The buoyancy ratio parameter Nr is showing its effect in Figure 12. An enhancement in the parameter Nr makes a sharp advancement in the nanofluid flow so velocity f( ζ ) becomes high. In Figure 13, the velocity f( ζ ) decreases for increasing quantities of bioconvection Rayleigh number Rb which characterizes the effect of up-swimming of microorganisms. It should be noted that there is no effect of Rb on the boundary layer thickness. For the case of pure natural convection, Rb = 0 which implies that microorganisms have the same density as water and, therefore, they do not contribute to density stratification (or which corresponds to the situation with zero concentration of microorganisms). Since the microorganisms swim on average in the upward direction, their presence produces a destabilizing effect and reduces the critical value of Rb. By increasing the average concentration of corresponding microorganisms in the suspension, Rb can be increased. Figure 14 depicts the role of Lewis number Le. Velocity f( ζ ) is found to be maximized with greater values of Le. Nanoparticles concentrations are consistently depressed with greater lewis number; nevertheless the over-riding influence is imposed by the lewis number on the diffusion of nanoparticles. The role of Schmidt number Sc is projected in Figure 15 which shows that velocity f( ζ ) amplifies with the effective diffusion of microorganisms and the acceleration of nanoparticles. Figure 16 demonstrates that the velocity f( ζ ) decelerates with the rising values of magnetic field parameter M. It is due to the fact that the magnetic field exerts a force, known as Lorentz force which suppresses the velocity. Figure 17 illustrates that the velocity f( ζ ) enhances by the high values of the Prandtl number Pr. It is evident that high quantities of Pr make enhancement in the mass transfer rate so f( ζ ) elevates due to gravity.

6.2. Temperature Profile

Temperature profile exhibits the interesting behavior for various parameters. Figure 18 shows the influence of second grade nanofluid parameter γ 1 . It demonstrates that the temperature θ ( ζ ) increases with the non-Newtonian effect of the nanofluid. The thermal slip effect is led by reduced heat transfer parameter γ 2 and is exhibited in Figure 19. The graph projects that with the elevation of reduced heat transfer parameter γ 2 , temperature θ ( ζ ) diminishes consequently the decrement comes into play in the thermal boundary layer thickness. Figure 20 reveals that as the porosity parameter γ 3 increases, the surface area of the porous media increases which drops down the temperature. Figure 21 depicts that as the values of inertial parameter γ 4 increases, the temperature diminishes in the boundary region. Figure 22 clearly shows that temperature increases with increasing chemical reaction parameter γ 5 . It is seen that temperature θ ( ζ ) increases with the sharpness of chemical reaction. Indeed, some chemical reactions like fissions and fusions make the temperature high. The thermophoresis parameter Nt gives rise in temperature with larger values of Nt as demonstrated by Figure 23. The thermophoresis parameter Nt is directly proportional to the heat transfer associated with the hot nanofluid. Thus an increment in the thermophoresis parameter Nt causes reduction in the heat transfer at the surface, consequently the temperature at the surface of the plate increases. Physically, it is due to the fact that with an enhancement of the thermophoresis parameter Nt, the random motion of the particles increases, resulting in an increase in the temperature profile. Due to the Brownian motion, the particles are much closer and so enhance the heat flow among the particles. Figure 24 depicts that the nondimensional temperature θ ( ζ ) decreases for the increasing values of magnetic field parameter M. It is due to the gravity because temperature rises at height and it falls down in tending towards ground. The influence of Prandtl number Pr has been projected in Figure 25 showing that the temperature θ ( ζ ) diminishes with the increasing values of Pr. The reason is that the increasing values of Prandtl number increase the viscosity of the fluid which absorbs the heat as a result, cooling phenomenon is generated.

6.3. Nanoparticle Concentration Profile

Nanoparticles having one or more dimensions of the order 100 nm or less have attracted great attention due to their unusual and fascinating properties and applications. Figure 26 displays that the concentration decreases with the non-Newtonian second grade nanofluid parameter γ 1 . The non-Newtonian second grade nanofluid is saturated so there is less opportunity for the diffusion of the nanoparticles, therefore, the concentration decreases. Figure 27 shows the effect of reduced heat parameter γ 2 on concentration ϕ ( ζ ). It is a common phenomenon that concentration depreciates with decreasing temperature so in the present case, the concentration ϕ ( ζ ) falls down with decreasing the reduced heat parameter γ 2 . Figure 28 focuses on the role of porosity parameter γ 3 . It provides that the concentration is enhanced with the rising values of porosity parameter γ 3 . The performance of inertial parameter γ 4 is shown in Figure 29 against the non-dimensional concentration profile ϕ ( ζ ). Concentration makes development with inertia since saturation is easily received. The highlights of the impact of chemical reaction parameter γ 5 is revealed in Figure 30. It is reported that for constructive (or generation) chemical reaction, the mass transfer shows an increasing behavior i.e., it is seen that the concentration profile increases significantly for the larger values of the chemical reaction parameter γ 5 . Figure 31 projects the role of Brownian motion parameter Nb. Due to Brownian motion, the random motion of the particles increases, which favors the non-dimensional nanoparticle concentration profile ϕ ( ζ ). Figure 32 illustrates that by enhancing the thermophoresis parameter Nt, the concentration profile ϕ ( ζ ) depreciates. The nanopartical concentration boundary layer thickness is an increasing function of the thermophoresis parameter Nt. This indicates that an increment in the thermophoresis parameter Nt induces resistance to the diffusion of the solute and this assists in the reduction of the concentration gradient at the surface. It is interesting to note that the effect of the thermophoresis parameter is more pronounced for the nanoparticle volume fraction than on the fluid velocity and temperature profiles. Figure 33 is specified for the effect of Lewis number Le on the concentration field ϕ ( ζ ) which reveals that the concentration becomes rich for the different values of the Lewis number Le. The fact is that the performance of the nanoparticles movement is developed, which enhances the nanoparticles concentration. An increment in the Lewis number Le causes the reduction in the mass diffusivity of the nanofluid which in turn reduces the solute concentration. Figure 34 shows that the nondimensional concentration field ϕ ( ζ ) increases as the magnetic field parameter M increases. It is due to the existence of Lorentz force which shows resistivity to the over all motion, consequently, nanoparticle concentration increases. The magnetic field parameter increases the friction between nanofluid layers so the concentration boundary layer thickness gets enhanced. The influence of Prandtl number Pr on the nanoparticle concentration profile ϕ ( ζ ) is depicted in Figure 35. Prandtl number is the fluid viscosity parameter so by increasing Prandtl number, the internal viscous forces of the fluid is strong and the concentration goes to peak which causes no special capability for further nanoparticles concentration. Therefore, concentration ϕ ( ζ ) becomes weak with larger values of Pr.

6.4. Gyrotactic Microorganism Concentration

There have been tremendous developments in the field of microorganism produced nanoparticles and their applications. The characteristics of the density motile microorganisms with respect to different parameters are discussed in the following. Figure 36 elucidates that the microorganism concentration profile Ω ( ζ ) grows up with increasing magnitude of second grade nanofluid parameter γ 1 . The reason is that increasing the quantity γ 1 increases the microorganism propulsion by providing the best platform for the growth of microorganism and hence explains the reason why the microorganism concentration boundary layer of second grade nanofluid enhances. Figure 37 shows the influence of reduced heat parameter γ 2 on the microorganism concentration Ω ( ζ ) profile. Intensity of temperature is not good for the survival of microorganism concentration. Therefore, in the present case, the microorganism concentration Ω ( ζ ) increases with moderate temperature. Figure 38 projects that the the microorganism concentration Ω ( ζ ) decreases with amplifying the porosity parameter γ 3 . In porous medium there is a resistance of holes to the flow of microorganisms so their progress is not significant in that case. Figure 39 illustrates that the non-dimensional microorganism concentration Ω ( ζ ) profile depreciates with the rising values of inertial parameter γ 4 . This is due to the fact that the inertial parameter affects the motion of the microorganisms since the environment is not clear to them. Figure 40 reveals that microorganism concentration field Ω ( ζ ) amplifies with the chemical reaction parameter γ 5 . Figure 41 shows that the microorganism concentration Ω ( ζ ) decreases against the Brownian motion parameter Nb because the propulsion of microorganisms is weak. Figure 42 demonstrates that by increasing thermophoresis parameter Nt, the microorganism concentration profile Ω ( ζ ) increases. This is due to the close relation of concentration with moderate temperature and it is expected that a higher thermophoresis parameter allows a deeper penetration of the microorganism concentration. This achieves enhancement in species in the nanofluid and elevates the motile microorganism concentration profile Ω ( ζ ). This indicates that a stronger thermophoresis effect induces a resistance in the mass diffusion. In Figure 43, the microorganism concentration field Ω ( ζ ) shows an upshot for the greater values of Lewis number Le. This happens because as the Lewis number Le increases, the viscous diffusion rate increases which in turn reduces the velocity at the surface and so elevates the density of the microorganisms. The effect of Schmidt number Sc on the microorganism concentration profile Ω ( ζ ) is projected in Figure 44. When the magnitude of Schmidt number Sc increases, the microorganism concentration profile Ω ( ζ ) decreases. This is attributed to the low propulsion of microorganisms. Figure 45 shows the effect of bioconvection Peclet number Pe on the microorganism concentration profile Ω ( ζ ). It indicates that the microorganism concentration profile Ω ( ζ ) grows large for different values of bioconvection Peclet number Pe. This is due to the fact that microorganisms diffuse more and more with larger values of bioconvection Peclet number Pe. The bioconvection Peclet number Pe helps to increase the speed of the microorganisms in respect of the nanofluid and so the density of the microorganisms is reduced near the surface. The behavior is highly influenced with an increase in the Lewis number Le which again helps in the reduction of the mass diffusivity. Figure 46 shows that by increasing magnetic field parameter M, the nondimensional motile density function profile Ω ( ζ ) decreases because the magnetic field opposes the flow due to Lorentz force.

7. Entropy Generation Analysis

It is necessary to study the irreversibility of the present system in terms of second-grade nanofluid, nanoparticles and microorganisms. Figure 47 shows that the irreversibility N G ( ζ ) of the system increases for high values of Reynolds number Re. Similarly, the values of Br make enhancement in entropy generation rate N G ( ζ ) which has the same increasing behavior in Figure 48. Different values of magnetic field parameter M in Figure 49 presents an increasing effect on the entropy generation rate N G ( ζ ). In Figure 50, the nondimensional diffusive constant parameter γ 6 assumes various high values (ranging from 0.50 to 1.20) with deep effects observations. γ 6 is related to the diffusivity. It shows that entropy generation rate N G ( ζ ) increases with increasing γ 6 where at the same time the temperature difference parameter θ w reveals a unique activity, including the depreciation of entropy generation rate N G ( ζ ) in Figure 51. The formation of curves of the entropy generation rate N G ( ζ ) in Figure 52 is due to the various values of the concentration difference parameter ϕ w , a reactive species making increment in N G ( ζ ). In Figure 53, it is shown that microorganisms concentration difference parameter Ω w favors the entropy generation rate.

8. Residual Errors

Residual errors are shown in the graphs to show the authentication of the problem. By fixing all the others parameters, different values are assigned to convergence control parameter h to show the effects in Figure 54, Figure 55, Figure 56 and Figure 57.

9. Conclusions

This article discusses the analytical solution of entropy generation and mixed convection in gravity-driven non-Newtonian second grade nanofluid flow containing both nanoparticles and gyrotactic microorganisms through a porous medium with a chemical reaction along a convectively heated vertical solid surface. The thin film considered in this work contains the mixture of nanoparticles and motile gyrotactic microorganisms [14]. The solution of the problem has been obtained by using analytical technique called Homotopy Analysis Method (HAM) for the velocity, temperature, concentration and microorganism concentration fields, which has been displayed in the diagrams. From Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20, Figure 21, Figure 22, Figure 23, Figure 24, Figure 25, Figure 26, Figure 27, Figure 28, Figure 29, Figure 30, Figure 31, Figure 32, Figure 33, Figure 34, Figure 35, Figure 36, Figure 37, Figure 38, Figure 39, Figure 40, Figure 41, Figure 42, Figure 43, Figure 44, Figure 45, Figure 46, Figure 47, Figure 48, Figure 49, Figure 50, Figure 51, Figure 52 and Figure 53, each figure displays the profound effect of the emerging parameters on the flow, heat transfer, nanoparticles concentration, microorganism concentration and entropy generation. Residual errors sketches are drawn in Figure 54, Figure 55, Figure 56 and Figure 57 to show the authentication of the HAM solution. The main findings of the study are summarized as following.
(i)
The velocity f( ζ ) depreciates for the porosity parameter γ 3 , inertial parameter γ 4 , bioconvection Rayleigh number Rb and magnetic field parameter M while it elevates for the second grade nanofluid parameter γ 1 , reduced heat transfer parameter γ 2 , chemical reaction parameter γ 5 , buoyancy parameter Gr, buoyancy ratio parameter Nr, Lewis number Le, Schmidt number Sc and Prandtl number Pr.
(ii)
The temperature θ ( ζ ) diminishes for the reduced heat transfer parameter γ 2 , porosity parameter γ 3 , inertial parameter γ 4 , magnetic field parameter M and Prandtl number Pr while it elevates for the second grade nanofluid parameter γ 1 , chemical reaction parameter γ 5 and thermophoresis parameter Nt.
(iii)
The nanoparticles concentration ϕ ( ζ ) diminishes for second grade nanofluid parameter γ 1 , reduced heat transfer parameter γ 2 , thermophoresis parameter Nt and Prandtl number Pr while it elevates for the porosity parameter γ 3 , inertial parameter γ 4 , chemical reaction parameter γ 5 , Brownian motion parameter Nb, Lewis number Le and magnetic field parameter M.
(iv)
The microorganism concentration Ω ( ζ ) diminishes for the porosity parameter γ 3 , inertial parameter γ 4 , Brownian motion parameter Nb, Schmidt number Sc and magnetic field parameter M while it elevates for the second grade nanofluid parameter γ 1 , reduced heat parameter γ 2 , chemical reaction parameter γ 5 , Lewis number Le, thermophoresis parameter Nt and bioconvection Peclet number Pe.
(v)
Entropy generation rate N G ( ζ ) diminishes with temperature difference parameter θ w while it elevates for Reynolds number Re, Brinkman number Br, magnetic field parameter M, diffusive constant parameter γ 6 , nanoparticles concentration difference parameter ϕ w and microorganism concentration difference parameter Ω w .
(vi)
Residual errors graphs are self explanatory for the efficiency of HAM solution.

Author Contributions

Conceptualization, N.S.K.; methodology, N.S.K.; software, Z.S.; validation, S.I.; formal analysis, I.K.; investigation, I.T.; resources, N.S.K.; data curation, I.T.; writing—original draft preparation, N.S.K.; writing—review and editing, T.A.A.; visualization, T.A.A.; supervision, N.S.K.

Funding

This research received no external funding.

Acknowledgments

The authors are thankful to the Higher Education Commission (HEC) Pakistan for providing the technical and financial support. The authors also express their deep and hearty thanks to the honorable reviewers for useful comments and information on the literature.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
xx-axis coordinate (m)
yy-axis coordinate (m)
uVelocity component along x-axis (m s 1 )
vVelocity component along y-axis (m s 1 )
v s Average swimming velocity (m s 1 )
UFree stream velocity (m s 1 )
W c Maximum cell swimming speed (m s 1 )
aConstant
bChemotaxis constant
B 0 Magnetic flux density (Tesla)
GrThermal Grashof number
MMagnetic field parameter
NrBuoyancy ratio parameter
RbBioconvection Rayleigh number
NtThermophoresis parameter
NbBrownian motion parameter
LeLewis number
ScSchmidt number
PeBioconvection Peclet number
PrPrandtl number
gGravitational acceleration (m s 2 )
PPressure (kg m 1 s 2 )
c p Specific heat at constant pressure (J kg 1 K 1 )
kThermal conductivity (W/m K)
TTemperature (K)
T f Convective surface temperature (K)
T Ambient fluid temperature (K)
h f Heat transfer coefficient
CNanoparticles concentration
C Ambient fluid concentration
NNumber density of motile microorganisms
N w Wall concentration of microorganisms
N Ambient concentration of microorganisms
N G Entropy generation number
DDiffusivity
D B Brownian diffusion coefficient
D T Thermophoretic diffusion coefficient
D n Diffusivity of microorganisms
f( ζ )Dimensionless velocity
LCharacteristic length (m)
RIdeal gas constant
ReReynolds number
BrBrinkman number
α 1 Normal stress moduli
σ Electrical conductivity (( Ω ·m) 1 )
μ Coefficient of viscosity (kg m 1 s 1 )
ρ Density (kg m 3 )
ν Kinematic viscosity (m 2 s 1 )
ψ Physical stream function (m 2 s 1 )
β Coefficient of volumetric volume expansion
Δ Difference operator
λ Thermal diffusivity of nanofluid (m 2 s 1 )
τ Ratio of the heat capacitances of nanoparticle and base fluid
ζ A scaled boundary layer coordinate
θ ( ζ )Dimensionless temperature
θ w Dimensionless temperature difference
ϕ ( ζ )Dimensionless concentration
ϕ w Dimensionless concentration difference
Ω ( ζ )Dimensionless microorganisms concentration
Ω w Dimensionless microorganisms concentration difference
γ a v Average volume of microorganisms (m 3 )
γ 1 Dimensionless second grade nanofluid parameter
γ 2 Reduced heat transfer parameter
γ 3 Porosity parameter
γ 4 Inertial parameter
γ 5 Chemical reaction parameter
γ 6 Diffusive constant parameter due to nanoparticles concentration
γ 7 Non-dimensional positive number
γ 8 Diffusive constant parameter due to microorganisms concentration
avAverage
BBrownian
cCell
pSolid particles
rReaction
nProperties related to microorganisms
fBase fluid
oOrigin
xLocal value
wProperties at the wall
Fluid properties at ambient condition
Superscripts
sSwimming
Differentiation w. r. t. ζ

References

  1. Khan, N.S.; Zuhra, Z.; Shah, Z.; Bonyah, E.; Khan, W.; Islam, S. Hall current and thermophoresis effects on magnetohydrodynamic mixed convective heat and mass transfer thin film flow. J. Phys. Commun. 2018. [Google Scholar] [CrossRef]
  2. Ishaq, M.; Ali, M.; Shah, Z.; Islam, S.; Muhammad, S. Entropy generation on nanofluid thin film flow of Eyring-Powell fluid with thermal radiation and MHD effect on an unsteady porous stretching sheet. Entropy 2018, 6. [Google Scholar] [CrossRef]
  3. Hayat, T.; Khan, M.I.; Qayyum, S.; Khan, M.I.; Alsaedi, A. Entropy generation for flow of Sisko fluid to rotating disk. J. Mol. Liquids 2018, 264, 375–385. [Google Scholar] [CrossRef]
  4. Khan, N.S.; Gul, T.; Khan, M.A.; Bonyah, E.; Islam, S. Mixed convection in gravity-driven thin film non-Newtonian nanofluids flow with gyrotactic microorganisms. Results Phys. 2017, 7, 4033–4049. [Google Scholar] [CrossRef]
  5. Zuhra, S.; Khan, N.S.; Islam, S. Magnetohydrodynamic second grade nanofluid flow containing nanoparticles and gyrotactic microorganisms. Comput. Appl. Math. 2018, 37, 6332–6358. [Google Scholar] [CrossRef]
  6. Khan, N.S. Bioconvection in second grade nanofluid flow containing nanoparticles and gyrotactic microorganisms. Braz. J. Phys. 2018, 43, 227–241. [Google Scholar] [CrossRef]
  7. Palwasha, Z.; Islam, S.; Khan, N.S.; Ayaz, H. Non-Newtonian nanoliquids thin film flow through a porous medium with magnetotactic microorganisms. Appl. Nanosci. 2018, 8, 1523–1544. [Google Scholar] [CrossRef]
  8. Raees, A.; Xu, H.; Liao, S.J. Unsteady mixed nano-bioconvection flow in a horizontal channel with its upper plate expanding or contracting. Int. J. Heat Mass Transf. 2015, 86, 174–182. [Google Scholar] [CrossRef]
  9. Zuhra, S.; Khan, N.S.; Shah, Z.; Islam, S.; Bonyah, E. Simulation of bioconvection in the suspension of second grade nanofluid containing nanoparticles and gyrotactic microorganisms. AIP Adv. 2018, 8, 105210. [Google Scholar] [CrossRef]
  10. Pedley, T.J. Instability of uniform microorganism suspensions revisited. J. Fluid Mech. 2010, 647, 335–359. [Google Scholar] [CrossRef]
  11. Xu, H.; Pop, I. Mixed convection flow of a nanofluid over a stretching surface with uniform free stream in the presence of both nanoparticles and gyrotactic microorganisms. Int. J. Heat Mass Transf. 2014, 75, 610–623. [Google Scholar] [CrossRef]
  12. Aziz, A.; Khan, W.A.; Pop, I. Free convection boundary layer flow past a horizontal flat plate embedded in porous medium filled by nanofluid containing gyrotactic microorganisms. Int. J. Therm. Sci. 2012, 56, 48–57. [Google Scholar] [CrossRef]
  13. Bhatti, M.M.; Mishra, S.R.; Abbas, T.; Rashidi, M.M. A mathematical model of MHD nanofluid flow having gyrotactic microorganisms with thermal radiation and chemical reaction effects. Neural Comput. Appl. 2016. [Google Scholar] [CrossRef]
  14. Raees, A.; Xu, H.; Sun, Q.; Pop, I. Mixed convection in gravity-driven nanoliquid film containing both nanoparticles and gyrotactic microorganisms. Appl. Math. Mech. 2015, 36, 163–178. [Google Scholar] [CrossRef]
  15. Ramzan, M.; Chung, J.D. Naeem Ullah, Radiative magnetohydrodynamic nanofluid flow due to gyrotactic microorganisms with chemical reaction and non-linear thermal radiation. Int. J. Mech. Sci. 2017. [Google Scholar] [CrossRef]
  16. Anjalidevi, S.P.; Kandasamy, R. Effects of chemical reaction, heat and mass transfer on laminar flow along a semi infinite horizontal plate. Heat Mass Transf. 1999, 35, 465–467. [Google Scholar] [CrossRef]
  17. Khan, N.S.; Gul, T.; Islam, S.; Khan, W. Thermophoresis and thermal radiation with heat and mass transfer in a magnetohydrodynamic thin film second-grade fluid of variable properties past a stretching sheet. Eur. Phys. J. Plus 2017, 132. [Google Scholar] [CrossRef]
  18. Palwasha, Z.; Khan, N.S.; Shah, Z.; Islam, S.; Bonyah, E. Study of two-dimensional boundary layer thin film fluid flow with variable thermophysical properties in three dimensions space. AIP Adv. 2018, 8, 105318. [Google Scholar] [CrossRef]
  19. Khan, N.S.; Zuhra, S.; Shah, Z.; Bonyah, E.; Khan, W.; Islam, S. Slip flow of Eyring-Powell nanoliquid film containing graphene nanoparticles. AIP Adv. 2018, 8, 115302. [Google Scholar] [CrossRef]
  20. Zuhra, S.; Khan, N.S.; Khan, M.A.; Islam, S.; Khan, W.; Bonyah, E. Flow and heat transfer in water based liquid film fluids dispensed with graphene nanoparticles. Results Phys. 2018, 8, 1143–1157. [Google Scholar] [CrossRef]
  21. Khan, N.S.; Gul, T.; Islam, S.; Khan, I.; Alqahtani, A.M.; Alshomrani, A.S. Magnetohydrodynamic nanoliquid thin film sprayed on a stretching cylinder with heat transfer. J. Appl. Sci. 2017, 7, 271. [Google Scholar] [CrossRef]
  22. Khan, N.S.; Gul, T.; Islam, S.; Khan, A.; Shah, Z. Brownian motion and thermophoresis effects on MHD mixed convective thin film second-grade nanofluid flow with Hall effect and heat transfer past a stretching sheet. J. Nanofluids 2017, 6, 812–829. [Google Scholar] [CrossRef]
  23. Khan, N.S.; Gul, T.; Islam, S.; Khan, W.; Khan, I.; Ali, L. Thin film flow of a second-grade fluid in a porous medium past a stretching sheet with heat transfer. Alex. Eng. J. 2017. [Google Scholar] [CrossRef]
  24. Wu, M.; Kuznetsov, A.V.; Jasper, W.J. Modeling of particle trajectories in an electro-statically charged channel. Phys. Fluids 2011, 22, 043301. [Google Scholar] [CrossRef]
  25. Kuznetsov, A.V. Emerging Topics in Heat and Mass Transfer in Porous Media, From Bioengineering and Microelectronics to Nanotechnology; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  26. Tham, L.; Nazar, R.; Pop, I. Steady mixed convection boundary layer flow on a horizontal circular cylinder with a constant surface temperature embedded in a porous medium saturated by a nanofluid containing both nanoparticles and gyrotactic microorganisms. J. Heat Transf. 2013, 135, 102601-1. [Google Scholar] [CrossRef]
  27. Liao, S.J. Homotopy Analysis Method in Nonlinear Differential Equations; Higher Education Press: Beijing, China; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
Figure 1. Geometry of the physical model and coordinates system.
Figure 1. Geometry of the physical model and coordinates system.
Entropy 21 00139 g001
Figure 2. The h curve of f( ζ ).
Figure 2. The h curve of f( ζ ).
Entropy 21 00139 g002
Figure 3. The h curve of θ ( ζ ).
Figure 3. The h curve of θ ( ζ ).
Entropy 21 00139 g003
Figure 4. The h curve of ϕ ( ζ ).
Figure 4. The h curve of ϕ ( ζ ).
Entropy 21 00139 g004
Figure 5. The h curve of Ω ( ζ ).
Figure 5. The h curve of Ω ( ζ ).
Entropy 21 00139 g005
Figure 6. Effect on velocity profile for h = −1.30, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Figure 6. Effect on velocity profile for h = −1.30, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Entropy 21 00139 g006
Figure 7. Effect on velocity profile for h = −1.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Figure 7. Effect on velocity profile for h = −1.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Entropy 21 00139 g007
Figure 8. Effect on velocity profile for h = 0.10, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Figure 8. Effect on velocity profile for h = 0.10, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Entropy 21 00139 g008
Figure 9. Effect on velocity profile for h = 0.10, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Figure 9. Effect on velocity profile for h = 0.10, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Entropy 21 00139 g009
Figure 10. Effect on velocity profile for h = 0.10, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 4 = 0.30, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Figure 10. Effect on velocity profile for h = 0.10, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 4 = 0.30, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Entropy 21 00139 g010
Figure 11. Effect on velocity profile for h = −1.10, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Gr.
Figure 11. Effect on velocity profile for h = −1.10, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Gr.
Entropy 21 00139 g011
Figure 12. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nr.
Figure 12. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nr.
Entropy 21 00139 g012
Figure 13. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Rb.
Figure 13. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Rb.
Entropy 21 00139 g013
Figure 14. Effect on velocity profile for h = −5.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Nr = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Le.
Figure 14. Effect on velocity profile for h = −5.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Nr = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Le.
Entropy 21 00139 g014
Figure 15. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Nr = 0.60, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Sc.
Figure 15. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Nr = 0.60, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Sc.
Entropy 21 00139 g015
Figure 16. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Nr = 0.60, Pe = 1.00, Pr = 10.00 and different values of M.
Figure 16. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Nr = 0.60, Pe = 1.00, Pr = 10.00 and different values of M.
Entropy 21 00139 g016
Figure 17. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, Nr = 0.60, M = 1.00 and different values of Pr.
Figure 17. Effect on velocity profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, Nr = 0.60, M = 1.00 and different values of Pr.
Entropy 21 00139 g017
Figure 18. Effect on temperature profile for h = −0.80, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Figure 18. Effect on temperature profile for h = −0.80, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Entropy 21 00139 g018
Figure 19. Effect on temperature profile for h = −0.80, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Figure 19. Effect on temperature profile for h = −0.80, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Entropy 21 00139 g019
Figure 20. Effect on temperature profile for h = −0.80, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Figure 20. Effect on temperature profile for h = −0.80, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Entropy 21 00139 g020
Figure 21. Effect on temperature profile for h = −0.80, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Figure 21. Effect on temperature profile for h = −0.80, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Entropy 21 00139 g021
Figure 22. Effect on temperature profile for h = −0.90, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 4 = 0.30, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Figure 22. Effect on temperature profile for h = −0.90, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 4 = 0.30, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Entropy 21 00139 g022
Figure 23. Effect on temperature profile for h = −0.90, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Gr = 0.50, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nt.
Figure 23. Effect on temperature profile for h = −0.90, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Gr = 0.50, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nt.
Entropy 21 00139 g023
Figure 24. Effect on temperature profile for h = −0.50, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Gr = 0.50, Sc = 0.70, Pe = 1.00, Pr = 10.00 and different values of M.
Figure 24. Effect on temperature profile for h = −0.50, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Gr = 0.50, Sc = 0.70, Pe = 1.00, Pr = 10.00 and different values of M.
Entropy 21 00139 g024
Figure 25. Effect on temperature profile for h = −0.50, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, Gr = 0.50, M = 1.00 and different values of Pr.
Figure 25. Effect on temperature profile for h = −0.50, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, Gr = 0.50, M = 1.00 and different values of Pr.
Entropy 21 00139 g025
Figure 26. Effect on concentration profile for h = −3.00, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Figure 26. Effect on concentration profile for h = −3.00, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Entropy 21 00139 g026
Figure 27. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Figure 27. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Entropy 21 00139 g027
Figure 28. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Figure 28. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Entropy 21 00139 g028
Figure 29. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Figure 29. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Entropy 21 00139 g029
Figure 30. Effect on concentration profile for h = −1.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 2 = 0.60, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Figure 30. Effect on concentration profile for h = −1.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 2 = 0.60, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Entropy 21 00139 g030
Figure 31. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nb.
Figure 31. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nb.
Entropy 21 00139 g031
Figure 32. Effect on concentration profile for h = −.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nt.
Figure 32. Effect on concentration profile for h = −.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nt.
Entropy 21 00139 g032
Figure 33. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Le.
Figure 33. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Le.
Entropy 21 00139 g033
Figure 34. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, Pr = 10.00 and different values of M.
Figure 34. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, Pr = 10.00 and different values of M.
Entropy 21 00139 g034
Figure 35. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00 and different values of Pr.
Figure 35. Effect on concentration profile for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00 and different values of Pr.
Entropy 21 00139 g035
Figure 36. Effect on microorganisms density function for h = −3.00, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Figure 36. Effect on microorganisms density function for h = −3.00, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 1 .
Entropy 21 00139 g036
Figure 37. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Figure 37. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 2 .
Entropy 21 00139 g037
Figure 38. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Figure 38. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 3 .
Entropy 21 00139 g038
Figure 39. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Figure 39. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 4 .
Entropy 21 00139 g039
Figure 40. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 4 = 0.30, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Figure 40. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 3 = 0.20, γ 2 = 0.60, γ 4 = 0.30, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of γ 5 .
Entropy 21 00139 g040
Figure 41. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nb.
Figure 41. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nb.
Entropy 21 00139 g041
Figure 42. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nt.
Figure 42. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nb = 0.80, Le = 0.60, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Nt.
Entropy 21 00139 g042
Figure 43. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Nb = 0.80, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Le.
Figure 43. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Nb = 0.80, Sc = 0.70, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Le.
Entropy 21 00139 g043
Figure 44. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Nb = 0.80, Nt = 0.90, Le = 0.60, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Sc.
Figure 44. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Nb = 0.80, Nt = 0.90, Le = 0.60, Pe = 1.00, M = 1.00, Pr = 10.00 and different values of Sc.
Entropy 21 00139 g044
Figure 45. Effect on microorganisms density function for h = −1.10, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Sc = 0.70, Nb = 0.80, M = 1.00, Pr = 10.00 and different values of Pe.
Figure 45. Effect on microorganisms density function for h = −1.10, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Sc = 0.70, Nb = 0.80, M = 1.00, Pr = 10.00 and different values of Pe.
Entropy 21 00139 g045
Figure 46. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Sc = 0.70, Nb = 0.80, Pe = 1.00, Pr = 10.00 and different values of M.
Figure 46. Effect on microorganisms density function for h = −3.00, γ 1 = 0.50, γ 2 = 0.60, γ 3 = 0.20, γ 4 = 0.30, γ 5 = 1.00, Gr = 0.50, Nr = 0.60, Rb = 0.70, Nt = 0.90, Sc = 0.70, Nb = 0.80, Pe = 1.00, Pr = 10.00 and different values of M.
Entropy 21 00139 g046
Figure 47. Effect on entropy generation rate for h = −1.50, Br = 0.40, M = 0.50, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of Re.
Figure 47. Effect on entropy generation rate for h = −1.50, Br = 0.40, M = 0.50, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of Re.
Entropy 21 00139 g047
Figure 48. Effect on entropy generation rate for h = −1.50, Re = 0.10, M = 0.50, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of Br.
Figure 48. Effect on entropy generation rate for h = −1.50, Re = 0.10, M = 0.50, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of Br.
Entropy 21 00139 g048
Figure 49. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of M.
Figure 49. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of M.
Entropy 21 00139 g049
Figure 50. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of γ 6 .
Figure 50. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, θ w = 0.70, ϕ w = 0.80, Ω w = 0.90 and different values of γ 6 .
Entropy 21 00139 g050
Figure 51. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, γ 6 = 0.60, ϕ w = 0.80, Ω w = 0.90 and different values of θ w .
Figure 51. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, γ 6 = 0.60, ϕ w = 0.80, Ω w = 0.90 and different values of θ w .
Entropy 21 00139 g051
Figure 52. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, γ 6 = 0.60, θ w = 0.70, Ω w = 0.90 and different values of ϕ w .
Figure 52. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, γ 6 = 0.60, θ w = 0.70, Ω w = 0.90 and different values of ϕ w .
Entropy 21 00139 g052
Figure 53. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80 and different values of Ω w .
Figure 53. Effect on entropy generation rate for h = −1.50, Re = 0.10, Br = 0.40, M = 0.50, γ 6 = 0.60, θ w = 0.70, ϕ w = 0.80 and different values of Ω w .
Entropy 21 00139 g053
Figure 54. Residual errors graph of velocity profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Figure 54. Residual errors graph of velocity profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Entropy 21 00139 g054
Figure 55. Residual errors graph of temperature profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Figure 55. Residual errors graph of temperature profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Entropy 21 00139 g055
Figure 56. Residual errors graph of concentration profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Figure 56. Residual errors graph of concentration profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Entropy 21 00139 g056
Figure 57. Residual errors graph of microorganism concentration profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Figure 57. Residual errors graph of microorganism concentration profile for γ = 0.30, γ 1 = 0.30, γ 2 = 0.30, γ 3 = 0.30, γ 4 = 0.30, Gr = 0.40, Nr = 0.40, Rb = 0.50, Nb = 0.50, Nt = 0.60, Le = 0.60, Sc = 0.40, Pe = 0.40, Pr = 0.70, M = 0.60 and different values of h.
Entropy 21 00139 g057
Table 1. Some recent developments related to the topic under investigation.
Table 1. Some recent developments related to the topic under investigation.
Author NamesAuthor WorksSome Outcomes
Khan et al. [1]Entropy generationEntropy generation increases with Reynolds number
Ishaq et al. [2]Entropy generationEntropy generation decreases with Eyring-Powell parameter
Hayat et al. [3]Entropy generationEntropy generation increases with Reynolds number
Khan et al. [4]BioconvectionBioconvection decreases with reduced heat transfer parameter
Zuhra et al. [5]BioconvectionGyrotactic microorganisms depreciates with magnetic field parameter
Khan [6]BioconvectionStratification increases with second-grade fluid parameter
Palwasha et al. [7]BioconvectionStratification increases with gravitational forces
Raees et al. [8]BioconvectionBioconvection depends on upper plate
Zuhra et al. [9]BioconvectionMicroorganisms decrease with increasing Brownian motion parameter
Pedley [10]InstabilityInstability of the system is due to microorganisms
Xu et al. [11]Mixed convectionNanoparticles favor the mixed convection
Aziz et al. [12]Free convectionInstability of the system increases due to microorganisms
Bhatti et al. [13]Chemical reactionMass transfer increases with chemical reaction
Raees et al. [14]Bioconvection‘Passively controlled nanofluid model provides better results
Ramzan et al. [15]Chemical reactionMass transfer increases with chemical reaction
Anjalidevi [16]Chemical reactionMass transfer increases with chemical reaction
Table 2. Various parameters in Equations (10)–(14).
Table 2. Various parameters in Equations (10)–(14).
Parameter NamesSymbols/NotationsDefined Values
Dimensionless second-grade nanofluid parameter γ 1 α 1 U 4 μ f x
Thermal Grashof numberGr 2 ( 1 C ) ρ f g β ( T f T ) 3 a ρ f
Buoyancy ratio parameterNr 2 ( ρ p ρ f ) C g 3 a ρ f
Bioconvection Rayleigh numberRb 2 ( N w N ) g γ a v Δ ρ 3 a ρ f
Prandtl numberPr ν f λ
Thermophoresis parameterNt τ D T ( T f T ) λ T
Brownian motion parameterNb τ D B C λ
Lewis numberLe ν f D B
Schmidt numberSc ν f D n
Bioconvection Peclet numberPe b W c D n
Reduced heat transfer parameter γ 2 h f k 1 4 x ν f 3 U
Porosity parameter γ 3 4 ν f 3 k 0 U
Inertial parameter γ 4 4 k 0 F 3 k 0
Chemical reaction parameter γ 5 4 k r ν x 3 b D B
Magnetic field parameterM 4 σ B 0 2 3 ρ f U x 2 a

Share and Cite

MDPI and ACS Style

Khan, N.S.; Shah, Z.; Islam, S.; Khan, I.; Alkanhal, T.A.; Tlili, I. Entropy Generation in MHD Mixed Convection Non-Newtonian Second-Grade Nanoliquid Thin Film Flow through a Porous Medium with Chemical Reaction and Stratification. Entropy 2019, 21, 139. https://doi.org/10.3390/e21020139

AMA Style

Khan NS, Shah Z, Islam S, Khan I, Alkanhal TA, Tlili I. Entropy Generation in MHD Mixed Convection Non-Newtonian Second-Grade Nanoliquid Thin Film Flow through a Porous Medium with Chemical Reaction and Stratification. Entropy. 2019; 21(2):139. https://doi.org/10.3390/e21020139

Chicago/Turabian Style

Khan, Noor Saeed, Zahir Shah, Saeed Islam, Ilyas Khan, Tawfeeq Abdullah Alkanhal, and Iskander Tlili. 2019. "Entropy Generation in MHD Mixed Convection Non-Newtonian Second-Grade Nanoliquid Thin Film Flow through a Porous Medium with Chemical Reaction and Stratification" Entropy 21, no. 2: 139. https://doi.org/10.3390/e21020139

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