Analytical Investigation of Laminar Viscoelastic Fluid Flow over a Wedge in the Presence of Buoyancy Force Effects

Abstract

An analytical strong method, the homotopy analysis method (HAM), is employed to study the mixed convective heat transfer in an incompressible steady two-dimensional viscoelastic fluid flow over a wedge in the presence of buoyancy effects. The two-dimensional boundary-layer governing partial differential equations (PDEs) are derived by the consideration of Boussinesq approximation. By the use of similarity transformation, we have obtained the ordinary differential nonlinear (ODE) forms of momentum and energy equations. The highly nonlinear forms of momentum and energy equations are solved analytically. The effects of different involved parameters such as viscoelastic parameter, Prandtl number, buoyancy parameter, and the wedge angle parameter, which is related to the exponent 𝑚 of the external velocity, on velocity and temperature distributions are plotted and discussed. An excellent agreement can be seen between the results and the previously published papers for 𝑓^′′(0) and 𝜃^′(0) in some of the tables and figures of the paper for velocity and temperature profiles for various values of viscoelastic parameter and Prandtl number. The effects of buoyancy parameter on the velocity and temperature distributions are completely illustrated in detail.

1. Introduction

Since the non-Newtonian fluid flow and heat transfer represents many important applications such as plastic films and artificial fibers, it is one of the most attractive fields in different aspects of engineering for the last few decades. One of the important studies is related to the convective heat transfer over a surface which can be observed vastly in engineering, agriculture, and petroleum industries [1, 2]. Hiemenz [3] was the first one who started the study of stagnation flow problem and obtained the ODE form of equations for the forced convective problem by introducing a similarity transformation. Dash and Behera [4] investigated laminar free convective viscoelastic fluid flow and heat transfer over an isothermal cylinder. Nazar et al. [5] studied the micropolar fluid flow over a stretching sheet in stagnation flow. Viscoelastic MHD flow and heat transfer over a stretching sheet was investigated by Abel et al. [6] considering viscous and ohmic dissipations. Nadeem and Akbar [7–10] solved different types of fluid flows such as non-Newtonian, Williamson, and tangent hyperbolic fluids in an endoscope analytically, numerically, and exactly. Presenting a numerical algorithm, Ariel [11] studied viscoelastic (second grade) fluid flow near a stagnation point. Finite difference method using Thomas algorithm was employed by Mahapatra and Gupta [12] to survey viscoelastic (Walters’ B liquid) fluid flow. An off-centered stagnation flow over a rotating disc was solved by Erfani et al. [13] by the modified differential transform method (MDTM). Vogel’s model of viscosity on the peristaltic flow of Jeffrey fluid was considered by Akbar et al. [14] in analytical and numerical forms. Ishak et al. [15] presented the results of stagnation-point flow in a permeable sheet numerically via an implicit finite difference scheme known as the Keller-Box method. Rashidi et al. [16] also analyzed the stagnation-point flow in a porous medium by DTM. Kasim et al. [17] considered heat generation in a free convective viscoelastic fluid flow over a horizontal circular cylinder at the lower stagnation point of the cylinder. Aman et al. [18] considered the slip condition in mixed convective boundary-layer flow numerically. The unsteady three-dimensional stagnation-point flow of a viscoelastic fluid was studied by Seshadri [19]. Turkyilmazoglu [20] presented exact solution for mixed convection over a stretching surface. Bhattacharyya and Layek [21] presented numerically the results of suction/blowing and thermal radiation effects over a porous shrinking sheet. Bachok et al. [22] and Layek et al. [23] employed Runge-Kutta-Fehlberg and fourth order classical Runge-Kutta methods in the study of stretching/shrinking sheets, respectively. Turkyilmazoglu [24–26] presented multiple solutions in viscoelastic MHD fluid flow and heat and mass transfer over stretching surfaces in slip condition.

