An empirical assessment of Tukey combined extended exponentially weighted moving average control chart

1.   Introduction

The suspension of ternary hybrid nanofluid is composed of three kinds of nanoparticles suspended in a base fluid (such as ethylene glycol) to improve its viscosity, thermal conductivity, and other thermophysical characteristics. These fluids are extensively researched in heat transfer applications, such as renewable energy systems, cooling systems, and biomedical applications. A finite element scheme is used to analyze the fluid flow, heat transfer, and other properties of ternary hybrid nanofluids. FEM is especially helpful for resolving non-linear problems and complex geometries. Sohail et al. ^[1] investigated the results of tetra-hybrid nanofluid towards a vertical disk. The induced magnetic field was utilized. They used non-Fourier's law and thermal radiation (non-linear). The numerical results have been done by finite element scheme. Shamshuddin et al. ^[2] used the HAM method to study a rotating disk system with tri-hybrid nanoparticles being affected by magnetic fields, viscous dissipation, and Hall currents. According to the results, tri-hybrid nanofluids have better velocity profiles than single (Cu) and hybrid (Cu-FeO₄) nanofluids. By adding nanoparticles, the base fluid's energy transfer efficiency is greatly increased. Shamshuddin et al. ^[3] studied the HAM to investigate how different chemical reactions affect the flow of hybrid nanofluids that incorporate Prandtl nanoparticles and Ohmic heating. Results show that lowering molecular diffusivity by increasing the chemical reaction parameter results in a decrease in concentration field profiles. This demonstrated that the interaction between reaction dynamics and fluid behavior is enhanced by nanoparticles. Shamshuddin et al. ^[4] worked on ternary hybrid nanofluids, like MWCNT-Au-Ag/blood mixtures. They included that THNF has various applications like atomic reactors, solar energy systems, heat pipes, cars, and biomedical engineering. The flow of such nanofluids over a bidirectional stretching sheet is investigated in this work, which is motivated by developments in nanotechnology. Ali et al. ^[5] focused on heat radiation and Soret effects over an inclined stretching sheet, this work explored nanofluid flow within a micropolar MHD framework. The LMB optimization technique is used to solve the governing ODEs. The analysis emphasizes how thermal and concentration gradients interact in sophisticated fluid mechanics. Khan et al. ^[6] studied on neural network-based optimization technique called the Levenberg–Marquardt (LM) backpropagation algorithm to numerically assess electromagnetohydrodynamic bioconvection in micropolar nanofluids. To analyze the flow behavior, the study takes into account stratification effects and thermal radiation. The method has described how effective sophisticated computational methods are in complex fluid dynamics. Alqahtani et al. ^[7] discussed energy enhancement in THB Carreau Yasuda nanofluid flow with a magnetic dipole across a vertical sheet. The fluid, was enriched with MWCNTs, Zn, and Cu in ethylene glycol using activation energy, heat source/sink, chemical reactions, and Darcy-Forchheimer effects. Practical effects like Joule heating and a magnetic field are incorporated into the flow analysis by Rahman et al. ^[8]. This investigates a hybrid nanofluid's two-dimensional flow in a porous medium enclosed by two parallel walls. Jan et al. ^[9] discussed the study of thermal radiation with a cross-model ternary hybrid nanofluid's laminar flow. The flow takes place in a porous medium under a magnetic field and is impacted by a stretching cylinder. Rahman et al. ^[10] examined Darcy-Brinkman flow across a stretched sheet while taking frictional heating and porous dissipation into account. The RK4 method was used to analyze the steady flow of a dusty fluid through a porous medium while accounting for slip effects. Rehman et al. ^[11] explored Darcy-Forchheimer flow across a porous curved stretching surface. To account for high-velocity effects in porous media, the model integrates Darcy's law with Forchheimer corrections. The governing equations are solved using HAM. Khan et al. ^[12] investigated chemically reactive exothermic fluid flow through a porous medium under gravity with MHD. The modified Darcy law was used to model the resistance of the porous medium. To solve the equations, a semi-implicit (FDA) has been used. Hayat et al. ^[13] estimated the peristaltic flow of MHD hybrid nanofluids in a symmetric channel. The Darcy–Forchheimer law is used to analyze the nonlinear behavior of porous media. They examined the effects of viscous dissipation and Joule heating. Adun et al. ^[14] studied a novel ternary hybrid nanofluid's thermophysical characteristics for use in energy applications. Its performance was optimized through experimental techniques. The results were intended to increase practical systems' energy efficiency. Lone et al. ^[15] utilized bvp4c method to solve the transformed equations. The ternary nanofluid's thermal efficiency can be increased by several factors, including thermophoresis, thermal radiation, Biot number, Eckert number, and magnetic field. The performance of heat transfer was greatly influenced by these factors. Khan et al. ^[16] examined the flow of ternary hybrid nanofluids over a stretching surface while considering porosity, magnetic effects, heat generation, and viscous dissipation. The ANN approach was used in a machine learning approach. Tanuja et al. ^[17] experienced the LTNE model used to investigate heat transfer in wetted rectangular porous fins. The Darcy model and Boussinesq approximation were used to analyze buoyancy with forces, and convection and radiation effects. It was estimated that the temperature profiles of the solid and ternary nanofluid phases were improved by lowering the surface-ambient radiation parameter. This suggests that heat transfer performance has been improved. Kim et al. ^[18] determined that heat transfer and magnetohydrodynamic natural convection in a ternary-hybrid nanofluid inside a partially heated porous square cavity are the main topics of this work. Additionally, the effects of thermal radiation and heat generation and absorption were examined. Mohanty et al. ^[19] estimated a cross-ternary hybrid nanofluid containing the nanoparticles SWCNT, MWCNT, and GO around a vertical cylinder. The thermophysical behavior of the nanoparticles was highlighted by the significant influence of Joule heating and radiation on thermal performance. Entropy generation is examined using the HAM. Ullah et al. ^[20] examined the thermal behavior of tri-hybrid nanofluid flow over an expanding disk, incorporating nanoparticles in a water-based solution. They investigate heat transfer and fluid dynamics under the influence of quadratic convection, ohmic heating, and a magnetic field. HAM is used to provide insights into enhanced thermal performance. Mohanty et al. ^[21] experienced the study of tri-hybrid nano-particles on a rotating disk with statistical analysis. They employed a finite element scheme for numerical outcomes. They also included that platelet-shaped nanofluid was more significant than other shapes of nanofluid for the case of heat transfer purposes. Faizan et al. ^[22] studied the mechanism of ternary hybrid nanofluid on the surface including Williamson fluid. They estimated entropy generation with thermal radiation and heat sink with convective boundary conditions. Ramzan et al. ^[23] discussed the role of Carreau fluid in multiple nanoscales under the influence of MHD with entropy generation. They estimated that maximum heat energy using ternary hybrid nano-fluid. Mishra et al. ^[24] discussed the novelty work of the machine learning approach on heat transfer inclusion of tri-hybrid nanofluid. They took three geometries considering heat sink, non-Fourier's law, and thermal radiation. They utilized bvp4c. Mishra ^[25] studied waste discharge with the hybrid nanofluid on a stretching/shrinking surface utilizing a chemical reaction. Mishra has worked on Darcy's law using the bvp4c approach. Mishra included that concentration profiles enhance with increment of unsteadiness parameter. Mishra ^[26] observed the significant role of hybrid nanofluid with base fluids (ethylene glycol and water) on a permeable plate. Thompson and Troian slip influences were considered. The numerical results have been obtained using the bvp4c approach. Mishra ^[27] discussed the study of gyrotactic microorganisms with Xue and Yamada-Ota models on an exponential sheet. Brownian motion, viscous dissipation, and radiation parameters are considered. Ramasekhar et al. ^[28] utilized a hybrid nanofluid in Casson fluid with a heat source and thermal radiation (non-linear). An approach called bvp5c has been utilized for numerical purposes. Basha et al. ^[29] discussed the characterizations of Casson fluid with thermal radiation. The numerical consequences have been attained using OHAM. Chabani and Mebarek-Oudina ^[30] estimated the role of a hybrid nanofluid utilizing a finite element approach. Mezaache et al. ^[31] visualized thermal energy in a heat exchanger channel. They discussed the influence of entropy generation. They studied Darcy–Forchheimer law in nanofluid with a thermal equilibrium approach.

