Fractional Nadeem trigonometric non-Newtonian (NTNN) fluid model based on Caputo-Fabrizio fractional derivative with heated boundaries

The fractional operator of Caputo-Fabrizio has significant advantages in various physical flow problems due to the implementations in manufacturing and engineering fields such as viscoelastic damping in polymer, image processing, wave propagation, and dielectric polymerization. The current study has the main objective of implementation of Caputo-Fabrizio fractional derivative on the flow phenomenon and heat transfer mechanism of trigonometric non-Newtonian fluid. The time-dependent flow mechanism is assumed to be developed through a vertical infinite plate. The thermal radiation’s effects are incorporated into the analysis of heat transfer. With the help of mathematical formulations, the physical flow system is expressed. The governing equations of the flow system acquire the dimensionless form through the involvement of the dimensionless variables. The application of Caputo-Fabrizio derivative is implemented to achieve the fractional model of the dimensionless system. An exact solution of the fractional-based dimensionless system of the equations is acquired through the technique of the Laplace transform. Physical interpretation of temperature and velocity distributions relative to the pertinent parameters is visualized via graphs. The current study concludes that the velocity distributions exhibit an accelerating nature corresponding to the increasing order of the fractional operator. Moreover, the graphical results are more significant corresponding to the greater time period.

Temperature at the surface V Velocity vector Fractional calculus deals with the study of differentiation of fractional order.The real data through a fractional derivative approach can be estimated more effectively as compared to an ordinary derivative.Many fields incorporate the concept of derivatives with fractional order.Examples of such disciplines are fluid dynamics, electrochemistry, dynamic complex systems, and biological systems 1 .The most used fractional derivatives are the Riemann-Liouville derivative having non-integer order and the Caputo derivative.The Riemann-Liouville derivative includes the singular kernel with the non-zero derivative of a constant.The concept of the Caputo derivative has overcome the problem of the Riemann-Liouville derivative by utilizing the power law kernel, which is still singular, but a constant has zero derivative 2 .After that, the concept of Caputo-Fabrizio fractional derivative involves a non-singular kernel of exponential form.With the implementation of numerous models of fractional derivatives, various non-Newtonian fluid models have been examined.An exploration of the fractional derivative on a time-dependent chemically reactive flow behavior of a non-Newtonian fluid in one direction was investigated by Siddique and Bukhari 3 .The authors compared the flow characteristics regarding classical and fractional derivatives and identified significant outcomes in the case of the fractional approach.Antonio Taneco-Hernández et al. 4 discussed the time-dependent flow mechanism through a bar by implementing different fractional derivatives.The authors scrutinized the viscoelastic properties of the considered fluid in the context of the various fractional operators.Nadeem et al. 5 discussed the time-dependent flow of a trigonometric non-Newtonian fluid by considering two heated surfaces.The authors observed that an increase in the fluid parameter results in a reduction of the fluid flow mechanism.Recently many researchers discussed flow problems by adopting the technique of fractional derivatives [6][7][8][9][10][11] .
In industries and engineering fields, non-Newtonian fluids have various practices.The implementations of such fluids include spray coating, production of paper, lubricants, gaseous diffusions, polymer production, etc. Non-Newtonian fluids have complicated nature which can be examined through numerous mathematical models.With the help of such models, numerous investigations on the behavior of non-Newtonian fluids have been conducted.Abbas et al. 12 assumed a Riga surface to observe the steady flow nature of a non-Newtonian fluid with physical boundary conditions.The authors explained the characteristics of the flow phenomenon relative to the emerging parameters.Nadeem et al. 13 studied the radiative flow phenomenon of a non-Newtonian fluid in the context of a stretchable Riga surface.Amplification in the flow field related to improved fluid parameter was identified in their observations.Yang et al. 14 discussed the model of a non-Newtonian fluid with the involvement of nanoparticles and explored the hydromagnetic flow mechanism on both shrinking and stretching surfaces.The significance of an inclined magnetic field in the flow behavior of a non-Newtonian fluid through a vertical surface was demonstrated by Ishtiaq and Nadeem 15 .An augmentation in the velocity distribution was indicated in their observations corresponding to the improved fluid parameter.Under the influence of physical impacts, the non-Newtonian fluid model with variable characteristics was demonstrated by Ahmad et al. 16 .Through a Riga plate, Nadeem et al. 17 study a non-Newtonian fluid with its axisymmetric flow nature influenced by variable properties.The authors investigated the flow properties affected by the fluid parameter.The flow characteristics of a non-Newtonian fluid and the impacts of double diffusion theory on the mechanism of heat transport were deliberated by Irfan et al. 18 .The authors numerically and analytically explore the flow phenomenon in relation to significant parameters.The two-phase model-based thermal analysis of a fluid incorporating nanoparticles with the involvement of a movable wedge was inspected by Baby et al. 19 .The study's findings revealed that an improved wedge parameter led to a diminished fluid flow behavior.][22][23][24] The current study has a significant aspect of implementing the Caputo-Fabrizio fractional derivative on the unsteady flow of a non-Newtonian fluid with trigonometric form.No one examined the fractional model of the trigonometric non-Newtonian fluid flow until now.The novelty of the ongoing study includes the fractionalbased unsteady flow mechanism of trigonometric non-Newtonian fluid proceed by a vertical infinite plate.The energy equation is exhibited by incorporating the impact of the thermal radiation which is further analyzed with the momentum equation by adopting the definition of the fractional derivative.To acquire an exact solution to the concerned problem, an effective Laplace transform methodology is accomplished.The flow velocity and temperature influenced by numerous pertinent parameters are graphically explained.The purpose of the current study is to address the following queries.
a. What are the advantages of using a fractional derivative approach on non-Newtonian fluid flow instead of an ordinary derivative approach?b.How fractional derivative-based equations tackled through the Laplace transform technique?c.What is the impact of the fractional parameter on the boundary layer thickness?d.What effect does the augmented fractional parameter have on the mechanism of heat transport?

