Study of the Carreau fluid over a linearly curved stretching sheet: A numerical approach

Abstract

Carreau fluid flow over a linearly curved stretching surface has significant applications in manufacturing industries. This study focuses on the impact of various parameters on the velocity, temperature, pressure skin friction, and heat transfer coefficient of a 2D Carreau flow. Both shear thinning (pseudoplastic) and shear thickening (dilatant) behaviors are considered. The parameters include suction, non-dimensional radius of curvature, Weissenberg number, magnetic parameter, power-law index, Eckert number (viscous dissipation parameter), and Prandtl number. The effects of viscous dissipation, magnetic field, convective boundary conditions, and suction are also taken into account. The study reveals that the flow velocity decreases with increasing magnetic parameter, curvature factor, and suction, but increases with the Weissenberg number and power law index. The heat transfer rate is reduced by the curvature factor, Weissenberg number, suction, power-law index, and Prandtl number, while the magnetic parameter has the opposite effect. The boundary value problem is simplified through appropriate similarity transformation and solved numerically using the shooting method and Matlab’s built-in function bvp4c. The results obtained from these two methods show good agreement. The findings also align well with existing literature in some limiting cases, confirming the validity of the study.

Introduction

Heat transfer flows due to moving boundaries have many important applications in industries. In manufacturing processes the stretching of sheets has significant effects on the quality of finished products. The final products are significantly influenced by stretching and heating/cooling; depending upon the skin friction or frictional drag and the rate of heat transfer. Under various stretching velocities, different real processes take place, viz; extrusion of polymer sheet, manufacturing of plastic films, continuous stretching, hot rolling, metal spinning, rubber sheet, or thin film and many more.

The boundary layer behavior prompted by the motion of a moving surface has taken attraction of researchers because of its applications in the manufacturing process, such as the flow generated in fiber spinning, glass blowing, hot rolling, paper production, and fibers production. Some of the famous applications of stretching process have been given by Fisher ¹ and Altan et al. ² Sakiadis ³ discussed the effects of constant velocity on Newtonian fluid caused by a moving plate. Erickson et al. ⁴ analyzed the impact of heat mass transfer past a moving plate. Tsou et al. ⁵ reported the influence of constant velocity over a stretching surface in the boundary layer together with temperature. Crane ⁶ carried forward the work of Sakiadis and discussed that in polymer industries it is sometimes essential to consider a stretching plastic sheet. The boundary layer analysis due to a stretchable porous plate of the heat and mass transfer is examined by Gupta and Gupta. ⁷ Recent studies on the curved surface, exponentially stretching curved surface and irregular surfaces emphasizes on the importance of its applications in biological and engineering industries.^8–12

Carreau fluids are a type of generalized Newtonian fluids whose model was put forth by Carreau in 1698. Power law models are enable to predict the viscosity for very high and low shear rate. Carreau fluid is ideally used to model the phenomena of non-Newtonian fluids like plasma, blood flow, and visco-elastics. This model was put forth by Carreau in 1698. Power law model are enable to predict the viscosity for very high and low shear rate. The Carreau Model is

$$\mathit{\tau }=[{\mu }_{\infty }+({\mu }_o-{\mu }_{\infty }){(1+{(\mathrm{\Gamma }\stackrel{·}{\mathit{\gamma }})}^2)}^{\frac{m-1}{2}}]{\stackrel{·}{\mathit{\gamma }}}_{\mathrm{ij}},$$

(1)

In the above equation, $\mathit{\tau }$ and $\stackrel{·}{\mathit{\gamma }}$ represents the extra stress tensor and generalized shear rate, ${\mu }_o$ is the zero shear rate viscosity, ${\mu }_{\infty }$ is the infinite shear rate viscosity, $\mathrm{\Gamma }$ is the material time, $n$ is the power law index. In terms of second invariant tensor $\Pi$ , the generalized shear rate can be expressed as

