Next Article in Journal
Fractal Analysis on Pore Structure and Hydration of Magnesium Oxysulfate Cements by First Principle, Thermodynamic and Microstructure-Based Methods
Next Article in Special Issue
Numerical Exploration via Least Squares Estimation on Three Dimensional MHD Yield Exhibiting Nanofluid Model with Porous Stretching Boundaries
Previous Article in Journal
A Study of Coupled Systems of ψ-Hilfer Type Sequential Fractional Differential Equations with Integro-Multipoint Boundary Conditions
Previous Article in Special Issue
Thermophysical Investigation of Oldroyd-B Fluid with Functional Effects of Permeability: Memory Effect Study Using Non-Singular Kernel Derivative Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Magnetic Field with Parabolic Motion on Fractional Second Grade Fluid

1
Department of Science & Humanities, Lahore Campus, National University of Computer and Emerging Sciences, Lahore 54000, Pakistan
2
Department of Mathematics, University of Management and Technology, Lahore 54770, Pakistan
3
Department of Automation, Biomechanics and Mechatronics, Lodz University of Technology, 1/15 Stefanowskiego St., 90-924 Lodz, Poland
4
Department of Mathematics, Art and Science Faculty, Siirt University, Siirt 56100, Turkey
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(4), 163; https://doi.org/10.3390/fractalfract5040163
Submission received: 10 August 2021 / Revised: 23 September 2021 / Accepted: 24 September 2021 / Published: 11 October 2021
(This article belongs to the Special Issue Recent Advances in Computational Physics with Fractional Application)

Abstract

:
This paper is an analysis of the flow of magnetohydrodynamics (MHD) second grade fluid (SGF) under the influence of chemical reaction, heat generation/absorption, ramped temperature and concentration and thermodiffusion. The fluid was made to flow through a porous medium. It has been proven in many already-published articles that heat and mass transfer do not always follow the classical mechanics process that is known as memoryless process. Therefore, the model using classical differentiation based on the rate of change cannot really replicate such a dynamical process very accurately; thus, a different concept of differentiation is needed to capture such a process. Very recently, new classes of differential operators were introduced and have been recognized to be efficient in capturing processes following the power law, the decay law and the crossover behaviors. For the study of heat and mass transfer, we applied the newly introduced differential operators to model such flow. The equations for heat, mass and momentum are established in the terms of Caputo (C), Caputo–Fabrizio (CF) and Atangana–Baleanu in Caputo sense (ABC) fractional derivatives. The Laplace transform, inversion algorithm and convolution theorem were used to derive the exact and semi-analytical solutions for all cases. The obtained analytical solutions were plotted for different values of existing parameters. It is concluded that the fluid velocity shows increasing behavior for κ , G r and G m , while velocity decreases for P r and M. For K r , both velocity and concentration curves show decreasing behavior. Fluid flow accelerates under the influence of S r and R. Temperature and concentration profiles increase for S r and R. Moreover, the ABC fractional operator presents a larger memory effect than C and CF fractional operators.

1. Introduction

Over the past thirty years, fractional derivatives have fascinated multiple investigators as compared to classical derivatives. Moreover, fractional derivatives are more credible in mathematical modeling of real-world problems [1,2,3,4].
In classical calculus, derivatives and integrals are uniquely computed. A similar situation exists in the case of fractional integrals. For example, Samko et al. [5], Podlubny [6], Oldham and Spanier [7], and Miller and Ross [8] used a similar definition to compute the fractional integrals. However, the circumstances are complicated in fractional order derivatives (FODs) because several different competing definitions exist in the literature. For instance, a few of those approaches include the Riemann-Liouville, the Caputo, the Hadamard, the Marchaud, the Granwald–Letinkov, the Erdelyi–Kober, the Riesz–Feller, the Caputo–Fabrizio and the Atangana–Baleanu approach. These definitions coincide only with some particular cases. Among all these approaches to define fractional differentiation and fractional integration, the approach of Riemann–Liouville is significant. However, the approach of Riemann–Liouville does not properly address the physics of some fractional derivative initial-boundary value problems. Furthermore, this definition can exhibit the derivative of a constant function other than zero. To overcome this problem, Caputo proposed an alternate definition of FOD in 1967 [9], and it was used in fluid dynamics to explain the theory of viscoelasticity [10]. Recently, Caputo and Fabrizio [11] provided a modern definition of a non-integer order derivative including exponential function and the Atangana and Baleanu [12], based on Mittatag–Leffler function, which is a generalization of the exponential function.
The FOD inherits a nonlocal nature, so it is an excellent tool to obtain a better understanding of the hereditary properties of different processes and materials. Possibly, the first utilization of non-integer calculus in physical problems was noticed due to the work of Abel [13] while finding the solution of integer order equation, known as tautochrone problem. In this problem, the curve of an object (frictionless wire), lying in a vertical plane, was determined by using the operator D 0 1 2 and assuming the dependence of the time position not on the starting point. Bagley [14] presented the first PhD thesis on the applications of FC in viscoelasticity models. Recently, the applications of FC have been observed in psychology to determine the time variation of the emotions of mankind [15,16]. The applications of FOD problems can be seen in dynamics and control systems [17], marine sciences and wave dynamics [18,19,20,21], diffusion processes [22,23,24], solid mechanics [25,26,27], medical sciences [28,29,30,31] and many more [32,33,34,35].
Convective flow is a self-sustained flow with a temperature gradient. Convective flow and magnetic effect combined play a very important role in real value problems. The SGF with warmth-transferring porous medium was discussed by Tan and Masuoka [36]. Mixed MHD convection on an upright plate with permeable space was analyzed by Aldose et al. [37]. Rashidi et al. [38] also investigated the difference between two normal kinds of liquid stream between clear liquid and permeable medium. The authors of [39] investigated the precarious MHD stream of turning SGF past a swaying plate. Khan et al. [40] discussed the precise answers for quickened flow behavior of a pivoting SGF in a permeable space. Bilal et al. [41] and Ali et al. [42] analyzed the SGF with a swaying plate under different conditions. The literature shows more interest developed in the numerical and approximate solutions on the convection flow of SGF [43,44,45,46,47,48,49]. Exact solutions for viscoelastic SGF fluid have been investigated by researchers [50,51,52]. In 2010, the SGF with Laplace transform methodology over a wavering plate was explored by Nazar et al. [53]. Ali et al. [54] explored the MHD fluid with permeable space. The procedure of heat transfer is expressed by using momentum, energy and continuity equations with stress and heat flux. The heat transfer phenomenon consists of the fins of the heat exchanger, the tabulator inside the tube or plates and the convective physical transport phenomenon [55]. MHD free convection radiative stream with different conditions has been considered by researchers [56,57,58,59]. Ali et al. [42] considered SGF with porous surface and determined the exact solutions. Some recent and useful work has also been done for SGF with fractional differential operators [60,61,62,63,64,65]. Heat and mass transfer phenomena in nanofluids with a porous medium have been investigated by [66,67,68]. Some significant work in the field of nanoparticles and carbon nanotubes via fractional derivatives have been done by researchers [69,70,71]. Recently, Rehman et al. [72] discussed heat and mass transfer of MHD unsteady SGF in the presence of ramped conditions. Song et al. [73] used the definition of ABC in order to study SGF with exponential heating and Darcy’s law. Moreover, Riaz et al. explored MHD SGF with ramped conditions via special functions [74].
In the present work, we propose the mathematical modeling of fractional SGF with the help of Laplace transform and fractional operators. Moreover, solutions are acquired for momentum, heat and mass profiles. An inversion algorithm is used for graphical interpretation of fractional models.

2. Mathematical Modeling