Research methodology
We consider the unsteady flow phenomenon of a non-Newtonian fluid with trigonometric type over a vertical infinite plate.In the setup of Cartesian coordinates, the static plate is placed toward the x-axis and in the perpendicular position of the plate, the y-axis is taken.Initially, both the plate and the considered fluid are in a static position with the fluid's uniform temperature T ∞ .As the plate is taken to be static, at time t = 0 + , the

Nadeem trigonometric non-Newtonian (NTNN) fluid model
The constitutive equation for the Nadeem trigonometric non-Newtonian (NTNN) fluid model is represented by the following expressions 5 , In Eq. ( 1), S exhibits the extra stress which is defined in the following way 25 Equation (2) converts into the following form after the first order expansion of sinh −1 1 2 γ ≈ 1 2 γ.
The expression of Eq. ( 3) from Eq. ( 2) can also be obtained if we use With the help of these trigonometric functions, the extra stress can be written as follows, By taking the first-order expansion of the above trigonometric functions, we have the following result (1) The expressions of extra stress tensor given in Eq. ( 4) yield the following single form after using Eq. ( 5) in Eq. ( 4). where For the ongoing problem, the direction of the time-dependent flow is taken towards the x-axis.So, we have the following form of velocity elements.After using Eq. ( 7) in Eqs.(1) and ( 6) and Boussinesq's approximation, the ongoing problem has the following form of momentum and energy equations 1, 25,26 With the Rosseland approximation of radiative heat flux q r = − 4σ * 3k 1 ∂T 4 ∂Y , Eq. ( 9) has the following expression 26 The physical boundary conditions of the concerned problem are described as 27 Now, familiarize the following non-dimensional variables 27 Equations ( 8) and (10) take the following form after the implementation of Eq. (12).
The dimensionless form of Eq. ( 11) has the following expressions

Fractional model
We implement the Caputo-Fabrizio time fractional derivative of order α ∈ [0, 1] to get the fractional model of the dimensionless system (13-14).In this fractional model, the time derivative is swapped with the Caputo-Fabrizio time fractional derivative which has the following form 28 In view of Eq. ( 17), Eqs. ( 13) and ( 14) take the following form

