Joule Heating and Viscosity-Ratio Effects on Dissipative Ternary Nanofluid Flow over a Permeable Surface
Abstract:
This study examines the effects of viscous dissipation, Joule heating, and coupled heat transfer on dissipative ternary nanofluid flow over a permeable surface. The ternary nanofluid is composed of $A{{l}_{2}}{{O}_{3}}$, $Si{{O}_{2}}$, and $Ti{{O}_{2}}$ nanoparticles dispersed in water as the base fluid. By introducing suitable similarity transformations, the governing partial differential equations are reduced to a coupled system of ordinary differential equations. The thermal field is analyzed for both prescribed surface temperature (PST) and prescribed heat flux (PHF) conditions, while a temperature-dependent heat source/sink term is incorporated to maintain energy balance within the fluid domain. The resulting energy equation is treated analytically with the aid of Kummer’s function and Laguerre polynomial techniques. The effects of the main controlling parameters, including the inverse Darcy number, magnetic parameter, viscosity-ratio parameter, and radiation parameter, are discussed with the support of graphical results. It is found that an increase in the magnetic parameter reduces the velocity by about 12% and raises the temperature by nearly 18%. These findings provide useful guidance for the design and thermal optimization of engineering systems involving complex nanofluids in porous media, including polymer extrusion and automotive cooling applications.
1. Introduction
Magnetic fields influence the movement of fluids in numerous industrial applications, where they serve as a novel, non-invasive tool to direct and control the flow of liquids. A typical example is the water treatment sector, where magnetic flow meters are chosen for their high accuracy in gauging the flow rate of water and wastewater. This is a pivotal technology for the management and distribution of water resources. Food and beverage industries rely on these meters to help maintain standard product quality by precisely measuring the flow of liquid ingredients. Likewise, the pharmaceutical industry turns to the accuracy of magnetic flow meters to carry out hygienic and exact measurements of sensitive fluids that are not in direct contact. These fluids allow for the adjustment of resistance and damping in different scenarios, thus improving the performance and safety of vehicles and machines. The scope of magnetic applications in fluid dynamics underlines their indispensable function in the evolution of industrial operations and product quality. Due to many applications, many authors examined magnetic flow, such as Dolui et al. [1], who examined the importance of ternary nanoparticle flow, particularly on the interactions at the Hiemenz stagnation point and the influence of a magnetic field on the flow behaviour as it moves over a surface. Turkyilmazoglu [2] explored the dynamics of porous media with heat transfer. Jan et al. [3] and Mishra et al. [4] studied the magnetohydrodynamics (MHD) 2-D ternary nanofluid flow based on polymers with different thermal conductivity and partial slip. The influences of the MHD and Cattaneo-Christov models on a Maxwell fluid stagnation-point movement are explored by Chen et al. [5].
The viscosity ratio plays a major role in various industrial and real-world applications. It is instrumental in the examination of fluid flow through porous media, which is a common scenario in sectors such as chemical engineering, environmental analysis, and biomedical engineering. For instance, in chemical engineering, the viscosity ratio is used to optimize the design of filtration mechanisms, ensuring efficient separation processes. In environmental implications, many researchers explored the viscosity ratio, such as Pavlov [6], who studied the two-dimensional MHD fluid movement. Zhang et al. [7] examined the electro-MHD characteristic in non-Newtonian fluids, particularly of the third grade, involving interactions between the fluid’s inherent characteristics and the external electric and magnetic fields. The Gauss-Seidel iteration method and quasi-linearization are used by Yao et al. [8] and Murali et al. [9] to study the 3-D Darcy-Brinkman-Forchheimer flow terms. The Brinkman impact on fluid flow was studied by Tso and Mahulikar [10] and Bèg et al. [11] by using numerical findings produced by a robust double-shot Runge-Kutta-Merson scheme and a finite-difference approach. Moreover, the viscosity ratio’s relevance extends to the biomedical sphere, particularly in drug delivery systems, where understanding fluid flow at the microscale can enhance the precision of targeted therapies. The adaptability of the viscosity ratio model to Newtonian and non-Newtonian fluids makes it a versatile method for researchers and engineers, allowing for precise predictions and improvements in fluid flow systems over a broad spectrum of applications.
Joule heating is used to keep precise temperatures in reactors and distillation columns, which is necessary for the uniformity and quality of the products. The plastics industry benefits from this effect in operations such as injection moulding and extrusion to adjust the viscosity of the material so that it flows easily and is shaped accurately. In the food and beverage industry, Joule heating is used for distillation, drying, cooking, etc., where consistent and uniform heat distribution is the key. Because of numerous applications, a number of researchers have studied Joule heating, such as Nayak et al. [12], who explored the MHD nano liquid flow via a permeable linear stretching surface with convective heating and a porous matrix. Leng et al. [13] and Nyundo et al. [14] explored the consequences of suction/injection and MHD Casson nanofluid flow for a nonlinearly expanding surface. Nazeer et al. [15] explored the Newtonian liquids Walters-B liquid transport mass and energy when they enter the flow via a penetrable surface that is semi-infinite. Khan et al. [16] researched the MHD fluid flow over a linear stretching surface with convective heating in the presence of a porous medium. Furthermore, in microfluidics, the heat generated by the Joule effect is used to open microvalves by the thermal expansion of the parts, thus showing its dual capability in both large-scale operations and small, highly precise applications.
Porous media is playing a main role in the direction and control of fluid flow in different industrial and real-world applications. The materials, which are permeable by nature, are allowing the controlled movement and spreading of fluids; porous media are very important in extraction processes. For example, it is essential to understand the multiphase flow within porous rocks for the efficient recovery of oil. Similarly, in environmental engineering, the study of fluid flow through soil and rock is fundamental for groundwater management and contamination mitigation. Sachhin et al. [17] studied the slanted magnetic field that affects fluid flow across a permeable, shrinking surface while allowing for heat transmission. Patnaik et al. [18] and Bég et al. [19] explored the convective boundary layer movement of a nanofluid via an extending sheet in a porous medium under unstable MHD conditions. A nanoparticle-enhanced porous medium embedded on a vertical flat surface with a steady mixed convection flow is investigated by Waini [20]. Joshi et al. [21] and Wang et al. [22] have examined the effects of heat radiation along with different viscosities on the flow of nanofluid over a rotating disk with permeable media. The development of microfluidic model porous media has made it possible to see and understand the complex fluid dynamics at the pore level, which has resulted in better designs and processes in industries such as food and wood production. Furthermore, the combination of data-driven methods and physics-informed machine learning is changing the traditional modeling and prediction of these flows, thus giving new insights and optimization strategies for these vital systems.
Thermal radiation fundamentally affects the communication of heat and thus energy transfer efficiency in fluid flows through industrial processes. From a thermal engineering perspective, figuring out how thermal radiation affects fluid flow is a must if one is to optimize industrial processes. For instance, the chemical industry uses radiation to keep reactors or distillation columns at the right temperature so the reaction is just efficient and the product quality is good. Likewise, in plastic industries, temperature control through radiation during injection molding and extrusion can help in obtaining the desired properties of materials. The food and beverage industry also uses radiation to provide the required heat in processes like distillation, drying, and cooking. Due to many applications, many researchers examined the radiation, such as Sachhin et al. [23], who explored how a slanted MHD influences fluid movement over an expanding/shrinking sheet while allowing for heat transmission. Yu and Wang [24] and Li et al. [25] explored the flow and temperature between two extending, spinning disks within a water flow. Bejawada and Nandeppanavar [26], Hosseinzadeh et al. [27], and Megahed [28] studied the micropolar liquid flow across a vertical permeable plate. Also, in high-temperature environments such as nuclear reactors or electronic cooling systems, thermal radiation helps to control heat transfer; thus, the safety and the life span of the equipment are improved. Incorporating the concepts of thermal radiation in fluid flow systems results in improvements in energy conservation and process efficiency, thereby highlighting its importance in industrial applications as well as in everyday life. Viscosity dissipation is taken into consideration by Ragueb and Mansouri [29] to study the laminar heat transfer inside a cross-section duct that is exposed to a uniform wall temperature. Schranner et al. [30] and Lelea and Cioabla [31] study the problem of heat transmission and the flow of incompressible fluid movements across various boundaries. Bataller [32], Tso et al. [33], and Megahed [28] study the viscous dissipation effects. Dadheech et al. [34], Crane [35], and Hamid et al. [36] study the heat transfer in stretching surfaces. Xiang et al. [37], Shanmugam and Maganti [38], and Agravat et al. [39] study the heat source/sink performance. Ullah et al. [40] and Flity et al. [41] study the effect of thermal conductivity on fluid flow. Bognár [42], Lone et al. [43], and Barna et al. [44] develop similarity solutions and analytical approaches for non-Newtonian and boundary layer flows. Klazly and Bognár [45] and Bognár and Hriczó [46] investigate empirical modeling of nanofluid properties and ferrofluid flow under magnetic effects. Usafzai et al. [47] comprehensively studies nanofluid and hybrid nanofluid transport phenomena over stretching and permeable surfaces.
Mahabaleshwar et al. [48], [49], [50] analyze porous media flow, convection, and magnetohydrodynamic effects in complex fluids. Khan et al. [51], Hatami and Ganji [52], Jenifer and Saikrishnan [53], Sachhin et al. [54], and Hamid et al. [36] present advanced studies on entropy generation and thermal transport mechanisms.
The present model is motivated by the thermal management of high-heat-flux electronic devices mounted on porous substrates, such as microprocessors, power electronics modules, and compact heat exchangers. In these systems, electrically conducting ternary nanofluids are circulated through porous heat sinks to enhance cooling performance. The presence of an external magnetic field enables control of the fluid motion, making MHD effects essential for accurate prediction. Since the coolant flows through porous structures (e.g., metal foams), the Darcy–Brinkman formulation is required to account for both permeability resistance and viscous diffusion within the medium. The electrical conductivity of the nanofluid leads to Joule heating under magnetic fields, while high shear near stretching surfaces produces viscous dissipation. At elevated temperatures, thermal radiation further contributes to heat transfer, necessitating the simultaneous inclusion of all these mechanisms.
The two thermal boundary conditions considered also reflect practical operating modes. The prescribed surface temperature (PST) condition represents actively temperature-controlled systems where the wall temperature is maintained constant, whereas the prescribed heat flux (PHF) condition corresponds to devices operating at fixed power dissipation, where the surface temperature adjusts according to the heat removal rate. Thus, the present formulation provides a realistic theoretical framework for magnetically controlled cooling of porous electronic systems using advanced ternary nanofluids.
The motivation of the present study comes from the real-life applications of thermal management and materials processing systems where electrically conducting nanofluids are made to flow over porous and permeable surfaces while being subjected to external magnetic fields. For instance, one may think of cooling the electronic components, which are mounted on porous substrates, using magnetic fields, or the polymer extrusion over porous molds that are electrically heated, or the energy systems where the nanoparticle‐enhanced coolants are circulated through the magnetized porous heat exchangers. In such scenarios, the coexistence of strong magnetic fields and electrically conducting hybrid ternary nanofluids leads to the generation of Joule heating, while the viscous dissipation, caused by the high shear rates in the vicinity of stretching or shrinking surfaces, results in thermal changes that are of significant magnitude. Besides that, at high operating temperatures, radiation heat transfer cannot be disregarded, especially in porous media where the radiative heat flux can combine with the conduction and convection mechanisms.
The porous medium model is used to describe the scenario accurately by including the contribution of viscous and permeability effects, which are inseparable in a porous structure and cannot be explained with the classical Darcy formulation, which refers only to the permeability effect. The choice of both surface temperature (PST) and heat flux (PHF) boundary conditions is made to indicate that, in practice, there exist two types of thermal control strategies, and both are to some extent standard industrial procedures, depending on whether it is surface temperature or heat input that gets regulated. Ignoring any of these mechanisms may cause one not to have a complete picture of the underlying physics and, in turn, may lead to the wrong evaluation of velocity, temperature, skin friction, and heat transfer rates. Therefore, the present work offers a consistent and physically reliable theoretical basis that integrates the co-action of hybrid ternary nanofluids, MHD, porous media, Joule heating, viscous dissipation, and thermal radiation into a single problem.
2. Mathematical Analysis and Solution
The current mathematical setup relies on several assumptions about physical phenomena. Figure 1 the physical representation of the problem. The flow is assumed steady, 2D, laminar, and incompressible. The nanofluid base with ternary nanoparticles uniformly distributed was treated as a single-phase medium with modified thermophysical properties. The porous medium is homogeneous and isotropic, and the porous medium’s effect is introduced by the Darcy-Brinkman model that allows for both permeability resistance and viscous diffusion effects. Thus, the magnetic field is tilted. Hall and ion-slip effects are not taken into account. The electrical conductivity of the fluid is considered constant. The pressure gradient along the surface is neglected, and body forces other than the Lorentz force are not taken into consideration. In the energy equation, viscous dissipation and Joule heating are kept since they are significant in high-shear and electrically conducting flows. Thermal radiation is accounted for by using the Rosseland diffusion approximation, which is suitable for optically thick media where the radiative heat flux acts like a diffusive process. Higher-degree temperature terms in the expansion of the radiative heat flux are discarded for the sake of mathematical convenience while still keeping the leading radiative part.

