# Study of Heat Transfer under the Impact of Thermal Radiation, Ramped Velocity, and Ramped Temperature on the MHD Oldroyd-B Fluid Subject to Noninteger Differentiable Operators

Authors: Syed Tauseef Saeed; Ilyas Khan; Muhammad Bilal Riaz; Syed Muhammad Husnine
Journal: Journal of Mathematics
 (2020)

---

## Abstract

This theoretical study explores the impact of heat generation/absorption with ramp wall velocity and ramp wall temperature on the magnetohydrodynamic (MHD) time-dependent Oldroyd-B fluid over an unbounded plate embedded in a porous surface. The mathematical analysis of fractional governing partial differential equations has been established using systematic and powerful techniques of Laplace transform with its numerical inversion algorithms. The fractionalized solutions have been traced out separately through all fractional differential operators. Nondimensional parameters along with Laplace transformation are used to find the solution of temperature and velocity profiles. Fractional time derivatives are used to analyze the impact of fractional parameters (memory effect) on the dynamics of the fluid. While making a comparison, it is observed that the fractional-order model is the best to explain the memory effect as compared to classical models. The obtained solutions are plotted graphically for different values of physical parameters. Our results suggest that the velocity profile decreases by increasing the effective Prandtl number. Furthermore, the existence of an effective Prandtl number may reflect the control of the thickness of momentum and enlargement of thermal conductivity.


---

## Body

## 1. Introduction

The interest in studying problems involving non-Newtonian fluid flows has considerably grown for their wide range of applications: from drilling oil and gas wells and well completion operations to industrial processes involving waste fluids, synthetic fibre foodstuffs, and the extrusion of molten plastics. The attributes of fluid flow trace the diversity of the physical structure for non-Newtonian fluid flow. In such a fluid, stress and rate of strain have a nonlinear relationship. Oldroyd-B fluids have become a significant model of rate-type fluid. The procedure for the flow of rate-type fluids was discussed by Oldroyd [1]. Viscoelasticity has important implications due to the characterization of viscoelastic parameters (relaxation and retardation phenomenon), elastic shearing strain, thermal relaxation, and other rheological properties [2, 3]. In this regard, the thermodynamical analysis for the constitutive model of thermoplastic, viscoelastic, and viscoplastic was observed by Krairi and Doghri [4] through the Cattaneo–Christov heat model. The temperature distribution and relaxation time of the heat flux were emphasized by the temperature equation.The technique of fractional calculus has been used to formulate mathematical modeling in various technological developments, engineering applications, and industrial sciences. Different valuable work has been discussed for modeling fluid dynamics, signal processing, viscoelasticity, electrochemistry, and biological structure through fractional time derivatives. This fractional differential operator found useful conclusions for experts to treat cancer cells with a suitable amount of heat source and has compared the results to see the memory effect of the temperature function. As compared to classical models, the memory effect is much stronger in fractional derivatives. From the past to the present, modeling of different processes is handled through various types of fractional derivatives and fractal-fractional differential operators, such as Caputo (power law), Atangana–Baleanu (Mittag–Leffler law), Caputo–Fabrizio (exponential law), Riemann–Liouville, and modified Riemann–Liouville (power law with boundaries) [5–16]. Ramped wall velocity and temperature with MHD fluid flow are gaining attention of many researchers. Physically, the implementation of ramped wall velocity with temperature in real-life problems has a significant role, but mathematically, it is difficult to handle such conditions. Ramped heating is used to control and increase the temperature with adiabatic conditions in an effective way. Firstly, Ahmed and Dutta [17] discussed the simultaneous use of ramped velocity and ramped temperature. Seth et al. [18–20] investigated heat and mass transfer phenomena with ramp temperature conditions. Recently, Tiwana et al. [21] and Anwar et al. [22] analyzed the MHD Oldroyd-B fluid under the effect of ramped temperature and velocity.In context with fractional differential operators, convective flow with ramped wall temperature for nonsingular kernel was analyzed by Riaz et al. [23]. Moreover, Riaz et al. [24] investigated the study of heat and mass transfer in the MHD Oldroyd-B fluid with ramped wall temperature using local and nonlocal differential operators. Additionally, the recent studies on modern fractional differential operators and viscoelastic fluids can be traced out in [25–38]. For this problem, the noninteger differentiable operator is chosen for the fractional MHD Oldroyd-B model which is developed under thermal radiation, ramp velocity, and ramp temperature associated with physical initial and boundary conditions. The model is solved via the Laplace transform technique and inversion algorithm. The required results are displayed in graphs with physical arguments.


