Non-similar forced convection analysis of Oldroyd-B fluid flow over an exponentially stretching surface

Abstract

In the study of boundary layer regions, it is in practice to dimensionalize the governing system and grouping variables together into dimensionless quantities in order to curtail the total number of variables. In similar flow phenomenon the physical parameters do not vary along the streamwise direction. However in non-similar flows the physical quantities change in the streamwise direction. In non-similar flows we are forced to non-dimensionalize the governing equations through non-similarity transformations. The forced flow of Oldroyd-B fluid is initiated as a result of stretching of a surface at an exponential rate. Flows over stretching surfaces are important because of their applications in extrusion processes. The forthright purpose of this study is to consider the non-similar aspects of forced convection from flat heated surface subjected to external viscoelastic fluid flow, described by the freely growing boundary layers enclosed by a region that involves without velocity and temperature gradients. The governing system of nonlinear partial differential equations (PDE’s) is transformed into dimensionless form by proposing new non-similar transformations. The dimensionless partial differential system is solved by using local non-similarity via bvp4c. Thermal transport analysis is conducted for distinct values of dimensionless numbers. It is revealed that heat shifting process expanded by the increase in the numerical values of Prandtl number and relaxation time. The dimensionless convective heat transfer coefficient results revealed that it is declining by expanding relaxation time constant $β_{1}$ and a boost is observed by enlarging the Pr and retardation time constant $β_{2}$ . A comparison of Nusselt number is presented.

Keywords

Exponential stretching sheet local non-similarity technique Oldroyd-B fluid non-similar modeling

Introduction

The analysis of boundary layer regions has attained an increasing interest and greatest importance in the theory of fluid mechanics in view of its popular applications. The boundary layer conception was inaugurated by Ludwig Prandtl in 1904. An assumption was made for two regions of fluid where viscosity consequences boosted inner part and insignificant for outer path of the boundary. At the continuous operating rigid sheet, a flow of two-dimensional, Newtonian fluid was examined by Sakiadis.¹ He took initiative of the boundary layer flow, which later on developed a huge interest in researchers working in this discipline. Crane² expanded the theory of Sakiadis.¹ He exhibited the exact solution of quiescent fluid above an unwavering expanding sheet. The stagnation location in flow above an expanding sheet has been considered by Chiam,³ where he inspected particular flow configuration.

The analysis of non-Newtonian fluids is on the ruling edge in research during last few decades because of its warm welcome in technological and scientific applications as compared to Newtonian fluids. There are many items such as food products (alcoholic beverages, ketchup, ice cream, yogurt and mayonnaise, etc.), chemical products (tooth pastes, paints, pharmaceutical chemicals, cosmetics, shampoos, etc.), and biological stuffs (blood, syrups, synovial fluid and vaccine, etc.) where Newtonian’s law of viscosity is not implemented, so they are distinguished as the non-Newtonian fluids.⁴ The industrial fluids differ in nature from viscous fluids due to their distinct characteristics, for example dying of paper and textiles plastic manufacturing, movement of biological fluids, and blade coating. The description of these materials is much more elaborated, due to this reason it is very difficult to compose a single constitutive relationship which assume the properties of such items. For non-Newtonian fluids, a lot of constitutive relations have been created to find their adequate conduct. Non-Newtonian fluids are organized as integral, rate, and differential.⁵ Oldroyd-B fluid falls in the category of rate type non-Newtonian fluid. This category interprets the relaxation and retardation time consequences. Tanner,⁶ Waters and King⁷ were the founders, they obtained exact solutions for Oldroyd-B fluids. Vieru et al.⁸ studied the flow of generic visco-elastic Oldroyd-B fluid for continuously accelerating plate. Fetecau et al.⁹ deliberated energetic balance for Rayleigh-Stokes problem of Oldroyd-B fluid. Unsteady helical flow of Oldroyd-B fluid presented by Jamil et al.¹⁰ Three-dimensional convective flow of Oldroyd-B fluid over expanding sheet was explored by Hayat et al.¹¹ Ramzan et al.¹² examined heating effects for three dimensional flow of Oldroyd-B fluid. Hayat et al.¹³ presented analysis for flow of Oldroyd-B fluid based on homogeneous-heterogeneous response and Cattaneo-Christov heat flux. Mixed convective three-dimensional radiative flow of Oldroyd-B fluid studied by Shehzad et al.¹⁴ Zheng et al.¹⁵ exhibited MHD flow of Oldroyd-B fluids with slip consequences.

