# Heat and Mass Transfer in Unsteady Boundary Layer Flow of Williamson Nanofluids

Authors: Tesfaye Kebede; Eshetu Haile; Gurju Awgichew; Tadesse Walelign
Journal: Journal of Applied Mathematics
 (2020)

---

## Abstract

In this paper, analytic approximation to the heat and mass transfer characteristics of a two-dimensional time-dependent flow of Williamson nanofluids over a permeable stretching sheet embedded in a porous medium has been presented by considering the effects of magnetic field, thermal radiation, and chemical reaction. The governing partial differential equations along with the boundary conditions were reduced to dimensionless forms by using suitable similarity transformation. The resulting system of ordinary differential equations with the corresponding boundary conditions was solved via the homotopy analysis method. The results of the study show that velocity, temperature, and concentration boundary layer thicknesses generally decrease as we move away from the surface of the stretching sheet and the Williamson parameter was found to retard the velocity but it enhances the temperature and concentration profiles near the surface. It was also found that increasing magnetic field strength, thermal radiation, or rate of chemical reaction speeds up the mass transfer but slows down the heat transfer rates in the boundary layer. The results of this study were compared with some previously published works under some restrictions, and they are found in excellent agreement.


---

## Body

## 1. Introduction

The term boundary layer flow refers to a kind of flow in a relatively narrow region near a solid surface where the effect of viscosity is significant. The theory of boundary layer flow was first introduced by the German scientist Prandtl in 1904. On the other hand, the field of magnetohydrodynamics (MHD) which deals with the interaction between the velocity of electrically conducting fluids and the electromagnetic fields was first introduced in 1942 by the Swedish physicist and electrical engineer Alfvén [1]. In 1961, Sakiadis [2] initiated the study of boundary layer flows over stretching surfaces and formulated the two-dimensional boundary layer equations. Following these pioneer works, many investigators have reported various useful study results. For instance, Kumar [3] used the Runge–Kutta fourth-order method along with shooting technique to analyze the effect of linear thermal stratification in stable stationary ambient fluid on steady MHD convective flow of a viscous incompressible electrically conducting fluid along a stretching sheet in the presence of mass transfer and magnetic effect. Also, Narasu et al. [4] used perturbation technique to obtain analytic solution for free convective unsteady fluid flow in the presence of thermal diffusion and chemical reaction past a vertical porous plate with heat source and slip effects.In order to enhance the transfer of heat in thermal systems, Choi and Eastmann [5] introduced the concept of nanofluids in 1995 and showed experimentally that the embedding of nanoparticles (particles having diameters less than 100 nm) with the conventional fluids such as water, oil, or ethylene glycol mixture dramatically increases the thermal conductivity and heat transfer of the fluids. The resulting uniform dispersion and stable suspension of nanoparticles in the base fluid is referred to as a nanofluid. Due to their improved heat transfer properties, nanofluids are working fluids that have been used in many energy saving activities of modern industries, cost saving production of technological devices, and advanced life-saving medical treatments. Consequently, several studies have been conducted to study the boundary layer flow behavior of nanofluids. For instance, Ishak et al. [6] anticipated the boundary layer flow of nanofluids past a porous shrinking sheet. Ibrahim and Shanker [7] reported a numerical study of unsteady laminar boundary layer flow and heat transfer of a viscous incompressible fluid over stretching sheet. Sekhar et al. [8] studied the chemical reaction effects on unsteady MHD oscillatory slip flow in an optically thin fluid through a planer channel filled with saturated porous medium in the presence of a temperature-dependent heat source. Haile and Shankar [9] presented the boundary layer flow of a nanofluid through a porous medium subjected to magnetic field, thermal radiation, viscous dissipation, and chemical reaction effects. Ali et al. [10] applied the idea of Caputo–Fabrizio time fractional derivatives to magnetohydrodynamics free convection flow of generalized Walters’-B fluid over a static vertical plate. The unsteady MHD free convection flows of alumina-water and single-walled carbon nanotube-water nanofluids within a symmetrical wavy trapezoidal enclosure was reported by Job et al. [11]. On the other hand, Sathish Kumar et al. [12] investigated the effect of nonlinear thermal radiation on unsteady MHD flow between parallel plates. Saqib et al. [13] used the Laplace transform technique to obtain exact solutions for natural convection flow model of a hybrid nanofluid in two vertical infinite parallel plates. Other recent studies [14–16] on the effects of various thermophysical parameters on the boundary layer flow of nanoparticles have also been reported.In 1929, Williamson [17] noticed that there are some fluid flows that show both viscous and elastic properties, now called pseudoplastic fluids. He then proposed an equation, called the Williamson model, and verified it experimentally to exhibit the shear thinning or pseudoplastic properties. The boundary layer flow of pesudoplastic fluids is of great interest due to its application in industry such as extrusion of polymer sheets, emulsion coated sheets like photographic films, solutions, melts of high molecular weight polymers, etc [18]. In 2013, Nadeem et al. [19] presented the modeling of a two-dimensional boundary layer equation for the flow of Williamson fluid past a linear and exponentially stretching sheet. Following these pioneer works, a number of researchers have been reporting useful results on boundary layer flow and transport mechanisms of Williamson nanofluids with various thermophysical effects. For instance, Nadeem and Hussain [20] analyzed the flow and heat transfer mechanisms of Williamson nanofluids. The effects of radiation and chemical reaction on the steady boundary layer flow of MHD Williamson fluid through porous medium toward a horizontal linearly stretching sheet in the presence of nanoparticles were investigated numerically by Krishnamurthy et al. [21]. They also studied these effects on Williamson nanofluid slip flow over a stretching sheet embedded in a porous medium [22]. Recently, Reddy et al. [23] examined the magnetohydrodynamic boundary layer flow with heat and mass transfer of Williamson nanofluid over a stretching sheet with variable thickness and variable thermal conductivity under the radiation effect. Kho et al. [24] analyzed the boundary layer flow of Williamson nanofluids past over a stretching sheet in the presence of thermal radiation effect. Shawky et al. [25] also examined the MHD flow with heat and mass transfer of Williamson nanofluids over stretching sheet through porous medium. Ibrahim and Gamachu [26] utilized the spectral Quasilinearization method to inspect a nonlinear convection flow of Williamson nanofluid past a radially stretching surface under the application of electric field.However, to the best of the authors’ knowledge, no study has been reported on the analytic solution of time-dependent boundary layer flow of Williamson nanofluids over a permeable stretching sheet embedded in a porous medium with the effects of magnetic field, thermal radiation, and chemical reaction. Thus, motivated by the aforementioned works, the present study attempts to fill the existing gaps in this area. To this end, the work of Bibi et al. [27] has been generalized in a way that the flow model has been extended by considering effects of thermal radiation, chemical reaction and porosity of the medium. Moreover, analytic solutions were obtained by employing an efficient method, namely, the homotopy analysis method.


