Orthogonal spline collocation scheme for multiterm fractional convection‐diffusion equation with variable coefficients

The orthogonal spline collocation (OSC) technique is an efficient way to solve a wide variety of problems that are modeled by ordinary and partial differential equations. In this article, by using OSC method in spatial direction and classical L1 approximation in temporal direction, a fully discrete scheme is established for a class of two‐dimensional multiterm fractional convection‐diffusion reaction equation with variable coefficients. The optimal estimates in H^j (j = 0, 1, 2) norms at each time step are derived. Also, estimate in space is provided. At last, we provide some numerical results to verify the accuracy and efficiency of the proposed algorithm.     

A backward euler orthogonal spline collocation method for the time‐fractional Fokker–Planck equation

We formulate and analyze a novel numerical method for solving a time‐fractional Fokker–Planck equation which models an anomalous subdiffusion process. In this method, orthogonal spline collocation is used for the spatial discretization and the time‐stepping is done using a backward Euler method based on the L1 approximation to the Caputo derivative. The stability and convergence of the method are considered, and the theoretical results are supported by numerical examples, which also exhibit superconvergence. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1534–1550, 2015     

An orthogonal spline collocation alternating direction implicit method for second-order hyperbolic problems

On a rectangular region, we consider a linear second-order hyperbolicinitial-boundary value problem involving a mixed derivativeterm, continuous variable coefficients and non-homogeneous Dirichletboundary conditions. In comparison to the alternating directionimplicit Laplace-modified method of Fernandes (1997), we formulateand analyse a new parameter-free alternating direction implicitscheme in which the standard central difference formula is usedfor the time approximation and orthogonal spline collocationis used for the spatial discretization. We establish unconditionalstability of the scheme, and its optimal order in the discretemaximum norm in time and the H¹ norm in space. Numerical experimentsindicate that the new scheme, which has the same order as themethod of Fernandes (1997, Numer. Math., 77, 223–241),is more accurate. We also show that the new scheme is easilygeneralized to the second-order hyperbolic problems on rectangularpolygons. Extensions of the scheme to problems with discontinuouscoefficients, nonlinear problems, and problems with other boundaryconditions are also discussed.     

Numerical simulation of a class of fractional subdiffusion equations via the alternating direction implicit method

In this article, a new numerical technique is proposed for solving the two‐dimensional time fractional subdiffusion equation with nonhomogeneous terms. After a transformation of the original problem, standard central difference approximation is used for the spatial discretization. For the time step, a new fractional alternating direction implicit (FADI) scheme based on the L₁ approximation is considered. This FADI scheme is constructed by adding a small term, so it is different from standard FADI methods. The solvability, unconditional stability and H¹ norm convergence are proved. Finally, numerical examples show the effectiveness and accuracy of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 531–547, 2016     

Alternating direction implicit OSC scheme for the two-dimensional fractional evolution equation with a weakly singular kernel

In this paper, a new kind of alternating direction implicit (ADI) Crank-Nicolson-type orthogonal spline collocation (OSC) method is formulated for the two-dimensional fractional evolution equation with a weakly singular kernel arising in the theory of linear viscoelasticity. The novel OSC method is used for the spatial discretization, and ADI Crank-Nicolson-type method combined with the second order fractional quadrature rule are considered for the temporal component. The stability of proposed scheme is rigourously established, and nearly optimal order error estimate is also derived. Numerical experiments are conducted to support the predicted convergence rates and also exhibit expected super-convergence phenomena.     

A second order finite difference-spectral method for space fractional diffusion equations

A high order finite difference-spectral method is derived for solving space fractional diffusion equations,by combining the second order finite difference method in time and the spectral Galerkin method in space.The stability and error estimates of the temporal semidiscrete scheme are rigorously discussed,and the convergence order of the proposed method is proved to be O(τ2+Nα-m)in L2-norm,whereτ,N,αand m are the time step size,polynomial degree,fractional derivative index and regularity of the exact solution,respectively.Numerical experiments are carried out to demonstrate the theoretical analysis.     

Superconvergence of the orthogonal spline collocation solution of Poisson's equation

Superconvergence phenomena have been observed numerically in the piecewise Hermite bicubic orthogonal spline collocation solution of Poisson's equation on a rectangle. The purpose of this article is to demonstrate theoretically the superconvergent fourth‐order accuracy in the first‐order partial derivatives of the collocation solution at the partition nodes. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 285–303, 1999     

A discrete collocation method for Fredholm integro-differential equations with weakly singular kernels

A fully discrete version of a piecewise polynomial collocation method is constructed to solve initial or boundary value problems of linear Fredholm integro-differential equations with weakly singular kernels. Using an integral equation reformulation and special graded grids, optimal global convergence estimates are derived. For special values of parameters an improvement of the convergence rate of elaborated numerical schemes is established. Some of our theoretical results are illustrated by numerical experiments.     

Numerical simulation for the three-dimension fractional sub-diffusion equation

