Abstract

The Coimbra concept of fractional order derivative is used to build a numerical approach using radial functions in this paper. The Coimbra derivative is capable of modelling a dynamic system with varying fractional order behaviour over time. The proposed scheme’s stability and convergence are investigated. In one and two space dimensions, the developed approach is validated for the given model. By applying a periodic boundary condition on a bounded domain, the model’s periodicity is shown statistically. The acquired findings demonstrate the new numerical scheme’s potency and, as a result, its high order accuracy.

1. Introduction

The Korteweg-De Vries (KdV) equation is first derived by Boussinesq in the year 1870. Later on in 1895, the same model was retrieved by Korteweg and de Vries [1] with the presumption of compact amplitude and huge wave length. In many nonlinear dispersive physical systems, the evolution of long wave can be expressed by the KdV type equation (see [25] and the references therein). In mathematical sciences and engineering, evolutionary nonlinear equations play a major role to model physical phenomena [2, 6]. In the theory of shallow water waves, the KdV equation is one of the most essential equations in nonlinear evolution developed in [4] and the references therein. Some of the important aspects of solutions of these dispersive equations discovered through observations are their long-time behaviour and known periodicity in time [7]. The important event of eventual periodicity has been presented previously in [8], and in more recent work [9, 10], a new solution is reestablished corresponding to the KdV equation. In addition, the forced oscillations and the stability of the KdV equation have been carried out in a very recent work [1114]. In applied mathematics, physics, and other related fields, a rich filed of research has evolved within the last century because of computational and analytic research on fractional and classical KdV equation [1519].

Both the theory and application of fractional calculus have advanced dramatically in the previous two decades. The nonlocal quality of fractional calculus and its effectiveness in reproducing anomalous diffusion that happens in transport dynamics in complex systems, such as fluid motion in viscoelastic medium, are the most important advantages [20], anomalous transfer in biology [21] and porous materials [22], etc. Control theory, entropy theory, image processing, and wave propagation phenomena all employing fractional calculus can be found in [2326]. The creation of tools to offer a mathematical structure for sophisticated physical systems and processes has been aided by breakthroughs in current variable order (VO) fractional calculus [27]. As a result of its appropriateness for modelling in a wide range of subjects, including science, engineering, and a variety of other disciplines, variable order fractional differential equations (VO-FDEs) have gained prominence [2831]. Physical modelling utilizing VO-FDE models has been the subject of a large-scale investigation. For example, Kobelev et al. [32] highlighted the dynamical and statistical systems with varying memory difficulties where the fractal dimension changes with coordinate and time. Coimbra et al. [33] used VO-fractional operators to investigate the viscoelasticity oscillator. Al-Mekhlafi and Sweilam [34] proposed a new multistrain TB model based on the VO-fractional derivative as a nonlinear ordinary differential equation extension. Due to the enormous number of applications, analytical and numerical techniques for solving variable order fractional order differential equations (VO-FDEs) have increased substantially in the last year. The analytical solution of VO-FDEs, on the other hand, is frequently difficult to obtain. Therefore, numerical approaches are used as sophisticated methods for numerical approximation of VO-FDEs in general [29, 3538].

The Caputo, the Liouville, the Marchaud, the Grunwald, and the Coimbra definitions are some of the recent variable order operator definitions suggested in the literature [33, 39]. Samko et al. [39] analyzed that the Riemann variable order definitions lost some features, meaning that the Marchaud operator is better than the Riemann-Liouville type operator. Ramirez et al. [40] also compared the variable order operators such as operators due to Riemann-Liouville, Marchaud, Caputo, and Coimbra using a simple criterion: the variable order operator must return the correct fractional derivative that corresponds to the argument of the functional order. Only the operator due to Coimbra and the Marchaud satisfy the aforementioned elementary condition [40], as well as the Coimbra variable order operator is more efficient numerically. Soon et al. [41] also demonstrated that the Coimbra variable order operator satisfies a mapping requirement and that it is the only formulation that returns the necessary derivatives as a function of for transitions between elastic and viscous regimes. Ramirez [40] demonstrated that the Coimbra concept is the most appropriate for physical modelling since it has essential properties that are desirable.

The numerical solution of the KdV problem of order and its eventual periodicity over confined domain is achieved using RBF with Coimbra variable order derivative. The following equations represent the proposed models in both one-dimensional and two-dimensional space:with the following initial conditionand the boundary conditions given bywhere .where , with the following boundary and initial conditions

