Next Article in Journal
Preparations of Silver/Montmorillonite Biocomposite Multilayers and Their Antifungal Activity
Next Article in Special Issue
Impact of Second-Order Slip and Double Stratification Coatings on 3D MHD Williamson Nanofluid Flow with Cattaneo–Christov Heat Flux
Previous Article in Journal
Synthesis and Characterization of Ni–W/Cr2O3 Nanocomposite Coatings Using Electrochemical Deposition Technique
Previous Article in Special Issue
Numerical Investigation of Multiple Solutions for Caputo Fractional-Order-Two Dimensional Magnetohydrodynamic Unsteady Flow of Generalized Viscous Fluid over a Shrinking Sheet Using the Adams-Type Predictor-Corrector Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Convective Heat Transfer and Magnetohydrodynamics across a Peristaltic Channel Coated with Nonlinear Nanofluid

1
Department of Mathematics, Division of Science and Technology, University of Education, Lahore 54770, Pakistan
2
Department of Mathematics, King Saud University, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Coatings 2019, 9(12), 816; https://doi.org/10.3390/coatings9120816
Submission received: 28 October 2019 / Revised: 24 November 2019 / Accepted: 28 November 2019 / Published: 2 December 2019
(This article belongs to the Special Issue Recent Trends in Coatings and Thin Film–Modeling and Application)

Abstract

:
The aim of the current study is to present an analytical and numerical treatment of a two-dimensional peristaltic channel along with the coating of laminar layers of nanoparticles with non-Newtonian (Williamson) base liquid. In addition to this, convective heat transfer and magnetic field effects also take into consideration. The geometry is considered as an asymmetric two dimensional channel experiencing sinusoidal waves propagating across the walls. The walls are supposed to have heat convection at the upper wall and the lower wall is having no temperature gradient. The problem is manufactured under the theory of lubrication approach. The mathematical models are evolved by using appropriate transformations. The obtained nonlinear differential equations are solved analytically. Graphical features are presented to find the influence of emerging physical parameters on the stream function, velocity of the nanofluid, heat transfer, nanoparticles concentration, pressure gradient, and pressure increase. It is found that the velocity decreases in the lower part while increasing in the upper side of the channel in the presence of nanoparticles. The temperature is becoming large with increasing amount of nanoparticles and heat convection at the boundaries. It is also observed that nanoparticle concentration is getting higher with Brownian motion parameter, but fluid becomes less thermal against thermophoresis parameter. The streamlines phenomenon clearly reflects the asymmetry of the channel. The characteristics of viscous fluid can be recovered by switching the Weissenbureg number (We) to zero.

1. Introduction

Nanofluids attract current predilection because of its heat conduction attributes. Changing the flow geometry, boundary conditions, or thermal conductivity of liquids can improve convective heat transfer. Over the years, researchers have tried to increase the thermal conductivity of liquids. For this purpose, with the idea of Maxwell [1] solid metal particles are introduced into the base liquid. The large micro-sized particles are used to make suspensions because the conductivity of solids is greater than that of liquids, but these particles tend to produce greater resistance to the flow of base fluid. Modern nanotechnology tends to take a new direction in this field. In 1995, Choi [2] proposed a liquid with nano-sized particles suspended in a base liquid to eliminate the disadvantages of micro-sized particles. These liquids have efficient convective heat transfer compared to pure liquids. Recently, the idea of nanofluid in peristalsis has been studied by some researchers [3,4,5,6,7,8,9,10].
Peristalsis is characterized as the extension and the arrival of a substance into a liquid that improves the formative waves that broaden the length of the conduit, blending and shipping the liquid toward the wave spread. It is a mechanism that is available in numerous organs of the human body. In some specific instruments—for example heart–lung machines, implantation gadgets, and other pumping apparatus—such types of processes are utilized. It is of specific significance in many species and especially in human body that the transportation of many tissues of the body under various conditions, for example, the sucking of blood by leaches, the heap from the kidneys to the bladder through filtration, transport of the spermatozoa to the male genital tract, the development of the bosom in the Fallopian tubes, vasomotion of little veins, just as the blending and transport of gastrointestinal entry material.
The use of heat is of particular importance in the field due to its wide scope in engineering and biomechanics. In addition, the common relationship of heat stress and peristalsis can be observed during the oxygenation process with the patient. The assessment of magnetic resonance in biological tissues has aroused great interest among researchers regarding physical problems such as blood.
The assessment of heat transfer is related to the conditions of convection used in processes such as thermal conductivity, mechanical properties, chemical reactions, and so on. Aziz [11] presented a similarity solution to incorporate the convective walls conditions for thermal boundary layer on a smooth plate. In another article, Makinde and Aziz [12] developed the MHD mixed model on a flat surface in a concise way in terms of compatibility. Makinde [13] also discussed the flow of the MHD component with the temperature and the mechanical evaluation of a plate on a flat surface with extended conditions. Merkin and Pop [14] considered the analysis of heat transfer by dynamically simulating the flow of a uniform current on a flat surface with a horizontal displacement. According to them, the heat flux near the main edge is dominated by the surface heat flux.
After knowing the significance of the above discussed phenomena, authors are keen to develop a series solution of peristaltic flow of nanofluid with Williamson fluid model as a base liquid with convective boundary conditions travelling through asymmetric channel. At least we know that this study has not been yet explored in the literature. This study will be a good base for the engineers to utilize the results in procedures like thermal energy storage, gas turbines, nuclear workshops, etc. The problem is modeled under the induction of lubrication approach. The series solutions of stream function, temperature distribution, and nanoparticle concentration are achieved by using a well-known converging method the homotopy perturbation method. The important features are analyzed more specifically by sketching graphs to estimate the impact of pertinent constant physical factors.

