Implementación del método de estabilización de petrov-galerkin en contracorriente para la solución numérica de problemas altamente advectivos

Abstract

El presente artículo describe el método Petrov-Galerkin en contracorriente (SUPG) como técnica de estabilización de las solucio- nes por elementos finitos de la ecuación diferencial de difusión –advección- reacción. En la primera parte del artículo se hace un breve análisis de la importancia de este tipo de ecuaciones diferenciales para el modelado de fenómenos físicos en múltiples campos. Posteriormente, se realiza una descripción unidimensional del método SUPG y se desarrolla la metodología para im- plementar la técnica en dos o tres dimensiones. Se presentan los resultados de un experimento numérico fuertemente advéctico y de elevada complejidad desde el punto de vista numérico. Los resultados muestran cómo la versión de la técnica SUPG imple- mentada permite aproximaciones estabilizadas en el espacio, aun para altos números de Peclet. Gráficas adicionales de los ex- perimentos numéricos aquí presentados pueden ser descargados de www.gnum.unal.edu.co.

This article describes the streamline upwind Petrov-Galerkin (SUPG) method as being a stabilisation technique for resolving the diffusion-advection-reaction equation by finite elements. The first part of this article has a short analysis of the importance of this type of differential equation in modelling physical phenomena in multiple fields. A one-dimensional description of the SUPG me- thod is then given to extend this basis to two and three dimensions. The outcome of a strongly advective and a high numerical complexity experiment is presented. The results show how the version of the implemented SUPG technique allowed stabilised approaches in space, even for high Peclet numbers. Additional graphs of the numerical experiments presented here can be downloaded from www.gnum.unal.edu.co.