The models in the above form are selected for the sake of comparison given in [42]. The Coimbra variable order derivative is defined by (7) in the next section.

1.1. The KdV Equation

Suppose that a wave propagates along a horizontal channel of the unform width along the positive direction of -axis alone. Let the depth of the channel be , be the time, and be the horizontal coordinate and let be the vertical distance of the fluid surface in equilibrium position. Let the amplitude of the wave be small enough, then the irrotational wave propagation can be modelled by the following equation known as the KdV equationwhere the first term represents the unform wave translation, and the other two terms and serve for the modification of the wave under the influences of nonlinear term and dispersive term , respectively.

1.2. Coimbra Variable Order Derivative

Modelling physical problems is better using the Coimbra concept. Variable order differentials are a useful tool for studying systems where the order changes with regard to one or more parameters, such as the management of a nonlinear viscoelasticity oscillator.

, , and the above operator require only one initial condition . The integer order derivative with respect to the variable is denoted by [33].

2. Analysis of RBF Approximation Method for Fractional Order KdV Equations

In the theory of multivariate approximation, the radial basis functions (RBF) method is the most extensively used tool. A generalized refinement of the multiquadric approach is RBF approximation methods. The MQ has a long history of application and theoretical research can be found in [43, 44]. The MQ approach is widely used in geology, geodesy, geophysics, and other domains, see [44]. Franke [45] conducted a comparative study in the field of MQ. Meanwhile, a key period in RBF history occurred, see for example [46], when Charles Micchelli refined the theory of the MQ method by establishing requirements that guarantee the system matrix nonsingularity for MQ methods. Schoenberg [47] is attribute with the results that generate the invertibility of the system matrix. Micchelli went on to say that Schoenberg’s constraints could be relaxed to allow many more functions to be included and that adequate conditions for functions could be applied to make the system matrix nonsingular. In 1990, physicist Kansa [48] discovery quickly disseminated the study, and RBF is used in a systematic approach for numerically solving partial differential equations and is meshless [49]. In numerous branches of applied areas [50, 51], a huge amount of mathematical applications of RBF are employed. In numerical techniques for solving PDEs with reasonable accuracy in multidimensions, Madych exposed the convergence rate of spectral order for MQ interpolation in [52]. In comparison to other state-of-the-art methodologies, these findings propelled RBF research forward swiftly, and the RBF methods drew appreciable attention in the literature as mesh-free approaches and their capacity to attain spectral accuracy for PDE numerical solutions on irregular domains [53]. In this work, a numerical scheme based on RBF and Coimbra derivative is constructed for fractional order KdV equations (1)–(5) defined in the following form:with the following boundary and initial conditionswith .

3. Variable Order Differential Operator Approximation

There are numerous definitions of varying order operators in the literature [54], but in the current study, we are using the definition due to Coimbra [33]. Because this variable order derivative has a great capability to model many complicated mechanical problems with accuracy, the Coimbra variable order operator has the capability to investigate and analyze the dynamics behaviour of many physical models, for example, the fractional forces which cannot be approximated accurately with constant order fractional operator or some other variable order derivatives. In the work [40], the authors performed a comparative study for solving a dynamical system and demonstrated that the Coimbra variable order derivative produced better results in many aspects than the nine definitions of variable order derivatives used in this study.

Now, for the numerical approximation of the Coimbra variable order derivative, we consider for , and let , where , then at time level , the Coimbra variable order derivative defined in (7) can be given by the following equation:

Let us denote the last term of this equation by , then we have

By using (11) in (10), we get the following form:where ; after further simplification, we getand by simplifying the integral involved, we have

Denoting the quantity by the , we getand assuming the value be denoted by , we get

Splitting the first term of this series and rewriting in the form, we have

Assuming , the approximation of Coimbra variable order derivative can be represented in the more simplified formwhich is the Coimbra variable order differential operator’s finite-difference approximation.

4. RBF Approximation Scheme

The RBF interpolant can be characterized as a linear combination of radial basis functions, as seen in the equation below. For a set of N scattered nodes ,where denotes the expansion coefficients at any time , denotes an RBF centred at , and denotes a distance norm in . It is possible to obtain the matrix form of (19) bywhere is a square matrix termed a system matrix, and the entries are . If is a spatial operator and is a boundary operator, then (19) is obtained.

