Numerical simulation of hybrid nanofluid flow with homogeneous and heterogeneous chemical reaction across an inclined permeable cylinder/plate

Abstract

The current study investigates the steady two-dimensional (2D) hybrid nanofluid (Hnf) flow over an inclined permeable plate/cylinder. The Hnf flow has been examined in the context of mixed convection, heterogeneous/homogenous chemical reaction, and permeable medium. The Hnf is prepared by dispersing silver (Ag), and iron ferrite (Fe₃O₄) nanoparticles (NPs) in water. The current research is motivated by the increasing demand for highly efficient cooling devices in a variety of industries and energy-associated operations. The energy transmission and fluid flow are mathematically specified by using a coupled nonlinear system of partial differential equations (PDEs). The system of PDEs is simplified into a dimensionless form of ODEs, which are then further numerically treated with the MATLAB package based on the finite difference method (bvp4c). It has been noticed that the permeability component develops the heat transfer curve while decreasing the flow rate of the fluid. The impact of heat source/sink increases the energy profile. Moreover, the plate surface demonstrates the dominant behavior of energy transportation than a cylinder with the variance of Ag-Fe₃O₄-NPs.

Keywords

Exponential heat source/sink hybrid nanofluid numerical solution (bvp4c)cylinder/plate homogeneous/heterogeneous chemical reaction

Introduction

Fluid flow over a plate or cylinder has a wide range of industrial applications. Some practical uses in industrial sectors involve heat exchangers, aerodynamics, cooling methods, and ecological engineering (Rafique, Mahmood, Adnan et al., 2024). Waqas et al. (2021) addressed the Darcy-Forchheimer bioconvective nanoliquid flow using a stretching plate/cylinder with modified mass and heat flux characteristics. It was revealed that the fluid velocity drops with the variation of the local inertia coefficient and mixed convection parameter. Xu et al. (2021) used a numerical method to assess the steady power-law nanofluids flow consisting of gyrotactic microbes through plates with the energy exchange and determined that the fluid velocity demonstrates an upsurge as the buoyancy convection factor increases. Bilal et al. (2021) calculated the mixed convection Darcy flow of a hybrid nanoliquid containing CNTs and Fe₃O₄-NPs across a prone expanding cylinder. The outcomes indicated that the hybrid nanoliquid is the effective source for enhancing heat conveyance and may be utilized as a cooling agent for industrial appliances. Algehyne et al. (2022) deliberated the hydromagnetic nanoliquid flow containing roaming microbes and NPs through a penetrable plate. The porosity parameter, buoyancy ratio, and inertial parameter have an adverse impact on the distribution of velocity. Alharbi et al. (2022) reported on the ternary nanoliquid flow with heat transmission involving metallic NPs across a prolonged cylindrical with the influence of magnetic induction using the bvp4c package. Salawu et al. (2023) used hybridized SWCNT-Ag and MWCNT-MoS4 magneto-nanomaterials to evaluate heat radiation and convection in engine oil throughout a cylinder. Bin Mizan et al. (2021) studied magnetic induction effects on an electroconductive incompressible nanofluid flowing over a cylindrical cylinder. Heat was distributed in a boundary layer using metallic nanoparticles or ferroparticles during flow. Raizah et al. (2023) investigated the 3D flow of Ag, MgO, and motile microbes’ water-based hybrid nanoliquids as they flowed across a circular cylinder sinusoidal in shape. It was noticed that the quantity of hybrid nanoparticles (Ag-MgO) enhances the energy fields. Shamshuddin et al. (2024) created a novel radiative nanoliquid flow mathematical model. Stretching a nanoparticle-immersed sheet with an inclined magnetic field creates this model. Some remarkable results related to fluid flow across cylinder/plate are presented by Ali, Jubair, Mahmood et al., 2024; Rafique, Mahmood, Ali et al., 2024; and Vijatha and Bala Anki Reddy, 2024.

The study of fluid flow over a cylinder or plate has several industrial and practical applications, such as aerodynamics, wind energy, heat exchangers, civil engineering, marine engineering, biomedical applications, and renewable energy (Islam, Mahmood, Khan, Ali et al., 2024a; Shah and Ali, 2019). In the current analysis, Ag and Fe₃O₄ NPs are considered in the based fluid water. Abbasi, Gul and Shehzad (2020) examined the effectiveness of Au, Cu, and Fe₃O₄ NPs on peristalsis flow hybrid nanofluids. Muhammad et al. (2021) highlighted the viscous dissipation and meting effect on the hybrid nanofluid consisting of CNTs and CuO-NPs. Larger values of the squeezing factor, melting constraint, and CNTs-CuO-NPs upsurge the fluid's velocity. Alqahtani et al. (2023) and Islam, Mahmood, Khan et al. (2024b) studied an incompressible magnetohydrodynamics nanofluid consisting of CoFe₂O₃-NPs flowing over a vertical plate. It was observed that the presence of CoFe₂O₃-NPs increased the mass conversation rate to 17.102% whereas increasing the energy transmission rate to 13.92%. Ahmad et al. (2024) studied the hybrid nanoliquid (Cu and Al₂O₃) flow surrounding the spinning sphere. It was observed that the hybrid nanoliquid (Cu and Al₂O₃/water) velocity is amplified by the influence of rotational and acceleration parameters. Adnan, Abbas, Said et al. (2024) scrutinized the thermal endurance of a water-based nanofluid containing oxide and metallic NPs (SiO₂, Au, and MoS₂). It was examined that the ranges for MoS₂ NPs concentrations (0.1–0.6%) and Eckert number efficiently improved the fluidity of the working fluids. Adnan, Abbas, Alqahtani et al. (2024) introduced a novel ternary nanoliquid model that describes the thermal radiations, mixed convection, heating sources, and generated magnetic fields. Yasmin et al. (2024) set out to investigate the effects of reactive species, heat sources and sinks, and hybrid nanofluids on mass and energy transmissions over a stretched cylinder. Further related results are presented by Ali et al. (2023) and Acharya (2024).

