The use of solutions on embedded grids for the approximation of singularly perturbed parabolic convection-diffusion equations on adapted grids

Abstract

The Dirichlet problem on a closed interval for a parabolic convection-diffusion equation is considered. The higher order derivative is multiplied by a parameter ɛ taking arbitrary values in the semi-open interval (0, 1]. For the boundary value problem, a finite difference scheme on a posteriori adapted grids is constructed. The classical approximations of the equation on uniform grids in the main domain are used; in some subdomains, these grids are subjected to refinement to improve the grid solution. The subdomains in which the grid should be refined are determined using the difference of the grid solutions of intermediate problems solved on embedded grids. Special schemes on a posteriori piecewise uniform grids are constructed that make it possible to obtain approximate solutions that converge almost ɛ-uniformly, i.e., with an error that weakly depends on the parameter ɛ: |u(x, t) − z(x, t)| ≤ M[N1−1 ln2N1 + N0−1 lnN0 + ɛ−1N1−K lnK−1N1], (x, t) e Ḡh, where N1 + 1 and N0 + 1 are the numbers of grid points in x and t, respectively; K is the number of refinement iterations (with respect to x) in the adapted grid; and M = M(K). Outside the σ-neighborhood of the outflow part of the boundary (in a neighborhood of the boundary layer), the scheme converges ɛ-uniformly at a rate O(N1−1 ln2N1 + N0−1 lnN0), where σ ≤ MN1−K + 1 lnK−1N1 for K ≥ 2.