Governing equations [1], [8], [12], [24]:
Boundary conditions are given as [2], [6], [17]:
PST and PHF case boundary conditions (B. Cs) are:
where, $D$ denotes a mass flux term, $S$ is the mass transpiration at the plate, $b$ and $l$ terms are constants.
Thermal radiative heat flux estimated by Rosseland’s approximation is written as [2], [6], [17]:
where, $ {{\sigma }^{*}}$ is the Stephen Boltzmann constant and $ {{k}^{*}}$ is the absorption parameter.
By using the Taylor series and solving $ {{T}^{4}}$ and ignoring the upper-order terms to get:
Hence, Eq. (3) becomes:
Similarity conversions are given as [2], [6], [17]:
Given the considerations above, the original system of PDEs in the description of the process is simplified to the standard boundary layer equations that allow similarity transformations to be carried out, thus obtaining the system of nonlinear ordinary differential equations (ODEs) in Eqs. (12)–(18):
PST case:
PHF case:
Modified B. Cs are obtained as:
PST case:
PHF case:
The governing nondimensional terms utilized in Eqs. (12) to (18) can be stated as:$ M\left( =\frac{{{\sigma }_{f}}{{B}_{0}}^{2}}{a{{\rho }_{f}}} \right)$ is the magnetic field, $ \Pr \left( =\frac{{{\mu }_{f}}{{\left( {{C}_{p}} \right)}_{f}}}{{{\kappa }_{f}}} \right)$ is Prandtl number, $ Nr\left( =\frac{16{{\sigma }^{*}}{{T}_{\infty }}^{3}}{3{{\kappa }_{f}}{{k}^{*}}} \right)$ is thermal radiation, $ Ni\left( =\frac{{{Q}_{0}}}{a{{\left( \rho {{C}_{p}} \right)}_{f}}} \right)$ represents heat source.
Now, ${{C}_{f}}$ and $ Nu$ is calculated by:
where,
Skin friction and Nusselt numbers as:
where, $ {{\operatorname{Re}}_{x}}=\frac{a\,{{x}^{2}}}{{{\nu }_{f}}}$ is the local Reynolds number.
By the B. Cs the solution obtained as:
where, $ \beta >0$ and:
Put Eqs. (24) and (25) in Eq. (12), then we get the following equation:
The roots of Eq. (26) are as:
where,
$ \begin{gathered} E_1=\Lambda \xi_2, E_2=\left(A_2 S \xi_2-\Lambda \xi_1\right), E_3=\left(A_2 S \xi_1-\xi_2 A_1 D a^{-1}-A_3 M \operatorname{Sin}^2(\tau) \xi_2\right), \\ E_4=\left(A_2 S+A_3 M \operatorname{Sin}^2(\tau) \xi_1+\xi_1 A_1 D a^{-1}\right) E_5=A_2 d+A_3 M \operatorname{Sin}^2(\tau)+A_1 D a^{-1} \\ \Delta_1=E_3^2-3 E_2 E_4+12 E_1 E_5, \Delta_2=2 E_3^2-9 E_2 E_3 E_4+27 E_1 E_4^2+E_2^2 E_5-72 E_1 E_3 E_5 \end{gathered} $
Reduced skin friction coefficient can be evaluated by:
To calculate the energy equation, let us take the new term t as:
The Eq. (13) calculated as:
where,
Modified B. Cs are:
By utilizing Eqs. (32) and (30) calculated as:
where, $ {{B}_{1}},$ $ {{B}_{2}},$ and $ {{B}_{3}}$ are defined as:
Now, the wall temperature gradient calculated as:
By Eqs. (32) and (35), then the PHF Eq. (14) has the same form as Eq. (33) and the parameters in Eq. (34) expect $ E{{c}_{3}}$ , which is to be replaced here by $ E{{c}_{4}}$ calculated as:
Modified B. Cs for PHF case written as:
Hence, the solution of the Eq. (14) is given by:
Further, the wall temperature is obtained as:
• Dynamic viscosity:
• Density:
• Thermal capacitance:
• Thermal conductivity:
• Electrical conductivity:
Table 1 represents the thermophysical properties of base fluid and nanoparticles. The dimensionless parameter ranges in Table 2 have been chosen in such a way that they correspond to MHD nanofluid flows in porous cooling systems and other similar situations in related applications where the literature has reported the values. The moderate magnetic parameter values limit the magnetic Reynolds number, thereby supporting the assumption of a negligible induced magnetic field. The inverse Darcy number is representative of porous materials such as metal foams and filtration media. The range of the Eckert number corresponds to the conditions of moderate velocity flows where the viscous dissipation effect is significant, but the boundary layer approximation is still valid. Likewise, the radiation parameter characterizes the condition of elevated temperatures at which radiative heat transfer takes place along with conduction and convection.
Properties | Water $\left(\boldsymbol{H_{2}O}\right)$ [53] | $\boldsymbol{Al_{2}O_{3}}$ [53], [54] | $\boldsymbol{SiO_{2}}$ [53] | $\boldsymbol{TiO_{2}}$ [54] |
|---|---|---|---|---|
$\rho$ | 997.1 | 3970 | 2270 | 4250 |
$C_{p}$ | 4179 | 765 | 730 | 686 |
$\kappa$ | 0.6013 | 40 | 1.4 | 8.85 |
$\sigma$ | 0.05 | 10$^{\text{-12}}$ | 3.5 $\times$ 10$^{\text{6}}$ | 10$^{\text{-10}}$ |
| Parameter | Physical Meaning | Engineering Relevance | Parameter Range |
|---|---|---|---|
| $\Lambda$ | Viscosity ratio | Polymer extrusion, microchannels | 0 $\le \Lambda \le$ 3 |
| $\phi$ | Nanoparticles volume fraction | Stable ternary nanofluids | 0 $\le \phi \le$ 0.02 |
| $M$ | Magnetic field parameter | MHD cooling, magnetic pumps | 0 $\le M \le$ 3 |
| $Da^{-1}$ | Inverse Darcy number | Porous heat exchangers, filtration | 0 $\le Da^{-1} \le$ 3 |
| $Nr$ | Thermal radiation parameter | High-temperature systems | 0 $\le Nr \le$ 3 |
| $N_i$ | Heat source/sink parameter | Reactive or electrically heated systems | -3 $\le N_i \le$ 3 |
| $d$ | Stretching/shrinking parameter | Sheet extrusion, coating | -3 $\le d \le$ 3 |
| $\xi_{1}$ | First-order velocity slip | Micro/nano flows | 0 $\le \xi_{1} \le$ 3 |
| $\xi_{2}$ | Second-order velocity slip | Higher-order rarefaction effects | 0 $\le \xi_{2} \le$ 3 |
| $S$ | Mass transpiration parameter | Boundary layer control | 0 $\le S \le$ 3 |
| $\Pr$ | Prandtl number | Water-based nanofluids | 6.21 |
| $Ec_{1}, Ec_{2}$ | Eckert numbers | High shear or high-speed flows | 0 $\le Ec_{1}, Ec_{2} \le$ 3 |
The volume fraction of nanoparticles is limited to very low concentrations, as is the case in stable ternary nanofluids; thus, a single-phase model and thermal equilibrium assumption can be adopted. As a result, the present model finds its use in the case of laminar boundary layer flows with low magnetic Reynolds numbers, homogeneous porous media, and well-dispersed nanoparticles. If one considers very high particle loading, strong magnetic induction effects, or turbulent regimes, then the use of advanced multiphase or fully coupled models would be necessary.
Table 3 demonstrates that the present analytical formulation reduces exactly to several well-established models available in the literature under appropriate limiting assumptions. The comparison confirms that the obtained momentum expressions and characteristic equations are consistent with previously reported results.
Related Works by Other Authors | Fluids | Value of Momentum Solution |
|---|---|---|
[35] | Newtonian | $\beta = 1$ |
[6] | Newtonian | $\beta = \sqrt{1+M}$ |
[36] | Newtonian | $f(\eta)=S+\frac{d}{\beta(1+b \beta)}\left(1-e^{-\beta \eta}\right)$ $(2+\lambda) b \beta^3+(2+2 \lambda-2 b S) \beta^2-2 S \beta-2 a=0$ |
[2] | Newtonian | $f(\eta)=S-\frac{1-e^{\beta \eta}}{\beta}$, |
[47] | Newtonian | $f(\eta)=S+\frac{a}{\left(\beta+b \beta^2\right)}\left(1-e^{-\beta \eta}\right)$, $(2+\varepsilon M) b \beta^3-(2+\varepsilon M-2 S b M) \beta^2-2 S M \lambda-2 a M=0$ |
Present work | Newtonian | $f(\eta)=D_1+D_2 e^{-\beta \eta}, D_1=S-D_2, D_2=\frac{d}{\xi_2 \beta^3-\xi_1 \beta^2-\beta}$ $\Lambda \beta^2-A_2 S \beta-A_3 M \operatorname{Sin}^2(\tau)-A_1 D a^{-1}=0$ |
• Excluding volume fraction, radiation, porous media and magnetic field term, viscous dissipation {Present results} $\to$ {results of Crane [35]}.
• Excluding permeable medium, ternary nanoparticles,heat dissipation, and including of nano fuid {Present results} $\to$ {results of Hamid et al. [36]}.
3. Results and Discussion
The current analysis examines the influence of viscous dissipation, Joule heating, and coupled heat transfer impact on ternary nanofluid flow, a ternary nanofluid formed by dissolving $ A{{l}_{2}}{{O}_{3}},Si{{O}_{2}},Ti{{O}_{2}}$ nanoparticles in water. The analysis delves into the formulation of complex PDEs into a more manageable set of ODEs via application of similarity terms. Further examined the energy equation by using Kummer’s and Laguerre polynomial methods by incorporating a novel variable; the analysis encompasses both PST and PHF cases, and many physical parameters are studied through graphs; Table 2 defines the physical meaning and practical ranges of parameters.
Figure 2 represents the axial momentum graph for the $D{{a}^{-1}}$, upsurging the $ D{{a}^{-1}}$ reduces the movement of the fluid. The inverse Darcy number represents the extent to which the porous medium hampers the fluid flow; thus, higher values correspond to lower permeability and more intense drag effects. Increasing this parameter obstructs the flow through the porous medium, thus lowering the velocity and increasing the thickness of the momentum boundary layer. The physical significance is in the interaction of the liquid with the solid matrix of the porous medium, which may result in various flow patterns ranging from laminar to turbulent. Knowledge of such interactions is vital for the accurate prediction and efficient operation of systems that use filtration, oil recovery, and groundwater movement. In porous media situations where better heat transfer performance is required, the ternary nanofluid outperforms the hybrid nanofluid due to its improved thermal transport capacity. However, this improvement causes a slight increase in flow resistance due to the increased effective viscosity generated by the presence of more nanoparticles.