## 2. Flow Analysis and Mathematical Formulation

In this study, the effects of magnetic field, thermal radiation, and chemical reaction on a two-dimensional unsteady boundary layer flow of an incompressible and viscous Williamson nanofluid over a heated permeable stretching sheet embedded in a porous medium were examined. It was assumed that the magnetic Reynolds number is small in liquid metals and partially ionized fluids; and the effect of polarization of charges is not taken into account. So, the induced magnetic field is negligible in comparison with the applied magnetic field. It was also assumed that the flow is generated from a slit by stretching of a uniformly permeable and semi-infinite flat plate with one end fixed at the slit and embedded in an optically thick porous medium as shown in Figure1.Figure 1
Physical configuration of the flow model.The Cartesian coordinate system has been used in a way that thex-axis is along the stretching sheet, y-axis is normal to the sheet, the origin is located at the slit, and the flow in the region y≥0 is considered. Now, assuming that the sheet starts stretching at t=0 and extends horizontally with nonuniform velocity, Uw=ax/1−ct, where a and c are positive constants with dimension (time)−1, then the basic equations for the balance of mass, momentum, energy, and nanoparticle volume fraction of the flow problem can be expressed in vector form as follows.


### 2.1. Continuity Equation

The continuity equation for conservation of mass becomes(1)∇⋅V=0,where ∇ is the differential operator and V is the flow velocity vector.


### 2.2. Conservation of Momentum Equation

The Navier–Stoke’s equation for the balance of linear momentum is given by(2)ρf∂V∂t+V⋅∇V=∇⋅S+μ▽2V+J×B−μK0V,where ρf is the density of the nanofluid; t is time; S is the Cauchy stress tensor; μ is the coefficient of dynamic viscosity; J is the current density; B is the external magnetic field; J×B is the Lorentz force produced by the interaction of the applied magnetic field with velocity of the fluid; and K0 is permeability of the porous medium.


### 2.3. Conservation of Energy Equation

In the absence of heat source/sink, viscous dissipation, and Joule heating effects, the conservation of energy for heat transfer is given by(3)ρCPf∂T∂t+V⋅∇T=κ▽2T+ρCPpDB∇C⋅∇T+DTT∞∇T⋅∇T−∇qr,where CP is specific heat at constant pressure; ρCPf and ρCPp are the specific heat capacity of the nanofluid and the nanoparticles, respectively; T is temperature of the nanofluid; κ is thermal conductivity; DB is the Brownian diffusion coefficient; DT is the thermophoresis diffusion coefficient; T∞ is the ambient temperature; and qr is the radiative heat flux.


### 2.4. Conservation of Nanoparticle Concentration Equation