We begin with SGF on an upright and unbounded plate, with the impact of a magnetic field having strength B 0 . The plate is perpendicular to η ˜ -axis and parallel to x ˜ -axis. At the start, the plate and fluid are not moving. T ˜ is fixed temperature and C ˜ is concentration at the surface. At the time τ ˜ > 0 , the temperature of the plate is either raised or lowered to T ˜ + ( T ˜ w T ˜ ) τ ˜ τ 0 , when τ ˜ τ 0 , and thereafter, for τ ˜ > t 0 , a constant temperature τ ˜ w is maintained and the level of mass transfer at the surface of the wall is either raised or lowered to C ˜ + ( C ˜ w C ˜ ) t ˜ τ 0 , when τ ˜ τ 0 , and thereafter, for τ ˜ > τ 0 is maintained at the constant surface concentration C ˜ w , respectively. The physical model of the problem can be given as follows in Figure 1 [75]. The detailed method to solve the problem is shown in Figure 2.
Governing equations for momentum, heat and mass are presented by [75]:
u ˜ τ ˜ = υ + γ ρ τ ˜ 2 u ˜ η ˜ 2 + g β T T ˜ g β T T ˜ + g β τ ˜ C ˜ g β C ˜ C ˜ σ B o 2 ρ u ˜ Φ k 1 ˜ υ + γ ρ τ ˜ u ˜ ,
ρ C p T ˜ τ ˜ = k 2 T ˜ η ˜ 2 q r ˜ η ˜ + Q 0 T ˜ Q 0 T ˜ ,
C ˜ η ˜ = D M 2 C ˜ η ˜ 2 + D T 2 T ˜ η ˜ 2 k 2 ˜ C ˜ k 2 ˜ C ˜ ,
and the imposed initial and boundary conditions are [75]:
τ ˜ 0 , u ˜ η ˜ , 0 = 0 , T ˜ η ˜ , 0 = T ˜ , C ˜ η ˜ , 0 = C ˜ , η 0 ,
τ ˜ > 0 , u ˜ η ˜ , τ ˜ = u o τ ˜ 2 , T ˜ ( 0 , τ ˜ ) = T ˜ + ( T ˜ w T ˜ ) τ ˜ τ 0 , 0 < τ ˜ τ 0 ; T ˜ w , τ ˜ > τ 0 ,
C ˜ ( 0 , τ ˜ ) = C ˜ + ( C ˜ w C ˜ ) τ ˜ τ 0 , 0 < τ ˜ τ 0 ; C ˜ w , τ ˜ > τ 0 , η ˜ = 0 ,
τ ˜ 0 , u ˜ η ˜ , τ ˜ 0 , T ˜ η ˜ , τ ˜ T ˜ , C ˜ η ˜ , τ ˜ C ˜ , η ˜ .
Figure 2. Flow chart of the method.
Figure 2. Flow chart of the method.
Fractalfract 05 00163 g002
For the simplification of Equations (1)–(7), we introduce the dimensionless variables given below:
ζ = u 0 τ ˜ 0 2 υ η ˜ , w = u ˜ τ ˜ 0 2 u 0 , θ = T ˜ T ˜ T ˜ w T ˜ , C = C ˜ C ˜ C ˜ w C ˜ , M 2 = σ B 0 2 τ 0 ρ , K r = τ 0 k 2 ˜ ,
S c = υ D M , G r = g β T ˜ T ˜ w T ˜ u ˜ 0 τ 0 , G m = g β C ˜ C ˜ w C ˜ u ˜ 0 τ 0 , τ 0 = υ u ˜ 0 2 1 5 , γ 1 = γ ρ υ τ 0 ,
S r = D T T ˜ w T ˜ υ C ˜ w C ˜ , P r = ρ υ C p k , R = 16 σ * T ˜ 3 3 k k * , 1 k 1 = u 0 υ τ 0 2 Φ k 1 ˜ , H = Q 0 υ τ 0 k ,
t = τ ˜ τ 0 2 h = 1 u 0 τ 0 k 1 ˜ , c = 1 + γ 1 h , b = M 2 + h .
Therefore, the dimensionless momentum, energy and mass equations are [75]
c w ζ , t t = 2 w ζ 2 + γ 1 3 w ζ , t t ζ 2 + G r θ ζ , t + G m C ζ , t b w ζ , t ,
P r θ ζ , t t = 1 + R 2 θ ζ , t ζ 2 + H θ ζ , t ,
C ζ , t t = 1 S c 2 C ζ , t ζ 2 + S r 2 θ ζ , t ζ 2 K r C ζ , t ,
with initial and boundary conditions
t 0 , w ζ , 0 = 0 , θ ζ , 0 = 0 , C ζ , 0 = 0 , ζ 0 ,
t > 0 , w ζ , t = t 2 , θ ( 0 , t ) = t 0 < t 1 ; 1 t > 1 , , C ( 0 , t ) = t 0 < t 1 ; 1 t > 1 , , ζ = 0 ,
t > 0 , w ζ , t 0 , θ ζ , t 0 , C ζ , t 0 , ζ .

3. Solution for Caputo Fractional Operator

The Caputo (C) fractional time derivative and its Laplace transform [76] are given below:
C D τ κ N ζ , τ = 1 Γ ( n κ ) 0 τ N ( n ) ( ξ ) ( τ ξ ) κ + 1 n d ξ ,
L C D τ κ N ζ , τ = s κ L N ( ζ , τ ) s κ 1 N ( ζ , 0 ) .

3.1. Temperature Field

The following Caputo derivative form of temperature Equation (13) is developed by using Equation (18):
P r C D t κ θ ( ζ , t ) = 1 + R 2 θ ζ 2 + H θ .
By implementing the Laplace transform on Equation (20), we obtain
2 θ ¯ c ( ζ , s ) ζ 2 P r 1 + R s κ H P r θ ¯ c ( ζ , s ) = 0 ,
where the homogenous solution of the above equation is
θ ¯ c ( ζ , s ) = c 1 e ζ P r 1 + R s κ H P r + c 2 e ζ P r 1 + R s κ H P r .
c 1 and c 2 can be determined by employing Equations (15)–(17), and the required solution is given below:
θ ¯ c ( ζ , s ) = 1 e s s 2 e ζ P r 1 + R s κ H P r .

3.2. Concentration Field

The following Caputo derivative form of concentration Equation (14) is developed by using Equation (18):
C D t κ C ( ζ , t ) = 1 S c 2 C ( ζ , t ) ζ 2 + S r 2 θ ( ζ , t ) ζ 2 K r C ( ζ , t ) .
By implementing the Laplace transform on Equation (24),
s κ + K r C ¯ c ( ζ , s ) = 1 S c 2 C ¯ c ( ζ , s ) ζ 2 + S r 2 θ ¯ c ( ζ , s ) ζ 2 ,
the homogenous solution of the above equation is
C ¯ h ( ζ , s ) = c 1 e ζ S c s κ ( 1 κ ) s + κ + K r + c 2 e ζ S c s κ ( 1 κ ) s + κ + K r ,
the particular solution of Equation (25) is
C ¯ p ( ζ , s ) = S c S r ( 1 e s ) P r H + H κ s κ H κ s 2 P r H ( 1 κ ) S c ( 1 + R ) S c K r ( 1 + R ) ( 1 κ ) s κ ( H κ + S c K r κ ( 1 + R ) ) e ζ P r 1 + R s κ ( 1 κ ) s + κ H P r ,
and the required solution is given below:
C ¯ c ( ζ , s ) = 1 e s s 2 A 1 e ζ S c s κ + K r + A 1 e ζ P r 1 + R s κ H P r .
where
A 1 = S c S r P r 1 + R s κ H P r 1 e s H 1 + R + S c K r × ( P r 1 + R S c ) s κ ( H 1 + R ) + S c K r + 1 s 2 1 ( H 1 + R + S c K r ) s κ P r 1 + R S c H 1 + R + S c K r .

3.3. Velocity Field

