# Effect of Local and Non-Local Kernels on Heat Transfer of Mixed Convection Flow of the Maxwell Fluid

Authors: M. B. Riaz; A. Atangana; Maryam Asgir; Muhammad Altaf Khan; Hafte Amsalu Kahsay
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The heat transfer study of mixed convection flow of the Maxwell fluid is carried out here. The fluid flow is demonstrated by the system of coupled partial differential equations in the dimensionless form firstly. Then, its fractional form is developed by using the new definition of the noninteger-order derivative with the singular kernel (Caputo/C) and nonsingular kernels (Caputo–Fabrizio/CF and Atangana–Baleanu (nonlocal)/ABC). The hybrid-form solutions are obtained by applying the Laplace transform, and for the inverse Laplace transform, the problem is tackled by the numerical algorithms of Stehfest and Tzou. The C, CF, and ABC solution comparison under the effects of considered different parameters is depicted. The physical aspects of the considered problem are well explained by C, CF, and ABC in comparison to the integer-order derivative due to its memory effects. Furthermore, the best fit model to explain the memory effects of velocity is CF. The solutions for the Newtonian fluid and ordinary Maxwell fluid are considered as a special case and found in the literature.


---

## Body

## 1. Introduction

The combination of natural and forced convection is called mixed convection. This phenomenon gained a lot of popularity in recent years due to its vast applications in many fields of engineering, such as nuclear, chemical, food, aerospace, electrical, fluid dynamics, and astrophysics. The combined convection process comes into reality when free or forced convection alone is not enough to describe the heat transfer process properly. Many diverse significant effects of mixed convection in the heat transfer phenomenon were seen in [1–4]. Integer-order derivatives are used widely to model real-world problems mathematically. The nonlocal nature of the noninteger-order derivative makes it a more efficient tool as compared to the integer-order derivative. The property of the noninteger derivative to get a better insight of the rheological and hereditary properties of the material under consideration makes it more valuable. It has become essential in modeling the description of many complex dynamics in the field of industries, engineering, technology, etc.Non-Newtonian fluids are more important as compared to the viscous and Newtonian fluids because they depict the complex phenomenon of real-life problems. Mostly, the problems that exist in engineering and industries exhibit a nonlinear relationship. Non-Newtonian fluids connect the shear stress and shear strain in a nonlinear relationship. This nonlinear relationship of the shear stress and shear strain makes the non-Newtonian fluid flow more subtle. Maxwell fluid is one of the widely studied non-Newtonian fluids due to its diversity and rheological properties [5]. The nature of the viscoelastic fluids is best described by the noninteger/fractional derivative due to the history of fluid flow. Bagley and Torvik [6] and Germnat [7] were pioneers who introduced the fractional approach in the viscoelastic fluids. The characteristics of the viscoelastic fluids were investigated by Mainardi and Spada [8] via a fractional approach, and they found that the results obtained were in a good settlement with the data given for the experiment. Caputo, Caputo–Fabrizio, and Atangana–Baleanu have introduced the definition of the fractional derivative which has been widely used in solving many problems of fluid dynamics and other engineering and science disciplines.The flow properties of the fractional Maxwell fluid between two parallel plates were investigated by Haito and Mingyu [9]. They obtained the solution of the pressure gradient by finite Fourier cosine and Laplace transform. Jamil [10] explored the analytical solutions of the fractionalized Maxwell fluid flowing past on a vertical plate with slip effects. He developed the results for the ordinary Maxwell and viscous fluid as a special case. The problem of unsteady MHD flow of the Maxwell fluid in a porous medium was solved by Khan et al. [11]. They studied the motion of the fluid induced due to the sine or cosine oscillation of the plate and revealed the interesting behaviors of velocity and shear stress. Vieru and Zafar [12] examined the Couette flow of the Maxwell fluid with the slip condition. The results for velocity and shear stress due to the constant and sinusoidal motion of the plate are obtained by the Laplace integral transform.The heat transfer phenomenon in the Maxwell fluid was studied by Imran et al. [13] and many others via a fractional approach. The analysis of heat transfer in the fractionalized Maxwell fluid under the effects of MHD and second-order slip effects was carried out by Aman et al. [14]. Semianalytical results for velocity and shear stress were obtained and compared graphically. Raza and Ullah [15] used the C and CF derivatives to study the heat transfer in the Maxwell fluid under the consideration of Newtonian heating. The comparison results revealed that the value of fluid temperature and velocity is lower in the C case. Merkin and Mahmood [16] studied the aspects of mixed convection of external flow of the surface which was chemical reactive and immersed in a porous medium. For vertical surfaces, such studies of mixed convective flow were conducted by Minto et al. [17]. The study of the mixed convection flow in a horizontal channel and under constant heat flux revealed interesting results which were demonstrated by Chou and Tsern [1]. Khan et al. [2] considered the exothermal and isothermal stretching disk with the phenomenon of mixed heat convection and presented the analytical results.The purpose of this communication is to consider the Maxwell fluid over an oscillating vertical plate with constant wall temperature and analyze the phenomenon of mixed convective heat transfer via the fractional derivative definition of C, CF, and ABC. Laplace transform is used to attain the solutions of the problem under deliberation, and these results of heat transfer analysis have much importance in the industry such as in cooling the circuits, radiator, heat exchanger, deicing system, solar panels, boilers, and many other home electric and industrial appliances.


