Abstract
A new approach for modeling real world problems called the “Eton Approach” was presented in this paper. The "Eton approach" combines both the concept of the variable order derivative together with Atangana derivative with memory derivative. The Atangana derivative with memory is used to account for the memory and fractional derivative for its filter effect. The approach was used to describe the potential energy field that is caused by a given charge or mass density distribution.We solve the modified model numerically and present supporting numerical simulations.
1 Introduction
The concept of memory effect has long been a problem within the modeling community. There are many physical problems that requested to be modeled taking into account the memory effect, however, the big challenge was to find a suitable operator that can efficiently describe the physically observed fact [1–6]. The approach was used to describe the potential energy field that is caused by a given charge or mass density distribution. Many researchers suggest that the memory effect could fully be described via fractional derivatives. Nonetheless, the idea of memory requests time and space components; for instance with the movement of pollution within a geological formations, one would like to know at what time and space the pollution was retarded in the aquifer. Indeed this cannot be described with fractional derivative. Therefore a new operator able to take care of this problem is needed. On the other hand, as shown in the literature, the local derivative could not explain many physical problems as compared with derivative with fractional order. It was recently suggested that since the derivative with fractional order are based upon the convolution design, the fractional derivative could therefore be considered as filter operators. These filters are able to get rid of all impurities introduced by the local derivative while modeling real world problems [7,8]. Thus in order to introduce the memory effect into the mathematical formulation of a given real world problem and also to filter all impurities introduced by local derivative, we propose in this paper a new approach that will combine the two derivatives.
2 New trends on fractional calculus
In this section we will present some useful information about the new derivative with fractional order. However, more information around this derivative can be found in the following works [7,8]
Definition 1. Let f ∈ H1 (a, b), b > a, a ∈ [0, 1] then, the new Caputo derivative of fractional derivative is defined as:
Where M(a) is a normalization function such that M(1) = 1[1]. However, if the function does not belongs to H1 (a, b) then, the derivative can be reformulated as
Remark 2. The authors remarked that, if
In addition,
Now after the introduction of a new derivative, the associate anti-derivative becomes important, the associated integral of the new Caputo derivative with fractional order was proposed by Nieto and Losada [2].
Definition 3. [2] Let 0 < a < 1. The fractional integral of order a of a function f is defined by
Remark [2]. Note that, according to the above definition, the fractional integral of Caputo type of function of order 0 < a < 1 is an average between function f and its integral of order one. This therefore imposes
The above expression yields an explicit formula for
Due to this, Nieto and Losada suggested that the new Caputo derivative of order 0 < α < 1 can be reformulated as
2.1 Some properties of the Caputo Fabrizio derivative
– For any natural number, n > 0, the following is obtained
– The Sumudu transform of a function f is given as:
Let f(t) be a function for which the Caputo-Fabrizio exists, then, the Sumudu transform of the Caputo-Frabrizio fractional derivative of f(t) is given as:
– If a function, f, is infinitely differentiable, the following result is obtained
2.2 Caputo-Fabrzio Derivative with fractional order in Riemann-Liouville(RL) sense
One problem faced by the Caputo-Fabrizio derivative with fractional order is that, when alpha goes to zero we do not recover the initial function but the function and addition term. To solve this problem Atangana and Goufo proposed [15] an alternative definition of fractional derivative with no singular kernel, the result was also obtained in the work done by Caputo and Fabrizio [16].
Definition 4. Let f be a function not necessary differential, let alpha be a real number such that 0 ≤ α ≤ 1, then the Caputo-Fabrizio derivative with order alpha in RL sense is given as:
If α = 0 is zero we have the following:
Using the argument by Caputo and Fabrizio, we also have that when α = 1 also goes to 1 we recover the first derivative.
3 Atangana derivative with memory
In mathematics, a dynamic scheme is a tuple (D, h, B) with D a manifold which can be a locally Banach space or Euclidean space, B the domain for time which is a set of non-negative real, and h is an evolution rule t → ft the range is of course a diffeomorphism of a manifold to itself [10].
Definition 5. Let D be a dynamic system with domain B (time or space), let u a positively defined function called uncertainty function of D, within the domain B, then if h ∈ D, the Atangana derivative with memory of a functionh denoted by Uu(f)is defined as:
Remark 6. If u = 1 we recover the first derivative (the the local derivative), if u = 0, we recover the initial function, this is conforms with the primary law of derivative.
3.1 Properties of Atangana derivative with memory
– Addition:
– Multiplication:
Proof.Using the definition of Atangana derivative with memory, we obtain the following:
This completes the proof. □
– Division:
– Lipchitz condition: Let f(t) and g(t) be two functions then
This proves that the Atangana derivative with memory possess the Lipchitz condition.
4 Poisson equation with Eton approach
In mathematics notably within applied mathematics, a partial differential equation is classified under elliptic type of equation and called Poisson’s equation is used in many fields of science for instance electrostatics, theoretical physics, and mechanical engineering. This equation is employed to describe the potential energy field caused by a given charge or mass density distribution. The aim of this section is to apply the new approach to the well-known Poisson equation, and present the difference between both models. We shall first recall that the one dimensional Poisson equation is given by:
Now, to apply the Eton approach of modeling real world problems, we apply first the Atangana derivative with memory as follows:
The above will first produce
Then using again the definition, we obtain:
For simplicity we have the following version
The second step is to replace the first order derivative with the variable order derivative to obtain
The above equation is the result of Poisson equation with Eton approach. In the next section, we shall present the existence and uniqueness of the above equation.
5 Existence and uniqueness of exact solution
In this section, we present in detail the analysis of existence of the solution of above equation. We start this investigation by converting Equation (25) into a Volterra equation by applying on both sides the anti-derivative of the C.F. derivative obtaining
The above equation can be reduced for simplicity to:
Since K is a linear operator, then we can find a positive constant, σ, such that
The above relation means the operator K satisfies Lipschitz condition. However assuming that, we can find two different solutions for Equation (27) then
Nonetheless using the Lipschitz property of K we obtain
Thus we have the following result
Nevertheless, we chose a specific, x0 > 0, such that
Therefore if
which implies that
This proves the uniqueness of the solution of the modified equation. We next present the existence of the solution, to do this we propose the following recursive formula
The solution of the time fractional with modified equation by the the "Eton approach" is given as
We shall now show that such solution. To achieve this, we shall first present the following preliminaries
Preliminaries: see([7,8,9,10]) Let us consider the following Banach space (X,||.|| ) be a Banach space and h a self-map of X. Let yn+1 = g(H, yn) be some iterative technique. Assuming that, F(H), the fixed point set of h has at least one element and that yn converges to a point p ∈ F(H). Let xn ⊆; X and define en = || xn+1 – g(H, xn)||. If
Theorem 7. (see [7,9]). Let (X,||.|| ) be a Banach space and h a self-map of X satisfying
|| Hx – Hy|| ≤ C||x – Hx|| + c||x – y||, for all x, y in X where 0 ≤ C, 0 ≤ α < 1. Suppose that H has fixed point p. Then, H is Picard H – Stable.
Theorem 8. Let h be a self-map that we chose to defined as follows
then the self-map is H-stable in L2( α, β) providing that the following inequality is satisfied,
Proof. For all n, m ∈ ℕ we have
The above can be converted using the triangular inequality property of the norm
Using the Lipschitz condition of K, we obtain the following
Rearranging, we obtain
Or we have
This shows that H-self-mapping has a fixed point and the proof is complete. □
In addition to this if we consider
We have the conditions of a hypothesis which holds for H. Thus, since all conditions of the hypothesis hold, then, H is Picard’s H-stable.
6 Numerical analysis
In this section, the new equation will be solved numerically using the well-known finite element method associate to the Caputo-Fabrizio derivative with fractional order. We will present the detail of stability and convergence of the numerical scheme. To accommodate readers that are not use to this new derivative, we present the first and second approximation of the Caputo-Fabrizio derivative [10]. The first approximation in space is given as:
where
With the second order approximation we derive the following:
At a particular point xm we have the following
Using the second order approximation of the ordinary derivative, we obtain
After integrating we obtained finally:
or in more simpler way:
Nevertheless replacing this in the modified equation obtained from the "Eton approach", we obtain
above recursive formula can therefore be used to generate numerical approximation of the new model of potential energy field caused by mass density distribution.
7 Numerical simulations
In this section we present some numerical simulations of the modified model for different function of uncertain derivative. For each function, we present the numerical representation for different values of alpha, the fractional order derivative. The uncertain order derivatives chosen here are
Additionally
8 Conclusion
The aim of this paper was to introduce a new approach of modeling real world problems. The presented approach combines the filter effect that is obtained by using the Caputo-Fabrizio derivative with fractional order and the memory effect, which is introduced by using the novel operator called uncertain derivative, this operator was proposed by Atangana. This approach is referred to as UFM or simply Eton approach. This novel approach was used to describe the potential energy field that is caused by a given charge or mass density distribution. The new approach help describe physical problems with memory.
Acknowledgement
The authors would like to extend their sincere appreciation to the Deanship of Scientific Research at King Saud University for funding this group No RG-1437-017.
References
1Coimbra C.F.M., Mechanics with variable-order differential operators, Annalen der Physik, 2003, 12, 692-703.10.1002/andp.200310032Search in Google Scholar
2Solomon T.H., Weeks E.R., Swinney H.L., Observation of anomalous diffusion and Lévy flights in a two-dimensional rotating flow, Phys. Rev. Lett., 1993, 71, 3975-3978.10.1103/PhysRevLett.71.3975Search in Google Scholar PubMed
3Magin R.L., Abdullah O., Baleanu D., Zhou X.J., Anomalous diffusion expressed through fractional order differential operators in the Bloch-Torrey equation, J. Magn. Reson., 2008, 190, 255-270.10.1016/j.jmr.2007.11.007Search in Google Scholar PubMed
4Sun H.G., Chen W., Chen Y.Q., Variable order fractional differential operators in anomalous diffusion modeling, Phys. A, 2009, 388, 4586-4592.10.1016/j.physa.2009.07.024Search in Google Scholar
5Azevedo E.N. et al., Concentration-dependent diffusivity and anomalous diffusion: a magnetic resonance imaging study of water ingress in porous zeolite, Phys. Rev. E, 2006, 73, Article ID 011204.10.1103/PhysRevE.73.011204Search in Google Scholar PubMed
6Pedro H.T.C., Kobayashi M.H., Pereira J.M.C., Coimbra C.F.M., Variable order modeling of diffusive-convective effects on the oscillatory flow past a sphere, J. Vib. Control, 2008, 14, 1659-1672.10.1177/1077546307087397Search in Google Scholar
7Caputo M., Fabrizio M., A new Definition of Fractional Derivative without Singular Kernel, Progr. Fract. Differ. Appl., 2015, 1, 73-85.Search in Google Scholar
8Losada J., Nieto J.J., Properties of a New Fractional Derivative without Singular Kernel, Progr. Fract. Differ. Appl., 2015, 1, 87-92.Search in Google Scholar
9Atangana A., On the new fractional derivative and application to nonlinear Fisher’s reaction-diffusion equation, Appl. Math. Comput., 2016, 273, 948-956.10.1016/j.amc.2015.10.021Search in Google Scholar
10Atangana A., Nieto J.J., Numerical solution for the model of RLC circuit via the fractional derivative without singular kernel, Adv. Mech. Eng., 2015, 7, 1-7.10.1177/1687814015613758Search in Google Scholar
11Odibat Z.M., Momani S., Application of variational iteration method to nonlinear differential equation of fractional order, Int. J. Nonlinear Sci. Numer. Simul., 2006, 7, 27-34.10.1515/IJNSNS.2006.7.1.27Search in Google Scholar
12Atangana A., Alkahtani B.S.T., Analysis of the Keller–Segel Model with a Fractional Derivative without Singular Kernel, Entropy, 2015, 17, 4439-4453.10.3390/e17064439Search in Google Scholar
13Benson D., Wheatcraft S., Meerschaert M., Application of a fractional advection-dispersion equation, Water Resources Res., 2015, 36, 1403-1412.10.1029/2000WR900031Search in Google Scholar
14Caputo M., Linear model of dissipation whose Q is almost frequency independent-II, Geophys. J. R. Astr. Soc., 1967, 13, 529-539.10.1111/j.1365-246X.1967.tb02303.xSearch in Google Scholar
15Doungmo Goufo E.M., Application of the Caputo-Fabrizio fractional derivative without singular kernel to Korteweg-de Vries-Bergers equation, Math. Model. Anal., 2015 (In press)10.3846/13926292.2016.1145607Search in Google Scholar
16Caputo M., Fabrizio M., Applications of new time and spatial fractional derivatives with exponential kernels, Progr. Fract. Differ. Appl., 2016, 2, 1-11.10.18576/pfda/020101Search in Google Scholar
© 2016 Badr Saad T. Alkahtani and A. Atangana, published by De Gruyter Open
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.