Figure 3 represents the axial velocity graph for the MHD; enhancing the magnetic term reduces the velocity of the fluid movement. When a $M$ term is applied to an electrically conducting fluid, it provides a Lorentz force that opposes the movement of the fluid, effectively altering its momentum. When a magnetic field is applied, a Lorentz force is produced, which acts against the movement of the electrically conducting fluid. Hence, by increasing the magnetic parameter, the fluid velocity is lowered and the momentum boundary layer thickness is also reduced. This interaction can be utilized to control and manipulate the flow in various applications, such as in cooling systems or magnetic pumps. Because of its higher effective thermal conductivity, the ternary nanofluid is more appropriate than the hybrid nanofluid under strong magnetic field situations that call for improved heat transfer. Because of the increased effective viscosity caused by the extra nanoparticles, this benefit results in a minor increase in flow resistance.

Figure 4 represents the axial velocity graph for the expanding/shrinking term; here solid dashed, and dotted lines denote ternary, hybrid, and nanofluid flow. Upsurging the stretching/shrinking term reduces the momentum; the stretching parameter qualitatively describes the extent to which the surface is stretched or shrunk; hence, it controls the growth of the boundary layer. Increasing the stretching parameter results in a higher surface-induced strain; this changes the velocity profile and causes a change in the thickness of the momentum boundary layer. The stretching/shrinking term significantly influences the momentum, particularly in boundary layer phenomena. When a surface stretches, it tends to accelerate the fluid particles near the surface, enhancing momentum due to the reduction of viscous effects. The ternary nanofluid is preferred to the hybrid nanofluid under higher stretching circumstances that need improved heat transfer due to its increased effective thermal conductivity.

