Mathematical Analysis of Two Phase Saturated Nanofluid Influenced by Magnetic Field Gradient

Abstract

Nanofluids are composed of nano-sized particles dispersed in a carrier liquid. The present investigation’s aim is to examine theoretically the magneto-thermomechanical coupling phenomena of a heated nanofluid on a stretched surface in the presence of magnetic dipole impact. Fourier’s law of heat conduction is used to evaluate the heat transmission rate of the carrier fluids ethylene glycol and water along with suspended nanoparticles of a cobalt–chromium–tungsten–nickel alloy and magnetite ferrite. A set of partial differential equations is transformed into a set of non-linear ordinary differential equations via a similarity approach. The computation is performed in Matlab by employing the shooting technique. The effect of the magneto-thermomechanical interaction on the velocity and temperature boundary layer profiles with the attendant effect on the skin friction and heat transfer is analyzed. The maximum and minimum thermal energy transfer rates are computed for the H₂O-Fe₃O₄ and C₂H₆O₂-CoCr₂₀W₁₅Ni magnetic nanofluids. Finally, the study’s results are compared with the previously available data and are found to be in good agreement.

1. Introduction

The problems of fluids’ flow and heat transfer for the rate of cooling and heating is very important in various industrial processes, such as polymer extrusion, drawing of plastic and sheets, fiberglass, the production of paper and many more. The transfer of heating or cooling in an object or between two distinct surfaces takes place due to temperature differences. A higher temperature difference causes higher thermal energy transmission [1,2,3]. The transmission of thermal energy can be computed theoretically using the well known law, i.e., Fourier’s law. This law uses the correlation between thermal conductivity and heat capacity to evaluate the transmission rate of thermal energy. The law also states that the transmission of thermal energy occurs with an inertial rate [4,5,6]. Later on, Cattaneo-Christov [7] added the thermal relaxation time into Fourier’s law. It plays an important role when experiments are difficult to conduct for heat transfer simulation. These laws are used for the theoretical computation of heat transfer [8,9,10]. The above discussed laws are used for evaluation of heat flux in two phase fluids as well as regular fluids.

The boundary layer flows on the stretching sheet has been examined by many researchers for the rate of heat transfer and cooling as it has promising applications in many industries. The stretching sheet velocity of the flow field has gained considerable attention as it has a certain effect on the excellence of the final products in the fabric and polymer industry. This stretching velocity can be linear, polynomial, hyperbolic or exponential. The pioneering work of Crane [11] has laid the foundation for theoretical investigations on the flow field over a stretching surface. He presented an exact analytical solution of the flow over a linearly stretching sheet. The boundary layer heat transport flow of multiphase magnetic fluids past a stretching sheet under the impact of a circular magnetic field was described in [12]. The general fluid model for the magnetohydrodynamic fluid flow and heat transfer was analyzed by Hatzikonstantinou and Vafeas [13]. The authors computed heat transfer along with the skin friction coefficient in the flow of a micropolar fluid. In this direction, Papadopoulos et al. [14] demonstrated the theoretical investigation of pipe flow in the presence of a magnetic field in a ferrofluid along with a cylindrical coil. The available literature in this direction includes the study of ferrite nanoparticles [15,16,17,18,19].

The literature consists of study of ferrite and regular particles, which are studied in different base fluids for heat transfer along with skin friction. These particles have been studied multiple times with different effects to compute which mathematical model better suits the physical situation. Thus, in the field of fluid dynamics the same problem is normally studied with different effects and parameters to predict the heat transfer that better suits the physical situation modeled. The fluid flow computed for heat transfer and skin friction is available in [20,21,22,23].

Ferromagnetic fluids were first synthesized by Rosensweig [24] in his pioneering work, in which he also first used the term ferrohydrodynamics. After that, the work has received increasing attention over the last few decades. In his later detailed work [25], he explored the basic equations for ferromagnetic fluids with internal rotation. The heat transfer is efficient in solids instead of liquids or gasses. Therefore, smart liquids called ferrofluids are synthesized by introducing ferromagnetic nanoparticles into a carrier fluid. The suspension of solid nanoparticles (1–100 nm) in a base fluid enhances the transport properties of the considered nanofluid and hence its heat transfer rate. The ferrofluids are famous due to the suspension of nano-sized ferrite particles instead of regular particles. These ferrite particles have the extra property of Curie temperature; thus, ferrofluids are of great interest and have been studied extensively by researchers after the work of Rosensweig [24]. To make electronic devices long-lasting, ferrofluids can evacuate heat from these devices, such as in speakers and Lenovo laptops [26]. Heat transfer in these devices takes place through the magnetic fluid, made up of a base fluid (water, oil, ethylene glycol) and ferrite nanoparticles. Later on, the flow problem discussed by Crane [11] was extended by Andreson and Valnes [27] with the addition of a magnetic field gradient on the stretched surface of an incompressible non-conducting ferrofluid. Odenbach [28] pointed out the behavior of magnet field impact over magnetite ferrofluid experimentally and discussed the practical applications in technical and medical fields.