## 2. Problem Statement

We discuss unsteady magnetohydrodynamic (MHD) fractional convective Oldroyd-B fluid flow under Boussinesq approximations over an infinite plate. Figure1 represents the flow geometry of the magnetized Oldroyd-B fluid. Under these presumptions, the governing equation for the Oldroyd-B fluid with appropriate conditions is defined as follows [22]:(1)1+λ1∂∂τ∂wη,τ∂τ=υ1+λ2∂∂τ∂2wη,τ∂η2−1+λ1∂∂τσB02ρwη,τ+gβT1+λ1∂∂τT−T∞−1+λ2∂∂τμϕρk∗wη,τ,(2)ρCp∂Tη,τ∂τ=κ∂2Tη,τ∂η2−∂qr∂η+QoT−T∞.Figure 1
Geometrical presentation of the Oldroyd-B model.The appropriate conditions are given as follows:(3)τ>0,wη,0=0,Tη,0=T∞,∂wη,0∂η=∂wη,0∂τ=0,η≥0,τ>0,w0,t=Uoττo,for0<τ<τo;Uo,forτ>τo,T0,τ=T∞+Tw−T∞ττ0,for0<τ<τo;Tw,forτ>τo,τ>0,wη,τ⟶0,Tη,τ⟶T∞,asη⟶∞.The dimensionless parameters in equations (1) to (3) are mentioned in the following:(4)w=VU0,η=ζU0υ,τ=tU02υ,τ0=υU02,θ=T−T∞Tw−T∞,λ1=λU02υ,λ2=λrU02υ,Pr=μCpk,Q=υ2QoρCpU02,M=σB02υρU02,1K=ϕυ2k∗U02.Applying (4) into (1)–(3), we required a set of dimensionless governing equations in the form of PDE’s system presented as follows:(5)a1+λ∂∂τ∂Vζ,t∂t=1+λr∂∂τ∂2Vζ,t∂ζ2+1+λ∂∂τGrθζ,t−a2Vζ,t,(6)∂θζ,t∂t=1Preff∂2θζ,t∂ζ2+Qθζ,t,where a1=1+λM+λr/Kanda2=M+1/K.The dimensionless corresponding conditions can be given as follows:(7)Vζ,0=θζ,0=0,Vtζ,0=Vζζ,0=0,forζ≥0,θ0,t=V0,t=t,for0<t<to,1,fort>to,Vζ,τ⟶0,θζ,τ⟶0,forζ⟶∞.


### 2.1. Governing Equations in terms of Singular and Nonsingular Kernels

We define Caputo time derivative with its Laplace transform defined in the following [39]:(8)DtκCgζ,t=1Γn−κ∫0tgnτt−τκ+1−ndτ,LDtκCgζ,t=sκLgζ,t−sκ−1gζ,0.The Caputo–Fabrizio fractional derivative and its Laplace transformation are defined as follows [40]:(9)DtκCFgζ,t=1n−κ∫0texpκ1−κ1−κ∂gζ,τ∂τdτ,LDtκCFgξ,t=sLgζ,t−gζ,01−κs+κ.The Atangana–Baleanu fractional derivative and its Laplace transformation are defined as follows [40]:(10)DtκABCgζ,t=1n−κ∫0tEκκt−κκ1−κ∂gζ,τ∂τdτ,LDtκABCgζ,t=sκLgζ,t−sκ−1gζ,01−κsκ+κ.


## 2.1. Governing Equations in terms of Singular and Nonsingular Kernels

