An efficient numerical method for 1D singularly perturbed parabolic convection-diffusion systems with repulsive interior turning points
In this work, we propose and study a numerical method to solve efficiently one-dimensional parabolic singularly perturbed systems of convection-diffusion type, for which the convection coefficient is zero at an interior point of the spatial domain. We focus our attention on the case of having the sa...
| Autores: | , |
|---|---|
| Tipo de recurso: | artículo |
| Estado: | Versión publicada |
| Fecha de publicación: | 2025 |
| País: | España |
| Institución: | Universidad Pública de Navarra |
| Repositorio: | Academica-e. Repositorio Institucional de la Universidad Pública de Navarra |
| OAI Identifier: | oai:academica-e.unavarra.es:2454/54351 |
| Acceso en línea: | https://hdl.handle.net/2454/54351 |
| Access Level: | acceso abierto |
| Palabra clave: | Singularly perturbed systems Turning points Boundary layers Fractional implicit Euler method Upwind scheme Shishkin meshes |
| Sumario: | In this work, we propose and study a numerical method to solve efficiently one-dimensional parabolic singularly perturbed systems of convection-diffusion type, for which the convection coefficient is zero at an interior point of the spatial domain. We focus our attention on the case of having the same diffusion parameter in both equations; as well we assume adequate signs on the convective coefficients in order to the interior turning point is of repulsive type. Under these conditions, if the data of the problem are composed by continuous functions, the exact evolutionary solution, in general, has regular boundary layers at the end points of the spatial domain. To solve this type of problems, we combine the fractional implicit Euler method and the classical upwind scheme, defined on a special mesh of Shishkin type. The resulting numerical method reach uniform convergence of first order in time and almost first order in space. Numerical results obtained for different test problems are shown which corroborate in practice the uniform convergence of the numerical algorithm and also their computational efficiency in comparison with classical numerical methods used for the same type of problems. |
|---|