# Numerical Simulation of Mixed Convection Squeezing Flow of a Hybrid Nanofluid Containing Magnetized Ferroparticles in 50%:50% of Ethylene Glycol–Water Mixture Base Fluids Between Two Disks With the Presence of a Non-linear Thermal Radiation Heat Flux

^{1}Department of Mathematics, College of Arts and Science, Prince Sattam Bin Abdulaziz University, Wadi Al Dawasir, Saudi Arabia^{2}Department of Mathematics and Social Sciences, Sukkur IBA University, Sukkur, Pakistan^{3}Department of Natural Sciences, The Begum Nusrat Bhutto Women University, Sukkur, Pakistan^{4}Department of Mathematical Sciences, Federal Urdu University of Arts, Science & Technology, Karachi, Pakistan^{5}Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City, Vietnam^{6}Department of Mathematics, Cankaya University, Ankara, Turkey^{7}Institute of Space Sciences, Măgurele, Romania^{8}Department of Medical Research, China Medical University Hospital, China Medical University, Taichung City, Taiwan

Ferroliquids are an example of a colloidal suspension of magnetic nanomaterials and regular liquids. These fluids have numerous applications in medical science such as cell separation, targeting of drugs, magnetic resonance imaging, etc. The hybrid nanofluid is composed by scattering the magnetic nanomaterial of more than one type nanoparticles suspended into the base fluid. It has different scientific applications such as heat dissipation, dynamic sealing, damping, etc. Owing to the vast ferrofluid applications, the time-dependent squeezed flow of hybrid ferroliquids under the impact of non-linear radiation and mixed convection within two disks was explored for the first time in this analysis. Here, the cobalt and magnetite ferrofluids are considered and scattered in a 50%:50% mixture of water–EG (ethylene glycol). The similarity technique is used to reduce the leading PDEs into coupled non-linear ODEs. The transmuted equations together with recommended boundary restrictions are numerically solved via Matlab solver bvp4c. The opposing and assisting flows are considered. The impacts of an emerging parameter on fluid velocity and temperature field of hybrid ferroliquids are examined through the different graphical aids. The results showed that the opposite trend is scrutinized due to the magnetic influence on the temperature and velocity in the case of assisting and opposing flows. The velocity augments due to the volume fraction of nanoparticles in the assisting flow and declines in the opposing flow, while the opposite direction is noticed in the temperature field.

## Introduction

The time-dependent squeezed flow is regularly addressed in realistic surroundings to describe the movement of fluid along an area of contracting of a recommended length. This flow can be set up via placing a preliminary sluggish liquid between a couple of rectangular parallel disks and plates to indicate the related mathematical models on specific planes of coordinates. The research interest in studying the squeezed flows is quickly developed owing to existing and rising applications in polymer processes of manufacturing and biological liquid transport, purification, lubrication, molding injection, filtration, etc. Ishizawa (1966) inspected the motion of squeezed flow by unstable distance in moving two disks and employed the perturbation technique to find the solution. An exact solution of a similar problem in a couple of elliptic plates with the multifold time-dependent series was obtained by Usha and Sridharan (1996). Duwairi et al. (2004) explored the time-dependent viscous squeezing flow between parallel two plates with a constant temperature. Domairry and Aziz (2009) applied the magnetic function on squeezing flow between two disks. The magnetic function on squeezing flow connecting two analogous permeable disks was inspected by Hayat et al. (2012a). Since then, many researchers have contributed to the progress of squeezed flow for distinct models of fluid (Hayat et al., 2012b; Qayyum et al., 2012; Sheikholeslami and Ganji, 2013; Sheikholeslami et al., 2013; Domairry and Hatami, 2014). In recent times, Hayat et al. (2018a,b) discussed the squeezed flow comprising carbon nanotubes through two parallel disks in a Darcy–Forchheimer medium with thermal radiation. Alzahrani et al. (2020) analyzed the squeezed 3D flow through a rotated in a channel with Dufour and Soret impacts. The impact of imposed magnetic forces on squeezed flow containing nanomaterials between two plates was examined by Babazadeh et al. (2020).

