Power Law Kernel Analysis of MHD Maxwell Fluid with Ramped Boundary Conditions: Transport Phenomena Solutions Based on Special Functions

In this paper, a new approach to find exact solutions is carried out for a generalized unsteady magnetohydrodynamic transport of a rate-type fluid near an unbounded upright plate, which is analyzed for ramped-wall temperature and velocity with constant concentration. The vertical plate is suspended in a porous medium and encounters the effects of radiation. An innovative definition of the time-fractional operator in power-law-kernel form is implemented to hypothesize the constitutive mass, energy, and momentum equations. The Laplace integral transformation technique is applied on a dimensionless form of governing partial differential equations by introducing some non-dimensional suitable parameters to establish the exact expressions in terms of special functions for ramped velocity, temperature, and constant-concentration fields. In order to validate the problem, the absence of the mass Grashof parameter led to the investigated solutions obtaining good agreement in existing literature. Additionally, several system parameters were used, such as as magnetic value M, Prandtl value Pr, Maxwell parameter λ, dimensionless time τ, Schmidt number “Sc”, fractional parameter α, andMass and Thermal Grashof numbers Gm and Gr, respectively, to examine their impacts on velocity, wall temperature, and constant concentration. Results are also discussed in detail and demonstrated graphically via Mathcad-15 software. A comprehensive comparative study between fractional and non-fractional models describes that the fractional model elucidate the memory effects more efficiently.