The following Caputo derivative form of velocity Equation (12) is developed by using Equation (18):
1 + h γ 1 ( C D t κ ) w ( ζ , t ) = 1 + γ 1 C D t κ 2 w ζ 2 + G r θ + G m C M 2 + h w .
By implementing the Laplace transform on Equation (30),
1 + γ 1 s κ 2 w ¯ c ( ζ , s ) ζ 2 c s κ + b w ¯ c ( ζ , s ) = G r θ ¯ c ( ζ , s ) + G m C ¯ c ( ζ , s ) ,
the homogenous solution of above equation is
w ¯ h ( ζ , s ) = c 1 e ζ 1 + γ 1 k 1 s κ ( 1 κ ) s + κ + M 2 + k 1 + c 2 e ζ 1 + γ 1 k 1 s κ ( 1 κ ) s + κ + M 2 + k 1 ,
the particular solution for Equation (31) is
w ¯ p ( ζ , s ) = G r ( 1 e s ) s κ + b 2 2 s 2 ( a 4 s a 5 ) ( a 12 s + b 2 ) ( a 13 s + a 14 ) ( s + b 2 ) e ζ P r 1 + R s κ ( 1 κ ) s + κ H P r G m ( 1 e s ) s κ + b 2 2 s 2 ( a 7 s κ a 8 ) ( a 12 s + b 2 ) ( a 13 s κ + a 14 ) ( s κ + b 2 ) × S c S r P r H + H κ s κ H κ s 2 P r H ( 1 κ ) S c ( 1 + R ) S c K r ( 1 + R ) ( 1 κ ) s κ ( H κ + S c K r κ ( 1 + R ) ) 1 e s s 2 e ζ S c s κ ( 1 κ ) s κ + κ + K r G m ( 1 e s ) s κ + b 2 2 s 2 ( a 7 s κ a 8 ) ( a 12 s κ + b 2 ) ( a 13 s κ + a 14 ) ( s κ + b 2 ) × e ζ S c s κ ( 1 κ ) s + κ + K r e ζ P r 1 + R s κ ( 1 κ ) s + κ H P r
and the required solution is given below:
w ¯ c ( ζ , s ) = 2 s 3 e ζ c s κ + b 1 + γ 1 s κ + G r ( 1 e s ) s 2 ( P r 1 + R ) ( s κ H P r ) ( 1 + γ 1 s κ ) ( c s κ + b ) × e ζ c s κ + b 1 + γ 1 s κ e ζ ( P r 1 + R ) ( s κ H P r ) + G m S c s κ + K r 1 + γ 1 s κ c s κ + b × 1 e s s 2 e ζ c s κ + b 1 + γ 1 s κ 1 e s s 2 e ζ S c s κ + K r + A 1 A 2 ,
where
A 2 = e ζ S c ( s κ + K r ) e ζ ( P r 1 + R ) ( s κ H P r ) , b 1 = 1 1 κ , b 2 = κ b 1 , a 1 = 1 + R P r , a 2 = H P r , a 3 = 1 a 1 , a 4 = a 3 b 1 a 3 a 2 , a 5 = a 3 a 2 b 2 , b 4 = a 3 S c , b 5 = a 3 a 2 + S c K r , a 6 = b 5 b 4 , a 7 = S c ( b 1 + K r ) , a 8 = S c K r b 2 , a 9 = a 4 a 7 , a 10 = a 5 + a 8 , a 11 = a 10 b 9 , a 12 = γ 1 κ + 1 , a 13 = c 1 κ + b , a 14 = b κ 1 κ + b .

4. Solution for Caputo–Fabrizio Fractional Operator

The CF derivative and its Laplace transform [77] are given below:
C F D τ κ N ζ , τ = 1 1 κ 0 τ e x p κ ( τ ξ ) 1 κ N / ( ξ ) d ξ , 0 < κ < 1 ,
L C F D τ κ N ζ , τ = s L N ( ζ , τ ) N ( ζ , 0 ) ( 1 κ ) s + κ .

4.1. Temperature Field

The following CF derivative form of temperature Equation (13) is developed by using Equation (36):
P r C F D t κ θ ( ζ , t ) = 1 + R 2 θ ζ 2 + H θ .
By implementing the Laplace transform on Equation (38),
2 θ ¯ c f ( ζ , s ) ζ 2 P r 1 + R s ( 1 κ ) s + κ H P r θ ¯ c f ( ζ , s ) = 0 ,
the homogenous solution of above equation is
θ ¯ c f ( ζ , s ) = c 1 e ζ P r 1 + R s ( 1 κ ) s + κ H P r + c 2 e ζ P r 1 + R s ( 1 κ ) s + κ H P r ,
and c 1 and c 2 can be determined by employing Equations (15)–(17). The required solution is given below:
θ ¯ c f ( ζ , s ) = 1 e s s 2 e ζ B 2 .

4.2. Concentration Field

The following CF derivative form of concentration Equation (13) is developed by using Equation (36):
C F D t κ C ( ζ , t ) = 1 S c 2 C ( ζ , t ) ζ 2 + S r 2 θ ( ζ , t ) ζ 2 K r C ( ζ , t ) .
By implementing the Laplace transform on Equation (42),
s ( 1 κ ) s + κ + K r C c f ( ζ , s ) = 1 S c 2 C ¯ c f ( ζ , s ) ζ 2 + S r 2 θ ¯ c f ( ζ , s ) ζ 2 ,
the homogenous solution of the above equation is
C ¯ h ( ζ , s ) = c 1 e ζ S c s ( 1 κ ) s + κ + K r + c 2 e ζ S c s ( 1 κ ) s + κ + K r ,
the particular solution of Equation (44) is
C ¯ p ( ζ , s ) = S c S r ( 1 e s ) P r H + H κ s H κ s 2 P r H ( 1 κ ) S c ( 1 + R ) S c K r ( 1 + R ) ( 1 κ ) s ( H κ + S c K r κ ( 1 + R ) ) e ζ P r 1 + R s ( 1 κ ) s + κ H P r ,
and the required solution is given below:
C ¯ c f ( ζ , s ) = 1 e s s 2 A 3 e ζ B 3 + A 3 e ζ B 2 ,
where
A 3 = S c S r P r 1 + R 1 1 κ H P r s H 1 + R κ 1 κ 1 e s H 1 + R + S c K r κ 1 κ × P r 1 + R 1 1 κ H P r S c 1 1 κ + K r s ( H 1 + R + S c K r ) κ 1 κ + 1 s 2 A 4 ,
A 4 = 1 ( H 1 + R + S c K r ) κ 1 κ s P r 1 + R 1 1 κ H P r H 1 + R + S c K r κ 1 κ ,
B 3 = S c 1 1 κ + K r s + S c K r κ 1 κ s + κ 1 κ ,
B 2 = P r 1 + R 1 1 κ H P r s H 1 + R κ 1 κ s + κ 1 κ .

4.3. Velocity Field