Figure 5a and Figure 5b portray the axial momentum graph for the velocity slips; here solid, dashed, and dotted lines denote ternary, hybrid, and nanofluid movement respectively. upsurging the velocity slips term reduces the momentum of the fluid. The velocity slip is a description of a situation where the fluid at the surface is not fully sticking to the surface but can move slightly in the direction parallel to the wall; thus, a finite tangential velocity is allowed at the boundary. When the slip parameter is increased, the wall shear stress is diminished, resulting in less momentum being extracted from the surface and transferred to the fluid. Hence, the velocity boundary layer reduces in thickness. This parameter accounts for the relative motion between the liquid and the boundary, which can alter the flow's momentum. The ternary nanofluid outperforms the hybrid nanofluid under higher velocity slip situations due to its superior heat transfer capabilities. However, because of the increased effective viscosity brought on by the extra nanoparticles, this advantage causes a slight rise in resistance to flow.


Figure 6 portrays the axial velocity plot for the mass transpiration; here solid, dashed and dotted lines portray the ternary, hybrid, and nanofluid movement. Upsurging the mass transpiration reduces the velocity. The suction parameter of mass refers to the removal of fluid from the permeable surface, which leads to the stabilization of the boundary layer. When suction is raised, it makes the velocity boundary layer thinner, and at the same time it halts the motion of fluid near the surface due to the increased mass extraction. The mass transpiration term is crucial in controlling the boundary layer characteristics and ensuring the stability of the flow. Under the higher suction conditions, the ternary nanofluid performs better than the hybrid nanofluid because of its enhanced thermal transport properties. This gain comes with a modest increase in flow resistance since the added nanoparticles increase the fluid’s effective viscosity.

