Numerical solvability of generalized Bagley–Torvik fractional models under Caputo–Fabrizio derivative

This paper deals with the generalized Bagley–Torvik equation based on the concept of the Caputo–Fabrizio fractional derivative using a modified reproducing kernel Hilbert space treatment. The generalized Bagley–Torvik equation is studied along with initial and boundary conditions to investigate numerical solution in the Caputo–Fabrizio sense. Regarding the generalized Bagley–Torvik equation with initial conditions, in order to have a better approach and lower cost, we reformulate the issue as a system of fractional differential equations while preserving the second type of these equations. Reproducing kernel functions are established to construct an orthogonal system used to formulate the analytical and approximate solutions of both equations in the appropriate Hilbert spaces. The feasibility of the proposed method and the effect of the novel derivative with the nonsingular kernel were verified by listing and treating several numerical examples with the required accuracy and speed. From a numerical point of view, the results obtained indicate the accuracy, efficiency, and reliability of the proposed method in solving various real life problems.


Introduction
In the last two decades, fractional calculus has become one of the most effective tools used in modeling of dynamical systems, to name a few, quantum, quantum physics, liquids, mechanics, optimization of biological models [1][2][3][4][5][6], where these phenomena are modeled and updated in equations that are closest to actually describing their condition. The first application for fractional calculus was probably what is referred to as a tautochrone problem that was given by Niels Henrik Abel in 1823 [2]. The great expansion in this field caused the attention and concern of researchers from various sciences, who in turn rushed to improve and correct some of the deficiencies and gaps in modeling.
However, there is no generally accepted definition of fractional operator, but some papers suggested mathematical properties that a fractional derivative must have. Hence, we can find many definitions each has advantages and disadvantages, see for example [5][6][7][8].
In 2015, the scientists Caputo and Fabrizio introduced a novel operator without a singular kernel [9], called Caputo-Fabrizio fractional derivative (CFFD). The advantage of this operator is that it provides an accurate update and description of some typical phenomena in which the most famous definitions often complain, similar to the Caputo, Riemann-Liouville, Riesz fractional operators, of limitations and deficiencies in the description [10][11][12][13][14][15], since all of these definitions have a singular kernel. The interest of the nonsingular kernel is due to the necessity of describing material heterogeneitics and structures with distinct scales which local theory fails to describe. For more information, the reader is kindly requested to read the papers [16][17][18][19][20][21].
Anyhow, Caputo-Fabrizio derivative has a nice characteristic; it has an equivalent representation with singular kernel which allows to study more general materials such as visco-plastic materials rather than only visco-elastic ones. The reader can refer to [22] in which the equivalence of CFFD with a model of singular kernel is presented. The Bagley-Torvik equation is one of the most important equations that occupies a leading status in the study of applied sciences and engineering applications. It was originally formulated in the eighties of the last century, where both Bagley and Torvik presented a study model viscoelasticity damped structure in the following equation [23]: along with the constants a i ∈ R, i = 1, 2, 3, h(ξ ) ∈ C[0, X] which is the space of continuous functions over the interval [0, X], and z(ξ ) is an unknown function to be determined. Moreover, the Bagley-Torvik equation (1.1) has been generalized as γ ∈ (0, 2), and it can be modeled from a generalized fractional vibration equation [24]. Most previously prepared works are mostly related to Riemann-Liouville and Caputo derivatives. Furthermore, the coefficients a i may also change parallel with the changes in fluid density and viscosity. That is, a i may be functions with respect to time. Exclusively, we generally consider the following generalized Bagley-Torvic equation (GBT) within CFFD: with initial conditions given by whereas the boundary conditions are of the form Physically, the fractional 1/2-order and 3/2-order derivatives are common to predict the frequency-dependent damping material. These derivatives are often used to model the motion of physical systems, as a study of an immersed plate in a Newtonian fluid was associated with the derivative 3/2 and gas in a fluid with the derivative 1/2. For more information, we refer to [25][26][27][28][29][30] and the references therein.
The value and importance of the Bagley-Torvic equation attracted the interest of researchers in finding approximate solutions using a range of iterative methods, among them, Adomian decomposition method, the generalized Taylor collocation method, homotopy analysis method, Chebyshev wavelet method, multistep methods, and predictorcorrector method of Adams type [31][32][33]. In this present paper, we are interested in studying and developing the modified reproducing kernel Hilbert space (MRKHS) under the influence of the CFFD, and making it suitable for solving the GBT under initial or boundary conditions, and that with greater accuracy and a lower time effort. The beginnings of the MRKHS method dates back to 1907, when Stanisław Zaremba presented research treating the boundary value problem of the harmonic functions. Thereafter, the process of developing this method proceeded until it has reached its current form. The MRKHS method occupies an important position among other numerical methods, as it is a very effective tool in many fields, such as machine learning, statistics, probability theory, economics, and the theory of integral operators [34][35][36][37][38]. In addition to that, this method is not limited to treating the well-posed problems, but also allows treating ill-posed problems.
The aim of this research paper can be summarized in several key points: • We undertake adjustments and improvements at the level of the MRKHS method, under the influence of the novel operator CFFD. • We use the CFFD properties to formulate GBT under initial conditions and convert it into an equivalent system of fractional differential equations. Then we apply the method given to solve this system. Finally, we get a relationship between the system solution and our equation. • We solve the GBT equation under the boundary conditions using the proposed method. • We prove some theorems related to the RKHS solution and its convergence to the exact solution. • We add mathematical simulations to determine the appropriateness and effectiveness of the accounts created.

