Blasius–Rayleigh–Stokes Flow of Hybrid Nanomaterial Liquid Past a Stretching Surface with Generalized Fourier’s and Fick’s Law

The effect of Stefan blowing on the Cattaneo–Christov characteristics of the Blasius–Rayleigh–Stokes flow of self-motive Ag-MgO/water hybrid nanofluids, with convective boundary conditions and a microorganism density, are examined in this study. Further, the impact of the transitive magnetic field, ablation/accretion, melting heat, and viscous dissipation effects are also discussed. By performing appropriate transformations, the mathematical models are turned into a couple of self-similarity equations. The bvp4c approach is used to solve the modified similarity equations numerically. The fluid flow, microorganism density, energy, and mass transfer features are investigated for dissimilar values of different variables including magnetic parameter, volume fraction parameter, Stefan blowing parameter, thermal and concentration Biot number, Eckert number, thermal and concentration relaxation parameter, bio-convection Lewis parameter, and Peclet number, to obtain a better understanding of the problem. The liquid velocity is improved for higher values of the volume fraction parameter and magnetic characteristic, due to the retardation effect. Further, a higher value of the Stefan blowing parameter improves the liquid momentum and velocity boundary layer thickness.


Introduction
At present, researchers are keen to study heat transfer applications involving nanofluids. Based on empirical results in various technical and medical fields, the heat transport, mass, and density during the flow have dynamic features. Many base liquids are not favorable for heat transfer applications due to poor thermal effectiveness. To overcome this problem and improve the heat transfer efficiency, we added nanoparticles to the base fluids. The nano-size particles are most effective in enhancing the heat transfer rate. These nanoparticles consist mostly of metals such as silver, steel, gold, copper, potassium, nitrides, and many others. Choi [1] was the first to introduce the theory of nanofluids. Later, Buongiorno [2] discussed the concept of natural convective heat transport in a nanofluid. Sheikholeslami and Chamkha [3] explored the Lorentz force effect on the nanofluid flow. Ahmad et al. [4] introduced the influence of carbon nanotube nanoparticles on the boundary layer flow with thermal radiation, double stratification, and heat generation. They found that the thermal stratification diminished the temperature distribution. The natural convection heat transfer and inclined magnetic field effects of a molybdenum disulfide (MoS2) and magnesium oxide (MgO) nanofluid were investigated by Hymavathi et al. [5] via a vertical stretched surface embedded in a porous medium with a non-uniform heat source or sink. Nadeem et al. [6] used a magnetic dipole with three different ferrite nanoparticles to assess the effects of heat transport phenomena in a ferrofluid. They discovered that the characteristic of magneto-thermomechanical cooperation reduced the movement of liquid molecules, and as a result, increased the coefficient of skin friction and the thermal transport rate at the surface. Together with the impacts of viscous dissipation and non-linear thermal radiative effects, Reddy et al. [7] investigated the role of gyrotactic microorganisms in the mass and heat transport for the time-dependent MHD flow of a cross liquid mediated through nanoparticles. Vajravelu et al. [8] evaluated the influence of the nanoparticle volume fraction on the flow and heat transfer characteristics of Ag-water and Cu-water nanofluids under the impact of internal heat absorption or generation and thermal buoyancy. Jamshed et al. [9] explored the Casson non-Newtonian Cu-water and TiO 2 -water nanofluid flows in terms of heat transport and entropy generation under the influence of solar thermal transfer and slip conditions. Many researchers (see [10][11][12]) have recently observed an improvement in the heat transport rate in the presence of various nanoparticles.
Many investigators have studied nanofluids to date, but studies related to hybrid nanofluids are attracting the attention of numerous scientists due to a wider range of heat conductivity applications. The combination of two or more nanoparticles in the base liquid generates a hybrid nanomaterial liquid. The reason that researchers are giving attention to this issue is the heat transport augmentation that can be accomplished by these hybrid nanofluids at a low cost of production. Mingzheng et al. [13] demonstrated the viscosity and thermal conductivity of various kinds of surfactant mixtures. Ahmad et al. [14] explored the thermal transport characteristics of a hybrid nanomaterial liquid with activation energy on a wedge, with the availability of variable velocity, Darcy-Forchheimer flow, and thermal slip. The influence of suction and buoyancy force effects on a Ag-MgO/water hybrid nanofluid flow across a stretching surface was introduced by Anuar et al. [15]. Maskeen et al. [16] explored the thermal transport properties of magneto-hydrodynamic hybrid aluminacopper/water nanomaterial liquid flow on an extending cylinder, with the influence of radiation and Joule heating. Ghalambaz et al. [17] demonstrated the thermal transport enhancement in a hybrid (Ag + MgO) nanofluid past a porous square cavity, using the local thermal non-equilibrium (LTNE) model. Manna et al. [18] proposed a new multibanded magnetic field approach for improving the controllability of convective transports. Applying four equal parts of a magnetic field horizontally over a heated system occupied by porous material saturated with Cu-Al 2 O 3 /water hybrid nano liquid demonstrates the multi-banding approach. Ahmad and Nadeem [19] examined the properties of mass and heat transport in the presence of hybrid SWCNT-MWCNT/water nanomaterial liquid with Hall slip, ion slip, and chemical reaction impacts. Esfe et al. [20] scrutinized the impact of the volume fraction of nanoparticles on the dynamic viscosity and thermal conductivity of Ag-MgO/water hybrid nanofluids with particle sizes of 25(Ag) and 40(MgO) nm and nanofluid volume fractions (50 percent Ag and 50 percent MgO by volume) ranging from 0 to 2 percent and found new interrelations. In a modified Buongiorno's model for the nanomaterial liquid magneto-transport phenomenon over an expanding cylinder in the vicinity of motile microorganism density, Rana et al. [21] integrated the Cattaneo-Christov mass flux (non-Fick's) and heat (non-Fourier's) ideas. Ma et al. [22] provided a (2D) numerical simulation of Ag-MgO nanomaterial forced convection and thermal transport in a channel with active coolers and heaters, to investigate the effect of a magnetic field on heat transport and nanomaterial liquid Ag-MgO forced convection.
Mass and heat transport phenomena occur due to the temperature and mass gradient respectively. Heat and mass transport have many industrial and engineering applications such as in heat exchangers, nuclear plants, heat pumps, heat conduction in tissues, refrigeration, heat transfer through materials, the diffusion of chemical impurities in rivers and oceans, the evaporation of water, etc. The convectional laws for the analysis of mass and heat transport used were Fourier's [23] and Fick's [24] laws. Fourier's and Fick's laws have a disadvantage; they give a parabolic type of energy and mass equation. To overcome this drawback, Cattaneo [25] modified Fick's and Fourier's laws by the addition of a time factor. Later, Christov [26] altered the Cattaneo law by replacing the time factor with the Oldroyd-B upper convective derivative. Han et al. [27] investigated the Maxwell liquid flow on a stretchable surface to present an evaluation of the Fourier heat flux and the Cattaneo-Christov model. With the addition of chemical reactions, a uniform heat source/sink, and thermal radiation, Venkateswarlu et al. [28] explored the impacts of a magnetic field on the flow of the Cattaneo-Christov heat flux model for MoS2 and MgO water-based nanofluids across a stretching sheet. Ali et al. [29] introduced the impact of the Cattaneo-Christov characteristics and bioconvection on self-motivated microorganisms in water-based nanoparticles with a leading-edge accretion/ablation and the Stefan blowing effect. The heat transport with respect to the Cattaneo-Christov concept with variable thermal relaxation time past a stretching sheet was discussed by Ahmad et al. [30]. Raju et al. [31] illustrated the heat, flow, and transport of mass attributes of a Maxwell nanofluid through a cylinder with a heat sink/source, using the Cattaneo-Christov premise. The Cattaneo-Christov premise was utilized by Malik et al. [32] to investigate the performance of Sisko fluid via a porous non-linearly stretched cylinder. Using a simplified mathematical model published by Jamshed and Aziz [33], the entropy production and heat transfer analysis of a thermal system containing hybrid nanomaterial with a Cattaneo-Christov flow model and thermal radiation impacts was explored. Garia et al. [34] used an extending surface with Joule heating and thermal radiation to solve the magneto-hydrodynamic flow of a SiO 2 −MoS 2 /water hybrid nanofluid past a wedge and a cone, with a generalized Fourier's model.
The transport of electric current through any conducting material produces Joule heating. The collision of moving particles is the reason behind this. As a result of this process, some kinetic energy is transformed to heat, and the temperature of the material rises. Scientists and engineers have been fascinated by the idea of improving the efficiency of numerous mechanical systems and industrial machinery in recent years. Such problems can be overcome by lowering the temperature induced by Ohmic dissipation or Joule heating. As a result, several academics studied the flow issues from diverse physical perspectives. Sahoo [35] explored the impacts of partial slip and Joule heating on the flow of second-grade MHD liquid on a stretched surface with thermal transmission. Shehzad et al. [36] reported the radiative MHD extended flow of a Jeffrey liquid with Joule heating. Waqas et al. [37] discussed the Carreau-Yasuda nanomaterial liquid flow past an extending/shrinking surface in the presence of Joule heating, motile microorganisms, and thermal radiation, under Robin's conditions. In the presence of thermal radiation and Joule heating, Kumar et al. [38] evaluated the 3D Oldroyd-B nanofluid flow across a stretching sheet. Hayat et al. [39] investigated the 3D steady second-grade nanomaterial flow via a rotating disc, with Joule heating and heat generation/absorption. For the constant wall heat flux condition, Chakraborty et al. [40] critically evaluated the transport-of-heat performance related to thermally developed coupled-electron magneto-hydrodynamic flows towards narrow flow conduits, taking electrokinetic impacts into account. Chen [41] conducted a study to investigate the impacts of Joule heating and viscous dissipation on the thermal transport and momentum, for magneto-hydrodynamic (MHD) flow towards an extending surface in the presence of radiation and surface suction/blowing. Technical engineers, designers, manufacturers, and developers are keen to study Blasius-Rayleigh-Stokes flow. To transform the time-dependent flow model's partial differential equations into a similar form, a new type of transformation was developed by Na [42], which consists of the Blasius-Rayleigh-Stokes variables. Fang [43] analyzed the heat transportation in an unsteady flow with the implementation of Blasius-Rayleigh-Stokes variables on a flat surface. Todd [44] was the first to establish the concept of leading-edge accretion/ablation, which is based on the unstable boundary layer model. The influence of radiation, Joule heating, and the Cattaneo-Christov theory applied to a Blasius-Rayleigh-Stokes flow towards a transitive magnetic field was discussed by Reddy et al. [45]. Several scholars have published studies on the Blasius-Rayleigh-Stokes variables in recent decades [46,47].
The unsteady MHD stagnation point Blasius-Rayleigh-Stokes hybrid nanomaterial liquid flow with Cattaneo-Christov theory, Joule heating, and convective conditions effect over an extending surface were investigated. To the best of our knowledge, no study has yet investigated the Blasius-Rayleigh-Stokes flow of hybrid nanomaterial liquid with convective boundaries. The Blasius-Rayleigh-Stokes variables were applied to transport the flow model into a couple of ODEs. The bvp4c MATLAB approach [48,49] was applied to solve these multiple equations numerically. The evolving parameters are discussed graphically along with the velocity, microorganism density, temperature, and concentration distribution. Further, a tabulated analysis of physical quantities is presented. To validate the problem, a comparison table is presented.

Mathematical Modeling
We consider the unsteady, laminar, incompressible, viscous, 2D, MHD boundary layer of the stagnation point Blasius-Rayleigh-Stokes flow of a hybrid nanofluid in the presence of viscous dissipation and microorganism density. The thermal and solutal energy transport analyses are presented with the influence of the Cattaneo-Christov theory and Joule heating. Furthermore, convective boundary conditions are applied on the surface of the stretching sheet and free-stream conditions are considered away from the surface. The transitive magnetic field is applied normally to the flow field. The free-stream and fluid velocities are U ∞ and U w , respectively, as shown in Figure 1. The fluid microorganism density, concentration, and temperature are N, C, and T, respectively. The wall temperature is T w , the wall concentration is C w , and the wall microorganism density is N w ; away from the wall these terms are T ∞ , C ∞ , and N ∞ , respectively. Using the boundary layer approximation of and the above supposition, the equations of momentum, mass, energy, concentration, and microorganism density become [11,21,45]: ∂v ∂y ∂N ∂t and the related boundary conditions are [29]: Figure 1. Problem configuration.
In the above equations, σ in f is the electrical conductivity of the nanofluid, the density of the nanomaterials liquid is ρ in f , the kinematic viscosity of the nanomaterials liquid is υ in f , the magnetic field is B, the thermal diffusion coefficient is α * in f , the microorganism concentration and diffusivity coefficient are D m and D n , respectively, h f , h g , and h n are the heat, mass, and microorganism transport coefficient, respectively, λ 1 and λ 2 are the thermal and concentration time relaxations, respectively, W c signifies the cell swimming speed, and b symbolizes the chemotaxis constant.

Similarity Analysis
Using the similarity transformation [29,44,45]: Applying Equation (7), the above nanoparticle phase equation transmutes to: g + PrL e 2 ζ cos(w * )g + f sin(w * )g + γ c sin(w * ) 3 f f g + ζ f g = 0, The transformed conditions are: The parameters ablation/accretion, thermal Biot number, Lewis number, concentration Biot number, microorganism Biot number, magnetic field parameter, bio-convection Lewis number, Stefan blowing parameter, Eckert number, thermal relaxation parameter, concentration relaxation parameter, Peclet number, and solid volume fraction of particles are symbolized by w * , B e , L e , B c , B n , M, L b , s, E c , γ e , γ c , P w and φ 1 , respectively. These parameters are defined as: The hypothetical relation is characterized as follows [20][21][22].

Quantities of Interest
The local friction drag and the local density number of motile microorganisms are: Using Equation (10), the dimensionless form is, The local Reynolds number is characterized as

Results and Discussion
In this section, the physical implications of prominent characteristics such as the magnetic characteristic (0 ≤ M ≤ 3), Stefan blowing parameter (0 ≤ s ≤ 6), solid volume fraction (0 ≤ φ 1 ≤ 0.03), thermal Biot number (0 ≤ B e ≤ 1.5), Eckert number , and microorganism Biot number (1 ≤ B n ≤ 3) against velocity, temperature, the volumetric concentration of nanoparticles, and motile microorganism number is reflected in Figures 3-16. Table 3 shows that higher values of φ 1 increase the skin friction, whereas higher estimates of M, α and s reduce the skin-friction coefficient. The effect of distinct values of L b , P e , B n and ω * along with the microorganism number is presented in Table 4. It was discovered that higher estimates of L b , P e and ω * increase the motile density transfer rate; however, higher values of B n have the reverse effect. Figure 2 describes the flowchart of bvp4c scheme. The influence of the magnetic parameter (M), volume fraction parameter (φ 1 ), and Stefan blowing parameter (s) on the liquid velocity is indicated in Figures 3-5. It can be noted from Figure 3 that a stronger estimate of M increases the velocity profile. The cause of this pattern is that a magnetic field in an electrically conducting liquid induces a Lorentz force and reduces the momentum boundary layer flow considerably. Figure 4 displays the influence of distinct values of the volume particle parameters versus velocity distribution, and we can perceive that the impact of solid volume fraction enhances thermal transmittance in both solutions, due to the increased number of nanoparticles in the fluid system, which increases resistance while also improving boundary layer thicknesses. Figure 5 displays the effect of the Stefan blowing characteristic on the fluid velocity. It can be noted that the fluid velocity is increased due to a higher estimate of s. Physically, buoyancy forces emerge as a result of the blowing effect, which intensify fluid movement and transport heated nanofluid away from the surface. The influence of the volume fraction parameter (φ 1 ), thermal Biot number (B e ), Eckert number (E c ), and thermal relaxation characteristic (γ e ) on the thermal field is presented in Figures 6-9. It can be seen from Figure 6 that larger values of φ 1 increase the fluid temperature because, due to a larger φ 1 , a retardation effect takes place, which boosts the temperature. The effect of the Biot number on the temperature is described in Figure 7. The Biot number, in reality, includes the heat transfer coefficient. As the Biot number increases, so does the heat transmission coefficient. The fluid temperature and the related boundary layer thickness increase as the heat transport coefficient increases. Figure 8 depicts the features of E c in the temperature distribution. High Eckert numbers cause the fluid temperature to rise faster because the fluid particles move faster, allowing them to collide more frequently and produce heat, causing the fluid temperature to rise faster. Figure 9 shows that an increase in the value of the thermal relaxation characteristic reduces the temperature. Physically, this is because the time it takes for material particles to transfer heat to adjacent particles reduces when γ e is increased. In other words, greater values of γ e cause the material to behave in a non-conducting manner, resulting in a drop in fluid temperature. Furthermore, it should be highlighted that for γ e = 0, heat is quickly transferred throughout the material. The temperature profile is greater for γ e = 0, i.e., in the Fourier's law case, than for the Cattaneo-Christov premise. The influence of the Stefan blowing parameter (s), concentration Biot number (B c ), Lewis number (L e ), and relaxation concentration characteristic (γ c ) on the concentration is described in Figures 10-13. Figure 10 shows that a larger estimate of s reduces the concentration. The effect of B c on the concentration is shown in Figure 11. Fluid particles require a longer time to diffuse across the material medium as B c is increased. As a result, the concentration distribution is reduced. Figure 12 shows the varying behavior of γ c along the concentration distribution. It is drawn showing that a larger assessment of γ c reduces the concentration. Physically, more time is required to transport the nanoparticles from one place to another. The impact of the Lewis number on the concentration is seen in Figure 13. It is interesting to note that L e is the ratio between thermal and mass diffusivity; as we increase L e , the mass diffusivity increases, and as a consequence the concentration profile declines. The influence of the bio-convection Lewis parameter (L b ), Peclet number (P e ), and microorganism Biot number (B n ) on the microorganism density is shown in Figures 14-16. It can be seen from Figures 14 and 15 that larger values of L b and P e reduce the microorganism density profile. Physically, with an increment in L b and P e , the microorganism diffusivity declines, and as a result both curves are reduced. Figure 16 shows the behavior of B n against the microorganism density profile. It shows that the microorganism density profile is enhanced due toa larger estimation of B n .

Concluding Remarks
The time-dependent 2D Blasius-Rayleigh-Stokes flow of hybrid nanomaterial liquid with magnetic field and microorganism effects is discussed. Convective conditions are applied at the surface of the extending surface. The main results of the current research are as follows: 1.
The velocity of the liquid is increased for higher estimates of the volume fraction parameter and magnetic parameter, due to the retardation effect.

2.
Higher values of the Stefan blowing parameter improve the velocity of the liquid and the momentum boundary layer thickness.

3.
A decaying trend occurs due to a higher thermal relaxation characteristic because particles have extra time to transport heat to nearby particles.

4.
Kinetic energy is transformed into heat energy due to the enhancement of the Eckert number. 5.
The nanoparticle concentration declines due to larger estimates of concentration relaxation and Stefan blowing parameter. 6.
A higher estimate of thermal and concentration Biot number improves the heat and mass transfer rates, respectively. 7.
A decreasing behavior occurs in the microorganism density profile due to larger values of L b and P e . 8.
The motile density transfer rate decays for larger values of B n . Funding: The author T. Abdeljawad would like to thank Prince Sultan University for paying the APC and support through the TAS research laboratory. In addition, this project was funded by Taif University Researchers Supporting Project (number TURSP-2020/252), Taif University, Taif, Saudi Arabia.
Data Availability Statement: All data supporting this study are available in the article.

Acknowledgments:
The author T. Abdeljawad would like to thank Prince Sultan University for paying the APC and support through the TAS research laboratory. In addition, the authors would like to acknowledge financial support from the Taif University Researchers Supporting Project (number TURSP-2020/252), Taif University, Taif, Saudi Arabia.

Conflicts of Interest:
The authors declare no conflict of interest. Heat, mass, and microorganism transport coefficients, respectively Subscripts w

Abbreviations
The boundary surface ∞ The ambient surface hn f Hybrid nanofluid n f Nanofluid