$$\stackrel{·}{\mathit{\gamma }}=\sqrt{\frac{1}{2}\Sigma \Sigma {\stackrel{·}{\mathit{\gamma }}}_{\mathrm{ij}}{\stackrel{·}{\mathit{\gamma }}}_{\mathrm{ij}}}=\sqrt{\frac{1}{2}\mathbf{\Pi }},$$

(2)

where $\Pi =\mathrm{trace}(\nabla U+{(\nabla U)}^T)$ , Here we consider the case for which ${\mu }_{\infty }=0$ . In the following form the extra stress tensor is rewritten as

$$\mathit{\tau }=\left[{\mu }_o{(1+{(\mathrm{\Gamma }\stackrel{·}{\mathit{\gamma }})}^2)}^{\frac{m-1}{2}}\right]{\stackrel{·}{\mathit{\gamma }}}_{\mathrm{ij}}$$

(3)

Hayat et al. ¹³ investigated the flow of Carreau fluid caused by linear stretching under the effects of Newtonian heating and chemical reactions. Effects of physical parameters on velocity, concentration, and temperature were considered. Khan et al. ¹⁴ implemented shooting scheme to study the behavior of Carreau nanofluids due to paraboloid upper horizontal surface. Carreau fluid over a linearly curved stretching surface is an important model in many engineering and industrial applications. After these initials, a lot of research has been done to compensate the applications of Carreau fluids in industries.^15–18

In numerical analysis, the word shooting for shooting method is originated artillery which describes its working in such a way that it places cannon at initial position and by continuously changing the angle it targets until it reaches and satisfies boundary conditions. Initial value problem is obtain by reducing boundary value problem and by applying appropriate similarity transformation, non-linear ordinary differential equations are obtain from PDE’s of differential equations. The equations are then solved by shooting method coupled with iv order Runge-Kutta integration scheme. This scheme is preferable to solve fluid dynamics problems because it takes advantage of speed and adaptivity over other numerical schemes. ¹⁹

The objectives of this study is concerned with the analysis of boundary driven flow of Carreau fluid over a curved stretching surface. Carreau fluid of both shear thinning (pseudoplastic) and shear thickening (dilatant) nature have been considered. The effects of viscous dissipation, suction and magnetic field are taken into account. Governing equations are reduced to self-similar equations via appropriate similarity transformation. The reduced system of equations are solved numerically through the implementation of shooting method and Matlab built in function bvp4c. The pertinent parameters in this communication are the suctions, non-dimensional radius of curvature, Weissenberg number, magnetic parameter, power-law index, Eckert number (viscous dissipation parameter), and Prandtl number.

Problem formulation

Carreau fluid flow over a curved surface is investigated with incompressible steady property. The surface is curved with radius $R$ and curvilinear natural coordinates are $s$ and $r$ . $R$ describes the distance from origin to the stretched surface. It defines the shape of the surface, for larger values of $R$ the surface goes from curved to flat. Linear stretching velocity is responsible for fluid flow that is $U_w=\mathrm{cs}$ along $s-$ direction. Therefore, the stream forms boundary layer regime through $r-$ direction and temperature distribution of the sheet is assumed $T_w=T_{o}s$ in which $c$ and $T_o$ show reference velocity and temperature. The geometry of flow is shown in Figure 1

Figure 1.

Flow geometry.

The governing equations for the flow are given in vector form as follows see Sajid et al. ²⁰ :

Continuity equation

$$\frac{\partial \rho }{\partial t}+\nabla ·(\rho \mathbf{U})=0,$$

(4)

Momentum equation

$$\rho \frac{D\mathbf{U}}{\mathrm{Dt}}=-\nabla p+\nabla ·\mathit{\tau }+\rho \mathbf{F},$$

(5)

Energy equation

$$\rho c_p\frac{\mathrm{DT}}{\mathrm{Dt}}=\nabla ·(k_o\nabla T)+\mathit{\varphi },$$

(6)