The following CF derivative form of velocity Equation (12) is developed by using Equation (36):
1 + h γ 1 ( C F D t κ ) w ( ζ , t ) = 1 + γ 1 C F D t κ 2 w ζ 2 + G r θ + G m C M 2 + h w .
By implementing the Laplace transform on Equation (51),
c 1 κ + b s + b κ 1 κ s + κ 1 κ w ¯ c f ( ζ , s ) = γ 1 1 κ + 1 s + κ 1 κ s + κ 1 κ 2 w ¯ c f ( ζ , s ) ζ 2 + G r θ ¯ c f ( ζ , s ) + G m C ¯ c f ( ζ , s ) ,
the homogenous solution of the above equation is
w ¯ h ( ζ , s ) = c 1 e ζ 1 + γ 1 k 1 s ( 1 κ ) s + κ + M 2 + k 1 + c 2 e ζ 1 + γ 1 k 1 s ( 1 κ ) s + κ + M 2 + k 1 .
The particular solution of Equation (52) is
w ¯ p ( ζ , s ) = G r ( 1 e s ) s + b 2 2 s 2 ( a 4 s a 5 ) ( a 12 s + b 2 ) ( a 13 s + a 14 ) ( s + b 2 ) e ζ P r 1 + R s ( 1 κ ) s + κ H P r G m ( 1 e s ) s + b 2 2 s 2 ( a 7 s a 8 ) ( a 12 s + b 2 ) ( a 13 s + a 14 ) ( s + b 2 ) × S c S r P r H + H κ s H κ s 2 P r H ( 1 κ ) S c ( 1 + R ) S c K r ( 1 + R ) ( 1 κ ) s ( H κ + S c K r κ ( 1 + R ) ) 1 e s s 2 e ζ S c s ( 1 κ ) s + κ + K r G m ( 1 e s ) s + b 2 2 s 2 ( a 7 s a 8 ) ( a 12 s + b 2 ) ( a 13 s + a 14 ) ( s + b 2 ) × e ζ S c s ( 1 κ ) s + κ + K r e ζ P r 1 + R s ( 1 κ ) s + κ H P r
and the required solution is given below:
w ¯ c f ( ζ , s ) = 2 s 3 e ζ B 1 + G r ( 1 e s ) s + κ 1 κ 2 s 2 A 7 e ζ B 1 e ζ B 2 + G m s + κ 1 κ 2 s 2 A 8 e ζ B 1 e ζ B 3 + A 5 e ζ B 1 e ζ B 2
where
A 7 = P r 1 + R 1 1 κ H P r s H 1 + R κ 1 κ γ 1 1 + κ + 1 s + κ 1 κ c 1 κ + b s + b κ 1 κ s + κ 1 κ
A 8 = S c 1 1 κ K r s + S c K r κ 1 κ γ 1 1 + κ + 1 s + κ 1 κ c 1 κ + b s + b κ 1 κ s + κ 1 κ
B 1 = c 1 κ + b s + b κ 1 κ γ 1 1 κ + 1 s + κ 1 κ

5. Solution for Atangana–Baleanu Fractional Operator

The ABC derivative and its Laplace transform [78] are given below:
A B C D τ α f ( ξ , τ ) = 1 1 α 0 τ E α α ( t τ ) α 1 α f ( ξ , τ ) τ d τ ,
L A B C D t κ f ξ , τ = s κ L f ( ξ , τ ) s κ 1 f ( ξ , 0 ) ( 1 κ ) s κ + κ .

5.1. Temperature Field

The following ABC derivative form of temperature Equation (13) is developed by using Equation (59):
P r A B C D t κ θ ( ζ , t ) = 1 + R 2 θ ζ 2 + H θ ,
The Laplace transform of Equation (61) is
2 θ ¯ a b c ( ζ , s ) ζ 2 P r 1 + R s κ ( 1 κ ) s κ + κ H P r θ ¯ a b c ( ζ , s ) = 0 ,
the homogenous solution of above equation is
θ ¯ a b c ( ζ , s ) = c 1 e ζ P r 1 + R s κ ( 1 κ ) s κ + κ H P r + c 2 e ζ P r 1 + R s κ ( 1 κ ) s κ + κ H P r .
and c 1 and c 2 can be determined by employing Equations (15)–(17). The required solution is given below:
θ ¯ a b c ( ζ , s ) = 1 e s s 2 e ζ B 2 .

5.2. Concentration Field

The following ABC derivative form of concentration Equation (14) is developed by using Equation (59):
A B C D t κ C ( ζ , t ) = 1 S c 2 C ( ζ , t ) ζ 2 + S r 2 θ ( ζ , t ) ζ 2 K r C ( ζ , t ) .
The Laplace transform of Equation (65) is
s κ ( 1 κ ) s κ + κ + K r C ¯ a b c ( ζ , s ) = 1 S c 2 C ¯ a b c ( ζ , s ) ζ 2 + S r 2 θ ¯ a b c ( ζ , s ) ζ 2 ,
and the homogenous solution of the above equation is
C ¯ h ( ζ , s ) = c 1 e ζ S c s κ ( 1 κ ) s κ + κ + K r + c 2 e ζ S c s κ ( 1 κ ) s κ + κ + K r .
The particular solution of Equation (66) is
C ¯ p ( ζ , s ) = S c S r ( 1 e s ) P r H + H κ s κ H κ s 2 P r H ( 1 κ ) S c ( 1 + R ) S c K r ( 1 + R ) ( 1 κ ) s κ ( H κ + S c K r κ ( 1 + R ) ) e ζ P r 1 + R s κ ( 1 κ ) s κ + κ H P r ,
and the required solution is given below:
C ¯ a b c ( ζ , s ) = 1 e s s 2 A 5 e ζ B 6 + A 5 e ζ B 5 ,
where
A 5 = S c S r P r 1 + R 1 1 κ H P r s κ H 1 + R κ 1 κ 1 e s H 1 + R + S c K r κ 1 κ × P r 1 + R 1 1 κ H P r S c 1 1 κ + K r s κ ( H 1 + R + S c K r ) κ 1 κ + 1 s 2 A 6 ,
A 6 = 1 ( H 1 + R + S c K r ) κ 1 κ s κ P r 1 + R 1 1 κ H P r H 1 + R + S c K r κ 1 κ ,
B 6 = S c 1 1 κ + K r s κ + S c K r κ 1 κ s κ + κ 1 κ ,
B 5 = P r 1 + R 1 1 κ H P r s κ H 1 + R κ 1 κ s κ + κ 1 κ .

5.3. Velocity Field

The following ABC derivative form of velocity Equation (12) is developed by using Equation (59):
1 + h γ 1 ( A B C D t κ ) w ( ζ , t ) = 1 + γ 1 A B C D t κ 2 w ζ 2 + G r θ + G m C M 2 + h w .
The Laplace transform of Equation (74) is
c 1 κ + b s κ + b κ 1 κ s κ + κ 1 κ w ¯ a b c ( ζ , s ) = γ 1 1 κ + 1 s κ + κ 1 κ s κ + κ 1 κ 2 w ¯ a b c ( ζ , s ) ζ 2 + G r θ ¯ a b c ( ζ , s ) + G m C ¯ a b c ( ζ , s ) ,
the homogenous solution of the above equation is
w ¯ h ( ζ , s ) = c 1 e ζ 1 + γ 1 k 1 s ( 1 κ ) s + κ + M 2 + k 1 + c 2 e ζ 1 + γ 1 k 1 s ( 1 κ ) s + κ + M 2 + k 1 ,
the particular solution of Equation (75) is
w ¯ p ( ζ , s ) = G r ( 1 e s ) s + b 2 2 s 2 ( a 4 s a 5 ) ( a 12 s + b 2 ) ( a 13 s + a 14 ) ( s + b 2 ) e ζ P r 1 + R s ( 1 κ ) s + κ H P r G m ( 1 e s ) s + b 2 2 s 2 ( a 7 s a 8 ) ( a 12 s + b 2 ) ( a 13 s + a 14 ) ( s + b 2 ) × S c S r P r H + H κ s H κ s 2 P r H ( 1 κ ) S c ( 1 + R ) S c K r ( 1 + R ) ( 1 κ ) s ( H κ + S c K r κ ( 1 + R ) ) 1 e s s 2 e ζ S c s ( 1 κ ) s + κ + K r G m ( 1 e s ) s + b 2 2 s 2 ( a 7 s a 8 ) ( a 12 s + b 2 ) ( a 13 s + a 14 ) ( s + b 2 ) × e ζ S c s ( 1 κ ) s + κ + K r e ζ P r 1 + R s ( 1 κ ) s + κ H P r
and the required solution is given below
w ¯ a b c ( ζ , s ) = 2 s 3 e ζ B 4 + G r ( 1 e s ) s κ + κ 1 κ 2 s 2 A 9 e ζ B 4 e ζ B 5 + G m s κ + κ 1 κ 2 s 2 A 10 e ζ B 4 e ζ B 6 + A 6 e ζ B 4 e ζ B 5 .
where
A 9 = P r 1 + R 1 1 κ H P r s κ H 1 + R κ 1 κ γ 1 1 + κ + 1 s κ + κ 1 κ c 1 κ + b s κ + b κ 1 κ s κ + κ 1 κ ,
A 10 = S c 1 1 κ K r s κ + S c K r κ 1 κ γ 1 1 + κ + 1 s κ + κ 1 κ c 1 κ + b s κ + b κ 1 κ s κ + κ 1 κ ,
B 4 = c 1 κ + b s κ + b κ 1 κ γ 1 1 κ + 1 s κ + κ 1 κ .
As κ 1 in Equations (23), (41) and (64) for temperature, Equations (28), (46) and (69) for concentration and Equations (34), (55) and (78) for velocity, we recuperate results for temperature, concentration and velocity profile for integer order shown in Kataria and Hari (Equations (13)–(15)) [75], respectively.
Stehfest’s formula [79] is one of the simplest algorithms we use to sort out the inverse Laplace transform.
w ( r , t ) = e 4.7 t 1 2 w ¯ r , 4.7 t + R e k = 1 N 1 ( 1 ) k w ¯ r , 4.7 + k π i t ,
where Re(.) is the real part, i is the imaginary unit and N 1 is a natural number.

