Conservative numerical method for a system of semilinear singularly perturbed parabolic reaction‐diffusion equations

Abstract

On a vertical strip, a Dirichlet problem is considered for a system of two semilinear singularly perturbed parabolic reaction‐diffusion equations connected only by terms that do not involve derivatives. The highest‐order derivatives in the equations, having divergent form, are multiplied by the perturbation parameter e ; ϵ ∈ (0,1]. When ϵ → 0, the parabolic boundary layer appears in a neighbourhood of the strip boundary.

Using the integro‐interpolational method, conservative nonlinear and linearized finite difference schemes are constructed on piecewise‐uniform meshes in the x ₁ ‐axis (orthogonal to the boundary) whose solutions converge ϵ‐uniformly at the rate O (N1⁻² ln² N ₁ + N ₂ ⁻² + N ₀ ⁻¹). Here N ₁ + 1 and N ₀ + 1 denote the number of nodes on the x ₁‐axis and t‐axis, respectively, and N ₂ + 1 is the number of nodes in the x ₂‐axis on per unit length.

First published online: 14 Oct 2010