Application of water based drilling clay-nanoparticles in heat transfer of fractional Maxwell fluid over an infinite flat surface

In the present paper, unsteady free convection flow of Maxwell fluid containing clay-nanoparticles is investigated. These particles are hanging in water, engine oil and kerosene. The values for nanofluids based on the Maxwell-Garnett and Brinkman models for effective thermal conductivity and viscosity are calculated numerically. The integer order governing equations are being extended to the novel non-integer order fractional derivative. Analytical solutions of temperature and velocity for Maxwell fluid are build using Laplace transform technique and expressed in such a way that they clearly satisfied the boundary conditions. To see the impact of different flow parameters on the velocity, we have drawn some graphs. As a result, we have seen that the fractional model is superior in narrate the decay property of field variables. Some limiting solutions are obtained and compared with the latest existing literature. Moreover, significant results can be observed for clay nanoparticles with different base fluids.

Since the last decades, in various disciplines, fractional calculus theory has brought greater attention for the researchers. In fact, it was found that use of fractional derivatives is very helpful in modifying many process related to thermal transport processes, engineering sciences, circuit analysis, Biotechnology and signal processing. There are so many other applications related to heat and mass transfer and fluid dynamics can be perceived in the references [1][2][3][4][5] . Different books related to fractional derivatives were written by Jagdev et al. 6 , Kolade and Atangana 7 and Baleanu et al. 8 and they also discuss the application of fractional derivatives and also about their operators.
In literature different concepts related to fractional derivatives can be find out. In physics fractional derivative plays an important role in designing different phenomena. Yet, it was mentioned in Caputo and Fabrizio 9 , that different circumstances belonging to material heterogeneities cannot be well-modeled using fractional derivatives introduced by Riemann-Liouville or Caputo. Due to this fact, new fractional derivative related to non-singular kernel were introduced by Caputo and Fabrizio 9 . It was important that, rather than power law and exponential decay function the kernel Mittag-Leffler function is more general. Therefore, both Riemann-Liouville and Caputo-Fabrizio are special cases of Atangana-Baleanu fractional operators [10][11][12][13] . In applied sciences, some new kinds of derivatives exist which are known as fractal derivative. Therefore there was a need to redefined the concept of expressing velocity in fractal media for example, scaling time in fractal (x, t α ) . Baleanu et al. 14 in 2020 introduced a new fractional operator using power law and are called hybrid fractional derivatives. This derivative is linear combination of constant proportional and Caputo type fractional derivative.
Ali 15 investigated the Atangana-Baleanu derivative with a novel approach. Circuits with fractional derivatives were developed by Hammouch and Mekkaoui 16 and they also discussed its behavioral dynamics. It is known that for description of price of opinion can be given by using Time Fractional Black Scholes Equation (TFBSE) with a time derivative of real order. The investigation of heat dissipation in transmission line of electrical circuit is given in Abro et al. 17 . An analysis of generalized Jeffery nanofluid in a rotating frame with non-singular fractional derivative is given in Ali et al. 18 . The behavior realted to heat transfer in different model with singular and non-singular is given in articles [19][20][21][22] . To control the entropy in the flow of heat is one of the main concerns of the industrial sectors. Since the problem had a large amplitude so a lot of researchers pay attention to solve it. They use different methods for different fluids to enhance the thermal conductivity 23 . Bejan et al. has came to the conclusion that viscous dissipation, heat transfer, mass transfer and chemical reaction are the main reasons are entropy enhancement in thermal systems [24][25][26][27] . To overcome the problem the first meaningful contribution was made by Choi, when he gave the idea of nanofluids 28 . This was the revolution in the flow of heat and has given the answers to the many unsolved problems. He just added nanosized particles of different type to solve the problem of entropy. Entropy in magnetohydrodynamics (MHD) by having exact analysis is investigated by Khan et al. 29 . A definition regarding the Bejan number which is useful to predict the power of magnetic field and entropy of fluid friction due to heat transfer is given by Awed 30,31 . Saouli and Aiboud-Saouli 32 has used an inclined plate and liquid film to analyzed entropy generation. Mahmud et al. 33 reported the experiment with added magnetic field influence. Selimefendigil et al. 34 has explored entropy for natural convection flow of a nanofluid. 35-39 made a useful contribution in this regard. 40,41 has discussed that how can nanofluids be used in solar energy system and evaporation. Brownian motion and thermophoresis effect in heat transfer was analyzed by Buongiorno 42 . Some valuable studies have been done in this regards by Biglarian et al. 43  The two-dimensional magnetohydrodynamics (MHD) flow of Casson fluid and heat transfer find out by Hamid et al. 53 . Mainly we focus the study is to survey the linear thermal radiation influence on dual solutions for both the steady and unsteady flow of Casson fluid under the effect of uniform magnetic field. The blood flow connecting nanoparticles through porous blood vessels in the occurence of magnetic field with the help of collocation and least squares techniques inspected by Usman et al. 54 . Blood is a non-Newtonian fluid having nanoparticles which are used for different models to find the viscosity of the nanofluids. Hamid et al. 55 disscussed the unsteady MHD flow of Williamson nanofluid between the permeable channel with heat source/sink. The effect of molybdenum disulfide ( MoS 2 ) nanoparticles forms on circling flow of nanofluid along an elastic stretched sheet. This nanofluid flow is measured in the existance of magnetic things, thermal radiation and variable thermal conductivity. Hamid et al. 56 discussed the different types of nanoparticles like Platelet, cylindrical and brick forms. Usman et al. 57 evaluate the flow of ethylene glycol and water based copper (Cu) nanoparticles between two squeezed parallel disks.
Newly, the industrialists are interested theoretically and experimentally studying the advancement of nanofluids. The industrialists are elaborated to find thermophysical properties (heat capacitance, thermal conductivity, density, electrical conductivity, thermal expansion, and viscosity) of unlike nanoparticles and base fluids using several procedures 58 because the next generations fluid is a nanofluid for heat transport which can deal additional thermal presentation in various industrial sectors like power generation, transportation, hyperthermia and air conditioning.
Inspiriting from the above-discussed literature, this study aims to focus on the applications of nanofluid in the drilling process. Different type of nanofluid is used in drilling activities such as oil-based drilling mud (OBM), water-based drilling mud. For this purpose we have used the Polymers nanoparticles and clay nanoparticles. It has been observed that for stability in temperature, least value of torque, prevention of fluid loss, to control the rheological possessions for scrubbing the hole and to filter the quality of cake, clay nanoparticles helps a lot. Khan et al. 59 has used three different base fluids to clean water. Nisar et al. 60 has studied the entropy of clay nanoparticles.
The main purpose of this paper is to extend the idea of Imran et al. 61 in which analytical solutions are obtained for viscous fluid. They used the Laplace transform method to obtain the solutions for temperature and velocity fields respectively. For the moment there is no such results regarding Maxwell fluid containing clay nanoparticles therefore, we have applied the most recent hybrid fractional operator for a Maxwell fluid of caly-water base nanofluids over an infinite vertical surface moving with constant velocity and obtained solutions with Laplace transform method. Some limiting solutions are also obtained and justified through graphical comparison and presented in the graphical section. Thermophysical properties of nanomaterials are defined in Table 1.