Figure 7 portrays the axial velocity plot for the viscosity ratio; here, solid, dashed, and dotted lines denote ternary, hybrid, and nanofluid flow. Upsurging the viscosity ratio term reduces the momentum of the fluid flow. The viscosity ratio reveals the share of the effective nanofluid viscosity in total viscosity and, hence, the change in flow resistance due to nanoparticle dispersion. It is well understood that a higher viscosity ratio increases the internal fluid friction that in turn lowers the velocity profile and makes the momentum diffusion effects more powerful. Physically, the ratio of viscous heat production to heat conduction in a moving viscous fluid is measured by the dimensionless viscosity ratio. It is frequently applied in polymer processing to measure the impact of heat transfer caused by viscous dissipation. Because of its superior thermal characteristics, the ternary nanofluid becomes more favorable than the hybrid nanofluid when higher thermal transfer performance is needed. However, the introduction of extra nanoparticles increases the fluid’s effective viscosity, resulting in a modest increase in flow resistance.

Figure 8a and Figure 8b represent the temperature plots for the Eckert term for both PST and PHF cases; here solid, dashed, and dotted lines portray ternary, hybrid, and nanofluid movement respectively. Eckert number quantifies the portion of kinetic energy that is transformed into internal energy by viscous dissipation. Under either PST or PHF, the fluid temperature is raised and the thermal boundary layer is thickened as the Eckert number increases because of the stronger internal heat generation; upsurging the Eckert number term raises the energy of the liquid movement. Physically, it denotes the ratio of the flow's kinetic energy to its enthalpy change, providing insight into the influence of fluid velocity on energy distribution. A high Eckert number denotes that the kinetic energy is significant compared to thermal energy, which can lead to notable temperature increases due to viscous dissipation. At higher Eckert number conditions, the ternary nanofluid exhibits superior thermal behavior compared with the hybrid nanofluid because it more effectively converts viscous dissipation into internal energy. Because of the increased effective viscosity brought on by the extra nanoparticles, this enhanced heat transfer is accompanied by a little increase in flow resistance.


Figure 9a and Figure 9b depict the temperature graphs for the $Nr$ for both PST and PHF cases; here solid, dashed, and dotted lines denote ternary, hybrid, and nanofluid flow respectively. Raising the thermal radiation raises the energy of the liquid movement. Figure 9c denotes the Nusselt number graph for thermal radiation. An increase the thermal radiation reduces the Nusselt number; a rise in the thermal radiation term typically results in a higher fluid temperature.