The exploration of nanofluids has acknowledged a prominent interest owing to their incredible range of applications involving the processing of nanomaterials, the coolant of automobiles, polymer coating, tribology of aerospace, and suspensions of medical sterilization. Choi (1995) coined the term “nanofluid” for the first time, which refers to the submicron scattering of solid material into a liquid, having greater thermal conductivity in a convectional liquid. It is perceptible that these nanomaterials are considered ultrafine (1–50 nm order of length); hence, the nanoliquids appeared to conduct further as a single-phase liquid than a suspension of solid–liquid. Comprehensive monographs regarding the nanofluids in conjunction with several applications can be seen in Tiwari and Das (2007), Chamkha et al. (2012), El-Kabeir et al. (2015), Wakif et al. (2018), Sheikholeslami (2019), Soomro et al. (2019), Mebarek-Oudina et al. (2020), and Raza et al. (2020). More precisely, the magnetic nanoliquid (ferroliquids), which is colloidally deferred of magnetic nanomaterials like ferrite, cobalt, and magnetite, is suspended in non-conducting fluids like heptanes, water, hydrocarbons, and kerosene. These ferrofluids have numerous applications in medical science like cell separation, targeting of drug medicine, imaging of magnetic resonance, etc. Neuringer (1996) inspected the thermal gradient and magnetic function on the flow of ferrofluids. The impact of uniform heat flux on slip flow with heat transport from a flat surface was analyzed by Khan et al. (2015). They considered three different ferrofluids (CoFe_{2}O_{4}, Fe_{3}O_{4}, and Mn-ZnFe_{2}O_{4}) with two different base fluids (water and kerosene). Rashad (2017) scrutinized the magnetic function on slip flow containing kerosene-based cobalt ferrofluid from a non-isothermal wedge with mixed convection and radiation effect. Zaib et al. (2019) inspected the entropy on mixed convective flow comprising the magnetite ferroliquid of micropolar fluid from a vertical plate. Recently, Ali et al. (2020) discussed the micropolar fluid with magnetic dipole impact comprising EG and water-based ferrofluids, namely, Fe and Fe_{3}O_{4} from a stretched sheet. On the other hand, hybrid nanoliquids are deliberate to catch superior properties of thermal as well as rheological via the mixing couple of distinct nanoparticles. A reason following the initiation of hybrid ferroliquids is to knob the heat transport efficiency in the fluid flow. It has different scientific applications as dynamic sealing, naval, damping, acoustics, microfluidics, etc. The impact of dissipation on time-dependent flow comprising hybrid nanoliquid through a rounded pipe was examined by Suresh et al. (2012). They accomplished that nanoliquid has a lower friction factor than hybrid nanoliquid. Madhesh and Kalaiselvam (2015) discussed the pressure drop, in which the Reynolds and Peclet numbers are investigated for the volume fraction of water-based Cu-TiO_{2} hybrid nanoliquid. Adriana (2017) presented the correlation of alumina hybrids and nanofluids through the date of the temperature gradient. The features of fluid flow of water-based CuO-Ag hybrid nanoliquids were inspected by Hayat and Nadeem (2017) under the impacts of chemical reaction and radiation. Mebarek-Oudina (2019) discussed the characteristics of thermal and hydrodynamic of Titania nanoliquids satisfying a cylindrical annulus. Mahanthesh et al. (2019) analyzed the influences of an exponential space-dependent heat source on magneto slip flow comprising carbon nanofluids from a stretchable rotating disk. Wakif et al. (2020) explored the influences of surface roughness and radiation containing the stability of Al_{2}O_{3}-CuO hybrid nanoliquid by employing the generalized Buongiorno model. Recently, Marzougui et al. (2020) investigated the entropy analysis through convective flow involving nanoliquid via MHD with chamfers in a cavity.

The scrutiny of the magnetic function has considerable solicitations in MHD bearings; geology, astrophysics, pumps, generators, medicine, control of boundary layer, etc. are few prominent magneto-hydrodynamics applications. Alshomrani and Gul (2017) examined the thin-film flow of water-based Al_{2}O_{3} and Cu nanofluid through a stretched cylinder under the impact of magnetic function. Sandeep and Malvandi (2016) investigated the features of magneto thermal transport containing the time-dependent flow of liquid movement of the nanofluid thin film of a starched surface. The features of hybrid nanoliquid flow with variable heat and drag forces were explored by Sandeep (2017). Ahmad and Nadeem (2017) scrutinized the magnetic impact on micropolar fluid comprising hybrid nanofluid with a heat sink/source and obtained the multiple results for hybrid nanofluid and micropolar fluid. Hamrelaine et al. (2019) investigated the magnetic impact of Jeffery-Hamel flow between non-parallel permeable walls or permeable plates. The magnetic effect on the radiative flow of the hybrid nanoliquid thin film with irregular heat sink/source was inspected by Kumar et al. (2020). Zaib et al. (2020) obtained the similarity of multiple results from magnetite ferroliquid conveying non-Newtonian blood flow with entropy generation. Wakif et al. (2019) estimated the impact of the magnetic field on flows of Stokes' second problem with entropy generation. Recently, Wakif (2020) examined the magneto flow of Casson blood flow from a horizontal stretched surface with rough geometry, thermal conductivity, and temperature viscosity.

As the literature discussed above, the non-linear radiative flow of hybrid ferroliquids from a squeezed mixed convective flow within two disks has not yet been explored. Thus, this investigation's core objective is to examine the non-linear radiative magneto flow of hybrid ferrofluids in a squeezed flow within two disks. Such sort of problems is created by biomedical enforcement and electric generators. The hybrid nanoparticles have constructive biomedical applications like cancer therapy, drug delivery, biosensing, etc. The liquid model is altered into ODEs through acceptable technique and investigated numerically through the Lobatto 3A formula.

## Problem Formulation

