Blog

Stability analysis of dual solutions for mixed convection and thermal radiation with hybrid nanofluid flow past shrinking/stretching curved surface | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  13, Article number: 21676 (2023 ) Cite this article Home Water Heater

Stability analysis of dual solutions for mixed convection and thermal radiation with hybrid nanofluid flow past shrinking/stretching curved surface | Scientific Reports

This study portrays the stability analysis and dual solutions of mixed convection and thermal radiation of hybrid nano-fluid flow past stretching/shrinking a curved surface in the presence of injection/suction conditions. A hybrid nano-fluid, in which water is used as the base fluid, copper and alumina are used as nano-particles, and the magnetic field is taken into account. The present study’s findings will provide fruitful implications for future research in the field of fluid dynamics. The bvp4c method using Matlab software is implemented to get the numerical solution of the nonlinear partial differential equation transformed into the ordinary differential equation. The behavior of the first and second solutions under the governing parameters on the curved surface of dimensionless velocity \(g'(\xi )\) , shear stress profile \(g''(\xi )\) , temperature profile \(\theta (\xi )\) , skin friction coefficient Cfs, and local Nusselt’s number Nus were visualized in figurative and tabular form. From this the following are investigated: as the values of \(\phi j\) increase, the velocity profile for the second solution decreases, and the opposite trend is observed for the first solution. For the values of K and \(\lambda 1\) , the shear stress profile increased for the first solution, and the opposite trend was observed for the second solution, though after some interval points, the inverse of this statement was observed.For the values of \(S,Pr,Rd,M, \;\text{and}\; \lambda\) , the upwind thermal boundary layer of the first solution is larger than the second solution.For the value of M uphill, the estimation of the absolute value of \(\lambda ci\) increases for both the skin friction coefficient and the local Nusselt number. In the second solution, increasing the values of \(\beta\) , Pr, S, Ec, and M has a similar effect on \(g''(0)\) , Cfs, -\(\theta '(0)\) , and Nus. In the first solution, increasing the values of Ec, S, and Pr on \(g''(0)\) and Cfs results in a decrease. The first solution has a positive eigenvalue, whereas the second solution has a negative eigenvalue. Agreement between the present analysis and literature is acceptable.

Nowadays, many researchers are paying attention to the examination of hybrid nano-fluid flow past a shrinking/stretching surface due to several implications in both industrial and manufacturing processes. These include polymer and wire drawing, hot rolling, rubber sheets, the performance of lubricants and paints, the aerodynamic extrusion of plastic sheets, metal spinning, and polymer extrusion. Hence, Yashkun et al.1 evaluated the characteristics of thermal radiation and suction influences on linear shrinking and stretching a sheet under magnetohydrodynamic hybrid nano-fluid for the existence of dual solutions. Waini et al.2 also discussed the thermal radiation effect on a vertically shrinking sheet with mixed convection magnetohydrodynamic flow in the presence of dust particles and \({{CuAl}}_{2} {{O}}_{3}\) nanoparticles. The stretching sheet in the presence of a heat sink/source, thermal radiation, and chemical reactions of first-order non-Newtonian nanofluid mixed convection flow was addressed by Hayat et al.3. Furthermore, Jamaludin et al.4 reported a numerical investigation of the vertical shrinking/stretching of a sheet under the influence of a heat source/sink, suction, and thermal radiation at the stagnation point in the presence of mixed convection flow and water as the base fluid, and with nanoparticles as titania \(TiO_{2}\) , copper Cu, and alumina Al2O3.

Searching for a new nanoparticle that enhances the heat transfer rate with minimum cost is the duty of noway researchers, but it is a challenging task for scientists in the field of fluid dynamics. The suspension of one type of nanoparticles into a nanofluid containing nanoparticles of other metals or non-metals can enhance the thermal properties to obtain a new kind of nanofluid, which is called hybrid nanofluid. In different industrial areas like nuclear system cooling, microelectronics, naval structures, manufacturing, biomedicine, and drug reduction, hybrid nanofluids have many applications. Yashkun et al.5 used water as the base fluid and copper and aluminum oxide as the nanoparticles, with Joule heating and mixed convection under an exponentially shrinking/stretching surface in the presence of heat transfer. A numerical investigation of unsteady mixed convection, magnetohydrodynamic flow of heat transfer, and hybrid nanofluid pasts of permeable vertical plates is reported by Wahid et al.6. Muhammad et al.7 computed the combined effects of mixed convection, slip, and dissipation due to the heating of hybrid nanofulids formed from Cu and MWCNTs past a surface with a curved-shaped body. Further, many researchers also investigated the effects of hybrid nanofluid flow8,9,10.

In varying practical applications such as transportation, electronics, etc., that is, a modern cutting-edge adhesive technologies, such as a flow past stretching a curved sheet have a wide range of applications. Some of the applications are growing crystals structures, plastic sheets preparation, manufacturing of electronic chips and materials, paper industry, cooling process. Due to this case, the curved sheet connected to the flow has come to the attention of more researchers in the present day. Researchers such as Rosca and Pop11 have reported on an electrical conduction fluid of incompressible, viscous, laminar, two-dimensional, unsteady flow with past shrinking sheets in the presence of a transversely identical magnetic field under a curved surface. Also, Revathi et al.12 investigated the influence of CH3OH + \(SiO_{2}\) + \(Al_{2}O_{3}\) , a methanol based hybrid nano-fluid, on a curved stretching surface in the presence of thermal radiation and cross diffusion due to increased activation energy and Brinkman number. The exception to the great output devices of entropy generation, activation energy due to forced natural convection flow over a curved surface of fully developed Darcy-Forchheimer is analyzed by Muhammad et al.13. Ahmed et al.14 developed the phenomena of heat and mass transport under a curving, expanding sheet of Cu-CuO dependent on sodium alginate, NaAlg hybrid nanofluid with boundary layer flow in the presence of nonlinear thermal radiation, chemical reactions, and a magnetic field. Ijaz et al.15 discussed the behavior of the entropy optimization rate of magnetohydrodynamic flow due to a curved stretching surface in the presence of Joule heating, heat generation, and viscous dissipation. Ibrahim and Gizewu16 reported the characteristic magnetohydrodynamic flow of a third-order slip boundary condition for a curved stretching sheet under gyrotactic microorganisms and bio-convective entropy generation in the presence of Dufour and Soret effects.

The phenomenon in which energy spreads from a heated surface to its absorption point in all directions in the form of electromagnetic waves is called thermal radiation. The thermal agitation of composite molecules in the body generates thermal radiation. Heating of the room due to an open fireplace is a common example of thermal radiation.The sun, the light bulb, and microwave radiation are typical examples. Furthermore, technological applications of it can be seen in solar power, nuclear power plants, combustion chambers, and chemical processes. The amount of heat energy emitted by the heat surface per unit area is directly proportional to the fourth power of the surface’s absolute temperature, as stated by Stefan Boltzmann’s law of radiative heat transfer. Wahid et al.17 evaluated the impacts of radiative mixed convection flow with a hybrid nano-fluid, copper-alumina/water, due to the existence of a magnetic field. Bejawada and Nandeppanavar18 analyzed the influence of thermal radiation in the presence of micro-polar fluid with magnetohydrodynamic heat transfer flow past a moving porous vertical plate. Lv et al.19 explored the behavior of activation energy and chemical reactions with two horizontal infinite plates without nano-fluid flow in the presence of variable thermal conductivity, the Hall effect, and thermal radiation past a permeable and stretchable sheet. Hussain et al.20 analyzed the impact of boundary layer flow on velocity slip, thermal radiation at heat absorption, and magnetohydrodynamics past a permeable exponentially stretching surface. The impact of magnetohydrodynamics, chemical reactions, and thermal radiation, in the presence of melting heat transfer, Brownian motion, and thermophoresis effects is due to a non-linear stretching surface, as discussed by Krishnamurthy et al.21. For more information, see22.

Mixed convection is a blend of free and forced convection in the fluid flow. When the effects of forced flow and free convection/buoyant forces in forced convection are present, a considerable mixed convection flow is produced. Jamaludin et al.23 presented a theoretical model for hybrid nano-fluids with mixed convection boundary layer flow due to an exponentially shrinking/stretching surface in the presence of injection/suction and viscous dissipation under the porous medium. Gohar et al.24 explored the influence of mixed convection flow in Casson hybrid nanofluid under a curved past stretching sheet on the contribution, according to Darcy-Forchheimer, of permeable media in an incompressible viscous fluid flow. Ibrahim and Gizewu25,26,27 examine mixed convection flow with modified Fick’s and Fourier’s diffusion theories of non-Newtonian tangent hyperbolic fluid in the presence of convective and slip boundary conditions past the non-uniform thickness, bi-directional stretching sheet and entropy generation under a thin film flow boundary condition.

Injection or Suction of fluid through the bounding surfaces, as a consequence, affect the rate of heat transfer from the bounding surfaces, as, for example, in mass transfer cooling, it can significantly change the flow field. In general, injection acts in the opposite manners whereas suction tends to increase the skin-friction and heat transfer coefficients . Injection/withdrawal of fluid through porous cooled or heated surface is of general interest in practical problems involving control of boundary layers, film cooling,etc. This can lead to enhanced cooling /heating of the system and can help from laminar flow to delay the transition28,29,30.

Existence of dual solutions in numerical computing have becoming an important topic in fluid dynamics problems for researchers to study its physical significant. Some of the researchers’ previous works on dual solutions are:- Shi et al.31 reported the impacts of the dual solution on Maxwell-fluid flow with mixed convection and thermal radiation in the presence of thermophoresis and Brownian motion past an exponentially shrinking sheet. Yahaya et al.32 evaluated the characteristics of a dual solution of two-dimensional unsteady mixed convection and thermal radiation flow for a Riga-plate with a stagnation point over convective boundary conditions. Nadeem et al.33 investigated the behavior of micro-polar hybrid nano-fluids with boundary layer flow past a shrinking/stretching surface by using the Runge–Kutta Butcher method along with the Nachtsheim–Swigert iteration with the existence of dual solutions. Naramgari et al.34 examined the characteristics of dual solutions in magnetohydrodynamics and mixed convection buoyancy of a non-isothermal stretching sheet in the presence of radiation, injection/suction, chemical reactions, a heat sink/source, and a magnetic field. Mousavi et al.35 and De et al.36 discuss the influence of a dual solution on the magnetic hydrodynamics of MgO–Ag/water hybrid nano-fluids of Casson flow past a shrinking/stretching surface in the presence of radiation, suction, and convective boundary conditions.