Introduction
The study of non-Newtonian fluids is a very interesting topic for many researchers due to its wide applications in modern technologies and many industrial sectors. Non-Newtonian fluids deals with problems related to chemical engineering, polymers, petroleum, paints, pharmaceutical, plastics, and food processing industries. Furthermore, these types of models also represent fluids such as biological materials, ketchup, polymeric liquids, motor oils, personal care products (shampoo, creams, toothpastes), complex mixtures, etc. However, in the literature, different types of fluid models exist, included rate-type, differential, and integral models. However, many researchers have addressed the rate-type fluid model, because from the research point of view, it is more applicable, due to its elastic effects and memory. Therefore, in this present work, rate-type models are recognized as Maxwell models. Further, Maxwell designed this model to estimate air's viscoelastic characteristics [1]. Furthermore, Jordan et al. [2] studied exact solutions for Stokes problems of the first type, including the Maxwell model. Moreover, Fetecau et al. [3] explored Stokes problems of second type by assuming the flow of fluid placed on a surface of infinite plate that is vacillating on the plane. Further, the model generated by Maxwell is known as an uncomplicated model designed for expecting the rheological performance of fluids, and it has both viscosity and elasticity as features. Fetecau [4] invetigated exact solution for an infinite vertical surface. The effects of thermophorsis and chemically reactive impacts on flow of Maxwell fluids were studied by Noor [5].
There is no denying the fact that when heat is extracted from a high-temperature wall through a fluid with a certain movement, this mechanism is referred to as heat convection (advection). This mechanism involves interactions between molecular diffusion and the motion of the fluid on the basis of natural as well as forced convection. This is because convective heat transfer transpires when the surface temperature varies from that of the encompassing fluid [6]. Solangi et al. [7] interrogated unique heat conduction properties for increasing concentrations. The focus of this work was to discuss the particle size control for heat exchange and mass concentration behavior on fluids. Soomro et al. [8] carried out the typical analysis of a stretching surface to develop the heat transmission for the flow of non-Newtonian nano-fluids. For physical aspects, numerical computation has been invoked to understand the governing equations based on the finite difference schemes that describe heat transfer phenomena. Shafiq et al. [9] explored an interesting study of the magnetohydrodynamic convective flow to explore the transfer rate of heat, motile microorganisms, and mass. They emphasized in their parametric study the problem for Brownian motion, buoyancy forces, thermophoretic, magnetic field, and Newtonian heating for temperature and concentration. Kashif et al. [10] investigated dual thermal analysis to understand the role of temperature dissimilarity versus temperature or time by using the suspension of nanoparticles. The core objective of this study was to highlight the problem through fractional operators and special functions. Heat transmission over the stretching sheet based on the magnetohydrodynamic two-dimensional Casson fluid flow was observed by Hamid et al. [11]. They examined linearly convected heat effects on twofold solutions in which numerical stability was discussed for the dual results from the governing equation of the problem. Abro et al. [12] suggested the fractional study for thermal radiation of Jeffery fluid and reduced the sub-solutions from fluid motion of the second grade with and without magnetic field. The thermal properties of the governing equations have been treated by integral transform approach. Sheikholeslami et al. [13] observed magnetizable hybrid fluid in the core of a circular cavity with two heaters in circular form subjected to carbon nanotubes. Their focus point in this work was to create magnetic strengths in the wires using electrical current. Abdelmalek et al. [14] applied a hybrid technique known as Control Volume Finite Element Scheme to a curvy circular heater with nano-sized particles using convective heat transmission. They suggested the finding and concluded that conformation of the curvy heater played a essential role in managing the heat emission rate and controlling the convectional flow within the enclosure. Kashif [15] applied surface modification technology to analyze the thermo-dissipation effects on time-dependent natural convectional flow of fluid. Finite Fourier sine transform, Laplace, and fractional techniques have been utilized in the governing equations for exhibiting typical and rheological properties of the problem [16][17][18]. Although studies on heat and transfer analysis can be continued, relevant studies can be observed in categorical formats such as heat transfer via analytical approaches [19][20][21][22][23][24][25][26], heat transfer via numerical approaches [27][28][29][30][31][32], heat transfer via fractional calculus approaches [33][34][35][36], and heat transfer via multi-dimensional approaches [37][38][39][40].
Fractional order calculus has been rapidly increasing in popularity these days because of its useful and exclusive features that are absent in non-fractional calculus. For example, classical models are unable to predict the earlier state of the process known as the memory effect, but fractional calculus is famous for having memory effects. It has been wellknown in recent decades for several real-world physical problems that fractional operators are more effective in differentiation compared to integral derivatives. Consequently, in modern times, a generalization of problems from the classical setting to the fractional setting is a subject of concern for many researchers. Generally, most circuit models are explained in terms of fractional differential equations. Therefore, fractional calculus is used in several fields, such as viscoelasticity, economics, relaxation processes, and fluid dynamics. It can be found in the literature that several fractional differential operators exist; for instance, Caputo with kernel (singular and local), Atangana-Baleanu with kernel (non-singular and non-local), Caputo-Fabrizio with kernel (non-singular and local), and a few others have been discussed [41][42][43][44]. For local and non-local kernels, convective flow with ramped conditions of temperature has been studied by Riaz et al. [45]. Furthermore, in a comparative study of MHD Maxwell fluid, the heat effect, with the application of local and non-local operators, is highlighted by Riaz et al. [46]. Some other fractional-associated references have been investigated [47][48][49][50][51], dealing with non-integer differential operators, MHD Jeffrey fluid movement, heat transport, and second-grade fluid. For more details, see [52][53][54][55].
Talha Anwar et al. [56] recently discussed the same problem with a classical approach, but without analyzing the diffusion equation, and executed the approximate result for the proposed problem by using the Laplace transformation technique and Durbin's numerical algorithm [57]. Based on the aforementioned literature, the object of this exploration is to develop the fractional model by using the modern interpretation of Caputo fractional operator and then derive the exact solution of the considered problem and make the comparison with results obtained by Talha Anwar et al. [56]. The consequences of different related physical parameters are explored, namely thermal Grashof number Gr, magnetic field M, dimensionless time τ, fractional parameter α, Schmidt number Sc, mass Grashof number Gm, and Prandtl number Pr, on non-dimensional velocity, concentration, and temperature. The results are discussed in detail and demonstrated graphically via Mathcad-15 software. Motivated by the above consideration, the main theme of this manuscript is to show the significance of fractional models corresponding to ordinary existing models for Maxwell fluid. In the end, solutions for some recently published work are compared with the current study, which endorses the authenticity of our solutions and proves that those investigations are limiting or special cases of the current problem.