In many different applications, including chemical reactors, gasoline and diesel engines, prevention of pollution, catalytic processes, atmospheric chemistry, wastewater treatment, electrochemistry, energy production, and green technologies, homogeneous and heterogeneous (HH) chemical reactions have a substantial impact on fluid flow (Damseh et al., 2009; Ghailan et al., 2023; Maraj et al., 2023). Haq et al. (2022) studied the upshot of a porous media, and HH chemical reactions on the hybrid nanoliquid flow over three different surfaces (conical shape, wedge, and plate). The decrease in HH reaction enhances the concentration profile. Pattnaik, Bhukta and Mishra (2022) inspected the effects of a nonuniform HH reaction on the Micropolar nano liquid flow via a vertically extending surface. Naveed, Imran and Gul (2023) and Bani-Fwaz et al. (2024) reported the heat exchange through nanoliquid in a channel created by extending and constricting walls, with the impact of thermal radiation, surface flexibility, and HH reaction. Waqas et al. (2024) calculated the nanoliquid flow through a vertical prolonging surface using numerical methods. It is evident that activating energy increases mass disseminating rate whereas chemical reactions lead it to decrease with the impact of magnetic field, HH reaction, and thermal radiation. Takhar, Chamkha and Nath (2000) studied mixed convection flow over a moving vertical thin cylinder, considering heat and mass diffusion. Some valuable results are recently delivered by Ali, Jubair and Siddiqui (2024); Ali, Jubair, Mahmood et al. (2024); and Xie et al. (2024).

This work is unique in using porous medium. Porous media enable fluids to pass through and interact with the solid matrix. Porous medium is essential in most academic fields. Computational fluid dynamics that account for spongy bone tissue's porosity may help predict drug distribution in medical applications. Limestone, beach sand and sandstone naturally store oil, gas and water. Some paper printing components, battery fuel cells, and baked product filters need porous media. Darcy pioneered porous media research to improve water distribution systems, according to Ali J Chamkha (1997). Darcy found that fluid flow through a porous medium is directly related to pressure difference and cross-sectional area and inversely proportional to viscosity. In porous media, conduction through the solid matrix and convection within the fluid inhabiting the pores interact to change heat transmission. Because of this complexity, managing heat transport within porous materials is challenging. To improve heat transfer processes in scientific, technological, and environmental sectors, these impacts must be understood (Mahmood, Rafique, Khan, 2024b). Because Darcy's rule is simple and robust, it may be used directly in numerical simulations for many purposes. Thus, the Darcy flow model in the boundary layer has been used in various recent porous medium studies. The effects of porous media on fluid dynamics and heat transport in Newtonian fluids are studied (Chamkha and Rashad, 2012; Chamkha, 1997).

In the previous paragraphs, the researcher addressed fluid concerns using typical methods. Nowadays, researchers use recent computer developments, such as optimization methods in machine learning models, have changed fluid mechanics. Artificial neural networks (ANNs) with the Levenberg-Marquardt (LM) strategy can solve fluid flow problems. The LM scheme optimizes nonlinear least squares problems and fits mathematical models that describe fluid flow behavior to experimental or simulated data. The LM technique minimizes sum of squared errors between observed data and model predictions to enhance fluid flow model parameter estimation. This helps us forecast flow patterns, turbulence, and other fluid phenomena. ANNs, which are affected by the human brain, can evaluate complex fluid flow patterns, forecast, and simulate turbulence and fluctuating flows. These approaches give a complete foundation for fluid mechanics’ complexity. Akbar, Zamir and Muhammad (2024) created an AI-driven neural network utilizing the LM approach to examine liquid mass and heat transport across a cone-disk device's conical gap. Akbar et al. (2024) used the LM technique, known for its innovative approach and convergent stability, to numerically model the evolution of a hybrid Cu-Al₂O₃/water nanofluid on a porous stretched sheet (HNF-PSS) using regression plots, state transition measures, histogram representations, and mean squared errors. Alghamdi et al. (2024) study a dual stratified common on the diversified convection barrier layer discharge of an Eyring-Powell fluid induced by a prone extended barrel. Concentration and temperature outside the barrel may be greater than in the moving stream. The flow equations are solved using a well-established numerically-based LM neural network scheme (LMNNS) aggregating solver.

In the current study, Hnf flow over an inclined permeable plate/cylinder. The Hnf flow has been observed under the effect of heterogeneous/homogeneous reaction, and permeable medium. The Hnf is prepared by dispersing Ag, and Fe₃O₄-NPs in water. The energy transmission and fluid flow are mathematically specified by using a coupled nonlinear system of PDEs. The system of PDEs is simplified into a dimensionless form of ODEs, which are then further numerically treated with the bvp4c. The main objectives of the present study are:

The current research is motivated by the increasing demand for highly efficient cooling devices in a variety of industries and energy-associated operations.

To investigate the steady 2D hybrid nanoliquid flow over an inclined absorbent plate/cylinder.

The Hnf flow is observed under the impact of mixed convection, heterogeneous/homogeneous chemical reaction, and permeable medium.

To obtain the numerical solution of the nonlinear system of PDEs using Matlab built-in package bvp4c.

To examine the influence of Ag-Fe₃O₄-NP on the velocity and energy transmission rate.

Applications of Ag-Fe₃O₄-NP: Nanoparticles of silver (Ag) and iron oxide (Fe₃O₄) improve heat transmission by increasing the base fluid's thermal conductivity. The nanofluid can be controlled with magnetic fields because Fe₃O₄ makes it magnetic. Applications requiring precision cooling targeting benefit from this. The hybrid nanofluid exhibits better stability and dispersion than single-component nanofluids, reducing particle agglomeration and sedimentation.

Mathematical formulation

The 2D laminar Hnf flow over a permeable inclined plate/cylinder is studied subject to the permeable medium, mixed convection, and heat source/sink. The silver (Ag) and iron oxide (Fe₃O₄) nanoparticles are used in water for the preparation of Hnf. x and r reveals the axial and radial directions, where u indicates the reference velocity as displayed in Figure 1. The ambient and surface temperature of the plate and cylinder is denoted by $T_{\infty}$ and $T_{w} .$ besides, external forces as well as pressure gradients are considered to have no impact on the nanofluid flow.