The combination of polymer and metal sheets develop the stretching surface through polymer extrusion process.⁴ Flow over stretching sheets have secured prominent value in the investigations of boundary layers because of their promising implementations in scientific, engineering medical, and industrial fields such as crystal extension, glass industry, hot rippling, silicon wafer process, polymer physics, space vehicles, paper manufacturing, spinning of fibers, and continuous casting of metals. Exponential stretching in surfaces is defined as expansion in surfaces at an exponential rate.¹⁶ It is a way to deal with flow or heat shifting by using distinct exteriors such that thermal transport over continuously expanding surfaces at an exponential rate, MHD boundary layer flow of nano-fluids above an exponentially expanding permeable surface and flow above an exponential expansion of sheet through porous medium. This field has drawn attention of many researchers and a significant progress is made in recent past. Convective heat transfer in fluid motion due to exponential stretching of surface is very important as it has influential role in most of the real life phenomenons. Some real world applications are manufacturing industry which includes artificial fibers, metal spinning, wire drawing, liquid films in condensation procedure, crystal growing, cable coating metallurgy, and polymers. Since Crane¹⁷ studied on the flows above expanding sheet, many scholars are involved in analyzing such tasks under distinct features. Magyari and Keller¹⁸ accomplished task by investigating thermo-physical features of viscous fluid above an exponentially stretch sheet. They utilized numerical and analytical techniques to illuminate the exponential alteration in temperature distribution. Sahoo and Poncet¹⁹ exhibited the influences of flow and heat sensation through utilization of partial slip boundary constraint above an exponential expanding sheet. Diverse assumptions on exponentially expanding sheet under distinct physical situations were spotted by Elbashbeshy,²⁰ Partha et al.,²¹ Mabood et al.,²² and El-Aziz.²³

Several problems in boundary-layer flow and heat transfer are non-similar. Non-similar modeling is very powerful tool. If we compare similar models versus non-similar models we will come to a conclusion that only a limited number of boundary layer equations can be tackled through similarity transformations. However a vast majority of boundary layer equations involving exponentially stretching surfaces, Darcy-Forchheimer models, non-Newtonian fluids such as Powell-Eyring fluid, Sisko fluid, Williamson fluid, Prandtl fluid, Burgers fluid, Casson fluid, etc. can only be tackled through appropriate use of non-similarity transformations. Theoretically and practically non-similar boundary layer studies are more important because of their vast applications. Three causes of non-similarities were considered by Sparrow and Yu²⁴ that is spatial variations in the free-stream velocity, transverse curvature, and surface mass transfer (suction/injection of fluid at surface). It may occur due to variation in wall temperature, variable wall thickness, buoyancy force effect, and inclination angle effects. The non-similarity of boundary layer can also arise from more than one factor. In-spite of this fact still very few publications in the field of non-similar flows exist as compared to similar flows. The non-similar PDE’s can be calculated by using the following methods which are implicit finite difference technique, homotopy analysis method.²⁴ and local non-similarity technique via suitable analytical and numerical method.²⁵

In this proposed non-similar model of Oldroyd-B fluid, local non-similarity along second level of truncation via bvp4c is employed to get numerical solutions. In the method of local similarity (first level truncation), PDE’s obtained after transformation have both similar and non-similar terms. Assuming the non-similarity part of the equation very small and neglecting these terms PDE’s become ODE’s. The main drawback of local similarity technique is that it neglects important non-similar expressions. The outcomes calculated by local similarity procedure cannot be precise. To overcome this drawback Minkowycz and Sparrow²⁷ presented a new method called local non-similarity (second level truncation) technique for obtaining solution for non-similar terms which are not zero. Later non-similar flows were studied by Massoudi.²⁸ Cui et al.³¹ exhibited non-similar flow for mixed convective flows of viscous fluids with nanoparticles.