6. Results and Discussion

This article shows the effect of heat and mass transfer in MHD SGF past a vertical plate. Three fractional models C, CF and ABC for flow, energy and mass equations are presented. Fractional derivatives and Laplace transform are applied to examine solutions for non-dimension fractional models. The limiting cases of fractional models are discussed. The impact of several parameters on momentum, heat and mass profiles are compared and studied by graphs.
Figure 3 highlights the behavior of momentum curves for κ . We see the velocity accelerates by raising κ . The reason is by an increment in κ , the thickness of the boundary layer will enhance, so the velocity increases. Moreover, velocity is highest for the ABC.
Figure 4 reveals the deviation in velocity distribution under the MHD condition. As M increases, the frictional force rises and hence fluid velocity decreases. For different values of M, an increase in Lorentz force effectively decreases flow accelerating forces; as a result, velocity is decelerated. Fluid velocity is maximum, moderate and minimum for ABC, CF and C models, respectively.
Figure 5 analyzes the influence of G r on momentum profile. Physically, large values respond to significant buoyancy force as it is related to strong convection currents. As G r increases, all buoyancy forces are dominant frictional forces and the hence momentum profile becomes amplified.
To highlight the velocity behavior for G m , we present Figure 6. Physically, the increment in buoyancy forces reduces the viscous force that leads to augmenting the flow raise with higher values of G m . The velocity curves show maximum behavior for the ABC model as compared to the other two models.
Figure 3. Velocity curves corresponding to C, CF and ABC with variable κ where P r = 7 , S c = 0.66 , G r = 10 , G m = 5 , H = 3 , R = 5 , S r = 3 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 3. Velocity curves corresponding to C, CF and ABC with variable κ where P r = 7 , S c = 0.66 , G r = 10 , G m = 5 , H = 3 , R = 5 , S r = 3 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g003
Figure 4. Velocity curves corresponding to C, CF and ABC with variable M where P r = 7 , S c = 0.66 , G r = 10 , G m = 5 , H = 3 , R = 5 , S r = 3 , K r = 2 , κ = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 4. Velocity curves corresponding to C, CF and ABC with variable M where P r = 7 , S c = 0.66 , G r = 10 , G m = 5 , H = 3 , R = 5 , S r = 3 , K r = 2 , κ = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g004
Figure 5. Velocity curves corresponding to C, CF and ABC with variable G r where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , R = 5 , S r = 3 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 5. Velocity curves corresponding to C, CF and ABC with variable G r where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , R = 5 , S r = 3 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g005
Figure 6. Velocity curves corresponding to C, CF and ABC with variable G m where P r = 7 , S c = 0.66 , κ = 0.5 , G r = 10 , H = 3 , R = 5 , S r = 3 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 6. Velocity curves corresponding to C, CF and ABC with variable G m where P r = 7 , S c = 0.66 , κ = 0.5 , G r = 10 , H = 3 , R = 5 , S r = 3 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g006
The significant impact of K r on momentum and mass profiles is shown in Figure 7. Both the fluid flow and the concentration decay with the rise in the K r . The presence of chemical reaction reduces the buoyancy effects, which decreases and hence weakens the flow field. Clearly, the ABC model shows the highest velocity and concentration.
Figure 8 describes the behavior of momentum and mass profile with increasing values of S r . As the velocity profile increases, the momentum boundary layer becomes thicker. Physically, large values of S r respond with a significant increase in mass buoyancy force; as a result, momentum and mass profile are raised. The effect is greatest for the ABC fractional MHD SGF model.
By raising thermal radiation parameter, momentum and heat profiles accelerate as shown in Figure 9. As the thermal radiation parameter increases, heat generation through flow increases, and as a result, bonds between fluid particles split which causes fluid to flow fast.
Figure 10 depicts velocity curves corresponding to C, CF and ABC for κ = 0.6 . Clearly, the ABC model shows significant behavior as compared to the other two curves. The reason is that Atangana and Baleanu propounded an advanced fractional operator by utilizing the generalized Mittag–Leffler function as a non-local and non-singular kernel. Comparison of Nusselt number with ref. [72] at Pr = 0.71 is given in Table 1.
Figure 9. Velocity and temperature curves corresponding to C, CF and ABC with variable R where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , S r = 3 , G r = 10 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 9. Velocity and temperature curves corresponding to C, CF and ABC with variable R where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , S r = 3 , G r = 10 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g009
Figure 10. Velocity curves corresponding to fractional models with κ = 0.6 and P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , S r = 3 , G r = 10 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 10. Velocity curves corresponding to fractional models with κ = 0.6 and P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , S r = 3 , G r = 10 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g010

7. Conclusions

This article is about the study of SGF with radiation and chemical reaction. Three fractional operators are applied to establish momentum, heat and mass profiles. Laplace transform and Stehfest’s formula are utilized for the solutions of the mentioned equations. Several graphs are presented to illustrate the impact of incipient parameters for the solutions. Some main results are given below:
  • Velocity curves are increasing for greater values of κ , G r and G m .
  • Fluid flow descends for P r and M.
  • Velocity and concentration curves show a decreasing behavior under the influence of K r .
  • Fluid velocity accelerates under the impact of S r and R.
  • Heat and mass profiles for S r and R are show an increasing behavior.
  • Curves show prominent behavior for ABC among C, CF and ABC.
Extending the work in this article as suggested below will be an interesting endeavor.
We have studied MHD SGF in the presence of ramped concentration and temperature with thermodiffusion. This study can be further carried out by considering more complex OBF models. Additionally, this study can be carried out by a rotational fluid model. Moreover, the present work can be extended by selecting different non-dimensional parameters and quantities.

Author Contributions

Conceptualization, N.I.; methodology, J.A.; software, A.A.; validation, M.B.R. and J.A.; formal analysis, A.A.; investigation, J.A.; data curation, M.B.R.; writing—original draft preparation, N.I. and J.A; writing—review and editing, N.I. and J.A; visualization, A.A.; supervision, A.A.; project administration, M.B.R.; Final checking, M.B.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the Polish National Science Centre under the grant OPUS 18 No. 2019/35/B/ST8/00980.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

My Manuscript has no associated data.

Acknowledgments

This work has been supported by the Polish National Science Centre under the grant OPUS 18 No. 2019/35/B/ST8/00980.

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.

Nomenclature

