Forced convection of non-darcy flow of ethylene glycol conveying copper(II) oxide and titanium dioxide nanoparticles subject to lorentz force on wedges: Non-newtonian casson model.

The topic of two-dimensional steady laminar MHD boundary layer flow across a wedge with non-Newtonian hybrid nanoliquid (CuO-TiO2/C2H6O2) with viscous dissipation and radiation is taken into consideration. The controlling partial differential equations have been converted to non-linear higher-order ordinary differential equations using the appropriate similarity transformations. It is demonstrated that a number of thermo-physical characteristics govern the transmuted model. The issue is then mathematically resolved. When the method's accuracy is compared to results that have already been published, an excellent agreement is found. While the thermal distribution increases with an increase in Eckert number, radiation and porosity parameters, the velocity distribution decreases as porosity increases.

1 Introduction

A solid-liquid dispersion of 1–100 nm sized nanoparticles or nanofibers makes up nanofluids. Due to the distinct physical and chemical properties of nanometer-sized particles, nanofluids have a wide range of commercial uses. Because of the improvement in their thermal properties, nanofluids have gained a lot of attention. This area has been the subject of extensive research. Fluids with low thermal conductivity include blends of water, oil, and ethylene glycol. These liquids serve as a cooling technique that boosts productivity and lowers operational expenses. Ahmadi et al. (2018) reviewed theoretically and experimentally the influential parameters on the heat conductivity of numerous nanofluids. According to their results, greater temperatures and nanoparticle concentrations typically result in nanofluids with better thermal conductivities. Philip et al. (2016) examined nitrate salts doped with CuO nanoparticles for thermal energy storage with better heat transmission problems. Shaw et al. (2022) studied boundary layer flow caused by radial stretching in TiO2 nanotubes-water nanofluid. They observed that the momentum transmission of that nanofluid is controlled while using the non-Newtonian fluid model. Aladdin and Bachok (2020) reported the range of dual solutions expand widely for suction and closely reduce for solid volume fraction of Al2O3 and TiO2 nanoparticles. Nishant and Shriram (2016) found that the total heat transfer coefficient and thermal conductivity of water and ethylene glycol-based nanofluids by the use of CuO and TiO2 nanoparticles are enhanced.

A hybrid nanofluid is made up of a base-fluid and two or more nanoparticles. A new class of such mixtures known as hybrid nanofluids has just been created. Compared to nanofluids, this class of materials has more effective heat transfer characteristics. Eshgarf et al. (1959) provided an excellent review of hybrid nanofluid properties, preparation, models, and stability. They found that the hybrid nanofluids had much better thermal conductivity than the traditional nanofluids (single particles). Additionally, it was shown that the choice of base fluid affects the thermal conductivity of the hybrid nanofluids. Mishra et al. (2022) observed that the Cu-Al2O3/EG hybrid nanofluid flowing through a rotating channel exhibits a considerable increase in the shear rate as well as the rate of heat transfer. Khashi’ie et al. (2022) investigated the thermal progression of a second-grade hybrid Cu-Al2O3/CH3OH nanofluid towards a permeable surface. It has been shown by Alhowaity et al. (2022) that the velocity and heat conduction rate through a stretching surface are greatly increased by the dispersion of copper and graphene nanoparticulates into the base fluid ethylene glycol using a power-law non-Newtonian model.