The Xue model is frequently used to predict the effective thermal conductivity of nanofluids that contain cylindrical nanoparticles. However, applications of the Xue model can be extended to various nanoparticles based on $CuO$, graphene oxide and ${Al}_{2}{O}_{3}$ in ethylene glycol. The thermal conductivity of $CuO$ and ${Al}_{2}{O}_{3}$ nanoparticles dissolved in ethylene glycol was examined in a study ^[32]. Sundar et al. ^[32] showed that both kinds of nanoparticles improved the base fluid's thermal conductivity, confirming the predictive value of models such as Xue's. Zhang et al. ^[33] studied the Xue and Maxwell models for hybrid nanofluids. They showed that the Xue model can be applied to different nanoparticle compositions and shapes, such as $CuO$ and ${Al}_{2}{O}_{3}$. Zhang et al. ^[34] discussed the application of the Xue model to predict the thermal behavior of nanofluids containing ethylene glycol as the base fluid. They observed that the Xue model has been effectively used to evaluate thermal conductivity in studies involving the peristaltic flow of hybrid nanofluids based on ethylene glycol. Tri-hybrid nanofluid is the suspension of three nanofluids in base fluid (ethylene glycol). Tri-hybrid nanofluid has the ability of heat transfer in comparison to base fluids like ethylene glycol, oil and water, etc. The various models of Xue, Yamada-Ota, and Hamilton-Crosser are studied for heat transfer enhancement. The Yamada Ota-model is used for predicting the volume fractions and particle size due to thermal conductivity. Xue model is called the two-phase model that is utilized in such processes where heat transfers in the solid phase and fluid phase. Both models have been used in cooling systems, solar thermal collectors, engine cooling fluids, drug delivery processes, photovoltaic cooling, and HVAC systems. Galal et al. ^[35] discussed the models (Xue and Yamada Ota) in Casson fluid considering various nanofluid over-stretching 3D surfaces via local thermal (non-equilibrium). Rafiq et al. ^[36] discussed comparison analysis between Yamad Ota and Xue, including hybrid nanofluid on disk via Lorentz force.