For a homogeneous chemical reaction, the concentration equation for nanoparticle volume fraction becomes(4)∂C∂t+V⋅∇C=DB▽2C+DTT∞▽2T−KrC−C∞,where C is nanoparticle volume fraction, and the rates Kr>0 and Kr<0 denote destructive and constructive reaction rates, respectively.In order to reduce the momentum (2), the Williamson constitutive equation is employed for the Cauchy stress tensor as defined in [28]:(5)S=−pI+τ1,where p is pressure; I is unit tensor; and τ1 is a deviatoric extra stress tensor given by(6)τ1=μ∞+μ0−μ∞1−Γγ˙A1,where μ0 and μ∞ are limiting viscosities at zero and infinite shear rates, respectively; Γ>0 is the material constant of the Williamson fluid; A1 is first Rivlin–Ericksen tensor; and γ˙=Π/2 is the shear rate with Π= trace A12. Here, for psedo-plastic fluids, we consider μ∞=0 and Γγ˙<1 which gives τ1=μ0/1−Γγ˙A1. Using binomial expansion, we get(7)τ1=μ01+Γγ˙A1.Details of the Williamson fluid model can be found in [17], [29], and [19].For an optically thick fluid, the energy (3) can also be reduced by applying the Rosseland diffusion approximation for the radiative heat flux given by(8)qr=−4σ∗3k∗∂T4∂y,where k* is the mean absorption constant and σ* is Stefan–Boltzmann constant. If the temperature difference within the flow is small, then we can expand T4 in Taylor series about T∞, which after neglecting higher order terms takes the form(9)T4≈4T∞3T−3T∞4.So, the heat flux can be approximated as(10)qr≈−16T∞3σ∗3k∗∂T∂y.Assume a time-dependent magnetic fieldBt=B0/1−ct is applied in the y-direction, where B0 is the initial magnetic field strength and ct<1. Then by using the boundary layer approximation, we rewrite the above conservation laws in two-dimensional Cartesian coordinates by using the equations given by Nadeem and Hussain [20] and Bibi et al. [27] as(11)∂u∂x+∂v∂y=0,(12)∂u∂t+u∂u∂x+v∂u∂y=υ∂2u∂y2+2υΓ∂u∂y∂2u∂y2−σB02ρf1−ct+υK0u,(13)∂T∂t+u∂T∂x+v∂T∂y=αf∂2T∂y2+τDB∂C∂y∂T∂y+DTT∞∂T∂y2+16σ∗T∞33ρCpfk∗∂2T∂y2,(14)∂C∂t+u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2−KrC−C∞,where u,v are the velocity components along x and y directions, respectively; υ=μ/ρf is the kinematic viscosity; σ is electric conductivity of the fluid; αf=κ/ρCpf is thermal diffusivity; and τ=ρCpp/ρCpf is ratio of effective heat capacities of nanoparticle and the ordinary fluid.The boundary conditions for the values of velocity, temperature, and nanoparticle concentration were taken as in Bibi et al. [27].Aty = 0,(15)u=Uwx,t=ax1−ct,v=Vwt=−V01−ct,T=Twx,t=T∞+T0Uwxυ1−ct,C=Cwx,t=C∞+C0Uwxυ1−ct,and as y⟶∞,(16)u⟶0,T⟶T∞,C⟶C∞,where Uw,Tw,Cw denote the surface velocity, temperature, and nanoparticle concentration, respectively; Vw is the mass transmission at the surface of the plate; a and a/1−ct are the initial and the effective stretching rates of the sheet, respectively; V0 is the constant value of velocity; T0,C0 are the constant values; and T∞,C∞ are the ambient values of temperature and nanoparticle concentration, respectively.In order to simplify the mathematical analysis, the following similarity transformations were introduced:(17)η=yUwυx,ψx,y,t=Uwυxfη,θη=T−T∞Tw−T∞,φη=C−C∞Cw−C∞,where η is the dimensionless similarity variable; fη is the dimensionless stream function; and θη and φη are the dimensionless temperature and nanoparticles volume fraction, respectively.If we choose the stream functionψx,y such that the velocity components are related as u=∂ψ/∂y and v=−∂ψ/∂x, then the continuity (11) is identically satisfied. Computing the required partial derivatives with respect to the new similarity variable η and substituting the values into equations (6)–(8), the governing system of partial differential equations are reduced to the following set of ordinary differential equations:(18)1+Wef″f‴+ff″−f′2−Af′+η2f″−M2+Kpf′=0,(19)1+43Rdθ″+Prfθ′−2f′θ−A23θ+ηθ′+Nbθ′φ′+Ntθ′2=0,(20)φ″+NtNbθ″+Scfφ′−2f′φ−A23φ+ηφ′−γφ=0,where the prime′ denotes differentiation with respect to η; We=Γx2a3/υ1−ct3 is the Weissenberg number representing the Williamson parameter; A=c/a is the unsteadiness parameter; M=σB02/aρf is the magnetic parameter; Kp=υx/UwK0 is the porosity parameter; Rd=4σ∗T∞3/κk∗ is the thermal radiation parameter; Pr=υ/αm is the Prandtl number; Nb=τDBCw−C∞/υ is the Brownian motion parameter; Nt=τDTTw−T∞/υT∞ is the thermophoresis parameter; Sc=υ/DB is the Schmidt number; and γ=Krx/Uw is the chemical reaction parameter with γ>0 and γ<0 denoting destructive and generative chemical reaction rates, respectively. Also using the similarity transformation in equation (17), the boundary conditions in equations (15) and (16) can be reduced as(21)f0=S,f′0=1,θ0=1,φ0=1,f′η⟶0,θη⟶0,φη⟶0,asη⟶∞,where S=V0/aυ is suction (for S<0) or injection (for s>0) parameter.From engineering point of view, it is useful to examine the impacts of skin friction coefficientCf, local Nusselt number Nux, and Sherwood number Shx on the boundary layer profiles of velocity, temperature, and concentration, respectively.The skin friction coefficientCf is defined as(22)cf=τw1/2ρfUw2,where τw is the shear stress at the permeable surface. For a Williamson fluid, it is given by(23)τw=−μ0∂u∂y+Γ2∂u∂y2y=0.Up on substitution, this gives(24)−12cfRex1/2=f″0+We2f″20,where Rex=ax2/υ1−ct is the local Reynolds number.The local Nusselt number is a dimensionless heat transfer coefficient defined as the ratio of convective and conductive heat transfer rates(25)Nux=xqwκTw−T∞,where qw=−κ+16σ∗T∞3/3k∗∂T/∂yy=0 is the heat flux that measures the heat transfer at the surface of a permeable stretching sheet. So, by substitution this gives(26)Nux=−1+43RdRex1/2θ′0.On the other hand, the local Sherwood numberShx is defined as the ratio of convective and diffusive mass transfer rates along a surface is which is given by(27)Shx=xJwDBCw−C∞,where Jw=−DB∂C/∂yy=0 is the mass flux at the surface. Thus, we get(28)Rex−1/2Shx=−φ′0.


