Nowadays, the study of nanofluids flow has the interest of researches because of its applications in medicine, biochemistry and industrial engineering. Nanofluids are moderately new category of fluids which consist of a base-fluid with nano-sized particles (1 - 100 nm) suspended within them. Choi [1] may be the first author to use the term “nanofluid”, where it was reported that one of the promising nanofluids applications in heat could transfer enhancement. In [2] , Choi et al. showed that the addition of a small amount (less than 1% by volume) of nanoparticles to conventional heat transfer liquids increase the thermal conductivity of the fluid up to approximately two times. An analysis of nanofluids examined by Buongiorno [3] was induced that this massive increase in the thermal conductivity occurs due to the presence of the Brownian diffusion and the thermophoretic diffusion of the nanoparticles. Recent articles on the nanofluid are cited in ReFS [4] [5] [6] [7] [8] .

Peristaltic motion in a channel or tube is considered as a type of flow that has great value in several physiological processes and industries. Peristalsis is a mechanism for mixing and transporting fluids through expansion and contraction of the wave propagation along the channel walls. This mechanism is seen in many biological system such as urine transport from kidney to bladder through the ureter, transport of lymph in the lymphatic vessels, swallowing food through the esophagus, the movement of chime in the gastrointestinal track, ovum movement in the fallopian tube, transport of spermatozoa vasomotion of small blood vessels such as venules and capillaries and blood flow in arteries, transport of corrosive fluids, sanitary fluid transport and blood pumps in heart lung machine etc. This analysis was first investigated by Latham [9] . He presented fluid motion in peristaltic pump and also discussed the characteristic of pressure rise verus flow rate. After the work of Latham [9] , Jaffrin and Shapiro [10] investigated the peristaltic pumping. They made the study under the assumption of long wavelength and low Reynolds number approximation. This work was further extended for Newtonian and non-Newtonian fluids with different flow geometries and boundary conditions as elaborated in the references [11] - [17] .

Recently, Eldabe et al. [18] have discussed peristaltic transport of a magneto non-Newtonian fluid through a porous medium in a horizontal finite channel. Safia Akram [19] has investigated the effects of nanofluid on peristaltic flow of a Carreau fluid model in an inclined magnetic field. Arshad et al. [20] have analyzed the peristaltic transport of a Carreau fluid in a compliant rectangular duct. Khalid Nowar [21] has studied peristaltic flow of a nanofluid under the effect of Hall current and porous medium. Noreen et al. [22] have examined numerical simulation of peristaltic flow of a Carreau naofluid in an asymmetric channel. Tasawar et al. [23] have discussed the radiative peristaltic flow of Jeffery nanofluid with slip conditions and joule heating.

Few attempts have been devoted to peristaltic flows in presence of heat and mass transfer; such investigations are of great importance, which is due to their extensive applications in medical and bio-engineering sciences, as it may be relevant in many processes in human body like oxygenation in lungs, hemodialysis and nutrients diffuse out of blood. Sohail and nadeem [24] have examined the effects of heat and mass transfer on peristaltic flow of Carreau fluid in a vertical annulus. Eldabe et al. [25] have studied the magneto-hydrodynamic flow and heat transfer for a peristaltic motion of Carreau fluid through a porous medium. Eldabe et al. [26] have studied the peristaltic motion of non-Newtonian fluid with heat and mass transfer through a porous medium in channel under uniform magnetic field. Sohail et al. [27] have discussed the effects of heat and mass transfer on peristaltic flow of a nanofluid between eccentric cylinders. Eldabe et al. [28] have discussed the peristaltic pumping of a conduction Sisko fluid through porous medium with heat and mass transfer. Ramesh and Devaker [29] have studied the effects of heat and mass transfer on the peristaltic transport of MHD couple stress fluid through porous medium in a vertical asymmetric channel.

The main aim of this work is to study the peristaltic motion of a Carreau nanofluid with heat and mass transfer through a porous medium in an asymmetric channel under the effects of Hall current, viscous dissipation and Joule heating. The analysis is performed under the well-established long wavelength and low Reynolds number approximations. A detailed mathematical formulation is presented and numerical solution graphically for velocity, temperature, nanoparticle phenomena and pressure gradient have been presented.