Exact solutions
To acquire the exact outcomes of the above Eqs.( 18) and ( 19) with fractional derivative, we utilized the technique of Laplace transform.The powerful methodology of Laplace transform is utilized to examine the transient linear systems.This methodology is basically consistent for linear systems.With the help of this technique, various linear models used in physical systems can be effectively handled.The initial conditions through this methodology can be accurately incorporated into the analysis.The stability analysis can be accomplished in the engineering control system via this method.Firstly, we execute the Laplace transform method on Eq. ( 19) and then obtain the solution of Eq. (18).The Laplace transforms of velocity, temperature, and Caputo-Fabrizio time fractional derivative are defined in the following equation 28 By using Eq. ( 20), Eq. ( 19) has the following form where γ = 1 1−α .The transformed boundary conditions of temperature are given as Equation ( 21) has the following solution corresponding to Eq. ( 22) where c 1 = 0 and c 2 = 1 s .Another form of Eq. ( 23) can be written as follow According to Eq. ( 25), the inverse Laplace transform of Eq. ( 24) has the following expression, To evaluate the above integral, we use the series representation of an exponential function.So, we have After evaluating the above integral, Eq. ( 27) can be written as We get the exact solution of the temperature distribution after using Eqs.( 27) and ( 28) as follows Now to acquire the solution of the velocity distribution given in Eq. ( 18), we apply the Laplace transform technique and get the following form The transformed conditions are Equation (30) has the following solution Using Eq. ( 31), we have the values of constants c 3 = 0, c 4 = − Gr(s+αγ ) s 2 γ −(1+N)Pr eff γ .Equation ( 32) can be written as ( 24) 4 . ( Vol.:(0123456789)The more suitable expression of Eq. ( 34) can be written as follows, The solution of velocity distribution is achieved after applying the inverse Laplace transformation on Eq. (36) as follows, where ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 can be evaluated in a same manner as ψ 2 (y, t, g, h) in Eq. (25).

Analysis of results
The fractional time derivative study of Nadeem trigonometric non-Newtonian fluid with an unsteady radiative flow over a vertical infinite plate has been examined in this article.The execution of adequate variables on the constitutive equations of the problem yields the dimensionless setup of equations.To get more significant results of the flow phenomenon, the fractional model of the Caputo-Fabrizio time derivative is implemented on the dimensionless system of coupled equations.The exact solutions of the temperature distribution and velocity distribution are obtained through an effective methodology of the Laplace transformation.For the correctness and confirmation of the ongoing problem, the temperature distribution values for specific parameters are compared with the previous study in Table 1.These values exhibit an excellent bonding with the previous findings, which indicate the validation of the present analysis.The physical interpretation of temperature and velocity fields relative to the pertinent parameters and time are graphically explored.The flow and thermal distributions of the considered fluid regarding different time periods and physical parameters are physically visualized.From the graphical representation of the temperature distribution and velocity distribution, it is perceived that the graphical results are more significant in the context of a larger period of time as compared to a smaller period of time.

Discussion of results
Figures 2a,b-5a,b are prepared to examine the influence of the pertinent parameters on the velocity field and temperature field regarding distinct periods of time.The purpose of organizing Fig. 2a,b is to disclose the nature of the velocity curve relative to the distinct time periods and escalation of the Grashof number.The increment of the Grashof number for both intervals of time develop an augmentation in the field of velocity.The reason for this graphical phenomenon is that the Grashof number characterized the connection of buoyancy and restraining forces.The reason for buoyancy forces is the variation of the fluid density across space.On the other hand, the viscosity of the fluid is the cause of the restraining forces.According to Shah et al. 29 , an increase in buoyancy (33) Table 1.Comparison of the values of the temperature distribution with Fetecau et al. 32 for t * = 0.5 results in a higher velocity of fluid flow, as described by Archimedes' principle.This principle exhibits that the buoyant force operating on an object submerged in a fluid is equivalent to the displaced fluid's weight.In a related study, Animasaun et al. 30 identified that when an object in a fluid experiences a greater buoyant force, it effectively reduces the net force acting against its motion, leading to an increase in velocity.Accordingly, the fluid accelerates objects with enhanced buoyancy more rapidly, resulting in a corresponding acceleration in the fluid velocity.This principle executes practical implementation in fields such as fluid dynamics, where the movement of objects in fluids is examined.Moreover, it plays a role in natural phenomena like the rising of hot air in a less dense and buoyant manner, contributing to atmospheric circulation patterns 31 .Figure 3a,b is sketched to observe the accelerating order of the fractional derivative at two different times on the curve of the fluid velocity.The larger magnitude of the fractional derivative yields an acceleration in the field of velocity.In Fig. 3b, at the large time period, the increasing effects of the fractional parameter are dominating near the plate as compared to the small period of time in Fig. 3a.According to the physical point of view, the thickness of the boundary layer becomes augmented due to the increment of fractional derivative order.Consequently, the flow field exhibits accelerating behavior.At both small and large time periods, Fig. 4a,b depicts the nature of temperature distribution corresponding to the escalating order of fractional derivative.The profile of the temperature is augmented with the increasing order of the fractional derivative for both time periods.The reason behind this phenomenon is that when the fractional operator's order is enhanced, the memory effect confined in the fractional operator becomes augmented which further affects the time period.Accordingly, with the increasing fractional parameter or fractional operator's order, the period of time is enhanced which escalates the temperature profile.The

Conclusion
The ongoing study focuses on the time-dependent flow phenomenon with the heat transfer mechanism in a trigonometric non-Newtonian fluid based on the fractional model.The fractional model of a trigonometric non-Newtonian fluid is obtained through the application of the Caputo-Fabrizio fractional derivative of order α .The thermal analysis of the problem is described with the significance of the thermal radiation.The suitable similarity variables yield the dimensionless setup of the equations.The implementation of the Laplace transform technique on the dimensionless setup of the equations provides the exact outcomes.This study has the following essential outcomes.
• Both the temperature and velocity distributions have dominating impacts related to the pertinent parameters at the large time period in comparison to the small time period.

Figure 1 .
Figure 1.Physical system of problem.

Figure 2 .
Figure 2. (a,b) Velocity distribution regarding Grashof number and time.

Figure 3 .
Figure 3. (a,b) Velocity distribution regarding fractional parameter and time.

Figure 4 .
Figure 4. (a,b) Temperature distribution regarding fractional parameter and time.

Figure 5 .
Figure 5. (a,b) Temperature distribution regarding effective Prandtl number and time.