SymbolQuantity
wVelocity of the fluid
θ Temperature of the fluid
CConcentration of the fluid
gAcceleration due to gravity
kThermal conductivity of the fluid
k 1 Permeability parameter
K r Parameter of chemical reaction
MParameter of magnetic field
Q 0 Coefficient of heat absorption/generation
P r Prandtl number
S c Schmidt number
G r Thermal Grashof number
G m Mass Grashof number
RParameter of thermal radiation
S r Soret number
D M Coefficient of mass diffusion
D T Coefficient of thermal diffusion
θ w Temperature of fluid at the plate
θ Temperature of fluid far away from the plate
C w Concentration level on the plate
C Concentration of the fluid far away from the plate
C p Specific heat at constant temperature
sLaplace transforms parameter
ρ Fluid density
κ Fractional parameter
γ One of the material modules of second grade fluids
γ 1 Second grade parameter
μ Dynamic viscosity
υ Kinematic viscosity
β T Volumetric coefficient of thermal expansion
β C Volumetric coefficient of expansion for mass concentration
Φ Porosity

References

  1. Bazhlekova, E.; Bazhlekov, I. Viscoelastic flows with fractional derivative models: Computational approach by convolutional calculus of Dimovski. Fract. Calc. Appl. Anal. 2014, 17, 954–976. [Google Scholar] [CrossRef]
  2. Duan, J.S.; Qiu, X. The periodic solution of Stokes’ second problem for viscoelastic fluids as characterized by a fractional constitutive equation. J. Non-Newton. Fluid Mech. 2014, 205, 11–15. [Google Scholar] [CrossRef]
  3. Friedrich, C.H.R. Relaxation and retardation functions of the Maxwell model with fractional derivatives. Rheo Acta 1991, 30, 151–158. [Google Scholar] [CrossRef]
  4. Meral, F.C.; Royston, T.J.; Magin, R. Fractional calculus in viscoelasticity: An experimental study. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 939–945. [Google Scholar] [CrossRef]
  5. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publish: Langhorne, PA, USA, 1993. [Google Scholar]
  6. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999; Volume 198. [Google Scholar]
  7. Oldham, K.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA; London, UK, 1974. [Google Scholar]
  8. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; John Wiley and Sons Inc.: New York, NY, USA, 1993. [Google Scholar]
  9. Caputo, M. Linear model of dissipation whose Q is almost frequency independent—II. Geophy. J. Int. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  10. Caputo, M. Elasticita e Dissipazione; Zanichelli: Bologna, Italy, 1969. [Google Scholar]
  11. Caputo, M.; Fabrizio, M. A new denition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 2015, 1, 1–13. [Google Scholar]
  12. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20. [Google Scholar] [CrossRef] [Green Version]
  13. Abel, N.H. Solution de Quelques Problemes Al’aide D’integrales Difinies, Oeuvres Completes; Grondahl: Christiania, Norway, 1881; Volume 1, pp. 16–18. [Google Scholar]
  14. Bagely, R.L. Applications of Generalized Derivatives to Viscoelasticity. Ph.D. Thesis, Air Force Institute of Technology, Kaduna, Nigeria, 1979. [Google Scholar]
  15. Ahmad, W.M.; El-Khazalib, R. Fractional order dynamical models of love. Chaos Solitons Fractals 2007, 33, 1367–1375. [Google Scholar] [CrossRef]
  16. Song, L.; Xu, S.; Yang, J. Dynamical models of happiness with fractional order. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 616–628. [Google Scholar] [CrossRef]
  17. Aksikas, I.; Fuxman, A.; Forbes, J.F.; Winkin, J. LQ control design of a class of hyperbolic PDE systems: Application to xed-bed reactor. Automatica 2009, 45, 1542–1548. [Google Scholar] [CrossRef]
  18. Arshad, M.S.; Mardan, S.A.; Riaz, M.B.; Altaf, S. Analysis of time-fractional semi-analytical solutions of strong interacting internal waves in rotating ocean. Punjab Univ. J. Math. 2020, 52, 99–111. [Google Scholar]
  19. Jhangeer, A.; Munawar, M.; Riaz, M.B.; Baleanu, D. Construction of traveling waves patterns of (1+n)-dimensional modified Zakharov-Kuznetsov equation in plasma physics. Results Phys. 2020, 19, 103330. [Google Scholar] [CrossRef]
  20. Lewis, P.A. A theory for the diffraction of the SH waves by randomly rough surfaces in two dimension. Q. J. Mech. Appl. Math. 1996, 49, 261–286. [Google Scholar] [CrossRef]
  21. Okrasinski, W. On a nonlinear convolution equation occuring in the theory of water percolation. Ann. Pol. Math. 1980, 37, 223–229. [Google Scholar] [CrossRef] [Green Version]
  22. Metzler, R.; Klafter, J. The random walks guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  23. Owolabi, K.M.; Atangana, A.; Akgul, A. Modelling and analysis of fractal fractional partial differential equations: Application to reaction-diffusion model. Alex. Eng. J. 2020, 59, 2477–2490. [Google Scholar] [CrossRef]
  24. Wyss, W. The fractional diffusion equation. J. Math. Phys. 1986, 27, 2782–2785. [Google Scholar] [CrossRef]
  25. Heibig, A.; Palade, L.I. On the rest state stability of an objective fractional derivative viscoelastic uid model. J. Math. Phys. 2008, 49, 043101. [Google Scholar] [CrossRef]
  26. Rosikin, Y.; Shitikova, M. Application of fractional calculus for dynamic problem of solid mechanics, novel trends and recent results. Appl. Mech. Rev. 2010, 63, 1–52. [Google Scholar] [CrossRef]
  27. Torvik, P.J.; Bagley, R.L. On the appearance of fractional derivatives in the behaviour of real materials. J. Appl. Mech. 1984, 51, 294–298. [Google Scholar] [CrossRef]
  28. Atangana, A. Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? Chaos Solitons Fractals 2020, 136, 109860. [Google Scholar] [CrossRef] [PubMed]
  29. Atangana, A.; Qureshi, S. Mathematical Modeling of an Autonomous Nonlinear Dynamical System for Malaria Transmission Using Caputo Derivative. Fract. Order Anal. Theory Methods Appl. 2020, 225–252. [Google Scholar] [CrossRef]
  30. Awrejcewicz, J.; Zafar, A.A.; Kudra, G.; Riaz, M.B. Theoretical study of the blood flow in arteries in the presence of magnetic particles and under periodic body acceleration. Chaos Solitons Fractals 2020, 140, 110204. [Google Scholar] [CrossRef]
  31. Sweilam, N.H.; Al-Mekhla, S.M.; Assiri, T.; Atangana, A. Optimal control for cancer treatment mathematical model using Atangana Baleanu Caputo fractional derivative. Adv. Diff. Equ. 2020, 2020, 334. [Google Scholar] [CrossRef]
  32. Abro, K.A.; Atangana, A. A comparative analysis of electromechanical model of piezoelectric actuator through Caputo Fabrizio and Atangana Baleanu fractional derivatives. Math. Methods Appl. Sci. 2020, 43, 9681–9691. [Google Scholar] [CrossRef]
  33. Arshad, M.S.; Baleanu, D.; Riaz, M.B.; Abbas, M. A Novel 2-Stage Fractional RungeKutta Method for a Time-Fractional Logistic Growth Model. Discret. Dyn. Nat. Soc. 2020, 2020, 1020472. [Google Scholar] [CrossRef]
  34. Asjad, M.I.; Aleem, M.; Riaz, M.B. Exact analysis of MHD Walters-B fluid flow with non-singular fractional derivatives of Caputo-Fabrizio in the presence of radiation and chemical reaction. J. Polym. Sci. Eng. 2018, 1, 599. [Google Scholar] [CrossRef]
  35. Zafar, A.A.; Kudra, G.; Awrejcewicz, J.; Abdeljawad, T.; Riaz, M.B. A comparative study of the fractional oscillators. Alex. Eng. J. 2020, 59, 2649–2676. [Google Scholar] [CrossRef]
  36. Tan, W.; Masuoka, T. Stoke’s first problem for a second grade fluid in a porous half-space with heated boundary. Int. J. Non-Linear Mech. 2005, 40, 515–522. [Google Scholar] [CrossRef]
  37. Aldoss, T.K.; Al-Nimr, M.A.; Jarrah, M.A.; Al-Shaer, B. Magnetohydrodynamics mixed convection from a vertical plate embedded in a porous medium. Numer. Heat Transf. Appl. 1995, 28, 635–645. [Google Scholar] [CrossRef]
  38. Rashidi, S.; Nouri-Borujerdi, A.; Valipour, M.S.; Ellahi, R.; Pop, I. Stress-jump and continuity interface conditions for a cylinder embedded in porous medium. Transp. Porous Media 2015, 107, 171–186. [Google Scholar] [CrossRef]
  39. Imran, M.A.; Imran, M.; Fetecau, C. MHD oscillating flows of rotating second grade fluid in a porous medium. Commun. Nonlinear Sci. Numer. Simul. 2014, 2014, 1–12. [Google Scholar] [CrossRef]
  40. Khan, I.; Farhad, A.; Norzieha, M. Exact solutions for accelerated flows of a rotating second grade fluid in a porous medium. World Appl. Sci. J. 2010, 9, 55–68. [Google Scholar]
  41. Riaz, M.B.; Zafar, A.A.; Vieru, D. On flows of generalized second grade fluids generated by an oscillating flat plate. In Matematica, Mecanica, Teoretica, Fizica; Buletinul Institutului Politehnic din Iaşi—Universitatea Tehnică: Iași, Romania, 2015; Volume 1. [Google Scholar]
  42. Riaz, M.B.; Awrejcewicz, J.; Rehman, A.U.; Akgül, A. Thermophysical Investigation of Oldroyd-B Fluid with Functional Effects of Permeability: Memory Effect Study Using Non-Singular Kernel Derivative Approach. Fractal Fract. 2021, 5, 124. [Google Scholar] [CrossRef]
  43. Hussain, M.; Hayat, T.; Asghar, S.; Fetecau, C. Oscillatory flows of second grade fluid in a porous space. Nonlinear Anal. Real World Appl. 2010, 11, 2403–2414. [Google Scholar] [CrossRef]
  44. Erdogan, M.E.; Imrak, C.E. An exact solution of the governing equation of a fluid of second grade for three dimensional vortex flow. Int. J. Eng. Sci. 2005, 43, 721–729. [Google Scholar] [CrossRef]
  45. Asghar, S.; Nadeem, S.; Hanif, K.; Hayat, T. Analytic solution of Stoke’s second problem for second gradefluid. Math Probl. Eng. 2006, 2006, 072468. [Google Scholar] [CrossRef]
  46. Tiwari, A.K.; Ravi, S.K. Analytical studies on transient rotating flow of a second grade fluid in a porous medium. Adv. Theor. Appl. Mech. 2009, 2, 33–41. [Google Scholar]
  47. Fetecau, C.; Fetecau, C. Starting solutions for the motion of second grade fluid. Int. J. Eng. Sci. 2005, 43, 781–789. [Google Scholar] [CrossRef]
  48. Khan, I.; Ellahi, R.; Fetecau, C. Some MHD flows of second grade fluid through the porous medium. J. Porous Media 2008, 11, 389–400. [Google Scholar]
  49. Fetecau, C.; Fetecau, C. Starting solutions for the motion of a second grade fluid due to longitudinal and torsional oscillations of a circular cylinder. Int. J. Eng. Sci. 2006, 44, 788–796. [Google Scholar] [CrossRef]
  50. Hsu, S.H.; Jamieson, A.M. Viscoelastic behavior at the thermal sol-gel transition of gelatin. Polymer 1993, 34, 2602–2608. [Google Scholar] [CrossRef]
  51. Bandelli, R. Unsteady unidirectional flows of second grade fluids in domains with heated boundaries. Int. J. Non-Linear Mech. 1995, 30, 263–269. [Google Scholar] [CrossRef]
  52. Damesh, R.A.; Shatnawi, A.S.; Chamkha, A.J.; Duwairi, H.M. Transient mixed convection flow of second grade viscoelastic fluid over a vertical surface. Nonlinear Anal. Model. Control. 2008, 13, 169–179. [Google Scholar] [CrossRef]
  53. Nazar, M.; Fetecau, C.; Vieru, D.; Fetecau, C. New exact solutions corresponding to the second problem of stokes’ for second grade fluids. Nonlinear Anal. Real World Appl. 2010, 11, 584–591. [Google Scholar] [CrossRef]
  54. Ali, F.; Norzieha, M.; Sharidan, S.; Khan, I.; Hayat, T. New exact solutions of stokes’ second problem for an MHD second grade fluid in a porous space. Int. J. Nonlinear Mech. 2012, 47, 521–525. [Google Scholar] [CrossRef]
  55. Makinde, O.D.; Khan, W.A.; Culham, J.R. MHD variable viscosity reacting flow over a convectively heated plate in a porous medium with thermophoresis and radiative heat transfer. Int. J. Heat Mass Transf. 2016, 93, 595–604. [Google Scholar] [CrossRef]
  56. Sheikholeslami, M.; Ganji, D.D.; Javed, M.Y.; Ellahi, R. Effect of thermal radiation on MHD nanofluid flow and heat transfer by means of two phase model. J. Magn. Magn. Mater. 2015, 374, 36–43. [Google Scholar] [CrossRef]
  57. Zhang, C.; Zheng, L.; Zhang, X.; Chen, G. MHD flow and radiation heat transfer of nanofluid in porous media with variable surface heat flux and chemical reaction. Appl. Math Model. 2015, 39, 165–181. [Google Scholar] [CrossRef]
  58. Rashidi, M.M.; Ali, M.; Freidoonimehr, N.; Rostami, B.; Hossain, M.A. Mixed convective heat transfer for MHD viscoelastic fluid flow over a porous wedge with thermal radiation. Adv. Mech. Eng. 2014, 6, 735939. [Google Scholar] [CrossRef] [Green Version]
  59. Dehghan, M.; Rahmani, Y.; Ganji, D.D.; Saedodin, S.; Valipour, M.S.; Rashidi, S. Convection-radiation heat transfer in solar heat exchangers filled with a porous medium: Homotopy perturbation method versus numerical analysis. Renew. Energy 2015, 74, 448–455. [Google Scholar] [CrossRef]
  60. Imran, M.A.; Shah, N.A.; Aleem, M.; Khan, I. Heat transfer analysis of fractional second-grade fluid subject to Newtonian heating with Caputo and Caputo-Fabrizio fractional derivatives: A comparison. Eur. Phys. J. Plus 2017, 132, 340. [Google Scholar]
  61. Tassaddiq, A. MHD flow of a fractional second grade fluid over an inclined heated plate. Chaos Solitons Fractals 2019, 123, 341–346. [Google Scholar] [CrossRef]
  62. Sene, N. Second-grade fluid model with Caputo–Liouville generalized fractional derivative. Chaos Solitons Fractals 2020, 133, 109631. [Google Scholar] [CrossRef]
  63. Haq, S.; Jan, S.; Jan, S.A.; Khan, I.; Singh, J. Heat and mass transfer of fractional second grade fluid with slippage and ramped wall temperature using Caputo-Fabrizio fractional derivative approach. AIMS Math. 2020, 5, 3056–3088. [Google Scholar] [CrossRef]
  64. Fatecau, C.; Zafar, A.A.; Vieru, D.; Awrejcewicz, J. Hydromagnetic flow over a moving plate of second grade fluids with time fractional derivatives having non-singular kernel. Chaos Solitons Fractals 2020, 130, 109454. [Google Scholar] [CrossRef]
  65. Siddique, I.; Tlili, I.; Bukhari, M.; Mahsud, Y. Heat transfer analysis in convective flows of fractional second grade fluids with Caputo–Fabrizio and Atangana–Baleanu derivative subject to Newtonion heating. Mech.-Time-Depend. Mater. 2019, 25, 291–311. [Google Scholar] [CrossRef]
  66. Mehryan, S.A.; Ghalambaz, M.; Vaezi, M.; Zadeh, S.M.; Sedaghatizadeh, N.; Younis, O.; Chamkha, A.J.; Abulkhair, H. Non-Newtonian phase change study of nano-enhanced n-octadecane comprising mesoporous silica in a porous medium. Appl. Math. Model. 2021, 97, 463–482. [Google Scholar] [CrossRef]
  67. Rana, B.M.; Arifuzzaman, S.M.; Islam, S.; Reza-E-Rabbi, S.; Al-Mamun, A.; Mazumder, M.; Roy, K.C.; Khan, M.S. Swimming of microbes in blood flow of nano-bioconvective Williamson fluid. Ther. Sci. Eng. Prog. 2021, 25, 101018. [Google Scholar] [CrossRef]
  68. Al-Mamun, A.; Arifuzzaman, S.M.; Reza-E-Rabbi, S.; Alam, U.S.; Islam, S.; Khan, M.S. Numerical simulation of periodic MHD casson nanofluid flow through porous stretching sheet. SN Appl. Sci. 2021, 3, 1–14. [Google Scholar] [CrossRef]
  69. Abro, K.A.; Gómez-Aguilar, J.F. Role of Fourier sine transform on the dynamical model of tensioned carbon nanotubes with fractional operator. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  70. Abro, K.A.; Laghari, M.H.; Gómez-Aguilar, J.F. Application of Atangana-Baleanu Fractional Derivative to Carbon Nanotubes Based Non-Newtonian Nanofluid: Applications in Nanotechnology. J. Appl. Comput. Mech. 2020, 6, 1260–1269. [Google Scholar] [CrossRef]
  71. Abro, K.A.; Khan, I.; Gomez-Aguilar, J.F. Heat transfer in magnetohydrodynamic free convection flow of generalized ferrofluid with magnetite nanoparticles. J. Therm. Anal. Calorim. 2021, 143, 3633–3642. [Google Scholar] [CrossRef]
  72. Rehman, A.U.; Riaz, M.B.; Akgül, A.; Saeed, S.T.; Baleanu, D. Heat and mass transport impact on MHD second-grade fluid: A comparative analysis of fractional operators. Heat Trans. 2021. [Google Scholar] [CrossRef]
  73. Song, Y.Q.; Raza, A.; Al-Khaled, K.; Farid, S.; Khan, M.I.; Khan, S.U.; Shi, Q.H.; Malik, M.Y.; Khan, M.I. Significances of exponential heating and Darcy’s law for second grade fluid flow over oscillating plate by using Atangana-Baleanu fractional derivatives. Case Stud. Therm. Eng. 2021, 27, 101266. [Google Scholar] [CrossRef]
  74. Riaz, M.B.; Abro, K.A.; Abualnaja, K.M.; Akgül, A.; Rehman, A.U.; Abbas, M.; Hamed, Y.S. Exact solutions involving special functions for unsteady convective flow of magnetohydrodynamic second grade fluid with ramped conditions. Adv. Diff. Equ. 2021, 2021, 408. [Google Scholar] [CrossRef]
  75. Kataria, H.R.; Patel, H.R. Effect of thermo-diffusion and parabolic motion on MHD Second grade fluid flow with ramped wall temperature and ramped surface concentration. Alex. Eng. J. 2018, 57, 173–185. [Google Scholar] [CrossRef] [Green Version]
  76. Riaz, M.B.; Atangana, A.; Iftikhar, N. Heat and mass transfer in Maxwell fluid in view of local and non-local differential operators. J. Therm. Anal. Calorim. 2020, 143, 4313–4329. [Google Scholar] [CrossRef]
  77. Iftikhar, N.; Baleanu, D.; Riaz, M.B.; Husnine, S.M. Heat and Mass Transfer of Natural Convective Flow with Slanted Magnetic Field via Fractional Operators. J. Appl. Comput. Mech. 2020, 7, 189–212. [Google Scholar] [CrossRef]
  78. Riaz, M.B.; Iftikhar, N. A comparative study of heat transfer analysis of MHD Maxwell fluid in view of local and non-local differential operators. Chaos Solitons Fractals 2020, 132, 109556. [Google Scholar] [CrossRef]
  79. Stehfest, H.A. Numerical inversion of Laplace transforms. In Proceedings of the Communications of the ACM, New York, NY, USA, 1 January 1970; Volume 13, pp. 9–47. [Google Scholar]
