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

2021 ◽  
Vol 6 (1) ◽  
pp. 23
Author(s):  
Dmitriy Tverdyi ◽  
Roman Parovik

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.

2021 ◽  
Vol 26 (3) ◽  
pp. 55
Author(s):  
Roman Parovik ◽  
Dmitriy Tverdyi

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.


2014 ◽  
Vol 2014 ◽  
pp. 1-22 ◽  
Author(s):  
Zongqi Liang ◽  
Yubin Yan ◽  
Guorong Cai

A standard Crank-Nicolson finite-difference scheme and a Dufort-Frankel finite-difference scheme are introduced to solve two-dimensional damped and undamped sine-Gordon equations. The stability and convergence of the numerical methods are considered. To avoid solving the nonlinear system, the predictor-corrector techniques are applied in the numerical methods. Numerical examples are given to show that the numerical results are consistent with the theoretical results.


Author(s):  
Fouad Mohammad Salama ◽  
Norhashidah Hj. Mohd Ali

In this paper, a hybrid method based on the Laplace transform and implicit finite difference scheme is applied to obtain the numerical solution of the two-dimensional time fractional advection-diffusion equation (2D-TFADE). Some of the major limitations in computing the numerical solution for fractional differential equations (FDEs) in multi-dimensional space are the huge computational cost and storage requirement, which are O(N^2) cost and O(MN) storage, N and M are the total number of time levels and space grid points, respectively. The proposed method reduced the computational complexity efficiently as it requires only O(N) computational cost and O(M) storage with reasonable accuracy when numerically solving the TFADE. The method’s stability and convergence are also investigated. The Results of numerical experiments of the proposed method are obtained and found to compare well with the results of existing standard finite difference scheme.


2001 ◽  
Vol 09 (01) ◽  
pp. 183-203 ◽  
Author(s):  
DMITRY MIKHIN

The energy conservation law and the flow reversal theorem are valid for underwater acoustic fields. In media at rest the theorem transforms into well-known reciprocity principle. The presented parabolic equation (PE) model strictly preserves these important physical properties in the numerical solution. The new PE is obtained from the one-way wave equation by Godin12 via Padé approximation of the square root operator and generalized to the case of moving media. The PE is range-dependent and explicitly includes range derivatives of the medium parameters. Implicit finite difference scheme solves the PE written in terms of energy flux. Such formalism inherently provides simple and exact energy-conserving boundary condition at vertical interfaces. The finite-difference operators, the discreet boundary conditions, and the self-starter are derived by discretization of the differential PE. Discreet energy conservation and flow reversal theorem are rigorously proved as mathematical properties of the finite-difference scheme and confirmed by numerical modeling. Numerical solution is shown to be reciprocal with accuracy of 10–12 decimal digits, which is the accuracy of round-off errors. Energy conservation and wide-angle capabilities of the model are illustrated by comparison with two-way normal mode solutions including the ASA benchmark wedge.


Sign in / Sign up

Export Citation Format

Share Document