Ref: RR254D
Revista Internacional Métodos numéricos para cálculo y diseño en ingeniería Vol25 Num: 4
ESTIMA DE ERROR RESIDUAL EXPLÍCITA PARA CANTIDADES DE INTERÉS UTILIZANDO FUNCIONES BURBUJA
En este trabajo se introduce un nuevo estimador de error residual explícito a posteriori para problemas
elípticos orientado a cantidades de interés. Se propone utilizar funciones burbuja sobre elementos y sobre
aristas. Se parte de la solución de elementos finitos del problema primal y de la de un problema adjunto (o
dual), asociado a una cantidad de interés dfinida por el usuario. Por ejemplo, la variación de temperatura
o el desplazamiento de un punto del dominio.
La estima se calcula en dos fases. Primero se aproxima el error en el interior de los elementos utilizando una
función burbuja. La aproximación del error interior sobre todo el dominio sería dada como una combinación
lineal de estas funciones. En segundo lugar, se introduce otra familia de funciones burbuja, asociada a las
aristas de la malla (lados de los elementos) y se complementa la primera fase de la estima añadiendo la
contribución al error que proviene de las aristas.
A diferencia de los estimadores del error explícitos habituales, el que se propone aquí no depende de constantes desconocidas. Esto es posible debido a que a pesar del carácter explícito del estimador, se está
aproximando la función de error asociada a uno de los problemas (el directo o el adjunto) en un espacio de
interpolación muy sencillo (el que dfinen las funciones burbuja). El cálculo explícito del residuo corresponde
sólo al otro de los dos problemas.
Summary
A new goal-oriented residual a posteriori error estimator is introduced for elliptic problems. The proposed
estimate is based on using bubble functions over both elements and edges. The error representation of
the quantity of interest prescribed by the user is performed using an adjoint problem. These quantities of
interest are either averaged values of the solution (mean temperature) or point values as, for instance, the
displacement of a given point.
The estimation procedure is organized in two phases. First, the error is approximated in the interior of
the elements projecting it into a given bubble function. The first contribution to the overall error (interior
estimate) is just the sum of all these projections. Second, a new family of bubble functions is considered,
each associated with one edge of the finite element mesh. The remaining part of the error is projected into
this new family of functions. This second error contribution (edge estimate) is added to the previous one to
obtain the complete estimate.
The present estimator is independent of the typical unknown constants appearing in explicit residual type
estimates. Despite of its explicit character, the present estimator gets rid of these constants by properly
combining the residuals in the adjoint and direct problems, although both are explicitly computed.
|