Thermal investigation into the Oldroyd-B hybrid nanofluid with the slip and Newtonian heating effect: Atangana–Baleanu fractional simulation

The significance of thermal conductivity, convection, and heat transportation of hybrid nanofluids (HNFs) based on different nanoparticles has enhanced an integral part in numerous industrial and natural processes. In this article, a fractionalized Oldroyd-B HNF along with other significant effects, such as Newtonian heating, constant concentration, and the wall slip condition on temperature close to an infinitely vertical flat plate, is examined. Aluminum oxide (Al2O3) and ferro-ferric oxide (Fe3O4) are the supposed nanoparticles, and water (H2O) and sodium alginate (C6H9NaO7) serve as the base fluids. For generalized memory effects, an innovative fractional model is developed based on the recently proposed Atangana–Baleanu time-fractional (AB) derivative through generalized Fourier and Fick’s law. This Laplace transform technique is used to solve the fractional governing equations of dimensionless temperature, velocity, and concentration profiles. The physical effects of diverse flow parameters are discussed and exhibited graphically by Mathcad software. We have considered 0.15 ≤ α ≤ 0.85 , 2 ≤ Pr ⁡ ≤ 9 , 5 ≤ G r ≤ 20 , 0.2 ≤ ϕ 1 , ϕ 2 ≤ 0.8 , 3.5 ≤ G m ≤ 8 , 0.1 ≤ Sc ≤ 0.8 , and 0.3 ≤ λ 1 , λ 2 ≤ 1.7 . Moreover, for validation of our present results, some limiting models, such as classical Maxwell and Newtonian fluid models, are recovered from the fractional Oldroyd-B fluid model. Furthermore, comparing the results between Oldroyd-B, Maxwell, and viscous fluid models for both classical and fractional cases, Stehfest and Tzou numerical methods are also employed to secure the validity of our solutions. Moreover, it is visualized that for a short time, temperature and momentum profiles are decayed for larger values of α , and this effect is reversed for a long time. Furthermore, the energy and velocity profiles are higher for water-based HNFs than those for the sodium alginate-based HNF.