## 2.1. Continuity Equation

The continuity equation for conservation of mass becomes(1)∇⋅V=0,where ∇ is the differential operator and V is the flow velocity vector.


## 2.2. Conservation of Momentum Equation

The Navier–Stoke’s equation for the balance of linear momentum is given by(2)ρf∂V∂t+V⋅∇V=∇⋅S+μ▽2V+J×B−μK0V,where ρf is the density of the nanofluid; t is time; S is the Cauchy stress tensor; μ is the coefficient of dynamic viscosity; J is the current density; B is the external magnetic field; J×B is the Lorentz force produced by the interaction of the applied magnetic field with velocity of the fluid; and K0 is permeability of the porous medium.


## 2.3. Conservation of Energy Equation

In the absence of heat source/sink, viscous dissipation, and Joule heating effects, the conservation of energy for heat transfer is given by(3)ρCPf∂T∂t+V⋅∇T=κ▽2T+ρCPpDB∇C⋅∇T+DTT∞∇T⋅∇T−∇qr,where CP is specific heat at constant pressure; ρCPf and ρCPp are the specific heat capacity of the nanofluid and the nanoparticles, respectively; T is temperature of the nanofluid; κ is thermal conductivity; DB is the Brownian diffusion coefficient; DT is the thermophoresis diffusion coefficient; T∞ is the ambient temperature; and qr is the radiative heat flux.


## 2.4. Conservation of Nanoparticle Concentration Equation

For a homogeneous chemical reaction, the concentration equation for nanoparticle volume fraction becomes(4)∂C∂t+V⋅∇C=DB▽2C+DTT∞▽2T−KrC−C∞,where C is nanoparticle volume fraction, and the rates Kr>0 and Kr<0 denote destructive and constructive reaction rates, respectively.In order to reduce the momentum (2), the Williamson constitutive equation is employed for the Cauchy stress tensor as defined in [28]:(5)S=−pI+τ1,where p is pressure; I is unit tensor; and τ1 is a deviatoric extra stress tensor given by(6)τ1=μ∞+μ0−μ∞1−Γγ˙A1,where μ0 and μ∞ are limiting viscosities at zero and infinite shear rates, respectively; Γ>0 is the material constant of the Williamson fluid; A1 is first Rivlin–Ericksen tensor; and γ˙=Π/2 is the shear rate with Π= trace A12. Here, for psedo-plastic fluids, we consider μ∞=0 and Γγ˙<1 which gives τ1=μ0/1−Γγ˙A1. Using binomial expansion, we get(7)τ1=μ01+Γγ˙A1.Details of the Williamson fluid model can be found in [17], [29], and [19].For an optically thick fluid, the energy (3) can also be reduced by applying the Rosseland diffusion approximation for the radiative heat flux given by(8)qr=−4σ∗3k∗∂T4∂y,where k* is the mean absorption constant and σ* is Stefan–Boltzmann constant. If the temperature difference within the flow is small, then we can expand T4 in Taylor series about T∞, which after neglecting higher order terms takes the form(9)T4≈4T∞3T−3T∞4.So, the heat flux can be approximated as(10)qr≈−16T∞3σ∗3k∗∂T∂y.Assume a time-dependent magnetic fieldBt=B0/1−ct is applied in the y-direction, where B0 is the initial magnetic field strength and ct<1. Then by using the boundary layer approximation, we rewrite the above conservation laws in two-dimensional Cartesian coordinates by using the equations given by Nadeem and Hussain [20] and Bibi et al. [27] as(11)∂u∂x+∂v∂y=0,(12)∂u∂t+u∂u∂x+v∂u∂y=υ∂2u∂y2+2υΓ∂u∂y∂2u∂y2−σB02ρf1−ct+υK0u,(13)∂T∂t+u∂T∂x+v∂T∂y=αf∂2T∂y2+τDB∂C∂y∂T∂y+DTT∞∂T∂y2+16σ∗T∞33ρCpfk∗∂2T∂y2,(14)∂C∂t+u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2−KrC−C∞,where u,v are the velocity components along x and y directions, respectively; υ=μ/ρf is the kinematic viscosity; σ is electric conductivity of the fluid; αf=κ/ρCpf is thermal diffusivity; and τ=ρCpp/ρCpf is ratio of effective heat capacities of nanoparticle and the ordinary fluid.The boundary conditions for the values of velocity, temperature, and nanoparticle concentration were taken as in Bibi et al. [27].Aty = 0,(15)u=Uwx,t=ax1−ct,v=Vwt=−V01−ct,T=Twx,t=T∞+T0Uwxυ1−ct,C=Cwx,t=C∞+C0Uwxυ1−ct,and as y⟶∞,(16)u⟶0,T⟶T∞,C⟶C∞,where Uw,Tw,Cw denote the surface velocity, temperature, and nanoparticle concentration, respectively; Vw is the mass transmission at the surface of the plate; a and a/1−ct are the initial and the effective stretching rates of the sheet, respectively; V0 is the constant value of velocity; T0,C0 are the constant values; and T∞,C∞ are the ambient values of temperature and nanoparticle concentration, respectively.In order to simplify the mathematical analysis, the following similarity transformations were introduced:(17)η=yUwυx,ψx,y,t=Uwυxfη,θη=T−T∞Tw−T∞,φη=C−C∞Cw−C∞,where η is the dimensionless similarity variable; fη is the dimensionless stream function; and θη and φη are the dimensionless temperature and nanoparticles volume fraction, respectively.If we choose the stream functionψx,y such that the velocity components are related as u=∂ψ/∂y and v=−∂ψ/∂x, then the continuity (11) is identically satisfied. Computing the required partial derivatives with respect to the new similarity variable η and substituting the values into equations (6)–(8), the governing system of partial differential equations are reduced to the following set of ordinary differential equations:(18)1+Wef″f‴+ff″−f′2−Af′+η2f″−M2+Kpf′=0,(19)1+43Rdθ″+Prfθ′−2f′θ−A23θ+ηθ′+Nbθ′φ′+Ntθ′2=0,(20)φ″+NtNbθ″+Scfφ′−2f′φ−A23φ+ηφ′−γφ=0,where the prime′ denotes differentiation with respect to η; We=Γx2a3/υ1−ct3 is the Weissenberg number representing the Williamson parameter; A=c/a is the unsteadiness parameter; M=σB02/aρf is the magnetic parameter; Kp=υx/UwK0 is the porosity parameter; Rd=4σ∗T∞3/κk∗ is the thermal radiation parameter; Pr=υ/αm is the Prandtl number; Nb=τDBCw−C∞/υ is the Brownian motion parameter; Nt=τDTTw−T∞/υT∞ is the thermophoresis parameter; Sc=υ/DB is the Schmidt number; and γ=Krx/Uw is the chemical reaction parameter with γ>0 and γ<0 denoting destructive and generative chemical reaction rates, respectively. Also using the similarity transformation in equation (17), the boundary conditions in equations (15) and (16) can be reduced as(21)f0=S,f′0=1,θ0=1,φ0=1,f′η⟶0,θη⟶0,φη⟶0,asη⟶∞,where S=V0/aυ is suction (for S<0) or injection (for s>0) parameter.From engineering point of view, it is useful to examine the impacts of skin friction coefficientCf, local Nusselt number Nux, and Sherwood number Shx on the boundary layer profiles of velocity, temperature, and concentration, respectively.The skin friction coefficientCf is defined as(22)cf=τw1/2ρfUw2,where τw is the shear stress at the permeable surface. For a Williamson fluid, it is given by(23)τw=−μ0∂u∂y+Γ2∂u∂y2y=0.Up on substitution, this gives(24)−12cfRex1/2=f″0+We2f″20,where Rex=ax2/υ1−ct is the local Reynolds number.The local Nusselt number is a dimensionless heat transfer coefficient defined as the ratio of convective and conductive heat transfer rates(25)Nux=xqwκTw−T∞,where qw=−κ+16σ∗T∞3/3k∗∂T/∂yy=0 is the heat flux that measures the heat transfer at the surface of a permeable stretching sheet. So, by substitution this gives(26)Nux=−1+43RdRex1/2θ′0.On the other hand, the local Sherwood numberShx is defined as the ratio of convective and diffusive mass transfer rates along a surface is which is given by(27)Shx=xJwDBCw−C∞,where Jw=−DB∂C/∂yy=0 is the mass flux at the surface. Thus, we get(28)Rex−1/2Shx=−φ′0.


