Dual solution of boundary-layer flow driven by variable plate and streaming-free velocity

Abstract

This article presents a numerical study to investigate boundary-layer heat transfer fluid associated with a moving flat body in cooperation of variable plate and streaming-free velocity along the boundary surface in the laminar flow. The thermal conductivity is supposed to vary linearly with temperature. Similarity transformations are applied to render the governing partial differential equations for mass, momentum and energy into a system of ordinary differential equations to reveal the possible existence of dual solutions. MATLAB package has been used to solve the boundary value problem numerically. We present the effects of various parameters such as velocity ratio, thermal conductivity and variable viscosity on velocity and temperature distribution. The analysis of the results concerning Skin friction and Nusselt number near the wall is also presented. It is focused on the detection and description of the dual solutions. The study reveals that the undertaken problem admits dual solutions in particular range of values of different physical parameters. It can be seen that for the first branch solution, the fluid velocity decreases near the sheet, but it increases far away from the sheet for velocity ratio parameter, whereas the opposite effect is induced for second branch solution. Skin friction coefficient and rate of heat transfer increase due to increase in thermal conductivity parameter.

Keywords

Convective flow variable plate dual solution flow profiles upper and lower branches

Introduction

Studies on boundary-layer flow over continuously moving surfaces are important due to their several applications in manufacturing processes such as glass-fibre production, the aerodynamic extrusion of plastic sheets, metal extrusion, materials-handling conveyors, and paper production.

The problem of a liquid over flat surface in uniform free stream was first elucidated by Blasius. This was done by introducing a new independent variable called the similarity variable. The similarity variable has been applied to solve thermal boundary layer for constant surface temperature along the plate. Steady free-convective and energy (heat) transfer through isothermal body was studied by Mahanti and Gaur¹ with a view to visualizing the effects of variable viscosity and conductivity properties. Magnetohydrodynamics (MHD) natural convection along with convective boundary surface condition past a flat surface has been reported by Rashidi et al.²

The boundary-layer flow and heat transfer over a continuous flat plate/surface moving in parallel or in reverse to the free stream was discussed by many researchers. The study of boundary-layer flow on a constant speed moving plate was first studied by Sakiadis.³ Abel et al.⁴ analysed the MHD, thermal radiation and buoyancy effects on a moving stretching surface. Elbashbeshy and Bazid⁵ reported the boundary-layer flow over a moving plate with the effect of force convection and temperature-dependent viscosity. Jhankal⁶ presented the two-dimensional MHD flow and heat transfer along an infinite porous continuous moving plate. The aim of this article is to study the combined effects of heat generation and chemical reaction on MHD natural convection flow over a moving plate embedded in a porous medium. Khan et al.⁷ studied the combined effects of heat generation and chemical reaction on MHD natural convection flow over a moving plate embedded in a porous medium. Sandeep and Sugunamma⁸ analysed the effects of radiation and rotation on unsteady hydromagnetic free convection flow of a viscous incompressible electrically conducting fluid past an impulsively moving vertical plate in a porous medium. Flow of an incompressible viscous fluid past a continuously moving semi-infinite plate is studied by Soundalgekar et al.,⁹ and they considered the variable viscosity and variable temperature. The effect of temperature dependent viscosity on a laminar mixed convection boundary-layer flow and heat transfer on continuous moving vertical surface is studied by Mohamed and Ali.¹⁰ Pop et al.¹¹ studied the flow and heat transfer problems, including the variation of fluid viscosity with temperature. Pantokratoras¹² analysed the boundary-layer flow along an isothermal, continuously moving plate, which is studied taking into account the variation of fluid viscosity with temperature.

Hassanien et al.¹³ presented the boundary-layer flow and heat transfer of an incompressible micropolar fluid over a plane moving surface in parallel or in reverse to the free stream and they considered the boundary condition is isothermal. Salleh et al.¹⁴ analysed the boundary-layer flow and heat transfer over a stretching sheet with Newtonian heating. Lin and Hung¹⁵ studied the flow and heat transfer of a plane surface moving in parallel and in reverse to the free stream.

