Effect of Thermal Radiation on Three-Dimensional Magnetized Rotating Flow of a Hybrid Nanofluid.

The effect of thermal radiation on the three-dimensional magnetized rotating flow of a hybrid nanofluid has been numerically investigated. Enhancing heat transmission is a contemporary engineering challenge in a range of sectors, including heat exchangers, electronics, chemical and biological reactors, and medical detectors. The main goal of the current study is to investigate the effect of magnetic parameter, solid volume fraction of copper, Eckert number, and radiation parameter on velocity and temperature distributions, and the consequence of solid volume fraction on declined skin friction and heat transfer against suction and a stretching/shrinking surface. A hybrid nanofluid is a contemporary type of nanofluid that is used to increase heat transfer performance. A linear similarity variable is-applied to convert the governing partial differential equations (PDEs) into corresponding ordinary differential equations (ODEs). Using the three-stage Labatto III-A method included in the MATLAB software's bvp4c solver, the ODE system is solved numerically. In certain ranges of involved parameters, two solutions are received. The temperature profile θη upsurges in both solutions with growing values of EC and Rd. Moreover, the conclusion is that solution duality exists when the suction parameter S≥Sci, while no flow of fluid is possible when S<Sci. Finally, stability analysis has been performed and it has been found that only the first solution is the stable one between both solutions.

1. Introduction

Fluid dynamics study has attracted the interest of experts, scholars, and researchers from numerous fields in current times owing to its multiple applications in engineering, science, and technology, as well as biopolymer-based detectors for medical sensing applications and diagnosis. Sakiadis [1] initially suggested the notion of a boundary layer steady flow on a stretching two-dimensional surface. Later, Crane [2] modernized the thought of Sakiadis. Then he applied it on steady-flow two-dimensional linear exponentially stretching surfaces. He claimed that the speed with which a sheet is stretched from a slit is related to the distance between them. Due to the fact that nanoparticles can pass through tissues and cells, chemical researchers state that nanotechnology can offer innovative solutions for treating stenosis, which is harmful and can cause death. It is noteworthy to mention the advanced evolvement of nanoparticles in drugs [3,4,5,6,7]. Choi [8] initiated research with a focus on nanoparticles by revealing their dynamic and abnormal attributes. Furthermore, he identified a nanofluid as a mixture generated by spreading nanoparticles in a fluid. He defined nanomaterials as heat transfer dispersions with better thermal performance than base fluids or conventional liquids. Nanofluids are composed of tiny quantities of solid particles measuring 100 nm or less in size. Nanofluids can be developed by evenly spread nanoparticles in a conventional liquid. Nanofluids often incorporate nanoparticles with unique chemical and physical properties, for instance, carbides, oxides, carbon nanotubes, and metals [9]. Water, organic liquids, polymeric solutions, oils, and other common fluids are examples of base fluids [10]. It has been discovered that nanofluids have better thermal efficiency than their basic fluids. Consequently, nanofluids have abundant real-world applications, including heat transfer. Agriculture, refrigerators, aerospace, and automobiles are just a few examples of their applications [11]. Khan and Pop [12] examined numerically the stretched surface in a nanofluid with a steady two-dimensional flow of the boundary layer. Miklavčič and Wang [13] are credited with being the first to explore the viscous steady three-dimensional flow across a shrinking surface using suction. Bachok et al. [14] examined a nanofluid with an unsteady two-dimensional boundary layer flow on a stretched/shrunk sheet. Moreover, Deroet al. [15] considered the two-dimensional model of a single phase of Tiwari–Das steady flow over a stretched and shrinking sheet during the analysis of a Casson-based nanofluid.