where $\frac{D}{\mathrm{Dt}}$ denotes the material derivative, $\mathbf{U}=(u,v)$ and $\mathbf{F}$ are the velocity vector and the body force, respectively. $p$ , $c_p$ , and $k_o$ are the pressure, specific heat and thermal conductivity. The function $\mathit{\varphi }$ is the viscous dissipation. In the presence of viscous dissipation the governing equations takes the following form:

Mass conservation equation

$$R{u}_s+[(r+R)v{]}_r=0,$$

(7)

Momentum equation

$$\frac{u^2}{r+R}=\frac{1}{\rho }{p}_r,$$

(8)

$$\begin{array}{c}v{u}_r+\frac{\mathrm{Ru}}{r+R}{u}_s+\frac{\mathrm{uv}}{r+R}=-\frac{1}{\rho }\frac{R}{r+R}{p}_s+\nu {[1+{\mathrm{\Gamma }}^2{(u_r-\frac{u}{r+R})}^2]}^{\frac{m-1}{2}}\hfill \\ (u_{\mathrm{rr}}+\frac{1}{r+R}{u}_r-\frac{u}{{(r+R)}^2})+{\mathrm{\Gamma }}^2\nu (m-1){[1+{\mathrm{\Gamma }}^2{(u_r-\frac{u}{r+R})}^2]}^{\frac{m-3}{2}}\hfill \\ {(u_r-\frac{u}{r+R})}^2(u_{\mathrm{rr}}-\frac{1}{r+R}{u}_r+\frac{u}{{(r+R)}^2})-\frac{\sigma B_o^2}{\rho }\frac{R^{2}u}{{(r+R)}^2},\hfill \end{array}$$

(9)

Energy equation

$$\begin{array}{c}\rho c_p(v{T}_r+\frac{\mathrm{Ru}}{r+R}{T}_s)=k_o(T_{\mathrm{rr}}+\frac{1}{r+R}{T}_r)\hfill \\ +\mu {[1+{\mathrm{\Gamma }}^2{(u_r-\frac{u}{r+R})}^2]}^{\frac{m-1}{2}}{(u_r-\frac{u}{r+R})}^2,\hfill \end{array}$$

(10)

boundary conditions for the above defined problems are (see Rosca and Pop ²¹ ):

$$u=u_w=\mathrm{cs},v=-v_w,T=T_w\mathrm{at}r=0,,$$

(11)

$$u\to 0,\frac{\partial u}{\partial r}\to 0,T\to T_{\infty }\mathrm{as}r\to \infty .$$

(12)

To formulate the system dimensionless, we utilized the following similarity variables

$$u=\mathrm{csf}\prime (\alpha ),v=-\frac{R}{r+R}\sqrt{c\nu }f(\alpha ),$$

(13)

$$\alpha =\sqrt{\frac{c}{\nu }}r,p=\rho c^2{s}^{2}P(\alpha ),$$

(14)

$$\theta =\frac{T-T_{\infty }}{T_w-T_{\infty }}.$$

(15)

Governing equations 7–10 with boundary conditions 11–12 are renewed into non-dimensional form as follows:

$$P\prime =\frac{{f\prime }^2}{\alpha +\delta },$$

(16)

$$\begin{array}{c}\frac{2\delta }{\alpha +\delta }P={[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-1}{2}}(f″\prime +\frac{f″}{\alpha +\delta }-\frac{f\prime }{{(\alpha +\delta )}^2})\hfill \\ +(m-1)\mathrm{We}{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-3}{2}}(f″\prime -\frac{f″}{\alpha +\delta }+\frac{f\prime }{{(\alpha +\delta )}^2})\hfill \\ {(f″-\frac{f\prime }{\alpha +\delta })}^2+\frac{\delta }{\alpha +\delta }\mathrm{ff}″+\frac{\delta }{{(\alpha +\delta )}^2}\mathrm{ff}\prime -\frac{\delta }{\alpha +\delta }{f\prime }^2\hfill \\ -\frac{M{\delta }^2}{{(\alpha +\delta )}^2}f\prime ,\hfill \end{array}$$

(17)