The mathematical analysis of the dual solutions of the boundary-layer flow over the moving surfaces has a practical impact in the engineering scenario. It gives the possibility to determine the most realistic, stable, and physically acceptable solutions. Ahmed et al.¹⁶ analysed the dual solution of steady-state, boundary-layer flow of a power law fluid along a moving flat plate in the presence of the thermal radiation, viscous dissipation, and internal heat generation or absorption. Afzal et al.¹⁷ presented the dual solutions of the boundary-layer velocity and temperature distributions. Mukhopadhyay and Golra¹⁸ generalized the dual solutions of boundary-layer flow of fluid along a moving surface with power law temperature. Furthermore, Mukhopadhyay¹⁹ also analysed the dual solution of boundary-layer flow of a moving fluid over a moving permeable surface in the presence of prescribed surface temperature and thermal radiation. Deswita et al.²⁰ reported the dual solutions for the power law fluid flow with suction and injection effects and neglected the heat transfer effects.

Our aim is to derive systematically the similarity transformation under similarity requirement for the governing equations, and then, we have transformed highly non-linear partial differential equations (PDEs) fluid flow of boundary-layer theory into ordinary differential equations (ODEs) by introducing the group of transformation. The reduced non-linear problem is numerically sketched by using bvp4c in MATLAB computational algorithm. The study reveals that there exist dual solutions in specific range of the vital parameters involved and that one of the two solutions is stable and physically realistic.

Model and transformations

We assume boundary-layer flow of a steady, two-dimensional, incompressible liquid past a moving flat surface. Density of the liquid is taken to be ρ. We consider a Cartesian coordinate system $(x, y)$ in which $x$ is measured along the body surface and y-axis is measured along the direction normal to the plate. We set the free stream of velocity $U_{e}$ and the flat body of velocity $U_{p}$ . It is presumed that $T_{w}$ is the temperature of the plate, while $T_{\infty}$ is the uniform temperature of the free stream.

Under the boundary-layer and the Boussinesq approximations, the governing equations for this physical problem are as follows^21,22

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

(1)

u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = U_{e} \frac{d U_{e}}{dx} + \frac{1}{ρ} \frac{\partial}{\partial y} (μ (t) \frac{\partial u}{\partial y})

(2)

ρ C_{p} (u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y}) = \frac{\partial}{\partial y} k (\frac{\partial T}{\partial y})

(3)

where $(uv)$ is the velocity components in $(x, y)$ directions, $C_{p}$ is the specific heat, $k$ is the variable thermal conductivity, and $μ (t)$ is the temperature-dependent dynamic viscosity.

We assume as in Sharma and Singh,²³ $U_{p} (x) = x^{n} U_{w}$ , $U_{e} (x) = x^{n} U_{\infty}$ , where $U_{w}$ and $U_{\infty}$ are constant reference velocities and $\bar{u} (x) = (U_{\infty} + U_{w}) x^{n} = B_{0} x^{n}$ , where $B_{0} = U_{\infty} + U_{w}$ .

Subject to the constraints

u = U_{w} (x), v = 0, T = T_{w} at y = 0

(4)

u = U_{e}, T = T_{\infty} at y \to \infty

(5)

The following transformation is applied to the system of governing equations

\begin{matrix} u = B_{0} x^{n} f' (η), η = \frac{y \sqrt{R e_{x}}}{x}, ξ = \frac{U_{\infty}}{U_{w} + U_{\infty}}, \\ θ = \frac{T - T_{\infty}}{T_{w} - T_{\infty}}, R e_{x} = \frac{B_{0} x^{n} * x}{ν}, k = k_{\infty} [1 + δ θ], \\ δ = \frac{k_{w} - k_{\infty}}{k_{\infty}}, μ (t) = \frac{μ_{\infty}}{1 + θ γ (T - T_{\infty})} \end{matrix}

(6)

Thus, the reduced form of our considered problem equations takes the form

\begin{matrix} n [{(f')}^{2} - ξ^{2}] - (\frac{n + 1}{2}) ff ″ \\ = (\frac{1}{1 + ε θ}) f ″' - (\frac{1}{{(1 + ε θ)}^{2}}) f ″ ε θ' \end{matrix}

(7)

(1 + δ θ) θ ″ + δ θ'^{2} + \Pr (\frac{n + 1}{2}) f θ' = 0

(8)

Subject to the constraints

\begin{matrix} f (0) = 0, f' (0) = 1 - ξ, θ (0) = 1, at η = 0 \\ f' (\infty) = ξ, θ (\infty) = 0, at η \to \infty \end{matrix}}

(9)

where $ε$ is the variable viscosity parameter, $δ$ is the thermal conductivity parameter, $ξ$ is the stretching ratio parameter, $\Pr$ is the Prandtl number, and $n$ is the pressure number.

Numerical technique

Now we solve the set of non-linear ordinary differential equations (7) and (8) with boundary conditions (9) numerically by using bvp4c function technique in MATLAB package. We consider $f = y_{1}, f' = y_{2}, f ″ = y_{3}, θ = y_{4}, θ' = y_{5}$ . Then, equations (7) and (8) are transformed into a system of first-order ordinary differential equations as given in the following

\begin{matrix} f' = y_{2} \\ f ″ = {y'}_{2} = y_{3} \\ f ″' = {y'}_{3} = n (y_{2}^{2} - ξ^{2}) \times (1 + ε \times y_{4}) - \frac{n + 1}{2} y_{1} y_{3} \times (1 + ε \times y_{4}) + \frac{ε}{(1 + ε \times y_{4})} y_{3} y_{5} \\ θ' = y_{5} \\ θ ″ = {y'}_{5} = - \Pr \times \frac{(n + 1)}{2} \times \frac{1}{1 + δ \times y_{4}} y_{1} y_{5} - \frac{δ}{1 + δ \times y_{4}} \times y_{5}^{2} \end{matrix}}

(10)

along with the initial boundary conditions

y_{1} (0) = 0, y_{2} (0) = 1 - ξ, y_{4} (0) = 1, y_{2} (\infty) = ξ, y_{4} (\infty) = 0

(11)

Equations (10) and (11) are integrated numerically as an initial value problem to a given terminal point. All these simplifications are made by using bvp4c function available in MATLAB software.

Results and discussion

The coupled systems subject to boundary conditions were solved numerically by using the bvp4c in MATLAB package for different values of $δ, ε, ξ, \Pr, n$ . The package was a boundary value procedure. The boundary transformed conditions (at infinity) of η: $η_{\max} = 20, f' (20) = ξ, θ (20) = 0$ .

The choice of $η_{max} = 20$ confirms that the far-field boundary conditions are satisfied precisely. Solutions for a range of problem parameters are, however, useful since they illustrate the main features of the boundary-layer flow. In this study, we have considered the values of viscosity variation parameter, $ε = 0.1, 0.2, 0.3, 0.4$ ; velocity ratio, $ξ = 0.05, 0.1, 0.15, 0.2, 0.30.4, 0.5, 0.6, 0.8$ ; pressure gradient $n = 0.3, 0.4, 0.5, 0.6, 1.0$ ; and thermal conductivity parameter, $δ = 0.1, 0.2, 0.3$ . At $0^{°} C$ , we have used Prandtl number 0.72 for air, 0.0288 for Mercury and 13.6 for water.

In order to test the accuracy of the present results, we have compared the obtained results with those of Hassanien et al.¹³ and Salleh et al.¹⁴ when we neglect the effects of some parameters. This comparison is shown in Table 1 and we notice an excellent agreement with those of Hassanien et al.¹³ and Salleh et al.¹⁴ This observation serves as a confirmation of the accuracy of the results that are performed in this communication.

Table 1.

Comparison results for $θ' (0)$ when $n = 1, ε = δ = ξ = 0$ with previously published results for different Prandtl number.

