2,766 research outputs found
Finite difference/local discontinuous Galerkin method for solving the fractional diffusion-wave equation
In this paper a finite difference/local discontinuous Galerkin method for the
fractional diffusion-wave equation is presented and analyzed. We first propose
a new finite difference method to approximate the time fractional derivatives,
and give a semidiscrete scheme in time with the truncation error , where is the time step size. Further we develop a fully
discrete scheme for the fractional diffusion-wave equation, and prove that the
method is unconditionally stable and convergent with order , where is the degree of piecewise polynomial. Extensive numerical
examples are carried out to confirm the theoretical convergence rates.Comment: 18 pages, 2 figure
Optimal Collocation Nodes for Fractional Derivative Operators
Spectral discretizations of fractional derivative operators are examined,
where the approximation basis is related to the set of Jacobi polynomials. The
pseudo-spectral method is implemented by assuming that the grid, used to
represent the function to be differentiated, may not be coincident with the
collocation grid. The new option opens the way to the analysis of alternative
techniques and the search of optimal distributions of collocation nodes, based
on the operator to be approximated. Once the initial representation grid has
been chosen, indications on how to recover the collocation grid are provided,
with the aim of enlarging the dimension of the approximation space. As a
results of this process, performances are improved. Applications to fractional
type advection-diffusion equations, and comparisons in terms of accuracy and
efficiency are made. As shown in the analysis, special choices of the nodes can
also suggest tricks to speed up computations
Domain Decomposition Methods for Space Fractional Partial Differential Equations
In this paper, a two-level additive Schwarz preconditioner is proposed for
solving the algebraic systems resulting from the finite element approximations
of space fractional partial differential equations (SFPDEs). It is shown that
the condition number of the preconditioned system is bounded by C(1+H/\delta),
where H is the maximum diameter of subdomains and \delta is the overlap size
among the subdomains. Numerical results are given to support our theoretical
findings.Comment: 19 pages, three figure
Strang-type preconditioners for solving fractional diffusion equations by boundary value methods
The finite difference scheme with the shifted Gr\"{u}nwarld formula is
employed to semi-discrete the fractional diffusion equations. This spatial
discretization can reduce to the large system of ordinary differential
equations (ODEs) with initial values. Recently, boundary value method (BVM) was
developed as a popular algorithm for solving large systems of ODEs. This method
requires the solutions of one or more nonsymmetric, large and sparse linear
systems. In this paper, the GMRES method with the block circulant
preconditioner is proposed for solving these linear systems. One of the main
results is that if an -stable boundary value method is used
for an m-by-m system of ODEs, then the preconditioner is invertible and the
preconditioned matrix can be decomposed as I+L, where I is the identity matrix
and the rank of L is at most . It means that when the GMRES
method is applied to solve the preconditioned linear systems, the method will
converge in at most iterations.Finally, extensive numerical
experiments are reported to illustrate the effectiveness of our methods for
solving the fractional diffusion equations.Comment: 19 pages,4 figure
An advanced meshless approach for the high-dimensional multi-term time-space-fractional PDEs on convex domains
In this article, an advanced differential quadrature (DQ) approach is
proposed for the high-dimensional multi-term time-space-fractional partial
differential equations (TSFPDEs) on convex domains. Firstly, a family of
high-order difference schemes is introduced to discretize the time-fractional
derivative and a semi-discrete scheme for the considered problems is presented.
We strictly prove its unconditional stability and error estimate. Further, we
derive a class of DQ formulas to evaluate the fractional derivatives, which
employs radial basis functions (RBFs) as test functions. Using these DQ
formulas in spatial discretization, a fully discrete DQ scheme is then
proposed. Our approach provides a flexible and high accurate alternative to
solve the high-dimensional multi-term TSFPDEs on convex domains and its actual
performance is illustrated by contrast to the other methods available in the
open literature. The numerical results confirm the theoretical analysis and the
capability of our proposed method finally.Comment: 22 pages, 26 figure
Multiresolution Galerkin method for solving the functional distribution of anomalous diffusion described by time-space fractional diffusion equation
The functional distributions of particle trajectories have wide applications,
including the occupation time in half-space, the first passage time, and the
maximal displacement, etc. The models discussed in this paper are for
characterizing the distribution of the functionals of the paths of anomalous
diffusion described by time-space fractional diffusion equation. This paper
focuses on providing effective computation methods for the models. Two kinds of
time stepping schemes are proposed for the fractional substantial derivative.
The multiresolution Galerkin method with wavelet B-spline is used for space
approximation. Compared with the finite element or spectral polynomial bases,
the wavelet B-spline bases have the advantage of keeping the Toeplitz structure
of the stiffness matrix, and being easy to generate the matrix elements and to
perform preconditioning. The unconditional stability and convergence of the
provided schemes are theoretically proved and numerically verified. Finally, we
also discuss the efficient implementations and some extensions of the schemes,
such as the wavelet preconditioning and the non-uniform time discretization.Comment: 31 pages, 1 figur
Discontinuous Galerkin methods and their adaptivity for the tempered fractional (convection) diffusion equations
This paper focuses on the adaptive discontinuous Galerkin (DG) methods for
the tempered fractional (convection) diffusion equations. The DG schemes with
interior penalty for the diffusion term and numerical flux for the convection
term are used to solve the equations, and the detailed stability and
convergence analyses are provided. Based on the derived posteriori error
estimates, the local error indicator is designed. The theoretical results and
the effectiveness of the adaptive DG methods are respectively verified and
displayed by the extensive numerical experiments. The strategy of designing
adaptive schemes presented in this paper works for the general PDEs with
fractional operators.Comment: 31 pages, 5 figure
Parallel-in-Time with Fully Finite Element Multigrid for 2-D Space-fractional Diffusion Equations
The paper investigates a non-intrusive parallel time integration with
multigrid for space-fractional diffusion equations in two spatial dimensions.
We firstly obtain a fully discrete scheme via using the linear finite element
method to discretize spatial and temporal derivatives to propagate solutions.
Next, we present a non-intrusive time-parallelization and its two-level
convergence analysis, where we algorithmically and theoretically generalize the
MGRIT to time-dependent fine time-grid propagators. Finally, numerical
illustrations show that the obtained numerical scheme possesses the saturation
error order, theoretical results of the two-level variant deliver good
predictions, and significant speedups can be achieved when compared to parareal
and the sequential time-stepping approach.Comment: 20 pages, 4 figures, 8 table
A spectral penalty method for two-sided fractional differential equations with general boundary conditions
We consider spectral approximations to the conservative form of the two-sided
Riemann-Liouville (R-L) and Caputo fractional differential equations (FDEs)
with nonhomogeneous Dirichlet (fractional and classical, respectively) and
Neumann (fractional) boundary conditions. In particular, we develop a spectral
penalty method (SPM) by using the Jacobi poly-fractonomial approximation for
the conservative R-L FDEs while using the polynomial approximation for the
conservative Caputo FDEs. We establish the well-posedness of the corresponding
weak problems and analyze sufficient conditions for the coercivity of the SPM
for different types of fractional boundary value problems. This analysis allows
us to estimate the proper values of the penalty parameters at boundary points.
We present several numerical examples to verify the theory and demonstrate the
high accuracy of SPM, both for stationary and time dependent FDEs. Moreover, we
compare the results against a Petrov-Galerkin spectral tau method (PGS-,
an extension of [Z. Mao, G.E. Karniadakis, SIAM J. Numer. Anal., 2018]) and
demonstrate the superior accuracy of SPM for all cases considered.Comment: 27 page
Operator-Based Uncertainty Quantification of Stochastic Fractional PDEs
Fractional calculus provides a rigorous mathematical framework to describe
anomalous stochastic processes by generalizing the notion of classical
differential equations to their fractional-order counterparts. By introducing
the fractional orders as uncertain variables, we develop an operator-based
uncertainty quantification framework in the context of stochastic fractional
partial differential equations (SFPDEs), subject to additive random noise. We
characterize different sources of uncertainty and then, propagate their
associated randomness to the system response by employing a probabilistic
collocation method (PCM). We develop a fast, stable, and convergent
Petrov-Galerkin spectral method in the physical domain in order to formulate
the forward solver in simulating each realization of random variables in the
sampling procedure
- …