In this study, the dynamics of heat transfer processes associated with the motion of viscous, incompressible fluids in the presence of a magnetic field and temperature gradient are analyzed. One of the striking features of the considered nanofluid is the dependence of the magnetization upon the temperature gradient, and this thermomagnetic coupling makes the study of these fluids demanding. Since the suspension of nano-sized particles in a viscous fluid enhances the transport properties of the flowing fluid, the transmission rate of the thermal energy is also enhanced. Therefore, the aim of the article is to describe the interaction between a magnetic field gradient and solid nano-sized particles and its influence on the convection, heat transfer and friction drag.

A cobalt–chromium–tungsten (wolfram)–nickel (CoCr₂₀W₁₅Ni) alloy and magnetite ferrite (Fe₃O₄) nano-sized particles were added into a base fluid of ethylene glycol (C₂H₆O₂) and water (H₂O). The graphical results are computed in the discussion section. The particles and the base fluid properties are taken under the assumption of isothermal equilibrium. The analysis is made in such a way that a single solid particle is suspended in the base fluid. The suspension is then computed for thermal energy transmission. The comparison is made for the ferrite (Fe₃O₄) and alloy (CoCr₂₀W₁₅Ni) two phase nanofluid flow. The results are discussed and presented for the base fluids ethylene glycol (C₂H₆O₂) and water (H₂O).

2. Mathematical Modeling

The steady state two phase incompressible laminar flow in the presence of an external magnetic dipole that is placed near the x surface at distance l is incorporated. The disturbance to the laminar flow is induced via the stretching x−surface sheet, and the stretching velocity is defined to be $u=Sx$. The temperature of the fluid at the surface is $T_w$, whereas the temperature of fluid far away from the sheet is taken to be a Curie temperature, i.e., $T=T_c$. The ferrite Fe₃O₄ and alloy CoCr₂₀W₁₅Ni solid nanoparticles satisfy the Curie temperature $T_c$. By contrast, the ferrite Fe₃O₄ and alloy CoCr₂₀W₁₅Ni along with the base fluids ethylene glycol C₂H₆O₂ and water H₂O are taken under the assumption of isothermal equilibrium and no slip occurs between them.

The vector form of the equation for the flow of the ferrofluid is given in [25]. In vector form the equations for the present problem are:

$$\nabla .\mathbf{V}=0,$$

(1)

$$\begin{array}{cc}& {\rho }_{nf}(\frac{\partial \mathbf{V}}{\partial t}+\mathbf{V}.\nabla \mathbf{V})=-\nabla P+{\mu }_{nf}{\nabla }^2\mathbf{V}+{\mu }_{0}M.\nabla \mathbf{H},\hfill \end{array}$$

(2)

$$\begin{array}{cc}& {\left(\rho c_p\right)}_{nf}(\frac{\partial T}{\partial t}+\mathbf{V}.\nabla T+{\mu }_{0}T\frac{\partial M}{\partial T}\mathbf{V}.\nabla \mathbf{H})=k_{nf}{\nabla }^{2}T,\hfill \end{array}$$

(3)

Physically, Equation (1) represents the conservation of mass. The left-hand side of Equation (2) represents momentum due to convection or inertial forces; the first term in Equation (2) on the right-hand side represents the pressure, whereas the second term appears due to surface forces or viscous forces, and the third term denotes the body forces due to the magnetic dipole moment. The first term on the left-hand side of Equation (3) represents heat transfer due to convection, and the second term represents the internal energy due to magnetic dipole. The right-hand side of Equation (3) shows the heat transfer due to conduction.

The mathematical equation for the steady state two-dimensional flow of an incompressible ferrofluid and boundary layer flow described by the dynamical equations governing convective flow and the heat transfer of the nanofluid can be written as:

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