We define Caputo time derivative with its Laplace transform defined in the following [39]:(8)DtκCgζ,t=1Γn−κ∫0tgnτt−τκ+1−ndτ,LDtκCgζ,t=sκLgζ,t−sκ−1gζ,0.The Caputo–Fabrizio fractional derivative and its Laplace transformation are defined as follows [40]:(9)DtκCFgζ,t=1n−κ∫0texpκ1−κ1−κ∂gζ,τ∂τdτ,LDtκCFgξ,t=sLgζ,t−gζ,01−κs+κ.The Atangana–Baleanu fractional derivative and its Laplace transformation are defined as follows [40]:(10)DtκABCgζ,t=1n−κ∫0tEκκt−κκ1−κ∂gζ,τ∂τdτ,LDtκABCgζ,t=sκLgζ,t−sκ−1gζ,01−κsκ+κ.


## 3. Temperature Profile via Fractional Operators

### 3.1. Temperature Profile via the Caputo Approach

Generating equation (6) for the fractional form, we imposed equation (8) on equation (6), and we have(11)DtκCθζ,t=1Preff∂2θζ,t∂ζ2+Qθζ,t.We prefer to apply Laplace transform given in (8) on equation (11). The resultant form of the above expression is(12)∂2θ¯ζ,p∂ζ2−Preffpκ−Qθ¯ζ,p=0.The solution of equation (12) is(13)θ¯ζ,p=c1eζPreffpκ−Q+c2e−ζPreffpκ−Q,and using equation (7), we find out the arbitrary parameter:(14)θ¯ζ,p=1−e−pp2e−ζPreffpκ−Q.


### 3.2. Temperature Profile via the Caputo–Fabrizio Approach

Generating equation (6) for the fractional form, we imposed equation (9) on equation (6), and we have(15)p1−κp+κθ¯ζ,p=1Preff∂2θ¯ζ,p∂ζ2+Qθ¯ζ,p,(16)∂2θ¯ζ,p∂ζ2−Preffb1p+b2b3p+b4θ¯ζ,p=0.The solution of equation (16) is(17)θ¯ζ,p=c1eζPreffb1p+b2/b3p+b4+c2e−ζPreffb1p+b2/b3p+b4,and using equation (7), we find out the arbitrary parameter:(18)θ¯ζ,p=1−e−pp2e−ζPreffb1p+b2/b3p+b4.


### 3.3. Temperature Profile via the Atangana–Baleanu Approach

Generating equation (6) for the fractional form, we imposed equation (10) on equation (6), and we have(19)pκ1−κpκ+κθ¯ζ,p=1Preff∂2θ¯ζ,p∂ζ2+Qθ¯ζ,p,(20)∂2θ¯ζ,p∂ζ2−Preffb1pκ+b2b3pκ+b4θ¯ζ,p=0.The solution of equation (20) is(21)θ¯ζ,p=c1eζPreffb1pκ+b2/b3pκ+b4+c2e−ζPreffb1pκ+b2/b3pκ+b4,and using equation (7), we find out the arbitrary parameter:(22)θ¯ζ,p=1−e−pp2e−ζPreffb1pκ+b2/b3pκ+b4,where b1=1−Q+κQ,b2=−κQ, and b3=1−κ,b4=κ.


## 3.1. Temperature Profile via the Caputo Approach

Generating equation (6) for the fractional form, we imposed equation (8) on equation (6), and we have(11)DtκCθζ,t=1Preff∂2θζ,t∂ζ2+Qθζ,t.We prefer to apply Laplace transform given in (8) on equation (11). The resultant form of the above expression is(12)∂2θ¯ζ,p∂ζ2−Preffpκ−Qθ¯ζ,p=0.The solution of equation (12) is(13)θ¯ζ,p=c1eζPreffpκ−Q+c2e−ζPreffpκ−Q,and using equation (7), we find out the arbitrary parameter:(14)θ¯ζ,p=1−e−pp2e−ζPreffpκ−Q.


## 3.2. Temperature Profile via the Caputo–Fabrizio Approach