Prandtl number	Present solution	Hassanien et al.¹³	Salleh et al.¹⁴
0.72	0.465014	0.46325	0.46317
1	0.582228	0.58198	0.58198
3	1.165218	1.16525	1.16522
5	1.568031		1.56806
7	1.895382		1.89548
10	2.307985	2.30801	2.30821
100	7.765803	7.74925	7.76249

In this section, we analysed the dual solutions of the boundary-layer velocity and temperature profiles. Dual solution refers to two different solutions that are obtained under the same condition by assuming several values of missing initial guesses. One of these two solutions is called the upper or first-branch solution and other one is called the lower or second-branch solution. In the following figures, let us specify that the solid lines represent the first or upper branch solution, while dotted lines represent the second or lower branch solution (Figure 1).

Figure 1.

Schematic diagram of the problem.

Figures 2 and 3 exhibit the effect of free-stream velocity ratio parameter $(ξ)$ on the non-dimensional velocity and temperature profiles. From these figures, we see that two solution branches exist for some particular governing parameter, which support the existence of dual solutions. For the case of first branch of solution, we also observe that near the sheet, the velocity is dropped with free-stream velocity ratio parameter, whereas opposite behaviour is shown for larger value of $η$ , that is, far away from the sheet. For the second branch of solution, velocity decreases with free-stream velocity ratio parameter $(ξ)$ for $0 \leq η \leq 10$ , whereas opposite behaviour is shown for $10 \leq η \leq 20$ . On the other hand, temperature distribution is lessened with the increase in free-stream velocity ratio parameter $(ξ)$ in the first solution, whereas reverse trend is shown in the second solution. From the numerical simulation, it is also noticed that the boundary-layer thickness for the first branch solution is smaller compared to the second branch solution.

Figure 2.

Velocity profile for various values of $ξ$ .

Figure 3.

Temperature profile for various values of $ξ$ .

The discrepancy of skin friction coefficient and rate of heat transfer with free-stream velocity ratio parameter $(ξ)$ for various values of velocity power index parameter (n) are shown in Figures 4 and 5.

Figure 4.

Variation of skin friction coefficient with $ξ$ for different values of n.

Figure 5.

Variation of rate of heat transfer coefficient with $ξ$ for different values of n.

According to the value of free-stream velocity ratio parameter, the following observation can possibly occur:

When $0 < ξ < 1,$ the plate and the fluid move in the same direction.

$ξ = 0$ , the plate is fixed.

$ξ < 0$ and $ξ > 1$ , the plate and the fluid move in the opposite direction.

$ξ = 1$ , the plate is moving at constant velocity and the fluid velocity is zero.

In the present analysis, we found that the dual solutions exist for a certain value of governing parameters. For velocity power index parameter $n = 0.5$ , the solution is unique for $ξ \geq 0.465$ and $ξ \leq 0.185$ . Dual solutions exist for $0.185 \leq ξ \leq 0.465$ , and for $n = 0.7$ , the solution is unique for $ξ \geq 0.285$ and $ξ \leq 0.185$ . Dual solutions exist for $0.185 \leq ξ \leq 0.285$ . Note that the range of the similarity solution is increased as velocity power index parameter becomes greater. We have also observed that skin friction coefficient increases with the rise in free-stream velocity ratio ( $ξ$ ). According to our observation for heat transfer rate from Figure 5, we can see that for the velocity power index parameter $n = 0.7$ , the solution is unique for $ξ \leq - 0.195$ and $ξ \geq 0.175$ and dual solutions exist for $- 0.195 \leq ξ \leq 0.175$ . It is also perceived from our study that the rate of heat transfer increases non-linearly with stretching ratio parameter; however, it decreases with velocity power index parameter.

