Implementation of LDG method for 3D unstructured meshes

Abstract

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.Keywords: Discontinuous finite element methods, high-order approximations, unstructured meshes, object-oriented programming.Mathematics Subject Classification: 65K05, 65N30, 65N55.

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.Palabras clave: Métodos de elemento finito discontinuos, aproximaciones de alto orden, mallas no estructuradas, programación orientada a objetos.Mathematics Subject Classification: 65K05, 65N30, 65N55.