Generating equation (6) for the fractional form, we imposed equation (9) on equation (6), and we have(15)p1−κp+κθ¯ζ,p=1Preff∂2θ¯ζ,p∂ζ2+Qθ¯ζ,p,(16)∂2θ¯ζ,p∂ζ2−Preffb1p+b2b3p+b4θ¯ζ,p=0.The solution of equation (16) is(17)θ¯ζ,p=c1eζPreffb1p+b2/b3p+b4+c2e−ζPreffb1p+b2/b3p+b4,and using equation (7), we find out the arbitrary parameter:(18)θ¯ζ,p=1−e−pp2e−ζPreffb1p+b2/b3p+b4.


## 3.3. Temperature Profile via the Atangana–Baleanu Approach

Generating equation (6) for the fractional form, we imposed equation (10) on equation (6), and we have(19)pκ1−κpκ+κθ¯ζ,p=1Preff∂2θ¯ζ,p∂ζ2+Qθ¯ζ,p,(20)∂2θ¯ζ,p∂ζ2−Preffb1pκ+b2b3pκ+b4θ¯ζ,p=0.The solution of equation (20) is(21)θ¯ζ,p=c1eζPreffb1pκ+b2/b3pκ+b4+c2e−ζPreffb1pκ+b2/b3pκ+b4,and using equation (7), we find out the arbitrary parameter:(22)θ¯ζ,p=1−e−pp2e−ζPreffb1pκ+b2/b3pκ+b4,where b1=1−Q+κQ,b2=−κQ, and b3=1−κ,b4=κ.


## 4. Velocity Profile via Fractional Operators

### 4.1. Velocity Profile via the Caputo Approach

We utilize Laplace transformation for the solutions of the velocity profile given by equation (5):(23)a1+λDtκC∂Vζ,t∂t=1+λrDtγC∂2Vζ,t∂ζ2+1+λDtκCGrTζ,t−a2Vζ,t.We prefer to apply Laplace transform given in (8) on equation (23). The resultant form of the above expression is(24)a1+λpκpV¯ζ,p=1+λrpγ∂2V¯ζ,p∂ζ2+1+λpκGrθ¯ζ,p−a2V¯ζ,p.The solution of the homogeneous part of (24) is(25)V¯ζ,p=c1eζa1+λpκp+a2/1+λrpγ+c2e−ζa1+λpκp+a2/1+λrpγ.The required solution can be written as(26)V¯ζ,p=c1eζa1+λpκp+a2/1+λrpγ+c2e−ζa1+λpκp+a2/1+λrpγ−Gr1+λpκ1−e−pe−ζPreffpκ−Qp2Preffpκ−Q1+λrpγ−a1+λpκp+a2,(27)V¯ζ,p=1−e−pp2e−ζa1+λpκp+a2/1+λrpγ−Gr1+λpκ1−e−pp2Preffpκ−Q1+λrpγ−a1+λpκp+a2×e−ζa1+λpκp+a2/1+λrpγ−e−ζPreffpκ−Q.


### 4.2. Velocity Profile via the Caputo–Fabrizio Approach

We utilize Laplace transformation for the solutions of the velocity profile given by equation (5). We prefer to apply Laplace transform (9) on equation (5). The resultant form of the above expression is(28)a1+λp1−κp+κpV¯ζ,p=1+λrp1−γp+γ∂2V¯ζ,p∂ζ2+1+λp1−κp+κGrθ¯ζ,p−a2V¯ζ,p.The simplest form of the above equation is as follows:(29)V¯ζ,p=c1eζp+a6a1p+a8p+a2/p+a4a7p+a6+c2e−ζp+a6a1p+a8p+a2/p+a4a7p+a6−Grp+a6a9p+a41−e−pp2Preffa3−Qp+Qa4a7p+a6−p+a6a1p+a8p+a2×e−ζPreffa3−Qp+Qa4/p+a4,and using (7), we find out the arbitrary parameter:(30)V¯ζ,p=1−e−pp2e−ζp+a6a1p+a8p+a2/p+a4a7p+a6−Grp+a6a9p+a41−e−pp2Preffa3−Qp+Qa4a7p+a6−p+a6a1p+a8p+a2×e−ζPreffa3−Qp+Qa4/p+a4−e−ζPreffa3−Qp+Qa4/p+a4,where a1=1+λM+λr/K,a2=1+M/K, a3=1/1−κ, a4=κ/1−κ, a5=1/1−γ,a6=γ/1−γ,a7=1+a5λr,a8=a1a4+λa3, and a9=1+λa3.