Unlike the number of studies related to stretching sheets, the articles on the convective flow over a wedge are very limited. A nonisothermal wedge in the presence of a heat source/sink was the base of investigation solved by Chamkha et al. [27] in the presence of thermal radiation effects by the finite-difference method. Kandasamy et al. [28] considered the effects of chemical reaction, variable viscosity, and thermophoresis over a porous wedge. Hossain et al. [29] presented the numerical results for unsteady flow over a symmetric wedge with variable surface temperature.

Homotopy analysis method (HAM) is one of the most well-known methods to solve highly nonlinear problems. The first one who employed HAM, Liao, offered a general analytical method for nonlinear problems [30, 31]. Rashidi et al. [32] used this method for mixed convective boundary-layer flow of a micropolar fluid towards a heated shrinking sheet. The stagnation-point flow of a nanofluid over a stretching sheet was investigated by Mustafa et al. [33] via HAM. Dinarvand et al. [34] used HAM to solve the unsteady laminar flow near stagnation point of a rotating and translating sphere. Hayat et al. [35] considered MHD flow of an upper-convected (UCM) fluid over a stretching surface by means of HAM. Currently, HAM has been used vastly by many researchers in different practical aspects of engineering and nonlinear problems [36–42].

In the present paper, we try to find the analytical solution for two-dimensional incompressible viscoelastic fluid flow over a wedge. Analytical solutions for the velocity and the temperature distributions are obtained using a powerful technique, namely, the HAM. The graphs are plotted and discussed for the variations of different involved parameters.

2. Flow Analysis

Consider a steady and laminar incompressible two-dimensional mixed convective heat transfer of a viscoelastic fluid flow over a wedge in the presence of buoyancy force effects. It is assumed that the external velocity is in the form of 𝑢_𝑒(𝑥)=𝑎𝑥^𝑚 where 𝑎 and 𝑚 are constants. The Cartesian coordinate system is supposed to help the solution in which the 𝑥-axis and the 𝑦-axis are along with and perpendicular to the wedge surface, respectively. The schematic diagram of the problem is plotted in Figure 1 in order to simplify the problem’s realization. Considering the Boussinesq, the boundary-layer approximations, and the above assumptions, we can derive the boundary-layer equations in the following format: 𝜕𝑢+𝜕𝑥𝜕𝑣𝜕𝑦=0,𝑢𝜕𝑢𝜕𝑥+𝑣𝜕𝑢𝜕𝑦=𝑢_𝑒𝑑𝑢_𝑒𝜕𝑑𝑥+𝜐²𝑢𝜕𝑦²+𝑘₀𝑢𝜕³𝑢𝜕𝑥𝜕𝑦²+𝜕𝑢𝜕𝜕𝑥²𝑢𝜕𝑦²+𝜕𝑢𝜕𝜕𝑦²𝑣𝜕𝑦²𝜕+𝑣³𝑢𝜕𝑦³+𝑔𝛽_𝑇𝑇−𝑇_∞ΩSin2,𝑢𝜕𝑇𝜕𝑥+𝑣𝜕𝑇=𝑘𝜕𝑦𝜌𝑐_𝑃𝜕²𝑇𝜕𝑦², (1)(2)(3) where 𝑢 and 𝑣 are velocity components in the directions of 𝑥 and 𝑦 (as shown in Figure 1). 𝜐 is the kinematic viscosity, 𝑘₀ is the viscoelasticity parameter, 𝑔 is the acceleration due to the gravity, 𝛽_𝑇 is the coefficient of thermal expansion, 𝑘 is the thermal conductivity, 𝜌 is the fluid density, 𝑐_𝑃 is the specific heat at constant pressure, and 𝑇 is the fluid temperature. Ω=𝜋𝛽 is the total angle of the wedge (for 𝛽=0 and 𝛽=1, the wedge takes the form of horizontal and vertical plate, resp.) and 𝛽 is the wedge angle parameter which is defined in the form of 𝛽=2𝑚/(𝑚+1).

