Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFast Algorithms and Error Analysis of Caputo Derivatives with Small Factional Orders
596
0
0.0
(
0
)
Added by
Zihang Zhang
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this paper, we investigate fast algorithms to approximate the Caputo derivative $^C_0D_t^alpha u(t)$ when $alpha$ is small. We focus on two fast algorithms, i.e. FIR and FIDR, both relying on the sum-of-exponential approximation to reduce the cost of evaluating the history part. FIR is the numerical scheme originally proposed in [16], and FIDR is an alternative scheme we propose in this work, and the latter shows superiority when $alpha$ is small. With quantitative estimates, we prove that given a certain error threshold, the computational cost of evaluating the history part of the Caputo derivative can be decreased as $alpha$ gets small. Hence, only minimal cost for the fast evaluation is required in the small $alpha$ regime, which matches prevailing protocols in engineering practice. We also present a stability and error analysis of FIDR for solving linear fractional diffusion equations. Finally, we carry out systematic numerical studies for the performances of both FIR and FIDR schemes, where we explore the trade-off between accuracy and efficiency when $alpha$ is small.
rate research
Read More
We consider a discretization of Caputo derivatives resulted from deconvolving a scheme for the corresponding Volterra integral. Properties of this discretization, including signs of the coefficients, comparison principles, and stability of the corresponding implicit schemes, are proved by its linkage to Volterra integrals with completely monotone kernels. We then apply the backward scheme corresponding to this discretization to two time fractional dissipative problems, and these implicit schemes are helpful for the analysis of the corresponding problems. In particular, we show that the overdamped generalized Langevin equation with fractional noise has a unique limiting measure for strongly convex potentials and establish the convergence of numerical solutions to the strong solutions of time fractional gradient flows. The proposed scheme and schemes derived using the same philosophy can be useful for many other applications as well.
We present a new fractional Taylor formula for singular functions whose Caputo fractional derivatives are of bounded variation. It bridges and ``interpolates the usual Taylor formulas with two consecutive integer orders. This enables us to obtain an analogous formula for the Legendre expansion coefficient of this type of singular functions, and further derive the optimal (weighted) $L^infty$-estimates and $L^2$-estimates of the Legendre polynomial approximations. This set of results can enrich the existing theory for $p$ and $hp$ methods for singular problems, and answer some open questions posed in some recent literature.
There are three main types of numerical computations for the Bessel function of the second kind: series expansion, continued fraction, and asymptotic expansion. In addition, they are combined in the appropriate domain for each. However, there are some regions where the combination of these types requires sufficient computation time to achieve sufficient accuracy, however, efficiency is significantly reduced when parallelized. In the proposed method, we adopt a simple numerical integration concept of integral representation. We coarsely refine the integration range beforehand, and stabilize the computation time by performing the integration calculation at a fixed number of intervals. Experiments demonstrate that the proposed method can achieve the same level of accuracy as existing methods in less than half the computation time.
This paper describes the analysis of Lagrange interpolation errors on tetrahedrons. In many textbooks, the error analysis of Lagrange interpolation is conducted under geometric assumptions such as shape regularity or the (generalized) maximum angle condition. In this paper, we present a new estimation in which the error is bounded in terms of the diameter and projected circumradius of the tetrahedron. Because we do not impose any geometric restrictions on the tetrahedron itself, our error estimation may be applied to any tetrahedralizations of domains including very thin tetrahedrons.
The purpose of this paper is to perform an error analysis of the variational integrators of mechanical systems subject to external forcing. Essentially, we prove that when a discretization of contact order $r$ of the Lagrangian and force are used, the integrator has the same contact order. Our analysis is performed first for discrete forced mechanical systems defined over $TQ$, where we study the existence of flows, the construction and properties of discrete exact systems and the contact order of the flows (variational integrators) in terms of the contact order of the original systems. Then we use those results to derive the corresponding analysis for the analogous forced systems defined over $Qtimes Q$.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