### 4.3. Velocity Profile via the Atangana–Baleanu Approach

We utilize Laplace transformation for the solutions of the velocity profile given by equation (5). We prefer to apply Laplace transform (10) on equation (5). The resultant form of the above expression is(31)a1+λpκ1−κpκ+κpV¯ζ,p=1+λrpγ1−γpγ+γ∂2V¯ζ,p∂ζ2+1+λpκ1−κpκ+κGrθ¯ζ,p−a2V¯ζ,p,V¯ζ,p=c1eζpγ+a6a1pκ+a8p+a2/pκ+a4a7pγ+a6+c2e−ζpγ+a6a1pκ+a8p+a2/pκ+a4a7pγ+a6−Grpκ+a6a9pγ+a41−e−pp2Preffa3−Qpκ+Qa4a7pγ+a6−pγ+a6a1pκ+a8p+a2×e−ζPreffa3−Qpκ+Qa4/pκ+a4,and using (7), we find out the arbitrary parameter:(32)V¯ζ,p=1−e−pp2e−ζpκ+a4a7pγ+a6/pκ+a4a7pγ+a6−Grpκ+a6a9pγ+a41−e−pp2Preffa3−Qpκ+Qa4a7pγ+a6−pγ+a6a1pκ+a8p+a2×e−ζpγ+a6a1pκ+a8p+a2/pκ+a4a7pγ+a6−e−ζPreffa3−Qpκ+Qa4/pκ+a4.Asκ⟶1, in the required velocity expressions (27), (30), and (32), we get the same result for the classical model as discussed in [22]. Furthermore, if we neglect λ1=0 and λ2=0, then the results are identical which were obtained by Riaz et al. [35]. This shows the validation of our obtained results. We use classical computational and numerical techniques such as Stehfest’s [41] and Tzou’s algorithms [42] for the inverse of Laplace transform. Tzou’s calculation for our numerical inverse Laplace is(33)vr,t=e4.7t12v¯r,4.7t+Re∑k=1N1−1κv¯r,4.7+kπit,where Re⋅ is the real part, i represents the imaginary part, and N1 is the natural number.


## 4.1. Velocity Profile via the Caputo Approach

We utilize Laplace transformation for the solutions of the velocity profile given by equation (5):(23)a1+λDtκC∂Vζ,t∂t=1+λrDtγC∂2Vζ,t∂ζ2+1+λDtκCGrTζ,t−a2Vζ,t.We prefer to apply Laplace transform given in (8) on equation (23). The resultant form of the above expression is(24)a1+λpκpV¯ζ,p=1+λrpγ∂2V¯ζ,p∂ζ2+1+λpκGrθ¯ζ,p−a2V¯ζ,p.The solution of the homogeneous part of (24) is(25)V¯ζ,p=c1eζa1+λpκp+a2/1+λrpγ+c2e−ζa1+λpκp+a2/1+λrpγ.The required solution can be written as(26)V¯ζ,p=c1eζa1+λpκp+a2/1+λrpγ+c2e−ζa1+λpκp+a2/1+λrpγ−Gr1+λpκ1−e−pe−ζPreffpκ−Qp2Preffpκ−Q1+λrpγ−a1+λpκp+a2,(27)V¯ζ,p=1−e−pp2e−ζa1+λpκp+a2/1+λrpγ−Gr1+λpκ1−e−pp2Preffpκ−Q1+λrpγ−a1+λpκp+a2×e−ζa1+λpκp+a2/1+λrpγ−e−ζPreffpκ−Q.


## 4.2. Velocity Profile via the Caputo–Fabrizio Approach

We utilize Laplace transformation for the solutions of the velocity profile given by equation (5). We prefer to apply Laplace transform (9) on equation (5). The resultant form of the above expression is(28)a1+λp1−κp+κpV¯ζ,p=1+λrp1−γp+γ∂2V¯ζ,p∂ζ2+1+λp1−κp+κGrθ¯ζ,p−a2V¯ζ,p.The simplest form of the above equation is as follows:(29)V¯ζ,p=c1eζp+a6a1p+a8p+a2/p+a4a7p+a6+c2e−ζp+a6a1p+a8p+a2/p+a4a7p+a6−Grp+a6a9p+a41−e−pp2Preffa3−Qp+Qa4a7p+a6−p+a6a1p+a8p+a2×e−ζPreffa3−Qp+Qa4/p+a4,and using (7), we find out the arbitrary parameter:(30)V¯ζ,p=1−e−pp2e−ζp+a6a1p+a8p+a2/p+a4a7p+a6−Grp+a6a9p+a41−e−pp2Preffa3−Qp+Qa4a7p+a6−p+a6a1p+a8p+a2×e−ζPreffa3−Qp+Qa4/p+a4−e−ζPreffa3−Qp+Qa4/p+a4,where a1=1+λM+λr/K,a2=1+M/K, a3=1/1−κ, a4=κ/1−κ, a5=1/1−γ,a6=γ/1−γ,a7=1+a5λr,a8=a1a4+λa3, and a9=1+λa3.


## 4.3. Velocity Profile via the Atangana–Baleanu Approach

We utilize Laplace transformation for the solutions of the velocity profile given by equation (5). We prefer to apply Laplace transform (10) on equation (5). The resultant form of the above expression is(31)a1+λpκ1−κpκ+κpV¯ζ,p=1+λrpγ1−γpγ+γ∂2V¯ζ,p∂ζ2+1+λpκ1−κpκ+κGrθ¯ζ,p−a2V¯ζ,p,V¯ζ,p=c1eζpγ+a6a1pκ+a8p+a2/pκ+a4a7pγ+a6+c2e−ζpγ+a6a1pκ+a8p+a2/pκ+a4a7pγ+a6−Grpκ+a6a9pγ+a41−e−pp2Preffa3−Qpκ+Qa4a7pγ+a6−pγ+a6a1pκ+a8p+a2×e−ζPreffa3−Qpκ+Qa4/pκ+a4,and using (7), we find out the arbitrary parameter:(32)V¯ζ,p=1−e−pp2e−ζpκ+a4a7pγ+a6/pκ+a4a7pγ+a6−Grpκ+a6a9pγ+a41−e−pp2Preffa3−Qpκ+Qa4a7pγ+a6−pγ+a6a1pκ+a8p+a2×e−ζpγ+a6a1pκ+a8p+a2/pκ+a4a7pγ+a6−e−ζPreffa3−Qpκ+Qa4/pκ+a4.Asκ⟶1, in the required velocity expressions (27), (30), and (32), we get the same result for the classical model as discussed in [22]. Furthermore, if we neglect λ1=0 and λ2=0, then the results are identical which were obtained by Riaz et al. [35]. This shows the validation of our obtained results. We use classical computational and numerical techniques such as Stehfest’s [41] and Tzou’s algorithms [42] for the inverse of Laplace transform. Tzou’s calculation for our numerical inverse Laplace is(33)vr,t=e4.7t12v¯r,4.7t+Re∑k=1N1−1κv¯r,4.7+kπit,where Re⋅ is the real part, i represents the imaginary part, and N1 is the natural number.


## 5. Results and Discussion