Figure 1 

The schematic diagram of the wedge.

The corresponding boundary conditions are as follows: 𝑢=0,𝑣=0,𝑇=𝑇_𝑤(𝑥),at𝑦=0,𝑢⟶𝑢_𝑒(𝑥),𝜕𝑢𝜕𝑦⟶0,𝑇⟶𝑇_∞as𝑦⟶∞. (4)

It is assumed that 𝑇_𝑤(𝑥)=𝑇_∞+𝑏𝑥^𝑚 where 𝑏 is constant value. After introducing the stream function 𝜓 and similarity variable 𝜂 and satisfying the continuity equation, we can derive the momentum and energy equations into the ordinary differential equations (ODEs) format as one can see in the following equations: 𝜂=𝑢_𝑒(𝑥)𝜐𝑥𝑦,𝜓=𝜐𝑥𝑢_𝑒(𝑥)𝑓𝜂,𝜃𝜂=𝑇−𝑇_∞𝑇_𝑤−𝑇_∞,𝑚𝑓^′2−1−𝑚+12𝑓𝑓^′′−𝑓^′′′−𝑘₁(3𝑚−1)𝑓^′𝑓^′′′−(3𝑚−1)2𝑓^2′′−(𝑚+1)2𝑓𝑓⁽⁴⁾−𝜆_𝑇𝑚Sin𝜋𝑚+1𝜃=0,𝜃^′′+Pr𝑚+12𝑓𝜃^′−𝑚𝑓^′𝜃=0, (5)(6)(7) where superscript ′ denotes the derivative with respect to 𝜂, 𝑘₁=𝑘₀𝑎𝑥^𝑚−1/𝜐 is the viscoelastic parameter (when 𝑚=1, the viscoelastic parameter takes the form of 𝑘₁=𝑘₀𝑎/𝜐 similar to the viscoelastic parameter obtained by Hayat et al. [44]), and 𝜆_𝑇=𝑔𝛽_𝑇(𝑇_𝑤−𝑇_∞)𝑥/(𝑎²𝑥^(2𝑚−1)𝑥)=Gr_𝑥/Re²_𝑥 is the buoyancy parameter, where Gr_𝑥=𝑔𝛽_𝑇(𝑇_𝑤−𝑇_∞)𝑥³/𝜐² is the Grashof number, Re_𝑥=𝑢_𝑒𝑥/𝜐=𝑎𝑥^𝑚+1/𝜐 is the Reynolds number, and Pr=𝜇𝑐_𝑃/𝑘 is the Prandtl number. The corresponding boundary conditions are as follows: 𝑓𝜂=0,𝑓^′𝜂𝜂𝑓=0,𝜃=1,at𝜂=0,_′𝜂=1,𝑓^′′𝜂=0,𝜃𝜂=0,as𝜂⟶∞. (8)

3. HAM Solution