Fractional sub-diffusion equations have been widely used to model sub-diffusive systems. Most algorithms are designed for one-dimensional problems due to the memory effect in fractional derivative. In this paper, the numerical simulation of the 3D fractional sub-diffusion equation with a time fractional derivative of order $\alpha$ $(0<\alpha <1)$ is considered. A fractional alternating direction implicit scheme (FADIS) is proposed. We prove that FADIS is uniquely solvable, unconditionally stable and convergent in $H^1$ norm by the energy method. A numerical example is given to demonstrate the efficiency of FADIS.     

An efficient nonpolynomial spline method for distributed order fractional subdiffusion equations

In this paper, we propose an efficient numerical method for a distributed order fractional subdiffusion problem using nonpolynomial spline approach. The solvability, stability, and convergence of the scheme are established rigorously, and it is shown that the spatial convergence order improves some previous work done. Simulation is then conducted to verify the accuracy of the proposed scheme as well as to compare with earlier work.     

A compact ADI scheme for the two dimensional time fractional diffusion‐wave equation in polar coordinates

In this article, we consider finite difference schemes for two dimensional time fractional diffusion‐wave equations on an annular domain. The problem is formulated in polar coordinates and, therefore, has variable coefficients. A compact alternating direction implicit scheme with accuracy order is derived, where τ, h₁, h₂ are the temporal and spatial step sizes, respectively. The stability and convergence of the proposed scheme are studied using its matrix form by the energy method. Numerical experiments are presented to support the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1692–1712, 2015     

A layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011     

High order convergent modified nodal bi-cubic spline collocation method for elliptic partial differential equation

A high order modified nodal bi-cubic spline collocation method is proposed for numerical solution of second-order elliptic partial differential equation subject to Dirichlet boundary conditions. The approximation is defined on a square mesh stencil using nine grid points. The solution of the method exists and is unique. Convergence analysis has been presented. Moreover, the superconvergent phenomena can be seen in proposed one step method. The numerical results clearly exhibit the superiority of the new approximation, in terms of both accuracy and computational efficiency.     

A second-order splitting combined with orthogonal cubic spline collocation method for the Rosenau equation

A second-order splitting method is applied to a KdV-like Rosenau equation in one space variable. Then an orthogonal cubic spline collocation procedure is employed to approximate the resulting system. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index 1. Error estimates in L² and L^∞ norms have been obtained for the semidiscrete approximations. For the temporal discretization, the time integrator RADAU5 is used for the resulting system. Some numerical experiments have been conducted to validate the theoretical results and to confirm the qualitative behaviors of the Rosenau equation. Finally, orthogonal cubic spline collocation method is directly applied to BBM (Benjamin–Bona–Mahony) and BBMB (Benjamin–Bona–Mahony–Burgers) equations and the well-known decay estimates are demonstrated for the computed solution. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 695–716, 1998     

Orthogonal spline collocation methods for the stream function‐vorticity formulation of the Navier–Stokes equations

An orthogonal spline collocation (OSC) spatial discretization is proposed for the solution of the fully coupled stream function‐vorticity formulation of the Navier–Stokes equations in two dimensions. For the time‐stepping, a three‐level leapfrog scheme is employed. This method is algebraically linear, and, at each time step, gives rise to a system of linear equations of the form arising in the OSC approximation of the biharmonic Dirichlet problem and can be solved by a fast direct method. Error estimates in the H^l–norm in space, l = 1,2, are derived for the semi‐discrete method and the fully‐discrete leapfrog scheme which is also shown to be second order accurate in time. Numerical results are presented which confirm the theoretical analysis and exhibit superconvergence phenomena, which provide superconvergent approximations to the components of the velocity. © John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Finite Difference/Collocation Method for Two-Dimensional Sub-Diffusion Equation with Generalized Time Fractional Derivative

In this paper, we propose a finite difference/collocation method for two-dimensional time fractional diffusion equation with generalized fractional operator. The main purpose of this paper is to design a high order numerical scheme for the new generalized time fractional diffusion equation. First, a finite difference approximation formula is derived for the generalized time fractional derivative, which is verified with order $2-\alpha$ $(0<\alpha<1)$. Then, collocation method is introduced for the two-dimensional space approximation. Unconditional stability of the scheme is proved. To make the method more efficient, the alternating direction implicit method is introduced to reduce the computational cost. At last, numerical experiments are carried out to verify the effectiveness of the scheme.     

A second-order numerical method for two-dimensional two-sided space fractional convection diffusion equation

Space fractional convection diffusion equation describes physical phenomena where particles or energy (or other physical quantities) are transferred inside a physical system due to two processes: convection and superdiffusion. In this paper, we discuss the practical alternating directions implicit method to solve the two-dimensional two-sided space fractional convection diffusion equation on a finite domain. We theoretically prove and numerically verify that the presented finite difference scheme is unconditionally von Neumann stable and second order convergent in both space and time directions.