2. Mathematical Modeling

The incompressible Williamson model is chosen as a base fluid for nanofluid in between an asymmetric channel experiencing heat convection at the peristaltic type surfaces. The width of the channel is taken as ( d 11 + d 12 ) . Flow is initiated due to the propagation of curved waves travelling with uniform speed c towards the flow. The exchange of heat is recognized by imposing temperatures T 0 and T 1 at the lower and upper areas, correspondingly. To discuss nano particle phenomenon, we have taken the nanoparticle concentration C 0 and on the lower side and upper one, accordingly (see Figure 1).
The magnetic field B 0 is exerted orthogonally. The wall surfaces are taken as
Y = H 1 = d 11 + a 11 cos [ 2 π λ X ]
Y = H 2 = d 12 b 11 cos [ 2 π λ X + φ ¯ ] ,   where   X = X c t
In upper defined equations, a 11 and b11 represent the wave amplitudes, λ gives the wavelength, t suggests the time, X depicts the wave’s direction, and Y is placed normally to X . The range of phase variance φ ¯ alters as 0 φ ¯ π . If φ ¯ = 0 , we meant that a symmetric dimensional channel is having waves located out of the phase and φ ¯ = π , suggest the waves within the phase. Moreover a 11 ,   b 11 ,   d 11 ,   d 12 and φ ¯ overcome the following relation
a 11 2 + b 11 2 + 2 a 11 b 11 cos φ ¯ ( d 11 + d 12 ) 2
The mathematical models of the considered problem given as
V = 0
ρ ( V t + V V ) = P ¯ + S + ρ f g α f ( T ˜ T 0 ) + ρ f g α f ( C ˜ C 0 ) + J × B
( ρ c ) f ( T ˜ t + V T ˜ ) = K T ˜ + S V + ( ρ c ) p ( D B ( C ˜ T ˜ ) + D T T m ( T ˜ T ˜ ) )
C ˜ t + V ^ C ˜ = ( D B C ˜ + D T T ˜ T o )
where g is the gravitational body force and α f represents the volumetric volume distension nanofluid’s coefficient. In above relations, ( ρ c ) f denotes the fluid’s heat capacity, ( ρ c ) p accounts for effective nanoparticles heat capacity, J = σ ( V × B ¯ ) reveals the current density, B ¯ = ( 0 , B 0 ) notifies the external magnetic field, and S is placed for the Cauchy stress tensor for Williamson fluid and is determined as
τ = ( μ + ( μ 0 + μ )   ( 1 Γ γ ) 1 )   γ ¯
where γ ¯ comprises the subsequent value
γ ¯ = 1 2 i j γ i j γ j i = 1 2 Π ¯
Here Π ¯ is the strain tensor. The velocity profile for the given problem is considered as V = ( U ,   V ) . Introducing a wavy frame we introduce the following transformations
x = X ,   y = Y ,   u = U c ,   v = V ,   p ( x ) = P ¯ ( X ,   t )
We suggest the following dimensionless parameters to be used in the above expressions
x ¯ = 2 π x λ ,   y ¯ = y d 11 ,   u ¯ = u c ,   v ¯ = v c δ ,   δ = d 11 λ ,   d = d 12 d 11 ,   p ¯ = d 11 2 p μ c λ ,   h 11 = H 11 d 11 ,   h 12 = H 12 d 12 , a 12 = a 11 d 11 ,   B r = E c Pr ,   b = b 11 d 11 ,   Re = ρ c d 11 μ ,   ψ ¯ = ψ c d 11 ,   θ = T ˜ T 0 T 1 T 0 ,   E c = c 2 c p ( T 1 T 0 ) , Pr = ρ ν c K ,   S ¯ = S d 11 μ c ,   W e = Γ c d 11 ,   φ = C ˜ C 0 C 1 C 0 ,   G r = ρ f g α f d 2 ( T 1 T 0 ) c μ ,   G c = ρ f g α f d 2 ( C 1 C 0 ) c μ , N b = τ D B ( C 1 C 0 ) ν ,   N t = τ D T ( T 1 T 0 ) τ T m , M = σ μ B 0 d 11 , W e = Γ c d 11
where M , W e , B r ,   Pr ,   N b ,   N t ,   G r , and G c represent the Hartman number, Weissenberg number, Brinkman number, Prandtl number, Brownian motion parameter, thermophoresis parameter, local temperature Grashof number, and local nanoparticle Grashof number, accordingly. After incorporating the above structured parameters and applying the conditions of large wavelength along with small Reynolds number in a wavy frame coordinates we have the final form of Equations (4)–(7)
u x + v y = 0 d p d x = y [ 2 ψ y 2 M 2 ψ + W e ( 2 ψ y 2 ) 2 ] + G r θ + G c φ  
or
2 y 2 [ 2 ψ y 2 + W e ( 2 ψ y 2 ) 2 M 2 ψ ] + G r θ y + G c φ y = 0
Pr [ N b φ y θ y + N t ( θ y ) 2 ] + 2 θ y 2 + B r [ ( 2 ψ y 2 ) 2 + W e ( 2 ψ y 2 ) 3 ] = 0
2 φ y 2 + N t N b 2 θ y 2 = 0
where ψ is stream function satisfying the relations u = ∂δ/∂y and v = δ ψ / x , The no-slip boundary conditions for velocity u and nanoparticles fraction φ and convective boundaries are taken into consideration for temperature θ which have the following dimensionless form in the wave frame [15]
ψ = F 2 ,     ψ y = 1 ,     at   y = h 11 ,     ψ = F 2 ,     ψ y = 1 ,     at   y = h 12 , θ y B i θ = B i   at   y = h 11   and   θ = 0   at   y = h 12 ) φ = 1               at   y = h 11 ,       and       φ = 0   at   y = h 12 , }
where h 11 = 1 + a 12 cos x and h 12 = d b cos ( x + φ ) . Also B i = h f d 11 / K is the Biot number, h f stands for the coefficient of convective thermal transport. The mean flow rate in dimensionless format is elaborated as
Q = F + 1 + d

