Two Finite Difference Schemes for Multi-Dimensional Fractional Wave Equations with Weakly Singular Solutions

Abstract A time-fractional initial-boundary value problem of wave type is considered, where the spatial domain is ( 0 , 1 ) d (0,1)^{d} for some d ∈ { 1 , 2 , 3 } d\in\{1,2,3\} . Regularity of the solution 𝑢 is discussed in detail. Typical solutions have a weak singularity at the initial time t = 0 t=0 : while 𝑢 and u t u_{t} are continuous at t = 0 t=0 , the second-order derivative u t ⁢ t u_{tt} blows up at t = 0 t=0 . To solve the problem numerically, a finite difference scheme is used on a mesh that is graded in time and uniform in space with the same mesh size ℎ in each coordinate direction. This scheme is generated through order reduction: one rewrites the differential equation as a system of two equations using the new variable v := u t v:=u_{t} ; then one uses a modified L1 scheme of Crank–Nicolson type for the driving equation. A fast variant of this finite difference scheme is also considered, using a sum-of-exponentials (SOE) approximation for the kernel function in the Caputo derivative. The stability and convergence of both difference schemes are analysed in detail. At each time level, the system of linear equations generated by the difference schemes is solved by a fast Poisson solver, thereby taking advantage of the fast difference scheme. Finally, numerical examples are presented to demonstrate the accuracy and efficiency of both numerical methods.

[1]  Jinye Shen.  Fast Finite Difference Schemes for Time-Fractional Diffusion Equations with a Weak Singularity at Initial Time , 2018, East Asian Journal on Applied Mathematics.

[2]  Ernst P. Stephan,et al.  Decompositions in Edge and Corner Singularities for the Solution of the Dirichlet Problem of the Laplacian in a Polyhedron , 1990 .

[3]  Zhi-Zhong Sun,et al.  Two Alternating Direction Implicit Difference Schemes for Solving the Two-Dimensional Time Distributed-Order Wave Equations , 2016, J. Sci. Comput..

[4]  William McLean,et al.  Superconvergence of a Discontinuous Galerkin Method for Fractional Diffusion and Wave Equations , 2012, SIAM J. Numer. Anal..

[5]  Martin Stynes,et al.  Too much regularity may force too much uniqueness , 2016, 1607.01955.

[6]  Zhi-Zhong Sun,et al.  An H2N2 Interpolation for Caputo Derivative with Order in (1, 2) and Its Application to Time-Fractional Wave Equations in More Than One Space Dimension , 2020, J. Sci. Comput..

[7]  Gene H. Golub,et al.  On direct methods for solving Poisson's equation , 1970, Milestones in Matrix Computation.

[8]  F. Mainardi.  Fractional Relaxation-Oscillation and Fractional Diffusion-Wave Phenomena , 1996 .

[9]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[10]  Changpin Li,et al.  A new Crank–Nicolson finite element method for the time-fractional subdiffusion equation☆ , 2017 .

[11]  Bangti Jin,et al.  Numerical methods for time-fractional evolution equations with nonsmooth data: A concise overview , 2018, Computer Methods in Applied Mechanics and Engineering.

[12]  Hong Sun,et al.  The Temporal Second Order Difference Schemes Based on the Interpolation Approximation for the Time Multi-term Fractional Wave Equation , 2019, J. Sci. Comput..

[13]  Masahiro Yamamoto,et al.  Initial value/boundary value problems for fractional diffusion-wave equations and applications to some inverse problems , 2011 .

[14]  Zhongqiang Zhang,et al.  Fast difference schemes for solving high-dimensional time-fractional subdiffusion equations , 2016, J. Comput. Phys..

[15]  Pin Lyu,et al.  A fast linearized finite difference method for the nonlinear multi-term time-fractional wave equation , 2019, Applied Numerical Mathematics.

[16]  Natalia Kopteva,et al.  Error analysis of the L1 method on graded and uniform meshes for a fractional-derivative problem in two and three dimensions , 2017, Math. Comput..

[17]  Zhi‐zhong Sun,et al.  A fully discrete difference scheme for a diffusion-wave system , 2006 .