$$\begin{array}{c}\theta ″+\frac{\theta \prime }{\alpha +\delta }+\frac{\delta \mathrm{Pr}}{\alpha +\delta }(f\theta \prime -f\prime \theta )+\mathrm{PrEc}{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-1}{2}}\hfill \\ {(f″-\frac{f\prime }{\alpha +\delta })}^2=0,\hfill \end{array}$$

(18)

and the pressure term from (17) refers to as:

$$\begin{array}{c}P=\frac{\alpha +\delta }{2\delta }{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-1}{2}}(f″\prime +\frac{f″}{\alpha +\delta }-\frac{f\prime }{{(\alpha +\delta )}^2})\hfill \\ +(m-1)\mathrm{We}{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-3}{2}}(f″\prime -\frac{f″}{\alpha +\delta }+\frac{f\prime }{{(\alpha +\delta )}^2})\hfill \\ {(f″-\frac{f\prime }{\alpha +\delta })}^2+\frac{\delta }{\alpha +\delta }\mathrm{ff}″+\frac{\delta }{{(\alpha +\delta )}^2}\mathrm{ff}\prime -\frac{\delta }{\alpha +\delta }{f\prime }^2\hfill \\ -\frac{M{\delta }^2}{{(\alpha +\delta )}^2}f\prime ,\hfill \end{array}$$

(19)

by utilizing equations (16) into (17), one can eliminate the pressure term

$$\begin{array}{c}(m-1)(m-3)W{e}^2{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-5}{2}}{(f″\prime -\frac{f″}{\alpha +\delta }+\frac{f\prime }{{(\alpha +\delta )}^2})}^2\hfill \\ {(f″-\frac{f\prime }{\alpha +\delta })}^3+(m-1)\mathrm{We}{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-3}{2}}{(f″-\frac{f\prime }{\alpha +\delta })}^2\hfill \\ (f^{\mathrm{iv}}+\frac{f″}{{(\alpha +\delta )}^2}-\frac{f\prime }{{(\alpha +\delta )}^3})+(m-1)\mathrm{We}{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-3}{2}}\hfill \\ (f″\prime -\frac{f″}{\alpha +\delta }+\frac{f\prime }{{(\alpha +\delta )}^2})(3f″\prime -\frac{f″}{\alpha +\delta }+\frac{f\prime }{{(\alpha +\delta )}^2})(f″-\frac{f\prime }{\alpha +\delta })\hfill \\ +{[1+\mathrm{We}{(f″-\frac{f\prime }{\alpha +\delta })}^2]}^{\frac{m-1}{2}}(f^{\mathrm{iv}}+\frac{2f″\prime }{\alpha +\delta }-\frac{f″}{{(\alpha +\delta )}^2}+\frac{f\prime }{{(\alpha +\delta )}^3})\hfill \\ +\frac{\delta }{(\alpha +\delta )}\mathrm{ff}″\prime +\frac{\delta }{{(\alpha +\delta )}^2}\mathrm{ff}″-\frac{\delta }{{(\alpha +\delta )}^3}\mathrm{ff}\prime -\frac{\delta }{{(\alpha +\delta )}^2}{f}^{\prime 2}-\frac{\delta }{(\alpha +\delta )}f\prime f″\hfill \\ -\frac{M{\delta }^2}{{(\alpha +\delta )}^2}f″-\frac{M{\delta }^2}{{(\alpha +\delta )}^3}f\prime =0,\hfill \end{array}$$

(20)

boundary conditions are:

$$f(0)=S,f\prime (0)=1,\theta (0)=1,$$

(21)

$$f\prime (\alpha )=0,f″(\alpha )=0,\theta (\alpha )=0\mathrm{at}\alpha \to \infty .$$

(22)

Where $S=-\frac{v_w}{\sqrt{c}\nu }$ $S>0$ shows suction, $\mathrm{We}={\mathrm{\Gamma }}^2{c}^{2}R{e}_s$ is the Weissenberg number, $\delta =\sqrt{\frac{c}{\nu }}R$ is the dimensionless radius of curvature and $M=\frac{\sigma B_o^2}{c\rho }$ is the magnetic parameter.