Figures 6 and 7 show the variation in the skin friction coefficient and rate of heat transfer with free-stream velocity ratio parameter $(ξ)$ for various viscosity parameter $(ε)$ . These figures demonstrate that there exist two solution branches. The first branch represents the stable solution, while the second branch denotes the unstable solution for each value of $(ξ)$ corresponding to a given value of viscosity parameter $(ε)$ . From Figure 6, we observe that the solution is unique when $ξ = ξ_{c}$ ; however, dual solutions exist when $ξ < ξ_{c}$ and no solutions exist when $ξ > ξ_{c}$ where the value of $ξ_{c} = 0.225, 0.065, - 0.445 .$ with specific values of $ε = 0.1, 0.2, 0.3 .$ $ξ_{c}$ is the critical value of $ξ$ , at which the two solution branches coincide and thus a unique solution is obtained.

Figure 6.

Variation of skin friction coefficient with $ξ$ for different values of $ε$ .

Figure 7.

variation of rate of heat transfer with $ξ$ for different values of $ε$ .

Variation of rate of heat transfer rate $θ' (0)$ on the wall with free-stream velocity parameter for various values of the viscosity parameter are shown in Figure 7. From this figure, it can be apprehended that for $ε = 0.1$ , the solution is unique when $ξ \leq - 0.195$ and $ξ \geq 0.175$ , while dual solutions exist when $- 0.195 \leq ξ \leq 0.175$ for some specific values of the governing parameter as $n = 0.7, \Pr = 0.72, δ = 0.1$ .

Furthermore, Figure 8 displays the variation of rate of heat transfer coefficient with $ξ$ for several values of $δ$ . From this graph, we observe that the dual solution exists. For $δ = 0.1$ , the solution is unique when $ξ \leq - 0.195$ and $ξ \geq 0.175$ , while dual solutions exist when $- 0.195 \leq ξ \leq 0.175$ for some specific values of the governing parameter as $n = 0.7, \Pr = 0.72, δ = 0.1 .$ We can inscribe here that the range of the similarity solution is decreased with the rise of thermal conductivity parameter $(δ)$ . It can also be seen that rate of heat transfer increases with $δ$ for both first- and second-branch solutions. We have also observed that rate of heat transfer is increased non-linearly with stretching ratio parameter. It is noticed that the second solution has a higher value collated to the first solution.

Figure 8.

Variation of rate of heat transfer with $ξ$ for different values of $δ$ .

Conclusion

The effect of variable viscosity and thermal conductivity with velocity power index parameter on the boundary-layer flow on a moving plate with free-stream velocity is analysed. Using suitable similarity analysis, the transformed non-dimensional boundary value problem has been solved with MATLAB symbolic software using shooting quadrature. The analysis reveals that the dual solution exists for some specific value of governing parameter. The current simulations have shown that

Fluid velocity decreases near the sheet but it increases far away from the sheet by enhancing velocity ratio parameter in the case of first branch of solution, whereas the opposite effect is induced with velocity ratio parameter for second-branch solution.

Fluid temperature reduces with the rise of velocity ratio parameter in the case of first solution. On the contrary, opposite behaviour is observed for the second solution.

Increasing the free-stream velocity ratio is found to consistently increase the Skin friction coefficient and rate of heat transfer near the surface wall.

A sustained decrease is generated in the skin friction coefficient and rate of heat transfer near the surface wall with an increase in velocity power index parameter.

The rise of the viscosity parameter had resulted in the reduction of the skin friction coefficient and rate of heat transfer at the surface wall.

Due to increase in thermal conductivity parameter, skin friction coefficient and rate of heat transfer also increase.

Similarity solution range of skin friction and rate of heat transfer coefficient drops with free-stream velocity ratio parameter with an increase in viscosity and thermal conductivity parameter.

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

M Ferdows

References

Mahanti

Gaur

. Effects of varying viscosity and thermal conductivity on steady free convective flow and heat transfer along and isothermal vertical plate in the presence of heat sink. J Appl Fluid Mech 2009; 2: 23–28.

Rashidi

Ferdows

Parsa

, et al. MHD natural convection with convective surface boundary condition over a flat plate. Abstract Appl Anal 2014; 2014: 1–10.