There are numerous thermal engineering applications where fluid flows across wedge-shaped bodies happen, including geothermal systems, crude oil extraction, thermal insulation, heat exchangers, and the storage of nuclear waste, etc., Holden (1970) done an experimental as well as the theoretical review on boundary layer flow on wedge, flat plate and sphere in 1970. Both the joined and detached interaction zones over the models were thoroughly examined by them. He resulted that the separation first happened on the compression surface downwind of the flat plate and that it was not always possible to determine early separation by looking for the first time that a pressure distribution inflexion point appeared. Romanov and Kolesnikova (2021) reviewed the boundary value problems past a wedge and briefly discussed on the properties and applications. Xiu et al. (2022) evaluated the forced convection flow of water-based ternary hybrid nanofluid on wedge surfaces at various temperatures with a focus on the comparative examination of large and small volume of nanoparticles. It was found that forced convection flow on wedges has higher friction linked with higher velocity as time grows long, regardless of whether the volume of nano-particles is large or tiny. Ali and Abdul Alim (2022) investigated above the wedge surface and finalised that, particularly mixed convection, thermophoresis and heat generation parameters, have no effect on skin friction. Haq et al. (2022) researched and discussed that the increment in angle of rotation enhances the local Nusselt number.

Convection is the movement of fluid that allows heat to move from one location to another. Convective heat transfer combines the processes of conduction and advection, although it is frequently mentioned as a separate type of heat transfer. Convective heat transfer plays a significant role in processes involving high temperatures. Convective boundary conditions are more useful in material drying, transpiration cooling operations and other industrial and engineering processes. In a forced convection mechanism or kind of transport, fluid motion is produced by an outside source. Many researchers have participated in the investigations into the impact of forced convection as a result of the outstanding concept. Khurana et al. (2017) reviewed on TiO2, CuO, Al2O3 forced convection thermal transfer and resulted that the researchers employed less than 3% volume concentration in their tests because nanofluids can improve heat transmission at low particle concentrations. Cui et al. (2022), Khan et al. (2022), Yasir et al. (2022), Ahmad et al. (2022), Hoi et al. (2022), Zhang et al. (2022) have studied the role that the effect of forced convection plays on different fluid kinds and surfaces.

Radiation is the release of energy in the form of electromagnetic waves or ionising particles, especially high-energy particles. Chamkha (2000) looked at and explored how thermal radiation and buoyancy impacted hydromagnetic flow over a permeable surface that was accelerating and had a heat source or sink. It was discovered that the wall heat transfer was reduced as a result of thermal radiation, positive wall mass transfer, magnetic field, or heat creation. In the presence of Soret and Dufour’s effects, Chamkha and Ben-Nakhi (2008) studied the MHD mixed convection-radiation interaction along a permeable surface submerged in a porous liquid. They came to the conclusion that thermal radiation decreased the local Nusselt number and increased the local Sherwood number, particularly when suction was present. Heat and mass transfer analysis of steady/unsteady hybrid nanofluid (MWCNT-Ag/water) flow across a stretching sheet with thermal radiation was covered by Sreedevi et al. (2020). They discovered that for both steady-state and unstable instances, the thermal boundary layer thickness increases with cumulating values of radiation. Numerous studies on flow and heat transfer over various geometries in the presence of radiation for mono/hybrid nanofluids have been published [see, for instance, (Chamkha et al., 2003; Al-Mdallal et al., 2020; Khan et al., 2020; Algehyne et al., 2022; Hamad et al., 2022)].

Heat and mass transmission in porous media have been the focus of several investigations. Due to their widespread occurrence in industrial and technological applications, geothermal reservoirs, drying of porous solids, thermal insulation, increased oil recovery, packed-bed catalytic reactors, and many other applications are examples of certain applications. Along with the size of the potential gradient passing through the porous medium, drag forces encountered by the fluid inside the porous medium also affect how fluid flows through it. The two basic categories of drag forces that oppose fluid movement are surface drag (friction) and drag caused by solid barriers [see Rasool et al. (2022)]. The Darcy equation simply takes into account surface drag in its simplest version. Darcy’s equation is therefore only applicable for Reynolds numbers lower than one. Solid barrier drag is comparable to surface drag at higher velocities, or for Reynolds numbers greater than one. The nonlinear term known as the Forchheimer term describes the drag pressure drop. When high velocities are present, non-Darcy behaviour is crucial for explaining fluid flow in porous media. The free and forced convection nanofluid flow across a non-Darcian porous surface, which is based on the laminar boundary-layer concept, has been the subject of extensive research. Consider, for example, Chamkhax (1996), Chamkha (1997), Chamkha et al. (2006), Modather et al. (2009), Kandasamy et al. (2014), Qawasmeh et al. (2019).

