A parameter‐uniform scheme for the parabolic singularly perturbed problem with a delay in time

In this paper, a parameter‐uniform numerical scheme for the solution of singularly perturbed parabolic convection–diffusion problems with a delay in time defined on a rectangular domain is suggested. The presence of the small diffusion parameter ? leads to a parabolic right boundary layer. A collocation method consisting of cubic B ‐spline basis functions on an appropriate piecewise‐uniform mesh is used to discretize the system of ordinary differential equations obtained by using Rothe's method on an equidistant mesh in the temporal direction. The parameter‐uniform convergence of the method is shown by establishing the theoretical error bounds. The numerical results of the test problems validate the theoretical error bounds.     

High‐order numerical methods for one‐dimensional parabolic singularly perturbed problems with regular layers

In this work we construct and analyze some finite difference schemes used to solve a class of time‐dependent one‐dimensional convection‐diffusion problems, which present only regular layers in their solution. We use the implicit Euler or the Crank‐Nicolson method to discretize the time variable and a HODIE finite difference scheme, defined on a piecewise uniform Shishkin mesh, to discretize the spatial variable. In both cases we prove that the numerical method is uniformly convergent with respect to the diffusion parameter, having order near two in space and order one or 3/2, depending on the method used, in time. We show some numerical examples which illustrate the theoretical results, in the case of using the Euler implicit method, and give better numerical behaviour than that predicted theoretically, showing order two in time and order N^?2log²N in space, if the Crank‐Nicolson scheme is used to discretize the time variable. Finally, we construct a numerical algorithm by combining a third order A‐stable SDIRK with two stages and a third‐order HODIE difference scheme, showing its uniformly convergent behavior, reaching order three, up to a logarithmic factor. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Numerical approximation of two‐dimensional convection‐diffusion equations with boundary layers

Our aim in this article is to show how one can improve the numerical solution of singularity perturbed problems involving boundary layers. Incorporating the structures of boundary layers into finite element spaces can improve the accuracy of approximate solutions and result in significant simplifications. In this article we discuss convection‐diffusion equations in the two‐dimensional space with a homogeneous Dirichlet boundary condition and a mixed boundary condition. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Multilevel numerical solutions of convection‐dominated diffusion problems by spline wavelets

In this article, we utilize spline wavelets to establish an adaptive multilevel numerical scheme for time‐dependent convection‐dominated diffusion problems within the frameworks of Galerkin formulation and Eulerian‐Lagrangian localized adjoint methods (ELLAM). In particular, we shall use linear Chui‐Quak semi‐orthogonal wavelets, which have explicit expressions and compact supports. Therefore, both the diffusion term and boundary conditions in the convection‐diffusion problems can be readily handled. Strategies for efficiently implementing the scheme are discussed and numerical results are interpreted from the viewpoint of nonlinear approximation. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Higher order numerical approximation for time dependent singularly perturbed differential‐difference convection‐diffusion equations

In this article, we develop a higher order numerical approximation for time dependent singularly perturbed differential‐difference convection‐diffusion equations. A priori bounds on the exact solution and its derivatives, which are useful for the error analysis of the numerical method are given. We approximate the retarded terms of the model problem using Taylor's series expansion and the resulting time‐dependent singularly perturbed problem is discretized by the implicit Euler scheme on uniform mesh in time direction and a special hybrid finite difference scheme on piecewise uniform Shishkin mesh in spatial direction. We first prove that the proposed numerical discretization is uniformly convergent of , where and denote the time step and number of mesh‐intervals in space, respectively. After that we design a Richardson extrapolation scheme to increase the order of convergence in time direction and then the new scheme is proved to be uniformly convergent of . Some numerical tests are performed to illustrate the high‐order accuracy and parameter uniform convergence obtained with the proposed numerical methods.     

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     

Graded Galerkin methods for the high‐order convection‐diffusion problem

We develop a Galerkin method using the Hermite spline on an admissible graded mesh for solving the high‐order singular perturbation problem of the convection‐diffusion type. We identify a special function class to which the solution of the convection‐diffusion problem belongs and characterize the approximation order of the Hermite spline for such a function class. The approximation order is then used to establish the optimal order of uniform convergence for the Galerkin method. Numerical results are presented to confirm the theoretical estimate.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A high‐order Godunov method for one‐dimensional convection‐diffusion‐reaction problems