In the method of local non-similarity the physical co-ordinates $(x, y)$ are transformed to ( $ξ, η$ ) dimensionless coordinates. If $ξ$ is zero or a constant, then the non-similar problem is converted into similar problem. In view of above mentioned discussion, the boundary layer flow of an Oldroyd-B fluid along an exponentially stretching surface is addressed. The relaxation and retardation effects are captured through fluid model. Local non-similarity technique is proposed by Sparrow et al.²⁹ and it is extensively used thereafter in many studies.^26,30 Local non-similarity method via bvp4c MATLAB based algorithm is used to compute coupled highly non-linear boundary value problems. The relaxation and retardation time and important physical parameter consequences are examined through graphs.

Mathematical formulation

The external forced flow of an incompressible, Oldroyd-B fluid is considered. Two-dimensional flow is produced by stretching of a sheet at an exponential rate in horizontal direction in Figure 1. The governing convection equations are modeled without incorporating the presence of viscous dissipation and thermal radiation effects.

\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0,

(1)

\begin{matrix} u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + Λ_{1} (u^{2} \frac{\partial^{2} u}{\partial x^{2}} + 2 uv \frac{\partial^{2} u}{\partial x \partial y} + v^{2} \frac{\partial^{2} u}{\partial y^{2}}) \\ = ν Λ_{2} (u \frac{\partial^{3} u}{\partial x \partial y^{2}} - \frac{\partial u}{\partial x} \frac{\partial^{2} u}{\partial y^{2}} - \frac{\partial u}{\partial x} \frac{\partial^{2} v}{\partial y^{2}} + v \frac{\partial^{3} u}{\partial y^{3}}) + ν (\frac{\partial^{2} u}{\partial y^{2}}), \end{matrix}

(2)

u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = α (\frac{\partial^{2} T}{\partial y^{2}}) .

(3)

Figure 1.

Physical configuration.

In this model, $u$ and $v$ are velocity components in direction of $x$ and $y . Λ_{1}$ is relaxation time, $T$ is temperature, and $Λ_{2}$ is retardation time, respectively. $ν$ is kinematic viscosity and $α$ is thermal diffusivity, respectively. Here boundary conditions are

\begin{matrix} At y = 0, u = U_{0} \exp (\frac{x}{l}), v = 0, T = T_{\infty} + T_{0} \exp (\frac{x}{2 l}), \\ As y \to \infty, u \to 0, v \to 0, T \to T_{\infty} . \end{matrix}

(4)

Here consider the non-similar flow by introducing new parameters $ξ (x)$ as non-similarity variable and $η (x, y)$ as pseudo-similarity variable written as

ξ = \exp (\frac{x}{l}), η = y \sqrt{\frac{U_{0}}{2 ν l}} \exp (\frac{x}{2 l}) .

(5)

Now by using chain-rule for connecting the parameters in ( $x, y)$ plane to those in $(ξ, η)$ plane, we get

\begin{matrix} u = U_{0} \exp (\frac{x}{l}) \frac{\partial f}{\partial η} (ξ, η), T = θ (ξ, η) (T_{I} - T_{\infty}) + T_{\infty}, \\ v = - \sqrt{\frac{ν U_{0}}{2 l}} \exp (\frac{x}{2 l}) (f (ξ, η) + η \frac{\partial f}{\partial η} + 2 ξ \frac{\partial f}{\partial ξ}) . \end{matrix}

(6)

Substitution of (6) in (1)–(3) we get the following momentum, energy, and nano-particle volume fraction equations as