Figure 1.

Fluid flow over an inclined cylinder/plate.

Furthermore, the nanofluid model including the heterogeneous and homogeneous reactions across two chemical components, A* and B* has been assumed. A cubic homogeneous (autocatalytic) reaction is considered in this model within the boundary layer flow. Similarly, a first-order heterogeneous reaction takes place on the catalyst surface, stated as Pooja, Narasimhamurthy and Anitha (2024):

\begin{aligned} A * + 2 B * \to 3 B *, rate = l_{1} a b^{2}, \\ A * \to B *, rate = l_{s} a . \end{aligned}

Here, a and b characterize the chemical species concentrations A* and B*, whereas

l_{1} and l_{s}

signify the constant rate.

The modeled equations based on the above suppositions are framed as Madhukesh et al. (2023) and Pooja, Narasimhamurthy and Anitha (2024):

\frac{\partial (r u)}{\partial x} + \frac{\partial (r v)}{\partial r} = 0,

(1)

u \frac{\partial (u)}{\partial x} + v \frac{\partial (u)}{\partial r} = v_{h n f} (\frac{\partial^{2} u}{\partial r^{2}} + \frac{1}{r} \frac{\partial u}{\partial r}) + \frac{{(ρ β)}_{h n f} g (T_{1} - T_{\infty}) \cos ς}{ρ_{h n f}} - \frac{v_{h n f}}{K_{1} *} u,

(2)

u \frac{\partial (T)}{\partial x} + v \frac{\partial (T)}{\partial r} = α_{h n f} (\frac{\partial^{2} T}{\partial r^{2}} + \frac{1}{r} \frac{\partial T}{\partial r}) - \frac{1}{{(ρ C_{p})}_{h n f}} \frac{\partial}{\partial y} (q_{r}) + \frac{Q_{e} *}{{(ρ C_{p})}_{h n f}} (T - T_{\infty}) e x p (- \sqrt{\frac{a}{υ_{f}} n y}),

(3)

u \frac{\partial a}{\partial x} + v \frac{\partial a}{\partial r} = D_{A} (\frac{\partial^{2} a}{\partial r^{2}} + \frac{1}{r} \frac{\partial a}{\partial r}) - l_{1} a b^{2},

(4)

u \frac{\partial b}{\partial x} + v \frac{\partial b}{\partial r} = D_{B} (\frac{\partial^{2} b}{\partial r^{2}} + \frac{1}{r} \frac{\partial b}{\partial r}) - l_{1} a b^{2},

(5)

In equation (3),

q_{r}

is defined as (Mahmood, Rafique, Khan, et al., 2024a):

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

(6)

The energy equation becomes:

\begin{aligned} u \frac{\partial (T)}{\partial x} + v \frac{\partial (T)}{\partial r} = α_{h n f} (\frac{\partial^{2} T}{\partial r^{2}} + \frac{1}{r} \frac{\partial T}{\partial r}) - \frac{1}{{(ρ C_{p})}_{h n f}} \frac{16 σ *}{3 k *} T_{h}^{3} \frac{\partial^{2} T}{\partial y^{2}} \\ + \frac{Q_{e} *}{{(ρ C_{p})}_{h n f}} (T - T_{\infty}) e x p (- \sqrt{\frac{a}{υ_{f}} n y}) . \end{aligned}

(7)

The boundary conditions (BCs) are (Madhukesh et al., 2023):

\begin{matrix} u = u_{w} = \frac{U_{0} * x}{l}, v = 0, - k \frac{\partial r}{\partial T} = h_{f} (T_{w} - T), D_{A} \frac{\partial a}{\partial r} = l_{s} a, D_{B} \frac{\partial b}{\partial r} = l_{s} a \\ at r = R u \to 0, T \to T_{\infty}, a \to a_{0}, b \to 0, r \to \infty \end{matrix}}

(8)

In equation (7) consider exponentially space-based heat source/sink term: The model needs an increasingly space-based heat source/sink to better reflect thermal systems in real life, where heat output and absorption vary by location. This method captures the uneven heat distribution, which is essential for understanding and improving thermal efficiency and disorder in nanofluid flows in real-world applications like electronics cooling systems and energy systems that require accurate thermal control.

The similarity variables are (Madhukesh et al., 2023):

\begin{aligned} ψ = \sqrt{u_{w} v_{f} x} R f (η), θ = \frac{T - T_{\infty}}{T_{w} - T_{\infty}}, u = \frac{1}{r} \frac{\partial ψ}{\partial r}, v = - \frac{1}{r} \frac{\partial ψ}{\partial x}, h (η) = \frac{a}{a_{0}}, g (η) = \frac{b}{a_{0}}, \\ η = \sqrt{\frac{u_{w}}{v_{f} x}} (\frac{r^{2} - r^{2}}{2 R}) . \end{aligned}

(9)

By using equation (9) in equations (1)–(5) and (8), the modeled equations become:

\frac{((1 + (2 δ_{1}) η) f^{″'} + (2 δ_{1}) f^{″})}{B_{1} B_{2}} - (f^{'})^{2} - \frac{P_{m}}{B_{1} B_{2}} f^{'} + f f^{″} + \frac{B_{3}}{B_{2}} γ_{1} θ \cos ζ = 0,

(10)

(1 + R d \frac{k_{T n f}}{k_{f} P r B_{4}}) ((1 + 2 δ_{1} η) θ^{″} + 2 δ_{1} θ^{'}) + f θ^{'} + Q_{e} \exp (- n η) .

(11)

(1 + 2 η δ_{1}) h^{″} + δ_{1} h^{'} + S c f h^{'} - S c L_{1} h g^{2},

(12)

δ * (1 + 2 η δ_{1}) g^{″} + δ * δ_{1} g^{'} + S c f g^{'} - S c L_{1} h g^{2},

(13)

The BCs for ODEs are:

\begin{matrix} f (η) = 0, f^{'} (η) = 1, θ^{'} (η) = - λ_{1} (1 + θ (η)), h^{'} (η) = L_{s} h (0), δ * g^{'} (0) = - L_{s} h (0) \\ at η = 0 f^{'} (\infty) = 0, θ (\infty) = 0, h (\infty) \to 1, g (\infty) \to 0 a s η \to \infty \end{matrix}}

(14)

The nondimensional constant and parameters are expressed in Table 1.

Table 1.

Physical parameter along with the symbols and expressions (Madhukesh et al., 2023; Pooja, Narasimhamurthy and Anitha, 2024).

Thermal Biot number	$λ_{1}$	$λ_{1} = h_{f} \sqrt{\frac{υ l}{U_{0} *}}$
Curvature factor	$δ_{1}$	$δ_{1} = \sqrt{\frac{v_{f} l}{U_{0} * r^{2}}}$
Cylinder surface	$δ_{1} > 0$
Plate surface	$δ_{1} = 0$
Mixed convection factor	$γ_{1}$	$γ_{1} = \frac{G r}{R e^{2}} = \frac{g β (T_{w} - T_{\infty}) l}{u_{w} U_{0} *}$
Porosity factor	$P_{m}$	$P_{m} = \frac{v_{f} l}{U_{0} * K *}$
Prandtl number	$P r$	$P r = \frac{v_{f}}{α_{f}}$
Heat source/sink constraint	$Q_{e}$	$Q_{e} = \frac{Q_{e} * l}{a {(ρ C_{p})}_{f}}$
Grashof number	$G r$	$G r = \frac{g β (T_{w} - T_{\infty}) l^{3}}{ν_{f}^{2}}$
Thermal radiation constant	Rd	$R d = \frac{16 σ * T_{h}^{3}}{3 k * k_{f}}$
Mass diffusion coefficient	$δ *$	$δ * = \frac{D_{A}}{D_{B}}$
Homogeneous reaction parameter	$L_{1}$	$L_{1} = \frac{l_{1} a_{0}^{2} l}{U}$
Heterogeneous reaction parameter	$L_{s}$	$L_{s} = \frac{l_{s}}{D_{A}} \sqrt{\frac{v l}{U}}$
Schmidt number	Sc	$S c = \frac{ν}{D_{A}}$

where

B_{1} = \frac{μ_{h n f}}{μ_{b f}}, B_{2} = \frac{ρ_{h n f}}{ρ_{b f}}, B_{3} = \frac{{(ρ β)}_{h n f}}{{(ρ β)}_{b f}}, B_{4} = \frac{{(ρ C_{p})}_{h n f}}{{(ρ C_{p})}_{b f}}, B_{5} = \frac{k_{h n f}}{k_{b f}}, B_{6} = \frac{σ_{h n f}}{σ_{b f}} .

Tables 2 and 3 show the numerical values and thermophysical properties of nanoparticles and base fluid.

Table 2.

The numerical values of water, Ag, and Fe₃O₄ (Zhang et al., 2021).

	$ρ (kg / m^{3})$	$C_{p} (j / kgK)$	$k (W / mK)$
Pure water H₂O	997.1	4179	0.613
Ag	10,500	235	429
Fe₃O₄	5200	670	6

Table 3.

Thermophysical properties of Hnf (Zhang et al., 2021).

Properties
Viscosity	$\frac{μ_{h n f}}{μ_{b f}} = \frac{1}{{(1 - ϕ_{F e_{3} O_{4}} - ϕ_{A g})}^{2}}$
Density	$\frac{ρ_{h n f}}{ρ_{b f}} = ϕ_{F e_{3} O_{4}} (\frac{ρ_{F e_{3} O_{4}}}{ρ_{b f}}) + ϕ_{A g} (\frac{ρ_{A g}}{ρ_{b f}}) + (1 - ϕ_{F e_{3} O_{4}} - ϕ_{A g})$
Thermal capacity	$\frac{{(ρ C_{p})}_{h n f}}{{(ρ C_{p})}_{b f}} = ϕ_{F e_{3} O_{4}} (\frac{{(ρ C_{p})}_{F e_{3} O_{4}}}{{(ρ C_{p})}_{b f}}) + ϕ_{A g} (\frac{{(ρ C_{p})}_{A g}}{{(ρ C_{p})}_{b f}}) + (1 - ϕ_{F e_{3} O_{4}} - ϕ_{A g})$
Thermal conductivity	$\frac{k_{h n f}}{k_{b f}} = [\frac{(\frac{ϕ_{F e_{3} O_{4}} k_{F e_{3} O_{4}} + ϕ_{A g} k_{A g}}{ϕ_{F e_{3} O_{4}} + ϕ_{F e_{2} O_{4}}}) + 2 k_{b f} + 2 (ϕ_{F e_{3} O_{4}} k_{F e_{3} O_{4}} + ϕ_{A g} k_{A g}) - 2 (ϕ_{F e_{3} O_{4}} + ϕ_{A g}) k_{b f}}{(\frac{ϕ_{F e_{3} O_{4}} k_{F e_{3} O_{4}} + ϕ_{A g} k_{A g}}{ϕ_{F e_{3} O_{4}} + ϕ_{A g}}) + 2 k_{b f} - 2 (k_{F e_{3} O_{4}} ϕ_{F e_{3} O_{4}} + k_{A g} ϕ_{A g}) + (ϕ_{F e_{3} O_{4}} + ϕ_{A g}) 2 k_{b f}}]$
Electrical conductivity	$\frac{σ_{h n f}}{σ_{b f}} = [\frac{(\frac{ϕ_{F e_{3} O_{4}} σ_{F e_{3} O_{4}} + σ_{A g} ϕ_{A g}}{ϕ_{F e_{2} O_{4}} + ϕ_{F e_{3} O_{4}}}) + 2 σ_{b f} + 2 (ϕ_{F e_{3} O_{4}} σ_{F e_{3} O_{4}} + ϕ_{A g} σ_{A g}) - 2 (ϕ_{F e_{3} O_{4}} + ϕ_{A g}) σ_{b f}}{(\frac{ϕ_{F e_{3} O_{4}} σ_{F e_{3} O_{4}} + ϕ_{A g} σ_{A g}}{ϕ_{F e_{3} O_{4}} + ϕ_{A g}}) + 2 σ_{b f} - (ϕ_{F e_{3} O_{4}} σ_{F e_{3} O_{4}} + ϕ_{A g} σ_{A g}) + (ϕ_{F e_{3} O_{4}} + ϕ_{A g}) σ_{b f}}]$

Hnf: hybrid nanofluid.

For skin friction and Nusselt number along with their transform form are:

\begin{aligned} C_{f} = {\frac{\partial u}{\partial r} |}_{r = R} (\frac{μ_{h n f}}{ρ_{f} u w^{2}}), N_{u} = {\frac{\partial T}{\partial r} |}_{r = R} (\frac{- x k_{h n f}}{k_{f} (T_{w} - T_{\infty})}), \\ C_{f} = \frac{f^{″} (0)}{B_{1} {(R e)}^{0.5}}, N_{u} = \frac{- (k_{h n f} + R d) {(R e)}^{0.5}}{k_{f}} θ^{'} (0) . \end{aligned}

(15)

In the coming section, equations (10)–(13) and (14) are numerically solved.

Numerical solution and validation

The MATLAB boundary value problem solver BVP4c is built around the finite difference approach. The concept of finite difference is in fact utilized in the framework of BVP4c to estimate differential equation solutions despite explicitly discretizing the problem. The iterative process is continued until a precise outcome of $10^{- 6}$ is achieved in every scenario. To solve equations (10)–(13) and (14), it will be reduced into the lowest order by employing the following variables.

\begin{aligned} ϖ_{1} (η) = f, ϖ_{2} (η) = f^{'}, ϖ_{3} (η) = f^{″}, ϖ_{4} = θ (η), ϖ_{5} = θ^{'} (η), ϖ_{6} = h, \\ ϖ_{7} = h^{'}, ϖ_{8} = g, ϖ_{9} = g^{'} . \end{aligned}

(16)

By inserting equation (13) in equations (9), (10) and (11), we get:

\begin{aligned} \frac{((1 + (2 δ_{1}) η) ϖ_{3}^{'} (η) + (2 δ_{1}) ϖ_{3} (η))}{B_{1} B_{2}} - (ϖ_{2} (η))^{2} - \frac{P_{m}}{B_{1} B_{2}} ϖ_{2} (η) \\ + ϖ (η)_{1} ϖ_{3} (η) + \frac{B_{3}}{B_{2}} γ_{1} θ \cos ζ = 0, \end{aligned}

(17)

(1 + R d \frac{k_{T n f}}{k_{f} P r B_{4}}) ((1 + 2 δ_{1} η) ϖ_{5}^{'} (η) + 2 δ ϖ_{5} (η)) + ϖ_{1} (η) ϖ_{5} (η) + Q_{e} \exp (- n η) .

(18)

(1 + 2 η δ_{1}) ϖ_{7}^{'} + δ_{1} ϖ_{7} + S c ϖ_{1} ϖ_{7} - S c L_{1} ϖ_{6} ϖ_{8}^{2},

(19)

δ * (1 + 2 η δ_{1}) ϖ_{9}^{'} + δ * δ_{1} ϖ_{9} + S c ϖ_{1} ϖ_{9} - S c L_{1} ϖ_{6} ϖ_{8}^{2},

(20)

The reduced BCs are:

\begin{matrix} ϖ_{1} (η) = 0, ϖ_{2} (η) = 1, ϖ_{5} (η) = - λ_{1} (1 + ϖ_{4} (η)), ϖ_{7} (η) = L_{s} ϖ_{6} (0), \\ δ * ϖ_{9} (0) = - L_{s} ϖ_{6} (0) at η = 0 \\ ϖ_{2} (\infty) = 0, ϖ_{4} (\infty) = 0, ϖ_{6} (\infty) \to 1, ϖ_{8} (\infty) \to 0 a s η \to \infty . \end{matrix}}

(21)

Results authentication

The reliability and accuracy of the results are numerically compared to the existing study in Table 4, demonstrating that the current outcomes are accurate and credible.

Table 4.

Results validation published literature for $- (f^{″} (0))$ .

P_m	Madhukesh et al. (2023)		Present work $- (f^{″} (0))$
P_m	Analytical	RK-45	bvp4c
1.0	1.4143135	1.414337	1.41433897
2.0	1.7321508	1.732151	1.73215371
3.0	2.4495897	2.449589	2.44959761
4.0	3.3167248	3.316724	3.31673920

Results and discussion

In this section, the physics behind the graphical results versus the variation of physical parameters is discussed.

Figure 2 reveals the influence of the porosity parameter $P_{m}$ on the $f^{'} (η) .$ The velocity curve of Hnf drops with the intensifying values of $P_{m}$ . Physically, the decreased permeation and fluid retention in open pores are the main causes of the drop in fluid velocity that occurs due to higher $P_{m}$ . These factors obstruct fluid flow and lower its velocity. Figure 3 exhibits that the influence of mixed convection $(γ_{1})$ enhances the velocity curve $f^{'} (η) .$ Physically, the rising effect of $γ_{1}$ determines that the buoyancy forces are more operative on the fluid flow. The velocity curve upsurges as a significance of $(γ_{1})$ . The $(γ_{1})$ drives the fluid motion, which results an enhancement in the velocity curve. Figure 4 reveals the effect of the inclination angle factor on $f^{'} (η) .$ It has been noticed that increasing the angle of cylinder or plate from $0 \circ$ to $90 \circ$ falls the velocity curve. The gravitational force is due to the rising values of $(ζ)$ , which opposing the fluid flow. $(ζ)$ is capable of changing the pressure gradient across the stream. The fluid velocity may decrease as a result of this changed pressure distribution as the flow meets different types of resistance in various directions. Figure 5 discloses the effect of Ag-Fe₃O₄-NPs on the $f^{'} (η) .$ It can be detected that the fluid velocity drops with the intensifying numbers of Ag-Fe₃O₄-NPs. Physically, the thermal conductivities of Ag-Fe₃O₄-NPs are generally higher than those of the base fluid. These NPs improve the nanofluid average thermal conductivity. A more effective exchange of heat across the Hnf is made possible by the resultant elevated thermal conductivity, which leads to more rapid cooling and drops fluid velocity.

Figure 2.

Hnf velocity $f^{'} (η)$ versus porosity parameter P_m.

Figure 3.

Hnf velocity $f^{'} (η)$ versus mixed convection $γ_{1}$ .

Figure 4.

Hnf velocity $f^{'} (η)$ versus inclination angle $ζ .$

Figure 5.

Hnf velocity $f^{'} (η)$ versus Ag-Fe₃O₄-NPs $ϕ_{1} = ϕ_{2}$ .

Figures 6 and 7 disclose the influence of Q_e and Rd on the fluid temperature. Physically, more heat gets injected into the fluid when the Q_e and Rd parameters rise. As the fluid takes in the heat from the source, this extra heat energy causes the temperature of the fluid to rise overall. Figure 8 reveals the influence of the thermal Biot number $λ_{1}$ on the temperature profile $θ (η) .$ Physically, the rate of energy distribution at a cylinder or plate's surface and the rate of conduction of heat inside the fluid are related by the Biot number. The fluid temperature upsurges due to the impact of thermal Biot number; because the slow rate of heat transferring over the plate/cylinder surface prevents the heat from quickly dissipating or absorbing, raising the temperature of the fluid.

Figure 6.

Fluid temperature $θ (η)$ versus heat source/sink parameter $Q_{e} .$

Figure 7.

Fluid temperature $θ (η)$ versus thermal radiation Rd.

Figure 8.

Fluid temperature $θ (η)$ versus thermal-Biot number $λ_{1}$ .

Figure 9 shows the influence of $γ_{1}$ (the mixed convection parameter) on the $θ (η) .$ It has been noticed that the upshot of $γ_{1}$ declines the $θ (η) .$ Physically, the rising effect of $γ_{1}$ determine that the buoyancy forces are more operative on the fluid temperature and flow. The velocity curve upsurges as a significance of the buoyancy force's ascendency. Buoyancy forces drive the fluid motion, which results in an enhancement in velocity and a dropping of energy profile $θ (η) .$ Figure 10 depicts that the energy field of Hnf increases with the mounting upshot of the $(ζ)$ . The variation in temperature of a fluid flowing through an inclined surface is influenced due to gravity. Potential energy is transformed into kinetic energy as the fluid propels upward, and friction and turbulence cause mechanical energy to be turned into thermal energy, which raises the fluid's temperature. Figure 11 exposes the effect of Ag-Fe₃O₄-NPs on the thermal profile of Hnf. Physically, the thermal conductivities of Ag-Fe₃O₄-NPs are generally higher than those of the base fluid. These NPs improve the nanofluid average thermal conductivity. A more effective exchange of heat across the Hnf is made possible by the resultant elevated thermal conductivity, which leads to more rapid cooling and a drop in temperature. In addition, the heat dissipation characteristics can be enhanced by the introduction of NPs. Heat can be absorbed and transferred by Ag-Fe₃O₄-NPs more proficiently. The enhanced capacity to dissipate heat results in a diminution in the thermal profile of the fluid.

Figure 9.

Fluid temperature $θ (η)$ versus mixed convection $γ_{1} .$

Figure 10.

Fluid temperature $θ (η)$ versus angle of inclination $ζ$ .

Figure 11.

Fluid temperature $θ (η)$ versus Ag-Fe₃O₄-NPs $ϕ_{1} = ϕ_{2}$ .

Figures 12–14 demonstrates the nature of mass profile $h (η)$ versus homogeneous reaction $L_{1}$ , heterogeneous reaction $L_{s}$ and Schmidt number Sc respectively. Every reactant in a homogeneous reaction is evenly incorporated at its atomic level. The amounts of both reactants and byproducts vary all over the fluid as the reaction proceeds. The concentration of such species goes down or up along the streams, depending on whether the reaction utilizes reactants or generates products. As the reaction progresses, this may cause the concentration curve $h (η)$ to decrease (Figure 12). A minimum of one reacting agent in a heterogeneous reaction is not in the same phase as the other reactants—solid, liquid, or gas. These reactions frequently take place at interfaces, like when reactants adhere to the outermost layer of a solid catalyst. The quantity of reactants in the larger amount of fluid falls as the reaction occurs on the surface because they have been utilized there. As the reaction progresses, this may cause the concentration profile $h (η)$ to decrease (Figure 13). Figure 14 also reveals that the mass curve drops with the variation of Sc. The molecular diffusivity falls off with the effect of Sc, whereas the kinetic viscosity of a fluid is enhanced, this results in the declination of the mass curve $h (η)$ .

Figure 12.

Mass profile $h (η)$ versus homogeneous reaction $L_{1}$ .

Figure 13.

Mass profile $h (η)$ versus heterogeneous reaction $L_{s}$ .

Figure 14.

Mass profile $h (η)$ versus Schmidt number Sc.

Figures 15 and 16 elucidate the graphical outcomes of Nusselt number $(N u)$ and Skin fraction $(C f)$ versus the Ag-Fe₃O₄-NPs. It can be perceived that the energy and velocity transference rate of hybrid nanofluid increases with the rising numbers of Ag- Fe₃O₄-NPs. This characteristic of hybrid nanofluid is of greater significance for applications in engineering and industry.

Figure 15.

Ag-Fe₃O₄-NPs versus Skin fraction (Cf).

Figure 16.

Ag-Fe₃O₄-NPs versus Nusselt number (Nu).

Conclusions

We have numerically evaluated the steady 2D Hnf flow over an inclined permeable plate/cylinder. The Hnf flow is examined under the impact of heterogeneous/homogeneous chemical reaction, and permeable medium. The Hnf is prepared by Ag and Fe₃O₄ NPs in water. The energy transmission and fluid flow are mathematically specified by using a coupled nonlinear system of PDEs. The system of PDEs is simplified into a dimensionless form of ODEs, which are then further numerically treated with the Matlab package bvp4c. The primary conclusions are:

The effect of the porosity parameter $P_{m}$ and inclination angle factor drops the velocity curve $f^{'} (η)$ in both cylinder and plate case.

The rising numbers of Ag-Fe₃O₄-NPs also decline the velocity field $f^{'} (η) .$

The influence of Q_e and heat radiation enhances the fluid temperature.

The fluid temperature upsurges due to the flourishing effect of thermal Biot number and angle inclination parameter.

The influence of mixed convection parameters $γ_{1}$ and rising numbers of Ag-Fe₃O₄-NPs drops the energy field of Hnf.

The significances of homogeneous reaction $L_{1}$ , heterogeneous reaction $L_{s}$ and Schmidt number Sc drop the mass profile $h (η)$ .

The energy transference rate of hybrid nanoliquid increases with the rising numbers of Ag-Fe₃O₄-NPs.

Footnotes

Data availability

All data used in this manuscript have been presented within the article.

Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The study was funded by Researchers Supporting Project number (RSPD2024R749), King Saud University, Riyadh, Saudi Arabia.

ORCID iDs

Khadija Rafique

Zafar Mahmood

Nomenclature

References

Abbasi

Gul

Shehzad

(2020) Effectiveness of temperature-dependent properties of Au, Ag, Fe₃O₄, Cu nanoparticles in peristalsis of nanofluids. International Communications in Heat and Mass Transfer 116: 104651.

Acharya

(2024) Hydrothermal scenario of buoyancy-driven magnetized multi-walled carbon nanotube-Fe₃O₄-water hybrid nanofluid flow within a discretely heated circular chamber fitted with fins. Journal of Magnetism and Magnetic Materials 589: 171612.

Adnan Abbas

Alqahtani

, et al. (2024) Numerical heat featuring in radiative convective ternary nanofluid under induced magnetic field and heat generating source. International Journal of Modern Physics B: 2550044.

Adnan Abbas

Said

, et al. (2024) Significance of coupled effects of resistive heating and perpendicular magnetic field on heat transfer process of mixed convective flow of ternary nanofluid. Journal of Thermal Analysis and Calorimetry 149(2): 879–892.

Ahmad

, et al. (2024) Energy and mass transmission through hybrid nanofluid flow passing over a spinning sphere with magnetic effect and heat source/sink. Nanotechnology Reviews 13(1): 20230194.

Akbar

Zamir

Muhammad

(2024) Levenberg-Marquardt technique analysis of thermal and concentration storage in cone-disk apparatus with neural network-enhancement. Thermal Science and Engineering Progress 50: 102529.

Akbar

, et al. (2024) Simulation of hybrid boiling nano fluid flow with convective boundary conditions through a porous stretching sheet through Levenberg-Marquardt artificial neural networks approach. International Journal of Heat and Mass Transfer 228: 125615.

Algehyne

, et al. (2022) Numerical simulation of bioconvective Darcy Forchhemier nanofluid flow with energy transition over a permeable vertical plate. Scientific Reports 12(1): 3228.

Alghamdi

, et al. (2024) Double layered combined convective heated flow of Eyring-Powell fluid across an elevated stretched cylinder using intelligent computing approach. Case Studies in Thermal Engineering 54: 104009.

10.

Alharbi

KAM

, et al. (2022) Computational valuation of Darcy ternary-hybrid nanofluid flow across an extending cylinder with induction effects. Micromachines 13(4): 588.

11.

Ali

Jubair

Mahmood

, et al. (2024) Analysis of unsteady MHD fluid flow across two parallel discs with uniform fluctuation subject to modified Hall and activation energy. International Journal of Modern Physics B: 2450424.

12.

Ali

Jubair

Siddiqui

MIH

. (2024) Numerical simulation of 3D Darcy-Forchheimer hybrid nanofluid flow with heat source/sink and partial slip effect across a spinning disc. Journal of Porous Media 27: 97–111.

13.

Ali

, et al. (2023) Exploring the impact of Hall and ion slip effects on mixed convective flow of Casson fluid model: A stochastic investigation through non-Fourier double diffusion theories using ANNs techniques. Thermal Science and Engineering Progress 46: 102237.

14.

Alqahtani

, et al. (2023) Energy transmission through Carreau Yasuda fluid influenced by ethylene glycol with activation energy and ternary hybrid nanocomposites by using a mathematical model. Heliyon 9(4).

15.

Bani-Fwaz

, et al. (2024) Computational investigation of thermal process in radiated nanofluid modulation influenced by nanoparticles (Al₂O₃) and molecular (H₂O) diameters. Journal of Computational Design and Engineering 11: 22–36.

16.

Bilal

, et al. (2021) Darcy-Forchheimer hybrid nano fluid flow with mixed convection past an inclined cylinder. CMC-Computers Materials & Continua 66(2): 2025–2039.

17.

Bin Mizan

, et al. (2021) Computation of ferromagnetic/nonmagnetic nanofluid flow over a stretching cylinder with induction and curvature effects. Heat Transfer 50(6): 5240–5266.

18.

Chamkha

(1997) Non-Darcy fully developed mixed convection in a porous medium channel with heat generation/absorption and hydromagnetic effects. Numerical Heat Transfer, Part A Applications 32(6): 653–675.

19.

Chamkha

Rashad

(2012) Natural convection from a vertical permeable cone in a nanofluid saturated porous media for uniform heat and nanoparticles volume fraction fluxes. International Journal of Numerical Methods for Heat & Fluid Flow 22(8): 1073–1085.

20.

Damseh

, et al. (2009) Combined effect of heat generation or absorption and first-order chemical reaction on micropolar fluid flows over a uniformly stretched permeable surface. International Journal of Thermal Sciences 48(8): 1658–1663.

21.

Ghailan

, et al. (2023) Biological analysis of emerging nanoparticles with blood through propagating flow along a plumb porous canal in the occurrence of energy and heat transfer. Surfaces and Interfaces 40: 103013.

22.

Haq

, et al. (2022) Impact of homogeneous and heterogeneous reactions in the presence of hybrid nanofluid flow on various geometries. Frontiers in Chemistry 10: 1032805.

23.

Islam

Mahmood

Khan

, et al. (2024a) Joule heating effects on triple diffusive free convective MHD flow over a convective surface: A lie-group transformation analysis. Modern Physics Letters B: 2450412.

24.

Islam

Mahmood

Khan

, et al. (2024b) Lie scaling transformations for the analysis of MHD flow with radiation, Soret, and Dufour effects as well as viscous dissipation across a convective surface during triple diffusion. International Journal of Heat and Fluid Flow 108: 109465.

25.

Madhukesh

, et al. (2023) Investigation of thermal performance of ternary hybrid nanofluid flow in a permeable inclined cylinder/plate. Energies 16(6): 2630.

26.

Mahmood

Rafique

Khan

, et al. (2024a) Heat transfer in radiative hybrid nanofluids over moving sheet with porous media and slip conditions: Numerical analysis of variable viscosity and thermal conductivity. Materials Today Communications 40: 109664.

27.

Mahmood

Rafique

Khan

, et al. (2024b) Importance of thermal conductivity models in analyzing heat transfer of radiative hybrid nanofluid across a stretching sheet using Darcy-Forchheimer flow. Journal of Porous Media 27(7): 1–24.

28.

Maraj

, et al. (2023) Thermally progressive particle-Cu/blood peristaltic transport with mass transfer in a non-uniform wavy channel: Closed-form exact solutions. Alexandria Engineering Journal 74: 453–466.

29.

Muhammad

, et al. (2021) Melting heat transfer in squeezing flow of basefluid (water), nanofluid (CNTs + water) and hybrid nanofluid (CNTs + CuO + water). Journal of Thermal Analysis and Calorimetry 143: 1157–1174.

30.

Naveed

Imran

Gul

(2023) Heat transfer analysis in hydromagnetic flow of couple stress fluid in presence of homogeneous and heterogeneous chemical reactions over a porous oscillatory stretchable sheet. Advances in Mechanical Engineering 15(2): 16878132231155824.

31.

Pattnaik

Bhukta

Mishra

(2022) Experience of non-uniform heat source/sink on the flow of micropolar nanofluid fluid with heterogeneous and homogenous chemical reaction. Waves in Random and Complex Media: 1–19.

32.

Pooja

Narasimhamurthy

Anitha

(2024) Numerical investigation of flow dynamics of Williamson fluid over an expanding cylinder/plate in presence of homogeneous/heterogeneous reactions. International Journal of Applied and Computational Mathematics 10(1): 10.

33.

Rafique

Mahmood

Adnan , et al. (2024) Numerical analysis of non-linear radiative Casson fluids containing CNTs having length and radius over permeable moving plate. Open Physics 22(1): 20240013.

34.

Rafique

Mahmood

Ali

, et al. (2024) Entropy analysis of Hall effect with variable viscosity and slip conditions on rotating hybrid nanofluid flow over nonlinear radiative surface. Materials Today Communications 109167.

35.

Raizah

, et al. (2023) Parametric simulation of stagnation point flow of motile microorganism hybrid nanofluid across a circular cylinder with sinusoidal radius. Open Physics 21(1): 20220205.

36.

Salawu

, et al. (2023) Elastic deformation of thermal radiative and convective hybrid SWCNT-Ag and MWCNT-MoS4 magneto-nanofluids flow in a cylinder. Results in Materials 17: 100380.

37.

Shah

Ali

(2019) Applications of hybrid nanofluids in solar energy, practical limitations and challenges: A critical review. Solar Energy 183: 173–203.

38.

Shamshuddin

, et al. (2024) Chebyshev spectral approach to an exponentially space-based heat generating single-phase nanofluid flowing on an elongated sheet with angled magnetic field. Numerical Heat Transfer, Part B: Fundamentals 85(2): 159–176.

39.

Takhar

Chamkha

Nath

(2000) Combined heat and mass transfer along a vertical moving cylinder with a free stream. Heat and Mass Transfer 36: 237–246.

40.

Vijatha

Bala Anki Reddy

(2024) ‘Second law analysis on magnetohydrodynamics flow of Eyring-Powell hybrid nanofluid over a cylinder/plate with non-linear thermal radiation and Cattaneo-Christov heat flux: Numerical and semi-analytical simulations. Numerical Heat Transfer, Part B: Fundamentals: 1–29.

41.

Waqas

, et al. (2021) ‘On bio-convection thermal radiation in Darcy–Forchheimer flow of nanofluid with gyrotactic motile microorganism under Wu’s slip over stretching cylinder/plate’. International Journal of Numerical Methods for Heat & Fluid Flow 31(5): 1520–1546.

42.

Waqas

, et al. (2024) Numerical analysis of MHD tangent hyperbolic nanofluid flow over a stretching surface subject to heat source/sink. Pramana 98(1): 27.

43.

Xie

, et al. (2024) Simulation of CO₂ dissolution reactions in saline aquifers using lattice Boltzmann method. Gas Science and Engineering 125: 205284.

44.

Y-J

, et al. (2021) Gyrotactic micro-organism flow of Maxwell nanofluid between two parallel plates. Scientific Reports 11(1): 15142.

45.

Yasmin

, et al. (2024) Optimized framework numerical solution for swirling hybrid nanofluid flow with silver/gold nanoparticles on a stretching cylinder with heat source/sink and reactive agents. Open Physics 22(1): 20230202.

46.

Zhang

X-H

, et al. (2021) The parametric study of hybrid nanofluid flow with heat transition characteristics over a fluctuating spinning disk. Plos one 16(8): e0254457.