## 2. Mathematical Description of the Problem

The incompressible Maxwell fluid unsteady mixed convective flow phenomenon on the oscillating vertical plate is considered here. Both fluid and plate are at rest at the start, and temperature isT∞.As the time starts to increase, the plate starts to move in sinusoidal motion with the velocityV=W0Htcosat (W0 is the amplitude of the motion, Ht is the unit step function, and a represents the frequency of the plate oscillation) which results in the start of the flow of the Maxwell fluid. The temperature of the plate increases/decreases to a constant value Tw. With the increase in time, the velocity reaches zero and temperature to free stream temperature T∞. The heat transfer phenomenon of Maxwell fluid flow due to mixed convection along shear stress is described by the following PDEs [3]. The list of parameters is given in Table 1.(1)ρ1+λ∂∂t∂wη,t∂t=μ∂2wη,t∂η2+1+λ∂∂tρgβT−T∞,1+λ∂∂tτη,t=μ∂wη,t∂η,ρcp∂Tη,t∂t=k∂2Tη,t∂η2,and the imposed initial and boundary conditions are(2)wη,t=0,Tη,t=T∞,att=0,wη,t=W0Htcosat,Tη,t=Tw,atη=0,wη,t=0,Tη,t=T∞,asη⟶∞.Table 1
Nomenclature.
SymbolQuantityρFluid densityλRelaxation timeμViscosity (dynamic)βVolumetric coefficientcpSpecific heatgAcceleration due to gravityThe dimensionless relations are defined as(3)w∗=wW0,η∗=ηW0υ,t∗=tW02υ,a∗=aυW02,τ∗=υτμW02,ϕ=T−T∞Tw−T∞,λ∗=λW02υ,Gr=νgβTw−T∞W03,Pr=μcpk.The dimensionless system of partial differential equations after dropping the “∗” notation is(4)1+λ∂∂t∂wη,t∂t=∂2wη,t∂η2+1+λ∂∂tGrϕη,t,1+λ∂∂tτη,t=∂wη,t∂η,Pr∂ϕη,t∂t=∂2ϕη,t∂η2,and initial and boundary conditions in the dimensionless form are(5)wη,t=0,ϕη,t=0,t=0,η>0,(6)wη,t=Htcosat,ϕη,t=1,t>0,η=0,(7)wη,t=0,ϕη,t=0,η⟶∞.


## 3. Preliminaries

Caputo (C) fractional time derivative is defined as(8)DCtαgη,t=1Γ1−α∫0tt−ξ−αg′η,ξdξ,α∈0,1,and its Laplace transform(9)LDCtαgη,t=sαLgη,t−sα−1gη,0.Caputo–Fabrizio (CF) fractional time derivative is(10)DCFtαgη,t=11−α∫0te−αt−ξ/1−αg′η,ξdξ,α∈0,1,with the Laplace transform defined as(11)LDCFtαgη,t=sLgη,t−gη,0s1−α+α.Atangana–Baleanu (ABC) time derivative is defined as(12)DCtαgη,t=11−α∫0tEα−αt−ξα1−αg′η,ξdξ,α∈0,1,and its Laplace transform is defined as(13)LDABCtαgη,t=sαLgη,t−sα−1gη,0sα1−α+α.


## 4. Fractional Formulation of the Problem

### 4.1. Caputo Time Fractional Model