(4)

$$\begin{array}{cc}& {\rho }_{nf}(u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y})=-\frac{\partial P}{\partial x}+{\mu }_{nf}\frac{{\partial }^{2}u}{\partial y^2}+{\mu }_{0}M\frac{\partial H}{\partial x},\hfill \end{array}$$

(5)

$$\begin{array}{cc}& {\left(\rho c_p\right)}_{nf}(u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y})+{\left(\rho c_p\right)}_{nf}{\mu }_{0}T\frac{\partial M}{\partial T}(u\frac{\partial H}{\partial x}+v\frac{\partial H}{\partial y})=k_{nf}\frac{{\partial }^{2}T}{\partial y^2}.\hfill \end{array}$$

(6)

The thermo-mechanical coupling at the wall and far from the wall is defined as:

$${\left.u\right|}_{y=0}=Sx,{\left.v\right|}_{y=0}=0,{\left.T\right|}_{y=0}=T_w,$$

(7)

$${\left.u\right|}_{y\to \infty }\to 0,{\left.T\right|}_{y\to \infty }\to T_c.$$

(8)

The term ${\mu }_{0}M\frac{\partial H}{\partial x}$ in Equation (5) represents the component of magnetic force per unit volume. This term depends on the existence of the magnetic field gradient along the corresponding x direction, called the Kelvin force, and is very well known in ferrohydrodynamics. The term ${\mu }_{0}T\frac{\partial M}{\partial T}(u\frac{\partial H}{\partial x}+v\frac{\partial H}{\partial y})$ in Equation (6) is due to the magneto caloric effect and represents the thermal power per unit of volume.

The relation of the ferrofluid between some of the physical characteristics of the solid nanoparticle and the considered carrier fluid were explained by Islam et al. [29] and is given in

Table 1. Mathematical expressions of thermo-physical properties.

Properties	Islam et al. [29]
Thermal Conductivity k	$\frac{k_{nf}}{k_f}=\frac{2k_f+k_s-2\mathsf{\Phi }(k_f-k_s)}{2k_f+k_s+\mathsf{\Phi }(k_f-k_s)}$
Viscosity $\mu$	${\mu }_{nf}={\mu }_f{(1-\mathsf{\Phi })}^{-25/10}$
Heat Capacity $\rho C_p$	${\left(\rho c_p\right)}_{nf}=(1-\mathsf{\Phi }){\left(\rho c_p\right)}_f+\mathsf{\Phi }{\left(\rho c_p\right)}_s$
Density $\rho$	${\rho }_{nf}=(1-\mathsf{\Phi }){\rho }_f+\mathsf{\Phi }{\rho }_s$

:

The thermo-physical properties for C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni are given in

Table 2. Thermo-physical values of C₂H₆O₂, H₂O, CoCr₂₀W₁₅Ni, and Fe₃O₄.

	k (W/mK)	$\mathit{\rho }$ (kg/m3)	${\mathit{C}}_{\mathit{p}}$ (J/kgK)
CoCr20W15Ni	15	9.60	400.0
Fe3O4	9.7	5180	670
C2H6O2	0.249	1116.6	2382
H2O	0.60	998.3	4182

.

2.1. Magnetic Dipole

Artificially synthesized nanofluids behave like normal fluids except that they experience a force due to magnetization. One of the striking features of nanofluids is the dependence of magnetization upon the temperature gradient, and this thermomagnetic coupling makes these fluids demanding in various applications. An external magnetic dipole is placed in the flowing fluid to magnetize the two phase C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni nanofluids. The introduction of magnetic forces into the magnetizable liquid gives rise to the effect known as ferrohydrodynamics. Now, since the nanoparticles are mechanically free to align with the fields of lines, the expression of the body force could be a good approximation for the flows of nanofluids. The magnetic scalar potential $\mathsf{\Omega }$ is stated below:

$$\mathsf{\Omega }=\frac{{\gamma }_1}{2\pi }\frac{x}{{(y+l)}^2+x^2},$$

(9)

here ${\gamma }_1$ symbolizes the dipole moment per unit length, while the x and y components for the magnetic field are:

$$H_x=-\frac{\partial \mathsf{\Omega }}{\partial x}=-\frac{{\gamma }_1}{2\pi }\frac{-{(y+d)}^2+x^2}{{\left(x^2+{(y+d)}^2\right)}^2},$$