The aim of this article is to investigate the stability analysis and dual solutions of mixed convection with thermal radiation flow of Cu–Al\(_{2}O_{3}/H_{2}O\) hybrid nano-fluid over a curved stretching/shrinking surface due to convective and slip conditions in the presence of suction/injection and magnetic field effects. The findings of the present article will provide fruitful results in the future research in the field of fluid dynamics. To the best of the author’s knowledge, these investigation have not been reported by any researcher yet. This article fills this huge gap in the existing literature.

In this article, 2D, incompressible, viscous boundary layer flow of hybrid nanofluid with the stability analysis for the existence of dual solutions of mixed convection with thermal radiation flow of Cu–Al\({_{2}} {{O}}_{3} {{/H}}_{2} {{O}}\) hybrid nano-fluid over a curved stretching/shrinking surface due to convection and slip conditions in the presence of suction/injection are considered. r and s are the coordinates of a radius of curvature R. (u, v) are the velocity components along the s and r directions, and a is the positive constant, where g is the acceleration due to gravity. The magnetic field Bo is taken into account in determining the direction of the flow. A hybrid nano-fluid, in which \(H_{2}O\) is used as the base fluid and Cu and \(Al_{2}O_{3}\) are used as solid nanoparticles, is taken into consideration. Figure 1 presents the physical flow model of the present paper. A system of differential equations under the above assumption for the boundary layers of the flow problem in unsteady form are govern as (see5,7,11,12,16):

Diagram of the present work.

The boundary conditions defined by11,16 are as the follows:

\(V_{w} = -\sqrt{a\nu _{f}}S\) , \(U_{w} = as\) and \(\lambda>0, \lambda< 0, S > 0, S<0, {\lambda _{2}}\) and \({\lambda _{3}}\) are representing stretching sheet, shrinking sheet, suction, injection, velocity slip and thermal jump parameters respectively.

The subsequent similarity variables employed in the assumed steady-state case of the governing problem in Eqs. (1)–(6) given by11 and16 as:

When the similarity transformation [see Eq. (19)] is adopted into Eq. (1) it is sufficed while the steady form of Eqs.(2)–(4) are altered into

where \(Pr = \frac{\nu _{f}}{\alpha _{f}}\) , \(\beta =\frac{\alpha }{a}\) , \(K=R\sqrt{\frac{a}{\nu _{f}}}\) , \(\lambda 1 = \frac{Gr}{(Res)^2}\) , \(Gr =\frac{g\beta _{f}(T_{w}-T_{\infty })s^3}{\nu ^{2}_{f}}\) ,

\(Rd = \frac{16\sigma ^* T^{3}_{\infty }}{3K*K_{f}}\) , \(Ec = \frac{\mu _{f} a^2s^2}{K_{f}(T_{w}-T_{\infty })}\) , \(M = \frac{\sigma Bo^{2}}{a\rho }\) are Prandtl number, parameter of unsteadiness, curvature parameter, mixed convection parameter, local Grashof number, radiation parameter, Eckert number and Hartman number parameter respectively.

When we compare Eqs. (9) and (8), we obtain:

The dimensionless forms of the associated boundary conditions are:

From Eq. (9) P would be written as the following:

The skin friction coefficient \(Cf_{r}\) and local Nusselt number \(Nu_{l}\) (see37) are the physical interests of engineering quantities expressed as follows:

are represented by the reduced local Nusselt’s number, reduced skin friction coefficient, and local Reynolds number, respectively.

The thermo-physical characteristics of nano-particles and base fluid as given by5 in the following Table 1:

The hybrid nano-fluid under consideration are given by5:

Researchers38,39 were reported the existence of dual solutions. To determine these solutions are realizable physically, stability analysis of equations (11)-(13) is required. In order to observe which solution is realizable the parameter \(\tau\) is utilized with an initial value problem where new variables given as the following:

One obtains after substitution of (19) into the unsteady form of (1) and (6) are:

Then, consider the following perturbation functions which is used to check the stability of the steady flow solutions \(g(\xi ) = g_{0}(\xi )\) and \(\theta (\xi ) = \theta _{0}(\xi )\) satisfying the boundary value problem (11)–(13) (see38,39)

where unknown eigenvalue parameter is \(\gamma\) and the functions \(G(\xi ,\tau )\) and \(\theta _{1}(\xi ,\tau )\) are relatively small compared to \(g_{0}(\xi )\) and \(\theta _{0}(\xi )\) respectively. The stability of the solutions are determined by the sign (positive or negative) of the eigenvalues. Putting of (23) into (20)–(22) gives the following linearized eigenvalue problem:

With their corresponding boundary conditions are:

by considering \(\tau = 0\) (24)–(26) are reduced as the following:

With their corresponding boundary conditions are:

