Heat Transport and Thermal Efficiency in Magnetohydrodynamics Ternary Hybrid Nanofluid Flow past a Vertical Deformable Surface with Viscous Dissipation and Joule Heating Effects

1 Introduction

In the recent decades, nanotechnology has drastically revolutionized the study of heat transport and computational fluid analysis. The traditional heat transport fluids like water, ethylene glycol, and engine oil have low thermal conductivity, which restricts their usage in contemporary heat exchange systems [1,2]. To bridge this research gap, nanofluids (NFs) with nanoscale particles of solid matter were proposed in order to increase the heat transfer characteristics via dispersing nanoparticles such as graphene, titanium dioxide and zinc oxide into a base liquid. Based on this advancement, hybrid nanofluids (HNFs) were established via mixing two various kinds of nanoparticles (NPs) to improve thermal conductivity, stability, and heat transport performance [3]. Over the past few years, scientists have taken this idea further to a more complex point, the so-called ternary hybrid nanofluids (THNFs) that are substances composed of three various kinds of NPs suspended in a supporting liquid [4]. The novel kind of engineered fluids has improved thermophysical properties than single and binary HNFs due to the synergistic influence of NPs. The co-existence of metals, metal oxides, and carbon-based nanomaterials assists the thermal conductivity, balancing of viscosity, and enhancement of convective thermal transport performance [5,6]. THNFs have therefore become a growing area of interest in thermal energy systems and industrial engineering applications. THNFs are also important as a result of their outstanding capacity to improve thermal transport mechanisms, and the fluid remains stable and energy is minimized. The importance of THNFs is due to their outstanding fitting in to improve the processes of thermal transport without disrupting fluid stability and decreasing energy losses [7]. Ref. [8] recently studied a theoretical examination on a single-phase boundary-layer model for a THNFs consisting of silicon dioxide, titanium dioxide, copper over an extending sheet. In their study, the authors developed new formulations for the effective thermophysical properties via extending the traditional Hamilton–Crosser and Brinkman models to accommodate three various kinds of NPs. Ref. [9] investigated the heat generation and Magnetohydrodynamic (MHD) flow framework for Cu, Fe2O3, SiO2 via including a novel machine-learning approach. Their results showed that increasing magnetic field intensity suppresses the velocity profile whereas improving the fluid temperature profile, mainly as a result of a Lorentz drag force. The THNFs showed superior thermal absorption with a 22% improvement in the local Nusselt number.

Magnetohydrodynamic (MHD) flow has received a remarkable interest in the recent year due to its application engineering, and industrial development. Magnetohydrodynamic flow describes the movement of an electrically conducting fluid including liquid metals, plasma, electrolyte solutions, and saltwater, when subjected to the influence of a magnetic field [10]. The relationship of fluid transport and electromagnetic fields leads to distinctive characteristics that differentiate Magnetohydrodynamic flows compared to normal hydrodynamic systems. The electric currents caused via the movement of charged particles interact with the imposed magnetic field in these flows to create Lorentz forces [11]. Such forces change the momentum transport characteristics, decrease turbulence and impact heat and mass transport processes in the hydrodynamic medium. The basic physics of magnetohydrodynamic flow arises from the connection between the Navier–Stokes and Maxwell’s equations [12,13]. When a magnetic field is implemented to a conductive flow field, it boosts electric currents that create Lorentz forces opposing the fluid’s motion. This occurrence allows effective control of flow operation in the absence of direct mechanical contact, making magnetohydrodynamic fundamentals greatly valuable for numerous modern technologies [14]. Significant applications such as propulsion (MHD thrusters), astrophysical magnetic field generation, and MRI machines. In practical engineering setups, the ability to manipulate fluid motion via electromagnetic fields allows greater control of thermal and velocity transport [15]. For example, heat exchange systems, an applied magnetic field can control the thickness of boundary layer, impact temperature gradients, and alter thermal conductivity, thereby increasing system productivity. Ref. [16] investigated the impacts of internal forced-convective thermal transport in Magnetohydrodynamic NFs, with emphasis on behavior improvement mechanisms for heat exchangers and compact devices. Ref. [17] examined the combined influence of Hall current and ion-slip on radiative Magnetohydrodynamic rotating flows. Their outcomes showed that the rate of heat transport upsurges with the strength of the magnetic field but declines as values of the ion-slip rises.