(10)

$$H_y=-\frac{\partial \mathsf{\Omega }}{\partial y}=\frac{{\gamma }_1}{2\pi }\frac{2x(y+d)}{{\left({(y+d)}^2+x^2\right)}^2}.$$

(11)

Since it is generally known that the magnetic body force is relative to the magnetic field gradient of H, we have:

$$H=\sqrt{{\left(\frac{\partial \mathsf{\Omega }}{\partial y}\right)}^2+{\left(\frac{\partial \mathsf{\Omega }}{\partial x}\right)}^2}.$$

(12)

Using Equations (10)–(12), and then differentiating the resultant equation with respect to x and y, we obtain the following equations:

$$\frac{\partial H}{\partial x}=-\frac{{\gamma }_1}{\pi }\frac{x}{{(y+d)}^4},$$

(13)

$$\frac{\partial H}{\partial y}=\frac{{\gamma }_1}{\pi }\left(\frac{2x^2}{{(y+d)}^5}-\frac{1}{{(y+d)}^3}\right).$$

(14)

The variation of the magnetization M with the magnetic field intensity H and temperature T can be fairly approximated by the linear relation shown below. The detailed explanation for the derivation of H is presented in [27]. The magnetization is defined as:

$$M=K_c(T-T_{\infty }).$$

(15)

The geometry for the two phase C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni nanofluids is evident in Figure 1.

2.2. Similarity Analysis

For the numerical solution of the governing boundary layer model, Equations (5) and (6), along with boundary conditions, need to be transformed into non-dimensional ordinary differential equations. The transformations is performed by introducing the similarity analysis. The similarity transformation for the two phase nanofluid introduced by Andersson and Valnes [27] is:

$$\begin{array}{cc}& \psi (\xi ,\eta )=\left(\frac{{\mu }_f}{{\rho }_f}\right)\eta f\left(\xi \right),\hfill \\ & \theta \left(\eta ,\xi \right)\equiv \frac{T-T_c}{T_w-Tc}={\theta }_1\left(\xi \right)+{\eta }^2{\theta }_2\left(\xi \right).\hfill \end{array}$$

(16)

The stream function $\psi (x,y)$ defined by the velocity components is $u=\frac{\partial \psi }{\partial y},v=-\frac{\partial \psi }{\partial x}$, which satisfies the conservation Equation (4) identically, while f and $\theta$ are the dimensionless functions.

$$\begin{array}{cc}& \xi =y(\frac{{\rho }_{f}S}{{\mu }_f}{)}^{\frac{1}{2}},u=Sx{f}^{\prime }\left(\xi \right),\hfill \\ & \eta =x(\frac{{\rho }_{f}S}{{\mu }_f}{)}^{\frac{1}{2}},v=-(S{\nu }_f{)}^{\frac{1}{2}}f\left(\xi \right).\hfill \end{array}$$

(17)

By using Equations (16) and (17) in the momentum and energy Equations (5) and (6), we obtain the transformed nonlinear ordinary differential equations of the momentum and temperature as follows:

$$\frac{A_1^{-1}}{{(1-\mathsf{\Phi })}^{\frac{25}{10}}}{f}^{‴}+f{f}^{″}-{\left(f^{\prime }\right)}^2-\frac{2A_1^{-1}\beta {\theta }_1}{{(\xi +\gamma )}^4}=0,$$

(18)

$$\frac{A_2^{-1}{k}_{nf}}{k_f}{\theta }_1^{″}+Pr(f{\theta }_1^{\prime }-2f^{\prime }{\theta }_1)-4A_2^{-1}{\sigma }_2{\left(f^{\prime }\right)}^2+\frac{2A_2^{-1}{\sigma }_2\beta f({\theta }_1-{\sigma }_1)}{{(\xi +\gamma )}^3}=0,$$

(19)

$$\begin{array}{cc}& \frac{A_2^{-1}{k}_{nf}}{k_f}{\theta }_2^{″}-Pr(4f^{\prime }{\theta }_2-f{\theta }_2^{\prime })+\frac{2A_2^{-1}{\sigma }_2\beta f{\theta }_2}{{(\xi +\gamma )}^3}\hfill \\ & -\frac{\beta ({\theta }_1-{\sigma }_1)}{{(\xi +\gamma )}^3}(\frac{4f}{{(\xi +\gamma )}^5}+\frac{2f^{\prime }}{{(\xi +\gamma )}^4})-A_2^{-1}{\sigma }_2{\left(f^{″}\right)}^2=0,\hfill \end{array}$$