Mathematical Model
Let us assume that the unsteady MHD, natural convection, in-compressible, time dependent, viscous motion of Maxwell fluid is very close to an infinite vertical plate embedded in a permeable medium with constant concentration and ramped conditions on temperature. In this case, consider the cartesian coordinates system (x,y); the plate is placed in a plane such as the vertically oriented x-axis and the y-axis in normal direction. Initially, for t = 0, the fluid and plate exhibit no motion and both possess uniform temperature T ∞ and concentration C ∞ . Moreover, the velocity can be stated piece-wise such that for 0 < t ≤ t 0 , it depends on time as u(0, t) = u 0 (t/t 0 ), and later for t > t 0 , the plate describes a constant motion with velocity u 0 . The physical quantities of the models' flow are described in Figure 1 The principle governing partial differential equations under usual Boussinesq's approximation with smaller Reynolds number and above mentioned assumptions are modeled as [43,44]: with initial and boundary conditions where and Introducing the set of dimensionless quantities are: After employing the dimensionless quantities, ignore the asterisk * notation; the following partial differential equations in dimensionless form are derived as: ∂θ(ψ, τ) ∂τ with conditions in dimensionless form being where

Preliminaries
Caputo fractional operator is described as and Laplace transformation of Equation (16) is obtained as where ℘ is named as a fractional parameter.

Solution of the Problem
To find the solution to the considered problem, the Laplace transformation technique is employed.

Exact Solution of Heat Profile with Caputo Time Fractional Derivative
The temperature Equation (10) in terms of Caputo fractional derivative is By applying the Laplace Transformation technique to write the solution to Equation (18) with conditions as Equations (13)-(15), we get Its solution is given byθ We applied boundary conditions for temperature given by Equation (20) to determined unknown constants, and we got which can be written asθ To get the required solution to Equation (23) we used Laplace inverse transformation, which is written as: where O(Φ) denotes a standard Heaviside function, Φ = τ − 1 and It is difficult to find θ c (ψ, τ) in exponential form, so we expressθ c (ψ, s) in its equivalent form asθ Employing inverse Laplace transformation, we get

Exact Solution of Mass Profile with Caputo Time Fractional Derivative
The fractional version of concentration Equation (12) in the form of Caputo fractional operator is derived as Solving Equation (28) using (13)- (15) and employing the Laplace Transformation technique, the resulting equations are written as: The solution in general form is The values of constants c 1 and c 2 , conditions are implemented given above Equation (30) for concentration, soC It is complicated to find C(ψ, τ) in exponential form, soC(ψ, s) in its equivalent form is written asC By employing inverse Laplace transformation, we obtained
Applying inverse Laplace transformation gives where O(τ − 1) is a standard Heaviside function. To calculate the value of u c (ψ, τ), we first expressū c (ψ, s) in summation form, which is equivalent to exponential form. Therefore, it can be written as Employing inverse Laplace transformation, we obtained .τ 1−αp−r .
Next, considerū 2 (ψ, s) to calculate u 2 (ψ, τ): This can be written in modified form as Taking inverse Laplace transformation yields Employing inverse Laplace transformation, we obtained similarly, the solution for u 3 (ψ, τ) by consideringū 3 (ψ, s) is derived, so we writē Finally, by applying Laplace inverse transformation, we get Employing inverse Laplace transformation, we obtained Furthermore, Employing inverse Laplace transformation, we obtained The function G a,b,c (., τ) used in the above expressions is known as the Generalized Lorenzo Hartly function, and its Laplace inverse is defined as: L −1

Solution of Shear Stress
The fractional version of Equation (11) in the form of a Caputo fractional operator is constructed as The solution to Equation (56), obtained by employing Laplace transform, is To find the value of dū(ψ,s) dψ , differentiating Equation (39) with respect to ψ yields where § = (1 + λs α )(M + s).
Substituting Equation (58) into Equation (57), which gives the solution of shear stress: Since the solutions given in Equation (59) contain complex terms of the Laplace parameter s, to derive the solution in a real time τ, we applied numerical inversion method known as Durbin Method [57].

Limiting Cases
As α → 1, we get the same expression for velocity profile as Talha Anwar et al. [56] obtained with ramped temperature, without considering the term Gm, i.e, Gm = 0. Furthermore, we derive the same result for velocity profile and temperature distribution of a viscous fluid when λ = 0 and Gm = 0 as Seth et al. [58]. This proves the authenticity of our newly established solution expressions.

