Non-similar analysis of heat generation and thermal radiation on Eyring-Powell Hybrid nanofluid flow across a stretching surface

Abstract

This paper portrays a two-dimensional mathematical model developed for the analysis of Eyring-Powell hybrid nanofluid flow, incorporating variable temperature and velocity profiles. This study investigates the flow through a porous media with a mixture of hybrid nanoparticles (copper dioxide, magnetite) integrated into ethylene glycol ( $EG$ ) across a stretching sheet. The model takes magnetic effects, heat generation, and thermal radiation into account. Additionally, nonsimilar partial differential equations (PDEs) are transformed into ordinary differential equations using the local nonsimilarity approach. These equations are then numerically solved using the built-in Matlab function bvp4c. To examine the effects of adjusting parameters on heat, mass transfer, and skin friction, the results are displayed graphically. With an increase in the given values of magnetic field $M$ , the velocity profile $f' (ξ, η)$ is seen to decay. Moreover, it is observed that the hybrid nanofluid density and dynamic viscosity increase as the volume fraction rises. It is estimated that the thermal conductivity and viscosity of the ${Fe}_{3} O 4 + Cu / EG$ hybrid nanofluid are expected to increase to 4.1% and 12.35%, respectively, and to 71.55% and 78.01%, respectively, for volume concentrations ranging from 0.02% to 0.05%. It is also concluded that the local skin-friction coefficient has the opposite trend while the local Nusselt number rises monotonically with both the slip velocity parameter and the surface convection parameter. This research has broad applications In the glass and polymer industries, metallic plate cooling, plastic sheet contractions, etc.

Keywords

Eyring-Powell hybrid nanofluid heat generation thermal radiation porous medium magnetic field

Introduction

Nanofluid is a kind of sophisticated manufactured fluid that is made up of a base fluid (usually a liquid, such ethylene glycol, water, or oil) suspended in nanoparticles, which particles smaller than a nanometer. Frequently, metals, metal oxides, carbon nanotubes, and other nanomaterials are used to create these nanoparticles. Engineering, physics, materials science, and other disciplines are interested in nanofluids because they have special and improved features that result from the addition of nanoparticles to the basic fluid. Thermal conductivity of nanofluids is typically much higher than that of the base fluid alone. Because of this characteristic, they can be applied to thermal management and heat transfer systems, including the creation of more effective heat exchangers and electronic cooling systems. Shafiq and Sindhu¹ conducted a statistical analysis of Williamson nanofluid on a porous stretching surface under the influence of magnetohydrodynamic (MHD) forces. Shekaramiz et al.² investigated the MHD boundary layer flow of a nanofluid over a stretched sheet. An investigation into Maxwell nanofluid flow with convective boundary conditions was carried out by Razzaq et al.³ Shah et al.⁴ considered heat radiation and viscous dissipation while examining the flow of nanofluid through a stretchable Riga wall. Khan et al.⁵ study was concerned with the Casson nanofluid flow above a nonlinear stretching plate. Farooq et al.⁶ examined the potential of MHD bioconvective micropolar nanofluids while taking into account the soret and dufour impacts. Ali et al.⁷ numerically investigated the transient MHD 3D rotating flow of Maxwell Cattaneo Christov heat flux model over a bidirectional stretching sheet. Khattak et al.⁸ explored the study of nanofluid flow over a porous and deformable plate, taking into account Joule heating, thermal radiation, as well as thermophoresis and Brownian motion effects within the framework of the Cattaeno–Christov heat flux model. Razzaq et al.⁹ examined the nanofluid flow for single-multilayer carbon nanotubes upon a vertical porous cone while considering the influences of MHD and thermal radiation.

Hybrid nanofluids represent a compelling solution for heat storage, exhibiting notably higher energy efficacious contrast to conventional nanofluids. The utilization of hybrid nanofluids, comprising manganese, zinc ferrite, copper, and silver nanocomposites on spinning discs has garnered increasing attention in both research and development due to their multifaceted applications. The MHD convection of $(A l_{2} O_{3} - Cu / Water$ ) with an internal heat sink/source in a porous media was investigated by Chamkha et al.¹⁰ It is concluded that raising the number of nanoparticles, the Rayleigh numbers reduce thermal reaction and increases heat transmission in the case of high-volume percentage. Cao et al.¹¹ investigated the effects of various nanomaterials under distinct degrees of partial slip using ternary-hybrid nanofluid models. Applying the concept of non-Fourier principles, Algehyne et al.¹² examined the flow of a fluid over a Riga plate, considering slip impacts and activation energy, to explore the influence of MHD on the movement of $Ag - MgO$ / $H_{2} O$ -embedded hybrid nanofluid containing motile microorganisms. Sarada et al.¹³ developed a heat flux model utilizing non-Fourier analysis to delve into the heat production implications within the hybridization flow of ternary nanofluids. Elnaqeeb et al.¹⁴ developed three-dimensional flow involving dual stretching and suction of ternary hybrid nano-fluid within a closed rectangular region, with various shapes and densities. Khazayinejad and Nourazar¹⁵ studied spatial fractional energy transmission and continuous laminar MHD boundary layer flow using a $Gr - F e_{3} O_{4} - H_{2} O$ hybrid nanofluid. The flow with heat transfer containing Eyring-Powell nanoparticles from a moving stretched/shrinking sheet was analyzed by Roşca and Pop.¹⁶ El-Zahar et al.¹⁷ used a radioactive rolling cylinder to study the mixed convection flow of magneto-hybrid nanofluids. It was discovered that the rigid volume fraction was lessened by the force of drag.