Figure 1 demonstrates the problem's physical background of the flow configuration schematically. Consider a 2D magneto-hydrodynamic (MHD), incompressible, and squeezing flow of a hybrid ferrofluid along with axisymmetric and time-dependent flow between the two infinite parallel disks alienated by a distance ${h}_{1}(t)={H}_{1}{(1-bt)}^{0.5},$ where *b* is called the characteristic parameter whose dimension is *t*^{−1}. Further, interpretation of this point shows that the squeezing movements of both the disk are occurring in the case *b* > 0 with velocity *w*_{1} = *dh*_{1}/*dt* until they touch each other at *bt* = 1 while the case *b* < 0 refers to the movement of the disk away to each other. The initial or starting point of the disk is represented by *H*_{1} at time *t* = 0. The regular liquid is caught to be 50%:50% (ethylene glycol + water). The time-dependent magnetic field is proportional to ${B}_{1}(t)={B}_{0}/{(1-bt)}^{0.5}$ and is enforced normally to the disks. Also, the field of induced magnetic field is ignored owing to the assumption of small Reynold's number. The influences of mixed convective and thermal radiation are invoked in the governing equations. *T*_{w} refers to the constant temperature of the lower disk *z*_{1} = 0, while *T*_{h1} called the constant temperature of the upper disk located at *z*_{1} = *h*_{1} (*t*). The upper disk located at *z*_{1} = *h*_{1} (*t*) is moving away or toward with the velocity $0.5b{H}_{1}{(1-bt)}^{0.5}$ from the stationary lower disk at *z*_{1} = 0. Based on these suppositions, the prevailing constitutive equations of the momentum and energy for the time-dependent viscous hybrid ferrofluid are (Hayat et al., 2012a, 2018a; Babazadeh et al., 2020).

with the boundary conditions (Babazadeh et al., 2020)

where the pressure and temperature of the flow of fluid are represented by the symbols *p* and *T*_{1}, respectively. The rest of the mathematical notions and letters of symbols in the Equations (1)–(6) are defined in the nomenclature, while the thermophysical properties of the normal fluid and the hybrid ferrofluid are shown in Table 1.

**Table 1**. Properties of thermophysical of normal liquid and ferroparticles (Rashad, 2017; Kumar et al., 2020).

In addition, the mathematical expressions of the studied hybrid ferrofluid are given by:

where ϕ is the total volume fraction of ferroparticles. The volume fractions of the individual ferroparticles are indicated by the symbol ϕ_{1} and ϕ_{2}, while the subscripts *hnf, f*, and *s* represent the hybrid ferrofluid, base fluid, and solid particles. Further, to ease this analysis, the similarity variables are introduced as the following

Now, using the aforementioned transformation along with Equation (6), the governing PDEs, Equations (2) to (4), along with the boundary condition (5) becomes

Eliminating the pressure term from Equations (8) and (9), the final form of the momentum and energy equation becomes

with the boundary restrictions

In which:

where the dimensional parameters contained in the equations above are the following: $S=\frac{b{{H}_{1}}^{2}}{2{\nu}_{f}},$ ${M}^{2}=\frac{{\sigma}_{f}{B}_{0}^{2}{{H}_{1}}^{2}}{{\rho}_{f}{\nu}_{f}},$ $\lambda =\frac{2{{H}_{1}}^{2}{(1-bt)}^{2}g{\beta}_{f}({T}_{w}-{T}_{{h}_{1}})}{b{r}_{1}{\nu}_{f}},$ ${R}_{d}=\frac{4{\sigma}^{*}{T}_{{h}_{1}}^{3}}{{k}^{*}{k}_{f}},$ ${\theta}_{{h}_{1}}=\frac{{T}_{w}}{{T}_{{h}_{1}}}$, $A=\frac{{w}_{0}}{b{H}_{1}}$, and $Pr=\frac{{\nu}_{f}}{{\alpha}_{f}}$ while these parameters are namely called the squeeze number, Hartmann number or magnetic, mixed convective parameter, radiation parameter, temperature ratio parameter, the blowing/suction parameter, and the Prandtl number, respectively.

## Engineering Quantities of Interest

The skin friction factor and local Nusselt number are signified as:

where

In terms of dimensional form, the local skin friction and Nusselt number are

where ${Re}_{{r}_{1}}=\frac{{r}_{1}b{H}_{1}\sqrt{1-bt}}{2{\nu}_{f}}$ is the local Reynolds number.

## Numerical Procedure

The higher-order Equations (11) and (12) along with the boundary restriction (13) are a highly non-linear dimensional form of ODEs in nature. The solution of these equations is very difficult to obtain analytically. Therefore, we utilized the numerical approach of bvp4c, which is based on the three-stage Lobatto 3A method. This process is fruitful to handle the one-point problem in the first order. In broad mode, Lobatto 3A techniques have been exercised for the two-point boundary value problem due to their very exceptional properties of stability and the fact that they contain accuracy of fourth order over the entire interval. To execute this approach, the equations in the dimensional form are exercised owing to new variables and converted into the first-order system of equations. To do the technique as outlined, let