(14)1+λDCtα∂wη,t∂t=∂2wη,t∂η2+1+λDCtαGrϕη,t,(15)1+λDCtατη,t=∂wη,t∂η,(16)PrDCtαϕη,t=∂2ϕη,t∂y2.


### 4.2. Caputo–Fabrizio Time Fractional Model

(17)1+λDCFtα∂wη,t∂t=∂2wη,t∂η2+1+λDCFtαGrϕη,t,(18)1+λDCFtατη,t=∂wη,t∂η,(19)PrDCFtαϕη,t=∂2ϕη,t∂η2.


### 4.3. Atangana–Baleanu Time Fractional Model

(20)1+λDABCtα∂wη,t∂t=∂2wη,t∂η2+1+λDABCtαGrϕη,t,(21)1+λDABCtατη,t=∂wη,t∂η,(22)PrDABCtαϕη,t=∂2ϕη,t∂η2.


## 4.1. Caputo Time Fractional Model

(14)1+λDCtα∂wη,t∂t=∂2wη,t∂η2+1+λDCtαGrϕη,t,(15)1+λDCtατη,t=∂wη,t∂η,(16)PrDCtαϕη,t=∂2ϕη,t∂y2.


## 4.2. Caputo–Fabrizio Time Fractional Model

(17)1+λDCFtα∂wη,t∂t=∂2wη,t∂η2+1+λDCFtαGrϕη,t,(18)1+λDCFtατη,t=∂wη,t∂η,(19)PrDCFtαϕη,t=∂2ϕη,t∂η2.


## 4.3. Atangana–Baleanu Time Fractional Model

(20)1+λDABCtα∂wη,t∂t=∂2wη,t∂η2+1+λDABCtαGrϕη,t,(21)1+λDABCtατη,t=∂wη,t∂η,(22)PrDABCtαϕη,t=∂2ϕη,t∂η2.


## 5. Solutions for Temperature

### 5.1. Caputo Sense

By implementing the Laplace transform on equation (16),(23)∂2ϕ¯η,s∂η2−Prsαϕ¯η,s=0,and the solution of the above equation is given as(24)ϕ¯η,s=1se−ηPrsα.


### 5.2. Caputo–Fabrizio Sense

Laplace transform of equation (19) is(25)∂2ϕ¯η,s∂η2−b0ss+γϕ¯η,s=0.The above differential equation solution after applying the transformed boundary conditions ofϕη,s from (6) and (7) is(26)ϕ¯η,s=1se−ηb0s/s+γ,where a0=1/1−α,b0=Pra0,andγ=αa0.


### 5.3. Atangana–Baleanu Sense

Applying the Laplace transform on equation (22),(27)∂2ϕ¯η,s∂η2−b0sαsα+γϕ¯η,s=0,and the Laplace solution of the above equation is(28)ϕ¯η,s=1se−ηb0sα/sα+γ,where a0=1/1−α,b0=Pra0,andγ=αa0.


## 5.1. Caputo Sense

By implementing the Laplace transform on equation (16),(23)∂2ϕ¯η,s∂η2−Prsαϕ¯η,s=0,and the solution of the above equation is given as(24)ϕ¯η,s=1se−ηPrsα.


## 5.2. Caputo–Fabrizio Sense

Laplace transform of equation (19) is(25)∂2ϕ¯η,s∂η2−b0ss+γϕ¯η,s=0.The above differential equation solution after applying the transformed boundary conditions ofϕη,s from (6) and (7) is(26)ϕ¯η,s=1se−ηb0s/s+γ,where a0=1/1−α,b0=Pra0,andγ=αa0.


## 5.3. Atangana–Baleanu Sense

Applying the Laplace transform on equation (22),(27)∂2ϕ¯η,s∂η2−b0sαsα+γϕ¯η,s=0,and the Laplace solution of the above equation is(28)ϕ¯η,s=1se−ηb0sα/sα+γ,where a0=1/1−α,b0=Pra0,andγ=αa0.


## 6. Solutions for Velocity

### 6.1. Caputo Sense

By taking the Laplace transform of equation (14) and substituting the value of ϕ¯η,s from equation (24),(29)∂2w¯η,s∂η2−1+λsαsw¯η,s=−Gr1+λsαse−ηPrsα.The solution of equation (29) after the use of the corresponding transformed boundary conditions from (6) and (7) is(30)w¯η,s=ss2+a2e−ηs1+λsα+Gr1+λsαPrsα+1−s21+λsαe−ηs1+λsα−e−ηPrsα.


### 6.2. Caputo–Fabrizio Sense

