Jungang Wang

Zongze Yang, Jungang Wang, Yan Li et al.

In this paper we investigate a sub-diffusion equation for simulating the anomalous diffusion phenomenon in real physical environment. Based on an equivalent transformation of the original sub-diffusion equation followed by the use of a smooth operator, we devise a high-order numerical scheme by combining the Nystrom method in temporal direction with the compact finite difference method and the spectral method in spatial direction. The distinct advantage of this approach in comparison with most current methods is its high convergence rate even though the solution of the anomalous sub-diffusion equation usually has lower regularity on the starting point. The effectiveness and efficiency of our proposed method are verified by several numerical experiments.

Jiahui Hu, Jungang Wang, Zhanbin Yuan et al.

In this paper, an alternating direction implicit (ADI) difference scheme for two-dimensional time-fractional wave equation of distributed-order with a nonlinear source term is presented. The unique solvability of the difference solution is discussed, and the unconditional stability and convergence order of the numerical scheme are analysed. Finally, numerical experiments are carried out to verify the effectiveness and accuracy of the algorithm.

Jiahui Hu, Jungang Wang, Zhanbin Yuan et al.

The fractional Feynman-Kac equations describe the distribution of functionals of non-Brownian motion, or anomalous diffusion, including two types called the forward and backward fractional Feynman-Kac equations, where the fractional substantial derivative is involved. This paper focuses on the more widely used backward version. Based on the discretized schemes for fractional substantial derivatives proposed recently, we construct compact finite difference schemes for the backward fractional Feynman-Kac equation, which has q-th (q=1, 2, 3, 4) order accuracy in temporal direction and fourth order accuracy in spatial direction, respectively. In the case q=1, the numerical stability and convergence of the difference scheme in the discrete L-infinity norm are proved strictly, where a new inner product is defined for the theoretical analysis. Finally, numerical examples are provided to verify the effectiveness and accuracy of the algorithms.