Plugging Equation (17) in the dimensional form of ODEs to get the first-order set of ODEs as follows:

Subject to corresponding initial conditions

To carry out this procedure, we have to choose the finite value of the boundary η → ∞, say η _{∞} and the computated value of this constraint depends on the other governing constraints involved in the system of the first-order ODEs. The aforementioned Equation (19) contains one unknown *A*, which is constant, and the value depends on the other constraints that influenced the system. In the current computation process, the value of the unknown constant *A* is fixed to be 0.2, while the other involved sundry constraints whose numerical fixed values are mentioned in the portion of Results and Discussion. Thus, to run the process of bvp4c, it required some initial early guesses that satisfied Equation (19). Moreover, the problem is bounded between 0 and 1. So, in the computation, the value of the step size is reserved as Δη = 0.01. The iterative process is repeated, waiting for the required tests to be obtained in compliance with the convergence criterion up to the accuracy point 10^{−10}.

## Results and Discussion

The physical trend of the engaged parameters on the temperature distribution and nanofluid velocity was investigated. During the computation, it allocated that the pertaining parameters be fixed as ϕ = 0.03, *M* = 0.3, *R*_{d} = 0.5, θ_{h1} = 0.5, *S* = 0.5, *A* = 0.2, and λ = 0.5, except for the values depicted in the tables and graphs. Table 2 illustrates the validation of the current outcomes of $\frac{{{H}_{1}}^{2}}{{{r}_{1}}^{2}}{Re}_{{r}_{1}}{C}_{F{r}_{1}}$ with Hayat et al.'s (2012a) existing results. It confirmed the outstanding harmony with available results. Tables 3, 4 are set to scrutinize the behavior of the physical parameters in the form of the values of $\frac{{{H}_{1}}^{2}}{{{r}_{1}}^{2}}{Re}_{{r}_{1}}{C}_{F{r}_{1}}$ and $\sqrt{1-bt}N{u}_{{r}_{1}}$ in the opposing and assisting flows. The friction factor and heat transport rate decline owing to *M* in the assisting flow and augment in the opposing flow. Whereas, the ferroparticle volume fraction tremendously augments the heat transport rate and friction factor in the assisting flow, a divergent performance is seen in the opposing flow. Physically, the greater amount of ϕ creates a larger energy conveyed via the fluid flow connected with the unbalanced progress of the ultrafine materials and thus generates a significant improvement in the friction factor and the heat transport rate. The squeezed parameter augments the friction factor and the heat transport rate in both flows.

**Table 2**. Justification of the outcomes for $\frac{{{H}_{1}}^{2}}{{{r}_{1}}^{2}}{Re}_{{r}_{1}}{C}_{F{r}_{1}}$ when *R*_{d} = θ_{h1} = λ = ϕ = 0.

**Table 3**. The numerical values of $\frac{{{H}_{1}}^{2}}{{{r}_{1}}^{2}}{Re}_{{r}_{1}}{C}_{F{r}_{1}}$ and $\sqrt{1-bt}N{u}_{{r}_{1}}$ for the assisting flow λ = 0.5 when *R*_{d} = θ_{h1} = 1, *A* = 0.2.

**Table 4**. The numerical values of $\frac{{{H}_{1}}^{2}}{{{r}_{1}}^{2}}{Re}_{{r}_{1}}{C}_{F{r}_{1}}$ and $\sqrt{1-bt}N{u}_{{r}_{1}}$ for the opposing flow λ = −0.5 when *R*_{d} = θ_{h1} = 1, *A* = 0.2.