When a magnetic field is applied to a fluid flow, an electromotive force is produced, altering the distribution of velocities. A magnetic field applied in the transverse direction of the flow, as opposed to one applied in the direction of the flow, directly affects the fluid’s velocity and may be more effective at controlling the flow. However, it requires more energy and increases drag forces [see Jha and Aina (2018)]. The interaction between the electrically conducting fluid and a magnetic field has an impact on several industrial equipment types, including boundary layer control, pumps, bearings, and MHD generators. Mittal et al. (2018) completed their review of the magnetohydrodynamic flow in nanofluids. They came to the conclusion that the magnetic field’s strength assisted the concentration rise as the velocity and temperature fell due to the Lorentz forces. The effects of Hall and ion slip on unstable MHD free convective rotating flow through a saturated porous media over an exponentially accelerating plate were discussed by Veera Krishna et al. (2020). It is reported that a magnetic field is said to prevent the stream from turning around. Several well-known authors (Chamkha, 1997), (Veera Krishna et al., 2018), (Ashraf et al., 2022), (Veera Krishna and Chamkha, 2020), (Veera Krishna and Chamkha, 2019), (Veera Krishna et al., 2021) addressed the MHD flow of electrically conducting liquids in various arrangements.

The viscosity of the fluid transforms some kinetic energy into thermal energy while fluid particles are in motion. Viscous dissipation is the term used to describe this irreversible process, which is induced by viscosity (Animasaun et al. (2022)). According to Koriko et al. (2021) an increase in the Eckert number implies an increase in the rate at which kinetic energy is being converted to internal energy. This is capable of boosting temperature distribution across the dynamics of blood-gold Carreau nanofluid and dusty fluid.

Any categorised fluid can attain density increases of two or three orders of magnitude, which are not insignificant in the polymer phase. Non-Newtonian fluids do not adhere to the law of viscosity. As a result, the idea of non-Newtonian dynamics was popularised by Isaac Newton. Most applications of non-Newtonian viscosity qualities are found in the field of biochemical engineering. According to Adewale et al. (2017), who tested the models for non-Newtonian liquids, the Casson physiochemical model accurately characterised in both low and high tensile conditions. Some other relevant studies are mentioned in (Lin and Lin, 1987; Chamkha and Ben-Nakhi, 2008; Iqbal et al., 2020a; Iqbal et al., 2020b; Kumar et al., 2020; Wakif et al., 2020; Alghamdi et al., 2021; Wakif et al., 2021; Zaydan et al., 2022).

Because of the importance of the above-mentioned factors, the main contribution of this research is an investigation of the non-Newtonian viscous dissipative mono nanoliquid TiO2/C2H6O2 and hybrid nanoliquid CuO - TiO2/C2H6O2 and subjected to an MHD forced convective flow past a non-Darcian wedge with radiation effects numerically. Therefore, this study is conducted to explore that what are the impacts of the non-Newtonian TiO2/C2H6O2 and CuO - TiO2/C2H6O2 nanoliquids boundary layer flow past the wedge with radiation and viscous dissipation effects?

2 Mathematical formulation

A steady 2-D MHD boundary layer flow of a Casson ethylene glycol-based hybrid CuO-TiO nanoliquid/TiO2nanoliquid over a wedge is studied at velocity u and reflects the atmospheric nanofluid velocity. The Cartesian coordinate system (x,y) is used, where the positive y-axis is taken normal to the wedge and the x-axis is picked along the wedge. The and represent the temperature of the wedge and the ambient fluid, respectively, as illustrated in Figure 1. The tensile state equation for Casson’s hybrid nanoliquid shape is as follows: [See (Xiu et al., 2022)] where is a component of deformation and is called the critical value, is the yield stress and is the plastic viscosity. The base liquid (Ethylene Glycol) is in thermal balance with the nanoparticles (cupric oxide CuO and titania TiO ) and there is no contact between the two solid components. Table 1 shows CuO, TiO and EG (nanoparticles and base liquid) thermophysical properties. It is assumed that the applied magnetic field is homogenous and that it is pointing in the direction normal to the surface. The maximum angle of the wedge is represented by , where is the Hartree pressure gradient. So, the set of governing equations of the hybrid nanoliquid movement are (Al-Mdallal et al., 2020), (Kandasamy et al., 2014) and (Qawasmeh et al., 2019)