Introduction
With the addition of nanometer-sized particles in various base fluids, thermophysical characteristics may improve in energy transfer schemes. This process signals an expansion in the thermal conductivity for base fluids, making it more reliable and ongoing. These significant fluids define nanofluids (NFs) with an extensive series of suggestions in several areas of science, as well as technology, with nuclear devices, heat exchangers, solar plates, vehicle heaters, and biotic and organic devices (Usman et al., 2018;Khan et al., 2022a;Khan et al., 2022b;Ahmed et al., 2022;Hassan et al., 2022;Khan et al., 2022c). First, Lee and Eastman presented the idea of NFs in 1995 (Lee et al., 1999). Numerous applications of NFs are discoursed by Kaufui et al. (Wong and Omar De Leon., 2010). Mahian et al. (2019) proposed important ideas and reflected novel innovations to completely explain the NFs. They were obsessed with innovative expansions in this field, comprehensive explanations of the thermophysical characteristics, and imitation of thermal transmission in NF flow. Waini et al. (2019) used a numerical scheme to discuss an unsteady thermal transmission flow past a shrinking sheet in an HNF. They presented different applications of NFs in numerous branches of science along with appreciated recommendations. NFs have achieved significant consideration from researchers due to their improved heat conversion characteristics. The rheological presentation of an NF using a revolving rheometer was proposed by Vallejo et al. (2019a). Different rheological characteristics of NFs are discussed in Vallejo et al. (2019b). Currently, NFs have been characterized as HNFs in several mechanisms (Rashad et al., 2018). HNFs are developed by mixing two dissimilar nanoparticles in the base liquid. Its main inspiration is to increase the thermal features of NFs. The variable thermal transmission of HNFs through magnetic influence was examined in Mohebbi et al. (2019). The heat transmission in the non-Newtonian HNF composed with entropy generation was discussed in Shahsavar et al., 2018). Furthermore, Farooq et al. (2018) deliberated on the entropy in the HNF flow in a stretching sheet. Asogwa et al. (2021) discussed chemical reactions and heat sinks over a ramped temperature. The analytical solution of governing equations was found with the Laplace transform. Asogwa et al. (2022a) used the Laplace approach to discuss a water-based NF containing aluminum oxide and copper in a moving plate and proved that thermal absorption causes a decline in aluminum oxide NF's thermal and momentum profiles with a copper NF. Shankar Goud et al. (2022) used the Keller-box scheme for the numerical solution along with thermal effects, momentum, and solutal slip on the thermal transmission with a description of the magnetohydrodynamic (MHD) flow of Casson fluid and an exponential porous surface with Dufour, chemical reaction, and Soret impacts. Khan et al. (2022d) studied a fractionalized electro-osmotic flow based on the Caputo operator of a Casson NF containing sodium alginate nanoparticles over a vertical microchannel with MHD effects. They proved that the inclination angle boosts the velocity. Asogwa et al. (2022b) and Asogwa et al. (2022c) considered the stimulation significance of the thermal transmission with the MHD flow of a NF through an extending sheet with MATLAB bvp4c. Furthermore, they investigated the radiative features of the MHD flow with collective heat transportation characteristics on a reactive stretching surface with the Casson NF numerically using MATLAB bvp4c. Goud et al. (2022) applied the bvp4c scheme to study the convection flow via an infinite porous plate on thermal transmission, as well as mass transmission. Asogwa et al. (2022d) discussed the influences of the movement of nanoparticles in NFs by an exponentially enhanced Riga plate. Reddy et al. (2022) calculated the effect of activation energy on a second-grade MHD NF flow over a convectively curved heated stretched surface by considering the Brownian motion and generation/absorption, and thermophoresis. They have shown that velocity and thermal profiles suggestively increase with the concurrent increasing estimation of the fluid parameter.
The fractional calculus (FC) has obtained substantial consideration from experts in previous decades. The important inventions have newly been presented in the application of the FC, where new derivatives, as well as integral operators, are hired (Awan et al., 2019). The new anticipated operators contain the generalized Mittag-Leffler function (MLF), and these features intensify the innovative constructions to achieve numerous attractive properties that are recognized in important outcomes. Subsequently, Atangana and Dumitru (2016) anticipated, the innovative and applicable time-fractional operator, which is expansively hired in numerous branches of science and engineering. It is exposed that the MLF is a more operative and vigorous screening apparatus than the exponential and power laws, constructing the AB-fractional operator, in terms of Caputo, an effective arithmetic procedure to simulate progressively perilous complex tasks. Due to their extensive implications, such fractional models are extensively identified for deriving fractional differential equations (FDEs) with no manufactured irregularities, as for Caputo, Riemann-Liouville (RL), and Caputo-Fabrizio (CF) derivatives, because of their characteristic non-orientation (Ali et al., 2021;Ali et al., 2022a;Raza et al., 2022;Zhang et al., 2022). We also perceived interest in these fractional derivatives on the topic of mathematical approaches, although scientifically approximating these operators' outcomes to compute different problems (Martyushev and Sheremet, 2012;Ali et al., 2022b). Batool et al. (2022) discussed the thermal and mass transmission processes of a micropolar NF under magnetic and buoyancy effects across an inclusion. Rasool et al. (2022a) examined the significance of the MHD Maxwell NF flow and obtained the solution to this problem by employing the homotopy analysis technique for diverse physical parameters. Moreover, they studied an electro-magnetohydrodynamic NF flow in a permeable medium with heating Frontiers in Materials frontiersin.org boundary conditions. Furthermore, they applied Buongiorno's method for the flow of radiating thixotropic NFs over a horizontal surface by considering the retardational effects of Lorentz forces and using the influence of Brownian and thermophoresis diffusions (Rasool et al., 2022b;Rasool et al., 2023). In this paper, a fractionalized Oldroyd-B HNF flow is examined by the recent definitions of the AB time-fractional derivative having a Mittage-Leffler kernel along with Newtonian heating, constant concentration, and the wall slip condition on temperature close to an infinite vertical flat plate. The AB fractional operator is introduced in the governing equations of temperature and diffusion by employing the generalized types of Fourier and Fick's law. The developed nondimensional fractional model is solved using the Laplace transform method. Graphical illustrations are used to depict the physical behavior of fractional derivatives and the consequence of diverse flow parameters on velocity, thermal, and concentration fields. Furthermore, for validation of our attained results, some limiting cases are considered to recover fractional derivatives, as well as classical models of Maxwell and Newtonian fluids. The impacts of diverse flow parameters on variable profiles are achieved and presented graphically with significant conclusions.