Sakiadis

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

Abel

Prasad

Mahaboob

. Buoyancy force and thermal radiation effects in MHD boundary layer visco-elastic fluid flow over continuously moving stretching surface. Int J Therm Sci 2005; 44: 465–476.

Elbashbeshy

EMA

Bazid

MAA

. The effect of temperature-dependent viscosity on heat transfer over a continuous moving surface. J Phys D: Appl Phys 2000; 33: 2716.

Jhankal

. MHD boundary layer flow and heat transfer along an infinite porous hot horizontal continuous moving plate. Indian J Chem Technol 2018; 25: 106–110.

Khan

, et al. Effects of relative magnetic field, chemical reaction, heat generation and Newtonian heating on convection flow of Casson fluid over a moving vertical plate embedded in a porous medium. Sci Rep 2019; 9: 400.

Sandeep

Sugunamma

. Radiation and inclined magnetic field effects on unsteady hydromagnetic free convection flow past an impulsively moving vertical plate in a porous medium. J Appl Fluid Mech 2014; 7: 275–286.

Soundalgekar

Takhar

Das

, et al. Effect of variable viscosity on boundary layer flow along a continuously moving plate with variable surface temperature. Heat Mass Trans 2004; 40: 421–424.

10.

Mohamed

and Ali

. The effect of variable viscosity on mixed convection heat transfer on a vertical moving surface. Int J Therm Sci 2006; 45: 60–69.

11.

Pop

Gorla

RSR

Rashidi

. The effect of variable viscosity on flow and heat transfer to a continuous moving flat plate. Int J Eng Sci 1992; 30: 1–6.

12.

Pantokratoras

. Further results on the variable viscosity on flow and heat transfer to a continuous moving flat plate. Int J Eng Sci 2004; 42: 17–18, 1891–1896.

13.

Hassanien

Shamardan

Moursy

, et al. Flow and heat transfer in the boundary layer of a micropolar fluid on a continuous moving surface. Int J Numer Method Heat Fluid Flow 1999; 9: 643–659.

14.

Salleh

Nazar

Pop

. Boundary layer flow and heat transfer over a stretching sheet with Newtonian heating. J Taiwan Inst Chem Eng 2010; 41: 651–655.

15.

Lin

Huang

. Flow and heat transfer of plane surfaces moving in parallel and reversely to the free stream. Int J Heat Mass Transfer 1994; 37: 333–336.

16.

Ahmed

Siddique

Sagheer

. Dual solutions in a boundary layer flow of a power law fluid over a moving permeable flat plate with thermal radiation, viscous dissipation and heat generation/absorption. Fluids 2018; 3: 6.

17.

Afzal

Badaruddin

Elgarvi

. Momentum and transport on a continuous flat surface moving in a parallel stream. Int J Heat Mass Transf 1993; 36: 3399–3403.

18.

Mukhopadhyay

Golra

RSR

. Dual solutions for boundary layer flow of moving fluid over a moving surface with power-law surface temperature. Int J Appl Mech Eng 2013; 18: 113–124.

19.

Mukhopadhyay

. Dual solutions in boundary layer flow of a moving fluid over a moving permeable surface in presence of prescribed surface temperature and thermal radiation. Chin Phys B 2014; 23: 014702.

20.

Deswita

Ishak

Nazar

. Power-law fluid flow on a moving wall with suction and injection effects. Aust J Basic Appl Sci 2010; 4: 2250–2256.

21.

Mwaonanji

. Boundary layer flow over a moving flat surface with temperature dependent viscosity. Dar es Salaam, Tanzania: University of Dar Es Salaam, 2012.

22.

Vyas

Rai

. Radiative flow with variable thermal conductivity over a non–isothermal stretching sheet. Int J Contemp Math Sciences 2010; 5: 2685–2698.

23.

Sharma

Singh

. Effects of variable thermal conductivity and heat source/sink on MHD flow near a stagnation point on a linearly stretching sheet. J Appl Fluid Mech 2009; 2: 13–21.