\begin{matrix} \frac{\partial^{3} f}{\partial η^{3}} + f \frac{\partial^{2} f}{\partial η^{2}} - 2 {(\frac{\partial f}{\partial η})}^{2} + \\ β_{2} (2 ξ \frac{\partial^{3} f}{\partial η^{3}} \frac{\partial f}{\partial η} - ξ f \frac{\partial^{4} f}{\partial η^{4}} + 3 ξ f {(\frac{\partial^{2} f}{\partial η^{2}})}^{2} - 2 ξ \frac{\partial f on is observed) . A concept of}{\partial ξ} \frac{\partial^{4} f}{\partial η^{4}}) + \\ β_{1} (\begin{matrix} - 4 ξ {(\frac{\partial f}{\partial η})}^{3} + η ξ {(\frac{\partial f}{\partial η})}^{2} \frac{\partial^{2} f}{\partial η^{2}} + 6 ξ f \frac{\partial f}{\partial η} \frac{\partial^{2} f}{\partial η^{2}} - 4 ξ^{3} {(\frac{\partial f}{\partial η})}^{2} \frac{\partial^{3} f}{\partial ξ^{2} \partial η} \\ - 12 ξ^{2} {(\frac{\partial f}{\partial η})}^{2} \frac{\partial^{2} f}{\partial ξ \partial η} - 4 ξ^{2} f \frac{\partial^{3} f}{\partial η^{3}} \frac{\partial f}{\partial ξ} - 4 ξ^{2} \frac{\partial^{3} f}{\partial η^{3}} {(\frac{\partial f}{\partial ξ})}^{2} + 4 ξ^{2} f \frac{\partial^{3} f}{\partial ξ \partial η^{2}} \frac{\partial f}{\partial η} \\ + 12 ξ^{2} \frac{\partial f}{\partial η} \frac{\partial f on is observed) . A concept of}{\partial ξ} \frac{\partial^{2} f}{\partial η^{2}} + 8 ξ^{3} \frac{\partial f}{\partial η} \frac{\partial f}{\partial ξ} \frac{\partial^{3} f}{\partial ξ \partial η^{2}} - ξ f^{2} \frac{\partial^{3} f}{\partial η^{3}} \end{matrix}) \\ = 2 ξ (\frac{\partial f}{\partial η} \frac{\partial^{2} f}{\partial ξ \partial η} - \frac{\partial^{2} f}{\partial η^{2}} \frac{\partial f}{\partial ξ}) . \end{matrix}

(7)

\frac{\partial^{2} θ}{\partial η^{2}} + \Pr (f \frac{\partial θ}{\partial η} - θ \frac{\partial f}{\partial η}) = 2 \Pr ξ (\frac{\partial f}{\partial η} \frac{\partial θ}{\partial ξ} - \frac{\partial θ}{\partial η} \frac{\partial f}{\partial ξ}) .

(8)

\begin{matrix} \frac{\partial f}{\partial η} (ξ, 0) = 1, \frac{\partial f}{\partial η} (ξ, \infty) = 0, θ (ξ, 0) = 1, \\ f (ξ, 0) + 2 ξ \frac{\partial f}{\partial ξ} = 0, θ (ξ, \infty) = 0 . \end{matrix}

(9)

where dimensionless relaxation time constant $β_{1}$ , retardation time constant $β_{2}$ , and the Prandtl number $\Pr$ are defined as

β_{1} = \frac{Λ_{1} U_{o}}{2 L}, β_{2} = \frac{Λ_{2} U_{o}}{2 L}, \Pr = \frac{ν}{α} .

The proportion between convective heat shifting to a conductive heat transfer of fluid bring up to Nusselt number $N u_{x}$ . It offers the heat shifting rate at surface, interpreted as

C_{f} = \frac{2 τ_{I}}{ρ U_{I}^{2}}, Nu = \frac{x q_{I}}{K (T_{\infty} - T_{I})} .

(10)

Dimensionless presentation mentioned below;