Mathematical formulation based on a hybrid nanofluid
Consider an unsteady and an incompressible Oldroyd-B HNF flow close to an infinite vertical flat plate. Initially, consider that the fluid and plate are at a relaxation position, with constant temperature T ∞ and concentration C ∞ . After some time, the plate is kept constant and the fluid begins to move with a temperature value T(0, t) − a 1 zT(0,t) zξ u 0 sin ωt, where u 0 is a constant that signifies the dimension of velocity. At that time, the plate obtains a temperature T w and concentration C w , which persist constantly. We supposed that velocity, temperature, and concentration profiles are the only functions of ξ and t. The configuration of the problem is shown in Figure 1.
By Boussinesq's estimation (Ali et al., 2021), the governing equations for an Oldroyd-B HNF are discussed by Martyushev and Sheremet (2012). The equation of motion is as follows: The energy balance equation is as follows (Awan et al., 2019): The Fourier law (Zhang et al., 2022) for thermal conduction is as follows: The diffusion equation (Awan et al., 2019) for The Fick law is as follows (Awan et al., 2019): The appropriate initial and boundary conditions are as follows: Table 1 shows the properties of thermal and under-conversation fluids and nanoparticles. Frontiers in Materials frontiersin.org The following are a set of non-dimensional parameters: By utilizing the aforementioned variables in Eqs. 1-8 and after dropping the * notation, we obtain where The properties of a HNF are defined by Zhang et al. (2022).

Fractional model based on a non-local kernel
Now, we develop a fractional Oldroyd-B HNF using Fourier and Fick's law based on the AB-fractional operator (Atangana and Dumitru, 2016), which is explained as the following expression for a function f(ξ, t) and the kernel Mittage-Leffler function E γ (τ) is defined by The Laplace transform is The governing equations for the AB-fractional derivative are obtained by substituting the ordinary derivative with the AB derivative operator AB D α η in Eqs. 11-15 as q ψ, η −Ω 6 zθ ψ, η zψ .
3 Solution of the problem

Momentum profile
Taking the Laplace transform on Eq. (24) with related conditions in Eqs. (15) 1 -(17) 1 , we have By using temperature values from Eq. (37) and concentration from Eq. (38) and with conditions of Eq. (40), we obtain the solution of the velocity field for Eq. (40) as  Simulation to explain the temperature for changing α when Pr 3 and ϕ 1 ϕ 2 0.4.

FIGURE 3
Simulation to explain the temperature for Pr when α 0.4.

Frontiers in Materials
frontiersin.org solution of dimensionless fractional equations of energy, concentration, and momentum is obtained with the Laplace method. To observe from the physical perception, the impacts of fractional derivatives and different flow parameters on concentration, velocity, and temperature are measured and shown in Figures 2-15 graphically. Figure 2 shows the influence of α on the temperature field. By setting other parameters constant and fluctuating the value of α, it is seen that for a small time, the temperature profile declined for larger values α and this effect is reversed for a greater time. We see that fluid characteristics can be measured by fractional parameters. For a different value of α, the temperature close to the plate is extreme. The temperature declines away from the plate and is asymptotic in the growing ξ direction, which satisfies our boundary conditions. Figure 3 shows the thermal behavior for Pr. For large estimations of Pr, the temperature declines. Substantially, the heat conductivity increasing the estimations of Pr, manufacturing the fluid thicker, sources the least thickness of the heat boundary layer. Figures 4, 5 show the temperature behavior with ϕ 1 and ϕ 2 . The temperature field represents an increasing function of ϕ 1 and ϕ 2 . As expected, with greater values of ϕ 1 and ϕ 2 , the capacity of the HNF expands to hold additional heat. Therefore, the heat conductivity of the NF increases and temperature increases at different times.
The fluid velocity declines as we increase α, as shown in Figure 6, when there is less time. For a long time, the velocity is enhanced. Physically, when α increases, the velocity and thermal boundary layer decline, and as a consequence, the velocity declines for a short time. Figure 7 shows the behavior of the