Basic definitions and concepts
In this part, we present some basic definitions and theories for the Caputo-Fabrizio fractional derivative (CFFD) used in our study. It should be noted that rewriting Eq. (1.1) in the sense of this operator is due to its property of having a nonsingular kernel that enables it to represent some phenomena that well brief the functions of kernels. This data is not presented in other definitions, which creates difficulties in modeling. The use of exponentially based kernel and the equivalence of CFFD with a model of singular kernel give high precision in realistic representation phenomena, especially for physical problems.

Definition 2.2 ([9]
) Let z be a smooth function over [0, X] and 0 ≤ α ≤ 1. Then the corresponding fractional integral of CFFD operator with order α of a function z is given by The following proposition will be used as an important hint for the higher-order derivations.

Proposition 2.4 ([10])
If α ∈ [0, 1] and n ∈ N , then the CFFD of order n + α is defined as Remark 2.5 In this present work, we take M(α) = 1 and a = 0, so we can reformulate Definition 2.1 as follows: where ( * ) is the convolution operator.
In addition to that, the CFFD of any constant is zero, that is, for all c ≥ 0, we have CF D α ξ c = 0. Anyhow, we end this section with the equivalent representation of the CFFD which is

Preliminaries of MRKHS method
In this section, some of the essential facts of reproducing kernel theory are presented to construct Hilbert spaces associated with the reproducing function of our method. For more details, please read the papers [39][40][41][42][43][44][45][46][47]. During this study, AC[0, T] denotes the absolutely continuous real functions.
for each v ∈ W is called reproducing-kernel function, whereas the property in (ii) is called "the reproducing property". This function possesses some important properties such as being of unique representation, conjugate symmetric, and positive-definite.
The inner product and the norm corresponding in H m [0, X] for m = 1, 2, 3 are given as follows:

Theorem 3.5 ([35])
The unique representation of the reproducing-kernel function associated with H 2 [0, X] can be written as as well as the inner product and the norm associated with W m [0, X] are built as follows:

.1 MRKHS solution of GBT equation along with boundary conditions
In this section, we present a brief description for the notations and preliminary definitions of the MRKHS theory. Additionally, we explain how to solve GBT with boundary conditions (1.2)-(1.4) using the MRKHS method. Accordingly, we construct orthonormal function systems of the space based on the process orthogonalization of Gram-Schmidt.
Consider the GBT equation within CFFD Now, to apply our technique to GBT equation (4.1) on the Hilbert space H 2 [0, X], we consider a linear differential operator defined as follows: Thus, using the simple transform z(ξ ) := (ω(ξ ) -(μ 2μ 1 )ξ + μ 1 ), the GBT with the boundary condition equation (4.1) can be equivalently converted to the form Proof In order to prove that L is a bounded operator, it is enough to find M > 0 such that By the definition of the inner product for m = 1 in (3.2) on the Hilbert space On the other hand, using the reproducing property of the MRKHS and by the Cauchy-Schwarz inequality, we can write Hence, We construct an orthonormal function system of H 3 [0, X] as follows: can be generated from the wellknown Gram-Schmidt orthogonalization process as follows: where σ ik > 0 are the orthogonalization coefficients.
Proof Note that Now, for each z(ξ ) ∈ H 3 [0, X], suppose z(ξ ), i (ξ ) = 0, thus we have and L is continuous, we get Lz(ξ i ) = 0. Using the existence of the inverse operator L -1 , we conclude that z(ξ ) = 0. So the proof is complete.
is a dense subset on [0, X] and the analytic solution z(ξ ) of (4.1) is unique, then the solution z(ξ ) can be represented in the following form: (4.10) Now, the approximate solution can be obtained by truncating the series in (4.9) as follows: In the following theorem, we prove that the error that results when approximating the solution in (4.9) by the form (4.11) is decreasing to zero.