Heat exchangers, nuclear power plants, jet engine cooling, electronic cooling, solar energy systems, drug delivery systems, and metal cutting/machining fluids are the applications of ternary hybrid nanofluids that are depicted in Figure 1.

Table 1 describes the comparative novelty analysis with published works ^{[37,38,39,40]}. It is noticed that tri-hybrid nanofluid, Darcy-Forchheimer law, heat source, and viscous dissipation on a vertical jet are not studied in the current model. Further, the finite element approach is used while the Box–Behnken design is utilized. It is claimed that novelty work is discussed in the section on the formulation of the problem. Applications for the study of tri-hybrid nanofluids in a vertical jet using the finite element method and Darcy-Forchheimer law include industrial heat exchangers, thermal energy storage, and sophisticated cooling systems. The Box-Behnken design enhances heat transfer efficiency, supporting microfluidic devices, biomedical cooling, and aerospace propulsion.

2.   Modeling of the problems

A 2D flow of Casson fluid with steady motion is considered on an oriented round jet along with axial symmetry. The coordinate system is in line with the direction normal to the jet and the axis of symmetry. Except for density, which drops exponentially with temperature, the fluid's characteristics are constant and it is incompressible. It is assumed that the centerline temperature is much higher than the fluid's ambient temperature. The effects of Casson fluid behavior and tri-hybrid nanofluid flow are considered. The Darcy-Forchheimer law, which takes into consideration the resistance of porous media, controls the jet's flow. To simulate energy dissipation, a heat sink is incorporated. The geometry is described in Figure 2. It is noticed that ethylene glycol is discovered as a non-Newtonian ^[41] and shear thinning ^[42]. The rheology of Casson fluid is non-Newtonian. Therefore, a combined study of Casson fluid with ethylene glycol is discussed in this model since such phenomena are used in energy-efficient thermal transport, biomedical fluids, and cooling systems. Thus, due to its superior thermal qualities such as its high thermal conductivity, low freezing point, and chemical stability, ethylene glycol is selected as the base fluid for heat transfer applications.

Figure 2 describes the diagram of a vertical jet, and the coordinate system is adjusted as the x-axis is taken as a normal direction while the y-axis is taken as a horizontal direction. The motion of the nanofluid on a vertical jet is produced using wall velocity. It is noticed that velocity components $u$ are along the y-axis and $v$ is taken along the x-axis. The flow is influenced by gravity ($G$). A cross-section area is denoted by tangential velocity $w$. The ambient temperature is assumed as ${t}_{\infty }.$ The centerline temperature of the jet is much higher than the ambient fluid temperature ${t}_{\infty }$, resulting in a thermal gradient that affects fluid motion and heat transfer as the jet expands vertically. Applications of vertical jets are cooling systems and fuel injection. It improves cooling in gas turbines, electronics, and material processing by causing a high-velocity jet to strike a surface. Optimizing industrial spraying, biomedical drug delivery, and propulsion systems through an understanding of velocity, pressure, and thermal boundary layers in such jets enhances performance and efficiency in engineering applications. The following assumptions are listed below.

➢ A steady flow of Casson fluid is considered;

➢ A vertical jet is used;

➢ Darcy's and Forchhiermer's law are discussed;

➢ An incompressible and two-dimensional flow are considered;

➢ Viscous dissipation and heat sink are utilized;

➢ Variable thermal conductivity is used;

➢ Xue and Yamada Ota hybrid nanofluid are inserted.

A nonlinear rheological model known as the Casson fluid model explains how viscoelastic fluids behave when they experience yield stress. The stress tensor of Casson fluid ^[43] is defined as

In Eq (1), ${\pi }_{e}$ is formulated using a component of deformation rate, ${\pi }_{e}$ is called critical value, $\mu$ is termed as plastic dynamic viscosity and ${P}_{y}$ is called yield stress. It is noticed that the stress tensor of Casson fluid for $n = 1$ is defined ^[43] as

Here, $A$ is known as the Rivlin-Ericksen tensor. It is noted that Eq (2) can be reduced into Casson fluid when ${P}_{y} = 0.$ The model of PDEs ^{[37,38,39,40]} is defined as

Equation (3) describe continuity equation for the case of steady and incompressible flow. Equation (4) describes the momentum equation in the presence of Casson fluid. The term related to $u\frac{\partial u}{\partial x}+u\frac{\partial u}{\partial x}$ is called convective acceleration, ${\mu }_{thb}\left(1+\frac{1}{\beta }\right)\frac{{\partial }^{2}u}{\partial {y}^{2}}$ is viscous dissipation in the presence of Casson fluid, $G\left({\rho }_{\infty }-{\rho }_{thb}\right)$ reveals the buoyancy force because of density variations, $-\frac{{\mu }_{thb}}{k}u-\frac{F}{{k}^{\frac{1}{2}}}{u}^{2}$ are formulated due to Darcy's Forchhiermer law. Equation (5) is called energy equation, $\left(u\frac{\partial t}{\partial x}+u\frac{\partial t}{\partial x}\right)$ is convective heat, ${\mathcal{K}}_{thb}\frac{1}{y}\frac{1}{\partial y}\left(y\frac{\partial t}{\partial y}\right)$ is diffusion of heat due variable thermal conductivity, $Q\left(t-{t}_{\infty }\right)$ is the heat source and $\frac{{\mu }_{thb}}{2}{\left(\frac{\partial u}{\partial y}\right)}^{2}$is the viscous dissipation term.

BCs ^[37] are defined as

Now, dimensionless parameters are defined ^[37] as

Dimensionless ODEs ^[37] are

Dimensionless BCs ^[37] are

The model for the ternary hybrid nanofluid ^[42] is provided as

The ternary hybrid Xue model ^[44] is described as

The Yamada Ota model (ternary hybrid) ^[44] is derived as

The variable density parameter, fluid's absorption number, Prandtl number, Eckert number, Darcy's number, Forchhiermer number, and heat source number are defined as

Skin friction coefficient for Newtonian fluid ^[37] is

Nusselt number ^[37] is

3.   Numerical solution

The above equations, which are created using primitive variables, can be solved using the Finite Element Method. A detailed description of how to use the finite element method to solve the system of PDEs (partial differential equations) is provided here. First, primitive variables are utilized for easy integration of the problem. Finite element approach is significant rather analytical approach for the case of handing complex problems.

➢ For complex geometries like porous media, spherical, cylindrical domains, and non-uniform surfaces, the finite element scheme works incredibly well;

➢ Such nonlinear PDEs are difficult for analytical approaches to solve, but the finite element scheme effectively manages them by discretizing the domain into discrete elements.

➢ By establishing conditions on the nodes of elements, the finite element scheme enables the flexible implementation of boundary conditions.

Here, primitive variables are defined as

Using a primitive variable in Eqs (8)–(10). Now, equations are

Weak Formulations: Integrate over the domain Ω after multiplying each PDE by a test function from the solution space. After multiplying each PDE by a test function from the solution space, integrate over domain Ω.

Discretization of domain: Segment the domain Ω into finite components, such as triangle elements.

Calculations of Stiffness elements: The calculations of stiffness elements are mentioned here

The Global stiffness matrix is derived as

Assembly Process: Use the weak form of the basic functions to calculate the element stiffness matrix and element load vector for each element. The global matrix is

Solve the system of equations: Use the proper numerical techniques to solve the linear or nonlinear algebraic system, such as the Gaussian elimination approach. Verifying that the maximum relative change between two successive iteration levels falls below 10 at each stage is known as the convergence of the solution. The flow chart of the numerical approach is mentioned in Figure 3. Table 3 describes the comparative results of heat transfer distribution with a published study ^[37] when ${Q}^{*} = 0, Fr = ϵ = 0, Ec = 0, {N}_{3} = {N}_{1} = 0$, and ${N}_{2} = 0.$

Error analysis: It is noticed that error analysis is derived as

and convergence condition is

4.   Results and discussion

The two-dimensional motion of steady Casson fluid is considered on a jet. The suspension of ternary hybrid nano-fluid is assumed. Darcy's Forchhiermer law and heat sink are considered. ANOVA of Box–Behnken Design is estimated for Nusselt number and divergent velocity. The chosen parameter ranges guarantee accurate modeling of fluid flow and heat transfer. The ranges of dimensionless parameters are $0.0\le D\le 4.0, 0.0\le Ec\le 5.0, -5.0\le {Q}^{*}\le 4.5$, $0.0\le D\le 7.0, 0.0\le m\ge \mathrm{2.5, 0.0} \le Fr\le 3$, $0.0\le ϵ\ge 7.0, 0.0\le {N}_{3}\le 0.006, {0.0\le N}_{1}\le 0.007$, ${0.0\le N}_{2}\le \mathrm{0.005, 0.0}\le \beta \le 3.0$, $5\le Pr\ge 12.$The detailed explanations are mentioned below.

In Figures 4–7, green lines are associated with THBN, HN is described by hybrid nanofluid, and NF describes the role of blue lines. Figure 4 estimates the variation of $Fr$ on $U$ in the presence of tri-hybrid nano-particles. With an increment of $Fr$, the velocity field declines. Darcy-Forchheimer law combines inertial and viscous effects to describe fluid flow through porous media. The Forchheimer term is added to Darcy's law, which applies to slow, viscous flow to account for higher flow velocities where inertial effects become noticeable. In porous media, inertial resistance is produced due to the Forchheimer parameter. It causes a rise in the nonlinear drag force, which declines fluid motion and lowers velocity, which lowers the flow rate. THBN on motion is higher than NF and HB on motion. The variation of $\mathrm{m}$ on the velocity field is captured in Figure 5. The curves of $\mathrm{U}$ are inclined when $\mathrm{m}$ is enhanced. Here, $\mathrm{m}$ is called the power law number. In such models, the power-law number m can describe how the resistance to flow varies with velocity or shear rate. Fluid density variations are amplified by increasing the variable density parameter, which raises inertia and motion resistance. This causes the velocity field to decrease and the flow to slow down due to a higher momentum diffusion rate. Figure 7 captures the behavior of $D$ on velocity curves. The velocity field increases when $D$ is enhanced. Here, $D$ is known as the fluid's absorption number. It is a dimensionless figure that expresses the quickly a fluid (or energy) is absorbed into a medium about other transport effects. The role of $Ec$ on the temperature profile is visualized in Figure 8 considering THBNF, HNF, and nanofluid. The concept of viscous dissipation is utilized for modeling of Eckert number. When the Eckert number is enhanced, more heat energy is dissipated. The temperature profile is enhanced. The amount of heat energy for nanofluid is less than the heat energy for THBNF and HNF. Figure 9 exhibits the trend of the heat sink on heat including THNF, HNF and nanofluid. The heat sink is denoted by ${Q}^{*}$ which is proportional to temperature difference. An increment of heat sink deals with an enhancement of heat energy. The heat source parameter $\left({Q}^{*}\right)$ controls both heat absorption (negative ${Q}^{*}$) and heat generation (positive ${Q}^{*}$). It has a major effect on heat transfer in a number of applications, such as industrial heat exchangers, electronic cooling, and nuclear reactors. The change of $D$ on heat energy is estimated in Figure 10 considering nanofluid, HNF and THNF. Thermal energy is boosted when $D$ is enhanced. Here, $D$ is called the fluid's absorption number. The production of heat energy for nanofluid is less than heat energy for HNF, and THNF. The variation of $m$ on the velocity curve is captured in Figure 11 including THNF, NF, and hybrid nanofluid. The velocity curve is enhanced with the change of $m.$ Convective heat transfer efficiency is decreased when fluid density variations are enhanced by an increase in the variable density parameter. As a result, the fluid retains more heat, raising the temperature field and thermal energy distribution. The variation of $Fr$ and $ϵ$ on divergent velocity is estimated in Figure 12 considering THNF, HNF, and nanofluid. Divergent velocity increases when $Fr$ and $ϵ$ are increased. Resistance to fluid expansion is decreased by increasing inertial and nonlinear effects while increasing $Fr$ and $ϵ.$ This encourages flow divergence, which raises divergent velocity because of stronger inertial forces. The variation of $D$ and $m$ on Cf is tackled in Figure 13. The skin friction coefficient is declined when $D$ and $m$ are increased. While increasing $m$ amplifies density variations and encourages flow expansion, increasing $D$ improves fluid absorption and lowers local resistance. Therefore, divergent velocity decreases. The production of tri-hybrid nanofluid is greater than the production of HNF and nanofluid. Figure 14 shows the roles of $m$ and ${Q}^{*}$ on the Nusselt number. It is observed that the Nusselt number enhances when $m$ and ${Q}^{*}$ are increased. The variation of $D$ and $M$ on the Nusselt number is observed in Figure 15 for the case of hybrid nanofluid and tri-hybrid nano-fluid. The heat transfer is improved when the buoyancy-driven convection is enhanced when $m$ is enhanced. An increasing in ${Q}^{*}$ supplies more heat energy, which enhances the temperature gradients. Therefore, Nusselt number increases. Figure 16 describes the change of $m$ and ${Q}^{*}$ on the Nusselt number for the Xue and Yamada Ota ternary hybrid nanofluid model. It is noticed that the heat transfer rate for Yamada Ota ternary hybrid nanofluid is higher than Xue-ternary hybrid nanofluid. It is used because of their precision in predicting thermal conductivity and effective viscosity in intricate multiphase flows, the Xue and Yamada-Ota hybrid nanofluid models were chosen. Yamada-Ota's model is perfect for high-performance cooling applications because it more accurately depicts the interactions between nanoparticles, which results in higher heat transfer rates. Despite its conservatism, Xue's model is still applicable to stability analysis. Their comparison ensures better thermal management by offering insights into how to optimize hybrid nanofluids for energy and industrial applications. The production of TNF is greater than HNF. Table 4 describes the variation of $Ec, m$ and $D$on Nusselt number for HNF and THNF. The production of heat transfer rate for THNF is less than HNF. Nusselt number declines when $Ec$ is enhanced. The rate of heat transfer decreases when $m$ is enhanced. Heat transfer rate decreases with the change of $D.$ The production of Nusselt number THNF is greater than HNF and nanofluid.

Box–Behnken design method: Nusselt number and divergent velocity are determined with the approach of the Box–Behnken design method. The Box–Behnken design method (see Figure 17) is a statistical technique for response variable optimization that is used to analyze the Nusselt number and divergent velocity. This approach systematically assesses the effects of several variables and how they interact. It effectively ascertains the ideal parameters for improved fluid system flow and heat transfer. The parameters and factor levels that are utilized in the Box–Behnken Design (BBD) for the analysis of skin friction coefficients (SKF) are listed in Table 5. The variables are the magnetic parameter (B, C) between 0.3 and 0.9, the power-law index (m, B) between 1 and 6, and the Forchheimer number (Fr, A) between 0.1 and 0.8. These ranges enable a thorough assessment of how they affect the responses. The parameters and factor ranges utilized in the Box–Behnken Design (BBD) for Nusselt number analysis are shown in Table 6. The variables include the diffusion parameter (D, C) between 0.1 and 0.8, the heat source parameter (Q*, B) between 2 and 8, and the Eckert number (Ec, A) between 0.3 and 0.9. These ranges aid in the investigation of their impact on heat transfer. BBD matrix for calculating the skin friction coefficient ($cf$) with factors $Fr$, $m$, and $B$ at different levels is displayed in Table 7. The 15 runs in the table each have unique parameter combinations and matching $Cf$ responses. The ANOVA results for the skin friction coefficient (Cf) are displayed in Table 8, highlighting the importance of the model terms. The linear term m (30.70%) and interaction effects like Fr×m (22.17%) are the major contributors to the model's 94.98% explanation of the variation, whereas the error only makes up 2.70 percent. BBD matrix for calculating the Nusselt number (Nu) with factors $Ec$, ${Q}^{*}$, and $D$ at different levels is shown in Table 9. It includes 15 runs that demonstrate particular parameter combinations and the Nu responses that correspond to them to evaluate the behavior of heat transfer. The Nusselt number (Nu) analysis is summarized in the ANOVA table, which also displays the significance of the model with an F-value of 3.22 and a P-value of 0.021. Quadratic terms ($Ec$: P = 0.029) and interactions such as Ec*D (P = 0.002) make significant contributions, demonstrating their influence on heat transfer. According to the ANOVA table 10, linear, square, and interaction terms all contribute to the model's significance (P = 0.021). Low P-values for important parameters like Q*, D, and their interactions show how strongly they affect the Nusselt number. Figure 19 depicts the residuals roughly following a straight line, which indicates normality, according to the Normal Probability Plot. Despite the residuals clustering for small fitted values, the Versus Fits plot shows no discernible pattern, indicating homoscedasticity. The model's assumptions are supported by the Histogram and Versus Order plots, which display a random residual distribution devoid of trends. The normality of the residuals is roughly indicated by the normality plot, with minor deviations at the extremes (see Figure 18). There are no discernible patterns in the Versus Fits and Versus Order plots, indicating constant variance and independence. The residuals' normal distribution is validated by the histogram, which supports symmetry. Figures 20a, b, c describe the role of contour plots of $Fr, M, B$ on the Skin friction coefficient. Figures 21a, b, c describe contour plots of $Ec, {Q}^{*}, D$ on the Skin friction coefficient. The change of ${Q}^{*}$ and $Ec$ on the Nusselt number is visualized in Figure 13. The Nusselt number is a decline function with the change of ${Q}^{*}$ and $Ec.$