Figure 1. Geometrical presentation of the problem.
Figure 1. Geometrical presentation of the problem.
Fractalfract 05 00163 g001
Figure 7. Velocity and concentration curves corresponding to C, CF and ABC with variable K r where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , R = 5 , S r = 3 , κ = 0.5 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 7. Velocity and concentration curves corresponding to C, CF and ABC with variable K r where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , R = 5 , S r = 3 , κ = 0.5 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g007
Figure 8. Velocity and concentration curves corresponding to C, CF and ABC with variable S r where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , R = 5 , G r = 10 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Figure 8. Velocity and concentration curves corresponding to C, CF and ABC with variable S r where P r = 7 , S c = 0.66 , κ = 0.5 , G m = 5 , H = 3 , R = 5 , G r = 10 , K r = 2 , M = 0.5 , γ 1 = 0.1 , h = 0.5 , c = 1.05 and b = 0.75 .
Fractalfract 05 00163 g008
Table 1. Comparison of Nusselt number with ref. [72] at P r = 0.71.
Table 1. Comparison of Nusselt number with ref. [72] at P r = 0.71.
R ϕ t Nu (Ref. [75]) for Ramped Temp Nu (C) for Ramped Temp Nu (CF) for Ramped Temp Nu (ABC) for Ramped Temp
230.30.38370.3840.3850.386
230.50.55830.5570.5580.559
230.70.72890.7270.7280.729
230.50.44980.4470.4480.449
230.50.55830.5570.5580.559
250.50.65210.6530.6540.655
230.50.55830.5570.5580.559
430.50.43240.4330.4340.435
630.50.36550.3660.3670.368
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Iftikhar, N.; Riaz, M.B.; Awrejcewicz, J.; Akgül, A. Effect of Magnetic Field with Parabolic Motion on Fractional Second Grade Fluid. Fractal Fract. 2021, 5, 163. https://doi.org/10.3390/fractalfract5040163

AMA Style

Iftikhar N, Riaz MB, Awrejcewicz J, Akgül A. Effect of Magnetic Field with Parabolic Motion on Fractional Second Grade Fluid. Fractal and Fractional. 2021; 5(4):163. https://doi.org/10.3390/fractalfract5040163

Chicago/Turabian Style

Iftikhar, Nazish, Muhammad Bilal Riaz, Jan Awrejcewicz, and Ali Akgül. 2021. "Effect of Magnetic Field with Parabolic Motion on Fractional Second Grade Fluid" Fractal and Fractional 5, no. 4: 163. https://doi.org/10.3390/fractalfract5040163

APA Style

Iftikhar, N., Riaz, M. B., Awrejcewicz, J., & Akgül, A. (2021). Effect of Magnetic Field with Parabolic Motion on Fractional Second Grade Fluid. Fractal and Fractional, 5(4), 163. https://doi.org/10.3390/fractalfract5040163

Article Metrics

Back to TopTop