The above two equations can be expressed in matrix form by

We obtain it in a more compressed formIn case of identity operators , the equation above can be written as

Model equations (1)–(5) can be approximated in the following way employing the weighted scheme, VODO finite-difference approximation, and RBF spatial operator approximation

Substituting the values of from (19) and (20), we get

This numerical strategy based on RBF can be solved at any point in time to acquire the value of , for which we can use (20). Using equation (20) to remove the value of , we get

The amplification matrix of the numerical scheme (27) is the following matrix :

In fact, the matrices and are identical, as shown by (20) and (24), because matrix is a special instance of matrix for identity operators. It is evident from the definitions of and that they are positive real integers, hence . As a result, (28) amplification matrix can be represented in a more basic form as follows:

Now, for , and denoting by , we obtain

Lemma 1. If is a square matrix of rank with negative eigenvalues, then the estimate for any isThis is also true for the Euclidean norm.

Proof:. Let the eigenvectors of the matrix be with the corresponding eigenvalues , thus we have the following equation:Suppose the vectors be orthonormal eigensystem for the matrix , since and all the eigenvalues of the matrix are negative, so we obtain

5. The Numerical Scheme’s Error Analysis

The VODO numerical scheme of order in time is , whereas RBF numerical scheme is mostly dependent on the RBF utilized for the derivation of other differentiation matrices and the RBF system matrix, as demonstrated in the previous work discussion. The order of convergence of several forms of RBF has been determined in [51]. Let the spatial numerical approximation corresponding to the present numerical scheme for a given RBF be of order , , and be the separation distance between the scattered nodes utilized for RBF interpolant. For the numerical scheme specified in (22), let be the approximate solution, be the precise solution, and be the error at time :

The above numerical technique’s amplification matrix, , is mostly determined by the type of RBF and the scale factor used. Assume that when the condition of Lemma 1 holds for a given optimal shape parameter value and optimal RBF option, thenis a criterion for the numerical scheme’s stability in (27). Assuming that both the initial solution value and the solution are sufficiently smooth along with ,where stands for a constant. At time , the error always fulfils the initial as well as the boundary condition via mathematical induction.when the condition (35) holds, then

This demonstrates that the current numerical method for VODO is convergent.

6. The Numerical Methodology for Variable Order Diffusion Models

Problem 1. We consider the KdV equation defined in (1)–(3) for the following value of function :with the following initial conditionwhere the boundary conditions can be extracted from the exact solutionwhere is the fractional order. The compactly supported radial basis function defined as with a support size of , is used to solve this issue across the spatial domain and the number of points are used. For various settings of order , step size , and collocation points , the results are presented in Figures 1 and 2 and Table 1, where the accuracy is quantified in terms of maximum error norm. When we gave a periodic boundary condition at , like with , the solution at each point of the domain is periodic in time for 1D fractional order KdV equation.

Problem 2. In this last example, we consider the following 2D KdV equation defined in (3)–(5) with the following value of the function with the following initial:and the boundary conditions can be taken from the following exact solution:, . The current RBF-based solution solves this problem, and the results are displayed in Figure 3 and Table 2, respectively. The current numerical technique appears to be convergent and stability attained when . This backs up the prior sections’ convergence and stability analysis of the current numerical system is achieved.

Problem 3. In the last example, we consider the irregular domain within the regular domain . The variable order which is used in this computation is defined by the function , . We used different number of nodes in the irregular domain. The radial basis function defined by is implemented in this problem and its corresponding shape parameter value changes solution accuracy, which is calculated using the formula [55]. The results are shown in graphical form and can be seen in Figure 4.

7. Conclusion

The numerical solution of the variable order KdV models in 1D and 2D is achieved using an RBF-based numerical approach. The RBF is used to approximate the spatial derivative, whereas the Coimbra derivative is used to approximate the variable order time differential operator. The numerical scheme’s stability and convergence are established. The current numerical technique is found to have a sensitivity in temporal integration. The periodicity of the KdV equation in 1D is explored, and it is demonstrated that the solution is periodic in time at each point of the domain for the fractional order KdV in 1D. The suggested numerical system provides the capacity to numerically approximate numerous complex mechanical problems with ease and precision. The Coimbra variable order operator can be used to examine and analyze the dynamics of a variety of physical models, such as fractional forces, which cannot be accurately modelled with a constant order fractional operator.

Data Availability

The date used to support this study are included in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.