Cette thèse porte sur la modélisation de l’écoulement et du transport en milieu poreux ;nous effectuons des simulations numériques et démontrons des résultats de convergence d’algorithmes.Au Chapitre 1, nous appliquons des méthodes de volumes finis pour la simulation d’écoulements à densité variable en milieu poreux ; il vient à résoudre une équation de convection diffusion parabolique pour la concentration couplée à une équation elliptique en pression.Nous nous appuyons sur la méthode des volumes finis standard pour le calcul des solutions de deux problèmes spécifiques : une interface en rotation entre eau salée et eau douce et le problème de Henry. Nous appliquons ensuite la méthode de volumes finis généralisés SUSHI pour la simulation des mêmes problèmes ainsi que celle d’un problème de bassin salé en dimension trois d’espace. Nous nous appuyons sur des maillages adaptatifs, basés sur des éléments de volume carrés ou cubiques.Au Chapitre 2, nous nous appuyons de nouveau sur la méthode de volumes finis généralisés SUSHI pour la discrétisation de l’équation de Richards, une équation elliptique parabolique pour le calcul d’écoulements en milieu poreux. Le terme de diffusion peut être anisotrope et hétérogène. Cette classe de méthodes localement conservatrices s’applique àune grande variété de mailles polyédriques non structurées qui peuvent ne pas se raccorder.La discrétisation en temps est totalement implicite. Nous obtenons un résultat de convergence basé sur des estimations a priori et sur l’application du théorème de compacité de Fréchet-Kolmogorov. Nous présentons aussi des tests numériques.Au Chapitre 3, nous discrétisons le problème de Signorini par un schéma de type gradient,qui s’écrit à l’aide d’une formulation variationnelle discrète et est basé sur des approximations indépendantes des fonctions et des gradients. On montre l’existence et l’unicité de la solution discrète ainsi que sa convergence vers la solution faible du problème continu. Nous présentons ensuite un schéma numérique basé sur la méthode SUSHI.Au Chapitre 4, nous appliquons un schéma semi-implicite en temps combiné avec la méthode SUSHI pour la résolution numérique d’un problème d’écoulements à densité variable ;il s’agit de résoudre des équations paraboliques de convection-diffusion pour la densité de soluté et le transport de la température ainsi que pour la pression. Nous simulons l’avance d’un front d’eau douce assez chaude et le transport de chaleur dans un aquifère captif qui est initialement chargé d’eau froide salée. Nous utilisons des maillages adaptatifs, basés sur des éléments de volume carrés. / This thesis bears on the modelling of groundwater flow and transport in porous media; we perform numerical simulations by means of finite volume methods and prove convergence results. In Chapter 1, we first apply a semi-implicit standard finite volume method and then the generalized finite volume method SUSHI for the numerical simulation of density driven flows in porous media; we solve a nonlinear convection-diffusion parabolic equation for the concentration coupled with an elliptic equation for the pressure. We apply the standard finite volume method to compute the solutions of a problem involving a rotating interface between salt and fresh water and of Henry's problem. We then apply the SUSHI scheme to the same problems as well as to a three dimensional saltpool problem. We use adaptive meshes, based upon square volume elements in space dimension two and cubic volume elements in space dimension three. In Chapter 2, we apply the generalized finite volume method SUSHI to the discretization of Richards equation, an elliptic-parabolic equation modeling groundwater flow, where the diffusion term can be anisotropic and heterogeneous. This class of locally conservative methods can be applied to a wide range of unstructured possibly non-matching polyhedral meshes in arbitrary space dimension. As is needed for Richards equation, the time discretization is fully implicit. We obtain a convergence result based upon a priori estimates and the application of the Fréchet-Kolmogorov compactness theorem. We implement the scheme and present numerical tests. In Chapter 3, we study a gradient scheme for the Signorini problem. Gradient schemes are nonconforming methods written in discrete variational formulation which are based on independent approximations of the functions and the gradients. We prove the existence and uniqueness of the discrete solution as well as its convergence to the weak solution of the Signorini problem. Finally we introduce a numerical scheme based upon the SUSHI discretization and present numerical results. In Chapter 4, we apply a semi-implicit scheme in time together with a generalized finite volume method for the numerical solution of density driven flows in porous media; it comes to solve nonlinear convection-diffusion parabolic equations for the solute and temperature transport as well as for the pressure. We compute the solutions for a specific problem which describes the advance of a warm fresh water front coupled to heat transfer in a confined aquifer which is initially charged with cold salt water. We use adaptive meshes, based upon square volume elements in space dimension two.