## 3. Method of Solution

The homotopy analysis method (HAM), first proposed by Liao in 1992, is one of the relatively recent and powerful analytic methods that has shown great efficiency over the past few years. Based on the concept of homotopy in topology, the method provides us extremely great freedom to choose different solution expressions from which we can choose the one that approximates our solution more efficiently. In 1997, Liao introduced a very useful parameter, called convergence-control parameter, denoted byℏ or c0, which provides a more convenient approach to ensure the convergence of HAM solutions and to adjust the rate and region of convergence of the solution.Now, to implement the HAM in this problem, we choose the initial approximations in the form(29)f0η=S+1−e−η,θ0η=e−η,φ0η=e−η.We then select the auxiliary linear operators as(30)Lf=d3fdη3+d2fdη2,Lθ=d2θdη2+dθdη,Lφ=d2φdη2+dφdη,with the property LfC1+C2η+C3e−η=0,LθC4+C5e−η=0,LφC6+C7e−η=0; where Cii=1−7 are integral constants to be determined by the boundary conditions.We also choose the nonzero auxiliary functions as(31)Hfη=Hθη=Hφη=e−η.Finally, based on the ordinary differential equations (12)–(14), the nonlinear operators can be defined as follows:(32)Nf=1+We∂2ϕf∂η2∂3ϕf∂η3+ϕf∂2ϕf∂η2−∂ϕf∂η2−A∂ϕf∂η+η2∂2ϕf∂η2−M2+Kp∂ϕf∂η,Nθ=1+43Rd∂2ϕθ∂η2+Prϕf∂ϕθ∂η−2ϕθ∂ϕf∂η−A23ϕθ+η∂ϕθ∂η+Nb∂ϕθ∂η∂ϕφ∂η+Nt∂ϕθ∂η2,Nφ=∂2ϕφ∂η2+NtNb∂2ϕθ∂η2+Scϕf∂ϕφ∂η−2ϕφ∂ϕf∂η−A23ϕφ+η∂ϕφ∂η−γϕφ,where ϕfη;q,ϕθη;q, and ϕφη;q are the homotopy approximations of f,θ, and φ, respectively, satisfying(33)ϕf0;q=S,ϕf′0;q=1,ϕθ0;q=1,ϕφ0;q=1,ϕf′η;q⟶0,ϕθη;q⟶0,ϕφη;q⟶0,asη⟶∞,with q∈0,1 is called the embedding parameter.According to Liao, the corresponding zeroth-order deformation equations can be constructed as(34)1−qLfϕf−f0η=ℏfHfNfϕf,ϕθ,ϕφ,(35)1−qLθϕθ−θ0η=ℏθHθNθϕf,ϕθ,ϕφ,(36)1−qLφϕφ−φ0η=ℏφHφNφϕf,ϕθ,ϕφ.Obviously, whenq=0, the solutions to equations (25)–(27) are ϕfη;0=f0η,ϕθη;0=θ0η, and ϕφη;0=φ0η. Similarly, when q=1, the solutions to equations (25)–(27) are given by ϕfη;1=fη,ϕθη;1=θη, and ϕφη;1=φη. That is, as the embedding parameter q increases from 0 to 1, the homotopy solutions ϕfη;q,ϕθη;q, and ϕφη;q vary continuously from the known initial approximations f0η,θ0η, and φ0η to the unknown exact solutions fη,θη, and φη, respectively. Now, if this continuous variation is smooth enough, then we expand ϕfη;q,ϕθη;q, and ϕφη;q in the Taylor series about q as follows:(37)ϕfη;q=ϕfη;0+∑m=1∞1m!∂mϕfη;q∂ηmq=0qm,(38)ϕθη;q=ϕθη;0+∑m=1∞1m!∂mϕθη;q∂ηmq=0qm,(39)ϕφη;q=ϕφη;0+∑m=1∞1m!∂mϕφη;q∂ηmq=0qm.As noted in Liao [30], the convergence of the above homotopy series depend on the auxiliary parameters ℏf,ℏθ, and ℏφ. Assuming that these parameters are properly chosen so that the series in equations (28)–(30) converge at q=1. Then using the deformations of ϕfη;q,ϕθη;q, and ϕφη;q, the homotopy series solutions become(40)fη=ϕfη;1≈f0η+∑m=1∞fmη,(41)θη=ϕθη;1≈θ0η+∑m=1∞θmη,(42)φη=ϕφη;1≈φ0η+∑m=1∞φmη,where fmη=∂mϕfη;q/∂ηmq=0,θmη=∂mϕθη;q/∂ηmq=0 and φmη=∂mϕφη;q/∂ηmq=0 are called the mth-order homotopy derivatives. Here, we require these derivatives to satisfy the boundary conditions:(43)fm0=S,fm′0=1,θm0=1,φm0=1,fm′η⟶0,θmη⟶0,φmη⟶0,asη⟶∞.It is important to note that the equations in (31)–(33) provide us with the relationship between the exact solutions and the initial approximations.In order to get the terms of the series in equations (28)–(30), we first define the vectors satisfying the boundary conditions (43) as(44)f⟶m−1η=f0η,f1η,f2η,…,fm−1η,θ⟶m−1η=θ0η,θ1η,θ2η,…θm−1η,φ⟶m−1η=φ0η,φ1η,φ2η,…,φm−1η.Differentiating the zeroth-order deformation equations (25)–(27) m times with respect to q, then setting q=0 and finally dividing the resulting equations by m!, we get the so-called mth-order deformation equations given by(45)Lffm−χmfm−1=ℏfHfRm−1ff⟶m−1,θ⟶m−1,φ⟶m−1,(46)Lθθm−χmθm−1=ℏθHθRm−1θf⟶m−1,θ⟶m−1,φ⟶m−1,(47)Lφφm−χmφm−1=ℏφHφRm−1φf⟶m−1,θ⟶m−1,φ⟶m−1,where χm=0ifm≤11ifm>1 is the unit step function,(48)Rm−1f=fm−1‴+We∑k=0m−1fm−1″fm−1−k‴+∑k=0m−1fkfm−1−k″−∑k=0m−1fk′fm−1−1′−Afm−1′+η2fm−1″−M2+Kpfm−1′,Rm−1θ=1+43Rdθm−1″+Pr∑k=0m−1fkθm−1−k′−2∑k=0m−1θkfm−1−k′−A23θm−1+ηθm−1′+Nb∑k=0m−1θm−1′fm−1−k′+PrNt∑k=0m−1θk′θm−1−k′, (49)Rm−1φ=φm−1″+NtNbθm−1″+Sc∑r=0m−1frφm−1−r′−2∑r=0m−1φrfm−1−r′−A23φm−1+ηφm−1′−γφm−1.Taking the inverse of the linear operators on both sides of the high-order deformation equations (35)–(37), we obtain the following iterative formula:(50)fmη=χmfm−1η+ℏfLf−1HfηRm−1ff⟶m−1η,θ⟶m−1η,φ⟶m−1η,θmη=χmθm−1η+ℏθLθ−1HθηRm−1θf⟶m−1η,θ⟶m−1η,φ⟶m−1η,φmη=χmφm−1η+ℏφLφ−1HφηRm−1φf⟶m−1η,θ⟶m−1η,φ⟶m−1η.These equations can be solved recursively one after the other in order as a function of all the previous terms of the series. That is, for eachm=1,2,… , one can generate the terms of the series and the corresponding partial sums as orders of successive approximations.Even if the proper values of the initial approximations, the auxiliary linear operators, and the auxiliary functions have been selected, the HAM series solutions still contain the auxiliary parametersℏf,ℏθ, and ℏφ that could influence the accuracy, region of convergence, and rate of approximation for the solutions. That is, we still have great freedom to choose the values of the auxiliary parameters. So it remains to choose appropriate values of ℏf,ℏθ, and ℏφ that will ensure the convergence of the HAM solutions. In 2003, Liao [30] suggested that the HAM solutions will converge to the exact solutions in regions for which the graphs of the partial sums of the functions or their first few derivatives at some parameters against ℏ, called the ℏ-curves, are essentially horizontal. That is, the interval of least change in the physical quantity signifies the best interval for the choice of the convergence-control parameters.To facilitate the computation, the HAM-based Mathematica package, BVPh 2.0 developed by Zhao and Liao in 2013 has been employed. This package has recently been used successfully by many researchers like Farooq and Zhi-Liang [31]; Demir et al. [32]; and Bano et al. [33]. For details regarding use of the package, we recommend the readers to refer Zhao and Liao [34].It can be seen from Figure2 that the ℏ-curves are nearly horizontal in the ranges(51)−1.2<ℏf<−0.3,−1.2<ℏθ<−0.2,−1.2<ℏφ<−0.2.Figure 2
ℏ-curves (a) and total squared residual error (b) for the 20th HAM approximation.That is, these intervals are the valid regions for which the HAM solutions converge to the exact solutions. So, any value from these intervals may be picked up to give the approximate solutions of our flow problem. By comparing the results at some of these possible values of the parameters, one may choose the optimal values of the parameters for which the approximate solutions will converge in the whole region ofη. More precisely, the optimal values of the convergence-control parameters can be obtained by minimizing the average squared residual errors of the governing equations at the 20th order of HAM approximation as indicated in Figure 2(b). Now, using the BVPh 2.0 package at 20th order HAM approximation, the optimal values for the convergence-control parameters were found to be(52)ℏf≈−1.03433,ℏθ≈−0.879735,ℏφ≈−0.767268.These values and the valuesM=0.2,We=0.2,A=0.3,Pr=0.72,Nb=0.2,Nt=0.1,Sc=1,S=0.3,Kp=1,Rd=0.3andγ=2 are utilized throughout this study unless stated otherwise.It is clear from Table1 that increasing the number of terms decreases the squared residual errors and leads the series to converge.Table 1
Convergence of some HAM solutions.
Order of HAM approximation
−
f
″
0
−
θ
′
0
−
φ
′
0
Squaredεf
Residualεθ
Errorsεφ
5
1.9352
1.0616
1.3121
10
1.9370
1.0531
1.3064
1.84
×
10
−
6
0.0000388
0.0000201
15
1.9381
1.0511
1.3059
20
1.9381
1.0499
1.3055
7.35
×
10
−
8
0.0000192
0.0000104
25
1.9381
1.0492
1.3053
30
1.9381
1.0487
1.3051
2.28
×
10
−
8
0.0000129
0.0000071
35
1.9381
1.0484
1.3050
40
1.9381
1.0481
1.3049
9.91
×
10
−
9
0.0000098
0.0000054


