Mathematical Modeling of Linear Fractional Oscillators

Application of the Explicit Euler Method for Numerical Analysis of a Nonlinear Fractional Oscillation Equation

In this paper, a numerical analysis of the oscillation equation with a derivative of a fractional variable Riemann–Liouville order in the dissipative term, which is responsible for viscous friction, is carried out. Using the theory of finite-difference schemes, an explicit finite-difference scheme (Euler’s method) was constructed on a uniform computational grid. For the first time, the issues of approximation, stability and convergence of the proposed explicit finite-difference scheme are considered. To compare the results, the Adams–Bashford–Moulton scheme was constructed as an experimental method. The theoretical results were confirmed using test examples, the computational accuracy of the method was evaluated, which is consistent with the theoretical one, and the simulation results were visualized. Using the example of a fractional Duffing oscillator, waveforms and phase trajectories, as well as its amplitude–frequency characteristics, were constructed using a finite-difference scheme. To identify chaotic regimes, the spectra of maximum Lyapunov exponents and Poincaré points were constructed. It is shown that an explicit finite-difference scheme can be acceptable under the condition of a step of the computational grid.

Application of the Fractional Riccati Equation for Mathematical Modeling of Dynamic Processes with Saturation and Memory Effect

In this study, the model Riccati equation with variable coefficients as functions, as well as a derivative of a fractional variable order (VO) of the Gerasimov-Caputo type, is used to approximate the data for some physical processes with saturation. In particular, the proposed model is applied to the description of solar activity (SA), namely the number of sunspots observed over the past 25 years. It is also used to describe data from Johns Hopkins University on coronavirus infection COVID-19, in particular data on the Russian Federation and the Republic of Uzbekistan. Finally, it is used to study issues related to seismic activity, in particular, the description of data on the volumetric activity of Radon (RVA). The Riccati equation used in the mathematical model was numerically solved by constructing an implicit finite difference scheme (IFDS) and its implementation by the modified Newton method (MNM). The calculated curves obtained in the study are compared with known experimental data. It is shown that if the model parameters are chosen appropriately, the model curves will give results that correlate well with real experimental data. Moreover, with other parameters of the model, it is possible to make some prediction about the possible course of the considered processes.

Investigation of Finite-Difference Schemes for the Numerical Solution of a Fractional Nonlinear Equation

The article discusses different schemes for the numerical solution of the fractional Riccati equation with variable coefficients and variable memory, where the fractional derivative is understood in the sense of Gerasimov-Caputo. For a nonlinear fractional equation, in the general case, theorems of approximation, stability, and convergence of a nonlocal implicit finite difference scheme (IFDS) are proved. For IFDS, it is shown that the scheme converges with the order corresponding to the estimate for approximating the Gerasimov-Caputo fractional operator. The IFDS scheme is solved by the modified Newton’s method (MNM), for which it is shown that the method is locally stable and converges with the first order of accuracy. In the case of the fractional Riccati equation, approximation, stability, and convergence theorems are proved for a nonlocal explicit finite difference scheme (EFDS). It is shown that EFDS conditionally converges with the first order of accuracy. On specific test examples, the computational accuracy of numerical methods was estimated according to Runge’s rule and compared with the exact solution. It is shown that the order of computational accuracy of numerical methods tends to the theoretical order of accuracy with increasing nodes of the computational grid.

Some Aspects of Numerical Analysis for a Model Nonlinear Fractional Variable Order Equation

The article proposes a nonlocal explicit finite-difference scheme for the numerical solution of a nonlinear, ordinary differential equation with a derivative of a fractional variable order of the Gerasimov–Caputo type. The questions of approximation, convergence, and stability of this scheme are studied. It is shown that the nonlocal finite-difference scheme is conditionally stable and converges to the first order. Using the fractional Riccati equation as an example, the computational accuracy of the numerical method is analyzed. It is shown that with an increase in the nodes of the computational grid, the order of computational accuracy tends to unity, i.e., to the theoretical value of the order of accuracy.

Exact Solutions of Bernoulli and Logistic Fractional Differential Equations with Power Law Coefficients

In this paper, we consider a nonlinear fractional differential equation. This equation takes the form of the Bernoulli differential equation, where we use the Caputo fractional derivative of non-integer order instead of the first-order derivative. The paper proposes an exact solution for this equation, in which coefficients are power law functions. We also give conditions for the existence of the exact solution for this non-linear fractional differential equation. The exact solution of the fractional logistic differential equation with power law coefficients is also proposed as a special case of the proposed solution for the Bernoulli fractional differential equation. Some applications of the Bernoulli fractional differential equation to describe dynamic processes with power law memory in physics and economics are suggested.

Mathematical Model of Fractional Duffing Oscillator with Variable Memory

The article investigates a mathematical model of the Duffing oscillator with a variable fractional order derivative of the Riemann–Liouville type. The study of the model is carried out using a numerical scheme based on the approximation of the fractional derivative of the Riemann–Liouville type by a discrete analog—the fractional derivative of Grunwald–Letnikov. The adequacy of the numerical scheme is verified using specific examples. Using a numerical algorithm, oscillograms and phase trajectories are constructed depending on the values of the model parameters. Chaotic regimes of the Duffing fractional oscillator are investigated using the Wolf–Bennetin algorithm. The forced oscillations of the Duffing fractional oscillator are investigated using the harmonic balance method. Analytical formulas for the amplitude-frequency, phase-frequency characteristics, and also the quality factor are obtained. It is shown that the fractional Duffing oscillator possesses different modes: regular, chaotic, multi-periodic. The relationship between the order of the fractional derivative and the quality factor of the oscillatory system is established.