The influence of the volume fraction of hybrid ferrofluid (ϕ) on the fluid velocity (*F*′(η)) and temperature (θ(η)) for the assisting flow (ASF) and the opposing flow (OPF) is revealed in Figures 2, 3, respectively. Figure 2 portrays that the fluid velocity increases for the ASF and decreases for the OPF as we augmented the values of (ϕ) simultaneously up to the approximate range of η_{max} = 0.3 and η_{max} = 0.58, respectively, while the reverse behavior of the velocity is seen in both the ASF and the OPF when the similarity variable is greater than the aforementioned ranges. The solution in the opposing flow is more dominant as compared to the assisting flow in the velocity profile. This is apparent from reality that a greater amount of (ϕ) communicates to increase the ferrofluid thermal conductivity, which inspires the diffusion of heat in order for the heat to impulsively disperse near the surface of the disks. Figure 3 demonstrates the temperature distribution for the ASF as well as the OPF against the similarity variable for the hybrid ferrofluid (ϕ). The temperature distribution diminishes in the ASF and upsurges in the OPF as we boost up the value of the hybrid ferrofluid (ϕ). The large amount of (ϕ) generates a great energy conveyed during the flow connected with the uneven progress of the ultrafine materials and hence produces a substantial improvement in the rate of heat transport that consequently augments the temperature. It is fascinating to see in the figure that the thermal boundary layer in the ASF is less extensive than the OPF. Figures 4, 5 show the Hartmann impact *M* on (*F*′(η)) and (θ(η)) against the similarity variable η for the phenomenon of both the ASF and OPF, respectively. In case of assisting flow, the fluid velocity reduces as we move along the horizontal axis within the range of η_{max} = 0.51, and in the same spotted range, the velocity is enhanced in the opposing flow for the continuous positive change in the values of *M*, while for the simultaneous cases, the opposite behavior is, respectively, shown as we move out of the aforementioned range. Further, it is scrutinized that the liquid velocity is maximum attained for the ASF than that of OPF. Figure 5 exhibits that the temperature distribution grows higher in the ASF and decreases in the OPF for the regular augmentation in the parameter *M*. More generally, the temperature behavior of the fluid is enhanced in the ASF due to the higher values of *M*. The constant augmentations in *M* will reduce the thickness of the thermal boundary layer of the fluid regularly. Physically, the escorting hybrid ferrofluid produces a sort of resistive force, that is to say, the Lorentz force. Together with this, the flow generates a certain amount of heat energy. It is worth mentioning that the hybrid ferroliquid is extremely prejudiced through the drag force in the ASF compared with OPF. The impact of the squeeze number *S* for both the cases of heating and cooled surfaces on the velocity and the distribution of temperature is plotted in Figures 6, 7, respectively. Here, we have taken the values of squeeze number to be greater than zero, which physically shows that the upper disk is moving away from the lower stationary disk, while in the case of opposite inequality of the squeeze number, corresponding to the lower disk moving away from the upper disk, this was not discussed graphically. More precisely, Figure 6 shows that the velocity is declining against the dimensionless variable η in the domain from 0 ≤ η ≤ 0.6 for the ASF when we improve the value of *S*, while the flow velocity is jumping up in the same domain for the OPF. When we altered the domain from 0.6 ≤ η ≤ 01, then the opposite behavior (like 0 ≤ η ≤ 0.6) for the flow velocity is accomplished. Despite this, the momentum boundary layer is lower in the OPF when we compared with the ASF. The variation in the temperature distribution is shown in Figure 7, where the distribution in the temperature of the fluid flow is consistently falling as we augment the value of *S* for both the ASF and OPF. It is interestingly seen in the pictures that both the thermal boundary layer and the distribution in the temperature are declining. The last profile of velocity is behaviorally exposed in Figure 8 for suction *A* > 0 and blowing *A* < 0. The value of the velocity decreases in both the ASF and OPF for the four selected values of (*A* = 0, 1, 1.5, 2) while the opposite trend of the velocity is achieved in the case of injection (*A* = 0, −1, −1.5, −2). In contrast to suction, the momentum boundary layer is more dominant in the injection case. The impact of the blowing/suction parameter *A* in the temperature profile against the new variable η for both the positive and negative values of λ is depicted in Figure 9. A very similar behavior is observed (such as in Figure 8), but the temperature and thermal boundary layer continue to fall due to the parameter value of (*A* = 0, 1, 1.5, 2) as compared to the parameter in the case of injection (*A* = 0, −1, −1.5, −2). Figures 10, 11 display the stimulus of the parameters *R*_{d} and θ_{h1} on (θ(η)) for the assisting and opposing flows, respectively. Figure 10 shows the increasing behavior of temperature in the ASF and in the OPF for the greater value of *R*_{d}. However, it is widely known that its radiation process is just the heat transfer phenomenon that emits the energy through fluid particles so that an extra heat is generated in the flow. Thus, it is identified as a development in the thermal boundary layer for greater *R*_{d}. Behavior-wise, one can get a similar temperature distribution (like in Figure 10), which is highlighted in Figure 11. We noticed very motivating results in that the heat curves θ(η) are further exaggerated in the impact of radiation parameter when we compared with the temperature ratio parameter. Finally, Figures 12, 13 are sketched to observe the behavior of streamlines for suction and blowing. These patterns suggest that streamlines' behavior in case of suction is smooth, while in the case of blowing, the streamlines are divided into two regions.

## Conclusion

This research succeeded in conferring the exploration of the mixed convective squeezed flow of hybrid ferroparticles within the two discs under the impacts of MHD and radiation. Water + ethylene 50%:50% is taken as a base liquid of the channel. The leading equations are resolved through the Lobatto 3A formula. To ensure the correctness of the outcomes, an assessment of existing results is made. The influences of numerous significant constraints on the fluid flow are consummated via graphs. The following significant outcomes are summarized:

• The velocity of hybrid nanofluid upsurges due to ϕ in ASF and decelerates in OPF, while the temperature declines in ASF and augments in OPF.

• The Hartmann number upsurges and shrinks the temperature, respectively, in ASF and OPF, while the opposing trend is seen in the velocity.

• The squeezed number declines the temperature in both flows and augments the velocity in the OPF and shrinks that in the ASF.

• The suction parameter augments the velocity and temperature, while the injection shrinks the velocity and temperature in both flows.

• The radiation and temperature ratio parameters upsurge the temperature in both flows.

