Natural convection flow of a fluid using Atangana and Baleanu fractional model

A modified fractional model for the magnetohydrodynamic (MHD) flow of a fluid is developed utilizing Atangana–Baleanu fractional derivative (ABFD). Natural convection and wall oscillation instigate the flow over a vertical plate positioned in a porous medium. The partial differential equations (PDEs) are transmuted to ordinary differential equations (ODEs). The Laplace transform method with its inversion is employed to accomplish the exact solutions of momentum and heat equations. The final solution is expressed in terms of gamma function, modified Bessel function, and Mittag-Leffler function. The previous definitions Caputo fractional and Riemann–Liouville are rarely used by the researchers now due to their limitations. The newly introduced ABFD has got significance nowadays due to its nonlocal and nonsingular kernel. This work focuses on the oscillating boundary conditions for the viscous model in terms of ABFD. The influence of involved parameters is interpreted through plots. The velocity profile is an increasing function of fractional parameter and jumps for a higher Grashof number due to buoyancy push. Furthermore, the Atangana–Baleanu (AB) model is compared with the ordinary derivative model for limiting case and analyzed in detail. It is noted that the ordinary fluid flows faster compared to the fractional fluid.


Introduction
The literature of natural convection is sufficiently rich for magnetohydrodynamic oscillatory flow with classical fluid models. These natural problems related to the engineering and sciences lack the memory effect description. The fractional derivatives can describe memory effects, rheology, viscoelastic effect, diffusive transport, and fluid flow. In sciences and engineering, such as fluid mechanics, biomedical engineering, earthquake engineering, chemical engineering, and cooling processes in industries, many problems are modeled in fractional differential equations. Finding an operator that can physically describe this fact was complicated. Several scientists proposed that fractional derivatives can fully describe the memory effects, but the time and space components are the vital factors in implementing this idea. Various real-world problems follow three mathematical functions, i.e., the power law function, the exponential decay function, and the generalized Mittag-Leffler function [1]. These mentioned functions are the basis for many mathematical definitions of the fractional differential operators such as Caputo and Riemann-Liouville fractional derivatives based on the power law functions [1,2]. These two definitions have been followed by researchers, but they have a limitation that the Caputo fractional derivative has a singular kernel and the derivative of a constant is not zero Riemann-Liouville fractional derivative [3]. These limitations can be overcome by Caputo-Fabrizio fractional derivative but due to the locality in its kernel this definition also exhibit limitations [4,5]. Shedding some light on these two definitions, some quality literature work is worth discussing here to understand their importance. Baleanu and Fernandez [6] established a new formula consisting of Mittag-Leffler kernel in the form of a series of Riemann-Liouville fractional integrals. Several varieties of fractional definitions can be found in the literature. Fernandez et al. [7] expressed the Prabhakar fractional model and its generalized forms as a series of Riemann-Liouville integrals. In a recent study, a generalization for the existing definitions of fractional derivatives and integrals was established by Fernandez et al. [8]. They introduced an integral operator with a general kernel and expressed it as an infinite series of the Riemann-Liouville integrals. Recently, Jarad and Abdeljawad [9] established generalized noninteger derivatives and the Laplace transform to solve dynamical systems in the fractional derivatives frame. Jarad et al. [10] developed a new class of fractional operators in the Reimann-Liouville and Caputo sense. Some relevant references to the fractional integrals and their applications can be found in [11][12][13][14][15] and the references therein.
Azhar et al. [16] used the Caputo-Fabrizio time fractional derivative to the problem of nanofluid flow over a moving vertical plate. Fetecau et al. [17] conducted an analysis on the flow of a nanofluid over an isothermal plate Caputo-time fractional derivative. The flow of a generalized second grade fluid between parallel plates with Riemann-Liouville fractional derivative model was investigated by Wenchang and Mingyu [18]. They acquired the exact analytical solution using the Laplace transform and the Fourier transform. The flow of a second order fluid induced by a plate moving impulsively with fractional anomalous diffusion was investigated by Mingyu and Wenchang [19]. The Rayleigh-Stokes problem for a fractional second grade fluid was studied by Shen et al. [20]. The fractional Laplace transform and Fourier sine transform were employed to obtain the exact solution. Exact analytical solution unsteady flow of a generalized Maxwell fluid between two circular cylinders was determined via Laplace and Hankel transforms [21]. Recently Shen et al. [22] studied a fractional Maxwell viscoelastic nanofluid for various particle shapes. Caputo time fractional derivative was implemented by Zhang et al. [23] to acquire the numerical and analytical solutions for the problem of 2D flow of Maxwell fluid under variable pressure gradient. They used the separation of variables method to acquire analytical solution while for numerical solution the finite difference method was used. Aman et al. [24] studied a fractional Maxwell fluid for exact analytical solution with second order slip. Now to tackle the above discussed issue of local kernel, Atangana and Baleanu [25] proposed a new fractional derivative with nonlocal and nonsingular kernels, which fulfills all the mathematical principles related to fractional calculus. The ABFD definition is based on the Mittag-Leffler function and promises an improved description of dynamics of the system with the memory effects [26]. ABFD in the Caputo sense is used to investigate the model of vertical transmission and to cure a vector-borne disease [27]. They found the model very effective in curing the disease with a decrease in transmission rate. Very limited research can be found on fluid mechanics problems reporting fluid flow, heat transfer, and other world problems in applied mathematics and mathematical physics. Gomez-Aguilar et al. [28] presented a three-dimensional cancer model using Caputo-Fabrizio-Caputo and Atangana-Baleanu-Caputo fractional order derivatives. This model was proposed to study the evolution of cancer growth by the interactions of tumor cells with healthy cells and activated immune system cells. A scheme via the Laplace transform, the Sumudu-Picard integration method, and the Adams-Moulton rule was used to acquire the solutions.
Jan et al. [29] determined the solution for Brinkman type nanofluid using an Atangana-Baleanu fractional model. Owolabi and Atangana [30] analyzed the numerical simulation of an Adams-Bashforth scheme using Atangana-Baleanu-Caputo fractional derivatives. Saad et al. [31] established numerical solutions for the fractional Fisher type equation using ABFD. They employed a spectral collocation method based on Chebyshev approximations. In this research work, the spectral collocation method was implemented for the first time to solve nonlinear equation with Atangana-Baleanu derivatives. The logarithmic-KdV equation with the AB fractional operator and Mittag-Leffler type kernel was investigated by Mustafa et al. [32] using the numerical iterative method. Some plenteous literature regarding ABFD, their analytical solution, and applications can be found in [33][34][35][36] and the references therein. A few previous studies with ABFD are referenced in [34,35,37,38], but mostly they are focused on non-Newtonian fluids, nanofluids, or some are formulated with constant wall velocity. The research work carried out on nanofluids connects to an entirely different formulation and outcomes as they emphasize more the thermophysical properties of the nanoparticles and fluid.
Recently research with huge practical significance was carried out on human immunodeficiency virus-tuberculosis (HIV-TB) coinfected model using ABFD and the Mittag-Leffler law by Khan et al. [39]. They carried out the simulations for the existence, uniqueness of solution, and Hyers-Ulam (HU) stability. Another interesting research was done using ABFD for human blood flow nanofluids by Khan [40]. Their target was the investigation of nanofluids with channel flow in terms of the ABFD model and via the Laplace transform technique. They considered constant wall velocity at the boundaries. There is a gap for these combinations of studies in terms of oscillating boundary conditions for the ABFD model. Thus, the proposed study will focus on the ABFD model of an oscillatory Newtonian fluid flow with magnetic and radiative heat flux.
The literature reveals the importance of ABFD in the fluid dynamics domain. There is limited literature that discusses the free convection fluid flow using an Atangana-Baleanu derivative model for oscillatory flow; therefore the authors are motivated to conduct this research. Exact solutions are to be established for the problem of natural convection flow of fluid over a vertical plate using ABFD. The analytical solution is traced out via the Laplace transform method with its inversion. The solution is particularized for limiting case, i.e., for an ordinary fluid model by making α = 1. Finally, a brief graphical interpretation has been made for various parameters and comparison with the ordinary fluid model.
The structure of the paper is arranged as follows. In Sect. 2 we derive the mathematical ABFD model for the problem under consideration. The analytical solution is given in Sect. 3, we also recover the limiting case here. The acquired results are explained briefly in Sect. 4 followed by the conclusion of the paper in Sect. 5.

Mathematical framework
Consider unsteady flow of a fluid over a vertical plate with heat flux taken into consideration. A uniform magnetic field is employed perpendicular to the flow of the fluid. The flow is induced by the natural buoyancy force and heat flux q r on the plate. The temperature on the surface of the plate is T w , while T ∞ is the temperature far away from the surface. The plate starts oscillations at time t > 0. The physical interpretation of the flow is shown in Fig. 1. The governing equations of momentum and energy are given as follows: Here, T and u denote temperature and velocity, ρ f , μ, k 1 , B 2 0 , σ f , β f , k f , q r are the density, dynamic viscosity, porosity parameter, magnetic parameter, electrical conductivity, thermal expansion coefficient, thermal conductivity, and radiative flux parameter. The radiative heat flux is defined as in [16,24]: The boundary conditions are where U 0 is the amplitude and ω is the frequency of oscillations. The similarity variables are as follows: (1)-(4) to dimensionless PDEs using Eq. (5), we get (asterisk* is omitted for convenience): with the boundary conditions where where k, M, Gr T , Pr, N , ν, α denote the porosity, magnetic, Grashof number, Prandtl number, radiation parameter, kinematic viscosity, and mean radiation absorption coefficient respectively. Here ν = μ ρ f . The ABFD model of the present problem is given by Pr · AB D α t θ = Atangana and Baleanu fractional derivative [25] is defined as where N(α) is a normalization function such that is a Mittag-Leffler function [41]. The Laplace transform of Eq. (13) is given as follows: where q denotes the Laplace transform operator.

