ترغب بنشر مسار تعليمي؟ اضغط هنا

Arbitrary order of convergence for Riesz fractional derivative via central difference method

80   0   0.0 ( 0 )
 نشر من قبل Pui Ho Lam
 تاريخ النشر 2021
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English




اسأل ChatGPT حول البحث

We propose a novel method to compute a finite difference stencil for Riesz derivative for artibitrary speed of convergence. This method is based on applying a pre-filter to the Grunwald-Letnikov type central difference stencil. The filter is obtained by solving for the inverse of a symmetric Vandemonde matrix and exploiting the relationship between the Taylors series coefficients and fast Fourier transform. The filter costs Oleft(N^{2}right) operations to evaluate for Oleft(h^{N}right) of convergence, where h is the sampling distance. The higher convergence speed should more than offset the overhead with the requirement of the number of nodal points for a desired error tolerance significantly reduced. The benefit of progressive generation of the stencil coefficients for adaptive grid size for dynamic problems with the Grunwald-Letnikov type difference scheme is also kept because of the application of filtering. The higher convergence rate is verified through numerical experiments.



قيم البحث

اقرأ أيضاً

201 - Jing Sun , Weihua Deng , Daxin Nie 2021
We make the split of the integral fractional Laplacian as $(-Delta)^s u=(-Delta)(-Delta)^{s-1}u$, where $sin(0,frac{1}{2})cup(frac{1}{2},1)$. Based on this splitting, we respectively discretize the one- and two-dimensional integral fractional Laplaci an with the inhomogeneous Dirichlet boundary condition and give the corresponding truncation errors with the help of the interpolation estimate. Moreover, the suitable corrections are proposed to guarantee the convergence in solving the inhomogeneous fractional Dirichlet problem and an $mathcal{O}(h^{1+alpha-2s})$ convergence rate is obtained when the solution $uin C^{1,alpha}(bar{Omega}^{delta}_{n})$, where $n$ is the dimension of the space, $alphain(max(0,2s-1),1]$, $delta$ is a fixed positive constant, and $h$ denotes mesh size. Finally, the performed numerical experiments confirm the theoretical results.
In this work, new finite difference schemes are presented for dealing with the upper-convected time derivative in the context of the generalized Lie derivative. The upper-convected time derivative, which is usually encountered in the constitutive equ ation of the popular viscoelastic models, is reformulated in order to obtain approximations of second-order in time for solving a simplified constitutive equation in one and two dimensions. The theoretical analysis of the truncation errors of the methods takes into account the linear and quadratic interpolation operators based on a Lagrangian framework. Numerical experiments illustrating the theoretical results for the model equation defined in one and two dimensions are included. Finally, the finite difference approximations of second-order in time are also applied for solving a two-dimensional Oldroyd-B constitutive equation subjected to a prescribed velocity field at different Weissenberg numbers.
393 - Vijay Kumar Patel 2021
The present article is devoting a numerical approach for solving a fractional partial differential equation (FPDE) arising from electromagnetic waves in dielectric media (EMWDM). The truncated Bernoulli and Hermite wavelets series with unknown coeffi cients have been used to approximate the solution in both the temporal and spatial variables. The basic idea for discretizing the FPDE is wavelet approximation based on the Bernoulli and Hermite wavelets operational matrices of integration and differentiation. The resulted system of a linear algebraic equation has been solved by the collocation method. Moreover, convergence and error analysis have been discussed. Finally, several numerical experiments with different fractional-order derivatives have been provided and compared with the exact analytical solutions to illustrate the accuracy and efficiency of the method.
The aim of this paper is to develop fast second-order accurate difference schemes for solving one- and two-dimensional time distributed-order and Riesz space fractional diffusion equations. We adopt the same measures for one- and two-dimensional prob lems as follows: we first transform the time distributed-order fractional diffusion problem into the multi-term time-space fractional diffusion problem with the composite trapezoid formula. Then, we propose a second-order accurate difference scheme based on the interpolation approximation on a special point to solve the resultant problem. Meanwhile, the unconditional stability and convergence of the new difference scheme in $L_2$-norm are proved. Furthermore, we find that the discretizations lead to a series of Toeplitz systems which can be efficiently solved by Krylov subspace methods with suitable circulant preconditioners. Finally, numerical results are presented to show the effectiveness of the proposed difference methods and demonstrate the fast convergence of our preconditioned Krylov subspace methods.
In this article, a numerical scheme is introduced for solving the fractional partial differential equation (FPDE) arising from electromagnetic waves in dielectric media (EMWDM) by using an efficient class of finite difference methods. The numerical s cheme is based on the Hermite formula. The Caputos fractional derivatives in time are discretized by a finite difference scheme of order $mathcal{O}(k^{(4-alpha)})$ & $mathcal{O}(k^{(4-beta)})$, $1<beta <alpha leq 2$. The stability and the convergence analysis of the proposed methods are given by a procedure similar to the standard von Neumann stability analysis under mild conditions. Also for FPDE, accuracy of order $mathcal{O}left( k^{(4-alpha)}+k^{(4-beta)}+h^2right) $ is investigated. Finally, several numerical experiments with different fractional-order derivatives are provided and compared with the exact solutions to illustrate the accuracy and efficiency of the scheme. A comparative numerical study is also done to demonstrate the efficiency of the proposed scheme.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا