Implementación del método LDG para mallas no estructuradas en 3D

Abstract

En este artículo se describe una implementación del método “Local Discontinuous Galerkin” (LDG) aplicado a problemas elípticos en 3D. Se discute la implementación de los principales operadores. En particular el uso de aproximaciones de alto orden y de mallas no estructuradas. Estructuras de datos eficientes que permiten un rápido ensamblado del sistema lineal en su formulación mixta son descritas en detalle.This paper describes an implementation of the Local Discontinuous Galerkin method (LDG) applied to elliptic problems in 3D. The implementation of the major operators is discussed. In particular the use of higher-order approximations and unstructured meshes. Efficient data structures that allow fast assembly of the linear system in the mixed formulation are described in detail.