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

Kebede, Tesfaye; Haile, Eshetu; Awgichew, Gurju; Walelign, Tadesse

doi:https://doi.org/10.1155/2020/1890972

Journal of Applied Mathematics

On this page

Abstract Introduction Results Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 1890972 | https://doi.org/10.1155/2020/1890972

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

Tesfaye Kebede,¹Eshetu Haile,¹Gurju Awgichew,¹and Tadesse Walelign¹

Academic Editor: Mustafa Inc

Received28 Jun 2019

Revised09 Sept 2019

Accepted03 Jan 2020

Published01 Feb 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.

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 Figure 1.

The Cartesian coordinate system has been used in a way that the x-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 is considered. Now, assuming that the sheet starts stretching at and extends horizontally with nonuniform velocity, , where a and c are positive constants with dimension (time)⁻¹, 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 becomeswhere 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 bywhere 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; is the Lorentz force produced by the interaction of the applied magnetic field with velocity of the fluid; and 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 bywhere is specific heat at constant pressure; and are the specific heat capacity of the nanofluid and the nanoparticles, respectively; T is temperature of the nanofluid; κ is thermal conductivity; is the Brownian diffusion coefficient; is the thermophoresis diffusion coefficient; is the ambient temperature; and is the radiative heat flux.

2.4. Conservation of Nanoparticle Concentration Equation

For a homogeneous chemical reaction, the concentration equation for nanoparticle volume fraction becomeswhere C is nanoparticle volume fraction, and the rates and 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]:where p is pressure; I is unit tensor; and is a deviatoric extra stress tensor given bywhere and are limiting viscosities at zero and infinite shear rates, respectively; is the material constant of the Williamson fluid; is first Rivlin–Ericksen tensor; and is the shear rate with trace . Here, for psedo-plastic fluids, we consider and which gives Using binomial expansion, we get

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 bywhere is the mean absorption constant and is Stefan–Boltzmann constant. If the temperature difference within the flow is small, then we can expand in Taylor series about , which after neglecting higher order terms takes the form

So, the heat flux can be approximated as

Assume a time-dependent magnetic field is applied in the y-direction, where is the initial magnetic field strength and . 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] aswhere are the velocity components along x and y directions, respectively; is the kinematic viscosity; σ is electric conductivity of the fluid; is thermal diffusivity; and 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].

At y = 0,and as ,where denote the surface velocity, temperature, and nanoparticle concentration, respectively; is the mass transmission at the surface of the plate; a and are the initial and the effective stretching rates of the sheet, respectively; is the constant value of velocity; are the constant values; and are the ambient values of temperature and nanoparticle concentration, respectively.

In order to simplify the mathematical analysis, the following similarity transformations were introduced:where η is the dimensionless similarity variable; is the dimensionless stream function; and and are the dimensionless temperature and nanoparticles volume fraction, respectively.

If we choose the stream function such that the velocity components are related as and , 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:where the prime′ denotes differentiation with respect to η; is the Weissenberg number representing the Williamson parameter; is the unsteadiness parameter; is the magnetic parameter; is the porosity parameter; is the thermal radiation parameter; is the Prandtl number; is the Brownian motion parameter; is the thermophoresis parameter; is the Schmidt number; and is the chemical reaction parameter with and 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 aswhere is suction (for ) or injection (for ) parameter.

From engineering point of view, it is useful to examine the impacts of skin friction coefficient , local Nusselt number , and Sherwood number on the boundary layer profiles of velocity, temperature, and concentration, respectively.

The skin friction coefficient is defined aswhere is the shear stress at the permeable surface. For a Williamson fluid, it is given by

Up on substitution, this giveswhere 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 rateswhere is the heat flux that measures the heat transfer at the surface of a permeable stretching sheet. So, by substitution this gives

On the other hand, the local Sherwood number is defined as the ratio of convective and diffusive mass transfer rates along a surface is which is given bywhere is the mass flux at the surface. Thus, we get

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 , 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

We then select the auxiliary linear operators aswith the property ; where are integral constants to be determined by the boundary conditions.

We also choose the nonzero auxiliary functions as

Finally, based on the ordinary differential equations (12)–(14), the nonlinear operators can be defined as follows:where , and are the homotopy approximations of , and φ, respectively, satisfyingwith is called the embedding parameter.

According to Liao, the corresponding zeroth-order deformation equations can be constructed as