(20)

The boundary conditions in (6) and (7) now become:

$$\begin{array}{cc}& f^{\prime }\left(\xi \right)=1,f\left(\xi \right)=0,{\theta }_1\left(\xi \right)=1,{\theta }_2\left(\xi \right)=0,at\xi =0,\hfill \\ & f^{\prime }\left(\xi \right)\to 0,{\theta }_1\left(\xi \right)\to 0,{\theta }_2\left(\xi \right)\to 0,at\xi \to \infty ,\hfill \end{array}$$

(21)

while $A_1$ and $A_2$ are:

$$A_1=1-\mathsf{\Phi }+\mathsf{\Phi }\frac{{\rho }_s}{{\rho }_f},A_2=1-\mathsf{\Phi }+\mathsf{\Phi }\frac{\rho c_p{}_s}{\rho c_p{}_f}.$$

(22)

The dimensionless parameters calculated in the problem include $\beta$ (ferrohydrodynamic interaction), ${\sigma }_2$ (viscous dissipation), $Pr$ (Prandtl number), and ${\sigma }_1$ (dimensionless Curie temperature). Mathematical expressions are presented below:

$$\begin{array}{cc}& {\sigma }_1=\frac{T_c}{T_w-T_c},\beta =\frac{{\gamma }_1}{2\pi }\frac{{\mu }_0{K}_c(T_w-T_c){\rho }_f}{{\mu }_f^2},\hfill \\ & {\sigma }_2=\frac{S{\mu }_f^2}{{\rho }_f{K}_c(T_w-T_c)},\gamma =\sqrt{\frac{S{\rho }_f{l}^2}{{\mu }_f}},Pr=\frac{{\nu }_f}{{\alpha }_f}.\hfill \end{array}$$

(23)

The physical interest parameters for the boundary layer flow problems are the skin friction coefficient, which is the rate of shear stress, and the Nusselt number, which is the rate of heat transfer. The wall shear stress characterized by the skin friction coefficient $C_f$, along with the Nusselt number, is:

$$\begin{array}{cc}& C_f=\frac{{\tau }_w}{\frac{1}{2}{\rho }_{nf}{S}^2},{\tau }_w={\mu }_{nf}{\left.\frac{\partial u}{\partial y}\right|}_{y=0},\hfill \\ & N{u}_x={\left.\frac{x{k}_{nf}}{k_f(T_w-T_c)}\frac{\partial T}{\partial y}\right|}_{y=0}.\hfill \end{array}$$

(24)

where $q_w=-k_{nf}{\left(\frac{\partial T}{\partial y}\right)}_{y=0}$ is the surface heat flux. The substitution of the non-dimensional similarity transformation considered in Equations (16) and (17) into Equation (24) transforms the skin friction and local Nusselt number into dimensionless form. Hence, Equation (24) takes the following form:

$$\begin{array}{cc}& \frac{1}{2}R{e}_x^{\frac{1}{2}}{C}_f=\frac{1}{{(1-\mathsf{\Phi })}^{\frac{25}{10}}}{f}^{″}\left(0\right),\hfill \\ & R{e}_x^{-\frac{1}{2}}N{u}_x=\frac{-k_{nf}}{k_f}\left({\theta }_1^{\prime }\left(0\right)+{\eta }^2{\theta }_2^{\prime }\left(0\right)\right).\hfill \end{array}$$

(25)

where the local Reynolds number $R{e}_x$ is defined as $R{e}_x=xS/{\nu }_f$, which is based on the stretching sheet velocity $S_{\left(x\right)}$.

3. Numerical Simulation

The non-linear ordinary differential momentum (Equation (18)) and energy equations (Equations (19) and (20)), together with the appropriate boundary conditions (Equations (21)) were solved numerically using the shooting technique for various values of physical parameters. The mathematical model for the ferromagnetic two phase C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni nanofluids was computed with the help of the shooting technique in Matlab. The analysis in the discussion section is based on the dimensionless parameters and numbers in the above equations. The results show the impact of the solid particles on the heat transfer rate and friction drag theoretically. The system of equations is transformed into the first order equations in order to solve via the shooting technique.