FIGURE 1

Physical Configuration.

TABLE 1

Thermophysical properties of nanoparticles and base liquid Khan et al. (2022).

Properties	ρ(kg/m^2)	C_p(J/kgK)	k(W/mK)	σ(S/m)	Pr
CuO	6,510	540	18	5.96 ×10^7
TiO2	4,250	686.2	8.9538	2.38 ×10^6
EG	1,113.5	2,430	0.252	5.5 ×10^−6	204

Since the frequency of electron-atom collisions is predicted to be minimal, Hall and ion slip currents can be disregarded. In dimensional patterns, the initial and boundary requirements with forced convection can be expressed as. at

Here and are velocity components along the and axes, respectively. is the potential flow velocity, where , is a constant, Qawasmeh et al. (2019), is the Casson liquid parameter, and are kinematic viscosity, density, thermal diffusivity, dynamic viscosity, specific heat and magnetic field capacity respectively. The thermophysical properties of nanoliquid and hybrid nanoliquid is tabulated in Table 2. The subscript denotes hybrid nanoliquid and denotes nanoliquid. The subscripts and denotes the partial differentiation with respect to and respectively. are permeability, porosity and Forchheimer coefficient of porous media respectively. Whereas, is the drag coefficient.

TABLE 2

Thermophysical properties of hybrid nanoliquid and nanoliquid Al-Mdallal et al. (2020).

Properties	Hybrid nanoliquid (TiO2-CuO/EG)
μ_hnl	μ_l(1−φ_1)^2.5(1−φ_2)^2.5
ρ_hnl	{(1−φ_2)[ρ_s_1+(1−φ_1)ρ_l]}+φ_2ρ_s_2
ρC_hnl	{(1−φ_2)[φ_1ρC_s_1+(1−φ_1)]ρC_l]+φ_2ρC_s_2
α_hnl	k_hnl(ρC_p)_hnl
σ_hnlσ_nl	σ_s_2+2σ_nl−2φ_2(σ_nl−σ_s_2)σ_s_2+2σ_nl+φ_2(σ_nl−σ_s_2)
	Where σ_nlσ_l=σ_s_1+2σ_l−2φ_1(σ_l−σ_s_1)σ_s_1+2σ_l+φ_1(σ_l−σ_s_1)
k_hnlk_nl	k_s_2+(n−1)k_nl−(n−1)φ_2(k_nl−k_s_2)k_s_2+(n−1)k_nl+φ_2(k_nl−k_s_2)
	Where k_nlk_l=k_s_1+(n−1)k_l−(n−1)φ_1(k_l−k_s_1)k_s_1+(n−1)k_l+φ_1(k_l−k_s_1)
Properties	Nanoliquid (TiO2/EG)
μ_nl	μ_l(1−φ)^2.5
ρ_nl	(1−φ)ρ_l+φρ_s
ρC_nl	(1−φ)(ρC_p)_l+φ(ρC_p)_s
α_nl	k_nl(ρC_p)_nl
σ_nlσ_l	σ_s+2σ_l−2φ(σ_l−σ_s)σ_s+2σ_l+φ(σ_l−σ_s)
k_nlk_l	k_s+(n−1)k_l−(n−1)φ(k_l−k_s)k_s+(n−1)k_l+φ_1(k_l−k_s)