In order to satisfy the boundary conditions, the initial approximations must be chosen in the appropriate form as follows: 𝑓₀𝜂=𝜂+𝑒^−𝜂,𝜃−1₀𝜂=𝑒^−𝜂. (9) The linear operators ℒ_𝑓(𝑓), ℒ_𝜃(𝜃) are defined as ℒ_𝑓𝑓=𝜕⁴𝑓𝜕𝜂⁴+𝜕³𝑓𝜕𝜂³,ℒ_𝜃(𝜃)=𝜕²𝜃𝜕𝜂²+𝜕𝜃,𝜕𝜂 (10) with the following properties: ℒ_𝑓𝑐₁+𝑐₂𝜂+𝑐₃𝜂²+𝑐₄𝑒^−𝜂ℒ=0,_𝜃𝑐₅+𝑐₆𝑒^−𝜂=0, (11) where 𝑐₁−𝑐₆ are arbitrary constants. The nonlinear operators are 𝒩_𝑓𝑓,𝜃𝜕𝜂;𝑞𝜂;𝑞=𝑚𝑓(𝜂;𝑞)𝜕𝜂²−−1𝑚+12𝑓𝜕𝜂;𝑞²𝑓𝜂;𝑞𝜕𝜂²−𝜕³𝑓𝜂;𝑞𝜕𝜂³−𝑘₁()𝜕𝑓3𝑚−1𝜂;𝑞𝜕𝜕𝜂³𝑓𝜂;𝑞𝜕𝜂³−𝑚+12𝑓𝜕𝜂;𝑞⁴𝑓𝜂;𝑞𝜕𝜂⁴−3𝑚−12𝜕²𝑓(𝜂;𝑞)𝜕𝜂²²−𝜆_𝑇𝑚Sin𝜋𝜃,𝒩𝑚+1𝜂;𝑞_𝜃𝑓,𝜃=𝜕𝜂;𝑞𝜂;𝑞₂𝜃𝜂;𝑞𝜕𝜂²+Pr𝑚+12𝑓𝜕𝜃𝜂;𝑞𝜂;𝑞𝜕𝑓𝜕𝜂−𝑚𝜂;𝑞𝜃.𝜕𝜂𝜂;𝑞 (12) The auxiliary functions are introduced as ℋ_𝑓𝜂=ℋ_𝜃𝜂=𝑒^−𝜂. (13) Now, the 𝑖th order deformation equations (14) are solved by the symbolic software MATHEMATICA: ℒ_𝑓𝑓_𝑖𝜂−𝜒_𝑖𝑓_𝑖−1𝜂=ℏℋ_𝑓𝜂𝑅_𝑓,𝑖𝜂,ℒ_𝜃𝜃_𝑖𝜂−𝜒_𝑖𝜃_𝑖−1𝜂=ℏℋ_𝜃𝜂𝑅_𝜃,𝑖𝜂, (14) where ℏ is the auxiliary nonzero parameter: 𝑅_𝑓,𝑖𝜂=^𝑖−1_𝑗=0𝑚𝜕𝑓_𝑗𝜂𝜕𝜂𝜕𝑓_{𝑖−1−𝑗}𝜂−𝜕𝜂𝑚+12𝑓_𝑗𝜂𝜕²𝑓_{𝑖−1−𝑗}𝜂𝜕𝜂²−𝜕³𝑓_𝑖−1𝜂𝜕𝜂³−𝑚1−𝜒_𝑖−^𝑖−1_𝑗=0𝑘₁(3𝑚−1)𝜕𝑓_𝑗𝜂𝜕𝜕𝜂³𝑓_{𝑖−1−𝑗}𝜂𝜕𝜂³−𝑚+12𝑓_𝑗𝜂𝜕⁴𝑓_{𝑖−1−𝑗}𝜂𝜕𝜂⁴−3𝑚−12𝜕²𝑓_𝑗𝜂𝜕𝜂²𝜕²𝑓_{𝑖−1−𝑗}𝜂𝜕𝜂²−𝜆_𝑇𝑚Sin𝜋𝜃𝑚+1^𝑖−1𝜂,𝑅_𝜃,𝑖𝜂=𝜕²𝜃_𝑖−1𝜂𝜕𝜂²+Pr^𝑖−1_𝑗=0𝑚+12𝑓_𝑗𝜂𝜕𝜃_{𝑖−1−𝑗}𝜂𝜕𝜂−𝑚𝜃_𝑗𝜂𝜕𝑓_{𝑖−1−𝑗}𝜂,𝜒𝜕𝜂_𝑖=0,𝑖≤1,1,𝑖>1. (15) For more information about the HAM solution, please see [30, 31].