Researchers have attempted to incorporate numerous solid nanoparticles through various types of base fluids in response to the increasing requirements for heat transfer in several branches of the industry, as previously mentioned. As a result, a hybrid nanofluid, the newest type of nanofluid, was developed. According to Devi and Devi [16], combining a metal nanoparticle in a small quantity/nanotube with an oxide/suspension of metal nanoparticles in a base fluid will dramatically boost thermal properties. Common fluids usually have more than one solid particle, such as glycol, paraffin oil, water, vegetable oil, engine oil, kerosene, and ethylene. According to Huminic and Huminic [17], there are several significant applications of hybrid nanofluids. For example, these liquids have been used for heat transfer in exchangers of plate heat, micro channels, air-conditioning systems, etc. Radiation’s effects on the steady flow of a two-dimensional hybrid nanofluid on a nonlinear stretched/shrunk surface were investigated by Waini et al. [18]. Moreover, Waini et al. [19] studied the influence of transpiration on the steady flow of a two-dimensional hybrid nanofluid and heat transfer through a stretched/shrunk sheet with a uniform shear flow.

Many fields use rotational flow, such as engineering processes (e.g., rotating machinery and computer storage devices), electronic devices, extrusion of plastic sheets, and the centrifugal filtration process. Hayat et al. [20] investigated the /water three-dimensional hybrid nanofluid’s steady flow in the presence of radiation, chemical reactions, and heat generation. It was discovered that a hybrid fluid has a superior rate of energy transfer compared toa regular nanofluid. When the rotation parameter is raised, the concentration profile increases. Furthermore, Hayat et al. [21] scrutinized a three-dimensional rotating hybrid nanofluid flow using boundary influence with radiation and partial slip. They came to the conclusion that the temperature profile improved as the rotation and radiation parameters were increased. Later, Anuaret al. [22] studied the hybrid nanofluid steady rotating flow for a stretched/shrunk surface under the influence of radiation. On hybrid nanofluids, numerous scholars have investigated the stretching/shrinking sheet and the rotating sheet in models with different considerations [23,24,25,26].

Magnetohydrodynamics (MHD) is a term that combines the words magneto (magnetic), hydro (fluids), and dynamics (motion). Magnetohydrodynamics is the study of the flow of electrically conducting fluids in the presence of a magnetic field. Since a magnetic field can occur anywhere in the universe, MHD phenomena can also occur when conducting fluids are present in natural phenomena. Engineers use MHD concepts in the project and design of an extensive range of applications in industries (e.g., heat exchanger generators and MHD pumps) [27]. So it is vital. Devi and Devi [16] numerically explored the two-dimensional magnetic –time-independent flow through a stretched sheet with suction and the Newtonian heating effect in a hybrid nanofluid. Furthermore, Devi and Devi [28], with Newtonian heating, expanded work to three-dimensional time-independent. They reported that a hybrid nanofluid transfers heat more quickly than a normal nanofluid. After that, Devi and Devi [29] examined the two-dimensional steady flow of a hybrid nanofluid over a stretched sheet. In a hybrid nanofluid, Aly and Pop [30] studied the two-dimensional steady flow over a stretched and shrunk plate using suction, the biot number, and MHD effects. Wainiet al. [31] explored the influence of hybrid nanoparticles on a fluid’s steady flow in an exponentially stretching/shrinking two-dimension allayer. It was found that suction effects and shrinking effects had dual solutions. Lund et al. [32] analyzed two-dimensional unsteady flow of a hybrid nanofluid on stretching and shrinking sheets under the influence of radiation and MHD. Teh and Asghar [33] studied the MHD hybrid nanofluid flow under the influence of Joule heating on a three-dimensional rotating stretching/shrinking surface. The act of creating heat by passing current through a conductor is known Joule heating, also known as ohmic heating, resistance heating, or resistive heating. Incandescent lightbulbs, resistance ovens, electric stoves, soldering irons, and cartridge heaters are all examples of Joule heating. Additionally, Khashi et al. [34] assessed the two-dimensional steady flow behavior of a hybrid /water nanofluid associated with a radially stretched/shrunk sheet under the impact of Joule heating, MHD, and suction. Moreover, Yan et al. [35] explored the impact of Joule heating on electrically conducting hybrid nanofluid flow moving over an exponential surface.