• The friction factor and heat transport rate decline due to *M* in the ASF and augment in the OPF.

• The ferroparticle volume fraction tremendously augments the friction factor and the heat transport rate in the ASF, while the opposing behavior is seen in the OPF.

In the end, this work can be further considered by taking hybrid nanoparticles with different shape factors along with mass transfer. Also, this problem may be extended into the two-phase model by taking into account the Buongiorno model.

## Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

## Author Contributions

UK did formulation and transformed the problem. AZ computed results. IK, KN, and DB wrote the manuscript. DB, IK, and KN revised the manuscript. All authors contributed to the article and approved the submitted version.

## Funding

This project was supported by the Deanship of Scientific Research at Prince Sattam Bin Abdulaziz University under the research project No. 2020/01/16436.

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

Adriana, M. A. (2017). Hybrid nanofluids based on Al_{2}O_{3}, TiO_{2} and SiO_{2}: numerical evaluation of different approaches. *Int. J. Heat Mass Transf*. 104, 852–860. doi: 10.1016/j.ijheatmasstransfer.2016.09.012

Ahmad, A., and Nadeem, S. (2017). Effects of magnetohydrodynamics and hybrid nanoparticles on a micropolar fluid with 6-types of stenosis. *Res. Phys*. 7, 4130–4139. doi: 10.1016/j.rinp.2017.10.032

Ali, L., Liu, X., Mujeed, S., Abdal, S., and Khan, S. A. (2020). Analysis of magnetic properties of nano-particles due to a magnetic dipole in micropolar fluid flow over a stretching sheet. *Coatings* 10:170. doi: 10.3390/coatings10020170

Alshomrani, A. S., and Gul, T. (2017). A convective study of Al_{2}O_{3}-H_{2}O and Cu–H_{2}O nanoliquid films sprayed over a stretching cylinder with viscous dissipation. *Eur. Phys. J. Plus* 132:495. doi: 10.1140/epjp/i2017-11740-1

Alzahrani, A. K., Malik, Z. U., and Muhammad, T. (2020). Numerical treatment for 3D squeezed flow in a rotating channel with Soret and Dufour effects. *Front. Phys*. 8:201. doi: 10.3389/fphy.2020.00201

Babazadeh, H., Muhammad, T., Shakeriaski, F., Ramzan, M., and Hajizadeh, M. R. (2020). Nanomaterial between two plates which are squeezed with impose magnetic force. *J. Therm. Anal. Calorim*. doi: 10.1007/s10973-020-09619-6

Chamkha, A. J., Abbasbandy, S., Rashad, A. M., and Vajravelu, K. (2012). Radiation effects on mixed convection over a wedge embedded in a porous medium filled with a nanofluid. *Trans. Porous Med*. 91, 261–279. doi: 10.1007/s11242-011-9843-5

Choi, S. U. S. (1995). “Enhancing thermal conductivity of fluids with nanoparticle,” in *Developments and Applications of Non-Newtonian Flows*, Vol. 231/MD 66, eds D. A. Siginer and H. P. Wang (New York, NY: ASME FED), 99–105.

Domairry, G., and Aziz, A. (2009). Approximate analysis of MHD Squeeze flow between two parallel disks with suction or injection by homotopy perturbation method. *Math. Prob. Eng*. 2009, 1–19. doi: 10.1155/2009/603916

Domairry, G., and Hatami, H. (2014). Squeezing Cu-water nanofluid flow analysis between parallel plates by DTM-Pade method. *J. Mol. Liq*. 193, 37–44. doi: 10.1016/j.molliq.2013.12.034

Duwairi, H. M., Tashtoush, B., and Domseh, R. A. (2004). On heat transfer effects of a viscous fluid squeezed and extruded between two parallel plates. *Heat Mass Transf*. 14, 112–117. doi: 10.1007/s00231-004-0525-5

El-Kabeir, S. M. M., Modather, M., and Rashad, A. M. (2015). Effect of thermal radiation on mixed convection flow of a nanofluid about a solid sphere in a saturated porous medium under convective boundary condition. *J. Porous Med*. 18, 569–584. doi: 10.1615/JPorMedia.v18.i6.20

Hamrelaine, S., Mebarek-Oudina, F., and Sari, M. R. (2019). Analysis of MHD Jeffery Hamel flow with suction/injection by homotopy analysis method. *J. Adv. Res. Fluid Mech. Therm. Sci*. 58, 173–186.

Hayat, T., Ayub, T., Muhammad, T., Alsaedi, A., and Mustafa, M. (2018a). On Darcy-Forchheimer squeezed flow of carbon nanotubes between two parallel disks. *Int. J. Numerical Methods Heat Fluid Flow* 28, 2784–2800. doi: 10.1108/HFF-06-2017-0250

Hayat, T., Haider, F., Muhammad, T., and Alsaedi, A. (2018b). Darcy–Forchheimer squeezed flow of carbon nanotubes with thermal radiation. *J. Phys. Chem. Solids* 120, 79–86. doi: 10.1016/j.jpcs.2018.04.016

