Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical solutions of electromagnetic wave model of fractional derivative using class of finite difference scheme
192
0
0.0
(
0
)
Added by
Vijay Kumar Patel
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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 scheme 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.
rate research
Read More
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 coefficients 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.
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.
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 equation 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.
In this work, we derive a nonstandard finite difference scheme for the SICA (Susceptible-Infected-Chronic-AIDS) model and analyze the dynamical properties of the discretized system. We prove that the discretized model is dynamically consistent with the continuous, maintaining the essential properties of the standard SICA model, namely, the positivity and boundedness of the solutions, equilibrium points, and their local and global stability.
The asymptotic stable region and long-time decay rate of solutions to linear homogeneous Caputo time fractional ordinary differential equations (F-ODEs) are known to be completely determined by the eigenvalues of the coefficient matrix. Very different from the exponential decay of solutions to classical ODEs, solutions of F-ODEs decay only polynomially, leading to the so-called Mittag-Leffler stability, which was already extended to semi-linear F-ODEs with small perturbations. This work is mainly devoted to the qualitative analysis of the long-time behavior of numerical solutions. By applying the singularity analysis of generating functions developed by Flajolet and Odlyzko (SIAM J. Disc. Math. 3 (1990), 216-240), we are able to prove that both $mathcal{L}$1 scheme and strong $A$-stable fractional linear multistep methods (F-LMMs) can preserve the numerical Mittag-Leffler stability for linear homogeneous F-ODEs exactly as in the continuous case. Through an improved estimate of the discrete fractional resolvent operator, we show that strong $A$-stable F-LMMs are also Mittag-Leffler stable for semi-linear F-ODEs under small perturbations. For the numerical schemes based on $alpha$-difference approximation to Caputo derivative, we establish the Mittag-Leffler stability for semi-linear problems by making use of properties of the Poisson transformation and the decay rate of the continuous fractional resolvent operator. Numerical experiments are presented for several typical time fractional evolutional equations, including time fractional sub-diffusion equations, fractional linear system and semi-linear F-ODEs. All the numerical results exhibit the typical long-time polynomial decay rate, which is fully consistent with our theoretical predictions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