Thermal radiation signifies a basic mode of thermal transport in which energy is released and transmitted as electromagnetic waves. In contrast to conduction and convection, which depend on a material medium, radiation occurs without of matter, such as in a vacuum [18]. All objects with temperatures above absolute zero release radiant energy exchange, with the strength and spectral distribution of this emission regulated via the material’s temperature and surface characteristics [19]. This performance is governed via the Stefan–Boltzmann principle, which explains that the radiation energy flux emitted via a surface is proportional to the fourth power of its absolute temperature. At raised temperatures, radiation becomes the predominant mode of thermal transport, playing an important function in determining the thermal productivity and behavior of engineering systems [20]. In several manufacturing processes and sustainable applications, thermal radiation energy flux plays a significant character in evaluating the energy balance. It is important in processes such as generating electricity, combustion, radiators of heat in cars, and medical treatments [21]. In aerospace and security technologies, for instance, radiation helps as the key mechanism for heat dissipation in space propulsion vehicle and re-entry vehicles operating within vacuum environments. Similarly, in sustainable energy production, operative control of radiative thermal transport is vital for preventing overheating and improving system optimization. When radiation interacts with fluid motion, it introduces an additional mechanism that significantly influences the overall heat transfer process. In flow of boundary layer, for instance, radiative heat transfer rate modifies the fluid temperature within the fluid, leading to changes in the thickness of thermal boundary layer and the overall rate of thermal transport [22]. Ref. [23] examined the 3D Magnetohydrodynamic (MHD) flow of nanofluids including unsteady, mixed convection effects, and radiation energy flux across an infinitely large rough disc. Their outcomes showed that thermal radiation significantly increase the fluid temperature. Ref. [24] study investigated Magnetohydrodynamic flow of a viscoelastic nanofluid near a stagnation point, explicitly including radiation energy flux and radiation energy flux at the surface. Outcomes indicate that thermal slip and radiation act together to increase the fluid temperature. Neuro-computational dynamics of an electromagnetically reactive copper-titania-water mixture within a quadratically accelerated Riga channel with graduated thermo-solutal conditions are study by [25]. Electromagnetic influences on fractional hemodynamics of trihybrid nanolayered nanoparticles in a diverging ciliated microtube: a machine learning-assisted analysis for biomedical applications are study by [26]. Buoyancy-driven magnetohydrodynamic nanofluid flow and heat transfer in a porous cavity with an exothermic reaction governed by Arrhenius kinetics are study by [27]. Fundamentals of nanofluids: Evolution, applications and new theory, international journal of biomathematics and systems biology are study by [28]. Magnetohydrodynamic boundary layer flow of nanofluid with variable chemical reaction in a radiative vertical plate are study by [29]. Latent heat-assisted thermal management of a metal hydride hydrogen storage reactor using novel shape-stabilized composite phase change material slurry is study by [30]. Photothermal-responsive form-stable phase change materials for enhanced solar absorption and evaporation in advanced desalination systems are study by [31].

This study presents a novel investigation of heat transport enhancement in a magnetohydrodynamic (MHD) stagnation-point flow of a ternary hybrid nanofluid composed of titanium dioxide (TiO2), cadmium telluride (CdTe), and molybdenum disulfide (MoS2) nanoparticles over a Darcy–Brinkman porous medium under thermal radiation. Unlike previous studies, which primarily focus on single or binary nanofluids, this work systematically explores the combined effects of three distinct nanoparticles, electromagnetic forces, viscous dissipation, and radiative heat transfer on both velocity and temperature distributions. The study further quantifies the influence of nanoparticle composition, magnetic field strength, and surface stretching/shrinking rates on the flow and thermal boundary layer, providing new insights into optimizing thermal performance in porous media. By integrating MHD effects with mixed convection in a ternary hybrid nanofluid context, this research offers a significant advancement in enhancing thermal transport efficiency, which has direct implications for advanced energy systems and engineering applications. Fig. 1 shows the schematic illustration of a (MHD) flow with NFs implementation.

Figure 1: Schematic illustration of a (MHD) flow with NFs implementation.

2 Mathematical Formulation of the Flow Problem

The study considers an incompressible flow, two-dimensional, radiative mixed convection flow of a viscous, through a vertical stretching/shrinking surface embedded in a Darcy–Brinkman porous medium. The flow and thermal transport are influenced via buoyancy-driven flow, radiative heat transfer, dissipation of viscous energy, and a heat source or sink within the stagnation-point flow. A Cartesian coordinate system ( x ˘ a , y ˘ a ) is used in such a way that the x ˘ a -axis lies along the surface (upward direction) while the y ˘ a -axis is taken normal to the surface (see Fig. 2a,b). The velocity components in the x ˘ a - and y ˘ a -directions are u ˘ a and v ˘ a , respectively. The wall (surface) velocity and temperature are denoted by U ˘ w a and T ˘ w a x ˘ a , whereas, the ambient (free-stream) velocity and temperature are denoted by the respective symbols U ˘ ∞ and T ˘ ∞ .

Figure 2: Demonstrates the geometry of the flow configuration together with the Cartesian coordinate system for two cases: (a) stretching sheet and (b) shrinking sheet.

The governing equations, after applying the aforementioned assumptions, are expressed as follows [13]:

The associated boundary conditions are:

In the above equations, u ⌣ a and v ⌣ a denote the velocity components of the (THNF) in the x ⌣ a - and y ⌣ a -directions, respectively, while x ⌣ a and y ⌣ a denote the corresponding Cartesian coordinates along and normal to the stretching surface. The parameter μ t h n f signifies the dynamic viscosity, and ρ t h n f is the effective density of the THNF. The term t h n f corresponds to the effective thermal conductivity, and ρ c p t h n f represents the effective heat capacity of the THNF. The symbol K a designates the permeability of the porous medium, whereas g denotes the gravitational acceleration. The quantity ( ρ β ) t h n f represents the product of the effective density and thermal expansion coefficient of the THNF. The term σ t h n f indicates the electrical conductivity of the thermally hybrid nanofluid, and σ f refers to the strength of the externally applied magnetic field. The parameter q r a d corresponds to the radiative heat flux in the transverse direction, while Q a represents the internal heat generation (or absorption) coefficient. The wall temperature distribution along the stretching surface is expressed by T w a x a , and the ambient temperature of the fluid far from the surface is denoted by T ⌣ ∞ . Similarly, U ⌣ w a and U ⌣ ∞ denote the wall stretching velocity and the free-stream velocity, respectively.

2.1 Radiative Heat Flux and Energy Equation Simplification

The radiative heat flux in the energy equation can be expressed in accordance with the Rosseland diffusion approximation as [23]: qrad=−4σ*3k*∂T⌣a4∂y⌣a(5) where σ * represents the Stefan–Boltzmann constant and k * denotes the mean absorption coefficient of the fluid medium. To simplify this nonlinear term, the temperature function T ⌣ a 4 is expanded about the free-stream temperature T ∞ using a Taylor series and neglecting higher-order terms.

Substituting Eq. (6) into Eq. (5) provides a linearized form of the radiative heat flux, which is then incorporated into the energy Eq. (3). Consequently, the modified energy equation can be written as

2.2 Thermophysical Properties of Based-Fluid and Nanoparticles

Table 1 presents the thermophysical attributes of TiO2 (titanium dioxide), CdTe (cadmium telluride), and MoS2 (molybdenum disulfide) and water-based fluid.

The properties of thermophysical properties of base fluid and nanoparticles are:

2.3 Nondimensionalization of Governing Equations

To simplify the governing equations and reduce the number of independent parameters, the following dimensionless variables are introduced [8]:

Substituting these transformations into the governing equations yields their nondimensional form, we obtain:

The associated boundary conditions are:

The dimensionless parameters and physical quantities are displayed in Table 2.

Here, l b is the characteristic length scale, U ⌣ ∞ denotes the free-stream velocity, and Re = U ⌣ ∞ l b / v f represents the Reynolds number of the thermally hybrid nanofluid. The term Δ T ⌣ = T ⌣ w a − T ⌣ ∞ defines the characteristic temperature difference between the wall and ambient fluid, while Q 0 and K 0 are the reference heat generation rate and characteristic permeability of the porous medium, respectively. Moreover, the condition γ d > 0 corresponds to an assisting flow (heated plate), γ d < 0 indicates an opposing flow (cooled plate), and γ d = 0 represents pure forced convection.

2.4 Similarity Transformations

In order to get the similarity solutions of the governing Eqs. (1), (2) and (6), it is necessary to introduce a proper transformation that results into a system of ordinary (similarity) equations. In this regard, the wall temperature, the heat source/sink and the permeability functions are determined as:

The stream function ψ is introduced in terms of the similarity variable ξ as: ψ=x⌣b12pξ,T⌣b=T⌣wbx⌣b=1x⌣bχξ,(14) in which the similarity variable is defined by:

The similarity transformations have been introduced in such a way that the independent variables are minimized by combining the spatial coordinates into one variable. As a result, the profiles of temperature and velocity generated will be geometrically the same at varying streamwise positions, despite the fact that the thickness of the boundary layer will be increasing with distance away the leading edge. This method is called the similarity transformation method and greatly simplifies the analysis since the same nonlinear partial differential equations are reduced to a simpler form of ordinary differential equations, which can be easily dealt with both analytically and numerically.

Subject to

2.5 Engineering Quantities of Interest

Two important quantities of engineering are the skin friction coefficient and the local Nusselt number which describe the surface drag and the rate of heat transfer respectively. Such parameters are given as: Cf=τwaρfU⌣∞2,Nu⌣xa=x⌣aqwakfΔT⌣(19) where τ w a denotes the wall shear stress and q w a represents the wall heat flux.

For nondimensional analysis, the shear stress and heat flux are expressed in terms of their dimensionless forms as:

By substituting the similarity variables and nondimensional quantities from Eqs. (8) and (14) into above equations, one obtains the dimensionless forms of the skin friction coefficient and Nusselt number as:

3 Numerical Solution

3.1 Application of Shifted Legendre Collocation Method (SLCM)

In this section, the step-by-step technique for applying the SLCM to the solution of the highly non-linear system of the differential equations resulting in the boundary layer analysis is discussed. The methodology involves formulating the problem into a finite interval, approximating using truncated Shifted Legendre series the unknown functions and reducing the system of ODEs and boundary conditions into a non-linear algebraic equation system which is solvable. The governing equations are:

Subject to

Step 1: Map the semi-infinite domain

We map the semi-infinite physical domain via a standard algebraic map ξ=Lt1−t,t=ξL+ξ(27) where L > 0 is a scaling parameter, under this map:

Let

Then

• First derivative:

• Second derivative:

Because a ′ = d a d t = − 2 1 − t L .

Third derivative (derived by applying d d ξ = a D t to the second derivative):

Step 2: Introduce shifted Legendre polynomials

Define the shifted Legendre polynomials P n * t on 0 , 1 by

Choose an approximation order N. Approximate p ξ and χ ξ by pξt≈pNt=∑k=0NakPk*t,χξt≈χNt=∑k=0NbkPk*tg(34) with unknown coefficients a k k = 0 N , b k k = 0 N . We will need derivatives D t p N , D t 2 p N , D t 3 p N

Similarly,

Step 3: Convert the ODEs to equations in t and substitute above expressions, we obtain

Step 4. Construct the collocation system

Choose collocation points t j in [0, 1]. A common and stable choice for polynomial collocation is the shifted Legendre–Gauss–Lobatto nodes in [0, 1] tj=1+xj2,j=0,…,N(39) where x j are the Legendre–Gauss–Lobatto nodes on [−1, 1]. residual functions R ( 1 ) ( t ) and R ( 2 ) ( t ) by substituting the expansions into the t-form of equations such as:

Enforce the collocation equations

R ( 1 ) ( t j ) = 0 , R 2 t j = 0 , j = 1 … , 2 . This yields 2 N − 1 nonlinear algebraic equations for the 2 N + 1 unknown coefficients a k k = 0 N , b k k = 0 N . The remaining equations come from the boundary conditions.

3.2 Residual Error and Convergence Analysis

Here, Table 3 illustrates spectral convergence, the exponential reduction in error with increasing degree of the polynomial. The following values are exemplary and demonstrate the approximate size and direction a well converged SLCM solution. Also, Table 4 displays the calculated values of p ″ 0 and χ ′ 0 as the polynomial degree N at N = 20 . The comparison of the considered code for the limiting cases are presented in Table 5. The numerical results show strong agreement with previously published data, confirming the reliability and precision of the computational model.

4 Result and Discussion

This study analyzes the steady, laminar, two-dimensional flow and heat transfer of ternary hybrid nanofluids, consisting of TiO2, CdTe, and MoS2 nanoparticles dispersed in water, over a vertical stretching or shrinking surface embedded in a Darcy–Brinkman porous medium. The study considers practical and engineering-relevant parameter ranges. The study considers practical and engineering-relevant parameter ranges, including permeability, δ d = 0.5, 0.7, 0.9, mixed convection parameters, γ d = 0.1, 0.3, 0.5, velocity ratio parameter, Σ d = 0.4, 0.6, 0.7, Eckert number, E c d = 0.2, 0.4, 0.6, heat source/sink parameter, Γ d = 0.2, 0.3, 0.5, and thermal radiation parameter, N r = 0.5, 0.7, 0.9. The present study analyzed the magnetohydrodynamic (MHD) flow of convection heat transfer in TiO2-CdTe-MoS2 nanoparticles with viscous dissipation and thermal radiation. Results obtained are discussed in terms of velocity and temperature profiles and Nusselt number and skin friction coefficient with different controlling parameters. From Table 6, increasing value of δ b upsurges the skin friction ( p ″ ( 0 ) ) but slightly declines the − χ ′ 0 . Increasing Pr leads to a thinner thermal boundary layer, drastically growing the − χ ′ 0 . It has no influence on p ″ ( 0 ) in this model structure. N r Improve the thermal transport via radiation, effectively thickening the thermal boundary layer, which declines the − χ ′ 0 . Increasing E c d added more heat into the fluid via friction, growing the internal energy and thus decreasing the − χ ′ 0 . δ d rises the velocity gradient at the wall, leading to a direct increase in p ″ 0 . It also slightly declines the − χ ′ 0 . A positive G r b adds energy, decreasing the temperature gradient at the wall and reducing the − χ ′ 0 .

Effect of permeability δ d and mixed convection parameters γ d on p ′ ξ : Fig. 3a,b, respectively, demonstrate the variations in the p ′ ξ due to the impact of the permeability δ d and mixed convection parameters γ d . It is noted that as the permeability parameter upsurges, the p ′ ξ of the region of the boundary layer decreases. This can be explained physically via the porous character of the medium which applies a resistance force to the movement of the fluid particles. The greater the permeability parameter the less permeable the medium is and hence the penetration of the fluid gets limited and consequently the intensity of the flow is suppressed. This drag effect is called the Darcy drag and occurs as a result of the frictional interaction between the solid matrix of the porous media and the fluid. It is shown in Fig. 3b that the γ d improves the velocity profile. The increase in the γ d amplifies the thermal buoyancy effects, which help the fluid movement along the surface. The temperature between the wall and the surrounding fluid causes buoyant forces that support the upward movement of the fluid and quickens the total velocity. Physically, this intensification is related with the transformation of the thermal energy into the kinetic energy inside the flow domain. This is a procedure that is significant in improving velocity and thermal distribution in systems where the thermal gradient is high. For both the Fig. 3a,b, the TiO2-CdTe-MoS2/water ternary hybrid nanofluid has a greater velocity increase than the conventional CdTe-TiO2/water hybrid nanofluid in both parameters δ d and γ d . This is because of the synergistic effect of the three different nanoparticles suspended in the base fluid, which together promote the thermal conductivity and alter the viscous characteristics of the suspension.

Figure 3: (a) Behavior of velocity distribution under the effect of δ d . (b) Variation of fluid flow with γ d .

Effect of magnetic parameter ( M ) and Grashof number ( G r b ) , and Σ d velocity ratio parameter on p ′ ξ : Fig. 4a,c displays the change in the p ′ ξ with the impact of the various physical parameters. It is noted that (M) has a strong retarding influence on the flow field. The resistive drag that suppresses the velocity profile by introducing a Lorentz force, who works in the opposite direction of fluid motion. Physically, this is due to the fact that when an electrically conducting fluid is moved in the presence of a magnetic field, currents develop which interact with the magnetic field to produce a new opposing electromagnetic force. The higher the M the greater this opposition, the kinetic energy is lost through heating and thus the p ′ ξ is slowed down. The Grashof number ( G r b ) is known to have an enhancing effect on the velocity profile, as revealed in Fig. 4b. An increase in this parameter causes the buoyancy forces due to temperature variations to overpower the viscous resistance and result in an increased upward or outward flow of the fluid. This upsurges the impact of natural convection and increases the rate of flow of the fluid. The improvement of the p ′ ξ with a growing G r b shows the significance of thermal buoyancy in convective transport in the CdTe-TiO2-MoS2/water ternary hybrid nanofluid systems. The heat convected to the nanoparticles enables a more intense fluid agitation, which consequently improves more intense p ′ ξ . An increase in Σ d means a higher rate of stretch of the surface compared to the ambient flow, which is inclined to carry a larger number of fluid particles along the surface, thereby making the velocity near the boundary layer region higher (Fig. 4c). The the CdTe-TiO2-MoS2/water ternary hybrid nanofluid exhibits a stronger improvement of the velocity magnitude than the the CdTe-TiO2-MoS2/water hybrid nanofluid in all the parameters investigated.

Figure 4: (a) Behavior of velocity distribution under the effect of M . (b) Variation of fluid flow velocity with G r b . (c) Impact of Σ d  on velocity field.

Fig. 5a,b shows that changes in the M strength and E c d can significantly impact the temperature field in both the CdTe-TiO2/water hybrid nanofluid and CdTe-TiO2-MoS2/water ternary hybrid nanofluid systems. The χ ξ is evidently growing with the increasing M strength and E c d , and this magnification is more prominent in the CdTe-TiO2-MoS2/water ternary hybrid nanofluid as compared to the traditional CdTe-TiO2/water hybrid nanofluid. This performance can be attributed to the intensified energy conversion processes and augmented internal frictional heating that take place within the fluid medium due to these parameters. The imposition of a magnetic field induces a Lorentz force that acts normal to the fluid flow. This electromagnetic opposition reduces the velocity gradients, thereby converting a portion of the kinetic energy into thermal energy. As a result, the χ ξ displays an upward trend as the M upsurges. In electrically conducting nanofluids, including those comprising metallic nanoparticles, the M induces extra Joule heating influence that further promote the χ ξ . The ternary hybrid nanofluid, including three different nanoparticles, displays greater thermal conductivity as a result of the synergistic relations among the nanoparticle constituents. The outcomes noticeably show that both the M strength and E c d number exert a positive effect on the χ ξ for CdTe-TiO2-MoS2/water ternary hybrid nanofluid flow.

Figure 5: (a) Behavior of χ ξ  under the effect of M .  (b): Variation of χ ξ  with E c d .

The variation in the χ ξ of both CdTe-TiO2/water hybrid nanofluid and CdTe-TiO2-MoS2/water ternary hybrid nanofluid under the combined effect of N r and Γ d is demonstrated in Fig. 6a,b. The outcomes noticeably display that an increase in the N r and Γ d leads to a substantial increase in the χ ξ within the boundary layer. This improvement is more pronounced in the ternary hybrid nanofluid compared to the hybrid nanofluid, showing the higher thermal transport potential of the ternary combination due to its improved thermal conductivity and synergistic nanoparticle interactions. Physically, the increase in the N r improve radiative heat flux within the fluid domain, which served as an extra mode of energy transfer apart from conduction and convection. Likewise, the increase in the Γ d parameter corresponds to a rise in volumetric energy discharge within the fluid system, thereby supplying to the internal thermal reservoir. Also, the effective surface area for energy absorption, thereby intensifying the impact of both N r and Γ d . Therefore, incorporating ternary CdTe-TiO2-MoS2/water ternary hybrid nanofluid in thermally intensive engineering applications can considerably improve productivity, reliability, and energy consumption.

Figure 6: (a) Behavior of χ ξ  under the effect of N r . (b): Variation of χ ξ  with Γ d .

5 Conclusion

In this study, heat transport enhancement in magnetohydrodynamic flow of nanofluids containing TiO2 (titanium dioxide), CdTe (cadmium telluride), and MoS2 (molybdenum disulfide) nanoparticles over a vertical stretching/shrinking surface with thermal radiation has been analyzed. The study revealed that the interaction between electromagnetic forces, dissipation of viscous energy, and radiation intensity strongly governs both flow and heat transfer boundary layers. These results emphasize the potential of THNFs for advanced regulation of thermal energy, energy conversion, and heat removal systems, providing a theoretical basis for optimizing MHD thermal energy transfer in mechanical systems. Based on the numerical findings, the significant final remarks are outlined as follows:

• The M parameter exerts a retarding effect on the fluid motion, leading to a noticeable reduction in the fluid flow due to the Lorentz force opposing the flow.
• Greater values of the Grb increase the fluid flow, confirming the dominance of buoyancy-driven forces that encourage stronger convective currents.
• The Σd parameter positively impacts the velocity distribution, increasing the speed near the boundary region.
• The M, Ecd, Nr, Γd parameters all contribute to advancement in the fluid temperature, increasing thermal energy within the boundary layer.
• Improved Ecd values reveal the conversion of kinetic energy into internal energy through viscous dissipation, resulting in greater thermal distribution.
• The rise in Nr improves thermal transport and increases the thermal boundary layer thickness.
• The ternary hybrid nanofluid consistently displays greater improvement in both fluid flow compared with the hybrid nanofluid, owing to enriched thermal conductivity and synergistic particle interaction.
• The combined effect of electromagnetic forces, buoyancy influence, and internal heat sources plays an important function in determining the speed and thermal distribution performance.
• Practical usage of this research includes thermal management in magnetohydrodynamic (MHD) generators, microelectronic cooling systems, heat exchangers solar collectors, and nuclear reactor.

Limitations and Future Research

Despite the valuable insights provided by this study, some limitations should be acknowledged. First, the analysis is based on a steady, laminar, and two-dimensional flow assumption, which may not fully capture transient effects, turbulence, or three-dimensional behaviors encountered in practical engineering applications. Second, the thermophysical properties of the ternary hybrid nanofluid are assumed to be constant, whereas in reality they may vary with temperature, nanoparticle concentration, agglomeration, and shear rate. In addition, nanoparticle shape, size distribution, and possible slip mechanisms such as Brownian motion and thermophoresis are neglected. The porous medium is modeled as homogeneous and isotropic using the Darcy–Brinkman formulation, which may oversimplify complex porous structures found in real systems. Moreover, the study is purely numerical and lacks experimental or industrial validation to corroborate the predicted trends.

Future research may extend the present work by considering unsteady, three-dimensional, and turbulent flow configurations, as well as temperature-dependent and concentration-dependent fluid properties. Incorporating nanoparticle interaction effects, variable magnetic fields, non-uniform heat sources, and more realistic porous media models would further enhance the physical relevance of the model. Experimental studies or high-fidelity numerical simulations could be conducted to validate the theoretical predictions. Additionally, coupling the present framework with optimization techniques or machine-learning approaches could provide deeper insight into optimal parameter ranges for advanced thermal management and energy systems.