In Figure 2  ℏ-curves are plotted for an especial case, obtained via 20th order of HAM solution. Choosing a proper value of auxiliary parameter from the valid region in straight line is very important to control the convergence of the approximation series in the so-called ℏ-curve. To check the accuracy of the method, we have presented Tables 1–4 and Figures 3 and 4. An excellent agreement can be found between our results and the published papers.

Table 1 

Comparison of 𝑓^′′(0) for various values of 𝑘₁ and Pr when 𝑚=1.0 and 𝜆_𝑇=0.2.

Table 2 

Comparison of −𝜃^′(0) for various values of 𝑘₁ and Pr when 𝑚=1.0 and 𝜆_𝑇=0.2.

Table 3 

Comparison of 𝑓^′′(0) for various values of 𝑘₁ and Pr when 𝑚=1.0 and 𝜆_𝑇=1.0.

Table 4 

Comparison of −𝜃^′(0) for various values of 𝑘₁ and Pr when 𝑚=1.0 and 𝜆_𝑇=1.0.

Figure 2 

The ℏ-curves of 𝑓′′(0) and 𝜃′(0) obtained by the 20th order approximation of the HAM solution when 𝑘1=1.0, 𝜆𝑇=0.2, Pr=1.0, and 𝑚=0.5.

Figure 3 

Verification of the results 𝑓′(𝜂) with the previously published paper by Li et al. when 𝑘1=0.3, 𝑚=1.0, and Pr=1.0.

Figure 4 

Verification of the results 𝜃(𝜂) with the previously published paper by Li et al. when 𝑘1=0.3, 𝑚=1.0, and Pr=1.0.

In order to choose the optimal value of auxiliary parameter ℏ, we have presented the average residual error as (see [40, 45–47], for more details) Δ_𝑓,𝑚=1𝐾^𝐾_𝑖=0Res_𝑓^𝑚_𝑗=0𝑓_𝑗(𝑖Δ𝑥)²,Δ_𝜃,𝑚=1𝐾^𝐾_𝑖=0Res_𝜃^𝑚_𝑗=0𝜃_𝑗(𝑖Δ𝑥)², (16) where Δ𝑥=10/𝐾 and 𝐾=20. For the given order of approximation 𝑚, the optimal value of ℏ is given by the minimum values of Δ_𝑓,𝑚 and Δ_𝜃,𝑚 corresponding to nonlinear algebraic equations: 𝑑Δ_𝑓,𝑚𝑑ℏ=0,𝑑Δ_𝜃,𝑚𝑑ℏ=0. (17)

For example, in order to find the optimal values of ℏ, the residual error for the HAM 20th order of solution is presented in Figure 5.

Figure 5 

The residual errors for temperature profile when 𝑘1=0.8, 𝜆𝑇=0.5, 𝑚=0.5, and Pr=5.0.

4. Results and Discussion

In this paper the mixed convection of a steady and incompressible two-dimensional viscoelastic fluid flow over a wedge surface is studied. To the best of authors’ knowledge, the current paper is the first paper which presents a similarity solution and includes the effects of buoyancy parameter in mixed convective flow over a horizontal wedge surrounded by the viscoelastic fluid. In this paper, the effects of involved parameters are taken into account and the velocity and temperature distributions are discussed and plotted analytically. To gain a vast understanding of the results, graphical representation is considered for different involved parameters. In Figures 6–13 the effects of different parameters such as viscoelastic parameter, Prandtl number, buoyancy parameter, and the wedge angle parameter are illustrated in detail. Figure 6 denotes the influence of 𝑘₁ on the velocity distribution clearly. Viscoelasticity parameter produces tensile stress, so the boundary-layer thickness decreases and makes it contract, transversely, and hence velocity decreases. As we anticipate and see in Figure 7, the increase in viscoelastic parameter causes the temperature to increase. The increase in 𝑘₁ leads to increase in wall temperature gradient and the nondimensional temperature distribution is enhanced.

Figure 6 

The effect of 𝑘1 on velocity profile when 𝜆𝑇=0.4, Pr=1.0, and 𝑚=0.5.