Consider the peristaltic flow of an incompressible viscous electrically conducting nanofluid which obeys Carreau model inside a two dimensional vertical asymmetric channel of width d 1 + d 2 through a porous medium. Asymmetry in the channel is produced by propagation of waves along the channel walls traveling with different amplitudes, phases but with constant speed. In the Cartesian coordinates system (X, Y), the right-hand side wall Y = H 1 and the left-hand side wall Y = H 2 are given by

Y = H 1 ( X , t ) = d 1 + a 1 cos ( 2 π λ ( X − c t ) ) (1)

Y = H 2 ( X , t ) = − d 2 − b 1 cos ( 2 π λ ( X − c t ) + ϕ ) (2)

where, a 1 and b 1 are the amplitudes of right and left walls respectively, λ is the wavelength, t is the time, and ϕ is the phase which varies in the range 0 ≤ ϕ ≤ π . When ϕ = 0 then symmetric channel with waves out of the phase can be described (i.e. both walls move outward or inward simultaneously), and for ϕ = π , the waves are in phase. Further d 1 , d 2 , a 1 , b 1 and ϕ satisfy the condition:

a 1 2 + b 1 2 + 2 a 1 b 1 cos ϕ ≤ ( d 1 + d 2 ) 2 , (3)

so that the walls will not intersect with each other.

A strong uniform magnetic field with magnetic flux density B = ( 0 , 0 , B 0 ) is applied and the Hall effects are taken into account. The induced magnetic field is neglected by assuming a very small magnetic Reynolds number ( R e m < < < 1 ) , also it is assumed that there is no applied polarization voltage so that the total electric field E = 0 . The expression for the current density J including the Hall effect and neglecting ion-slip and thermoelectric effects [The generalized Ohm’s law] is given by [30] :

J = σ [ E + V × B − 1 e n e ( J × B ) ] (4)

where σ is the electrical conductivity of the fluid, V is the velocity vector, e is the electric charge of electrons, n e is the number density of electrons. Equation (1) can be solved in J to yield the Lorentz force vector in the form:

J × B = − σ B 0 2 1 + m 2 [ ( U − m V ) i ^ + ( m U + V ) j ^ ] (5)

where U and V are the X and Y components of the velocity vector, m = σ B 0 e n e is

the hall parameter. The heat transfer and nanopatricle processes are maintained by considering temperature T 0 , T 1 and nanoparticle phenomena C 0 , C 1 to the right and left sides wall, respectively.

The constitutive equation for a carreau fluid is given by [24] :

η − η ∞ η 0 − η ∞ = [ 1 + ( Γ γ ˙ ) 2 ] n − 1 2 , (6)

τ i j = η 0 [ 1 + n − 1 2 ( Γ γ ˙ ) 2 ] γ ˙ i j , (7)

in which τ i j is the extra stress tensor, η ∞ is the infinite shear rate viscosity, η 0 is the zero shear rate viscosity, Γ is the time constant, n is the power law index, and γ ˙ is defined as:

γ ˙ = 1 2 ∑ i ∑ j γ ˙ i j γ ˙ j i = 1 2 Π (8)

where Π = t r a c ( g r a d V + ( g r a d V ) T ) 2 is the second invariant strain tensor.

The governing equations for the present problem are described as:

Continuity equation

∂ U ∂ X + ∂ V ∂ Y = 0 (9)

Equations of motion

ρ f ( ∂ U ∂ t + U ∂ U ∂ X + V ∂ U ∂ Y ) = − ∂ P ∂ X + ∂ τ X X ∂ X + ∂ τ X Y ∂ Y − η 0 k 1 U + σ B 0 2 1 + m 2 ( m V − U ) + ρ f g α ( T − T 0 ) + + ρ f g α ( C − C 0 ) (10)