In this article we develop a high‐order Godunov method for one‐dimensional convection‐diffusion‐reaction problems where convection dominates diffusion. The heart of this method comes from incorporating the diffusion term via the slope of the linear representation (recovery) of the solution on each grid cell. The method is conservative and explicit. Therefore, it is efficient in computing time. For constant coefficient linear convection, diffusion, and Lipschitz‐type reaction, the properties of the total variation stability and monotonicity preservation are proved. An error estimation is derived. Computational examples are presented and compared with the exact solutions. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 495–512, 2000     

Parameter uniform approximations for time‐dependent reaction‐diffusion problems

Using discrete Green's functions techniques, we present a classification of fitted mesh methods for time‐dependent reaction diffusion problems, based on the analyses of Linß (Linß, Numer Algor 40 (2005), 23–32) for the analogous steady‐state problem and of Kopteva (Kopteva, Computing 66 (2001), 179–197) for time‐dependent convection‐diffusion problems. As examples of how to apply the analysis, we derive error estimates for the fitted meshes of Shishkin and Bakhvalov, and provide supporting numerical results. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

One‐level Newton–Krylov–Schwarz algorithm for unsteady non‐linear radiation diffusion problem

In this paper, we present a parallel Newton–Krylov–Schwarz (NKS)‐based non‐linearly implicit algorithm for the numerical solution of the unsteady non‐linear multimaterial radiation diffusion problem in two‐dimensional space. A robust solver technology is required for handling the high non‐linearity and large jumps in material coefficients typically associated with simulations of radiation diffusion phenomena. We show numerically that NKS converges well even with rather large inflow flux boundary conditions. We observe that the approach is non‐linearly scalable, but not linearly scalable in terms of iteration numbers. However, CPU time is more important than the iteration numbers, and our numerical experiments show that the algorithm is CPU‐time‐scalable even without a coarse space given that the mesh is fine enough. This makes the algorithm potentially more attractive than multilevel methods, especially on unstructured grids, where course grids are often not easy to construct. Copyright © 2004 John Wiley & Sons, Ltd.     

Superconvergence of the direct discontinuous Galerkin method for a time‐fractional initial‐boundary value problem

A time‐fractional reaction–diffusion initial‐boundary value problem with periodic boundary condition is considered on Q ? Ω × [0, T] , where Ω is the interval [0, l] . Typical solutions of such problem have a weak singularity at the initial time t = 0. The numerical method of the paper uses a direct discontinuous Galerkin (DDG) finite element method in space on a uniform mesh, with piecewise polynomials of degree k ≥ 2 . In the temporal direction we use the L1 approximation of the Caputo derivative on a suitably graded mesh. We prove that at each time level of the mesh, our L1‐DDG solution is superconvergent of order k + 2 in L²(Ω) to a particular projection of the exact solution. Moreover, the L1‐DDG solution achieves superconvergence of order (k + 2) in a discrete L²(Q) norm computed at the Lobatto points, and order (k + 1) superconvergence in a discrete H¹(Q) seminorm at the Gauss points; numerical results show that these estimates are sharp.     

Analysis of a decoupled time‐stepping algorithm for reduced MHD system modeling magneto‐convection

In this article, we analyze the stability and error estimate of a decoupled algorithm for a magneto‐convection problem. Magneto‐convection is assumed to be modeled by a coupled system of reduced magneto‐hydrodynamic (RMHD) equations and convection‐diffusion equation. The proposed algorithm applies the second‐order backward difference formula in time and finite element in space. To obtain a noniterative decouple algorithm from the fully discrete nonlinear system, we use a second‐order extrapolation in time to the nonlinear terms such that their skew symmetry properties are preserved. We prove the stability of the algorithm and derive error estimates without assuming any stability conditions. The algorithm is unconditionally stable and requires the solution of one RMHD problem and one convection‐diffusion equation per time step. Numerical test is presented that illustrates the accuracy and efficiency of the algorithm.     

An IMEX‐BDF2 compact scheme for pricing options under regime‐switching jump‐diffusion models

In this paper, an implicit‐explicit two‐step backward differentiation formula (IMEX‐BDF2) together with finite difference compact scheme is developed for the numerical pricing of European and American options whose asset price dynamics follow the regime‐switching jump‐diffusion process. It is shown that IMEX‐BDF2 method for solving this system of coupled partial integro‐differential equations is stable with the second‐order accuracy in time. On the basis of IMEX‐BDF2 time semi‐discrete method, we derive a fourth‐order compact (FOC) finite difference scheme for spatial discretization. Since the payoff function of the option at the strike price is not differentiable, the results show only second‐order accuracy in space. To remedy this, a local mesh refinement strategy is used near the strike price so that the accuracy achieves fourth order. Numerical results illustrate the effectiveness of the proposed method for European and American options under regime‐switching jump‐diffusion models.     

