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

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.

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

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;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 singlephase 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 timedependent 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. Figure 1 demonstrates the problem's physical background of the flow configuration schematically. Consider a 2D magnetohydrodynamic (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.

PROBLEM FORMULATION
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 h 1 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.5bH 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(Hayat et al., , 2018aBabazadeh 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.
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) 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: bH 1 , and Pr = ν f α 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 τ r 1 z 1 = µ hnf ∂u 1 ∂z 1 + ∂w 1 ∂r 1 z 1 =h 1 (t) , . (15) In terms of dimensional form, the local skin friction and Nusselt number are where Re r 1 = r 1 bH 1 √ 1−bt 2ν 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 onepoint 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: Frontiers in Chemistry | www.frontiersin.org 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 . FIGURE 2 | Impact of φ 1 , φ 2 on F ′ (η).

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, θ h 1 = 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 H 1 2 r 1 2 Re r 1 C Fr 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 H 1 2 r 1 2 Re r 1 C Fr 1 and √ 1 − btNu 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 FIGURE 3 | Impact of φ 1 , φ 2 on θ(η). 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.
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 θ h 1 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.