Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems on Anisotropically Refined Meshes

In this paper we consider the a posteriori and a priori error analysis of discontinuous Galerkin interior penalty methods for second-order partial differential equations with nonnegative characteristic form on anisotropically refined computational meshes. In particular, we discuss the question of error estimation for linear target functionals, such as the outflow flux and the local average of the solution. Based on our a posteriori error bound we design and implement the corresponding adaptive algorithm to ensure reliable and efficient control of the error in the prescribed functional to within a given tolerance. This involves exploiting both local isotropic and anisotropic mesh refinement. The theoretical results are illustrated by a series of numerical experiments

Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems on Anisotropically Refined Meshes

In this paper we consider the a posteriori and a priori error analysis of discontinuous Galerkin interior penalty methods for second-order partial differential equations with nonnegative characteristic form on anisotropically refined computational meshes. In particular, we discuss the question of error estimation for linear target functionals, such as the outflow flux and the local average of the solution. Based on our a posteriori error bound we design and implement the corresponding adaptive algorithm to ensure reliable and efficient control of the error in the prescribed functional to within a given tolerance. This involves exploiting both local isotropic and anisotropic mesh refinement. The theoretical results are illustrated by a series of numerical experiments

A Pre-processing Moving Mesh Method for Discontinuous Galerkin Approximations of Advection-Diffusion-Reaction Problems

We propose a pre-processing mesh re-distribution algorithm based upon harmonic maps employed in conjunction with discontinuous Galerkin approximations of advection-diffusion-reaction problems. Extensive two-dimensional numerical experiments with different choices of monitor functions, including monitor functions derived from goal-oriented a posteriori error indicators are presented. The examples presented clearly demonstrate the capabilities and the benefits of combining our pre-processing mesh movement algorithm with both uniform, as well as, adaptive isotropic and anisotropic mesh refinement

hp-Version Composite Discontinuous Galerkin Methods for Elliptic Problems on Complicated Domains

In this paper we introduce the hp-version discontinuous Galerkin composite finite element method for the discretization of second-order elliptic partial differential equations. This class of methods allows for the approximation of problems posed on computational domains which may contain a huge number of local geometrical features, or microstructures. While standard numerical methods can be devised for such problems, the computational effort may be extremely high, as the minimal number of elements needed to represent the underlying domain can be very large. In contrast, the minimal dimension of the underlying composite finite element space is independent of the number of geometric features. The key idea in the construction of this latter class of methods is that the computational domain Ω is no longer resolved by the mesh; instead, the finite element basis (or shape) functions are adapted to the geometric details present in Ω. In this paper, we extend these ideas to the discontinuous Galerkin setting, based on employing the hp-version of the finite element method. Numerical experiments highlighting the practical application of the proposed numerical scheme will be presented

Exponentially Fitted Discontinuous Galerkin Schemes for Singularly Perturbed Problems

New discontinuous Galerkin schemes in mixed form are introduced for symmetric elliptic problems of second order. They exhibit reduced connectivity with respect to the standard ones. The modifications in the choice of the approximation spaces and in the stabilization term do not spoil the error estimates. These methods are then used for designing new exponentially fitted schemes for advection dominated equations. The presented numerical tests show the good performances of the proposed schemes.Fil: Lombardi, Ariel Luis. Universidad Nacional de General Sarmiento. Instituto de Ciencias; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Matemática; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Pietra, Paola. Consiglio Nazionale delle Ricerche. Istituto di Matematica Applicata e Tecnologie Informatiche; Itali