The radiative heat flux [see (Chamkha and Ben-Nakhi, 2008), (Sreedevi et al., 2020), (Hamad et al., 2022)] is now simplified using the Rosseland approximation for radiation:

The Stefen–Boltzmann constant and the mean absorption coefficient, respectively, are and . The temperature differences within the flow are assumed to be such that the term may be written as a linear function of temperature. As a result, extending in a Taylor series around while ignoring the higher-order variables yields:

Using transformation of similarity Eqs 3–5, 7 becomes with boundary conditions,

Initially,

In Eqs 8–10 the superscripts denote differentiation with respect to . Here is the first order porous resistance parameter, is the second order porous resistance parameter, is the magnetic field parameter, is the Prandtl number, is the Eckert number and is the radiation parameter.

Physical quantities skin friction coefficient which have now been described as

Using transformation of similarity Eqs 6, 10 becomes and local Nusselt number Nu is described as

Using Eq 6, Eq 12 becomes

3 Numerical simulation methodology

The numerical results for cupric oxide—titania/ethylene glycol and titania/ethylene glycol are obtained using MATLAB’s bvp4c function. The work findings are presented graphically, with a focus on the model’s mathematical key components and their effects on velocity, temperature, and engineering interest quantities. The first step is to make a first-order scheme of ODEs out of the third order ODEs (8) to (10). with limiting equations

The function bvp4c requires a first try at the answer because the tolerance for the current issue is set to 10–8. The assumption we made must adhere to the solution’s behaviour and satisfy the boundary conditions (26) and (27). Eqs 19–25 are numerically designed in order to get the original value argument to a conclusion. The bvp4c function in the MATLAB software is used to make both of these simplifications.

4 Results and discussion

It should be noted that in the current work, ethylene glycol (C2H6O2) is first mixed with nanosized cupric oxide (CuO) particles and then fixed to create a non-Newtonian nanoliquid. To create the hybrid nanoliquid TiO2 - CuO/C2H6O2, titania TiO2 nanoparticles are added to the nanoliquid and the thermophysical properties of these nanoparticles are mentioned in Table 1. Table 2 compares the findings of Qawasmeh et al. (2019), Chamkha, (1997), Chamkha and Ben-Nakhi (2008), for a given Prandtl number with m = 0, Ec = 0, , , R = 0 when , for Newtonian liquid. A fairly strong consensus can be seen from this table. The Prandtl number is fixed at 204 due to the base fluid being ethylene glycol. The physical parameter values are fixed as m = 0.0909, Ec = 0.01, , , R = 0.5 and throughout the study unless otherwise indicated in the graph.

Figure 2 depicts the velocity variation with magnetic parameter M. It is evident that as the magnetic parameter increases, the velocity increases as well for the TiO /C H O and TiO -CuO/C H O hybrid nanoliquids. This is due to the Lorentz force, which results in a retarded force on the velocity profile when a magnetic field is present. Moreover, as observed from the figure, the increment is similar for both TiO2/C H O and TiO2-CuO/C2H6O2.

FIGURE 2

Effects of velocity profiles for various .

Figure 3 depicts how the porosity parameter affects velocity curves. It is impressive that when rises, the nanofluid velocity drops on porous surfaces and its boundary layer widens. The velocity profile gets very close at increased porosity, which is seen in Figure 3 for TiO2/C2H6O2 and TiO2-CuO/C2H6O2 nanofluids similarly. This pattern is explained by the relatively low permeability and large porous inertia factors.

FIGURE 3

Effects of velocity profiles for various .

The effect of the magnetic parameter M on the thermal distribution is seen in Figure 4. The graph shows that when magnetic parameter values increase, temperature profiles decrease for the TiO -CuO/C H O hybrid nanoliquid and the TiO /C H O nanoliquid. The strong Lorentz force that causes the fluid flow to slow down as a result of the increased development of resistance temperature is caused up to after that lowering the temperature profile.

FIGURE 4

Effects of temperature profiles for various .