Two‐grid methods for time‐harmonic Maxwell equations

In this paper, we develop several two‐grid methods for the Nédélec edge finite element approximation of the time‐harmonic Maxwell equations. We first present a two‐grid method that uses a coarse space to solve the original problem and then use a fine space to solve a corresponding symmetric positive definite problem. Then, we present two types of iterative two‐grid methods, one is to add the kernel of the curl ‐operator in the fine space to a coarse mesh space to solve the original problem and the other is to use an inner iterative method for dealing with the kernel of the curl ‐operator in the fine space and the coarse space, separately. We provide the error estimates for the first two methods and present numerical experiments to show the efficiency of our methods.Copyright © 2012 John Wiley & Sons, Ltd.     

A high‐order compact ADI method for solving three‐dimensional unsteady convection‐diffusion problems

We derive a high‐order compact alternating direction implicit (ADI) method for solving three‐dimentional unsteady convection‐diffusion problems. The method is fourth‐order in space and second‐order in time. It permits multiple uses of the one‐dimensional tridiagonal algorithm with a considerable saving in computing time and results in a very efficient solver. It is shown through a discrete Fourier analysis that the method is unconditionally stable in the diffusion case. Numerical experiments are conducted to test its high order and to compare it with the standard second‐order Douglas‐Gunn ADI method and the spatial fourth‐order compact scheme by Karaa. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Parallel characteristic finite element method for time‐dependent convection–diffusion problem

Based on the overlapping‐domain decomposition and parallel subspace correction method, a new parallel algorithm is established for solving time‐dependent convection–diffusion problem with characteristic finite element scheme. The algorithm is fully parallel. We analyze the convergence of this algorithm, and study the dependence of the convergent rate on the spacial mesh size, time increment, iteration times and sub‐domains overlapping degree. Both theoretical analysis and numerical results suggest that only one or two iterations are needed to reach to optimal accuracy at each time step. Copyright © 2010 John Wiley & Sons, Ltd.     

A higher order difference method for singularly perturbed parabolic partial differential equations

This paper studies a higher order numerical method for the singularly perturbed parabolic convection-diffusion problems where the diffusion term is multiplied by a small perturbation parameter. In general, the solutions of these type of problems have a boundary layer. Here, we generate a spatial adaptive mesh based on the equidistribution of a positive monitor function. Implicit Euler method is used to discretize the time variable and an upwind scheme is considered in space direction. A higher order convergent solution with respect to space and time is obtained using the postprocessing based extrapolation approach. It is observed that the convergence is independent of perturbation parameter. This technique enhances the order of accuracy from first order uniform convergence to second order uniform convergence in space as well as in time. Comparative study with the existed meshes show the highly effective behavior of the present method.     

A space‐time spectral method for one‐dimensional time fractional convection diffusion equations

In this paper, we propose a space‐time spectral method for solving a class of time fractional convection diffusion equations. Because both fractional derivative and spectral method have global characteristics in bounded domains, we propose a space‐time spectral‐Galerkin method. The convergence result of the method is proved by providing a priori error estimate. Numerical results further confirm the expected convergence rate and illustrate the versatility of our method. Copyright © 2016 John Wiley & Sons, Ltd.     

Solving time‐periodic fractional diffusion equations via diagonalization technique and multigrid

This paper addresses numerical computation of time‐periodic diffusion equations with fractional Laplacian. Time‐periodic differential equations present fundamental challenges for numerical computation because we have to consider all the discrete solutions once in all instead of one by one. An idea based on the diagonalization technique is proposed, which yields a direct parallel‐in‐time computation for all the discrete solutions. The major computation cost is therefore reduced to solve a series of independent linear algebraic systems with complex coefficients, for which we apply a multigrid method using the damped Richardson iteration as the smoother. Such a linear solver possesses mesh‐independent convergence factor, and we make an optimization for the damping parameter to minimize such a constant convergence factor. Numerical results are provided to support our theoretical analysis.     

High-order finite difference schemes for the solution of second-order BVPs

We introduce new methods in the class of boundary value methods (BVMs) to solve boundary value problems (BVPs) for a second-order ODE. These formulae correspond to the high-order generalizations of classical finite difference schemes for the first and second derivatives. In this research, we carry out the analysis of the conditioning and of the time-reversal symmetry of the discrete solution for a linear convection–diffusion ODE problem. We present numerical examples emphasizing the good convergence behavior of the new schemes. Finally, we show how these methods can be applied in several space dimensions on a uniform mesh.