Mathematical formulation and solution
Let water, kerosene and engine oils are base fluids to carve-up in the flow of clay nanoparticles. The flow of the fluid is in the region y 1 > 0 , close to a heated flat vertical plate. The plate is normal to y-axis and is fixed. In the beginning it is assumed that the fluid is at rest on the plate having surrounding temperature T ∞ . This ambient temperature changes from T ∞ to T w in no time causing the motion in the plate with velocity V 0 forcing the fluid to move in x-direction as shown in Fig. 1. The governing equations are given as follows [62][63][64][65] . Table 1. Thermophysical properties of nanofluids 59 .

Material
Symbol where -p is the pressure, I is the identity matrix, µ is the viscosity, and A 1 is the first Rivlin-Eriksen tensor defined as where D Dt is the material time derivative, and L is gradient of the velocity. Fractional stress tensor for Maxwell fluid with constant proportional Caputo time fractional derivative 14 where τ = S xy is the non zero component of extra stress tensor and ∂α ∂t α is the constant proportional Caputo derivative of non integer order α . The innovative fractional derivative is defined and given in 14 .

Numerical outcomes and analysis
This paper deals with the investigation of Clay nanoparticles in free convection of a Maxwell fluid. The analytical solutions satisfies the initial and boundary conditions. The solutions are obtained with the application of novel fractional derivative and Laplace transformation. The influence of nanoparticles as well fractional parameter are discussed through some graphs.
(−y 1 √ c 1 ) n 1 (c 2 ) n 2 ( ) n 3 (K 1 (α)) n 4 t 1− n 1 2 −αn 3 +n 4 n 1 !n 3 !n 4 !(b 0 ) n 2 (K 0 (α)) n 4 −n 2 −n 3 Ŵ(n 2 + n 3 + 1)Ŵ( n 1 2 − n 2 + 1) Ŵ(n 2 + n 3 + 1 − n 4 )Ŵ( n 1 2 − n 2 + 1 − n 3 )Ŵ(2 − n 1 2 − αn 3 + n 4 ) .     Fig. 3. Further, fractional parameter can be used to control the momentum boundary layer thickness. The relationship of velocity of clay nanoparticles and volume fraction φ is discussed in Fig. 4. For large values of φ fluid velocity will decrease as Table 2. Effect of fractional parameter on dimensionless velocity.     www.nature.com/scientificreports/ viscous forces became stronger with increasing φ . Figure 5 high lights the influence of Grashof number Gr . Fluid velocity increases as we enhance the value of Gr . Gr describes the influence of the thermal buoyancy force to the viscous force. If Gr equal to zero, then there is no free convection current, if Gr is greater than zero, plate is outwardly chilled and if Gr is less than zero, plate is outside frenzied. It means larger plate is outwardly chilled increasing the velocity and effect is reverse of Pr . Figure 6 represents the contrast of velocity profile for three different types of base fluids (water, kerosene, and engine oil). It is observed that the velocity of water-based clay nanofluid fluid is larger than kerosene oil and engine oil-based clay nanofluids, respectively. Because the  www.nature.com/scientificreports/ thermal conductivity of water is larger than that of kerosene oil and engine oil then the velocity of water-based clay nanofluid is larger than the others. Figure 7 our obtained results by taking = 0 are compared with the results of Imran et al. 61 while keeping other parametric values constant. Table 2 represents the velocity comparison of present paper when = 0 and Imran et al. 61 for different values of parameter α . We have seen that both Fig. 7 and Table 2 are in good agreement with each other. Figure 8 is plotted to see the velocity comparison of hybrid fractional derivative and Khan et al. 59 . Since the velocity obtained in 59 is for viscous fluid and in the present study for Maxwell fluid. This figure shows the viscous fluid is swiftest than Maxwell fluid. The reason is that Maxwell fluid is non-Newtonian one and more thicker than viscous. Table 3 shows the velocity comparison of current paper and Khan et al. 59 for different values of parameter α . In both cases we have found that velocity is smaller for CPC fractional model. Physically, fractional operator is responsible for the history of the model and can have better control for momentum boundary layer. Figure 9 shows the velocity comparison between present result and Danish et al. 69 with the solutions obtained with the same fractional operator for viscous fluid and we see that velocity is smaller for constant proportional Caputo model of Maxwell fluid over viscous fluid. Due to less viscosity of viscous fluid, the it flows faster than Maxwell fluid. The velocity contrast of present result and Danish et al. 69 is shown in Table 4 and velocity is minimum for constant proportional Caputo model. The velocity comparison between present result when =0, α = 1 , Imran et al. 61 when α = 1 and Khan et al. 59 is shown in Fig. 10. We see that these results are in good agreement. Table 5 represents the velocity contrast of present result and published results and we see that results are same. The influence of fractional parameter β on Nusselt number is studied numerically in Table 6. Nusselt number is decreasing function of fractional parameter β . The influence of fractional parameter α on Skin friction is evaluated numerically in Table 7. Skin friction is decreasing function of fractional parameter α.