$$\begin{array}{cc}& W_1=f\left(\xi \right),W_2=f^{\prime }\left(\xi \right),W_3=f^{″}\left(\xi \right),W_3^{\prime }=f^{‴}\left(\xi \right),W_4={\theta }_1\left(\xi \right),\hfill \\ & W_5={\theta }_1^{\prime }\left(\xi \right),W_5^{\prime }={\theta }_1^{″}\left(\xi \right),W_6={\theta }_2\left(\xi \right),W_7={\theta }_2^{\prime }\left(\xi \right),W_7^{\prime }={\theta }_2^{″}\left(\xi \right).\hfill \end{array}$$

(26)

Using the transformation given in Equation (26), it reduces the Equations (18)–(21):

$$\begin{array}{cc}& W_3^{\prime }={(1-\mathsf{\Phi })}^{\frac{25}{10}}{A}_1(W_2*W_2-W_3*W_1+\frac{2\beta W_4}{A_1{\left(\xi +\gamma \right)}^4}),\hfill \\ & W_5^{\prime }=-\frac{k_f{A}_2}{k_{nf}}(Pr(W_5{W}_1-2W_4{W}_2)+\frac{2\lambda \beta W_4(W_1-\epsilon )}{A_2{(\xi +\gamma )}^3}-\frac{4\lambda }{A_2}{W}_2^2),\hfill \\ & W_7^{\prime }=\frac{k_f{A}_2}{k_{nf}}(Pr(4W_6{W}_2-W_7{W}_1)-\frac{2\lambda \beta W_1{W}_6}{A_2{(\xi +\gamma )}^3}+\frac{\lambda }{A_2}{W}_2^2\hfill \\ & +\frac{\lambda \beta (W_4-\epsilon )}{{(\xi +\gamma )}^3}(\frac{2W_2}{{(\xi +\gamma )}^4}+\frac{4W_1}{{(\xi +\gamma )}^5})),\hfill \end{array}$$

(27)

$$\begin{array}{cc}& W_1=1,W_2=0,W_4=1,W_6=0,at\xi =0,\hfill \\ & W_2\to 0,W_4\to 0,W_6\to 0,at\xi \to \infty .\hfill \end{array}$$

(28)

4. Discussion

The two phase C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄O-CoCr₂₀W₁₅Ni magnetic nanofluids were computed via the shooting technique in Matlab. The influence of physical parameters and dimensionless numbers on the flowing magnetic nanofluid is discussed in this section. Two base fluids along with two distinct solid particles are examined. The ferrite and alloy nanoparticles are entertained in the analysis. The magnetic two phase C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni nanofluids are simulated. The flows of C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni nanofluids were computed and then compared with each other. The variation of velocity field was noticed to be in an order, H₂O-Fe₃O₄, H₂O-CoCr₂₀W₁₅Ni, C₂H₆O₂-Fe₃O₄, and C₂H₆O₂-CoCr₂₀W₁₅Ni. Therefore, we conclude that the maximum velocity field is determined for the H₂O-Fe₃O₄ nanofluid and the minimum distribution of velocity is determined for the magnetic C₂H₆O₂-CoCr₂₀W₁₅Ni nanofluid. This variation occurs due to the thermophysical values of the transport properties of the solid particles and base fluids. The comparison is shown in Figure 2. The relation is compared for the temperature field in Figure 3. The variation in temperature distribution is noticed in the order C₂H₆O₂-CoCr₂₀W₁₅Ni, C₂H₆O₂-Fe₃O₄, H₂O-CoCr₂₀W₁₅Ni, and H₂O-Fe₃O₄. This means that the H₂O-Fe₃O₄ nanofluid leads to the lowest temperature, whereas the C₂H₆O₂-CoCr₂₀W₁₅Ni magnetic nanofluid exhibits the higher temperature field. The ferrite nanoparticles considered in the problem have higher thermal properties, which help in the enhancement of heat transfer; as a result, the lower temperature field occurs for the magnetite ferrite base magnetic nanofluid. Cobalt particles, however, have the property of Curie temperature, but that is not sufficient for a higher temperature field. Thus, the lower temperature field is observed for the alloy-based magnetic nanofluid.