Thermal radiation has a significant influence at high operating temperatures that cannot be overlooked. Radiation sensitivity is crucial in the design of a suitable method since several industrial processes take place at extremely high temperatures. It also plays a substantial role in a variety of industrial applications, such as glass processing, furnace construction, internal combustion engines, plasma physics, and spacecraft [36].Sreedevi et al. [37] investigated the two-dimensional unsteady flow of mass and heat transfer of a hybrid nanofluid on a stretched surface by radiation. Recently, Dero et al. [38] studied the numerical analysis of the water hybrid nanofluid in terms of cross-flow and stream wise under the influence of thermal radiation. Some study has been done on flow and thermal radiation. Some current publications can be found in these studies [39,40,41].

In the two-dimensional case, research on hybrid nanofluids under MHD, radiation, and Joule heating has been rigorous (see Waini et al. [18], Khashi’ie et al. [34], and Yan et al. [35]). According to the literature survey on previous studies, it is critical to address such a problem as three-dimensional flow since this setting enables more extensive and realistic real-world applications. The present study expanded on the work of Anuar et al. [22] by incorporating the influence of radiation, MHD, and Joule heating using the research of Tiwari and Das [42]. Therefore, for this study, a novel physical model of the influence of thermal radiation on the three-dimensional magnetized rotating flow of a hybrid nanofluid with Joule heating is developed. To reach high convective heat transfer efficiency, in the hybrid nanofluid measured in this study, alumina () and copper () nanoparticles are used. By suspending two nanoparticles, a hybrid nanofluid (/water) is designed. Profiles of temperature and velocity for several values of volume fraction copper, magnetic, Eckert number, and radiation parameters of solid nanoparticles are examined in this study. Furthermore, the effect of solid nanoparticle copper against stretched/shrunk and suction parameters in terms of declined heat transfer and change in skin friction are also included in this study. The present numerical findings are compared to the results of prior investigations for comparison purposes. To the best of our knowledge, this model is different and new and no related article has been found in the literature.

2. Materials and Methods

2.1. Mathematical Formulation

Three-dimensional hybrid nanofluid steady flow with heat transfer in stretched and shrunk sheets is shown in Figure 1, which has Cartesian coordinates with the x-axis and the y-axis measured in the plane. The fluid that occupies half of the area is at . In the x-direction, the surface is deemed to be stretched/shrunk only when velocity and temperature is [43]. The entire system rotates at a constant velocity along a perpendicular to the surface z-axis. In the z-axis, there is a variable magnetic field and radiation . The surface of the hybrid nanofluid is constant and has ambient temperature . is the wall temperature, and is the characteristic temperature. The mathematical governing equations help understand the physical model by indicating that the flow is incompressible, Newtonian, and laminar. Figure 1 shows that fluid stages and nanoparticles are considered in a thermal equilibrium state. They are distinctly tiny and uniform in shape, allowing the slip velocity among phases to be discounted.

Figure 1

The physical model of a stretching/shrinking surface.

The governing equation for a hybrid nanofluid is as follows [22]:

The boundary conditions are

The hybrid nanofluid velocity components and are along the x-axis, the y-axis, and the z-axis, respectively. is the stretched/shrunk parameter such that shows the shrinking surface, represents the static surface, and indicates the stretching sheet. represents temperature. Hybrid nanofluid thermophysical properties were used to determine the above equation, as mentioned in Table 1 and Table 2. Alumina () nanoparticles are denoted by , and copper () nanoparticles are indicated by . It should be noted that the current model is used for only spherical nanoparticles because it gives better heat transfer performance. Moreover, , and correspond with specific heat capacity, density, dynamic viscosity, electrical conductivity, and thermal conductivity, respectively. Fluid, nanofluid, hybrid nanofluid, solid nanoparticles ( and solid nanoparticles () are denoted by subscripts and respectively.

Table 1

Thermophysical properties of a hybrid nanofluid [33].

Names	Properties
Dynamic viscosity	μ_hnf=μ_f1−ϕ_1^521−ϕ_2^52
Density	ρ_hnf=1−ϕ_21−ϕ_1ρ_f+ϕ_1ρ_s1+ϕ_2ρ_s2
Thermal conductivity	k_hnf=k_s2+2k_nf−2ϕ_2k_nf−k_s2k_s2+2k_nf+ϕ_2k_nf−k_s2×k_nf where k_nf=k_s1+2k_f−2ϕ_1k_f−k_s1k_s1+2k_f+ϕ_1k_f−k_s1×k_f
Heat capacity	ρc_p_hnf=1−ϕ_21−ϕ_1ρc_p_f+ϕ_1ρc_p_s1+ϕ_2ρc_p_s2
Electrical conductivity	σ_hnf=σ_2+2σ_nf−2ϕ_2σ_nf−σ_2σ_2+2σ_nf+ϕ_2σ_nf−σ_2×σ_nf where σ_nf=σ_1+2σ_f−2ϕ_1σ_f−σ_1σ_1+2σ_f+ϕ_1σ_f−σ_1×σ_f

Table 2

Solid nanoparticles and base fluid(water) thermophysical properties.

Fluids	ρ (kg/m3)	σ Sm	c_p (J/kg K)	k (W/m K)	Pr
Copper (Cu)	8933	5.96 × 10^7	385	400	-
Alumina (Al_2O_3)	3970	3.69 × 10^7	765	40	-
Water (H_2O)	997.1	0.05	4179	0.613	6.2

The similarity factors indicated below are employed in this example [44].

Here, prime shows the differentiation for , is a stretching constant, and is the kinematic viscosity. Although , stands for the injection/suction parameter. When the flow is injection, and when the , the flow is suction.

Equation (1) is identically satisfied by inserting Equation (6) in Equations (2)–(5). Subsequently, Equations (2)–(5) are transformed into the following ODEs.

Subject to conditions: where is the rotation, the magnetic parameters, is the Eckert number, is the radiative heat flux, and is Prandtl number, which is described by:

The skin friction coefficients and show the x-axis and the y-axis, respectively. The local Nusselt number is specified as:

By using Equations (6) and (12), we get:

Here, is the local Reynolds number.

The physical properties of the solid nanoparticles and the base fluid(water) are presented in Table 2 [27,33].

2.2. Stability Analysis

The system of Equations (7)–(9) indicates that more than one numerical solution exists in a specific array of stretched/shrunk and suction parameters. Therefore, an analysis of stability is needed to determine which of the solutions is stable [45,46]. That being stated, the transformation requires the introduction of new similarity variables, which are defined as where represents the dimensionless time (). When Equation (14) is substituted into the unsteady Equations (2)–(4), the following equations can be acquired.

With the new boundary condition

We use perturbed equations to verify the steady flow stability of , and , which satisfy Equations (7)–(9). where , and are modest relative to and , respectively, and is the eigenvalue parameter. The following linearized eigenvalue issue is generated by inserting Equation (19) into Equations (15)–(17), fixing and simplifying the equations. subject to boundary conditions

The boundary conditions stated in Equation (23) must be changed prior to resolving the problem of eigenvalue. Harris et al. [47] proposed relaxing and replacing one of the far-field conditions. So, in this research, relaxing is impossible since it will be substituted by , which can influence the other numerical solutions and boundary conditions. As a result, we selected to be loosened and substituted with . This allows the stability of the solution to be determined, where indicates that the solution is stable, whereas shows that the solution is unstable owing to the development of the disturbance of the solution.

3. Results and Discussion

The solution’s duality has been achieved in the figures by using various initial guesses for , , and , with the outcome that both velocities and temperature profiles satisfied the boundary condition asymptotically. By using the software MATLAB bv4pc solver, we numerically solved nonlinear ODEs (Equations (7)–(9)) with the boundary equation (Equation (10)). The bvp4c solver was created by Jacek Kierzenka and Lawrence F. Shampine of Southern Methodist University, Texas [48]. The bvp4c solver is a finite difference method that uses the three-stage Lobatto IIIA algorithm to provide numerical solutions with fourth-order precision. Before starting to discuss the results of the current study, we have compared the coding of a numerical method to make sure that our computer code is working properly. First, to validate the coding of a numerical scheme in this study, the reduced skin friction and for pure waterwhen ,, and is compared with values issued by Nazar et al. [49] and Anuaret al. [22], as in Table 3, for stretching sheet = 1. The findings are consistent with those of previous studies. To yield the wanted hybrid nanofluid, i.e., /water, first, we show that alumina is isolated in the base fluid (water). After that, copper () is dispersed into the nanofluid /water. Here, represents alumina ( and represents copper ().

Table 3

The comparison results of rotation parameter at different values when,,, and .

	Nazar et al. [49]		Anuar et al. [22]		Current Results
ω	f^″0	g^′0	f^″0	g^′0	f^″0	g^′0
0.5	−1.1384	−0.5128	−1.13838	−0.51276	−1.138374	−0.512760
1.0	−1.3250	−0.8371	−1.35503	−0.83710	−1.325028	−0.837098
2.0	−1.6523	−1.2873	−1.65235	−1.28726	−1.65235	−1.287258

The suction parameter ( ) of various values with the variation in , , and with three values of solid nanoparticle volume fraction copper ispresented in Figure 2a–c for (shrinking sheet) with the occurrence of different parameter values and . The range of solid nanoparticles is 0 to 0.06. Here, shows the critical point at the suction parameter where both solutions meet each other. It can be shown that when , no solution exists. Furthermore, the values of = 1.9868, 1.8810, and 1.7681 are relative critical points of and 0.06, respectively. It is worth revealing here that when , it is purely water-based nanofluid and . After that, 2% of is added and is achieved. Moreover, the value of appears to rise as 6% of the solid volume fraction of is added to the hybrid nanofluid.The boundary layer separation extends as the value of increases. From Figure 2a, in the first solution, when values of increase, the reduced skin friction increases and on the other side, the values of decrease in the second solution. In Figure 2b, the first solution declines and the second solution of is improved when the values of increase. As can be seen in Figure 2c, when the value of rises, the magnitude of is also decreased in both first and second solutions.

Figure 2

(a–c): Influence of against suction on , , and .

Due to the existence of multiple solutions, we investigated which of the governing parameters can play a role in the formation of dual solutions. Against varying values of solid nanoparticles, Figure 3a–c demonstrates volume fraction along with reduced skin friction and and reduced heat transfer rate through the presence of parameters and . For shrinking sheet , we can understand that the presence of non-unique solutions is notable. When , represents the critical point at which both solutions meet at the shrinking sheet parameter. It can be shown that no solution exists when but a unique solution is observed when . Beyond these critical values, the boundary layer estimate is no longer defensible. Furthermore, the values of ,, and are the corresponding critical points of and 0.06, respectively. This is proof that increasing the value of can cause the boundary layer separation to be delayed. Figure 3a shows that in the first solution, the value of increases when and decreases when , while the second solution declines after the increment of . In Figure 3b, when the influence of improves, the value of decreases in the first and increases in the second solution. In Figure 3c, declines in both solutions when increases. Similar findings can be seen in Anuar et al. [22,26].

Figure 3

(a–c): Influence of against stretching/shrinking on , , and .

Figure 4a–c depicts different values of = 0, 1, and 3 when the velocity and temperature profiles are and respectively, under different parameters, such as and . Figure 4a illustrates that profiles decrease with an upsurged value of in both (first and second) solutions due to the high drag force created by the magnetic number. The effect of on profiles is shown in Figure 4b. It is observed that profiles reduce with an upsurged value of in both solutions. Figure 4c shows that increases as the value of increases in both solutions. As a consequence, the influence of heat transfer is increased as increases. The Lorentz force, which is caused by the magnetic field, makes the rate of transfer less resistant [49].

Figure 4

(a–c) Profiles of velocity , , and temperature for different .

The plots of the velocity profile and the temperature profile against solid nanoparticle volume fractions and of copper are portrayed in Figure 5a–c, respectively. Figure 5a shows that in thefirst and second solutions, the fluid velocity decreases when the value of is enhanced. From Figure 5b, the profile in both solutions declines with a growing value of . Figure 5c show that in the first and second solutions, the temperature profile increases when the value of is enhanced.

Figure 5

(a–c): Profiles of velocity , , and temperature for different .

Figure 6 illustrates the value of at a different point for the temperature profile . In the first and second solutions, the values of increase, causing to rise as well. The intensity of heat transfer rises as the value of grows due to the increasing heat created by Joule heating. It is worth noting that the influence of the increasing values of on is not important because, as seen in Table 4, the values remain the same for the increasing values of . Some variations in the value of are provided in Table 4.

Figure 6

Profile of for different parameters.

Table 4

The results of the Eckert number when , ,, , =, and .

Ec	−θ^′0
0	8.247357 (8.476272)
0.3	8.247162 (8.476187)
0.6	8.246967 (8.476103)

Note: The values inside the brackets represent the values of the second solution.

Figure 7 shows the various values of radiation parameter and for the temperature profile . The radiation parameter only exists in Equation (9). As a result, the values of have no influence on the velocity profile since they are uncoupled from the momentum Equations (7) and (8).The thermal boundary layer thickness constantly rises in the first and second solutions with a growing value of , which means that higher values result in a reduced temperature gradient at the surface.Due to the existence of high radiation, a big quantity of heat energy is produced in the system, which causes the temperature to increase and implies that the temperature of the fluid rises.

Figure 7

Profile of for different parameters.

Because dual solutions are obvious in the data, we ran a stability analysis to assess the solution’s practicability. Based on our findings, we determined that only the first option is stable and viable for use as guidance in real-world applications. This claim is supported by tabulation (see Table 5), which shows that the solution is stable when the executed minimum eigenvalues are positive and non-stable otherwise.

Table 5

The results of the suction parameter at different values of when

		ε_1
ϕ_2	S	First Solution	Second Solution
0	3	0.2042	−0.2233
	2.5	0.1584	−0.2093
	2	0.0806	−0.1003
	1.9868	0.0071	−0.0251
0.02	3	0.8668	−0.8978
	2.8	0.7704	−0.7929
	2.6	0.5954	−0.5988
	2.4	0.3723	−0.4012
	2.2	0.1912	−0.0853
	1.8810	0.0108	−0.0071
0.06	3	1.1380	−1.0836
	2.8	1.0230	−0.9651
	2.6	0.8186	−0.7554
	2.4	0.6591	−0.5691
	2.2	0.3532	−0.2856
	2	0.0993	−0.1003
	1.7681	0.0080	−0.0006

4. Conclusions

In this study, a three-dimensional /water MHD hybrid nanofluid with heat transfer rotating flow on a linear stretched/shrunk surface underthe effect of radiation and Joule heating has been investigated through a bvp4c solver on the MATLAB computing platform. This research concentrates on the behavior of , and underthe effect of suction/injection, MHD, radiation, and Joule heating on the above-stated hybrid nanofluid flow. The following are the key conclusions of the present investigation:

The outcomes of Equations (2)–(5) are non-unique when and unique when . No solution exists when is noticed for a shrinking sheet.

The presence of solutions is based on the values of the suction parameter for assumed parameter values.

The intended model displays good heat transfer performance when nanoparticle sizes are less than 10%.

The temperature and velocity profiles of the boundary layer can be changed by applying nanoparticles in the hybrid nanofluid.

The profile is reduced when the magnetic parameter is enhanced.

and values are directly proportional to temperature profiles in both solutions.

The first solution is a stable solution.