Within the field of fluid dynamics, materials’ behavior is divided into two categories: Newtonian and non-Newtonian fluids. Newton’s law of viscosity governs Newtonian fluids like water and air, but non-Newtonian fluids behave more complexly, when viscosity varying with time, stress, and shear rate, among other factors. Eyring-Powell nanofluid model is an example of non-Newtonian fluids. It is applied to explain the rheological behavior of colloidal suspensions of nanoparticles in a base fluid. Cui et al.¹⁸ used nonsimilar modeling to study the forced convection of magnetized Eyring-Powell nanofluid flow in a porous media above a Riga sheet. Mushtaq et al.¹⁹ examined the effects of an exponentially growing sheet under the assumption of Eyring-Powell fluid. They examined the opposite impacts of viscoelastic properties on heat-velocity boundary layer. Rahimi et al.²⁰ utilized the collocation technique to estimate the boundary layer of the Powell-Eyring fluid. Akbar et al.²¹ presented numerical solutions for the Powell-Eyring magnetofluid above an expanding sheet. Hassan et al.²² conducted research focusing on examining the flow characteristics of the Powell-Eyring fluid and its associated entropy generation concerning time scale and viscosity parameters. Ghadikolaei et al.²³ showcased the time-dependent attribute of a MHD flow of Eyring-Powell fluid constrained by a channel, influenced by Joule dissipation.

Thermal radiation is a ubiquitous phenomenon that plays a vital role in numerous natural and technological processes. It is invaluable in fields ranging from engineering and physics to astronomy and environmental science because of its distinctive qualities and applications. Determining the properties of thermal radiation is crucial for developing effective heat transfer systems, utilizing renewable energy sources, and solving cosmic mysteries. Jan et al.²⁴ studied the influences of thermal generation-absorption and applied magnetic field on Sisko fluid model. Cui et al.²⁵ used a model of Casson fluid with radiative nanofluids serving as the magnetized boundary layer to simulate a stretched sheet of fluid. Abrar et al.²⁶ explored the entropy generation within cilia-driven nanofluid transport, considering the impacts of thermal radiation, internal heating, and viscous dissipation. Kumar et al.²⁷ developed a computational model for the radiative heat transfer in the presence of magnetic field influences and viscous dissipation over an infinite vertical sheet. In addition, the author highlights the application of the under-study problem in the chemical and metallurgy sector. Agbaje et al.²⁸ performed a computational study to determine the unstable flux in the presence of thermal radiation from Powell-Eyring nanofluid and heat obstetrics. The findings demonstrated that the fluid’s temperature and the thermal boundary layer’s thickness were elevated by the characteristics of the heat source and thermal radiation. Oreyeni et al.²⁹ made analysis on mixed hybrid magnetite particles ( $F e_{3} O_{4}$ ) with silver (Ag) to investigate the effects of heat radiation, magnetic fields, and viscous dissipation on the thermal characterization function in solar-powered ships. The analysis methodology applied was the Galerkin weighted residual method, and thermal boundary layer’s heat flow was assessed utilizing the Cattaneo-Christov model.

The method of local non-similarity was proposed by Sparrow and Yu³⁰ as a solution to the serious flaws and errors in the similarity method. One frequent issue with using a local similarity method is that the system is not fully converted into ordinary differential equations (ODEs). In comparison to the conventional similarity strategy, this yields a more consistent and accurate result. The local non-similarity approach addresses this constraint by considering non-similar transformations in the problem’s immediate region. In a non-Darcy porous medium, Farooq et al.³¹ used a non-similar boundary layer model to study forced convection in the steady flow of Casson nanofluid. Similarly, non-similar solution for Maxwell hybrid nanofluid in the presence of MHD discussed by Sagheer et al.³² The MHD convective transport of a nanofluid across a stretchy surface was investigated by Riaz et al.³³ using the bvp4c approach, accounting for several effects. Using a non-similarity technique, Sagheer et al.³⁴ investigated increased thermal utilization for EMHD nanofluid with variable thermal flux along the stretching surface.

This study focuses on the selection of copper oxide $(Cu),$ magnetite $(F e_{3} O_{4})$ , and ethylene glycol $(EG)$ based on their distinct thermal characteristics. Metallic oxide nanoparticles $F e_{3} O_{4}$ and $Cu$ are easily soluble in the base fluid, while $EG$ exhibits strong thermal features. Pure metallic nanoparticles $(Cu$ , $F e_{3} O_{4}) and EG$ are recognized for their larger aspect ratio, superior heat conductivity, and low specific gravity. Although nanofluids tend to be unstable over short durations, the addition of $F e_{3} O_{4} + Cu$ / $EG$ ensures long-term stability and enhances thermal properties significantly. According to existing literature, it is observed that the influence of viscous dissipation on MHD heat and mass flow of hybrid-nano fluid ( $Cu$ + $F e_{3} O_{4}$ ) across a stretched sheet through a permeable media using non similarity transformation remains unexplored. Utilizing the local non similarity method alongside the bvp4c algorithm up to its second truncation level enables the computation of numerical solutions to nonlinear-couple partial differential equations. Through adjustment of pertinent parameters, relevant findings on the subject matter are elucidated and visually presented using tables and diagrams.

Mathematical modeling of the problem

The following assumptions are performed under the considered flow.

• Considering a steady, two-dimensional flow of Eyring-Powell nanofluid above a linearly stretching sheet, characterized by a velocity profile of $u = xc$ .

• The Cartesian coordinate system aligns the x-axis with the extending surface, and the y-axis is normal to it.

• The stretching surface undergoes expansion and is subjected to reheating, starting at an initial temperature $T_{w}$ , while $T_{\infty}$ represents the ambient temperature.

• The flow is generated by the action of two equal and opposite forces along the x-axis and the sheet is assumed to be linearly stretching.

• The velocity components are denoted as $v$ and $u$ in the $y$ and $x$ directions, respectively. Figure 1 illustrates the flow coordinate system and provides a detailed description of the model assumptions involved in the setup.

Figure 1.

Flow configuration.

The Cauchy stress tensor for Eyring-Powell fluid can be expressed as:

A = - PI + τ,

(1)

τ = [μ + \frac{1}{δ γ} \sin h^{- 1} (\frac{1}{m_{1}} γ)] B_{1},

(2)