3. Solution of the Problem

The above obtained Equations (12)–(15) display the nonlinear ordinary differential equations in which ψ , θ , and φ are mutually dependent. Such types of problems cannot be handled by exact techniques. Therefore, we chose a more appropriate solution procedure, the homotopy perturbation method (HPM) [16,17] to solve the current highly complicated boundary value problems. The deformation equations for ψ , θ , and φ can be constructed as
( 1 q ) 1 ( ψ ^ ψ 0 ) + q [ 2 y 2 ( 2 ψ ^ y 2 + W e ( 2 ψ ^ y 2 ) 2 M 2 ψ ^ ) + G r θ ^ y + G c φ ^ y ] =   0 ,
( 1 q ) 2 ( θ ^ θ 0 ) + q [ N b θ ^ y φ ^ y + N t ( θ ^ y ) 2 + 2 θ ^ y 2 + B r ( 2 ψ ^ y 2 ) 2 + W e ( 2 ψ ^ y 2 ) 3 ] =   0 ,
( 1 q ) 2 ( φ ^ φ 0 ) + q [ 2 φ ^ y 2 + N t N b 2 θ ^ y 2 ] =   0 ,
where £ 1 and £ 2 are linear operators which are picked as £
£ 1 = 4 y 4   and   £ 2 = 2 y 2
and ψ 0 , θ 0 , and φ 0 are the initial approximations which must satisfy the boundary conditions as well as differential operator. The initial approximations for ψ ,   θ ,   and   φ are elected as
ψ 0 = ( h 11 h 12 2 y )   ( 2 ( h 11 h 12 )   ( h 11 y )   ( h 12 y ) ) 2 ( h 11 h 12 ) 2 + F ( h 11 2 4 h 11 h 12 + h 12 2 + 2 ( h 11 + h 12 )   y 2 y 2 2 ( h 11 h 12 ) 2 )
θ 0 = B i h 12 B i y 1 B i h 11 + B i h 12
φ 0 = h 12 + y h 11 h 12
Applying perturbation on small embedding parameters F [ 0 ,   1 ] , we suggest the following series solutions
ψ = ψ 0 + q ψ 1 + q 2 ψ 2  
θ = θ 0 + q θ 1 + q 2 θ 2
φ = φ 0 + q φ 1 + q 2 φ 2  
After substituting the above series solutions in Equations (18)–(20), we get the two systems for ψ , θ , and φ .
  • Zeroth Order System
    £ 1 [ ψ 0 ψ 0 ] = 0 , ψ 0 = F 2 ,   ψ 0 y = 1 ,   at   y = h 1 ,   ψ 0 = F 2 ,   ψ 0 y = 1 ,   at   y = h 2 , }
    £ 2 [ θ 0 θ 0 ] = 0 , θ 0 ( h 1 ) B i θ 0 ( h 1 ) = B i   at   y = h i   and   θ 0 = 0   at   y = h 2 }
    £ 2 [ φ 0 φ 0 ] = 0 , φ 0 = 1 ,     at   y = h 1   and   φ 0 = 0   at   y = h 2 , }
  • First Order System
    £ 1 [ ψ ] + 2 y 2 [ 2 ψ 0 y 2 + W e ( 2 ψ 0 y 2 ) 2 M 2 ψ 0 ] + G r θ 0 y + G c φ 0 y = 0 , ψ 1 = 0 ,   ψ 1 y = 0 ,     at   y = h 1   and   ψ 1 = 0 ,   ψ 1 y = 0 ,     at   y = h 2 , }
    £ 2 ( θ 1 ) + N b [ φ 0 y θ 0 y ] + N t ( θ 0 y ) 2 + 2 θ 0 y 2 + B r [ ( 2 ψ 0 y 2 ) 2 + W e ( 2 ψ 0 y 2 ) 3 ] = 0 , θ 1 ( h 1 ) B i θ 1 ( h 1 ) = 0   at   y = h 1   and   θ 1 = 0   at   y = h 2 }
    £ 2 ( φ 1 ) + 2 φ 0 y 2 + N t N b 2 θ 0 y 2 = 0 , φ 1 = 0 ,     a t   y = h 1 ,   φ 1 = 0   a t   y = h 2 , }
  • Zeroth Order Solutions
    By solving zeroth order systems by built-in technique in mathematical software, we obtain
    ψ 0 = ψ 0 = ( h 11 + h 12 2 y )   ( 2 ( h 11 h 12 )   ( h 11 y )   h 12 y ) 2 ( h 11 h 12 ) + F ( h 11 2 4 h 11 h 12 + h 12 2 + 2 ( h 11 + h 12 )   y ) 2 y 2 ( h 11 h 12 ) 2  
    θ 0 = θ 0 = B i h 2 B i y 1 B i h 1 + B i h 2
    φ 0 = φ 0 = h 12 + y h 11 h 12
  • First Order Solutions
    The first order system has acquired the following general solutions
    ψ 1 = 1 6 ( 1 + B i ( h 11 h 12 ) 7 ) [ 1 / 4 G c ( 1 + B i ( h 11 h 12 ) 6 6 ( F + h 11 h 12 ) ( h 11 M 2 + 2 h 11 3 h 12 M 2 2 h 11 h 12 3 M 2 + h 12 4 M 2 + 48 W e ( F + h 11 h 12 ) + B i ( h 11 h 12 ) ( G r ( h 11 h 12 ) 5 + 6 ( F + h 11 h 12 )   ( h 11 4 M 2 + 2 h 11 3 h 12 M 2 + + h 12 4 M 2 + 48 W e ( F + h 11 h 12 ) ) ) ) y 4 + 3 5 ( 1 + B i ( h 11 h 12 ) 4 ( F + h 11 h 12 )   M 2 y 5 ) ] + L 11 + y L 12 + y 2 L 13 + y 3 L 14
    θ 1 = 1 ( 1 + B i ( h 11 h 12 ) ) 2 ( h 11 h 12 ) 9 [ 1 / 2 B i ( h 11 h 12 ) 8 ( ( 1 + B i ( h 11 h 12 ) )   N b + B i ( h 11 h 12 )   N t ) Pr + 36 B r ( 1 + B i ( h 11 h 12 ) ) 2 ( F + h 11 h 12 ) 2 ( h 11 h 12 ) 2 ( h 11 3 3 h 12 2 ( h 12 2 W e ) + 3 h 11 ( h 12 2 + 2 F W e ) h 12 ( h 12 2 6 W e + 6 h 12 W e ) ) ) y 2 24 B r ( 1 + B i ( h 11 h 12 ) ) 2 ( F + h 11 h 12 ) 2 ( h 11 + h 12 )   h 11 3 3 h 11 2 ( h 12 3 W e ) + 3 h 12 ( h 12 2 + 3 F W e ) ( h 12 ( h 12 2 9 F W e + 9 h 12 W e ) ) y 3 + 12 B r ( 1 + B i ( h 11 h 12 ) ) 2 ( F + h 11 h 12 ) 2 ( h 11 3 3 h 11 2 ( h 12 6 W e ) + 3 h 11 ( h 12 2 + 6 F W e ) h 12 ( h 12 2 18 F W e + 18 h 12 W e ) ) y 4 ) ] 432 5 B r ( 1 + B i ( h 11 h 12 )   ( F + h 11 h 12 ) 2 ) 3 W e y 5 + L 15 + y L 16  
    φ 1 = 1 ( 1 + B i ( h 11 h 12 ) ) 2 ( h 11 h 12 ) 9 N b N t [ 1 / 3 B i ( h 11 h 12 ) 8 ( ( 1 + B i ( h 11 h 12 ) )   N b + B i ( h 11 h 12 )   N t ) Pr + 36 B r ( 1 + B i ( h 11 h 12 ) ) 2 ( F + h 11 h 12 ) 2 ( h 11 h 12 ) 2 ( h 11 3 3 h 11 2 ( h 12 2 W e + 3 h 1 ( h 12 2 + 2 F W e ) h 12 ( h 12 2 6 F W e + 6 h 12 W e ) ) ) Y 2 24 B r ( 1 + B i ( h 11 h 12 ) ) 2 ( F + h 11 h 12 ) 2 ( h 11 + h 12 )   h 11 3 3 h 12 2 ( h 12 3 W e ) + h 11 ( h 12 2 + 3 F W e )   ( h 12 ( h 12 2 9 F W e ) + 9 h 2 W e )   Y 3 + 12 G c ( 1 + B i ( h 11 h 12 ) 2 ( F + h 11 h 12 ) 2 ( h 11 3 3 h 11 2 ( h 12 6 W e ) + 3 h 11 2 ( h 12 6 W e ) + 3 h 11 ( h 11 2 + 6 F W e ) h 12 ( h 12 2 18 F W e ) )   y 4 ] ( 432 5 B r ( 1 + B i ( h 11 h 12 ) ) 2 ( F + h 11 h 12 ) 2 ) 3 W e y 5 + L 17 + y L 18
The final solutions according to the concept of HPM are given by using q → 1 in Equations (2)–(27).
ψ = ψ 0 + ψ 1 +
θ = θ 0 + θ 1 +
φ = φ 0 + φ 1 +
where constants L i j , i = 1 ,   j = 1 8 can be found by routine calculation. The complete solutions of ψ , θ , and φ can be obtained by supposed solutions. The solution for pressure gradient d p / d x can be found by simply substituting the values in Equation (12). The mathematical formula for the pressure increase function p can been visualized in next equation that has been solved numerically by built-in technique numerical integration on Mathematica.
p = 0 1 ( d p d x )   d x

4. Results and Discussion

This portion comprises of graphical results and discussion of obtained results for velocity, temperature, nanoparticles, pressure gradient, and stream functions. The numerical data of the pressure rise function p is also sketched against the domain of flow rate and found the effects of physical parameters separately. Figure 2 and Figure 3 are sketched for the velocity profile with varying the values of ( G r ) and ( G c ) , respectively in corresponding order. From Figure 2, it is clearly visible that velocity is decreasing in lower part and increasing in upper part of the channel and enhances its maximum peak at the center under the effect of G r . One can see the similar behavior by taking increasing values of G c but here the difference is that the velocity is not varying much under the effect of G c in Figure 3. Figure 4, Figure 5 and Figure 6 contain correspondingly the alteration of temperature profile θ with the variability of Biot number ( B i ) , Brinkman number ( B r ) , and the Prandtl number ( Pr ) . From Figure 4, one can notice that the temperature profile is stretched vertically with the increase in magnitudes of B i . It depicts that heat convection at the boundaries enhances the temperature of the Williamson nanofluid. It is also notable here that the temperature is maximum at lower wall and minimum at the lower surface and there is much variation in temperature level at upper region as compared to lower side. Figure 5 reflects the observation that temperature is an increasing function of B r and the temperature gradients are prominent at the lower portions as equated with the upper ones, but the extent of heat is similar at both the surfaces as was observed for B i . It can be received from Figure 6 that temperature profile is increasing in linear fashion for numerically increasing magnitudes of P r but the change in heat is calculated more significantly in the central parts of the enclosure which is the totally different result than we have achieved in Figure 4 and Figure 5. Figure 7 and Figure 8 are presented to see the behavior of nanoparticles volume fraction φ with increasing magnitudes of ( N b ) and ( N t ) . Figure 7 shows that φ is getting higher when someone increases N b . It is also explicit here that nanoparticles are dispersed in the region between the lower and upper surfaces. On the other hand, Figure 8 revels different story, the increase in N t decreases the nanoparticles concentration. Figure 9 is plotted for pressure gradient d p / d x for N b . It is seen that d p / d x is increasing as we increase N b and gets maximum height at the center of the domain, i.e., x = 0.5 . From Figure 10, we can see that pressure gradient is varying quite opposite manner for the parameter N t . It can also be noticed from Figure 9 and Figure 10 that pressure gradient gets positive values only in the central part and remains negative at the corners. Figure 11 and Figure 12 are displaced to see the effects of parameters M and W e on pressure rise p . Here the whole area is broken into three zones, namely Region I–III. The Region I is recognized by the portion where Q > 0 ,   p > 0 . Region II is named the place where Q > 0 and p < 0 while Region III is composed of the part Q ,   p < 0 . Figure 11 shows that p curves are increasing in Region I and II while decreasing in Region III with the variation of M . Also, the free pumping exists at Q 1.5 . In Figure 12, it is observed that in Region I and II, p is increasing and in Region III, it is decreasing. Also, the peristaltic pumping occurs in Regions I and II between the interval ( 1.7 ,   0.5 ) . The streamlines are drawn in Figure 13, Figure 14 and Figure 15 for the parameters G c , W e , and M , respectively. From Figure 13, it is clear that the number of boluses is increasing, but size of the trapped bolus is decreasing in lower part of the channel, while in upper portion, the situation is totally reflected in opposite ways. Figure 14 gives the streamlines variation under the different values of W e . It is attained here that, in the lower part, the number of boluses is increasing but size is changing randomly. The stream function for M has been sketched in Figure 15 and it is noted in both the lower and upper parts, the size of bolus in increasing while number is decreasing. It is also admitted by Figure 13, Figure 14 and Figure 15 that trapped boluses are displaced towards left from upper to lower side due to asymmetric dimensions of the channel which can be made symmetric by imposing φ ¯ = 0 .

5. Conclusions

In this article, the authors have discovered the mathematical treatment of the peristaltic flow of Williamson nanofluid coated with the walls of an asymmetric heated channel. The flow has been studied analytically and graphically through variation of some pertinent parameters. From the above discussion, the main findings are given below:
(1)
The velocity of nanofluid is decreasing in the lower part while increasing in the upper side with local temperature Grashof number and local nanoparticle Grashof number.
(2)
The temperature is becoming large with an increase in Biot number, Brinkman number, and Prandtl number.
(3)
The nano particle concentration is getting higher when we increase Brownian motion parameter, but diminishes with thermophoresis parameter.
(4)
The pressure gradient is increasing with Brownian motion parameter, but lessening for thermophoresis parameter.
(5)
The peristaltic pumping fasten up with Hartman number and Weissenberg number.
(6)
In the upper portion, the size of the trapped bolus is decreasing, but increasing in lower portion when we increase local nanoparticle Grashof numbers and Weissenberg numbers, but it varies in a random manner with Hartman numbers.
(7)
It is important to notice that boluses are trapped by their position in lower and upper corners of the channel due to its asymmetric structure. We can recover the results of symmetric channel by neglecting the phase difference.
(8)
The study of viscous nanofluid can be approached by neglecting Weissenburg number.

Author Contributions

Methodology, formal Analysis and writing—original draft preparation, A.R. (Arshad Riaz); writing—review and editing, A.R. (Abdul Razaq); funding acquisition, H.A.

Funding

This research project was supported by a grant from the Research Center of the Center for Female Scientific and Medical Colleges, Deanship of Scientific Research, King Saud University, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maxwell, J.C. A Treatise on Electricity and Magnetism, 2nd ed.; Clarendan Press: Oxford, UK, 1881. [Google Scholar]
  2. Choi, S.U.S.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Argonne National Lab: Du Page County, IL, USA, 1995. [Google Scholar]
  3. Hassan, M.; Marin, M.; Alsharif, A.; Ellahi, R. Convective heat transfer flow of nanofluid in a porous medium over wavy surface. Phys. Lett. A 2018, 382, 2749–2753. [Google Scholar] [CrossRef]
  4. Safaei, M.R.; Togun, H.; Vafai, K.; Kazi, S.N.; Badarudin, A. Investigation of heat transfer enhancement in a forward-facing contracting channel using FMWCNT nanofluids. Numer. Heat Transf. Part A Appl. 2014, 66, 1321–1340. [Google Scholar] [CrossRef]
  5. Alsagri, A.S.; Nasir, S.; Gul, T.; Islam, S.; Nisar, K.S.; Shah, Z.; Khan, I. MHD Thin Film Flow and Thermal Analysis of Blood with CNTs Nanofluid. Coatings 2019, 9, 175. [Google Scholar] [CrossRef]
  6. B’eg, O.A.; Tripathi, D. Mathematica simulation of peristaltic pumping with double-diffusive convection in nanofluids: A bio-nano-engineering model. J. Nanoeng. Nanosyst. 2011, 225, 99–114. [Google Scholar]
  7. Akbar, N.S.; Maraj, E.N.; Butt, A.W. Copper nanoparticles impinging on a curved channel with compliant walls and peristalsis. Eur. Phys. J. Plus. 2014, 129, 183. [Google Scholar] [CrossRef]
  8. Bhatti, M.M.; Rashidi, M.M. Effects of thermo-diffusion and thermal radiation on Williamson nanofluid over a porous shrinking/stretching sheet. J. Mol. Liq. 2016, 221, 567–573. [Google Scholar] [CrossRef]
  9. Bhatti, M.M.; Sheikholeslami, M.; Zeeshan, A. Entropy analysis on electro kinetically modulated peristaltic propulsion of magnetized nanofluid flow through a microchannel. Entropy 2017, 19, 481. [Google Scholar] [CrossRef]
  10. Prakash, J.; Tripathi, D.; Tiwari, A.K.; Sait, S.M.; Ellahi, R. Peristaltic Pumping of Nanofluids through a Tapered Channel in a Porous Environment: Applications in Blood Flow. Symmetry 2019, 11, 868. [Google Scholar] [CrossRef]
  11. Aziz, A. A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 1064–1068. [Google Scholar] [CrossRef]
  12. Makinde, O.D.; Aziz, A. MHD mixed convection from a vertical plate embedded in a porous medium with a convective boundary condition. Int. J. Therm. Sci. 2010, 49, 1813–1820. [Google Scholar] [CrossRef]
  13. Makinde, O.D. Similarity solution of hydromagnetic heat and mass transfer over a vertical plate with a convective surface boundary condition. Int. J. Phys. Sci. 2010, 5, 700–710. [Google Scholar]
  14. Merkin, J.H.; Pop, I. The forced convection flow of a uniform stream over a flat surface with a convective surface boundary condition. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 602–3609. [Google Scholar] [CrossRef]
  15. Akbar, N.S. Non-Newtonian fluid flow in an asymmetric channel with convective surface boundary condition: A note. J. Power Technol. 2014, 94, 34–41. [Google Scholar]
  16. He, J.H. Homo topy perturbation method for solving boundary value problems. Phys. Lett. A 2006, 350, 87–88. [Google Scholar] [CrossRef]
  17. He, J.H. Homotopy perturbation method: A new nonlinear analytical technique. Appl. Math. Comput. 2003, 135, 73–79. [Google Scholar] [CrossRef]
Figure 1. Geometry of the channel.
Figure 1. Geometry of the channel.
Coatings 09 00816 g001
Figure 2. Modification of velocity profile against G r for x = 1 ,   F = 2 ,   a = 0.2 ,   b = 0.1 ,   d = 1.5 ,   φ = 1.5 ,   G c = 0.3 ,   W e = 0.01 ,   M = 0.1 ,   B i = 0.5 .
Figure 2. Modification of velocity profile against G r for x = 1 ,   F = 2 ,   a = 0.2 ,   b = 0.1 ,   d = 1.5 ,   φ = 1.5 ,   G c = 0.3 ,   W e = 0.01 ,   M = 0.1 ,   B i = 0.5 .
Coatings 09 00816 g002
Figure 3. Modification of velocity profile against G c for x = 1 ,   F = 2 ,   a = 0.2 ,   b = 0.1 ,   d = 1.5 ,   φ = 1.5 ,   G r = 0.1 ,   W e = 0.01 ,   M = 0.1 ,   B i = 0.5 .
Figure 3. Modification of velocity profile against G c for x = 1 ,   F = 2 ,   a = 0.2 ,   b = 0.1 ,   d = 1.5 ,   φ = 1.5 ,   G r = 0.1 ,   W e = 0.01 ,   M = 0.1 ,   B i = 0.5 .
Coatings 09 00816 g003
Figure 4. Modification of temperature profile against B i for, n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W = 0.01 ,   M = 0.1 ,   G c = 0.9 ,   G r = 4 ,   G c = 0.3 ,   N b = 0.5 ,   Pr = 0.4 ,   N t = 0.2 .
Figure 4. Modification of temperature profile against B i for, n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W = 0.01 ,   M = 0.1 ,   G c = 0.9 ,   G r = 4 ,   G c = 0.3 ,   N b = 0.5 ,   Pr = 0.4 ,   N t = 0.2 .
Coatings 09 00816 g004
Figure 5. Modification of temperature profile against B r for n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.1 ,   φ = 0.1 ,   W e = 0.01 ,   G c = 0.9 ,   G r = 1 ,   B i     = 10 ,   N b = 0.5 ,   Pr = 0.4 ,   N t = 0.2 .
Figure 5. Modification of temperature profile against B r for n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.1 ,   φ = 0.1 ,   W e = 0.01 ,   G c = 0.9 ,   G r = 1 ,   B i     = 10 ,   N b = 0.5 ,   Pr = 0.4 ,   N t = 0.2 .
Coatings 09 00816 g005
Figure 6. Modification of temperature profile against Pr for n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.1 ,   W e = 0.01 ,   G c = 0.9 ,   G r = 1 ,   B i     = 5 ,   N b = 0.5 ,   G c = 0.01 ,   N t = 0.2 .
Figure 6. Modification of temperature profile against Pr for n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.1 ,   W e = 0.01 ,   G c = 0.9 ,   G r = 1 ,   B i     = 5 ,   N b = 0.5 ,   G c = 0.01 ,   N t = 0.2 .
Coatings 09 00816 g006
Figure 7. Modification of nanoparticles concentration against N b for n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W e = 0.01 ,   G c = 0.9 ,   G r = 4 ,   B i     = 0.5 ,   Pr = 0.4 ,   G c = 0.1 ,   N t = 0.2 .
Figure 7. Modification of nanoparticles concentration against N b for n = 2 ,   x = 0.1 ,   F = 5 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W e = 0.01 ,   G c = 0.9 ,   G r = 4 ,   B i     = 0.5 ,   Pr = 0.4 ,   G c = 0.1 ,   N t = 0.2 .
Coatings 09 00816 g007
Figure 8. Modification of nanoparticles concentration against N t for, n = 2 ,   x = 0.1 ,   F = 2 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W e = 0.01 , , G c = 0.9 ,   G r = 4 ,   B i = 0.5 ,   Pr = 0.4 ,   G c = 0.1 ,   N t = 0.2 .
Figure 8. Modification of nanoparticles concentration against N t for, n = 2 ,   x = 0.1 ,   F = 2 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W e = 0.01 , , G c = 0.9 ,   G r = 4 ,   B i = 0.5 ,   Pr = 0.4 ,   G c = 0.1 ,   N t = 0.2 .
Coatings 09 00816 g008
Figure 9. Modification of pressure gradient against N b for n = 2 ,   y = 0.1 ,   F = 10 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 , W e = 0.1 ,   G c = 0.9 ,   G r = 4 ,   B i     = 0.5 ,   Pr = 0.4 ,   G c = 0.1 ,   M = 1.5 ,   N t = 0.2 .
Figure 9. Modification of pressure gradient against N b for n = 2 ,   y = 0.1 ,   F = 10 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 , W e = 0.1 ,   G c = 0.9 ,   G r = 4 ,   B i     = 0.5 ,   Pr = 0.4 ,   G c = 0.1 ,   M = 1.5 ,   N t = 0.2 .
Coatings 09 00816 g009
Figure 10. Modification of pressure gradient against N t for, n = 2 ,   y = 0.1 ,   F = 10 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W e = 0.1 , , G c = 0.9 ,   G r = 4 ,   B i     = 0.09 ,   Pr = 0.4 ,   G c = 0.3 ,   M = 1.5 ,   N b = 0.1 .
Figure 10. Modification of pressure gradient against N t for, n = 2 ,   y = 0.1 ,   F = 10 ,   a = 0.2 ,   b = 0.1 ,   d = 0.51 ,   φ = 0.01 ,   W e = 0.1 , , G c = 0.9 ,   G r = 4 ,   B i     = 0.09 ,   Pr = 0.4 ,   G c = 0.3 ,   M = 1.5 ,   N b = 0.1 .
Coatings 09 00816 g010
Figure 11. Modification of pressure rise against M for, n = 2 ,   y = 0.1 ,   F = 10 ,   a = 0.2 ,   b = 0.3 ,   d = 0.5 ,   φ = 0.01 ,   M = 1.3 ,   G c = 0.3 ,   G r = 0.1 ,   B i     = 0.3 ,   Pr = 0.4 ,   G c = 0.3 ,   N b = 0.3 ,   N t = 0 .
Figure 11. Modification of pressure rise against M for, n = 2 ,   y = 0.1 ,   F = 10 ,   a = 0.2 ,   b = 0.3 ,   d = 0.5 ,   φ = 0.01 ,   M = 1.3 ,   G c = 0.3 ,   G r = 0.1 ,   B i     = 0.3 ,   Pr = 0.4 ,   G c = 0.3 ,   N b = 0.3 ,   N t = 0 .
Coatings 09 00816 g011
Figure 12. Modification of pressure rise against W e for, n = 0.1 ,   y = 0.1 ,   a = 0.2 ,   b = 0.3 ,   d = 0.5 ,   φ = 0.01 ,   M = 1.3 ,   G c = 0.3 ,   G r = 0.1 ,   B i     = 0.3 ,   Pr = 0.4 ,   G c = 0.3 ,   N b = 0.3 ,   N t = 0.2 .
Figure 12. Modification of pressure rise against W e for, n = 0.1 ,   y = 0.1 ,   a = 0.2 ,   b = 0.3 ,   d = 0.5 ,   φ = 0.01 ,   M = 1.3 ,   G c = 0.3 ,   G r = 0.1 ,   B i     = 0.3 ,   Pr = 0.4 ,   G c = 0.3 ,   N b = 0.3 ,   N t = 0.2 .
Coatings 09 00816 g012
Figure 13. Modification of streamlines for G c = { 0.1 , 0.5 , 0.9 } when n = 0.1 ,   y = 0.1 ,   F = 5 ,   a = 0.3 ,   b = 0.2 ,   d = 0.1 ,   φ = 0.01 ,   M = 0.1 ,   W e = 0.1 ,   G r = 4 ,   B i     = 0.09 ,   Pr = 0.4 ,   B r = 0.9 ,   N b = 0.5 ,   N t = 0.2 .
Figure 13. Modification of streamlines for G c = { 0.1 , 0.5 , 0.9 } when n = 0.1 ,   y = 0.1 ,   F = 5 ,   a = 0.3 ,   b = 0.2 ,   d = 0.1 ,   φ = 0.01 ,   M = 0.1 ,   W e = 0.1 ,   G r = 4 ,   B i     = 0.09 ,   Pr = 0.4 ,   B r = 0.9 ,   N b = 0.5 ,   N t = 0.2 .
Coatings 09 00816 g013
Figure 14. Modification of streamlines for W e = { 0.1 , 0.2 , 0.3 } when n = 0.1 ,   y = 0.1 ,   a = 0.2 ,   b = 0.3 ,   d = 0.5 ,   φ = 0.01 ,   M = 0.1 ,   B i     = 0.09 ,   Pr = 0.4 ,   G c = 0.9 ,   G r = 4 ,   N b = 0.5 ,   N t = 0.2 .
Figure 14. Modification of streamlines for W e = { 0.1 , 0.2 , 0.3 } when n = 0.1 ,   y = 0.1 ,   a = 0.2 ,   b = 0.3 ,   d = 0.5 ,   φ = 0.01 ,   M = 0.1 ,   B i     = 0.09 ,   Pr = 0.4 ,   G c = 0.9 ,   G r = 4 ,   N b = 0.5 ,   N t = 0.2 .
Coatings 09 00816 g014
Figure 15. Modification of streamlines for M = { 0.1 , 0.9 , 1.7 } when n = 2 ,   F = 5 ,   a = 0.3 ,   b = 0.2 ,   d = 1 ,   φ = 0.01 ,   W e = 0.1 ,   G r = 4 ,   G c = 0.1 ,   B i     = 0.09 ,   Pr = 0.4 ,   G c = 0.9 ,   N b = 0.5 ,   N t = 0.2 .
Figure 15. Modification of streamlines for M = { 0.1 , 0.9 , 1.7 } when n = 2 ,   F = 5 ,   a = 0.3 ,   b = 0.2 ,   d = 1 ,   φ = 0.01 ,   W e = 0.1 ,   G r = 4 ,   G c = 0.1 ,   B i     = 0.09 ,   Pr = 0.4 ,   G c = 0.9 ,   N b = 0.5 ,   N t = 0.2 .
Coatings 09 00816 g015

Share and Cite

MDPI and ACS Style

Riaz, A.; Alolaiyan, H.; Razaq, A. Convective Heat Transfer and Magnetohydrodynamics across a Peristaltic Channel Coated with Nonlinear Nanofluid. Coatings 2019, 9, 816. https://doi.org/10.3390/coatings9120816

AMA Style

Riaz A, Alolaiyan H, Razaq A. Convective Heat Transfer and Magnetohydrodynamics across a Peristaltic Channel Coated with Nonlinear Nanofluid. Coatings. 2019; 9(12):816. https://doi.org/10.3390/coatings9120816

Chicago/Turabian Style

Riaz, Arshad, Hanan Alolaiyan, and Abdul Razaq. 2019. "Convective Heat Transfer and Magnetohydrodynamics across a Peristaltic Channel Coated with Nonlinear Nanofluid" Coatings 9, no. 12: 816. https://doi.org/10.3390/coatings9120816

APA Style

Riaz, A., Alolaiyan, H., & Razaq, A. (2019). Convective Heat Transfer and Magnetohydrodynamics across a Peristaltic Channel Coated with Nonlinear Nanofluid. Coatings, 9(12), 816. https://doi.org/10.3390/coatings9120816

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