Frontiers in Materials
frontiersin.org velocity with Pr. The velocity field also decreases with increasing Pr. Enhancement in Pr decreases the thermal conductivity and increases the viscosity of the fluid because of which the momentum profile declines with Pr. Figure 8 shows the influence of Gr on the momentum profile. By increasing Gr, the velocity profile is enhanced. Since Gr exhibits the buoyancy force that increases the natural convection, therefore the velocity grows. Figure 9 shows the impact of Gm on the velocity by considering the changing Gm with time. The ratio of the buoyant force and viscous force is named the mass Grashof number that sources unrestricted convection. Figure 9 shows that velocity is enhanced for enhancing Gm. Figures 10, 11 show the effect of ϕ 1 and ϕ 2 on velocity.
The velocity decreases with increase in ϕ 1 and ϕ 2 . This means that with the addition of nanoparticles to the base liquids, the resulting HNF becomes denser, so they become more viscous than the regular fluid. Also, the boundary layer of regular fluids is thinner than that of the HNF, and as a result, the velocity shows a declining behavior with increasing values of ϕ 1 and ϕ 2 . Moreover, the impact of a water-based HNF has more progressive values as compared to that of the sodium alginatebased HNF on the profiles of energy and velocity. Figure 12 shows a comparison of different fluid models. It is observed that the solutions of Maxwell nanofluids for both ordinary and fractional cases have developed curves as compared to Oldroyd-B and viscous nanofluids. Figure 13 shows the velocity for the slip

Frontiers in Materials
frontiersin.org and no-slip conditions. It can be seen that the slip condition shows a lesser profile for velocity than the no-slip conditions. Figure 14 shows the temperature and velocity behaviors for the comparison of

FIGURE 14
Simulation to explain the concentration, temperature, and velocity for the comparison of different inversion numerical algorithms.

FIGURE 15
Comparison of our results with the results by Chen et al. (2022) for validation.  (Stehfest and Tzou's algorithm). The overlapping of profiles shows that these algorithms are strongly validated with each other. Figure 15 shows the validation of our results with Chen et al. (2022). By overlapping both curves, it is observed from these graphs that our achieved results match those developed by Chen et al. (2022). The numerical comparison of energy, concentration, and velocity profiles by different numerical methods is shown in Table 2. Table 3 shows the numerical results of the Nusselt number, Sherwood number, and skin friction. The comparison of the momentum profile with the work of Chen et al. (2022) is shown in Table 4.

Conclusion
This article examines the investigations of the unsteady, convective flow of the Oldroyd-B HNF flowing over a flat plate with wall slip conditions on temperature and constant concentration. The model is developed using the AB-fractional operator and solved with the Laplace transform method. The Laplace inversion is computed with the wellknown Stehfest and Tzou numerical schemes. Finally, the effect of diverse flow parameters is planned to estimate the physical clarification of the achieved results of governed equations. The main results from the previous section are summarized in the following: ❖ For a short time, the temperature and momentum profile decayed for a larger value of α, and this effect for both profiles is reversed for a longer time. ❖ By increasing Pr, the temperature and velocity show a decreasing behavior. ❖ By increasing Gm and Gr, the velocity profile is improved. ❖ The velocity decreases with increasing ϕ 1 and ϕ 2 . ❖ The energy and velocity profiles are larger for a water-based HNF than those of the sodium alginate-based HNF. ❖ The graphs of Maxwell nanofluids for both classical and fractional models have more advanced curves than Oldroyd-B and viscous nanofluids.
❖ The slip condition shows a lower profile for velocity than the no-slip condition. ❖ The comparison of diverse numerical algorithms (Stehfest and Tzou) strongly validated our study's solutions. ❖ Chen et al. (2022), the overlapping of both curves validate the achieved results of our study.

Future recommendation
For extension of this fractional problem examined in this article, we idolized the following proposal based on investigation, approaches, extensions, and geometries, as demarcated in the following: • The same problem can also be considered over a horizontal plate by using Prabhakar's time-fractional approach with an MHD effect in a porous medium. • A comparative study of this study can be solved by the natural and Laplace transform methods. • The same problem may be discussed by the Keller-box scheme.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.