Figure 7 

The effect of 𝑘1 on temperature profile when 𝜆𝑇=0.4, Pr=1.0, and 𝑚=0.5.

Figure 8 

The effect of 𝜆𝑇 on velocity profile when 𝑘1=0.6, Pr=1.0, and 𝑚=0.5.

Figure 9 

The effect of 𝜆𝑇 on temperature profile when 𝑘1=0.6, Pr=1.0, and 𝑚=0.5.

Figure 10 

The effect of Pr on velocity profile when 𝑘1=0.8, 𝜆𝑇=0.5, and 𝑚=0.5.

Figure 11 

The effect of Pr on temperature profile when 𝑘1=0.8, 𝜆𝑇=0.5, and 𝑚=0.5.

Figure 12 

The effect of 𝑚 on velocity profile when 𝑘1=1.0, 𝜆𝑇=0.2, and Pr=1.0.

Figure 13 

The effect of 𝑚 on temperature profile when 𝑘1=1.0, 𝜆𝑇=0.2, and Pr=1.0.

The influence of 𝜆_𝑇 on the velocity and temperature profiles is plotted in Figures 8 and 9. The only term in momentum equation which causes the energy and momentum equations to be coupled is the buoyancy parameter. When 𝜆_𝑇 increases, the effect of the temperature field on the velocity field will be invigorated and consequently both the velocity and the boundary-layer thickness augment, as shown in Figure 8. This parameter indicates the effect of free convection and buoyancy parameter in equation systems and so reinforces the influence of thermal variations on the velocity component. The buoyancy parameter increases the effect of convection on velocity control. It is worth mentioning that the buoyancy force and favorable pressure gradient have the same impact on velocity in the boundary-layer region and accelerate the fluid. Figure 9 illustrates the effect of 𝜆_𝑇 on temperature distribution. It is clearly depicted that the thermal boundary-layer thickness decreases with the increase in buoyancy parameter because of increase in heat transfer rate. If we want to illustrate the effect of buoyancy parameter in detail, we should go through this parameter. Based on its definition, the thermal buoyancy parameter is the ratio of buoyancy force to viscous force, so the increase in its value suggests a notable increase in the flow velocity. The Grashof number accelerates the fluid, so the velocity and the boundary-layer thickness increase with the increase in 𝜆_𝑇, as shown in Figure 8 due to this fact that the equations are coupled together only by the buoyancy parameter. In fact the buoyancy force acts like a favorable pressure gradient and accelerates the fluid, so the velocity and the boundary-layer thickness increase with the increase in Grashof number and more production occurs. The buoyancy force leads to the increase in temperature gradient and heat transfer rate and the temperature decreases (Figure 9). When Pr increases, the fluid heat capacity is enhanced and the impact of the thermal expansion decreases, so the boundary-layer thickness and velocity distribution decrease (Figure 10). In Figure 11, the effect of Prandtl number on temperature distribution is plotted. Due to the definition of Prandtl number, the ratio of momentum diffusion to thermal diffusion, the thermal diffusion decreases and the fluid thermal capacity increases with the increase in Pr, so the thermal boundary layer becomes thinner and temperature decreases. In Figure 10, the decreasing behavior of velocity component with Pr is not very vivid when Prandtl number varies between 3 and 7 and can be neglected, but the variation in temperature profiles is very distinct and can be easily understood (Figure 11). The effect of parameter 𝑚 is depicted in Figures 12 and 13. Since the increase in 𝑚 causes β to increase, this parameter represents the influence of wedge angle parameter. As 𝑚 increases, boundary-layer thicknesses, velocity, and temperature distributions decrease. The variation of velocity distribution with the wedge angle parameter is the same as its behavior with Pr and is not vivid, either.