By taking the value ofϕ¯η,s from equation (26), the transformed velocity equation becomes(31)∂2w¯η,s∂η2−c0s+γs+γsw¯η,s=−Grc0s+γs+γe−ηb0s/s+γs.The solution of the above equation is given as(32)w¯η,s=ss2+a2e−ηsc0s+γ/s+γ+Grc0s+γb0s2−s2c0s+γe−ηsc0s+γ/s+γ−e−ηb0s/s+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.


### 6.3. Atangana–Baleanu Sense

Imposing the Laplace transform on equation (20), with the substitution of the value of ϕ¯η,s from equation (28), we get(33)∂2w¯η,s∂η2−c0sα+γsα+γsw¯η,s=−Grc0sα+γsα+γe−ηb0sα/sα+γs.Solving the above equation, we get the solution(34)w¯η,s=ss2+a2e−ηsc0sα+γ/sα+γ+Grc0sα+γb0sα+1−s2c0sα+γe−ηsc0sα+γ/sα+γ−e−ηb0sα/sα+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.


## 6.1. Caputo Sense

By taking the Laplace transform of equation (14) and substituting the value of ϕ¯η,s from equation (24),(29)∂2w¯η,s∂η2−1+λsαsw¯η,s=−Gr1+λsαse−ηPrsα.The solution of equation (29) after the use of the corresponding transformed boundary conditions from (6) and (7) is(30)w¯η,s=ss2+a2e−ηs1+λsα+Gr1+λsαPrsα+1−s21+λsαe−ηs1+λsα−e−ηPrsα.


## 6.2. Caputo–Fabrizio Sense

By taking the value ofϕ¯η,s from equation (26), the transformed velocity equation becomes(31)∂2w¯η,s∂η2−c0s+γs+γsw¯η,s=−Grc0s+γs+γe−ηb0s/s+γs.The solution of the above equation is given as(32)w¯η,s=ss2+a2e−ηsc0s+γ/s+γ+Grc0s+γb0s2−s2c0s+γe−ηsc0s+γ/s+γ−e−ηb0s/s+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.


## 6.3. Atangana–Baleanu Sense

Imposing the Laplace transform on equation (20), with the substitution of the value of ϕ¯η,s from equation (28), we get(33)∂2w¯η,s∂η2−c0sα+γsα+γsw¯η,s=−Grc0sα+γsα+γe−ηb0sα/sα+γs.Solving the above equation, we get the solution(34)w¯η,s=ss2+a2e−ηsc0sα+γ/sα+γ+Grc0sα+γb0sα+1−s2c0sα+γe−ηsc0sα+γ/sα+γ−e−ηb0sα/sα+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.


## 7. Solutions for Shear Stress

### 7.1. Caputo Sense

The shear stress equation after applying the Laplace transform and simplification is(35)τ¯η,s=11+λsα∂w¯η,s∂η.Taking the derivative of velocity equation (30) with respect to η and putting the value in equation (35), we get the expression for shear stress:(36)w¯η,s=ss2+a2−s1+λsαe−ηs1+λsα+GrPrsα+1−s21+λsαPrsαe−ηPrsα−s1+λsαe−ηs1+λsα.


### 7.2. Caputo–Fabrizio Sense

The equation of the shear stress after employing the Laplace transform on equation (18) is(37)τ¯η,s=s+γc0s+γ∂w¯η,s∂η.Taking the derivative of velocity equation (32) with respect to η and putting the value in equation (37), we get the expression for shear stress:(38)w¯η,s=ss2+a2−ss+γc0s+γe−ηsc0s+γ/s+γ+Grs+γb0s2−s2c0s+γb0ss+γe−ηb0s/s+γ−sc0s+γs+γe−ηsc0s+γ/s+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.


### 7.3. Atangana–Baleanu Sense

Applying the Laplace transform on equation (21) and simplifying give(39)τ¯η,s=sα+γc0sα+γ∂w¯η,s∂η.The solution of the shear stress after substituting the derivative of velocity equation (34) in the above equation is(40)w¯η,s=ss2+a2−ssα+γc0sα+γe−ηsc0sα+γ/sα+γ+Grc0sα+γb0sα+1−s2c0sα+γb0sαsα+γe−ηb0sα/sα+γ−sc0sα+γsα+γe−ηsc0sα+γ/sα+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.Similarly, solutions for velocity, shear stress, and temperature for the ABC model were given by Abro et al. [18] (M = 0 and Φ=0 in Equation (2.8)). For the CF solution of velocity, temperature, and shear stress, we found the results in [19, 20] (α2=0 in Equation (7)).The Laplace inverse transform for the solution of velocity and temperature is attained by applying Stehfest and Tzou numerical algorithms [21, 22].Stehfest’s algorithm is defined as(41)ωΨ,τ=ln2τ∑z=12ndzω¯Ψ,zln2τ,where n is a positive integer and(42)dz=−1z+n∑i=z+1/2minz,nin2i!n−i!i!i−1!z−i!2i−z!,where the real number integer part is described by [.]. Tzou’s numerical algorithm is defined as(43)w¯γ,t=e4.7t12w¯γ,4.7t+Re∑ι=1N1−1ιw¯γ,4.7+zπit,where Re(.) is the real part, the imaginary part is i, and N1 >> 1 represents the natural number.


## 7.1. Caputo Sense

The shear stress equation after applying the Laplace transform and simplification is(35)τ¯η,s=11+λsα∂w¯η,s∂η.Taking the derivative of velocity equation (30) with respect to η and putting the value in equation (35), we get the expression for shear stress:(36)w¯η,s=ss2+a2−s1+λsαe−ηs1+λsα+GrPrsα+1−s21+λsαPrsαe−ηPrsα−s1+λsαe−ηs1+λsα.


## 7.2. Caputo–Fabrizio Sense

The equation of the shear stress after employing the Laplace transform on equation (18) is(37)τ¯η,s=s+γc0s+γ∂w¯η,s∂η.Taking the derivative of velocity equation (32) with respect to η and putting the value in equation (37), we get the expression for shear stress:(38)w¯η,s=ss2+a2−ss+γc0s+γe−ηsc0s+γ/s+γ+Grs+γb0s2−s2c0s+γb0ss+γe−ηb0s/s+γ−sc0s+γs+γe−ηsc0s+γ/s+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.


## 7.3. Atangana–Baleanu Sense

Applying the Laplace transform on equation (21) and simplifying give(39)τ¯η,s=sα+γc0sα+γ∂w¯η,s∂η.The solution of the shear stress after substituting the derivative of velocity equation (34) in the above equation is(40)w¯η,s=ss2+a2−ssα+γc0sα+γe−ηsc0sα+γ/sα+γ+Grc0sα+γb0sα+1−s2c0sα+γb0sαsα+γe−ηb0sα/sα+γ−sc0sα+γsα+γe−ηsc0sα+γ/sα+γ,where a0=1/1−α,b0=Pra0,c0=1+λa0,andγ=αa0.Similarly, solutions for velocity, shear stress, and temperature for the ABC model were given by Abro et al. [18] (M = 0 and Φ=0 in Equation (2.8)). For the CF solution of velocity, temperature, and shear stress, we found the results in [19, 20] (α2=0 in Equation (7)).The Laplace inverse transform for the solution of velocity and temperature is attained by applying Stehfest and Tzou numerical algorithms [21, 22].Stehfest’s algorithm is defined as(41)ωΨ,τ=ln2τ∑z=12ndzω¯Ψ,zln2τ,where n is a positive integer and(42)dz=−1z+n∑i=z+1/2minz,nin2i!n−i!i!i−1!z−i!2i−z!,where the real number integer part is described by [.]. Tzou’s numerical algorithm is defined as(43)w¯γ,t=e4.7t12w¯γ,4.7t+Re∑ι=1N1−1ιw¯γ,4.7+zπit,where Re(.) is the real part, the imaginary part is i, and N1 >> 1 represents the natural number.


## 8. Results and Discussion