Physical quantities

The quantity of physical interest in this problem is the skin friction coefficient $C_{\mathrm{fs}}$ and heat transfer coefficient $N{u}_s$ which is defined as

$$C_{\mathrm{fs}}=\frac{{\tau }_{\mathrm{rs}}}{\rho u_w^2},N{u}_s=\frac{s{q}_w}{k_o(T_w-T_{\infty })},$$

(23)

where

$${\tau }_{\mathrm{rs}}={\mu }_o{\left[{(1+{\mathrm{\Gamma }}^2{(\frac{\partial u}{\partial r}-\frac{u}{r+R})}^2)}^{\frac{m-1}{2}}(\frac{\partial u}{\partial r}-\frac{u}{r+R})\right]}_{r=0},$$

(24)

$$q_w=-k_o{\left(\frac{\partial T}{\partial r}\right)}_{r=0},$$

(25)

using equations (13)–(15) in (35), we obtain;

$${\mathrm{Re}}_s^{\frac{1}{2}}{C}_{\mathrm{fs}}={[1+\mathrm{We}{(f″(0)-\frac{f\prime (0)}{\delta })}^2]}^{\frac{m-1}{2}}(f″(0)-\frac{f\prime (0)}{\delta })$$

(26)

$${\mathrm{Re}}_s^{-\frac{1}{2}}N{u}_s=-\theta \prime (0).$$

(27)

The Reynolds number is $R{e}_s=\frac{c{s}^2}{\nu }$ .

Numerical scheme

Solutions to boundary value problems related to fluid flow can be determined using appropriate numerical techniques. The shooting method, one of the oldest and most efficient numerical methods, is often used. This method initially transforms the boundary value problem (BVP) into a system of first-order initial value problems (IVP). Then, the missing initial conditions at the starting point are estimated. This transformed IVP is solved using a highly efficient sixth-order Runge-Kutta method. The accuracy of the resulting solution is verified by comparing it with the given values at the end point. If the solution’s accuracy meets the desired level, it is accepted. Otherwise, the estimated values are iteratively refined using the Newton-Raphson method until the boundary conditions are satisfied within a specified tolerance limit at infinity.

Shooting method

The non-linear ODE for heat and momentum equations that is, subject to the boundary conditions are solved numerically using sixth order Runge-Kutta method along with shooting technique.

$$f\prime =u,f″=v,f″\prime =w,\theta \prime =P,$$

(28)

$$u_1=\frac{1}{\alpha +\delta },u_2=\frac{1}{{(\alpha +\delta )}^2},u_3=\frac{1}{{(\alpha +\delta )}^3},u_4=\frac{\delta }{\alpha +\delta },$$

(29)

$$u_5=\frac{\delta }{{(\alpha +\delta )}^2},u_6=\frac{\delta }{{(\alpha +\delta )}^3},u_7=\frac{M{\delta }^2}{{(\alpha +\delta )}^2},u_8=\frac{M{\delta }^2}{{(\alpha +\delta )}^3},$$

(30)

$$u_9=\frac{\delta \mathrm{Pr}}{\alpha +\delta },u_{10}=\mathrm{PrEc}$$

(31)

$$\begin{array}{c}(m-1)(m-3)\mathrm{We}{[1+\mathrm{We}{(v-u{u}_1)}^2]}^{\frac{m-5}{2}}{(w-v{u}_1+u{u}_2)}^2{(v-u{u}_1)}^3\hfill \\ +(m-1)\mathrm{We}{[1+\mathrm{We}{(v-u{u}_1)}^2]}^{\frac{m-3}{2}}{(v-u{u}_1)}^2(w\prime +v{u}_2+u{u}_3)+\hfill \\ (m-1)\mathrm{We}{[1+\mathrm{We}{(v-u{u}_1)}^2]}^{\frac{m-3}{2}}(w-v{u}_1+u{u}_2)(3w-v{u}_1+u{u}_2)\hfill \\ (v-u{u}_1){[1+\mathrm{We}{(v-u{u}_1)}^2]}^{\frac{m-1}{2}}\hfill \\ [w\prime +2w{u}_1-v{u}_2+u{u}_3]+u_4\mathrm{fw}+u_5\mathrm{fv}-u_5{u}^2-u_4\mathrm{uv}-u_{7}v-u_{8}u=0,\hfill \end{array}$$