Hayat, T., and Nadeem, S. (2017). Heat transfer enhancement with Ag–CuO/ water hybrid nanofluid. *Res. Phys*. 7, 2317–2324. doi: 10.1016/j.rinp.2017.06.034

Hayat, T., Yousaf, A., Mustafa, M., and Asghar, S. (2012a). Influence of heat transfer in the squeezing flow between parallel disks. *Chem. Eng. Comm*. 199, 1044–1062. doi: 10.1080/00986445.2011.631203

Hayat, T., Yousaf, A., Mustafa, M., and Obaidat, S. (2012b). MHD squeezing flow of second-grade fluid between two parallel disks. *Int. J. Numer. Methods Fluids* 69, 399–410. doi: 10.1002/fld.2565

Ishizawa, S. (1966). The unsteady flow between two parallel discs with arbitrary varying gap width. *Bull. Jpn. Soc. Mech. Eng*. 9, 533–550. doi: 10.1299/jsme1958.9.533

Khan, W. A., Khan, Z. H., and Haq, R. U. (2015). Flow and heat transfer of ferrofluids over a flat plate with uniform heat flux. *Eur. Phys. J. Plus* 130:86. doi: 10.1140/epjp/i2015-15086-4

Kumar, K. A., Sandeep, N., Sugunamma, V., and Animasaun, I. L. (2020). Effect of irregular heat source/sink on the radiative thin film flow of MHD hybrid ferrofluid. *J. Therm. Anal. Cal.* 139, 2145–2153. doi: 10.1007/s10973-019-08628-4

Madhesh, D., and Kalaiselvam, S. (2015). Experimental study on heat transfer and rheological characteristics of hybrid nanofluids for cooling applications. *J. Exp. Nanosci*. 10, 1194–1213. doi: 10.1080/17458080.2014.989551

Mahanthesh, B., Gireesha, B. J., Animasaun, I. L., Muhammad, T., and Shashikumar, N. S. (2019). MHD flow of SWCNT and MWCNT nanoliquids past a rotating stretchable disk with thermal and exponential space dependent heat source. *Phys. Scr*. 94:085214. doi: 10.1088/1402-4896/ab18ba

Marzougui, S., Mebarek-Oudina, F., Assia, A., Magherbi, M., Shah, Z., and Ramesh, K. (2020). Entropy generation on magneto convective flow of copper–water nanofluid in a cavity with chamfers. *J. Therm. Anal. Calorim*. doi: 10.1007/s10973-020-09662-3

Mebarek-Oudina, F. (2019). Convective heat transfer of Titania nanofluids of different base fluids in cylindrical annulus with discrete heat source. *Heat Transf. Asian Res*. 48, 135–147. doi: 10.1002/htj.21375

Mebarek-Oudina, F., Aissa, A., Mahanthesh, B., and Öztop, H. F. (2020). Heat transport of magnetized Newtonian nanoliquids in an annular space between porous vertical cylinders with discrete heat source. *Int. Commun. Heat Mass Transf*. 117:104737. doi: 10.1016/j.icheatmasstransfer.2020.104737

Neuringer, J. L. (1996). Some viscous flows of a saturated ferrofluid under the combined influence of thermal and magnetic field gradients. *Int. J. Non Linear Mech*. 1, 123–137. doi: 10.1016/0020-7462(66)90025-4

Qayyum, A., Awais, M., Alsaedi, A., and Hayat, T. (2012). Unsteady squeezing flow of Jeffery fluid between two parallel disks. *Chin. Phys. Lett*. 29:034701. doi: 10.1088/0256-307X/29/3/034701

Rashad, A. M. (2017). Impact of thermal radiation on MHD slip flow of a ferrofluid over a non-isothermal wedge. *J. Mag. Magnetic Mat*. 422, 25–31. doi: 10.1016/j.jmmm.2016.08.056

Raza, J., Mebarek-Oudina, F., Ram, P., and Makinde, O. D. (2020). MHD flow of non-Newtonian molybdenum disulfide nanofluid in a converging/diverging channel with Rosseland radiation. *Defect Diffus. Forum* 401, 92–106. doi: 10.4028/www.scientific.net/DDF.401.92

Sandeep, N. (2017). Effect of aligned magnetic field on liquid film flow of magnetic-nanofluid embedded with graphene nanoparticles. *Adv. Powder Technol*. 28, 865–875. doi: 10.1016/j.apt.2016.12.012

Sandeep, N., and Malvandi, A. (2016). Enhanced heat transfer in liquid thin film flow of non-Newtonian nanofluids embedded with grapheme nanoparticles. *Adv. Powder Technol*. 27, 2448–2456. doi: 10.1016/j.apt.2016.08.023

Sheikholeslami, M. (2019). Numerical approach for MHD Al_{2}O_{3}-water nanofluid transportation inside a permeable medium using innovative computer method. *Comput. Methods Appl. Mech. Eng*. 344, 306–318. doi: 10.1016/j.cma.2018.09.042