Theorem 4.4
Let E n = zz n H 3 , where z, z n are respectively the exact and the approximate solution of (4.1) represented in (4.9) and (4.11), then the error E n decreases monotonously in the sense of · H 3 , and E n → 0 as n → ∞.
Proof We have Note that E 2 n-1 > E 2 n , so we conclude that the error E n is monotone decreasing in the sense of · 2 H 3 , and because ∞ i=1 z(ξ ), i (ξ ) i (ξ ) is convergent, E n → 0 as n → ∞. Hence, the proof is complete.

MRKHS solution of GBT equation along with initial conditions
Now, we give some notations and preliminary definitions of the MRKHS theory. We then explain how to reformulate the GBT under the initial conditions (1.2)-(1.3) into an equivalent system of first-order fractional differential equations and how to implement our method to solve this system, with highlighting the relationship between the solution of the system and the solution of GBT. Accordingly, we construct an orthonormal function system of the space based on the Gram-Schmidt orthogonalization process.
Consequently, the initial value GBT equation can be converted into the form Applying the Schwarz inequality and using the fact that z, a 1 , a 2 are continuous over [0, X] and CF D 1 2 ξ V {2} ξ is continuous and uniformly bounded, we get and similarly, where ϒ ij ∈ + , ∀i, j = 1, 2. Hence, where M = max{M 1 , M 2 }. So, the proof is complete.
Next, we create an orthonormal function system of W 2 [0, X] as follows: Let ij (·) = V {1} ξ i (·) and ij (·) = L ij (·) for each i = 1, 2, . . . , and j = 1, 2, where L is the adjoint operator of L and {ξ i } ∞ i=1 is a countable dense subset of [0, X]. Moreover, using the properties of the reproducing kernel, we find To build the representative form of the MRKHS solutions of SFDE (4.14) equipped with the initial conditions (4.15) in the space W 2 [0, X], we use the well-known Gram-Schmidt process that outputs an orthonormal function , (4.18) where σ j ik , j = 1, 2, are the orthogonalization coefficients.
is a dense subset on [0, X] and the analytic solution Z(ξ ) of SFDE (4.14) is unique, then the analytic solution Z(ξ ) can be represented in the following form: .

(4.21)
Moreover, if we take finitely many terms in the series representation for the analytic solution Z(ξ ), we get directly the approximate solution SFDE (4.14), and it is given as the form .
Proof From Theorem 4.6, the proof is direct. Proof We have Note that E 2 n-1 > E 2 n . So, we conclude that the error E n is monotone decreasing in the sense of · 2 W 2 , and because ∞ i=1 Z(ξ ), i (ξ ) i (ξ ) is convergent, then E n → 0 as n → ∞. So, the proof is complete.

Numerical experiments
In order to demonstrate the effectiveness and eligibility of the proposed method for solving generalized Bagley-Torvik equations under the effect of CFFD derivative, we provide some illustrative examples and then solve them using MRKHS. In each example, we compare their results with the exact solution. All numerical computations are carried out using Mathematica 12.0 software package.
which can be rewritten along with homogeneous ICs as follows: where the exact solution is z(ξ ) = ξ (1ξ ).
In order to apply the MRKHS method with high efficiency, we make the substitution z 1 (ξ ) = v(ξ ) and z 2 (ξ ) = v (ξ ) to convert it into the equivalent system: Now, we apply the proposed method with n = 70 to get some numerical results as shown in Table 1. A comparison between the exact and the approximate curves is displayed in Fig. 1, which shows the accuracy of our method.
Applying the methodology described in this paper, with n = 90, some tabulated and graphical results are shown in Table 3 and Fig. 3.
Example 5.4 Consider the GBT equation equipped with the following ICs: whose the exact solution is z(ξ ) = 2(1ξ ).
The comparison between the exact and the approximate curves for n = 5 is displayed in Fig. 4, which shows high accuracy of this method.
(5.6)  In this example, we assume different values of the fractional CFFD orders to see the effect of this nonsingular kernel derivative to the GBT equations. Hence, a comparison between the approximate solutions for α ∈ {1, 2, 1.9, 1.8, 1.7, 1.6, 1.5} is shown in Table 4. By applying the methodology described in this paper for n = 55, some tabulated and graphical results are shown in Table 5 and Fig. 6.

Conclusion
In this paper, solutions of generalized Bagely-Torvik equations under the fractional derivative of a nonsingular kernel, the Caputo-Fabrizio derivative, are meaningfully discussed. These types of differential equations were solved under appropriate initial or boundary conditions. In the case of ICs, we prefer to convert the problem into a system of fractional differential equations seeking more efficiency and simplicity of the RKHS technique. Some theories related to the proposed method solutions, convergence, and error estimation have been proven. Numerical examples have been presented to reveal the