Yufeng Nie

Hao Dong, Junzhi Cui, Yufeng Nie et al.

This study develops a novel multiscale computational method for heat conduction problems of composite structures with diverse periodic configurations in different subdomains. Firstly, the second-order two-scale (SOTS) solutions for these multiscale problems are successfully obtained based on asymptotic homogenization method. Then, the error analysis in the pointwise sense is given to illustrate the importance of developing SOTS solutions. Furthermore, the error estimates for the SOTS approximate solutions in the integral sense is presented. In addition, a SOTS numerical algorithm is proposed to effectively solve these problems based on finite element method. Finally, some numerical examples verify the feasibility and effectiveness of the SOTS numerical algorithm we proposed.

Xiangrong Li, Nan Qi, Yufeng Nie et al.

The properties and applications of superconvergence on size-guaranteed Delaunay triangulation generated by bubble placement method (BPM), are studied in this paper. First, we derive a mesh condition that the difference between the actual side length and the desired length $h$ is as small as ${\cal O}(h^{1+α})$ $(α>0)$. Second, the superconvergence estimations are analyzed on linear and quadratic finite element for elliptic boundary value problem based on the above mesh condition. In particular, the mesh condition is suitable for many known superconvergence estimations of different equations. Numerical tests are provided to verify the theoretical findings and to exhibit the superconvergence property on BPM-based grids.

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.

Dan Li, Yufeng Nie, Chunmei Wang

A superconvergence error estimate for the gradient approximation of the second order elliptic problem in three dimensions is analyzed by using weak Galerkin finite element scheme on the uniform and non-uniform cubic partitions. Due to the loss of the symmetric property from two dimensions to three dimensions, this superconvergence result in three dimensions is not a trivial extension of the recent superconvergence result in two dimensions \cite{sup_LWW2018} from rectangular partitions to cubic partitions. The error estimate for the numerical gradient in the $L^{2}$-norm arrives at a superconvergence order of ${\cal O}(h^r) (1.5 \leq r\leq 2)$ when the lowest order weak Galerkin finite elements consisting of piecewise linear polynomials in the interior of the elements and piecewise constants on the faces of the elements are employed. A series of numerical experiments are illustrated to confirm the established superconvergence theory in three dimensions.

Xiaogang Zhu, Yufeng Nie

In this work, a new collocation approach using a combination of a wavelet operational matrix method and the exponential spline interpolation is proposed to solve the time-fractional convection-diffusion equation with variable coefficients. The operational matrix of fractional order integration is first derived based on sine-cosine wavelet functions, which helps to convert the underlying equation into a linear algebraic system. Then, an exponential B-spline method is introduced in spatial direction. On selecting a set of proper collocation points, the method in presence is evaluated on several test problems and the numerical results finally illustrate its validity and applicability.

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.