(32)

$$P\prime +P{u}_1+u_9(f_p-u_{\theta })+u_{10}{[1+\mathrm{We}(v-u{u}_1)]}^{\frac{m-1}{2}}$$

(33)

The boundary conditions are defined as:

$$f(0)=S,U(0)=1,\theta (0)=1,$$

(34)

$$U(\infty )=0,V(\infty )=0,\theta (\infty )=0.$$

(35)

Results

Discussion

The current study examines the behavior of Carreau fluid over a stretching sheet. It explores how various physical parameters, such as the radius of curvature, Weissenberg number, magnetic parameter, and suction parameter, influence the velocity, temperature, and pressure of the fluid. This is considered for both shear thinning $(m<1)$ and shear thickening $(m>1)$ scenarios. The impact of an increasing radius of curvature on velocity is depicted in Figures 2 and 3. It’s observed that for both shear thinning and shear thickening cases, the velocity profile decreases. This is because an increase in the radius of curvature reduces the centripetal force responsible for circular motion, leading to a decrease in velocity. The Weissenberg number, a viscoelastic parameter of the fluid, signifies the degree of deformation. As the Weissenberg number increases, the fluid’s thickness also increases, resulting in a decline in the velocity profile. This is because the Weissenberg number measures the fluid’s relaxation time. Figures 4 and 5 illustrate the effect of the Weissenberg number for shear thinning and shear thickening cases. In the shear thickening case, an increase in the Weissenberg number leads to an increase in the particles’ relaxation time, which in turn increases the resistance for molecules, resulting in a decrease in velocity. However, the opposite behavior is observed in the shear thinning case. A magnetic field applied perpendicular to the particle’s velocity directly affects the particle’s motion direction. An increase in the magnetic force reduces the radius of the circular path, which in turn decreases the particle’s velocity. Figures 6 and 7 demonstrate the decrease in the velocity profile with an increase in the magnetic field parameter M for both $(m<1)$ and $(m>1)$ cases.

Figure 2.

Radius of curvature $\delta$ impact on velocity profile.

Figure 3.

Radius of curvature $\delta$ impact on velocity profile.

Figure 4.

Weissenberg number $\mathrm{We}$ impact on velocity profile.

Figure 5.

Weissenberg number $\mathrm{We}$ impact on velocity profile.

Figure 6.

Magnetic parameter $M$ impact on velocity profile.

Figure 7.

Magnetic parameter $M$ impact on velocity profile.

Suction parameter represents the amount of fluid being removed from the system. It means that suction parameter reduces the amount of fluid from system which results in decrease of velocity profile of fluid. Figures 8 and 9 deliberates the effect of suction parameter by decrease in velocity profile for both $M$ for $(m<1)$ and $(m>1)$ .

Figure 8.

Suction parameter $S$ impact on velocity profile.

Figure 9.

Suction parameter $S$ impact on velocity profile.

Figures 10 and 11 are sketch to know the behavior of radius of curvature $\delta$ on velocity profile for shear thinning and shear thickening fluid. It is seen in both cases that velocity decreases for increasing values of $\delta$ . Increase in Weissenberg number is directly related to increase in turbulence of fluid. The turbulence in return causes the increase in heat transfer. Whereas increase in magnetic field hinders the motion of particles of fluid which results in decrease of motion of particles. With increase in Weissenberg number the thermal field shows up surd for smaller value of magnetic field for $(M=0.5)$ Figures 12 and 13. However, an opposite behavior is noticed for greater value of magnetic parameter with $(M=1.5)$ . Figures 14 and 15 have been provided to describe the features of temperature against the magnetic parameter $M$ for $(m<1)$ and $(m>1)$ .