In equation (1), the term “ $I$ ” represents the identity matrix, and “ $P$ ” signifies the pressure component. The symbol τ denotes the extra stress tensor, where $m_{1}$ and δ represent fluid parameters, μ denotes the dynamic viscosity, with $m_{1}$ having dimensions of $(time)^{- 1} .$ Additionally, $B_{1}$ represents the first Rivlin-Erickson tensor and is expressed as follows:

B_{1} = gradV + grad V^{T} and γ = \sqrt{\frac{1}{2} tr {B_{1}}^{2}}

(3)

Using the second-order approximations of the $\sinh - 1$ function and ignoring the third and so on terms² we obtain

\sin h^{- 1} (\frac{1}{m_{1}} γ) \approx \frac{γ}{m_{1}} - \frac{1}{6} {(\frac{γ}{m_{1}})}^{3}, with \frac{1}{6} {(\frac{γ}{m_{1}})}^{5} < < < 1 .

(4)

The governing equations related to flow problem are³⁵:

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

(5)

\begin{matrix} u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = (ν_{hnf} + \frac{1}{δ m_{1} ρ_{hnf}}) \frac{\partial^{2} u}{\partial y^{2}} - \frac{1}{2 δ ρ_{hnf} {m_{1}}^{3}} \frac{\partial^{2} u}{\partial y^{2}} {(\frac{\partial u}{\partial y})}^{2} \\ - \frac{σ_{hnf} B_{o}^{2}}{hnf} u - \frac{μ_{hnf}}{ρ_{hnf}} \frac{u}{K_{1}}, \end{matrix}

(6)

\begin{matrix} u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = \frac{k_{hnf}}{{(ρ C_{p})}_{hnf}} (\frac{\partial^{2} T}{\partial y^{2}}) - \frac{1}{{(ρ C_{p})}_{hnf}} \frac{\partial q_{r}}{\partial y} \\ + \frac{Q_{0}}{{(ρ C_{p})}_{hnf}} (T - T_{\infty}), \end{matrix}

(7)

The boundary conditions for a given problem are:

u = cx, v = 0, T = T_{\infty} at y = 0, as y \to u = 0, T = T_{\infty} .

(8)

The radiative heat flux $q_{r}$ is defined as :

q_{r} = - \frac{4 σ °}{3 k °} ({\frac{\partial T}{\partial y}}^{4}) .

(9)

Utilization of Taylor series expansion:

T^{4} = T_{\infty}^{4} + 4 T_{\infty}^{3} (T - T_{\infty}) + 6 T_{\infty}^{2} {(T - T_{\infty})}^{2} + \dots

Substitution of equation (10) in equation (7), we obtain equation (11) as given below:

q_{r} = - \frac{16 T_{\infty}^{3} σ °}{3 k °} (\frac{\partial T}{\partial y}),

(10)

\begin{matrix} u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = \frac{k_{hnf}}{{(ρ C_{p})}_{hnf}} (\frac{\partial^{2} T}{\partial y^{2}}) - \frac{1}{{(ρ C_{p})}_{hnf}} \\ \times \frac{16 T_{\infty}^{3} σ °}{3 K °} \frac{\partial^{2} T}{\partial y^{2}} + \frac{Q_{o}}{ρ C_{p}} (T - T_{\infty}) . \end{matrix}

(11)

Nonsimilar system

In the context of nonsimilar flow, we will introduce new parameters denoted as $ξ (x)$ for nonsimilarity variable, and $η (x, y)$ to represent a pseudo-similarity variable.^36,37

\begin{matrix} ξ = \frac{x}{l}, η = \sqrt{\frac{a}{ν}} y, u = ax \frac{\partial f}{\partial η} (ξ, η), \\ v = - \sqrt{a ν} (ξ \frac{\partial f}{\partial ξ} + f (ξ, η)), θ (ξ, η) = \frac{T - T_{\infty}}{(T_{w} - T_{\infty})} \end{matrix}

(12)

Putting equation (12) in equations (6)–(8), we get the following dimensionless system.

\begin{matrix} \frac{\partial^{3} f}{\partial η^{3}} (B_{1} ε + B_{2} - ξ^{2} ε ω {(\frac{\partial^{2} f}{\partial η^{2}})}^{2}) - M B_{3} \frac{\partial f}{\partial η} - B_{4} K \frac{\partial f}{\partial η} \\ - {(\frac{\partial f}{\partial η})}^{2} + f \frac{\partial^{2} f}{\partial η^{2}} = ξ (\frac{\partial f}{\partial η} \frac{\partial^{2} f}{\partial ξ \partial η} - \frac{\partial^{2} f}{\partial η^{2}} \frac{\partial f}{\partial ξ}), \end{matrix}

(13)

(B_{5} + \frac{4}{3} R B_{6}) \frac{\partial^{2} θ}{\partial η^{2}} + \Pr (f \frac{\partial θ}{\partial η} + B_{6} Q θ) = \Pr ξ (\frac{\partial f}{\partial η} \frac{\partial θ}{\partial ξ} - \frac{\partial θ}{\partial η} \frac{\partial f}{\partial ξ}) .

(14)

Subjected BCs;

\begin{matrix} \frac{\partial f (0, ξ)}{\partial η} = 1, f (0, ξ) = 0, θ (ξ, 0) = 1, at y = 0, \frac{\partial f (ξ, \infty)}{\partial η} \\ = 0, θ (ξ, \infty) = 0, at y \to \infty \end{matrix}

