Some efficient numerical schemes are proposed for solving one-dimensional (1D) and two-dimensional (2D) multi-term time fractional sub-diffusion equations, combining the compact difference approach for the spatial discretisation and L1 approximation for the multi-term time Caputo fractional derivatives. The stability and convergence of these difference schemes are theoretically established. Several numerical examples are implemented, testifying to their efficiency and confirming their convergence order.

Mean curvature-based image registration model firstly proposed by Chumchob-Chen-Brito (2011) offered a better regularizer technique for both smooth and nonsmooth deformation fields. However, it is extremely challenging to solve efficiently this model and the existing methods are slow or become efficient only with strong assumptions on the smoothing parameter β. In this paper, we take a different solution approach. Firstly, we discretize the joint energy functional, following an idea of relaxed fixed point is implemented and combine with Gauss-Newton scheme with Armijo's Linear Search for solving the discretized mean curvature model and further to combine with a multilevel method to achieve fast convergence. Numerical experiments not only confirm that our proposed method is efficient and stable, but also it can give more satisfying registration results according to image quality.

The numerical solution of the time-fractional sub-diffusion equation on an unbounded domain in two-dimensional space is considered, where a circular artificial boundary is introduced to divide the unbounded domain into a bounded computational domain and an unbounded exterior domain. The local artificial boundary conditions for the fractional sub-diffusion equation are designed on the circular artificial boundary by a joint Laplace transform and Fourier series expansion, and some auxiliary variables are introduced to circumvent high-order derivatives in the artificial boundary conditions. The original problem defined on the unbounded domain is thus reduced to an initial boundary value problem on a bounded computational domain. A finite difference and L1 approximation are applied for the space variables and the Caputo time-fractional derivative, respectively. Two numerical examples demonstrate the performance of the proposed method.

A class of finite volume methods is developed for pricing either European or American options under jump-diffusion models based on a linear finite element space. An easy to implement linear interpolation technique is derived to evaluate the integral term involved, and numerical analyses show that the full discrete system matrices are M-matrices. For European option pricing, the resulting dense linear systems are solved by the generalised minimal residual (GMRES) method; while for American options the resulting linear complementarity problems (LCP) are solved using the modulus-based successive overrelaxation (MSOR) method, where the H+-matrix property of the system matrix guarantees convergence. Numerical results are presented to demonstrate the accuracy, efficiency and robustness of these methods.

We review recent theoretical progress in the physical understanding of far-from-equilibrium phenomena seen experimentally in epitaxial growth and erosion on crystal surfaces. The formation and dynamics of various interface structures (pyramids, ripples, etc.), and also kinetic phase transitions observed between these structures, can all be understood within a simple continuum model based on the mass conservation law and respecting the symmetries of the growing crystal surface. In particular, theoretical predictions and experimental results are compared for (001), (110) and (111) crystal surfaces.

Traffic flow is treated as a continuum governed by the kinematic LWR model and the Greenshield flux function. The model is modified to describe traffic flow on a road with traffic lights or a roundabout. Parameters introduced determine the traffic flow behaviour and queue formation, and numerical simulations based on the Godunov method are carried out. The numerical procedure is shown to converge, and produces results consistent with previous analytic results.

In this paper, we consider a two-point boundary value problem with Caputo fractional derivative, where the second order derivative of the exact solution is unbounded. Based on the equivalent form of the main equation, a finite difference scheme is derived. The L∞ convergence of the difference system is discussed rigorously. The convergence rate in general improves previous results. Numerical examples are provided to demonstrate the theoretical results.

This article is devoted to the study of some high-order difference schemes for the distributed-order time-fractional equations in both one and two space dimensions. Based on the composite Simpson formula and Lubich second-order operator, a difference scheme is constructed with convergence in the L1(L∞)-norm for the one-dimensional case, where τ,h and σ are the respective step sizes in time, space and distributed-order. Unconditional stability and convergence are proven. An ADI difference scheme is also derived for the two-dimensional case, and proven to be unconditionally stable and convergent in the L1(L∞)-norm, where h1 and h2 are the spatial step sizes. Some numerical examples are also given to demonstrate our theoretical results.

We propose amixed spectral method for heat transfer in unbounded domains, using generalised Hermite functions and Legendre polynomials. Some basic results on the mixed generalised Hermite-Legendre orthogonal approximation are established, which plays important roles in spectral methods for various problems defined on unbounded domains. As an example, the mixed generalised Hermite-Legendre spectral scheme is constructed for anisotropic heat transfer. Its convergence is proven, and some numerical results demonstrate the spectral accuracy of this approach.

High-dimensional two-sided space fractional diffusion equations with variable diffusion coefficients are discussed. The problems can be solved by an implicit finite difference scheme that is proven to be uniquely solvable, unconditionally stable and first-order convergent in the infinity norm. A nonsingular multilevel circulant pre-conditoner is proposed to accelerate the convergence rate of the Krylov subspace linear system solver efficiently. The preconditoned matrix for fast convergence is a sum of the identity matrix, a matrix with small norm, and a matrix with low rank under certain conditions. Moreover, the preconditioner is practical, with an O(N logN) operation cost and O(N) memory requirement. Illustrative numerical examples are also presented.