The thermal radiation parameter characterizes the impact of radiative heat transfer in comparison to conductive heat transfer within the fluid. When this parameter increases, energy transport by radiation is enhanced, which consequently leads to an increase in fluid temperature and a thicker thermal boundary layer. This is due to the enhanced energy absorption by the fluid, which can be particularly pronounced when there is a large energy difference between the liquid and its surroundings. Under strong thermal radiation conditions, the ternary nanofluid performs better than the hybrid nanofluid because of its higher effective thermal conductivity and improved heat transfer capability. Because of the increased effective viscosity brought on by the extra nanoparticles, this improvement is accompanied by a little increase in flow resistance.
Figure 10a and Figure 10b portray the energy plots for MHD on both PST and PHF cases; here solid, dashed, and dotted lines depict ternary, hybrid, and nanofluid, respectively. Raising the magnetic field raises the temperature of the fluid flow. In both heat boundary conditions, i.e., with PST and PHF, the magnetic field induces a resistive Lorentz force that converts the velocity energy into thermal energy via Joule heating. Hence, increasing the magnetic parameter leads to an increase in the temperature of the fluid and, consequently, to the enhancement of the thermal boundary layer thickness in both thermal boundary cases; magnetic fields can raise heat transfer in microscale flows between parallel plates, affecting the energy distribution. Furthermore, because of the combined thermal conductivity of many nanoparticles, the ternary nanofluid exhibits higher temperature enhancement under strong magnetic field circumstances than hybrid and mono nanofluids. However, because the ternary nanofluid’s effective viscosity is comparatively greater than the hybrid nanofluid’s, this enhanced thermal performance is accompanied by a little increase in flow resistance.


Figure 11a and Figure 11b portray the energy plots for the heat source/sink for both PST and PHF cases, enhancing the heat source/sink upsurge of the energy of the fluid. Figure 11c denotes the Nusselt number plot for the heat source/sink. As the heat source/sink increases, the Nusselt term decreases, which denotes the ratio of convective to conductive heat exchanging at the boundary of a fluid movement. In both PST and PHF scenarios, the presence of a positive heat source parameter results in additional heat inside the system; thus, the temperature of the fluid is raised. When the rate of heat generation increases, the temperature layer near the surface becomes more extended, whereas a heat sink results in the local cooling effect. An upsurge in the heat source intensity typically raises the fluid temperature and enhances the convective heat transfer, tending toward an upper Nusselt term. This signifies a more efficient heat transfer from the surface to the liquid, which is essential in designing effective cooling systems and understanding the thermal behavior of fluids. Because of its enhanced thermal conductivity and heat transfer ability, the ternary nanofluid is more appropriate than the hybrid nanofluid in situations of intense heat production and high thermal gradients. Nevertheless, this enhancement slightly increases the flow resistance because more nanoparticles increase the fluid’s effective viscosity.



Figure 12a and Figure 12b portray the skin friction graphs for the $D{{a}^{-1}}$ and MHD; here solid, dashed, and dotted lines represent ternary, hybrid, and nanofluids, upsurging the magnetic field and $D{{a}^{-1}}$ decreasing the skin friction of the fluid flow. Enhancing the magnetic field causes the Lorentz force to become stronger, which is a force that acts to resist the fluid motion and, as a result, changes the wall shear stress, thus causing an increase in skin friction. On the other hand, a bigger inverse Darcy number makes porous medium resistance more dominant, hence the momentum getting suppressed near the surface and the skin friction coefficient becoming larger. From a practical perspective, the ternary nanofluid provides better thermal performance than the hybrid nanofluid, although it leads to a modest rise in skin friction. Therefore, despite the minor increase in flow resistance, the ternary nanofluid is more appropriate for situations where improved heat transfer is preferred.


Figure 13 portrays the streamlined plots for various choices of viscosity ratio, MHD, and inverse Darcy term; here solid, dashed, and dotted lines denote ternary, hybrid, and nanofluid movement, respectively. Increasing the magnetic parameter will strengthen the Lorentz force, which, by enhancing the resistance to fluid motion, will change the wall shear stress and thus increase skin friction. Similarly, a larger inverse Darcy number will increase the level of resistance of the porous medium, which will elevate momentum suppression near the surface and thus increase the skin friction coefficient. The streamline snapshots illustrate that an increase in the viscosity ratio causes a significant reduction in the fluid circulation because of stronger internal friction, thus resulting in more compact and quite weakly flowing structures near the surface. The Lorentz force, which is the result of an increase in the magnetic parameter, acts against the primary flow direction; thus, the strength of the vortical motion is decreased, and the streamlines are pushed to the boundary. A higher inverse Darcy number similarly means that the porous medium offers a higher degree of resistance; this results in the suppression of fluid penetration and the diminishing of recirculation zones. A larger viscosity ratio, a stronger magnetic field, and an increased porous resistance are the factors that together exert such influence that the flow configuration gets a major damping, features of which include a reduced streamline spacing as well as lower momentum transport in general.