ρ f ( ∂ V ∂ t + U ∂ V ∂ X + V ∂ V ∂ Y ) = − ∂ P ∂ Y + ∂ τ Y X ∂ X + ∂ τ Y Y ∂ Y − η 0 k 1 V + σ B 0 2 1 + m 2 ( m U + V ) (11)

Energy equation

( ρ c ) f ( ∂ T ∂ t + U ∂ T ∂ X + V ∂ T ∂ Y ) = K ( ∂ 2 T ∂ X 2 + ∂ 2 T ∂ Y 2 ) + [ τ X X ∂ U ∂ X + τ X Y ∂ U ∂ Y + τ Y X ∂ V ∂ X + τ Y Y ∂ V ∂ Y ] + σ B 0 2 1 + m 2 ( U 2 + V 2 ) + ( ρ c ) p { D B [ ∂ C ∂ X ∂ T ∂ X + ∂ C ∂ Y ∂ T ∂ Y ] + D T T 0 [ ( ∂ T ∂ X ) 2 + ( ∂ T ∂ Y ) 2 ] } (12)

Concentration equation

∂ C ∂ t + U ∂ C ∂ X + V ∂ C ∂ Y = D B ( ∂ 2 C ∂ X 2 + ∂ 2 C ∂ Y 2 ) + D T T 0 ( ∂ 2 T ∂ X 2 + ∂ 2 T ∂ Y 2 ) (13)

where ρ f is the density of the fluid, P is the pressure, k 1 is the permeability of the porous medium, g is the acceleration due to gravity, α is the volume expansion coefficient, T is the temperature of the fluid, C is the nanoprticle concentration, ( ρ c ) f is the heat capacity of the fluid, K is the thermal conductivity, ( ρ c ) p is the effective heat capacity of the nanoparticle material, D B is the Brownian diffusion coefficient and D T is the thermophoretic diffusion.

Introducing a wave frame ( x , y ) moving with the velocity c away from the fixed frame ( X , Y ) by the transformation

x = X − c t , y = Y , u = U − c , p ( x ) = P ( X , t ) (14)

in which ( x , y ) , ( u , v ) and p are the coordinates, velocity components and pressure in the wave frame.

Defining the following non-dimensional quantities

x ¯ = x λ ,   y ¯ = y d 1 ,   u ¯ = u c ,   v ¯ = v c ,   t ¯ = t λ ,   δ = d 1 λ ,   d = d 2 d 1 ,   a = a 1 d 1 , b = b 1 d 1 ,   h 1 = H 1 d 1 ,   h 2 = H 2 d 2 ,   p ¯ = d 1 2 p η 0 c λ ,   Re = ρ f c d 1 η 0 ,   ψ ¯ = ψ c d 1 , τ ¯ x ¯ x ¯ = λ η 0 c τ x x ,   τ ¯ x ¯ y ¯ = d 1 η 0 c τ x y ,   τ ¯ y ¯ y ¯ = d 1 η 0 c τ y y ,   W e = Γ c d 1 ,   γ ˙ ¯ = d 1 γ ˙ c , s = k 1 d 1 2 ,   M = σ η 0 B 0 d 1 ,   θ ¯ = T − T 0 T 1 − T 0 ,   Ω ¯ = C − C 0 C 1 − C 0 ,   Pr = μ ρ f α , G r = ρ f g α d 1 2 ( T 1 − T 0 ) c η 0 ,   B r = ρ f g α d 1 2 ( C 1 − C 0 ) c η 0 ,   E c = c 1 2 c f ( T 1 − T 0 ) , S c = μ ρ f D B ,   N b = τ D B ( C 1 − C 0 ) α ,   N t = τ D B ( T 1 − T 0 ) α ,   α = K ( ρ c ) f ,   τ = ( ρ c ) p ( ρ c ) f . (15)

where δ is the dimensionless wave number, R e is the Reynolds number, W e is the Weissenberg number, s is the porosity parameter, M is the magnetic parameter, G r is the local temperature Grashof number, B r is the local nanoparticle Grashof number, P r is the Prandtl number, E c is the Eckert number, S c is the Schmidt number, N b is the Brownain motion parameter and N t is the thermophoresis parameter.