## 4. Results and Discussions

In this section, we present and discuss the main results of our mathematical analysis for the unsteady boundary layer flow of Williamson nanofluids induced by a permeable stretching plate embedded in a porous medium in the presence of magnetic field, thermal radiation, and chemical reaction. In particular, the effects of various thermophysical parameters on velocity, temperature, and concentration profiles were examined and presented as follows.Since there are many theoretically and practically important flow field and transport phenomena that can be unsteady due to sudden movement of surfaces, we begin our result analysis by presenting the effects of unsteadiness parameter on our flow problem.The result in Figure3 explains that the velocity gradient, temperature, and concentration profiles are decreasing functions of the unsteadiness parameter in the boundary layer. This is due to the fact that as the unsteadiness parameter A=c/a increases, the velocity Uw=ax/1−ct of the stretching sheet also decreases which causes the transfer of less amount of heat and mass from the sheet to the nanofluid in the boundary layer region.Figure 3
Impacts of unsteadiness parameter on velocity, temperature, and concentration profiles.The influence of magnetic field on the unknown functions of the given problem is illustrated in Figure4. It reveals that the increase in magnetic parameter increases the temperature and concentration profiles but reduces the velocity gradient.Figure 4
Effects of magnetic parameter on velocity gradient, temperature, and concentration profiles.This result agrees well with the fact that the increase in magnetic field strength tends to increase a dragging force which resists the motion of the fluid particles and produces heat in the boundary layer.The effect of Williamson viscoelastic parameter in the flow field can be described in terms of the Weissenberg numberWe as illustrated in Figure 5.Figure 5
Effects of Weissenberg numberWe on velocity, temperature, and concentration profiles.The result in Figure5 reveals that the velocity gradient decreases near the stretching sheet as the Weissenberg number We increases. The temperature and concentration profiles are increasing with increasing of We in the boundary layer.The effect of porosityKp in the flow problem has been investigated as presented in Figure 6.Figure 6
Impacts of porosity parameterKp on velocity gradient, temperature, and concentration profiles.Figure6 indicates that as porosity of the medium increases, the temperature and concentration profiles also increase in the boundary layer. On the other hand, the velocity gradient is found to decrease near the stretching surface and increase away from the surface with increasing porosity of the medium. This is due to the fact that increasing the value of Kp has the tendency to resist the flow.The result in Figure7 displays that the rise in thermal radiation increases the temperature but reduces the concentration profiles. This is because the increase in thermal radiation provides more heat to the nanofluid which causes the increase in temperature and thermal boundary layer thickness. Furthermore, this parameter has no significant effect on the velocity gradient of the fluid flow.Figure 7
Effects of thermal radiation parameterRd on velocity gradient, temperature, and concentration profiles.The chemical reaction that takes place in nanofluid flow results the conversion of species in the boundary layer. The effects of chemical reaction parameter on the concentration profile of our problem are illustrated in Figure8.The result in Figure8 shows that the concentration profile increases with the generative chemical reaction and decreases with the destructive chemical reaction. This is true because as generative and destructive chemical reactions take place, the amount of nanoparticles within the fluid is getting larger and smaller, respectively.Figure 8
Effects of chemical reaction parameterγ on concentration profile.We see from Figure9 that the velocity gradient, temperature profile, and concentration of nanoparticles are decreasing functions of the injection parameter.Figure 9
Effects of surface permeability parameterS on velocity gradient and temperature profile.It can be observed from Table2 that all the coefficients of skin friction, Nusselt number, and Sherwood number are increasing functions of A. Also, the skin friction coefficient can be maximized by increasing M,Kp, or γ; and the coefficient of Nusselt number is increased by decreasing Kp,Rd, or the generative chemical reaction parameter γ. It is also shown that increasing Rd or γ or decreasing Kp raises the coefficient of Sherwood number.Table 2
Coefficients of skin friction, Nusselt number, and Sherwood number.
A
M
K
p
R
d
γ
−
f
″
0
−
θ
′
0
−
φ
′
0
0.3
1.93811
1.04995
1.30552
0.6
2.04031
1.16875
1.42042
0.9
0.2
2.14202
1.27181
1.52403
0.8
2.21440
1.02165
1.27691
1.4
2.91874
1.21627
1.45894
1
2.44662
1.25545
1.50555
4
3.05318
1.21715
1.46277
7
3.33303
1.46380
1.78493
0.1
3.41649
1.32197
1.60133
0.4
3.33312
1.20886
1.88957
0.7
0.2
3.33277
1.03745
1.95725
0.5
3.33281
1.04622
2.06716
0.8
3.33283
1.04553
2.17143
−0.2
3.33275
1.03952
1.78885
−0.5
3.33271
1.04150
1.64232
−0.8
3.33263
1.04410
1.46839We further validate our solutions by comparing them with previously published works in some limited cases. IfWe=M=Rd=Kp=Rd=Sc=S=γ=0, then our analytic approximation to the skin friction coefficient, −f″0, with some values of the unsteadiness parameter A can be compared with the results of some previous studies as shown in Table 3.The comparison presented in Table3 reveals that for the selected values of the parameters, the values of −f″0 determined in this study are in excellent agreement with some of the previously published works.Table 3
Comparisons of results on−f″0 for We=M=Kp=Rd=Sc=S=γ=0,Pr=0.72,Nb=0.2,Nt=0.1 against some values of A.
A
Sharidan et al. [35]
Chamkha et al. [36]
Mukhopadhyay and Gorla. [37]
Khan and Azam [38]
Bibi et al. [27]
Present study
0.0
—
—
—
1.00000
1.0005
1.00000
0.2
—
—
—
1.06801
1.0685
1.06874
0.4
—
—
—
1.13469
1.1349
1.13521
0.6
—
—
—
1.19912
1.1992
1.19930
0.8
1.26104
1.261512
1.261479
1.26104
—
1.26099
1.2
1.37772
1.378052
1.377850
1.37772
—
1.37755
2.0
—
—
—
1.58737
—
1.58740


