In this paper, we propose a numerical scheme for a system of two linear singularly perturbed parabolic convection-diffusion equations. The presented numerical scheme consists of a classical backward-Euler scheme on a uniform mesh for the time discretization and an upwind finite difference scheme on an arbitrary nonuniform mesh for the spatial discretization. Then, for the time semidiscretization scheme, an a priori and an a posteriori error estimations in the maximum norm are obtained. It should be pointed out that the a posteriori error bound is suitable to design an adaptive algorithm, which is used to generate an adaptive spatial grid. It is proved that the method converges uniformly in the discrete maximum norm with first-order time and spatial accuracy, respectively, for the fully discrete scheme. At last, some numerical results are given to validate the theoretical results.
In this paper, an adaptive grid method for a singularly perturbed Volterra integro-differential equation is studied. Firstly, this problem is discretized by a new second-order finite difference scheme, for which a truncation error analysis is conducted. Then, based on this truncation error bound and the mesh equidistribution principle, we show that there is a mesh that provides an optimal error bound of O(N−2), which is robust with respect to the perturbation parameter. Finally, based on an approximation monitor function, an adaptive grid generation algorithm is constructed and some numerical results are given to support our theoretical results.
A nonlinear fractional differential equation with a Caputo derivative of order is studied. This problem is discretized by using the L1 scheme on an arbitrary nonuniform mesh. By utilizing the Taylor expansion with integral remainder term, an optimal local truncation error estimation of L1 scheme is proved. Based on this truncation error estimation and the mesh equidistribution principle, a new monitor function is constructed to construct an adaptive grid generation algorithm. Numerical experiments are performed to confirm the accuracy of our new adaptive grid algorithm.