It should be reported that for various values of Pr the stability of the steady flow solution \(g_{0}(\xi )\) is determined by the smallest eigenvalue \(\gamma\) . According to38,39, the range of possible eigenvalues can be determined by relaxing a boundary condition on \(G_{0}(\xi )\) . Thus, we relax the condition that \(G_{0}''(\xi ) \rightarrow 0\) as \(\xi \rightarrow \infty\) and for a fixed value of \(\gamma\) (27)–(28) along with the new boundary condition \(G_{0}'''(0) = 1\) has to be solved.

To find the numerical solution of the non-linear ordinary differential equations subject to the boundary conditions of the flow problem, bvp4c solver using Matlab software was employed. To solve the non-linear ordinary differential equations (11) and (12), with their respective boundary conditions (13), we form a system of first order ordinary differential equations by letting in the form of \(g = f1\) ; \(g' = f2\) ; \(g'' = f3\) ;\(g''' = f4\) ; \(g^{iv} = f4'\) ; \(\theta = f5\) ; \(\theta ' = f6\) ; \(\theta '' = f6'\) ; \(gi(a) = fi(a)\) , \(gi(b) = fi(b)\) , \(\theta i(a) = fi(a)\) , \(\theta i(b) = fi(b)\) where \(i = 1,2,3,\cdots 6\) , \(a= 0\) and \(b \rightarrow \infty\) .

Then Eqs. 11–13 are reduced to the following form:

The following values were used to analyze the effects of current governing parameters: \(M= 0.9, Ec = 1, S = 5, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) . The characteristics of stability analysis and dual solutions for the dimensionless parameters of the present study are given as follows: The effect of M; Ec; \(\beta\) ; K; Rd; S; \(\lambda 1\) ; \(\lambda 2\) ; \(\lambda 3\) ; Pr; \(\lambda\) ; \(\phi i\) and \(\phi j\) on the velocity profile, shear stress profile, temperature profile, skin friction coefficient profile, and local Nusselt number profile are shown in Fig. 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24 and 25 and Tables 2, 3 and 4.

Portrays the existence of dual solutions on \(g'(\xi )\) with variation of \(\lambda\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g'(\xi )\) with variation of S when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g'(\xi )\) with variation of\(\lambda 1\) on \(g'(\xi )\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g'(\xi )\) with variation of \(\lambda 2\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g'(\xi )\) with variation of \(\phi j\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of \(\beta\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10\) and \(\phi i = 0.05\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of K when \(M= 0.9, Ec = 1, \xi = 0.4, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of \(\lambda 1\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of \(\lambda 2\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50, Rd = 0.2, S= 5, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of \(\lambda\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of \(\phi i\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of \(\phi j\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10\) and \(\phi i = 0.05\) .

Portrays the existence of dual solutions on \(g''(\xi )\) with variation of S when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of \(\beta\) when \(M= 0.9, Ec = 1, \xi = 0.4, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of Ec when \(M= 0.9, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of K when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of \(\lambda\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2,\phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of M when \(Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of \(\phi j\) when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10\) and \(, \phi i = 0.05\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of Pr when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of Rd when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50, S= 5, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(\theta (\xi )\) with variation of S when \(M= 0.9, Ec = 1, \xi = 0.4, \beta = -1, K= 50,Rd = 0.2, \lambda 1 = 0.1, \lambda 2 = 0.2, \lambda 3 = 0.2, Pr = 6.2, \lambda = -10, \phi i = 0.05\) and \(\phi j = 0.02\) .

Portrays the existence of dual solutions on \(Re_s^{1/2}C_f\) with variation of M.

Portrays the existence of dual solutions on \(Re_s^{-1/2}Nu_s\) with variation of M.

Figures 2, 3, 4, 5 and Fig. 6 depict the velocity profile \(g'(\xi )\) with variations of \(\lambda\) , S, \(\lambda 1\) , and \(\lambda 2\) . Figures 2 and 3 illustrate the velocity profile \(g'(\xi )\) with variations of \(\lambda\) and S, respectively. In these figures, boundary layer thickness of the first solution is greater than the second solution. As the values of both \(\lambda\) and S increase, the velocity profile for the first solution decreases and the opposite trend is observed for the second solution. Figure 4 portrays the velocity profile \(g'(\xi )\) with a variation of \(\lambda 1\) . In this figure, the boundary layer thickness of the first solution is smaller than that of second solution. Physically, this is due to the fact that an enhancement in \(\lambda 1\) increases the buoyancy force, which results in an increase in the velocity profile for both solutions. Figure 5 labels the velocity profile \(g'(\xi )\) with a variation of \(\lambda 2\) . In this figure, the boundary layer thickness of the first solution is grater than the second solution. It is physically noted that the velocity profile for first solution is increasing when \(\lambda 2\) increases. As a result, the opposite trend was observed for both solutions when various values of \(\lambda 2\) were presented. Figure 6 portrays the velocity profile \(g'(\xi )\) with variation of \(\phi j\) .In this figure, the boundary layer thickness of the first solution is greater than the second solution. As the values of \(\phi j\) increase, the velocity profile for the second solution decreases, and the opposite trend is observed for the first solution.

Figures 7, 8, 9, 10, 11, 12, 13 and 14 depict the shear stress profile \(g''(\xi )\) with variations of \(\beta\) , K, \(\lambda 1\) , \(\lambda 2\) , \(\lambda\) , and S. Figure 7 illustrates the shear stress profile \(g''(\xi )\) with a variation of \(\beta\) . In this figure, the boundary layer thickness of first solution is greater than the second solution. It can be observed that the thickness of the shear stress increases or decreases for the second or first solution. As a result, the shear stress profile increased for the second solution, and the opposite trend was observed for the first solution. Figures 8 and 9 delineate the shear stress profile \(g''(\xi )\) with variations of K and \(\lambda 1\) , respectively. It can be observed from both figures that the thickness of the shear stress profile for the first solution is larger than that of the second solution, even though after some interval point of \(\xi\) the opposite trend was observed. As a result, the shear stress profile increased for the first solution, and the opposite trend was observed for the second solution, though after some interval, the inverse of this statement was observed. Figures 10, 11 and 12 label the shear stress profile \(g''(\xi )\) with variations of \(\lambda 2\) ,\(\lambda\) and \(\phi i\) . It can be observed that from all Figs. 10, 11 and 12 the shear stress profile increased for the first solution at some interval point of \(\xi\) , and after this interval, inverse of this statement was observed. Even though the shear stress profile decreases for the second solution in all Figs. 10, 11 and 12. Figures 13 and 14 illustrate the shear stress profile \(g''(\xi )\) with variation of the volume fraction of \(Al_{2}O_{3}\) \(\phi j\) and mass suction parameter S. From Fig. 13 and 14, we came to understand that in both figures, the thickness of the shear stress profile for first solution is smaller than second solution even though after some interval point of \(\xi\) the opposite trend was observed. As a result, the shear stress profile decrease for the first solution , and the opposite trend was observed for the second solution, even though some interval points of \(\xi\) inverse of this statement were observed.

Figures 15, 16, 17, 18, 19, 20, 21, 22 and 23 show the effects of the temperature profile \(\theta (\xi )\) with various values of \(\beta\) , Ec, K, \(\lambda\) , M, \(\phi j\) , Pr, Rd, and S. Figure 15 sketches the temperature profile \(\theta (\xi )\) with a variation of \(\beta\) . It is observed from this figure that the thickness of the thermal boundary layer for the second solution is smaller than the first solution. As the values of \(\beta\) increase, the temperature profile for the first solution increases and the opposite trend is observed for the second solution. Figure 16 sketches the temperature profile \(\theta (\xi )\) with variation of Ec. It is observed from this figure that the thickness of the thermal boundary layer for the second solution is smaller than the first solution. Physically,the ratio of the dynamic temperature to the temperature or the kinetic energy to the enthalpy driving force for heat transfer is Ec. As a result, as the values of Ec increase, the temperature profiles for both solutions increase. Figure 17 sketches the temperature profile \(\theta (\xi )\) with a variation of K. The thickness of the thermal boundary layer becomes thinner as the K impact is increased, even though for the second solution is larger than the first solution. As the values of K increase the temperature profile of both the first and second solutions decreases. Figure 18 depicts the temperature profile \(\theta (\xi )\) with a variation of \(\lambda\) . It is seen from this figure that the thickness of the thermal boundary layer for the second solution is smaller than the first solution. As a result, The second solution is a decrease; the opposite trend was observed. Figure 19 sketches the temperature profile \(\theta (\xi )\) with variations of M. It is seen from this figure that the thermal boundary layer thickness for the first solution is greater than that of the second solution. Physically, the Lorentz force due to the transverse magnetic field has the property of relaxing the fluid velocity and temperature distributions. Accordingly, the temperature distribution of both solutions enlarges as the values of the magnetic parameter M increase. Figure 20 sketches the temperature profile \(\theta (\xi )\) with variation of \(\phi j\) . It can be observed that the thermal boundary layer thickness increases or decreases for the second or first solution. As a result, the temperature profile increased for the first solution, and the opposite trend was observed for the second solution via \(\phi j\) . Figure 21 sketches the temperature profile \(\theta (\xi )\) with variation of Pr. It is seen from this figure that the thermal boundary layer thickness for the first solution is larger than the second solution. physically, the ratio of momentum diffusivity to thermal diffusivity or kinematic viscosity to thermal diffusivity. From this fact, for various values of the Pr, temperature profile have opposite characteristics for both solutions. Figure 22 delineates the temperature profile \(\theta (\xi )\) with variation of Rd. The fact that the relative contribution of conduction heat transfer to thermal radiation transfer is a radiation parameter. As it can be seen from the figure, the thermal boundary layer thickness for the first solution is larger than that of the second solution. This shows that as the rate of radiation increases, the temperature profile in both solutions also increases. Figure 23 sketches the temperature profile \(\theta (\xi )\) with variation of S. It is seen from this figure that the thermal boundary layer thickness for first solution is greater than that of the second solution. As a result, for both solutions, the opposite trend in the suction parameter is increased.

The impact of the dual solution and the critical value for various values of M on the skin friction coefficient and the local Nusselt number are portrayed in Figs. 24 and 25 respectively. Figure 24 portrays the behavior of the dual solution and the critical values for M on the skin friction coefficient. Critical value for M = 12.0, 12.1, and 12.2 are \(\lambda c1 = -1.141414, \lambda c2 = -1.292929\) , and \(\lambda c3 = -1.494949\) respectively, where \(\lambda < 0\) up to which the solution exists. Dual solution exist between \(\lambda ci < \lambda\) where i = 1,2 and 3. For the value of M uphill, of the absolute value of \(\lambda ci\) increase for both the skin friction coefficient and the local Nusselt number. The same property, which can be easily understood is registered in Fig. 25, delineates the characteristic of dual solution and critical value for various value of M on the local Nusselt number.

Table 2 summarizes the current result and compares it with previous work in the open literature for various values of K on -Cfs. The result confirms good agreement, this shows that the present method gives an acceptable result. Tables 3 and 4 show the effect of dual solutions with varying M, S, Br, and Pr values on the \(g''(0)\) , Cfs, and Nus of \(Cu-Al_{2}O_{3}/H_{2}O\) , a hybrid nano-fluid. The results present the properties of the non-dimensional parameters M, S, Br, and Pr on \(g''(0)\) , Cfs, and Nus in an easy-to-understand manner, as shown in Tables 3 and 4. For different estimations of Pr the smallest eigenvalue are shown in Table 5. From this Table 5 we generalize that the first solution is stable (physically realizable) whereas the second solution is unstable (not physically realizable).

In this study, two-dimensional, incompressible, viscous boundary layer flow of hybrid nanofluid with stability analysis and a dual solution of mixed convection and thermal radiation with hybrid nano-fluid over curved surfaces past a stretching/shrinking surface are considered under slip conditions. A hybrid nanofluid,in which water is used as the base fluid and copper and alumina are used as nano-particles, and magnetic field are taken into account. The bvp4c method is used to solve the numerical solution of the governing equations of the present study. The main outputs of this study are:

For the values of both \(\lambda\) and S, the velocity profile for the first solution decreased, and the opposite trend was observed for the second solution.

As the values of \(\phi j\) increase, the velocity profile for the second solution decreases, and the opposite trend is observed for the first solution.

As the values of \(\beta\) increase, the shear stress profile increases for the second solution, and the opposite trend is observed for the first solution.

For the values of K and \(\lambda 1\) , the shear stress profile increased for the first solution, and the opposite trend was observed for the second solution, though after some interval points, the inverse of this statement was observed.

For the values of K, the temperature profiles increase, and for both the first and second solutions, they decrease.

As the values of Ec increase, the temperature profiles for both solutions increase.

For the values of \(S,Pr,Rd,M, \;\text{and}\; \lambda\) , the upwind thermal boundary layer of the first solution is larger than the second solution.

For the value of M uphill, the estimation of the absolute value of \(\lambda ci\) increases for both the skin friction coefficient and the local Nusselt number.

In the second solution, increasing the values of \(\beta\) , Pr, S, Ec, and M has a similar effect on \(g''(0)\) , Cfs, -\(\theta '(0)\) ,and Nus.

In the first solution, increasing the values of Ec, S, and Pr on \(g''(0)\) and Cfs result in a decrease.

The first solution has a positive eigenvalue, whereas the second solution has a negative eigenvalue.

All data generated or analysed during this study are included in this published article.

Thermal conductivity of hybrid nanofluid (\(\hbox {Wm}^{-1}\;\hbox {K}^{-1}\) )

Thermal expansion coefficient of hybrid nanofluid (\(\hbox {K}^{-1}\) )

Dimensionless scale in radial direction (−)

Electrical conductivity of hybrid nanofluid (S/m)

Dynamic viscosity of hybrid nanofluid (\(\hbox {kgm}^{-1}\;\hbox {s}^{-1}\) )

Heat capacity of hybrid nanofluid (\(\hbox {J}/(\hbox {kg}\;\hbox {K}))\) )

Density of the hybrid nanofluid (\(\hbox {kgm}^{-3}\) )

Electrical conductivity of nanofluid (S/m)

Dynamic viscosity of nanofluid (\(\hbox {kgm}^{-1}\;\hbox {s}^{-1}\) )

Heat capacity of nanofluid (\(\hbox {J}/(\hbox {kg}\;\hbox {K}))\) )

Density of the nanofluid (\(\hbox {kgm}^{-3}\) )

Thermal conductivity of nanofluid (\(\hbox {Wm}^{-1}\;\hbox {K}^{-1}\) )

Thermal expansion coefficient of nanofluid (\(\hbox {K}^{-1}\) )

Dynamic viscosity of the fluid (\(\hbox {kgm}^{-1}\;\hbox {s}^{-1}\) )

Specific heat at constant pressure (\(\hbox {J}\,\hbox {kg}^{-1}\;\hbox {K}^{-1}\) )

Density of the fluid (\(\hbox {kgm}^{-3}\) )

Thermal conductivity of the fluid (\(\hbox {Wm}^{-1}\;\hbox {K}^{-1}\) )

Yashkun, U., Zaimi, K., Bakar, N. A. A., Ishak, A. & Pop, I. MHD hybrid nanofluid flow over a permeable stretching/shrinking sheet with thermal radiation effect. Int. J. Numer. Methods Heat Fluid Flow 31(3), 1014–31 (2020).

Waini, I., Ishak, A. & Pop, I. Magnetohydrodynamic flow past a shrinking vertical sheet in a dusty hybrid nanofluid with thermal radiation. Appl. Math. Mech. 1–14 (2022).

Hayat, T., Bilal Ashraf, M., Shehzad, S. & Alsaedi, A. Mixed convection flow of Casson nanofluid over a stretching sheet with convectively heated chemical reaction and heat source/sink. J. Appl. Fluid Mech. 8(4), 803–13 (2015).

Jamaludin, A., Nazar, R. & Pop, I. Mixed convection stagnation-point flow of a nanofluid past a permeable stretching/shrinking sheet in the presence of thermal radiation and heat source/sink. Energies 12(5), 788 (2019).

Yashkun, U., Zaimi, K., Ishak, A., Pop, I. & Sidaoui, R. Hybrid nanofluid flow through an exponentially stretching/shrinking sheet with mixed convection and Joule heating. Int. J. Numer. Methods Heat Fluid Flow 31(6), 1930–50 (2021).

Wahid, N. S. et al. Unsteady MHD mixed convection flow of a hybrid nanofluid with thermal radiation and convective boundary condition. Chin. J. Phys. 77, 378–92 (2022).

Article  MathSciNet  CAS  Google Scholar 

Muhammad, K., Hayat, T., Alsaedi, A., Ahmad, B. & Momani, S. Mixed convective slip flow of hybrid nanofluid (MWCNTs+ Cu+ Water), nanofluid (MWCNTs+ Water) and base fluid (Water): A comparative investigation. J. Therm. Anal. Calorim. 143, 1523–36 (2021).

Khan, U. et al. Stagnation point flow of a water-based graphene-oxide over a stretching/shrinking sheet under an induced magnetic field with homogeneous-heterogeneous chemical reaction. J. Magn. Magn. Mater. 565, 170287 (2023).

Khan, U. et al. Features of hybridized AA7072 and AA7075 alloys nanomaterials with melting heat transfer past a movable cylinder with Thompson and Troian slip effect. Arab. J. Chem. 16(2), 104503 (2023).

Xu, H. Mixed convective flow of a hybrid nanofluid between two parallel inclined plates under wall-slip condition. Appl. Math. Mech. 43(1), 113–26 (2022).

Article  MathSciNet  MATH  Google Scholar 

Roşca, N. C. & Pop, I. Unsteady boundary layer flow over a permeable curved stretching/shrinking surface. Eur. J. Mech. B Fluids 51, 61–7 (2015).

Article  ADS  MathSciNet  MATH  Google Scholar 

Revathi, G. et al. Entropy optimization in hybrid radiative nanofluid (CH 3 OH+ SiO 2+ Al 2 O 3) flow by a curved stretching sheet with cross-diffusion effects. Appl. Nanosci. 1–15 (2021).

Muhammad, R., Khan, M. I., Jameel, M. & Khan, N. B. Fully developed Darcy–Forchheimer mixed convective flow over a curved surface with activation energy and entropy generation. Comput. Methods Progr. Biomed. 188, 105298 (2020).

Ahmed, N. et al. Nonlinear thermal radiation and chemical reaction effects on a (Cu-CuO)/NaAlg hybrid nanofluid flow past a stretching curved surface. Processes 7(12), 962 (2019).

Ijaz Khan, M., Khan, S. A., Hayat, T., Qayyum, S. & Alsaedi, A. Entropy generation analysis in MHD flow of viscous fluid by a curved stretching surface with cubic autocatalysis chemical reaction. Eur. Phys. J. Plus 135(2), 1–17 (2020).

Ibrahim, W. & Gizewu, T. Analysis of entropy generation of bio-convective on curved stretching surface with gyrotactic micro-organisms and third order slip flow. Int. J. Thermofluids 100277 (2023).

Wahid, N. S. et al. MHD mixed convection flow of a hybrid nanofluid past a permeable vertical flat plate with thermal radiation effect. Alex. Eng. J. 61(4), 3323–33 (2022).

Bejawada, S. G. & Nandeppanavar, M. M. Effect of thermal radiation on magnetohydrodynamics heat transfer micropolar fluid flow over a vertical moving porous plate. Exp. Comput. Multiphase Flow 5(2), 149–58 (2023).

Lv, Y. P. et al. Chemical reaction and thermal radiation impact on a nanofluid flow in a rotating channel with Hall current. Sci. Rep. 11(1), 19747 (2021).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Hussain, S. M., Sharma, R., Seth, G. S. & Mishra, M. R. Thermal radiation impact on boundary layer dissipative flow of magneto-nanofluid over an exponentially stretching sheet. Int. J. Heat Technol. 36(4), 1163–73 (2018).

Krishnamurthy, M., Gireesha, B., Prasannakumara, B. & Gorla, R. S. R. Thermal radiation and chemical reaction effects on boundary layer slip flow and melting heat transfer of nanofluid induced by a nonlinear stretching sheet. Nonlinear Eng. 5(3), 147–59 (2016).

Saleem, S. et al. Heat transfer in a permeable cavity filled with a ferrofluid under electric force and radiation effects. AIP Adv 9(9) (2019).

Jamaludin, A., Nazar, R., Naganthran, K. & Pop, I. Mixed convection hybrid nanofluid flow over an exponentially accelerating surface in a porous media. Neural Comput. Appl. 33(22), 15719–29 (2021).

Gohar, S. K. T., Khan, I., Gul, T. & Bilal, M. Mixed convection and thermally radiative hybrid nanofluid flow over a curved surface. Adv. Mech. Eng.https://doi.org/10.1177/1687813222108284 (2022).

Ibrahim, W. & Gizewu, T. Nonlinear mixed convection flow of a tangent hyperbolic fluid with activation energy. Heat Transf. 49(5), 2427–48 (2020).

Ibrahim, W. & Gizewu, T. Tangent hyperbolic nanofluid with mixed convection flow: An application of improved Fourier and Fick’s diffusion model. Heat Transf. Asian Res. 48(8), 4217–39 (2019).

Ibrahim, W. & Gizewu, T. Thin film flow of tangent hyperbolic fluid with nonlinear mixed convection flow and entropy generation. Math. Probl. Eng. 2021, 1–16 (2021).

Chaudhary, M. & Merkin, J. The effects of blowing and suction on free convection boundary layers on vertical surfaces with prescribed heat flux. J. Eng. Math. 27(3), 265–92 (1993).

Article  MathSciNet  MATH  Google Scholar 

Hussain, M., Ali, A., Inc, M., Sene, N. & Hussan, M. Impacts of chemical reaction and suction/injection on the mixed convective Williamson fluid past a penetrable porous wedge. J. Math. (2022).

Uwanta, I. & Hamza, M. Effect of suction/injection on unsteady hydromagnetic convective flow of reactive viscous fluid between vertical porous plates with thermal diffusion. Int. Sch. Res. Not. 2, 14 (2014).

Shi, Q. H. et al. Dual solution framework for mixed convection flow of Maxwell nanofluid instigated by exponentially shrinking surface with thermal radiation. Sci. Rep. 11(1), 15944 (2021).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Yahaya, R. I., Arifin, N. M., Pop, I., Ali, F. M. & Isa, S. S. P. M. Dual solutions of unsteady mixed convection hybrid nanofluid flow past a vertical Riga plate with radiation effect. Mathematics. 11(1), 215 (2023).

Nadeem, S., Israr-ur Rehman, M., Saleem, S. & Bonyah, E. Dual solutions in MHD stagnation point flow of nanofluid induced by porous stretching/shrinking sheet with anisotropic slip. AIP Adv. 10(6), 065207 (2020).

Article  ADS  CAS  Google Scholar 

Naramgari, S. & Sulochana, C. Dual solutions for MHD stagnation-point flow of a nanofluid over a stretching surface with induced magneticfield. Int. J. Sci. Eng. 9(1), 1–8 (2015).

Mousavi, S. M. et al. Dual solutions for Casson hybrid nanofluid flow due to a stretching/shrinking sheet: A new combination of theoretical and experimental models. Chin. J. Phys. 71, 574–88 (2021).

Article  MathSciNet  CAS  Google Scholar 

De, P., Mondal, H. & Bera, U. K. Dual solutions of heat and mass transfer of nanofluid over a stretching/shrinking sheet with thermal radiation. Meccanica 51, 117–24 (2016).

Article  MathSciNet  MATH  Google Scholar 

Wahid, N. S. et al. Flow and heat transfer of hybrid nanofluid induced by an exponentially stretching/shrinking curved surface. Case Stud. Therm. Eng. 25, 100982 (2021).

Roşca, A. V. & Pop, I. Flow and heat transfer over a vertical permeable stretching/shrinking sheet with a second order slip. Int. J. Heat Mass Transf. 60, 355–64 (2013).

Weidman, P., Kubitschek, D. & Davis, A. The effect of transpiration on self-similar boundary layer flow over moving surfaces. Int. J. Eng. Sci. 44(11–12), 730–7 (2006).

Article  CAS  MATH  Google Scholar 

Department of Mathematics, Ambo University, Ambo, Ethiopia

Wubshet Ibrahim & Tezera Gizeu

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

Both authors contributed to the study conception and design. Derivation of governing equations, Numerical computations, Matlab code writing and analysis were performed by both authors. The first draft of the manuscript was written by T.G. and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Ibrahim, W., Gizewu, T. Stability analysis of dual solutions for mixed convection and thermal radiation with hybrid nanofluid flow past shrinking/stretching curved surface. Sci Rep 13, 21676 (2023). https://doi.org/10.1038/s41598-023-48728-8

DOI: https://doi.org/10.1038/s41598-023-48728-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Stability analysis of dual solutions for mixed convection and thermal radiation with hybrid nanofluid flow past shrinking/stretching curved surface | Scientific Reports

Electric Heating Tube Sign up for the Nature Briefing: AI and Robotics newsletter — what matters in AI and robotics research, free to your inbox weekly.