{(2 Re)}^{\frac{1}{2}} C_{f} = \sqrt{2} f^{″} (ξ, 0), R e^{- \frac{1}{2}} Nu = - \sqrt{\frac{1}{2}} θ^{'} (ξ, 0) \ln (ξ) .

(11)

To interpret several non-similar boundary layer problems many researchers utilize the local non-similarity technique.^31,32 The local similarity principle is usually assumed for consequences of non-similar boundary layer. With the help of local non-similar technique, we have approximated non-similar PDE’s through ODE’s. A technique for finding locally non-similar boundary layer is discussed and demonstrated in detail by Farooq et al.^33,34

Results and discussion

The prime purpose of this section is to express the outcomes for dimensionless parameters of velocity $f^{'} (η)$ and temperature $θ (η)$ . Numerical solutions have been calculated for equations (7) and (8) subject to boundary conditions (9) by using MATLAB based algorithm bvp4c. Figures 2 and 3 presents the influences of relaxation and retardation time constant $β_{1}, β_{2}$ respectively on the flow field $f^{'} (η)$ . The consequences of $β_{1}$ on velocity function $f^{'}$ are demonstrated in Figure 2. The velocity frame $f^{'} (η)$ and the boundary layer thickness declined for wider $β_{1}$ . The reason behind this certainty that lazy recovery technique is noted for greater relaxation time, which is a source for boundary layer width to raise at delaying rate.

Figure 2.

Graph of $f^{'} (η)$ for distinct $β_{1}$ .

Figure 3.

Graph of $f^{'} (η)$ for distinct $β_{2}$ .

Figure 3 shows the consequences of retardation time $β_{2}$ on velocity $f^{'} (η) .$ Present, by raising $β_{2}$ the enlargement in boundary layer thickness and fluid flow is acquired.

Figure 4 illustrates the distribution of velocity profile $f^{'} (η)$ versus dimensionless streamwise coordinate $ξ$ for the stretching sheet with forced convection. The velocity profile $f^{'} (η)$ is highest at prime edge and its decaying continuously and approaches to zero with increasing $η$ . Increasing points of $ξ$ revealed an expansion in the velocity.

Figure 4.

Graph of $f^{'} (η)$ for distinct $ξ .$

Figure 5 portray the impact of Prandtl number, by raising Prandtl number, temperature profile reduced. It is because of bigger Prandtl number, diffusion of energy will be lesser. Hence, enhancing Prandtl number concluded a powerful reduction in temperature figure of fluid which becomes reason for thinner boundary layer.

Figure 5.

Graph of $θ (η)$ for distinct $\Pr .$

Figures 6 and 7 are plotted for the effects of $β_{1}$ , $β_{2}$ which are relaxation and retardation time constants on the temperature θ $(η)$ . Influence of $β_{1}$ on θ $(η)$ can be clearly viewed in Figure 6, kinetic energy of molecules within substance increases due to rise in temperature. So, velocity of the molecules will be increased and collisions will occur more frequently and hence the relaxation time will be decreased.

Figure 6.

Graph of $θ (η)$ for distinct $β_{1} .$

Figure 7.

Graph of $θ (η)$ for distinct $β_{2} .$

Figure 7 clearly exhibits the influences of $β_{2}$ against θ $(η)$ . It is remarked that an expansion in $β_{2}$ lowered the temperature portrait θ $(η)$ . Thus, it deduced that $β_{1}$ and $β_{2}$ has opposite outcomes on $f^{'} (η),$ θ $(η)$ because of relaxation and retardation times.

Table 1 shows the range table for distinct parameters. Table 2 exhibits the comparison of heat transfer coefficient for various values of $\Pr$ . Table 3 presents the skin friction and local Nusselt number, it also shows the comparison of present outcomes with published work⁴ for the local Nusselt number, it has shown the outcomes for various parameters of surface heat shifting rate. These tabular values reveal that Nusselt number declined by expanding the Prandtl number Pr and retardation time constant $β_{2}$ , while it boost by increasing relaxation time constant $β_{1}$ .

Table 1.

The range of distinct parameters for a stable solution.

$β_{1}$	$β_{2}$	$\Pr$
0.1–20	0.3	3
0.1	0.1–2.0	3
0.3	0.3	1–30

Table 2.

Comparison of heat transfer coefficient, for various values of $\Pr = 1, 2, 3$ when $β_{1} = β_{2} = ξ = 0$ with Mukhopadhyay,³⁵ Bidin and Nazar,³⁶ and Khan et al.³⁷

$\Pr$	Mukhopadhyay³⁵	Khan et al.³⁷	Bidin and Nazar³⁶	Present results
	Shooting method	Shooting method	Keller-Box method	bvp4c
1	0.9547	0.95473	0.9547	0.9547
2	1.4714	1.47151	1.4714	1.4714
3	1.8961	1.8961	1.8961	1.8961

Table 3.

Comparison of local similar and non-similar $Cf, Nu$ (local skin friction, local Nusselt number) for distinct values of $\Pr$ , $β_{1} and β_{2}$ when $ξ = 1.0$ .

Parameters			Present	Present	Hayat et al.⁴	Difference (%)
$β_{1}$	$β_{2}$	$\Pr$	$Cf ({Re}_{x}^{\frac{1}{2}})$	$Nu ({Re}_{x}^{- \frac{1}{2}})$	Hayat et al.⁴	Difference (%)
0.1	0.1	1	−0.3752841074	−1.1505686009	−0.96468	19
0.2			−0.4197797111	−1.1301294995	−0.95470	18
0.3			−0.4565666995	−1.116288424	−0.94533	18
0.4			−0.4878735710	−1.1070892517	−0.93621	18
0.5			−0.5149484768	−1.0946129530	−0.92758	18
0.1	0.2	1	−0.3547247166	−1.6770848486	−0.98989	69
	0.3		−0.2986584919	−1.6882361267	−1.01010	67
	0.4		−0.2723104359	−1.6848666221	−1.02626	64
	0.5		−0.2696032243	−1.6888764297	−1.04050	62
0.1	0.1	1.1	−0.3752841073	−1.2124937391	−1.02560	18
		1.2	−0.3752841080	−1.2715955504	−1.08390	17
		1.5	−0.3752841076	−1.4351676038	−1.24580	15
		1.8	−0.3752841076	−1.5827626912	−1.39240	14

Conclusion

In this research the important outcomes are as follows:

The velocity portrait increases for greater $ξ .$

The influences of $β_{2}$ are against to those of $β_{1}$ on velocity figure $f^{'}$ .

For large values of Prandtl number, the temperature and thermal boundary surface width reduces.

The dissimilarity of $β_{1}$ and $β_{2}$ are qualitatively close on the temperature view.

Nusselt number expanded for greater $β_{1}$ and Pr while bring down by raising $β_{2}$ .

Footnotes

Handling Editor: James Baldwin

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Umer Farooq

References

Sakiadis

. Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. Am Inst Chem Eng J 1961; 7: 26–28.

Crane

LJ.

Flow past a stretching plate. J Appl Math Phys 1970; 21: 645–647.

Chiam

TC.

Stagnation point flow towards a stretching plate. J Phys Soc Jpn 1994; 63: 2443–2444.

Hayat

Imtiaz

Alsaedi

Boundary layer flow of Oldroyd-B fluid by exponentially stretching sheet. Appl Math Mech 2016; 37: 573–582.

Farooq

Munir

, et al. MHD flow of Maxwell fluid with nanomaterials due to an exponentially stretching surface. Sci Rep 2019; 9: 7312.

Tanner

RY.

Note of the Rayleigh problem for a visco-elastic fluid. Z Angew Math Phys 1962; 13: 573–580.

Waters

N D

King

MJ.

Unsteady flow of an elastico-viscous liquid. Rheol Acta 1970; 9: 345–355.

Vieru

Fetecau

Flow of a generalized Oldroyd-B fluid due to a constantly accelerating plate. Appl Math Comput 2008; 201: 834–842.

Fetecau

Hayat

Zierep

, et al. Energetic balance for the Rayleigh-Stokes problem of an Oldroyd-B fluid. Nonlinear Anal Real World Appl 2011; 12: 1–13.

10.

Jamil

Fetecau

Imran

Unsteady helical flows of Oldroyd-B fluids. Commun Nonlinear Sci Numer Simul 2011; 16: 1378–1386.

11.

Hayat

Shehzad

Alsaedi

, et al. Three-dimensional flow of Oldroyd-B fluid over surface with convective boundary conditions. Appl Math Mech 2013; 34: 489–500.

12.

Ramzan

Farooq

Alhothuali

, et al. Three dimensional flow of an Oldroyd-B fluid with Newtonian heating. Int J Numer Methods Heat Fluid Flow 2015; 25: 68–85.

13.

Hayat

Imtiaz

Alsaedi

, et al. On Cattaneo-Christov heat flux in MHD flow of Oldroyd-B fluid with homogeneous-heterogeneous reactions. J Magn Magn Mater 2016; 401: 296–303.

14.

Shehzad

Alsaedi

Hayat

, et al. Thermophoresis particle deposition in mixed convection three-dimensional radiative flow of an Oldroyd-B fluid. J Taiwan Inst Chem Eng 2014; 45: 787–794.

15.

Zheng

Liu

Zhang

Exact solutions for MHD flow of generalized Oldroyd-B fluid due to an infinite accelerating plate. Math Comput Model 2011; 54: 780–788.

16.

Kishan

Kalyani

Chenna Krishna Reddy

MHD boundary layer flow of a nanofluid over an exponentially permeable stretching sheet with radiation and heat Source/Sink. Nano Micro Scale Sci Technol 2016; 4: 44–51.

17.

Crane

LJ.

Flow past a stretching plate. Z Angew Math Phys 1961: 7: 21–28.

18.

Magyari

Keller

Heat and mass transfer in the boundary layers on an exponentially stretching continuous surface. J Phys D Appl Phys 1999; 32: 577–585.

19.

Sahoo

Poncet

Flow and heat transfer of a third grade fluid past an exponentially stretching sheet with partial slip boundary condition. Int J Heat Mass Transf 2011; 54: 5010–5019.

20.

Elbashbeshy

EMA

. Heat transfer over an exponentially stretching continuous surface with suction. Arch Mech 2001; 53: 643–651.

21.

Partha

Murthy

Rajasekhar

GP.

Effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface. Int J Heat Mass Transf 2005; 41: 360–366.

22.

Mabood

Khan

Ismail

AM.

MHD flow over exponential radiating stretching sheet using homotopy analysis method. J King Saud Univ Eng Sci 2017; 29: 68–74.

23.

El-Aziz

MA.

Viscous dissipation effect on mixed convection flow of a micropolar fluid over an exponentially stretching sheet. Can J Phys 2009; 87: 359–368.

24.

Sparrow

HS.

Local non-similarity thermal boundary-layer solutions. J Heat Transf 1971; 93: 328–334.

25.

Farooq

Hayat

Alsaedi

, et al. Series solutions of non-similarity boundary layer flows of nano-fluids over stretching surfaces. Numer Algorithms 2015; 70: 43–59.

26.

Farooq

Ijaz

Khan

, et al. Modeling and non-similar analysis for Darch-Forchheimer-Brinkman model of Cassson fluid in a prous media. Int Commun Heat Mass Transf 2020; 119: 104955.

27.

Minkowycz

Sparrow

EM.

Local nonsimilar solutions for natural convection on a vertical cylinder. J Heat Transf 1974; 96: 178–184.

28.

Massoudi

Local non-similarity solutions for the flow of a non-Newtonian fluid over a wedge. Int J Non Linear Mech 2001; 36: 961–976.

29.

Sparrow

Quack

Boerner

CJ.

Local non-similarity boundary layer solutions. Am Inst Aeronaut Astronaut J 1970; 8: 1936–1942.

30.

Mushtaq

Asghar

Hossain

MA.

Mixed convection flow of second grade fluid along a vertical stretching flat surface with variable surface temperature. Heat Mass Transf 2007; 43: 1049–1061.

31.

Cui

Farooq

Razzaq

, et al. Closure to “Computational analysis for mixed convective flows of viscous fluids with nanoparticles” (Farooq, U, Lu, D C, Ahmed, S, and Ramzan, M. ASME J. Therm. Sci. Eng. Appl. 2019; 11(2), p. 021013). ASME J Thermal Sci Eng Appl 2021; 13: 066001.

32.

Raees

Farooq

Hussain

, et al. Non-similar mixed convection analysis for magnetic flow of second-grade nanofluid over a vertically stretching sheet. Commun Theor Phys 2021; 73: 065801.

33.

Farooq

Razzaq

Khan

, et al. Modeling and numerical computation of non-similar forced convective flow of viscous material towards an exponential surface. Int J Mod Phys B 2021; 15: 2150118.

34.

Farooq

Hussain

Ijaz

, et al. Impact of non-similar modeling on Darcy-Forchheimer-Brinkman model for forced convection of Casson nano-fluid in non-Darcy porous media, Int Commun Heat Mass Transf 2021; 125: 105312.

35.

Mukhopadhyay

MHD boundary layer flow and heat transfer over an exponentially stretching sheet embedded in a thermally stratified medium. Alex Eng J 2013; 52: 259–265.

36.

Bidin

Nazar

Numerical solution of the boundary layer flow over an exponentially stretching sheet with thermal radiation. Eur J Sci Res 2009; 33: 710–717.

37.

Khan

Hayat

Alsaedi

, et al. Melting heat transportation in radiative flow of nanomaterials with irreversibility analysis. Renew Sustain Energy Rev 2021; 140: 110739.