The comparative profiles of hybrid nanofluid ($A{{l}_{2}}{{O}_{3}}$-$Si{{O}_{2}}/{{H}_{2}}O$) and ternary nanofluid ($A{{l}_{2}}{{O}_{3}}$-$Si{{O}_{2}}$-$Ti{{O}_{2}}/\\{{H}_{2}}O$) show that the ternary suspension gives a higher thermal performance due to the synergistic effect of the three different types of nanoparticles. As seen also in Figure 9b and Figure 9c, the ternary nanofluid shows higher temperature levels and a different behavior of Nusselt number with the usage of thermal radiation; in other words, it is capable of stronger energy absorption and heat transport than the hybrid one.
Nevertheless, this thermal improvement goes hand-in-hand with the increased difficulty of the flow. Figure 12a and Figure 12b illustrate that the ternary nanofluid gives larger skin friction coefficients than the hybrid nanofluid at the same inverse Darcy number; the reason for this is the higher effective viscosity of the ternary mixture. As a consequence, a more powerful pump will be needed to maintain the flow. So, the findings disclose the existence of a trade-off in practice: on the one hand, ternary nanofluids possess outstanding heat transfer abilities; on the other hand, they cause a rise in hydrodynamic drag. Such fluids find their greatest use in a situation where the biggest heat rejection is the primary concern rather than the smallest flow resistance, e.g., in the cooling of very powerful electronics and in compact thermal systems.
4. Comparative Analysis of Prescribed Surface Temperature and Prescribed Heat Flux Boundary Conditions
Although both PST and PHF conditions depict the thermal transport at the wall, they portray fundamentally different thermal control mechanisms. In the PST case, the wall temperature is fixed, and the thermal boundary layer adjusts by altering the temperature gradient at the sheet. In contrast, under PHF cases, the wall heat flux is specified, and the sheet temperature becomes part of the solution.
This distinction leads to various thermal sensitivities. In the PST case, variations in parameters such as the Eckert number, radiation term, and heat source primarily impact the wall temperature gradient and thus the Nusselt number. However, in the PHF configuration, these terms influence both the wall temperature and the thermal gradient, producing a different quantitative energy transfer response. Therefore, the choice of boundary condition plays an important role in determining thermal boundary layer characteristics in magnetized ternary nanofluid flow.
A quantitative analysis of the temperature profiles reveals that the PHF arrangement changes its temperature more drastically with changes in the thermal parameters than the PST case. Thus, the radiation parameter change increases the temperature by only about 12% under PST, while the same parameter change increases the temperature under PHF by almost 18%. The same kind of trend is evidenced by the Eckert number, where viscous dissipation increases temperature by around 15% for PST but about 22% for PHF conditions.
The main reason for this difference is that, in the PST case, the wall temperature is kept constant, so the potential for thermal variation is limited, whereas in the PHF configuration, the wall temperature is allowed to change in accordance with the heat flux, thereby amplifying the overall thermal response. Thus, PHF boundary conditions result in thermal components being more sensitive to slight changes, while PST conditions offer temperature regulation with less fluctuation.
Table 4 shows a quantitative comparison of the percentage variations in the temperature distribution for both PST and PHF conditions, as well as the changes in the skin friction coefficient and Nusselt number for specific governing parameters. The magnetic parameter $M$ improves thermal responsiveness by 20% in the PST case and 25% in the PHF case, while also increasing skin friction coefficient by roughly 30% owing to resistive Lorentz force. The Eckert number $Ec$ which represents viscous dissipation, leads to a considerable increase in temperature in both thermal boundary conditions. Furthermore, the inverse Darcy number raises the skin friction coefficient by around 32%, suggesting higher resistance provided by the porous medium. The thermal radiation parameter $Nr$ and heat source/sink parameter $Ni$ have significant influence on the temperature field and Nusselt number, illustrating the sensitivity of heat transfer characteristics to these factors.
Parameter | Temperature (%) | Skin Friction (%) | Nusselt Number (%) | |
PST | PHF | |||
$M$ (Magnetic parameter) | 20% | 25% | 30% | - |
$Ec$ (Eckert number) | 15% | 22% | - | - |
$Da^{-1}$ (Inverse Darcy number) | - | - | 32% | - |
$Nr$ (Thermal radiation) | 12% | 18% | - | 9% |
$N_i$ (Heat source/sink) | 27% | 24% | - | 26% |
Although closed-form analytical expressions were gained for the momentum and temperature fields, the evaluation of these solutions across a wide parameter range requires systematic computational implementation. The roots of the momentum equation were first computed for each term set to ensure physically admissible reducing solutions. The analytical solutions involving exponential functions and generalized Laguerre polynomials were then evaluated analytically.
Parameter sweeps were performed by changing one parameter at a time while keeping others fixed within physically realistic ranges, as summarized in Section 4.1. For each term combination, the momentum profile, energy distribution, reduced skin friction coefficient, and Nusselt number were calculated. Graphical outcomes were produced by discretizing the similarity variable across a sufficiently large domain to ensure asymptotic convergence of the boundary layer.
Standard analytical methods were used for the root’s evaluation and special functions implementation in all calculations. It was ensured that the temperature solution converges for various values of the radiation and dissipation terms. This efficient computational protocol makes the analytical platform accessible for scientific-type parametric analysis.
5. Conclusion
The present study examined the behavior of dissipative ternary nanofluid flow over a permeable surface in the presence of viscous dissipation, Joule heating, and coupled heat-transfer effects. The ternary nanofluid, composed of $ A{{l}_{2}}{{O}_{3}}$ , $ Si{{O}_{2}}$ , and $ Ti{{O}_{2}}$ nanoparticles dispersed in water, was analysed under both PST and PHF conditions. By means of suitable similarity transformations, the governing equations were reduced to a tractable form, which made it possible to assess the influence of the main physical parameters on the velocity and thermal fields.
The principal findings may be summarised as follows:
• An increase in the magnetic parameter suppresses the fluid velocity while increasing the skin-friction coefficient and the thermal field in both PST and PHF cases.
• A stronger porous-medium resistance reduces the momentum boundary layer and increases wall shear.
• Higher mass transpiration and stronger stretching/shrinking effects reduce the fluid momentum.
• An increase in the viscosity-ratio parameter weakens the velocity boundary layer.
• Larger values of the Eckert number enhance the temperature field under both PST and PHF boundary conditions.
• Thermal radiation increases the temperature distribution and reduces the Nusselt number.
• A stronger heat source parameter raises the fluid temperature.
Overall, the results confirm that viscous dissipation and Joule heating have a pronounced influence on the thermal response of the system, whereas magnetic and porous-medium effects play an important role in controlling the flow structure. The combined action of these parameters determines the balance between heat-transfer enhancement and flow resistance in ternary nanofluid transport over permeable surfaces.
The present results may be useful in the design of magnetically controlled thermal systems involving ternary nanofluids in porous structures. In particular, the analysis shows that parameters such as viscous dissipation and thermal radiation can significantly increase the fluid temperature, which may be beneficial in processes where enhanced thermal transport is required. At the same time, stronger magnetic effects and greater porous resistance provide additional means of flow control, although these advantages are accompanied by reduced velocity and increased wall shear, which may lead to a higher pumping-power requirement.
The ternary nanofluid considered here is expected to offer better thermal performance than conventional or hybrid nanofluids, but this improvement is not without cost. The gain in heat-transfer capability is associated with increased hydrodynamic resistance, and this trade-off must be considered in practical design. From an engineering point of view, the best operating condition is therefore not obtained simply by maximising thermal enhancement but by selecting parameter combinations that improve heat removal without causing excessive flow resistance or undesirable thermal buildup near the surface.
For this reason, the present model may serve as a useful predictive framework for the preliminary design and optimisation of advanced cooling technologies, porous heat exchangers, and other thermal-management systems in which coupled magnetic, porous, and dissipative effects are relevant.
The present analysis is limited to a steady, laminar, two-dimensional boundary-layer flow model based on a single-phase representation of the ternary nanofluid. The thermophysical properties are assumed to be constant, and the induced magnetic field is neglected. These assumptions make the problem analytically tractable, but they also restrict the direct applicability of the model to more complex operating conditions.
In particular, when the nanoparticle concentration becomes high, when the magnetic field is very strong, or when non-Darcy porous effects become significant, the present formulation may no longer provide an adequate physical description. Under such conditions, a fully coupled numerical treatment beyond the boundary-layer approximation may be required. In addition, the single-phase effective-property approach does not account for particle interaction, agglomeration, slip mechanisms, or local thermal non-equilibrium effects, all of which may become important at higher particle loadings.
Future work may therefore consider unsteady flow, variable thermophysical properties, non-Newtonian base fluids, entropy-generation analysis, and experimental validation. Such extensions would provide a broader basis for assessing the practical relevance of the present model and for improving its applicability to real engineering systems.
The analytical solutions obtained in this study are valid within the framework of laminar boundary-layer theory and for parameter ranges consistent with the assumptions used in the derivation. They are therefore most appropriate for physically moderate regimes in which the flow remains steady, two-dimensional, and sufficiently close to the boundary-layer structure adopted in the analysis.
If the magnetic field becomes excessively strong, or if viscous dissipation and Joule heating become dominant, the coupled nonlinear effects may become too pronounced for the present analytical framework to remain fully reliable. In such cases, the flow may depart from the assumed boundary-layer behaviour, and a fully numerical approach would be more appropriate. Similarly, the Darcy-Brinkman porous-medium model adopted here assumes homogeneous permeability and cannot represent more complex situations involving anisotropic, heterogeneous, or strongly non-Darcy porous structures.
It should also be noted that the ternary nanofluid has been modeled through an effective-property single-phase approach. At higher particle concentrations, additional physical mechanisms such as inter-particle interaction, agglomeration, slip velocity, and thermal non-equilibrium may become important and would require more advanced multiphase or two-component formulations. Accordingly, the present analytical solutions should be interpreted as valid for the class of physical situations described above and not as universally applicable to all ternary nanofluid systems.
Conceptualization, K.N. and S.M.S.; methodology, S. M.S. and U.S.M.; validation, U.S.M. and L.M.P.; formal analysis, S.M.S.; investigation, S.M.S.; resources, S.M.S and K.N.; data curation, S.M.S and K.N.; writing—original draft preparation, S.M.S.; writing—review and editing, U.S.M. and G.L.; visualization, K.N.; supervision, U.S.M.; project administration, G.L. and U.S.M. All authors were actively involved in discussing the findings and refining the final manuscript.
The data used to support the findings of this study are available from the corresponding author upon request.
The authors declare that they have no conflicts of interest.
$A$ | First order slip coefficient |
$A_{1}, A_{2}, A_{3}, A_{4}, A_{5}$ | Constants |
$\mathrm{Al_{2}O_{3}}$ | Aluminium oxide |
$B$ | Second order slip coefficient |
$\mathrm{SiO_{2}}$ | Silicon dioxide |
$\mathrm{TiO_{2}}$ | Titanium dioxide |
$C_{p}$ | Specific heat |
$D$ | Mass flux |
$Da^{-1}$ | Inverse Darcy number |
$a$ | Stretching rate |
$c$ | Acceleration rate |
$B_{0}$ | Magnetic strength |
$b, l$ | Constants |
$Ec_{1}, Ec_{2}$ | Eckert numbers |
$K^{*}$ | Permeability of porous medium |
$M$ | Magnetic field parameter |
$Nr$ | Thermal radiation parameter |
$N_i$ | Heat source/sink parameter |
$\Lambda$ | Viscosity ratio |
$\Pr$ | Prandtl number |
$\operatorname{Re}_{x}$ | Reynolds number |
$q_{r}$ | Radiative heat flux |
$S$ | Mass suction/injection |
$d$ | Stretching/shrinking parameter |
$f$ | Velocity function |
$Q_{0}$ | Heat source/sink coefficient |
$T_{w}$ | Surface temperature |
$T_{\infty}$ | Ambient temperature |
$T$ | Fluid temperature |
$u, v$ | Velocity components |
$u_{w}, v_{w}$ | Wall velocities |
$x, y$ | Coordinate axes |
Greek Symbols
$\beta$ | Solution term |
$\eta$ | Similarity variable |
$\rho$ | Fluid density |
$\psi$ | Stream function |
$\kappa$ | Thermal conductivity |
$k^{*}$ | Absorption coefficient |
$\mu$ | Dynamic viscosity |
$\nu$ | Kinematic viscosity |
$\sigma$ | Electrical conductivity |
$\sigma^{*}$ | Stefan--Boltzmann constant |
$\theta$ | PST temperature |
$\varphi$ | PHF temperature |
$\phi$ | Nanoparticles volume fraction |
$\xi_{1}$ | First order velocity slip parameter |
$\xi_{2}$ | Second order velocity slip parameter |
Subscripts
HNF | Hybrid nanofluid |
B.Cs | Boundary conditions |
PDE | Partial differential equation |
ODE | Ordinary differential equation |
MHD | Magnetohydrodynamics |
PHF | Prescribed heat flux |
PST | Prescribed surface temperature |
