Grid approximation of singularly perturbed parabolic convection-diffusion equations subject to a piecewise smooth initial condition

Abstract

A boundary value problem for a singularly perturbed parabolic convection-diffusion equation on an interval is considered. The higher order derivative in the equation is multiplied by a parameter ɛ that can take arbitrary values in the half-open interval (0, 1]. The first derivative of the initial function has a discontinuity of the first kind at the point x0. For small values of ɛ, a boundary layer with the typical width of ɛ appears in a neighborhood of the part of the boundary through which the convective flow leaves the domain; in a neighborhood of the characteristic of the reduced equation outgoing from the point (x0, 0), a transient (moving in time) layer with the typical width of ɛ1/2 appears. Using the method of special grids that condense in a neighborhood of the boundary layer and the method of additive separation of the singularity of the transient layer, special difference schemes are designed that make it possible to approximate the solution of the boundary value problem ɛ-uniformly on the entire set \(\bar G\), approximate the diffusion flow (i.e., the product ɛ(ϱ/ϱx)u(x, t)) on the set \(\bar G^ * = \bar G\backslash \{ (x_0 ,0)\} \), and approximate the derivative (ϱ/ϱx)u(x, t) on the same set outside the m-neighborhood of the boundary layer. The approximation of the derivatives ɛ2(ϱ2/ϱx2)u(x, t) and (ϱ/ϱt)u(x, t) on the set \(\bar G^ * \) is also examined.