Using Equation (14) and the above set of non-dimensional quantities (15) into Equations (10)-(13), the resulting equations in terms of stream function

ψ ( u = ∂ ψ ∂ x , v = − δ ∂ ψ ∂ y ) can be written after dropping bars in the following

Non-dimensionless form as:

Re δ [ ( ∂ ψ ∂ y ∂ ∂ x − ∂ ψ ∂ x ∂ ∂ y ) ∂ ψ ∂ y ] = − ∂ p ∂ x + δ 2 ∂ ∂ x ( τ x x ) + ∂ ∂ y ( τ x y ) − δ m M 2 1 + m 2 ( ∂ ψ ∂ x ) − M 2 1 + m 2 ( ∂ ψ ∂ y + 1 ) − 1 s ( ∂ ψ ∂ y + 1 ) + G r θ + B r Ω (16)

Re δ 3 [ ( − ∂ ψ ∂ y ∂ ∂ x + ∂ ψ ∂ x ∂ ∂ y ) ∂ ψ ∂ x ] = − ∂ p ∂ y + δ 2 ∂ ∂ x ( τ y x ) + δ ∂ ∂ y ( τ y y ) − δ m M 2 1 + m 2 ( ∂ ψ ∂ y − 1 ) − δ 2 M 2 1 + m 2 ( ∂ ψ ∂ y + 1 ) − δ 1 s ( ∂ ψ ∂ x ) (17)

Re δ [ ∂ ψ ∂ y ∂ θ ∂ x − ∂ ψ ∂ x ∂ θ ∂ y ] = 1 P r ( ∂ 2 θ ∂ y 2 + δ 2 ∂ 2 θ ∂ x 2 ) + E c [ ( δ 2 τ x x ∂ 2 ψ ∂ x ∂ y + τ x y ∂ 2 ψ ∂ y 2 − δ 2 τ y x ∂ 2 ψ ∂ x 2 − δ τ y y ∂ 2 ψ ∂ y ∂ x ) ] + E c M 2 1 + m 2 [ δ 2 ( ∂ ψ ∂ x ) 2 + ( ∂ ψ ∂ y + 1 ) 2 ] E c M 2 1 + m 2 [ δ 2 ( ∂ ψ ∂ x ) 2 + ( ∂ ψ ∂ y + 1 ) 2 ] + N t Pr [ δ 2 ( ∂ θ ∂ x ) 2 + ( ∂ θ ∂ y ) 2 ] (18)

Re δ S c [ ∂ ψ ∂ y ∂ Ω ∂ x − ∂ ψ ∂ x ∂ Ω ∂ y ] = ( ∂ 2 Ω ∂ y 2 + δ 2 ∂ 2 Ω ∂ x 2 ) + N t N b ( ∂ 2 θ ∂ y 2 + δ 2 ∂ 2 θ ∂ x 2 ) (19)

where

τ x x = 2 [ 1 + n − 1 2 W e 2 γ ˙ 2 ] ∂ 2 ψ ∂ x ∂ y , (20)

τ x y = [ 1 + n − 1 2 W e 2 γ ˙ 2 ] ( ∂ 2 ψ ∂ y 2 − δ 2 ∂ 2 ψ ∂ x 2 ) = τ y x , (21)

τ y y = 2 δ [ 1 + n − 1 2 W e 2 γ ˙ 2 ] ∂ 2 ψ ∂ x ∂ y , (22)

and,

γ ˙ = [ 2 δ 2 ( ∂ 2 ψ ∂ x ∂ y ) 2 + ( ∂ 2 ψ ∂ y 2 − δ 2 ∂ 2 ψ ∂ x 2 ) 2 + 2 δ 2 ( ∂ 2 ψ ∂ y ∂ x ) 2 ] 1 2 (23)

Under the assumption of long wavelength δ ≪ 1 and low Reynolds number, Equations (16)-(19) become

