Seminario del Departamento de
Ecuaciones Diferenciales y Análisis Numérico

Fecha : 12 de junio de 2014
Hora  : 10:30
Lugar : Seminario del Departamento (Fac. de Matemáticas, 3a. planta, módulo 34)

Giordano Tierra Chica

(Mathematical Institute of the Charles University, Praga, República Checa)

Resumen

In this seminar I will present the a priori and a posteriori error analyses of a mixed finite element method for Darcy's equations with porosity depending exponentially on the pressure. A simple change of variable for this unknown allows to transform the original nonlinear problem into a linear one whose dual-mixed variational formulation falls into the frameworks of the generalized linear saddle point problems and the fixed point equations satisfied by an affne mapping. According to the latter, we are able to show the well-posedness of both the continuous and discrete schemes, as well as the associated Cea estimate, by simply applying a suitable combination of the classical Babuska-Brezzi theory and the Banach fixed point Theorem. In particular, given any integer k ≥ 0, the stability of the Galerkin scheme is guaranteed by employing Raviart-Thomas elements of order k for the velocity, piecewise polynomials of degree k for the pressure, and continuous piecewise polynomials of degree k + 1 for an additional Lagrange multiplier given by the trace of the pressure on the Neumann boundary. Note that the two ways of writing the continuous formulation suggest accordingly two different methods for solving the discrete schemes. Next, we derive a reliable and efficient residual-based a posteriori error estimator for this problem. The global inf-sup condition satisfied by the continuous formulation, Helmholtz decompositions, and the local approximation properties of the Raviart-Thomas and Clément interpolation operators are the main tools for proving the reliability. In turn, inverse and discrete inequalities, and the localization technique based on triangle-bubble and edge-bubble functions are utilized to show the efficiency. Finally, several numerical results illustrating the good performance of both methods, confirrming the aforementioned properties of the estimator, and showing the behaviour of the associated adaptive algorithm, will be reported.
This contribution is based on recent joint work with Gabriel N. Gatica (Universidad de Concepción, Chile) and Ricardo Ruiz-Baier (University of Lausanne, Switzerland).