Obviously, when , the solutions to equations (25)–(27) are , and . Similarly, when , the solutions to equations (25)–(27) are given by , and . That is, as the embedding parameter q increases from 0 to 1, the homotopy solutions , and vary continuously from the known initial approximations , and to the unknown exact solutions , and , respectively. Now, if this continuous variation is smooth enough, then we expand , and in the Taylor series about q as follows:

As noted in Liao [30], the convergence of the above homotopy series depend on the auxiliary parameters , and . Assuming that these parameters are properly chosen so that the series in equations (28)–(30) converge at . Then using the deformations of , and , the homotopy series solutions becomewhere and are called the mth-order homotopy derivatives. Here, we require these derivatives to satisfy the boundary conditions:

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

Differentiating the zeroth-order deformation equations (25)–(27) m times with respect to q, then setting and finally dividing the resulting equations by , we get the so-called mth-order deformation equations given bywhere is the unit step function,

Taking the inverse of the linear operators on both sides of the high-order deformation equations (35)–(37), we obtain the following iterative formula:

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 each , 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 , 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 , 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 Figure 2 that the -curves are nearly horizontal in the ranges

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

These values and the values are utilized throughout this study unless stated otherwise.

It is clear from Table 1 that increasing the number of terms decreases the squared residual errors and leads the series to converge.

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 Figure 3 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 increases, the velocity 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.

The influence of magnetic field on the unknown functions of the given problem is illustrated in Figure 4. It reveals that the increase in magnetic parameter increases the temperature and concentration profiles but reduces the velocity gradient.

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 number as illustrated in Figure 5.

The result in Figure 5 reveals that the velocity gradient decreases near the stretching sheet as the Weissenberg number increases. The temperature and concentration profiles are increasing with increasing of in the boundary layer.

The effect of porosity in the flow problem has been investigated as presented in Figure 6.

Figure 6 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 has the tendency to resist the flow.

The result in Figure 7 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.

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 Figure 8.

The result in Figure 8 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.

We see from Figure 9 that the velocity gradient, temperature profile, and concentration of nanoparticles are decreasing functions of the injection parameter.

It can be observed from Table 2 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 , or γ; and the coefficient of Nusselt number is increased by decreasing , or the generative chemical reaction parameter γ. It is also shown that increasing or γ or decreasing raises the coefficient of Sherwood number.

We further validate our solutions by comparing them with previously published works in some limited cases. If , then our analytic approximation to the skin friction coefficient, , 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 Table 3 reveals that for the selected values of the parameters, the values of determined in this study are in excellent agreement with some of the previously published works.

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.

Data Availability

The parametric variation data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