In Tables 1–4 the values of 𝑓^′′(0) and −𝜃^′(0) obtained by the 20th order of HAM solution are compared to other different authors’ published results. Through these tables, an excellent agreement can be seen between the present results and the previously published papers’ results. In Tables 1 and 3 the values of 𝑓^′′(0) decrease with the increase in Prandtl number (assisting flow) and the viscoelastic parameter (both assisting and opposing flow regions). Analytical values of wall temperature gradient −𝜃^′(0) for different values of Pr and 𝑘₁ are computed and presented in Tables 2 and 4. Analysis of the tabular data shows that the values of viscoelastic parameter reduce the local heat transfer, while the Prandtl number has the reverse effect. The last term of (2) refers to the buoyancy force. The positive case represents buoyancy assisting and negative case is for the buoyancy opposing flow regions. Buoyancy force assists the upper half of the flow field and opposes the lower half of the flow field. Practically in “+” case the external velocity is opposite to the gravity acceleration and in “−” case the external velocity and the acceleration due to gravity have the same direction (for more details, please see [48]).

In this paper, only assisting flow region, the upper half of the flow field, has been taken into account, but in tables, the results of the both cases have been presented for verification purpose.

5. Conclusion

In the present paper, a steady incompressible viscoelastic fluid flow over a wedge in the presence of buoyancy force effects has been studied analytically by HAM. This analytical solution shows excellent agreement with the data available in the literature (Tables 1–4 and Figures 3 and 4). As one can easily understand, the values of 𝑓^′′(0) decrease with the increase in Prandtl number in assisting flow and vice versa for opposing flow region (Tables 1 and 3). These tables clearly show that the skin friction coefficient decreases with the increase in the value of the viscoelastic parameter for both assisting and opposing flows. Analytical values of wall temperature gradient −𝜃^′(0) for different values of Prandtl number Pr and viscoelastic parameter 𝑘₁ are compared and presented in Tables 2 and 4. Analysis of the tabular data shows that the viscoelastic parameter reduces the Nusselt number. On the contrary, the effect of Prandtl number is to augment the rate of heat transfer. The dimensionless velocity profiles decrease with the increase in 𝑘₁ and the temperature distribution shows increasing behavior. The effect of increasing the buoyancy parameter is to reduce the thermal boundary-layer thicknesses, but the opposite behavior can be seen for velocity component. The impact of Prandtl number and wedge angle parameter on velocity and temperature profiles is the same. In both cases the decreasing behavior of thermal boundary layer is notable, but the effects of Pr and 𝑚 on velocity distribution cannot be distinguished distinctly.

Nomenclature

𝑎,𝑏:	Constant values
𝑐𝑖:	Arbitrary constant
𝑐𝑃:	Specific heat at constant pressure
𝑔:	Acceleration due to gravity
Gr𝑥:	Grashof number
𝑘:	Thermal conductivity
𝑘1:	Viscoelastic parameter
𝑚:	Constant value
Pr:	Prandtl number
Re𝑥:	Reynolds number
𝑇:	Temperature
𝑢,𝑣:	Velocity components along and perpendicular to the wedge
𝑥,𝑦:	Distance along and perpendicular to the wedge, respectively.

Greek Letters

𝛽:	Wedge angle parameter
𝛽𝑇:	Coefficient of thermal expansion
Ω:	Total angle of the wedge
ℏ:	Auxiliary nonzero parameter
ℋ:	Auxiliary function
ℒ:	Auxiliary linear operator
𝒩:	Nonlinear operator
𝜂:	Similarity variable
𝜃:	Dimensionless fluid temperature
𝜌:	Density
𝜆𝑇:	Buoyancy parameter
𝜐:	Fluid kinematic viscosity
𝜓:	Stream function.

Subscripts

𝑒:	External condition
𝑚:	Mean condition
𝑤:	Wall condition
∞:	Infinity condition.

Superscript

′:	Differentiation with respect to 𝜂.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors acknowledge the support from the National Research Foundation of South Africa.