Effiziente und zuverlässige a posteriori Fehlerabschätzungen sind eine Hauptzutat für die effiziente numerische Berechnung von Lösungen zu Variationsungleichungen durch die Finite-Elemente-Methode. Die vorliegende Arbeit untersucht zuverlässige und effiziente Fehlerabschätzungen für beliebige Finite-Elemente-Methoden und drei Variationsungleichungen, nämlich dem Hindernisproblem, dem Signorini Problem und dem Bingham Problem in zwei Raumdimensionen. Die Fehlerabschätzungen hängen vom zum Problem gehörenden Lagrange Multiplikator ab, der eine Verbindung zwischen der Variationsungleichung und dem zugehörigen linearen Problem darstellt. Effizienz und Zuverlässigkeit werden bezüglich eines totalen Fehlers gezeigt. Die Fehleranschätzungen fordern minimale Regularität. Die Approximation der exakten Lösung erfüllt die Dirichlet Randbedingungen und die Approximation des Lagrange Multiplikators ist nicht-positiv im Falle des Hindernis- und Signoriniproblems, und hat Betrag kleiner gleich 1 für das Bingham Problem. Dieses allgemeine Vorgehen ermöglicht das Einbinden nicht-exakter diskreter Lösungen, welche im Kontext dieser Ungleichungen auftreten. Aus dem Blickwinkel der Anwendungen ist Effizienz und Zuverlässigkeit im Bezug auf den Fehler der primalen Variablen in der Energienorm von großem Interesse. Solche Abschätzungen hängen von der Wahl eines effizienten diskreten Lagrange Multiplikators ab. Im Falle des Hindernis- und Signorini Problems werden postive Beispiele für drei Finite-Elemente Methoden, der konformen Courant Methode, der nicht-konformen Crouzeix-Raviart Methode und der gemischten Raviart-Thomas Methode niedrigster Ordnung hergeleitet. Partielle Resultate liegen im Fall des Bingham Problems vor. Numerischer Experimente heben die theoretischen Ergebnisse hervor und zeigen Effizienz und Zuverlässigkeit. Die numerischen Tests legen nahe, dass der aus den Abschätzungen resultierende adaptive Algorithmus mit optimaler Konvergenzrate konvergiert. / Efficient and reliable a posteriori error estimates are a key ingredient for the efficient numerical computation of solutions for variational inequalities by the finite element method. This thesis studies such reliable and efficient error estimates for arbitrary finite element methods and three representative variational inequalities, namely the obstacle problem, the Signorini problem, and the Bingham problem in two space dimensions. The error estimates rely on a problem connected Lagrange multiplier, which presents a connection between the variational inequality and the corresponding linear problem. Reliability and efficiency are shown with respect to some total error. Reliability and efficiency are shown under minimal regularity assumptions. The approximation to the exact solution satisfies the Dirichlet boundary conditions, and an approximation of the Lagrange multiplier is non-positive in the case of the obstacle and Signorini problem and has an absolute value smaller than 1 for the Bingham flow problem. These general assumptions allow for reliable and efficient a posteriori error analysis even in the presence of inexact solve, which naturally occurs in the context of variational inequalities. From the point of view of the applications, reliability and efficiency with respect to the error of the primal variable in the energy norm is of great interest. Such estimates depend on the efficient design of a discrete Lagrange multiplier. Affirmative examples of discrete Lagrange multipliers are presented for the obstacle and Signorini problem and three different first-order finite element methods, namely the conforming Courant, the non-conforming Crouzeix-Raviart, and the mixed Raviart-Thomas FEM. Partial results exist for the Bingham flow problem. Numerical experiments highlight the theoretical results, and show efficiency and reliability. The numerical tests suggest that the resulting adaptive algorithms converge with optimal convergence rates.

The main emphasis of this thesis is to study a posteriori error analysis of discontinuous Galerkin (DG) methods for the elliptic variational inequalities. The DG methods have become very pop-ular in the last two decades due to its nature of handling complex geometries, allowing irregular meshes with hanging nodes and different degrees of polynomial approximation on different ele-ments. Moreover they are high order accurate and stable methods. Adaptive algorithms refine the mesh locally in the region where the solution exhibits irregular behaviour and a posteriori error estimates are the main ingredients to steer the adaptive mesh refinement. The solution of linear elliptic problem exhibits singularities due to change in boundary con-ditions, irregularity of coefficients and reentrant corners in the domain. Apart from this, the solu-tion of variational inequality exhibits additional irregular behaviour due to occurrence of the free boundary (the part of the domain which is a priori unknown and must be found as a component of the solution). In the lack of full elliptic regularity of the solution, uniform refinement is inefficient and it does not yield optimal convergence rate. But adaptive refinement, which is based on the residuals ( or a posteriori error estimator) of the problem, enhance the efficiency by refining the mesh locally and provides the optimal convergence. In this thesis, we derive a posteriori error estimates of the DG methods for the elliptic variational inequalities of the first kind and the second kind. This thesis contains seven chapters including an introductory chapter and a concluding chap-ter. In the introductory chapter, we review some fundamental preliminary results which will be used in the subsequent analysis. In Chapter 2, a posteriori error estimates for a class of DG meth-ods have been derived for the second order elliptic obstacle problem, which is a prototype for elliptic variational inequalities of the first kind. The analysis of Chapter 2 is carried out for the general obstacle function therefore the error estimator obtained therein involves the min/max func-tion and hence the computation of the error estimator becomes a bit complicated. With a mild assumption on the trace of the obstacle, we have derived a significantly simple and easily com-putable error estimator in Chapter 3. Numerical experiments illustrates that this error estimator indeed behaves better than the error estimator derived in Chapter 2. In Chapter 4, we have carried out a posteriori analysis of DG methods for the Signorini problem which arises from the study of the frictionless contact problems. A nonlinear smoothing map from the DG finite element space to conforming finite element space has been constructed and used extensively, in the analysis of Chapter 2, Chapter 3 and Chapter 4. Also, a common property shared by all DG methods allows us to carry out the analysis in unified setting. In Chapter 5, we study the C0 interior penalty method for the plate frictional contact problem, which is a fourth order variational inequality of the second kind. In this chapter, we have also established the medius analysis along with a posteriori analy-sis. Numerical results have been presented at the end of every chapter to illustrate the theoretical results derived in respective chapters. We discuss the possible extension and future proposal of the work presented in the Chapter 6. In the last chapter, we have documented the FEM codes used in the numerical experiments.