Sheikholeslami, M., and Ganji, D. D. (2013). Heat transfer of Cu-water nanofluid flow between parallel plates. *Powder Technol*. 235, 873–879. doi: 10.1016/j.powtec.2012.11.030

Sheikholeslami, M., Ganji, D. D., and Ashorynejad, H. R. (2013). Investigation of squeezing unsteady nanofluid flow using ADM. *Powder Technol*. 239, 259–265. doi: 10.1016/j.powtec.2013.02.006

Soomro, F. A., Zaib, A., Haq, R. U., and Sheikholeslami, M. (2019). Dual nature solution of water functionalized copper nanoparticles along a permeable shrinking cylinder: FDM approach. *Int. J. Heat Mass Transf*. 129, 1242–1249. doi: 10.1016/j.ijheatmasstransfer.2018.10.051

Suresh, S., Venkitaraj, K. P., Selvakumar, P., and Chandrasekar, M. (2012). Effect of Al_{2}O_{3}-Cu/ water hybrid nanofluid in heat transfer. *Exp. Therm. Fluid Sci*. 38, 54–60. doi: 10.1016/j.expthermflusci.2011.11.007

Tiwari, R. K., and Das, M. K. (2007). Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. *Int. J. Heat Mass Transf*. 50, 2002–2018. doi: 10.1016/j.ijheatmasstransfer.2006.09.034

Usha, R., and Sridharan, R. (1996). Arbitrary squeezing of a viscous fluid between elliptic plates. *Fluid Dyn. Res*. 18, 35–51. doi: 10.1016/0169-5983(96)00002-0

Wakif, A. (2020). A novel numerical procedure for simulating steady MHD convective flows of radiative Casson fluids over a horizontal stretching sheet with irregular geometry under the combined influence of temperature-dependent viscosity and thermal conductivity. *Math. Prob. Eng*. 2020, 1–20. doi: 10.1155/2020/1675350

Wakif, A., Boulahia, Z., Mishra, S. R., Rashidi, M. M., and Sehaqui, R. (2018). Influence of a uniform transverse magnetic field on the thermo-hydrodynamic stability in water-based nanofluids with metallic nanoparticles using the generalized Buongiorno's mathematical model. *Eur. Phys. J. Plus* 133:181. doi: 10.1140/epjp/i2018-12037-7

Wakif, A., Chamkha, A. J., Thumma, T., Animasaun, I. L., and Sehaqui, R. (2020). Thermal radiation and surface roughness effects on the thermo-magneto-hydrodynamic stability of alumina–copper oxide hybrid nanofluids utilizing the generalized Buongiorno's nanofluid model. *J. Therm. Anal. Calorim.* doi: 10.1007/s10973-020-09488-z

Wakif, A., Qasim, M., Afridi, M. I., Saleem, S., and Al-Qarni, M. M. (2019). Numerical examination of the entropic energy harvesting in a magnetohydrodynamic dissipative flow of Stokes' second problem: utilization of the gear-generalized differential quadrature method. *J. Non Equil. Thermo*. 44, 385–403. doi: 10.1515/jnet-2018-0099

Zaib, A., Khan, U., Shah, Z., Kumam, P., and Thounthong, P. (2019). Optimization of entropy generation in flow of micropolar mixed convective magnetite (Fe_{3}O_{4}) ferroparticle over a vertical plate. *Alexandria Eng. J*. 58, 1461–1470. doi: 10.1016/j.aej.2019.11.019

Zaib, A., Khan, U., Wakif, A., and Zaydan, M. (2020). Numerical entropic analysis of mixed MHD convective flows from a non-isothermal vertical flat plate for radiative tangent hyperbolic blood biofluids conveying magnetite ferroparticles: dual similarity solutions. *Arabian J. Sci. Eng*. 45, 5311–5330. doi: 10.1007/s13369-020-04393-x

## Nomenclature

Keywords: squeeze flow, mixed convection, magnetize hybrid ferrofluids, non-linear radiation, disks

Citation: Nisar KS, Khan U, Zaib A, Khan I and Baleanu D (2020) Numerical Simulation of Mixed Convection Squeezing Flow of a Hybrid Nanofluid Containing Magnetized Ferroparticles in 50%:50% of Ethylene Glycol–Water Mixture Base Fluids Between Two Disks With the Presence of a Non-linear Thermal Radiation Heat Flux. *Front. Chem.* 8:792. doi: 10.3389/fchem.2020.00792

Received: 28 May 2020; Accepted: 29 July 2020;

Published: 23 September 2020.

Edited by:

Abderrahim Wakif, University of Hassan II Casablanca, MoroccoReviewed by:

Taseer Muhammad, King Khalid University, Saudi ArabiaMostafa Zaydan, University of Hassan II Casablanca, Morocco

Fateh Mebarek-Oudina, University of Skikda, Algeria

Copyright © 2020 Nisar, Khan, Zaib, Khan and Baleanu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ilyas Khan, ilyaskhan@tdtu.edu.vn