This section is dedicated to present physical interpretation of the obtained results via C, CF, and AB differential operators on the MHD fractional Oldroyd-B fluid over an infinite vertical plate on the porous medium. Analytical results are investigated via Laplace transformation with the inversion algorithm for velocity and energy profiles. The graphical representations are depicted for showing the influences of different physical parameters on velocity and temperature using the package of Mathcad-15. The time influence on all fractional derivative operators is analyzed in Figure2. It clearly shows that, for the variation of time, the behavior of the velocity profile is the same. The resultant velocity of the ABC model is larger than the other fractional models.Effect ofκ: the influence of fraction parameter κ on velocity can been seen through Figure 3. Clearly, fluid velocity reduces with the increase in the fractional parameter for small and large time. It is worth mentioning that profiles for these are best to explain the history (memory) of the fluids. While making comparison, velocity for the Atangana–Baleanu model is greatest because it has a nonlocal kernel. Velocity for CF is greater than C. This is because CF has a nonsingular kernel that imitates C with the singular kernel.Effect ofM: Figure 4 investigates the impact of the magnetic force on all fractional operators. This graphical representation indicates that, with an increase in the magnetic field, the velocity reduces due to Lorentz force. By increasing the parameter of M, the Lorentz force also increases. Fluid flow on the boundary layer slows down due to this force.Effect ofGr: Figure 5 shows the impact on Gr for the velocity field versus time. It can be seen that the velocity field enhances by increasing Gr. It is supported by the physical fact that Gr is the fraction of buoyancy and viscous forces. An increase in Gr means that the buoyancy force gets stronger near the plate such that it overcomes the viscous force and that the fluid gets accelerated.Effect ofλ: Figure 6 shows the impact on the velocity field for λ. As λ increases, the thickness of the momentum boundary layer reduces which results in the deceleration of the fluid. As a relaxation time increment implies that the fluid will take extra time to calm, it readily justifies the decrease in velocity. It is quite a reverse behavior as compared to λr.Effect ofλr: Figure 7 shows the behavior of velocity curves for λr. It is observed that velocity enhances with the increase in λr for all fractional models. The velocity behavior is also observed for the variation of time.Effect ofPreff: Figure 8 discusses the effect of Preff using C, CF, and ABC models with the variation of time. Specific heat boundary thickness depends on Preff. The thickness of the momentum and boundary layer is controlled by Preff. It is seen from the graph that the decrease in the velocity is observed by the increase in the value of Preff.Effect ofκ: Figure 9 highlights the effect of the fractional parameter on the temperature profile for fractional models. With the increase in κ, the resultant temperature decreases. Temperature for CF and ABC is more as compared to C in all cases. Moreover, as κ⟶1, temperature curves for noninteger order approach integer order.Figure 2
Comparison of the velocity profile for C, CF, and ABC with the variation of time.
(a)(b)Figure 3
Plot via C, CF, and AB approaches for velocity with different values ofκ.
(a)(b)(c)(d)(e)(f)Figure 4
Plot via C, CF, and AB approaches for velocity with different values ofM.
(a)(b)(c)(d)(e)(f)Figure 5
Plot via C, CF, and AB approaches for velocity with different values ofGr.
(a)(b)(c)(d)(e)(f)Figure 6
Plot via C, CF, and AB approaches for velocity with different values ofλ.
(a)(b)(c)(d)(e)(f)Figure 7
Plot via C, CF, and AB approaches for velocity with different values ofλr.
(a)(b)(c)(d)(e)(f)Figure 8
Plot via C, CF, and AB approaches for velocity with different values ofPreff.
(a)(b)(c)(d)(e)(f)Figure 9
Plot via C, CF, and AB approaches for the temperature profile with different values ofκ.
(a)(b)(c)(d)(e)(f)


## 6. Conclusion

This paper studies MHD Oldroyd-B fluid flow with ramped wall temperature and velocity under the influence of the thermal radiation in a porous medium. Fractional derivative operators with the inversion algorithm are used to acquire the solution of velocity and temperature. The significant remarks for this article are as follows:(1)
Velocity curves show decreasing behavior for fractional parametersκ and M. The velocity field decreases by increasing the value of Pr.(2)
Velocity increases asGr increases for all fractional models.(3)
Velocity profile is a decreasing function ofPreff for all fractional operators.(4)
Velocity profile shows an opposite behavior forλ1 and λr for all fractional operators.(5)
Temperature decreases by enhancing the value of the fractional parameter.(6)
ABC fractional operator is more considerable as compared to all the other fractional operators.



---

*Source: 8890820-2020-12-07.xml*