A comparative study of heat transfer in mixed convective flow of the Maxwell fluid with the noninteger-order approaches C, CF, and ABC is elaborated. To obtain the solutions of velocity, temperature, and shear stress via C, CF, and ABC fractional approaches, the Laplace transform is employed, and numerical algorithms are used to gain the solutions of the inverse Laplace transform. Effects of a fractional parameter, Prandtl number, and Grashof number on velocity, temperature, and shear stress are depicted graphically, and a comparison is made between three approaches to reveal the interesting results.Figure1 shows the physical insight of the fractional parameter on the temperature profile. Temperature increases with an increment in fractional parameter values in the C, CF, and ABC approach along with the variation of time. We see that, with the increase in time, the temperature increases and reaches a steady state rapidly in CF and ABC. A slight difference is noticed in the result of CF and ABC.Figure 1
Temperature profile comparison via C, CF, and ABC with the variation of fractional parameterα, with Pr = 5.
(a)(b)(c)(d)(e)(f)The Prandtl numberPr controls the thickness of the thermal boundary layer, resulting in a decrease in temperature with time variation as seen in Figure 2. We see a sharp decline in the temperature profile in the CF model, and it tends to the steady state slowly in comparison with C and ABC.Figure 2
Temperature profile comparison via C, CF, and ABC with the Prandtl number variation, withα = 0.7.
(a)(b)(c)(d)(e)(f)The fractional parameterα controls the velocity of the fluid which is presented in Figure 3. As time varies from small to large, we see that the velocity becomes the increasing function of α. The reason, α increases, boundary layer thickness increases, resulting in acceleration of the velocity. We see the highest velocity for CF. For α⟶1, we will recover the solution for integer-order derivatives which are already present in [3].Figure 3
Velocity profile comparison via C, CF, and ABC with the variation of fractional parameterα, and Pr = 5 and Gr = 5.
(a)(b)(c)(d)(e)(f)Figure4 depicts the consequences of the Grashof number Gr on the velocity profile at varying time scale.Figure 4
Velocity profile comparison via C, CF, and ABC with the variation of the Grashof number Gr, and Pr = 5,a = 2, and λ=0.7.
(a)(b)(c)(d)(e)(f)The graphical illustration depicts the velocity of fluid flow rises with the increase in Gr. The buoyancy forces take control over the viscous forces and result in accelerating the velocity. CF shows the highest velocity compared to C and ABC.The significant impact ofPr on the velocity profile is depicted in Figure 5. The fluid flow decreases with the Pr rise. Pr increases, and the thermal conductivity decreases, resulting in an increase in viscosity which slows down the fluid flow. Clearly, the CF model shows the highest velocity.Figure 5
Velocity profile comparison via C, CF, and ABC with the variation of the Prandtl number Pr, and Gr = 5,a = 2, and λ=0.7.
(a)(b)(c)(d)(e)(f)Figure6 analyzes the influence of the Grashof number Gr on shear stress at different time scales. Shear stress increases with the increase in the Grashof number Gr near the boundary, and then there comes a critical point in y where it starts to decline. The critical point moves further for CF as compared to C, but it moves closer to the boundary in case of ABC.Figure 6
Shear stress comparison via C, CF, and ABC with the variation of the Grashof number Gr, and Pr = 0.3,a = 2, and λ=0.5.
(a)(b)(c)(d)(e)(f)To highlight thePr behavior on shear stress, we observed Figure 7. The quite opposite behavior of Pr is observed compared to Gr. The velocity decreases with an increase in time and Pr, and at the critical point, it changes its behavior in the opposite direction.Figure 7
Shear stress comparison via C, CF, and ABC with the variation of the Prandtl number Pr, and Gr = 10,a = 2, and λ=0.3.
(a)(b)(c)(d)(e)(f)Forλ⟶0, the results for Newtonian fluids are recovered and found in [23]. The results of the comparison of Maxwell and Newtonian fluids for C, CF, and ABC are illustrated in Figure 8.Figure 8
Velocity comparison via C, CF, and ABC for Maxwell and Newtonian fluids.
(a)(b)


## 9. Conclusion

The problem of the heat transfer due to mixed convective fluid flow in the Maxwell fluid is brought into consideration here. The noninteger-order derivative definitions of C, CF, and ABC were employed to formulate the problem. The semianalytical solutions were obtained by the Laplace transform and numerical inverse Laplace transform. The comparative analysis of the behavior of the solutions of the velocity, temperature, and shear stress under the influence of different pertinent parameters is brought to light. The following results are concluded from this analysis:(1)
Temperature decreases with the increase inPr for the variation in time(2)
Velocity increases with the rise of Gr values and declines with the increment in Pr(3)
Shear stress shows the increasing behavior up to a critical point and then decreases under the influence of Gr(4)
Shear stress shows the decreasing behavior up to a critical point and then increases under the impact of Pr(5)
While comparing between the noninteger-order derivative approaches (C, CF, and ABC), we noticed that the significance of CF is more due to the nonsingular kernel(6)
For the validation of our results, the limiting cases ofα⟶1andλ⟶0 were considered, and obtained results were already present in the literature



---

*Source: 9914408-2021-06-12.xml*