H. Alfvén, “Existence of electromagnetic-hydrodynamic waves,” Nature, vol. 150, no. 3805, pp. 405-406, 1942.
View at: Publisher Site | Google Scholar
B. C. Sakiadis, “Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow,” AIChE Journal, vol. 7, no. 1, pp. 26–28, 1961.
View at: Publisher Site | Google Scholar
R. Kumar, “MHD boundary layer flow on heat and mass transfer over a stretching sheet with slip effect,” Journal of Naval Architecture and Marine Engineering, vol. 10, no. 2, pp. 109–118, 2013.
View at: Publisher Site | Google Scholar
N. S. Kumar, R. Kumar, and A. G. Vijaya Kumar, “Thermal diffusion and chemical reaction effects on unsteady flow past a vertical porous plate with heat source dependent in slip flow regime,” Journal of Naval Architecture and Marine Engineering, vol. 13, no. 1, pp. 51–62, 2016.
View at: Publisher Site | Google Scholar
S. U. S. Choi and J. A. Eastman, “Enhancing thermal conductivity of fluids with nanoparticles,” in Proceedings of the 1995 International Mechanical Engineering Congress and Exhibition, San Francisco, CA, USA, November 1995.
View at: Google Scholar
A. Ishak, Y. Y. Lok, and I. Pop, “Stagnation-point flow over a shrinking sheet in a micropolar fluid,” Chemical Engineering Communications, vol. 197, no. 11, pp. 1417–1427, 2010.
View at: Publisher Site | Google Scholar
W. Ibrahim and B. Shanker, “Unsteady boundary layer flow and heat transfer due to a stretching sheet by quasilinearization technique,” World Journal of Mechanics, vol. 1, no. 6, pp. 288–293, 2011.
View at: Publisher Site | Google Scholar
K. Raja Sekhar, G. V. Ramana Reddy, and B. D. C. N. Prasad, “Chemically reacting on MHD oscillatory slip flow in a planer channel with varying temperature and concentration,” Advances in Applied Science Research, vol. 3, no. 5, pp. 2652–2659, 2012.
View at: Google Scholar
E. Haile and B. Shankar, “Heat and mass transfer through a porous media of MHD flow of nanofluids with thermal radiation, viscous dissipation and chemical reaction effects,” American Chemical Science Journal, vol. 4, no. 6, pp. 828–846, 2014.
View at: Publisher Site | Google Scholar
F. Ali, M. Saqib, I. Khan, and N. Ahmad Sheikh, “Application of caputo-fabrizio derivatives to MHD free convection flow of generalized Walters’-B fluid model,” The European Physical Journal Plus, vol. 131, no. 10, 2016.
View at: Publisher Site | Google Scholar
V. M. Job, S. R. Gunakala, B. Rushi Kumar, and R. Sivaraj, “Time-dependent hydromagnetic free convection nanofluid flows within a wavy trapezoidal enclosure,” Applied Thermal Engineering, vol. 115, pp. 363–377, 2016.
View at: Publisher Site | Google Scholar
M. Sathish Kumar, N. Sandeep, and B. Rushi Kumar, “Effect of nonlinear thermal radiation on unsteady MHD flow between parallel plates,” Global Journal of Pure and Applied Mathematics, vol. 12, p. 1, 2016.
View at: Google Scholar
M. Saqib, I. Khan, and S. Shafie, “Application of fractional differential equations to heat transfer in hybrid nanofluid: modeling and solution via integral transforms,” Advances in Differential Equations, vol. 2019, no. 1, 2019.
View at: Publisher Site | Google Scholar
M. Saqib, F. Ali, I. Khan, N. A. Sheikh, and A. Khan, “Entropy generation in different types of fractionalized nanofluids,” Arabian Journal for Science and Engineering, vol. 44, no. 1, pp. 531–540, 2019.
View at: Publisher Site | Google Scholar
M. Saqib, I. Khan, and S. Shafie, “Natural convection channel flow of CMC-based CNTs nanofluid,” The European Physical Journal Plus, vol. 133, p. 549, 2018.
View at: Publisher Site | Google Scholar
M. Saqib, F. Ali, I. Khan, N. Ahmad Sheikh, and S. Bin Shafie, “Convection in ethylene glycol-based molybdenum disulfide nanofluid: atangana–baleanu fractional derivative approach,” Journal of Thermal Analysis and Calorimetry, vol. 135, no. 1, pp. 523–532, 2019.
View at: Publisher Site | Google Scholar
R. V. Williamson, “The flow of pseudoplastic materials,” Industrial & Engineering Chemistry, vol. 21, no. 11, pp. 1108–1111, 1929.
View at: Publisher Site | Google Scholar
N. C. Sacheti, T. E. Bashir, and P. Chandran, “Perturbation analysis of 2-dimensional boundary layer flow of an inelastic fluid using Williamson model,” International Journal of Applied Engineering Research, vol. 22, pp. 12728–12734, 2017.
View at: Google Scholar
S. Nadeem, S. T. Hussain, and C. Lee, “Flow of a Williamson fluid over a stretching sheet,” Brazilian Journal of Chemical Engineering, vol. 30, no. 3, pp. 619–625, 2013.
View at: Publisher Site | Google Scholar
S. Nadeem and S. T. Hussain, “Flow and heat transfer analysis of Williamson nanofluid,” Applied Nanoscience, vol. 4, no. 8, pp. 1005–1012, 2014.
View at: Publisher Site | Google Scholar
M. R. Krishnamurthy, B. C. Prasannakumara, B. J. Gireesha, and R. S. R. Gorla, “Effect of chemical reaction on MHD boundary layer flow and melting heat transfer of Williamson nanofluid in porous medium,” Engineering Science and Technology, an International Journal, vol. 19, no. 1, pp. 53–61, 2016.
View at: Publisher Site | Google Scholar
B. C. Prasannakumara, B. J. Gireesha, R. S. R. Gorla, and M. R. Krishnamurthy, “Effects of chemical reaction and nonlinear thermal radiation on Williamson nanofluid slip flow over a stretching sheet embedded in a porous medium,” Journal of Aerospace Engineering, vol. 29, no. 5, Article ID 04016019, 2016.
View at: Publisher Site | Google Scholar
S. Reddy, K. Naikoti, and M. M. Rashidi, “MHD flow and heat transfer characteristics of Williamson nanofluid over a stretching sheet with variable thickness and variable thermal conductivity,” Transactions of A. Razmadze Mathematical Institute, vol. 171, no. 2, pp. 195–211, 2017.
View at: Publisher Site | Google Scholar
Y. B. Kho, A. Hussanan, M. K. A. Mohamed, N. M. Sarif, Z. Ismail, and M. Z. Salleh, “Thermal radiation effect on MHD flow and heat transfer analysis of Williamson nanofluid past over a stretching sheet with constant wall temperature,” Journal of Physics: Conference Series, vol. 890, Article ID 012034, 2017.
View at: Publisher Site | Google Scholar
H. M. Shawky, N. T. M. Eldabe, K. A. Kamel, and E. A. Abd-Aziz, MHD Flow with Heat and Mass Transfer of Williamson Nanofluid over Stretching Sheet through Porous Medium, Springer, Berlin, Germany, 2018.
W. Ibrahima and D. Gamachu, “Nonlinear convection flow of Williamson nanofluid past a radially stretching surface,” AIP Advances, vol. 9, no. 8, Article ID 085026, 2019.
View at: Publisher Site | Google Scholar
M. Bibi, K. U. Rehman, M. Y. Malik, and M. Tahir, “Numerical study of Unsteady Williamson fluid flow and heat transfer in the presence of MHD through a permeable stretching surface,” The European Physical Journal Plus, vol. 133, no. 4, p. 154, 2018.
View at: Publisher Site | Google Scholar
I. Dapra and G. Scarpi, “Perturbation solution for pulsatile flow of a non-Newtonian Williamson fluid in a rock fracture,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 2, pp. 271–278, 2007.
View at: Publisher Site | Google Scholar
C. Vasudev, U. R. Rao, M. V. S. Reddy, and G. P. Rao, “Peristaltic pumping of Williamson fluid through a porous medium in a horizontal channel with heat transfer,” American Journal of Scientific and Industrial Research, vol. 1, no. 3, pp. 656–666, 2010.
View at: Publisher Site | Google Scholar
S. J. Liao, Beyond Perturbation: Introduction to Homotopy Analysis Method, Chapman and Hall/CRC, Boca Raton, FL, USA, 2003.
U. Farooq and L. Z. Liang, “Nonlinear heat transfer in a two-layer flow with nanofluids by OHAM,” Journal of Heat Transfer, vol. 136, no. 2, 2014.
View at: Publisher Site | Google Scholar
M. Ş. Demir, M. S. Dokuz, and S. Barış, “Semi-analytical solution for the MHD three dimensional stagnation-point flow of a Newtonian fluid towards a uniformly heated and moving vertical plate,” Journal of Environmental Science, Computer Science and Engineering & Technology, vol. 5, pp. 39–51, 2016.
View at: Google Scholar
N. Bano, B. B. Singh, and S. R. Sayyed, “Homotopy analysis for MHD hiemenz flow in a porous medium with thermal radiation, velocity and thermal slips effects,” Frontiers in Heat and Mass Transfer, vol. 10, 2018.
View at: Publisher Site | Google Scholar
Y. Zhao and S. J. Liao, “HAM-based mathematica package BVPh 2.0 for nonlinear boundary value problems,” Advances in the Homotopy Analysis Method, vol. 361, 2013.
View at: Google Scholar
S. Sharidan, T. Mahmood, and I. Pop, “Similarity solutions for the unsteady boundary layer flow and heat transfer due to a stretching sheet,” International Journal of Applied Mechanics and Engineering, vol. 11, pp. 647–654, 2006.
View at: Google Scholar
A. J. Chamkha, A. M. Aly, and M. A. Mansour, “Similarity solution for unsteady heat and mass transfer from a stretching surface embedded in a porous medium with suction/injection and chemical reaction effects,” Chemical Engineering Communications, vol. 197, no. 6, pp. 846–858, 2010.
View at: Publisher Site | Google Scholar
S. Mukhopadhyay and R. S. R. Gorla, “Unsteady MHD boundary layer flow of an upper convected Maxwell fluid past a stretching sheet with first order constructive/destructive chemical reaction,” Journal of Naval Architecture and Marine Engineering, vol. 9, no. 2, pp. 123–133, 2012.
View at: Publisher Site | Google Scholar
M. Khan and M. Azam, “Unsteady heat and mass transfer mechanisms in MHD Carreau nanofluid flow,” Journal of Molecular Liquids, vol. 225, pp. 554–562, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Tesfaye Kebede et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

2588

Downloads

1610

Citations