An efficient numerical method to solve 2D parabolic singularly perturbed coupled systems of convection-diffusion type with multi-parameters on a Bakhvalov–Shishkin mesh

Shiromani, Ram ; Clavero, Carmelo (Universidad de Zaragoza)

Resumen: This study addresses the efficient solution of a class of 2D parabolic singularly perturbed weakly coupled systems of convection-diffusion type. In the model problem, small positive parameters appear in both the diffusion and the convection terms. We assume that the diffusion parameters can be distinct, but the convection parameter remains the same for both equations. Then, for sufficiently small values of the parameters, overlapping boundary layers appear on the boundary of the spatial domain. To solve the problem, a numerical method is employed that combines the implicit Euler scheme, defined on a uniform mesh, with the upwind scheme for spatial discretization. Then, if the spatial discretization is carried out on an adequate nonuniform Bakhvalov–Shishkin (BS) mesh, the fully discrete scheme attains uniform convergence, with respect to all perturbation parameters; moreover, it has first-order accuracy in both temporal and spatial variables. Note that the construction of the BS mesh depends on the value and the ratio between the diffusion and the convection parameters, and special generating functions are needed to construct them. Numerical experiments illustrating the performance of the algorithm for some test problems are showed, which corroborate the uniform convergence of the method in agreement with the theoretical results.

Idioma: Inglés

DOI: 10.3934/math.2026076

Año: 2026

Publicado en: AIMS Mathematics 11, 1 (2026), 1820-1856
ISSN: 2473-6988

Financiación: info:eu-repo/grantAgreement/ES/DGA-FSE/E24-17R
Financiación: info:eu-repo/grantAgreement/ES/MCINN/PID2022-136441NB-I00

Tipo y forma: Article (Published version)

Área (Departamento): Área Matemática Aplicada (Dpto. Matemática Aplicada)