Results and Discussion
In this paper, we investigated the effects of constant concentration, temperature, and velocity with ramped conditions on time-dependent MHD convective flow of Maxwell fluid. We derived exact solution expressions in terms of special functions for dimensionless velocity, concentration, and ramped wall temperature via a Caputo fractional time operator by using the method of Laplace transformation as in the existing literature. To analyze the physical importance of the proposed problem, estimated the results with the help of graphically illustration for several values of various parameters involved in the problem such as Schmidt number Sc, Prandtl number Pr, magnetic number M, mass Grashof number Gm , fractional parameter α, and thermal Grashof number Gr. In Figures 2-8; behavior of temperature, concentration, and velocity are discussed for two different values of α. In Figures 9-15, we compare the classical result with the Caputo fractional time derivative for temperature, concentration, and velocity for involved parameters in the considered problem, and Figures 16-18 present graphs for varying values of α. Figure 2 describes the temperature for distinct values of Pr with α = 0.1 and α = 0.9; when the value of Pr increases, the result is falling the temperature. Figure 3 illustrates the behavior of the concentration curve for various values of the Schmidt number Sc parameter with α = 0.1 and α = 0.9; it is deceptive that for maximum values of Sc, the concentration curve decays. The involvement of the concentration factor of fluid velocity on the fluid flow is very important and cannot be overlooked. Figure 4 describes the effects of Gm when α = 0.1 and α = 0.9. The mass Grashof number is generally defined as the ratio of mass buoyancy force to viscous force, which causes unrestricted convection. It is noticed that the velocity is enhanced in the case of high values of Gm. Figure 5 presents the impacts of Gr with α = 0.1 and α = 0.9. It is remarked that the velocity contour is elevated due to the large values of Gr. Thermal Grashof number is the proportion of thermal buoyancy force to viscous force, which causes free convection. Figure 6 illustrates the variation in velocity field for different values of M with α = 0.1 and α = 0.9; it can be observed that the fluid velocity declines when M is maximum, causing an increase in the Lorentz force, which act as a dragging force along with allied forces, thereby leading to a fall in the fluid velocity. Furthermore, this force becomes weaker when this is far away from the plate surface, and eventually fluid flow stops. In Figure 7, the impact of Pr on the velocity with α = 0.1 and α = 0.9 is analyzed. It causes a reduced fluid velocity for every increase in the values of Pr. In Figure 8, the effect of Sc with α = 0.1 and α = 0.9 on velocity profile is analyzed: the flow decreases as Sc increases. Schmidt number Sc signifies concentration, and velocity in free convection flow regimes of fluids are by the relative effectiveness momentum.
In Figures 9-15, we observe that temperature, concentration and velocity profile with the Caputo fractional time derivative displays more decay than the classical approach presented by Talha Anwar et al. [56]. Therefore, it is concluded that the Caputo fractional derivative used in the present paper efficiently shows the decay in temperature field, concentration, and velocity profile compared to the classical approach. Furthermore, in Figures 16-18 , it is noticed that we recovered the classical result as α → 1.

Conclusions
In this paper, we analyzed the effects of ramped temperature and velocity with constant concentration on MHD, time-dependent, and convective flows of Maxwell fluids. The governing partial differential equation is inscribed into a dimensionless form. To develop the fractional model by using the modern interpretation of the Caputo fractional time-derivative operator and accomplish the comparison with ordinary models. The Laplace transformation technique is applied to establish the analytical solution for velocity and the concentration and energy equations in terms of special functions for the proposed problem. Furthermore, conferred dissimilar parameters, i.e., Prandtl number Pr, Maxwell parameter λ, fractional parameter α, Schmidt number Sc, Magnetic field M, and Mass and Thermal Grashof numbers Gm and Gr, respectively, were used to analyze the impacts of all these parameters on fluid velocity field, constant concentration, and ramped wall temperature with the help of graphical illustrations. Some noteworthy remarks and concluding results from this work are as follows: • The temperature field declines with larger values of P r . Furthermore, we noted a reduced concentration for increasing values of Sc. Data Availability Statement: During the current study, no data sets were developed or investigated. Therefore, no data sharing is applicable to this article.

Conflicts of Interest:
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.