\begin{array}{l} B_{1} = \frac{1}{[(1 - ϕ_{2}) {(1 - ϕ_{1}) + ϕ_{1} ρ_{s 1}} + ϕ_{2} ρ_{s 2}]}, \\ B_{2} = [\frac{1}{{(1 - ϕ_{2})}^{2.5} {(1 - ϕ_{1})}^{2.5}}], \\ B_{3} = [1 + \frac{3 (S_{f} - 1) (σ_{s 2} - 1) ϕ_{2}}{(σ_{s 2} + 2) - (σ_{s 2} - 1) (S_{f} - 1) ϕ_{2}} \\ + \frac{3 (S_{f} - 1) (σ_{s 1} - 1) ϕ_{1}}{(σ_{s 1} + 2) - (σ_{s 1} - 1) (S_{f} - 1) ϕ_{1}}], \\ B_{4} = [(1 - ϕ_{2}) {(1 - ϕ_{1}) + ϕ_{1} ρ_{s 1}} + ϕ_{2} ρ_{s 2}] \\ \frac{1}{{(1 - ϕ_{2})}^{2.5} {(1 - ϕ_{1})}^{2.5}}, \\ B_{5} = [(1 - ϕ_{2}) {(1 - ϕ_{1}) + ϕ_{1} ρ_{s 1}} + ϕ_{2} ρ_{s 2}] \\ \times \frac{k_{s 2} + (s f - 1) k_{b f} - (s f - 1) ϕ_{2} (k_{b f} - k_{s 2})}{k_{s 2} + (s f - 1) k_{b f} + ϕ_{2} (k_{b f} - k_{s 1})} \cdot \\ \frac{k_{s 1} + (s f - 1) k_{f} - (s f - 1) ϕ_{1} (k_{f} - k_{s 1})}{k_{s 1} + (s f - 1) k_{f} + ϕ_{1} (k_{f} - k_{s 1})}, \\ B_{6} = (1 - ϕ_{2}) {ϕ_{1} {(ρ C_{p})}_{s 1} + (1 - ϕ_{1}) {(ρ C_{p})}_{f}} \\ + ϕ_{2} {(ρ C_{p})}_{s 2}, A_{1} = \frac{μ_{h n f}}{μ_{f}}, A_{2} = \frac{k_{h n f}}{k_{f}} \end{array}

(15)

Now the magnetic parameter is expressed as $M = \frac{σ_{f} B_{o}^{2}}{ρ_{fc}}, R = \frac{16 T_{\infty}^{3} σ °}{3 k_{f} K °}$ is thermal radiation parameter, $\Pr = \frac{ν_{f} {(ρ c_{p})}_{f}}{k_{f}}$ is Prandtl number, $w = \frac{c^{3} l^{2}}{2 {m_{1}}^{2} ν_{f}} and ε = \frac{1}{m_{1} ν_{f} ρ_{f δ}}$ represents the material constants of Eyring–Powell fluid. The permeability of porous medium is denoted by $K = \frac{ν_{f}}{c K_{1}}$ and $Q = \frac{Q_{o}}{c ρ_{f} {C_{p}}_{f}}$ represents thermal sink/source. The skin friction and Nusselt number expressed as:

\begin{matrix} R e^{\frac{1}{2}} C_{f} = A_{1} + ε f ″ (ξ, 0) - \frac{ε ω}{3} ξ^{2} f^{''^{3}} (ξ, 0), \\ R e^{\frac{- 1}{2}} Nu = (A_{1} + \frac{4}{3} R) θ' (ξ, 0) . \end{matrix}

(16)

Where $Re = \frac{c l^{2}}{ν_{f}} .$

The specific assumption for above stated physical parameters are settled as

M = K = 0.5, ε = 0.1, w = 0.2, ϕ_{1} = 0.02, ϕ_{1} = 0.04

When dealing with nonsimilar boundary layer problems, local nonsimilarity is an often-employed method. Before moving on to local nonsimilar solutions, it is highly beneficial to assess the local similarity technique. This method presupposes a boundary layer equation structure with associated conditions and determines that the right-hand side, denoted as term $(\partial (.)) / \partial ξ)$ of equations (13)–(15) to be adequately small, thereby allowing it to be approximated as zero.

\begin{matrix} \frac{\partial^{3} f}{\partial η^{3}} (B_{1} ε + B_{2} - ξ^{2} ε ω {(\frac{\partial^{2} f}{\partial η^{2}})}^{2}) - M B_{3} \frac{\partial f}{\partial η} - B_{4} K \frac{\partial f}{\partial η} \\ - {(\frac{\partial f}{\partial η})}^{2} + f \frac{\partial^{2} f}{\partial η^{2}} = 0, \end{matrix}

(17)

(B_{5} + \frac{4}{3} R B_{6}) \frac{\partial^{2} θ}{\partial η^{2}} + \Pr (f \frac{\partial θ}{\partial η} + B_{6} Q θ) = 0 .

(18)

\begin{matrix} \frac{\partial f (ξ, 0)}{\partial η} = 1, f (ξ, 0) = 0, \frac{\partial θ (ξ, 0)}{\partial η} = 0, \\ \partial f (ξ, \infty) / \partial η \to 0, θ (ξ, \infty) \to 0 \end{matrix}

(19)

The value of $ξ$ can be regarded as a constant parameter at any stream-wise point. Therefore, equations (19)–(21) mentioned above can be regarded as ordinary differential equations. It would be mentioned that some results and accuracy are questionable when terms $(\partial (.) / \partial ξ)$ are not comprehended. In order to get around these difficulties, Sparrow and Yu³⁰ and Minkowycz and Sparrow³⁸ suggested local non-similarity technique to determine the implications of the non-similar boundary layer equations. Now let’s introduce functions using a local non-similarity approach as.

\begin{matrix} \frac{\partial f}{\partial η} = f', \frac{\partial^{2} f}{\partial η^{2}} = f ″, \frac{\partial^{3} f}{\partial η^{3}} = f ″', \frac{\partial θ}{\partial η} = θ', \frac{\partial^{2} θ}{\partial η^{2}} = θ ″, \\ g (ξ, η) = \frac{\partial f (ξ, η)}{\partial ξ}, p (ξ, η) = \frac{\partial θ (ξ, η)}{\partial ξ} . \end{matrix}

(20)

Now by using (20) in (17)–(19) the equations are transformed as

\begin{matrix} f ″' (B_{1} ε + B_{2} - ξ^{2} ε ω {f ″}^{2}) + M B_{3} f' - B_{4} K f' - {f ″}^{2} + f f ″ \\ = ξ ({f'}^{2} g - f ″ g) . \end{matrix}

(21)

(B_{5} + \frac{4}{3} R B_{6}) θ ″ + \Pr (B_{6} Q θ + ff') = \Pr ξ (f' p - θ' g) .

(22)

f' (ξ, 0) = 1, f (ξ, 0) = 0, θ' (ξ, 0) = 0 .

f' (ξ, \infty) \to 0, θ (ξ, \infty) \to 0 .

(23)

We get (24)–(26) by taking partial derivatives of (21)–(23) with respect to $ξ$ :

\begin{matrix} g ″' (B_{1} ε + B_{2} - ξ^{2} ε ω {f ″}^{2}) - 2 ξ ε ω f ″' ({f ″}^{2} + ξ g ″ f ″) \\ - B_{3} M g' - B_{4} K g' + g ″ + 2 g f ″ + 2 f ″ g + g ″ f \\ - gf''^{2} = ξ g (2 f' g' - g ″) . \end{matrix}

(24)

\begin{matrix} (B_{5} + \frac{4}{3} R B_{6}) p ″ + \Pr (B_{6} Qp + f p' + g θ' - f' p + p' g) \\ = \Pr ξ (g' p - p' g) . \end{matrix}

(25)

g (ξ, 0) = 0, g' (ξ, 0) = 0, p' (ξ, 0) = 0, g' (ξ, \infty) \to 0, p (ξ, \infty) \to 0 .

(26)

Numerical solution

The precise solution to the overall problem presents significant challenges due to the highly non-linear nature of the governing equations, as addressed by the local non-similar method. To tackle this, the local non-similar equations are initially reformulated into a system of first-order differential equations by introducing new variables. These reformulated equations are then processed using the bvp4c solver, which employs the Lobatto IIIa formula-a three-stage method. This MATLAB routine ensures a solution with accuracy up to the third order over the entire interval. Consequently, the proposed problem must be translated into a system of first-order equations to facilitate the approximation of the solution, following the outlined procedure

\begin{matrix} (f = f (1), f' = f (2), f ″ = f (3)), (θ = f (4), θ' = f (5)) \\ (g = f (6), g' = f (7), g ″ = f (8)), (p = f (9), P' = f (10)) \end{matrix}

Result and discussion

This section illustrates the Cu and $F e_{3} O_{4}$ hybrid nanofluid flow embedded in ethylene glycol with expanding surfaces. Figures and tables illustrate the significance of flow limitations in relation to the energy and velocity profiles. Table 3 showed excellent consistency between the current investigation and the literature under limited settings.^39,40 Numerical solutions for the examined fluid flow issues are generated using the MATLAB-based tool bvp4c, which considers physical factors in both scenarios, particularly $F e_{3} O_{4}$ /EG (nanofluid) and $F e_{3} O_{4} + Cu$ /EG used as “hybrid-nanofluid.” In Table 1 thermo-physical properties for hybrid nanofluids and nanofluids are discussed by Rashad et al.³⁵ and Table 2 represents the thermal proprieties for Fe₃O₄, Cu, and EG.^36,37 The result is verified by contrasting it with the body of knowledge, as indicated in Table 3 below. Table 4 shows the impact of governing factors on skin friction factor and temperature gradient at boundary. Computations are performed for $F e_{3} O_{4}$ − $Cu / EG$ hybrid nanofluids. At the boundary, the skin friction factor assumes an increasing trend. Furthermore, the hybrid nanofluid exhibits a thicker momentum boundary layer at a fixed value of Eyring-Powell fluid parameter in comparison to the conventional nanofluid. Figure 2 illustrate how the dimensionless magnetic parameter affects velocity distribution. With an increase in the prescribed magnetic field $M$ , $f' (ξ, η)$ is seen to decay. The magnetic parameter $M$ gave rise to a drag force nominated the Lorentz force. The Lorentz force affects the fluid’s flow patterns by suppressing velocity in specific directions. It usually resists motion caused by velocity variations due to stretching surface, acting as a damping force. The influences of $K$ on the velocity portray is depicted in Figure 3, which also demonstrates that with higher levels of permeability, $f' (ξ, η)$ and the momentum boundary layer thickness decreases substantially. Permeability is the property of fluid motion resistance that decreases with increasing $K$ values, which in turn causes a decrease in velocity distribution. Improved heat transfer across the porous medium is made possible by stronger fiber packing, which offsets the hybrid nanoparticles in the fluid’s inability to function physically at higher K values.

Table 1.

Thermal properties of the Fe₃ O₄, Cu, and EG are discussed by Rashad et al.³⁵

Materials	$C_{p} (\frac{J}{Kkg})$	$ρ (\frac{Kg}{m^{3}})$	$k (\frac{W}{mK})$	$σ (\frac{S}{m})$
$F e_{3} O_{4}$	2090	5180	0.15	2500
$Cu$	385	8933	400	5.96 × 10⁷
EG	2430	1115	0.253	$1.07 \times 10^{- 11}$

Table 2.

Thermophysical properties for hybrid nanofluids and nanofluids are presented by Razzaq and Farooq, Razzaq et al.^36,37

Density $(Kg (m)^{- 3})$	$ρ_{hnf} = [ϕ_{2} ρ_{s 2} + (1 - ϕ_{2}) {ϕ_{1} ρ_{s 1} + (1 - ϕ_{1}) ρ_{f}}]$
Dynamic viscosity $(kg (ms)^{- 1})$	$μ_{hnf} = \frac{μ_{f}}{{(1 - ϕ_{2})}^{2.5} {(1 - ϕ_{1})}^{2.5}}$
Heat capacity $(J (kg . K)^{- 1})$	${(ρ C_{p})}_{hnf} = (1 - ϕ_{2}) {ϕ_{1} {(ρ C_{p})}_{s 1} + (1 - ϕ_{1}) {(ρ C_{p})}_{f}} + ϕ_{21} {(ρ C_{p})}_{s 2}$
Thermal conductivity $(W (mK)^{- 1})$	$\frac{k_{hnf}}{k_{bf}} L = \frac{k_{s 2} + (sf - 1) k_{bf} - (sf - 1) ϕ_{2} l (k_{bf} - k_{s 2})}{k_{s 1} + (sf - 1) k_{bf} + ϕ l (k_{bf} - k_{s 1})},$ $\frac{k_{bf}}{k_{f}} = \frac{k_{s 1} + (sf - 1) k_{f} - L (sf - 1) ϕ_{1} (k_{f} - k_{s 1})}{k_{s 1} + (sf - 1) k_{f} + L ϕ (k_{bf} - k_{s 1})},$
Electrical conductivity $(S m^{- 1})$	$\frac{σ_{hnf}}{σ_{bf}} = \frac{3 ϕ (σ_{1} ϕ_{1} + σ_{2} ϕ_{2 -} σ_{bf} (ϕ_{1} + ϕ_{2}))}{(σ_{1} ϕ_{1} + σ_{2} ϕ_{2} + 2 ϕ σ_{bf}) - ϕ σ_{bf} ((σ_{1} ϕ_{1} + σ_{2} ϕ_{2}) - σ_{bf} (ϕ_{1} + ϕ_{2}))} + 1,$ where $\frac{σ_{nf}}{σ_{f}} = 1 + \frac{3 (\frac{σ_{s}}{σ_{f}} - 1) ϕ}{(\frac{σ_{s}}{σ_{f}} + 2) - (\frac{σ_{s}}{σ_{f}} - 1) ϕ}$ ,

Table 3.

$- θ' (0)$ values are compared for various $\Pr$ inputs when the remaining parameters are zero.

$\Pr$	Ishak et al.³⁹	Alsaedi et al.⁴⁰	Present results
2	1.3349	1.3333	1.4236
5	3.2927	3.3165	3.8953
10	4.7742	4.7969	4.9523
15	5.9097	5.9320	5.9377

Table 4.

Results for skin friction for some values of $M, K, ε, ω and R$ .

$M$	$K$	$ε$	$ω$	$R$	$ϕ$	$- C_{fx} {(R e_{x})}^{\frac{1}{2}}$	$- Nu {(R e_{x})}^{\frac{1}{2}}$
0.2	0.5	0.2	0.2	0.5	0.04	1.8543	0.1736
0.5						2.1017	0.1766
1.0						2.3240	0.1798
0.5	0.1					1.9425	0.1766
	0.3					1.9838	0.1777
	0.5					2.0240	0.1792
	0.1	0.1				1.5029	0.17363
		0.2				1.4672	0.1688
		0.3				1.3341	0.1624
			0.1			1.6749	0.1762
			0.3			1.6812	0.1769
			0.5			1.6929	0.1798
				0.2		1.5016	0.1793
				0.5		1.5763	0.1889
				0.8		1.5765	0.1926
					0.02	0.1362	0.1453
					0.03	0.1455	0.1597
					0.05	0.1673	0.1761

Figure 2.

Variation of $M$ versus $f' (ξ, η)$ .

Figure 3.

Variation of $K$ versus $f' (ξ, η)$ .

As seen in Figures 4 and 5, the function with $R$ is in opposition to $f' (ξ, η), θ (ξ, η),$ and the thickness of the significant boundary layer. This is because increasing $R$ causes thermal diffusivity to lag, which in turn causes the area of the thermal boundary layer’s thickness to decrease. In Figures 6 and 7 the behavior of heat generation parameter $(Q)$ is represented with respect to velocity distribution $f' (ξ, η)$ and thermal distribution $θ (ξ, η)$ . The interpretation of Figures 8 and 9 clearly points to an increase in $Q$ values, $f' (ξ, η)$ and $θ (ξ, η)$ as a result of an improvement in thermal boundary layer thickness due to internal heat energy growth in the ascending heat source parameter circumstance. Figures 8 and 9 illustrate the velocity diagram for varying volume fractions of $F e_{3} O_{4}$ and $Cu,$ when $0 \leq ξ \leq 1$ . The hybrid nanofluid density and dynamic viscosity increase as the volume fraction rises. As shown in these figures, this subsequently modifies $f' (ξ, η)$ and results in decreased momentum generation. As can be shown in Figure 10, the temperature directly depends on the $Cu$ nanoparticles. This pattern is consistent with the fact that copper nanoparticles have higher thermal conductivity than basic fluids. Figure 11 presents the temperature account influenced by the momentum to mass diffusivity ratio; a decreasing trend is noted toward increased values of the parameter that is employed, namely $\Pr$ . The main reason for the temperature drop is that when $\Pr$ increases, thermal diffusivity falls. Because fluids with higher Prandtl numbers have relatively poorer thermal conductivity, conduction, and the thickness of the thermal boundary layer decreases, lowering temperature. Figures 12 and 13 display the computed estimates of skin friction and Nusselt number for various values of $M$ versus nonsimilarity variable $(ξ)$ . The figures indicate that an increase in $M$ leads to an improvement in both heat transmission and skin friction coefficient because of the Lorentz force’s resistance; nevertheless, an opposite pattern is observed as $ξ$ increases. With respect to the mixed convection parameter, and for various permeability values, Figures 14 and 15 illustrate that the permeability parameter exhibits the same behavior as the magnetic factor on the Nusselt number and friction coefficient. The Nusselt number and skin friction profiles against $ξ$ with different inputs of $R$ are shown in Figures 16 and 17. When the $ξ$ and R parameters rise, the drag force and Nusselt factor also rise. The reason for this is that when the flux improves and absorbs more heat, the momentum boundary layer thickness increases.

Figure 4.

Variation of $R$ versus $f' (ξ, η)$ .

Figure 5.

Variation of $R$ versus $θ (ξ, η)$ .

Figure 6.

Variation of $Q$ versus $f' (ξ, η)$ .

Figure 7.

Variation of $Q$ versus $θ (ξ, η)$ .

Figure 8.

Variation of $ϕ_{1}$ versus $θ (ξ, η)$ .

Figure 9.

Variation of $ϕ_{2}$ versus $f' (ξ, η)$ .

Figure 10.

Variation of $ϕ_{2}$ versus $θ (ξ, η)$ .

Figure 11.

Variation of $\Pr$ versus $θ (ξ, η)$ .

Figure 12.

Variation of $M$ versus $R e^{1 / 2} Cf$ .

Figure 13.

Variation of $M$ versus $R e^{- 1 / 2} Nu$ .

Figure 14.

Variation of $K$ versus $R e^{1 / 2} Cf$ .

Figure 15.

Variation of $K$ versus $R e^{- 1 / 2} Nu$ .

Figure 16.

Variation of $R$ versus $R e^{1 / 2} Cf$ .

Figure 17.

Variation of $R$ versus $R e^{- 1 / 2} Nu$ .

Conclusion

This article examines how heat generation and radiative heat affect the mixed convection flux in a porous medium containing a magnetic Eyring–Powell hybrid nanofluid. With the aid of non similarity method which solved by using the built-in Matlab algorithm bvp4c, the results have been obtained. The main conclusions are.

• With an increase in the parameter values of $M, K, R, ϕ_{1}$ , and $ϕ_{2}$ , $f' (ξ, η)$ is dropped. However, it is the heat sink/source that exhibits the contrary behavior.

• Raising the values of $Q$ and $ϕ_{2}$ parameters cause an increase in $θ (ξ, η) .$

• Thermal enhancement is caused by an increase in thermal radiation and $F e_{3} O_{4}$ nanoparticle concentrations.

• Drag force decreases with increasing magnetic field, porous material permeability, and heat radiation.

• As the magnetic field, the permeability of porous media, and the heat source/sink increase, heat transfer decreases.

• It is estimated that the thermal conductivity and viscosity of the ${Fe}_{3} O_{4} + Cu / EG$ hybrid nanofluid are expected to increase to 4.1% and 12.35%, respectively, and to 71.55% and 78.01%, respectively, for volume concentrations ranging from 0.02% to 0.05 and temperatures between 20°C and 60°C.

The methodologies and findings of this study have practical implications for various engineering applications, including aerospace, automotive, and environmental engineering. Researchers in these fields can leverage this investigation insight to improve the design and performance of systems involving boundary layer flows. By contributing these advancements, the study serves as a valuable resource for researchers looking to solve complex boundary layer problems more effectively and efficiently.

Footnotes

Appendix

Handling Editor: Sharmili Pandian

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 iDs

Saiqa Sagheer

Umer Farooq

References

Shafiq

Sindhu

TN.

Statistical study of hydromagnetic boundary layer flow of Williamson fluid regarding a radiative surface. Results Phys 2018; 7: 3059–3067.

Shekaramiz

Fathi

Ataabadi

, et al. MHD nanofluid free convection inside the wavy triangular cavity considering periodic temperature boundary condition and velocity slip mechanisms. Int J Therm Sci 2021; 170: 107179.

Razzaq

Farooq

Mirza

HR.

Nonsimilar forced convection analysis of Maxwell nanofluid flow over an exponentially stretching sheet with convective boundary conditions. Z Angew Math Mech 2023; 103: e202200623.

Shah

Hayat

Alsaedi

Entropy optimization in a fourth-grade nanofluid flow over a stretchable Riga wall with thermal radiation and viscous dissipation. Int Commun Heat Mass Transf 2021; 127: 105398.

Khan

Abrar

Uddin

, et al. Entropy generation due to micro-rotating Casson’s nanofluid flow over a nonlinear stretching plate: numerical treatment. Waves Random Complex Media. Epub ahead of print 19 May 2022. DOI: 10.1080/17455030.2022.2067376.

Farooq

Liu

Farooq

, et al. Non-similar analysis of bioconvection MHD micropolar nanofluid on a stretching sheet with the influences of Soret and Dufour effects. Appl Water Sci 2024; 14: 116.

Ali

Thumma

Habib

, et al. Finite element analysis on transient MHD 3D rotating flow of Maxwell and tangent hyperbolic nanofluid past a bidirectional stretching sheet with Cattaneo Christov heat flux model. Therm Sci Eng Progress 2022; 28: 101089.

Khattak

Ahmed

Abrar

, et al. Numerical simulation of Cattaneo–Christov heat flux model in a porous media past a stretching sheet. Waves in Random and Complex Media. Epub ahead of print 9 February 2022. DOI: 10.1080/17455030.2022.2030503.

Razzaq

Farooq

Aldandani

Nonsimilar convection analysis of single and multilayer carbon nanotubes based nanofluid flow over a vertical cone in a complex porous media subjected to thermal radiations and chemical reaction. J Magn Magn Mater 2023; 572: 170583.

10.

Chamkha

Armaghani

Mansour

, et al MHD convection of an Al2O3Cu/water hybrid nanofluid in an inclined porous cavity with internal heat generation/absorption. Iran J Chem Chem Eng 2021; 41: 936–956.

11.

Cao

Animasaun

Yook

, et al. Simulation of the dynamics of a colloidal mixture of water with various nanoparticles at different levels of partial slip: ternary-hybrid nanofluid. Int Commun Heat Mass Transf 2022; 135: 106069.

12.

Algehyne

El-Zahar

Sohail

, et al. Thermal improvement in pseudo-plastic material using ternary hybrid nanoparticles via non-Fourier’s law over porous heated surface. Energies 2021; 14: 8115.

13.

Sarada

Gamaoun

Abdulrahman

, et al. Impact of the exponential form of internal heat generation on water-based ternary hybrid nanofluid flow by capitalizing non-Fourier heat flux model. Case Stud Therm Eng 2022; 38: 102332.

14.

Elnaqeeb

Animasaun

Shah

NA.

Ternary-hybrid nanofluids: significance of suction and dual-stretching on three-dimensional flow of water conveying nanoparticles with various shapes and densities. Z Nat Forsch A 2021; 6: 231–243.

15.

Khazayinejad

Nourazar

SS.

On the effect of spatial fractional heat conduction in MHD boundary layer flow using Gr-Fe₃O₄–H₂O hybrid nanofluid. Int J Therm Sci 2022; 172: 107265.

16.

Roşca

Pop

Flow and heat transfer of Powell–Eyring fluid over a shrinking surface in a parallel free stream. Int J Heat Mass Transf 2014; 71: 321–327.

17.

El-Zahar

Rashad

Saad

, et al. Magneto-hybrid nanofluids flow via mixed convection past a radiative circular cylinder. Sci Rep 2020; 10: 10494.

18.

Cui

Razzaq

Azam

, et al. Nonsimilar forced convection analysis of chemically reactive magnetized Eyring-Powell nanofluid flow in a porous medium over a stretched riga surface. J Porous Media 2022; 25: 67–81.

19.

Mushtaq

Mustafa

Hayat

, et al. Exponentially stretching sheet in a Powell-Eyring fluid: numerical and series solutions. Z Nat Forsch A 2013; 68: 791–798.

20.

Rahimi

Ganji

Khaki

, et al. Solution of the boundary layer flow of an Eyring-Powell non-Newtonian fluid over a linear stretching sheet by collocation method. Alex Eng J 2017; 56: 621–627.

21.

Akbar

Ebaid

Khan

ZH.

Numerical analysis of magnetic field effects on Eyring-Powell fluid flow towards a stretching sheet. J Magn Magn Mater 2015; 382: 355–358.

22.

Hassan

Ahsan Usman

, et al. Entropy generation and flow characteristics of Powell Eyring fluid under effects of time sale and viscosities parameters. Sci Rep 2023; 13: 8376.

23.

Ghadikolaei

Hosseinzadeh

Ganji

DD.

Analysis of unsteady MHD Eyring-Powell squeezing flow in stretching channel with considering thermal radiation and Joule heating effect using AGM. Case Stud Therm Eng 2017; 10: 579–594.

24.

Jan

Mushtaq

Farooq

, et al. Nonsimilar analysis of magnetized Sisko nanofluid flow subjected to heat generation/absorption and viscous dissipation. J Magn Magn Mater 2022; 564: 170153.

25.

Cui

Jan

Farooq

, et al. Thermal analysis of radiative Darcy–Forchheimer nanofluid flow across an inclined stretching surface. Nanomaterials 2022; 12: 4291.

26.

Abrar

Ul Haq

Awais

, et al. Entropy analysis in a cilia transport of nanofluid under the influence of magnetic field. Nucl Eng Technol 2017; 49: 1680–1688.

27.

Kumar

Reddy

Rao

, et al. Thermal radiation impact on MHD heat transfer natural convective nanofluid flow over an impulsively started vertical plate. Case Stud Therm Eng 2021; 24: 100826

28.

Agbaje

Mondal

Motsa

, et al. A numerical study of unsteady non-Newtonian Powell-Eyring nanofluid flow over a shrinking sheet with heat generation and thermal radiation. Alex Eng J 2017; 56: 81–91.

29.

Oreyeni

Akindele

Obalalu

, et al. Thermal performance of radiative magnetohydrodynamic Oldroyd-B hybrid nanofluid with Cattaneo-Christov heat flux model: solar-powered ship application. Numeri Heat Transf A Appl 2024; 85: 1954–1972.

30.

Sparrow

HS.

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

31.

Farooq

Hussain

Ijaz

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

32.

Sagheer

Farooq

Hussain

Non-similar investigation of enhanced thermal efficiency of Maxwell based hybrid nanofluid (MoS₂+ZnO) across a linearly stretched surface. J Magn Magn Mater 2023; 565: 170285.

33.

Riaz

Naheed

Farooq

, et al. Non-similar investigation of magnetized boundary layer flow of nanofluid with the effects of Joule heating, viscous dissipation and heat source/sink. J Magn Magn Mater 2023; 574: 170707.

34.

Sagheer

Razzaq

Vafai

Local non-similar solutions for EMHD nanofluid flow with radiation and variable heat flux along a stretching sheet. Numeri Heat Transf A Appl. Epub ahead of print 3 June 2024. DOI: 10.1080/10407782.2024.2360082.

35.

Rashad

Nafe

Eisa

DA.

Heat generation and thermal radiation impacts on flow of magnetic Eyring–Powell hybrid nanofluid in a porous medium. Arab J Sci Eng 2023; 48: 939–952.

36.

Razzaq

Farooq

Non-similar analysis of MHD hybrid nanofluid flow over an exponentially stretching/shrinking sheet with the influences of thermal radiation and viscous dissipation. Numeri Heat Transf B Fundamentals. Epub ahead of print 13 February 2024. DOI: 10.1080/10407790.2024.2312958.

37.

Razzaq

Sagheer

Farooq

Local non-similar solutions of magnetohydrodynamic casson nanofluid flow over a non-linear inclined surface with thermal radiation and heat generation effects: a utilization of up to third truncation. J Porous Media 2024; 27: 67–83.

38.

Minkowycz

Sparrow

EM.

Local nonsimilar solutions for natural convection on a vertical cylinder. J Heat Transfer 1974; 1: 178–183.

39.

Ishak

Nazar

Pop

Mixed convection on the stagnation point flow toward a vertical, continuously stretching sheet. J Heat Transfer 2007; 129: 1087–1090.

40.

Alsaedi

Hayat

Qayyum

, et al. Eyring–Powell nanofluid flow with nonlinear mixed convection: entropy generation minimization. Comput Methods Programs Biomed 2020; 186: 105183.