Figure 10.

Radius of curvature d impact on temperature profile with n = 0.5.

Figure 11.

Radius of curvature d impact on temperature profile with m = 0.5.

Figure 12.

Weissenberg number We impact on temperature profile with m = 0.5.

Figure 13.

Weissenberg number We impact on temperature profile with m = 1.5.

Figure 14.

Magnetic parameter M impact on temperature profile with m = 0.5.

Figure 15.

Magnetic parameter M impact on temperature profile with m = 1.5.

Suction parameter refers to the amount of fluid being removed from system. When the amount of fluid reduces, the internal energy of the particles also reduces which effects the temperature of system. Figures 16 and 17 portrays the effect of increase in suction parameter on temperature profile which shows decline with increases in suction parameter.

Figure 16.

Suction parameter S impact on temperature profile with m = 0.5.

Figure 17.

Suction parameter S impact on temperature profile with m = 1.5.

Figures 18 and 19 interpret the influence of Prandtl number on the thermal field. Both cases showed decline in temperature profile.

Figure 18.

Prandtl number Pr impact on temperature profile with m = 0.5.

Figure 19.

Prandtl number Pr impact on temperature with m = 1.5.

Increases in Eckert number produces heat in fluid due to frictional heating. Eckert number is fraction of kinetic energy to specific heat enthalpy. So when the Eckert number increases, it works against fluid stresses to transform kinetic energy into internal energy. Figures 20 and 21 display the increase in temperature profile in both cases by increasing Eckert number $\mathrm{Ec}$ .

Figure 20.

Eckert number Ec impact on temperature profile with m = 0.5.

Figure 21.

Eckert number Ec impact on temperature profile with m = 1.5.

When the radius of curvature keep increasing the surface become flat hence causes a decrease in pressure of fluid. The effect of increasing radius of curvature and Magnetic parameter on the pressure is plotted in Figures 22 and 23. For large values of $\delta$ the velocity decreases, while the pressure tends to zero for the flat surface. Figure 23 shows the effect of radius of curvature on pressure. Figure 24 displays the behavior of powerlaw index on pressure profile. It is observed that pressure decreases for large values of $m$ . With increases in power law index parameter, the fluid become less viscous, and it can flow more conveniently. So, fluid will require less pressure to flow. Therefore, in Figure 24 we noticed a decrease in pressure with increases in value of power law index parameter $m$ .

Figure 22.

Radius of curvature $\delta$ impact on pressure profile.

Figure 23.

Magnetic parameter $M$ impact on temperature profile.

Figure 24.

Power-law index $m$ impact on temperature profile.

Table 1 is presented to show the validity of the present results with the previous published work. The results for all values of $\delta$ are align with previous published works. The distinction in the skin friction coefficient and the heat transfer rate owing to the alteration of $\mathrm{We}$ , $M$ , $\delta$ , $\mathrm{Ec}$ , and $\mathrm{Pr}$ are computed to tabulate the numerical values for shear thinning $(m<1)$ and shear thickening $(m>1)$ case in Tables 2 and 3. Both the skin friction coefficient and local heat transfer rate experience a decrease in magnitude with an increase in $\mathrm{We}$ for shear thinning $(m<1)$ and it increases for shear thickening $(m>1)$ case. The skin friction coefficient increases with an increase in Magnetic parameter $M$ and suction $S$ , but decreases with an increase in radius of curvature parameter $\delta$ . An increase in $M$ , $\delta$ and $\mathrm{Ec}$ results in a decrease in local heat transfer rate, while an increase in $\mathrm{Pr}$ leads to an increase in local heat transfer rate.

Table 1.

Comparison of present values of the local friction coefficient $-{\mathrm{Re}}_s^{\frac{1}{2}}{C}_{\mathrm{fs}}$ with Rosca and Pop ²¹ for several values of $\delta$ .