The Regression equation of Nusselt number and divergent velocity is defined as

The limitations of BBD are listed below.

➢ According to BBD, the output responses and input parameters have a quadratic relationship. For highly non-linear problems, BBD might not be able to capture intricate relationships.

➢ BBD needs many simulations as the number of input factors rises;

The computational cost rises dramatically for problems with multiple parameters.

5.   Conclusions

Two-dimensional Casson rheology with Darcy's Forchhiermer law is utilized in Jet. The correlation of THNF is utilized with heat sink and viscous dissipation. The finite difference approach is used for numerical outcomes. The key points are listed below.

➢ The motion of nanofluid is enhanced when Forchiermer number and $B$ are enhanced. However, the motion of the nanofluid is decreased when $m$ and $D$ are enhanced;

➢ The motion of THNF is greater than the motion of HNF and nanofluid;

➢ Thickness related to TBLT for THNF is less than the thickness of TBLT of HNF and nanofluid;

➢ The temperature profile is boosted when heat sink and Eckert number, $D$, and the power law index number are enhanced;

➢ THNF has remarkable achievement in heat transfer rate on the divergent velocity rather than HNF on the divergent velocity;

➢ According to the ANOVA, the model's reliability is confirmed by the 94.98% contribution, which is significant and dominated by linear terms (51.72%), especially mmm (30.70%), and significant two-way interactions. Error contributes only 2.70%.

➢ The choice of $Go, CuO$, and ${Al}_{2}{O}_{3}$ with Casson fluid and ethylene glycol is founded on stability and thermal and rheological properties. Such a choice is made to determine the heat transfer and fluid stability for engineering and industrial applications. Further, the combination of three nanofluids in a base fluid is better than a single nanofluid.

Limitations

➢ The mesh density and quality have an impact on the accuracy of the finite element method. Slow convergence, errors, and instabilities can result from a bad mesh.

➢ Research has indicated that the Yamada–Ota model may overestimate in some situations because it frequently predicts higher temperatures than the Xue model.

➢ When dealing with three or four independent variables, the Box-Behnken method works well; however, as the number of input parameters increases, it becomes computationally costly and ineffective.

➢ When turbulence occurs in high-velocity or highly inertial flows, the Darcy-Forchheimer model may not hold because it assumes a particular type of resistance in porous media.

Future directions: The current model is based on Casson fluid in a Vertical jet using a ternary hybrid nanofluid with BBD. It is mentioned that this model can be solved in different non-Newtonian fluids like Jeffery fluid, Maxwell fluid, hyperbolic tangent fluid, Second-grade fluid, micropolar fluid, and Williamson fluid. In the future, the correlations of tetra hybrid nanofluid will be used for the measurement of thermal energy. Also, non-Fourier's law will be implemented. Such models can be used in biomedical engineering, power generation, automotive, and thermal management process.

Author contributions

Umar Nazir: Conceptualization, Methodology, Software, Formal analysis, Writing—original draft preparation, Project administration; Abdelaziz Nasr: Conceptualization, Methodology, Software, Validation, Formal analysis, Investigation, Resources, Writing—review and editing, Supervision, Funding acquisition. All authors have read and agreed to the published version of the manuscript.

Use of Generative-AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Funding

This research work was funded by Umm Al-Qura University, Saudi Arabia under grant number 25UQU4340531GSSR01.

Acknowledgments

The authors extend their appreciation to Umm Al-Qura University, Saudi Arabia for funding this research work through grant number 25UQU4340531GSSR01.

Conflicts of Interest

The authors declare no conflicts of interest.