Una implementación computacional del método Vem Mixto para el problema de Brinkman en 2D

Abstract

En este artículo se describen algunos aspectos específicos sobre una implementación computacional para la formulación mixta de elementos virtuales (mixed-VEM, por sus siglas en inglés) del problema lineal de Brinkman en dos dimensiones, con condiciones de frontera de Dirichlet no homogéneas. La formulación empleada fue originalmente propuesta y analizada en [18]. La implementación planteada aquí considera cualquier grado polinomial k ≥ 0 de manera natural al construir diversas matrices locales de bajo tamaño. Además, se propone un algoritmo para el ensamblaje del sistema lineal global asociado, que garantiza la continuidad de la componente normal en la solución discreta.In this paper we describe some specific aspects on the computational implementation of the a mixed virtual element method (mixed-VEM) for the two-dimensional linear Brinkman model with non-homogeneous Dirichlet boundary conditions. The formulation used below was originally proposed and analysed in [18]. The implementation presented here consider any polynomial degree k ≥ 0 in a natural way by building several local matrices of small size. In addition, an algorithm is proposed for the assembly of the associated global linear system, which guarantees the continuity of the normal component in the discrete solution.