0 = − ∂ p ∂ x + ∂ ∂ y [ 1 + n − 1 2 W e 2 ( ∂ 2 ψ ∂ y 2 ) 2 ] ∂ 2 ψ ∂ y 2 − ( M 2 1 + m 2 + 1 s ) [ ∂ ψ ∂ y + 1 ] + G r θ + B r Ω (24)

0 = − ∂ p ∂ y (25)

0 = 1 Pr ∂ 2 θ ∂ y 2 + E c [ ( ∂ 2 ψ ∂ y 2 ) 2 + n − 1 2 W e 2 ( ∂ 2 ψ ∂ y 2 ) 4 ] + E c M 2 1 + m 2 ( ∂ ψ ∂ y + 1 ) 2 + N b Pr ( ∂ Ω ∂ y ∂ θ ∂ y ) + N t Pr ( ∂ θ ∂ y ) 2 (26)

0 = ∂ 2 Ω ∂ y 2 + N t N b ∂ 2 θ ∂ y 2 (27)

Eliminating the pressure from Equations (24) and (25) yields

∂ 2 ∂ y 2 [ ( 1 + n − 1 2 W e 2 ( ∂ 2 ψ ∂ y 2 ) 2 ) ∂ 2 ψ ∂ y 2 − ( M 2 1 + m 2 + 1 s ) ψ ] + G r ∂ θ ∂ y + B r ∂ Ω ∂ y = 0 (28)

The non-dimensional boundaries will take the form

h 1 = 1 + a cos ( 2 π x ) , h 2 = − d − b cos ( 2 π x + ϕ ) (29)

where a , b , ϕ , d satisfies the relation

a 2 + b 2 + 2 a b cos ϕ ≤ ( 1 + d ) 2 (30)

The corresponding boundary conditions are

ψ = F 2 , ∂ ψ ∂ y = − 1 , θ = 0 , Ω = 0 , at y = h 1 ψ = − F 2 , ∂ ψ ∂ y = − 1 , θ = 1 , Ω = 1 , at y = h 2 (31)

where F is the dimensionless average flux in the wave frame defined by

F = ∫ h 2 ( x ) h 1 ( x ) ∂ ψ ∂ y d y (32)

The time mean Q in the wave frame is defined by

Q = F + 1 + d (33)

In this work, a program was designed by Mathematica Software (version 10) simulate the application of parametric ND Solve package to find the numerical behaviors of our dimensionless system. Graphical results to the velocity u, the temperature θ , the nano particle concentration Ω and the pressure gradient

d p d x are obtained under the impact of emerging parameters at moving boundaries

of the fluid.

The peristaltic flow of a Carreau nanofluid through a porous medium with heat and mass transfer in the presence of Hall current, Joule heating and viscous dissipation is studied under assumption of long wavelength and low Reynolds number. The main results are summarized as follow:

・ The velocity of Carraeu nanofluid increases at the neighborhood of the channel walls and decreases near the center of the channel by increasing the Weissenber number We and the power law index n.

・ It is observed that the magnetic parameter M and the thermophoresis parameter n_t have opposite effects on the velocity to that of we and n.

・ The temperature profile increases with an increase in the Brownian motion parameter Nb and decreases with an increase in the thermophoresis parameter n_t.

・ The concetration profile decreases with an increase in the Brownian motion parameter Nb and increases with an increase in the thermophoresis parameter n_t.

・ The pressure gradient decreases with the increase of the magnetic parameter M, the power law index n, the Weissenber number we and the local nanoparticle Grashof number br, while it increases with the increase of the local temperature Grashof number g_r and the thermophoresis parameter n_t.

Eldabe, N.T.M., Abo-Seida, O.M., Abo-Seliem, A.A.S., El- Shakhipy, A.A. and Hegazy, N. (2017) Peristaltic Transport of Magnetohydrodynamic Carreau Nanofluid with Heat and Mass Transfer inside Asymmetric Channel. American Journal of Computational Mathematics, 7, 1-20. https://doi.org/10.4236/acjm.2017.71001