## 5. Concluding Remarks

In this study, the homotopy analysis method (HAM) has been successfully applied to obtain analytical series solutions for the unsteady boundary layer flow of Williamson nanofluids induced by a permeable stretching sheet embedded in a porous medium in the presence of magnetic field, thermal radiation, and chemical reaction. The convergence of the HAM solutions is ensured by using successive iterations,ℏ-curves, and the squared residual error analysis. The validity of the results are also verified by comparing them with previous results under some restricted assumptions and found to be in excellent agreement. The effects of various pertinent thermophysical parameters on the dimensionless velocity, temperature, and concentration profiles are examined, and the results are summarized as follows:(i)
Since the fluid flow is caused only by the stretching sheet, the velocity, temperature, and concentration boundary layer thickness generally decrease as we move away from the surface of the sheet.(ii)
The velocity gradient of the flow problem is found to be a decreasing function of the unsteadiness parameter, magnetic field strength, and injection parameter. It also decreases near the stretching surface with the increase in porosity of the medium and Weissenberg number.(iii)
The temperature profile can be increased in the boundary layer by increasing magnetic field strength, porosity of the medium, thermal radiation, or Weissenberg number. This profile can also be increased by decreasing the injection parameter or the unsteadiness parameter.(iv)
The concentration profile of nanoparticles can be increased in the boundary layer by increasing the magnetic field strength, Weissenberg number, porosity of the medium, or generative chemical reaction. This profile can also be increased by decreasing the unsteadiness parameter, destructive chemical reaction, thermal radiation or injection parameter;(v)
The rate of momentum transfer can be enhanced mainly by decreasing magnetic field strength, porosity of the medium, chemical reaction, or unsteadiness parameter; it can also be decreased by reducing the effect of thermal radiation.(vi)
Heat transfer rate can be speeded up by increasing the unsteadiness parameter or by decreasing the porosity of the medium, thermal radiation, or destructive chemical reaction.(vii)
Mass transfer rate can also be enhanced by increasing the thermal radiation, constructive chemical reaction, or unsteadiness parameter. It can also be facilitated by decreasing the effect of thermal radiation or destructive chemical reaction.



---

*Source: 1890972-2020-02-01.xml*