Present values
$\delta$	Rosca and Pop 21	SM	bvp4c	Percentage error
5	1.1516	1.1516	1.1516	0
10	1.0734	1.0719	1.0719	8.405
20	1.0350	1.0350	1.0350	0
30	1.0231	1.0231	1.0231	0
40	1.0171	1.0173	1.0173	0.0196
50	1.0318	1.0138	1.0138	0
100	1.0068	1.0068	1.0068	0
200	1.0034	1.0034	1.0034	0
1000	1.0006	1.0006	1.0006	0
$\infty$	1.0000	1.0000	1.0000	0

Table 2.

Numerical values of the local friction coefficient $-{\mathrm{Re}}_s^{\frac{1}{2}}{C}_{\mathrm{fs}}$ for different physical parameters.

$\mathrm{We}$	$M$	$\delta$	$S$	$n=0.5$		$n=1.5$
				SM	bvp4c	SM	bvp4c
0	0.5	5	0.5	1.6405	1.6405	1.6405	1.6405
1				1.4678	1.4678	1.7720	1.7720
2				1.3779	1.3779	1.8388	1.8388
0.6	0.2	5	0.5	1.4092	1.4092	1.5811	1.5811
	0.4			1.4860	1.4860	1.6849	1.6849
	0.6			1.5558	1.5558	1.7817	1.7817
0.6	0.5	5	0.5	1.5217	1.5217	1.7341	1.7341
		10		1.4648	1.4648	1.6467	1.6467
		1000		1.4128	1.4128	1.5690	1.5690
0.6	0.5	5	0	1.2978	1.2978	1.4785	1.4785
			0.4	1.4730	1.4730	1.6790	1.6790
			0.8	1.6795	1.6795	1.9109	1.9109

Table 3.

Numerical values of local heat transfer coefficient for different physical parameters.

$\mathrm{We}$	$M$	$\delta$	$\mathrm{Pr}$	$\mathrm{Ec}$	$n=0.5$		$n=1.5$
					SM	bvp4c	SM	bvp4c
0	0.5	5	0.7		0.8896	0.8896	0.8896	0.8896
1					0.8492	0.8492	0.9157	0.9157
2					0.8261	0.8261	0.8261	0.9271
0.6	0.2	5	0.7		0.9039	0.9039	0.9432	0.9432
	0.4				0.8748	0.8748	0.91918	0.9191
	0.6				0.8503	0.8503	0.8985	0.8985
0.6	0.5	5	0.7		0.8623	0.8623	0.9087	0.9087
		10			0.8507	0.8507	0.8945	0.8945
		1000			0.8376	0.8376	0.8789	0.8789
0.6	0.5	5	0.7		0.8448	0.8448	0.8905	0.8905
			1		1.0985	1.0985	1.1494	1.1494
			0.2		1.3297	1.3297	1.3857	1.3857
0.6	0.5	5		0	0.9739	0.9739	1.02896	1.0289
				0.2	0.8623	0.8623	0.9087	0.9087
				0.4	0.7507	0.7507	0.7884	0.7884

The above promising results gives future directions for this model to utilize for further industrial approaches by considering the effect of entropy generation. Moreover, for other non-Newtonian fluids this work can be carry forward.

Conclusions

Numerical investigations have been carried out on the flow of the Carreau fluid along with a linear stretchable curved wall. On the basis of our investigation, the key points of the pre-sent analysis are given below.

• For shear thinning fluids, increases in Weissenberg number intensify the velocity of fluid while it mitigates the velocity for shear thickening fluids.

• An increase in the magnetic parameter, curvature factor and suction parameter results in a decrease in the flow velocity. On the other hand, fluid velocity increases with an increase in the power law index.

• For shear thinning fluids, the temperature profile increases while it decreases for shear thickening fluids when the Weissenberg number is large.

• The heat transfer rate decreases with an increase in the suction parameter, curvature factor, power-law index and Prandtl number. Conversely, the opposite behavior is observed for the magnetic parameter.