The influence of ferrohydrodynamic interaction parameter $\beta$ on the temperature and the axial velocity is displayed in Figure 4 and Figure 5. The magnetic property of the solid particles, i.e., alloys and ferrite, makes the two phase nanofluid magnetize in the presence of a magnetic dipole placed near the surface. The dipole attracts the alloy and the ferrite particles until the alloy and ferrite reach their Curie temperature. Thus, the axial velocity for varying $\beta$ declines, whereas the temperature profile for $\beta$ shows the reverse nature for all of the magnetic nanofluids C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni. The impact of $\beta$ on the axial velocity is determined in Figure 4, whereas its influence on the temperature field is evident in Figure 5. The interaction between the nanoparticles of the ferrite and alloys with the flowing fluid leads to resisting the flow and the internal energy of the fluid. The higher resistance to the flowing fluid is induced by the alloy as well as the ferrite nanoparticles; as a result, the velocity field declines, as presented in Figure 4. On the other hand, the internal energy in the presence of the alloy and ferrite nanoparticles arises when $\beta$ is enhanced; thus, the temperature field enhancement is evident, as shown in Figure 5.

Generally, thermal engineers are interested in the reduction of wall shear stresses, whilst they need to enhance the heat transfer rate of the system. For this purpose, the mathematician constructs different mathematical models that interpret the heat transfer and the fluid flow to examine which model better predicts the friction drag and Nusselt number. Regular fluids have low thermal properties; thus, they are not fast in transmission of thermal energy, whereas solids have higher transport properties, and hence their dispersion in base fluids make them more efficient for heat transfer. Therefore, C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni magnetic nanofluids were incorporated in this work. A higher resistance induced by the surface to the flowing fluid was noticed for the C₂H₆O₂-CoCr₂₀W₁₅Ni magnetic fluid, which indicates that the alloy particles can cause more resistance as compared to the ferrite nanoparticles, as evident in Figure 6. The higher resistance has a better impact on the corresponding velocity. This means that the corresponding velocity will be higher for the considered solid particle and the base fluid. Figure 7 demonstrates the Nusselt number for the C₂H₆O₂-Fe₃O₄, C₂H₆O₂-CoCr₂₀W₁₅Ni, H₂O-Fe₃O₄, and H₂O-CoCr₂₀W₁₅Ni magnetic nanofluids. The maximum transmission of heat is observed for the H₂O-Fe₃O₄ magnetic nanofluid. The thermal properties are higher for the base fluid as well as for the solid ferrite nanoparticles. Thus, the maximum heat transfer is noticed for the H₂O-Fe₃O₄ magnetic nanofluid, while the temperature behavior is lowered. A comparison of the present results was made with the available results examined by Ishak et al. [30]. The comparison is made for the Nusselt number as presented in

Table 3. Comparison of Nusselt number for various values of Pr (Prandtl number) in the flow of C₂H₆O₂, H₂O, CoCr₂₀W₁₅Ni, and Fe₃O₄.

Pr	Ishak et al. [30]	Present Results
0.01	0.0197	0.0155
0.72	0.8086	0.8088
1.0	1.0000	1.0000
3.0	1.9237	1.9298
7.0	3.0723	3.0765
10.0	3.7207	3.7200
100.0	5.2941	5.2965

.

5. Concluding Remarks

The phenomena of heat transfer analysis in ferromagnetic nanofluids on a stretching surface was explored with nanoparticles of magnetite ferrite and alloy in water and ethylene glycol (EG) as a carrier fluid. The main objective of the present theoretical study was to reduce the friction drag in the boundary layer flow and to intensify the efficiency of the considered nanofluid. The characteristics of nanoparticles together with the magnetic dipole impact enhanced the thermal conductivity of the engineered nanofluid and hence their influence on the heat flow. An appropriate transformation was employed on the mathematical model to convert the described system of partial differential equations into nonlinear ordinary differential equations and then computations were performed using the shooting method.

Combinations of nanoparticles of ferrite magnetite with a base fluid of water and ethylene glycol, and then of a cobalt–chromium–tungsten–nickel alloy in water and EG were examined. Moreover, the impact of various physical parameters involved in the flow model, such as ferrohydrodynamics interaction and others, were also analyzed on the velocity and temperature profiles to predict heat transfer rate. Finally, a comparison of different values of the Prandtl number indicated excellent agreement with previous data available in the literature. In conclusion, the nanofluid H₂O-Fe₃O₄ was found to exhibit a maximum heat transfer rate and C₂H₆O₂-CoCr₂₀W₁₅Ni with higher resistance as compared to the rest of the nanofluid combinations considered.