Limiting case
Taking α = 1 in Eq. (14), we acquire the limiting case by recovering the solution for the classical model of the present problem: Thus, the solution we recovered is given as follows:

Results and discussion
This  Fig. 2(a), the temperature profile for N = 0.5 is low, while in Fig. 2(b) the temperature shows a unique rising behavior at α = 0.7 and α = 0.9 for greater radiation parameter, i.e., N = 2. Figure 3 Fig. 4(a) with variation of α. A noteworthy increasing behavior can be seen for Gr = 5.5 in Fig. 4(b) due to a huge push of buoyancy forces. All these observations are accomplished at t = 2.   Fig. 6(a) with variation of k. The increasing behavior for α = 0.9 in Fig. 6(b) is considerable compared to that in case of smaller α due to a huge push of buoyancy forces. All these observations are carried out at t = 2.
In Fig. 7  Fluid velocity tends to decrease gradually in Figs. 7(a) and 7(b), but the effect is more significant with higher velocity range in case of a greater value of α.
In Fig. 8 In Fig. 9, the temperature profile of the current fractional problem is compared with the ordinary model of the problem by making α = 1 in the heat equation (Eq. (2)). The deviation at some points is clear at t = 1 with constant parameters Pr = 5, N = 1.5. It is worth noting that both ordinary and fractional models have qualitatively identical behavior of fluid temperature. Figure 10 depicts the comparison for velocity profile of fractional derivative    Table 1, corresponding to the results in Figs. 9-10. These numerical values are evaluated for different distances of fluid from the plate, i.e., 0 ≤ y ≤ 2. In consistency with these plots, the numerical values signify the abrupt changes when α = 1, i.e., the fractional fluid model is replaced by the ordinary one. For the fractional fluid, the highest value of velocity u(y, t) is 0.422, while for the ordinary fluid it is 0.173.

Conclusion
In this section, the concluding remarks for the analysis of a fluid under free convection and radiation are emphasized. The ABFD model has been extended for oscillatory flow and heat transfer model for which the analytical solution has been established in terms of Mittag-Leffler function via the attractive technique "Laplace transform method". Limiting case is recovered for an ordinary fluid when fractional parameter is taken to be one, i.e., α = 1. The main outcomes and concluding remarks are given as follows. Fluid velocity increases with increasing effect of the radiation parameter. But for the same constant parameters and varying radiation parameter, the velocity acquires a significant increase for greater value of the Grashof number. Fluid flows faster when the value of porosity parameter increases. The more permeable the medium is, the faster the fluid flows. The maximizing range of fractional parameter gives a unique jump to the velocity profile when the Grashof number is as high as 5.5 or above. A similar behavior is found in case of temperature for greater value of the radiation parameter. For the ordinary fluid, the temperature profile has identical behavior to the fractional model. The velocity profile is higher for the ordinary fluid compared to the fractional fluid. × exp -trur α cos απ sin ur α sin απ dr du, ( A 4 ) , L -1 q q 2 + a 2 = cos(at). (A5)