
In this paper, we propose a timefractional molecular beam epitaxy (MBE)
model with slope selection and its efficient, accurate, full discrete, linear
numerical approximation. The numerical scheme utilizes the fast algorithm for
the Caputo fractional derivative operator in time discretization and Fourier
spectral method in spatial discretization. Refinement tests are conducted to
verify the $2\alpha$ order of time convergence, with $\alpha \in (0, 1]$ the
fractional order of derivative. Several numerical simulations are presented to
demonstrate the accuracy and efficiency of our newly proposed scheme. By
exploring the fast algorithm calculating the Caputo fractional derivative, our
numerical scheme makes it practice for long time simulation of MBE coarsening,
which is essential for MBE model in practice. With the proposed fractional MBE
model, we observe that the scaling law for the energy decays as $
O(t^{\frac{\alpha}{3}})$ and the roughness increases as
$O(t^{\frac{\alpha}{3}})$, during the coarsening dynamics with random initial
condition. That is to say, the coarsening rate of MBE model could be
manipulated by the fractional order $\alpha$, and it is linearly proportional
to $\alpha$. This is the first time in literature to report/discover such
scaling correlation. It provides a potential application field for fractional
differential equations. Besides, the numerical approximation strategy proposed
in this paper can be readily applied to study many classes of timefractional
and high dimensional phase field models.

In this paper, we introduce a class of high order immersed finite volume
methods (IFVM) for onedimensional interface problems. We show the optimal
convergence of IFVM in H1 and L2 norms. We also prove some superconvergence
results of IFVM. To be more precise, the IFVM solution is superconvergent of
order p+2 at the roots of generalized Lobatto polynomials, and the flux is
superconvergent of order p+1 at generalized Gauss points on each element
including the interface element. Furthermore, for diffusion interface problems,
the convergence rates for IFVM solution at the mesh points and the flux at
generalized Gaussian points can both be raised to 2p. These superconvergence
results are consistent with those for the standard finite volume methods.
Numerical examples are provided to confirm our theoretical analysis.

In this article, we study superconvergence properties of immersed finite
element methods for the one dimensional elliptic interface problem. Due to low
global regularity of the solution, classical superconvergence phenomenon for
finite element methods disappears unless the discontinuity of the coefficient
is resolved by partition. We show that immersed finite element solutions
inherit all desired superconvergence properties from standard finite element
methods without requiring the mesh to be aligned with the interface. In
particular, on interface elements, superconvergence occurs at roots of
generalized orthogonal polynomials that satisfy both orthogonality and
interface jump conditions.

In this paper, we study superconvergence properties of the local
discontinuous Galerkin method for onedimensional linear parabolic equations
when alternating fluxes are used. We prove, for any polynomial degree $k$, that
the numerical fluxes converge at a rate of $2k+1$ (or $2k+1/2$) for all mesh
nodes and the domain average under some suitable initial discretization.
We further prove a $k+1$th superconvergence rate for the derivative
approximation and a $k+2$th superconvergence rate for the function value
approximation at the Radau points.
Numerical experiments demonstrate that in most cases, our error estimates are
optimal, i.e., the error bounds are sharp.

This paper is concerned with superconvergence properties of a class of finite
volume methods of arbitrary order over rectangular meshes. Our main result is
to prove {\it 2kconjecture}: at each vertex of the underlying rectangular
mesh, the bi$k$ degree finite volume solution approximates the exact solution
with an order
$ O(h^{2k})$, where $h$ is the mesh size. As byproducts, superconvergence
properties for finite volume discretization errors at Lobatto and Gauss points
are also obtained. All theoretical findings are confirmed by numerical
experiments.

We present and analyze a finite volume scheme of arbitrary order for elliptic
equations in the onedimensional setting. In this scheme, the control volumes
are constructed by using the Gauss points in subintervals of the underlying
mesh. We provide a unified proof for the infsup condition, and show that our
finite volume scheme has optimal convergence rate under the energy and $L^2$
norms of the approximate error. Furthermore, we prove that the derivative error
is superconvergent at all Gauss points and in some special case, the
convergence rate can reach $h^{2r}$, where $r$ is the polynomial degree of the
trial space. All theoretical results are justified by numerical tests.