Figure 5 shows how the Eckert number Ec affects the temperature profile. The Eckert number illustrates the function-versus-viscous fluid tension conversion of external energy into internal energy, which causes the process to be irreversible. It has been established that a rise in the viscous dissipation parameter is accompanied by a rise in temperature. From Figure 5. It is noted that for higher values of the Eckert number, the thermal distribution is a little higher for TiO2-CuO/C2H6O2 nanoliquid than for TiO /C2H6O2 nanoliquid.

FIGURE 5

Effects of temperature profiles for various .

On the other hand, Figure 6. Reveals that the increasing values of radiation parameter (R) result in the increasing temperature distribution of the TiO -CuO/C2H6O2 hybrid nanoliquid and the TiO /C2H6O2 nanoliquid flow region. This is because there would be an increase in the thermal boundary layer thickness with the increase of radiation parameter R.

FIGURE 6

Effects of temperature profiles for various .

Figure 7 illustrates how the porosity affects the non-dimensional temperature profile. The temperature initially drops as the porosity increases, then it rises. The thickness of the thermal boundary layer also increases with porosity. This behaviour is the same for the mono nanofluid and hybrid nanofluid.

FIGURE 7

Effects of temperature profiles for various .

Figures 8–10 show hybrid nanoliquid and nanoliquid skin friction past a wedge. It is revealed from Figures 8–10 that the increasing values of porosity and Casson parameter result in a decreasing skin friction coefficient. In the meantime, the increasing values of m and M result in increasing the skin friction coefficient. Figures 8–10. Illustrate the common trend that the TiO -CuO/C2H6O2 hybrid nanoliquid has heavier skin friction values than the TiO2/C2H6O2 mono nanoliquid. This is due to the fact that the friction made by hybrid nanoparticles TiO and CuO is greater than mono TiO nanoparticles.

FIGURE 8

Effects of skin friction for certain and .

FIGURE 10

Effects of skin friction for certain and .

Figures 11, 12 illustrate the TiO -CuO/EG hybrid nanoliquid and the TiO2/EG nanoliquid Nusselt number over the wedge for the increment of R, and Ec. It is clearly observed that the increasing values of R, and m escalate Nusselt number for the TiO /EG nanoliquid than the TiO -CuO/EG hybrid nanoliquid. Simultaneously, the increasing values of Ec drop the Nusselt number. While relating Figures 11, 12, it is revealed that the TiO -CuO/EG hybrid nanoliquid has a lesser Nusselt number than the TiO2/EG mono nanoliquid.

FIGURE 11

Effects of Nusselt number for certain and .

FIGURE 12

Effects of Nusselt number for certain and .

5 Conclusion

The momentum and thermal behaviour of viscous dissipative mono and hybrid nanofluids and ) past a wedge embedded in the non-Darcy porous medium with radiation and magnetic field are examined numerically. From the computational results, the following conclusions are given (Table 3).

TABLE 3

Amount of for certain Prandtl number with for Newtonian liquid when .

Pr	Current study	Qawasmeh et al. (2019)	Chamkha (1997b)	Alghamdi et al. (2021)
10^0	0.33226941	0.332054	0.332173	0.332058
10^1	0.72830834	0.728136	0.72831	0.728148
10^2	1.57217714	1.571821	1.57218	1.57186
10^3	3.38788227	3.387073	3.38809	3.38710
10^4	7.29934030	7.297260	7.30080	7.29742

1) The velocity increases when M is increased and decreases when is increased for both the nanofluid and hybrid nanofluid.

2) For and , temperature rises when and Ec are enhanced and declines when M is raised.

3) An increase in the values of and decreases and the raising values of M increases the skin friction respectively.

4) A local Nusselt number is an increasing function for the increasing values of R and

5) Both the skin friction and Nusselt number increase with the effect of m.

6) Greater values of Ec lower the Nusselt number for both and .

7) While comparing the and nanoliquids, the skin friction is a little greater for and the Nusselt number is a little higher for nanoliquid.