Zusammenfassung Diese Arbeit ist ein Versuch, verschiedene klassische und neuere Methodender glatten bzw. nichtglatten Optimierung zu verallgemeinern und in ihrem Zusammenhang darzustellen. Als Hauptinstrument erweist sich dabei die sogenannte verallgemeinerte Kojima-Funktion. Neben reichlichen Beispielen setzen wir einen besonderen Akzent auf die Betrachtung von Variationsungleichungen, Komplementaritaetsaufgaben und der Standartaufgabeder mathematischen Programmierung. Unter natuerlichen Voraussetzungen an diese Probleme kann man u.a. Barriere-, Straf- und SQP-Typ-Methoden, die auf Newton-Verfahrenbasieren, aber auch Modelle, die sogenannte NCP-Funktionen benutzen, mittelsspezieller Stoerungen der Kojima-Funktion exakt modellieren. Daneben werdendurch explizite und natuerliche Wahl der Stoerungsparameter auch neue Methoden dieser Arten vorgeschlagen. Die Vorteile solcher Modellierungsind ueberzeugend vor allem wegen der direkt moeglichen (auf Stabilitaetseigenschaften der Kojima-Gleichung beruhendenden)Loesungsabschaetzungen und weil die entsprechenden Nullstellen ziemlich einfach als Loesungen bekannter Ersatzprobleme interpretiert werden koennen. Ein weiterer Aspekt der Arbeit besteht in der genaueren Untersuchungder "nichtglatten Faelle". Hier wird die Theorie von verschiedenen verallgemeinerten Ableitungen und dadurch entstehenden verallgemeinerten Newton-Verfahren, die im Buch "Nonsmooth Equations in Optimization" von B. Kummer und D. Klatte vorgeschlagen und untersucht wurde, intensiv benutzt. Entscheidend ist dabei, dass die benutzten verallgemeinerten Ableitungen auch praktisch angewandt werden koennen, da man sie exakt ausrechnen kann. / This work attempts to generalize various classical and new methods of smooth or nonsmooth optimization and to show them in their interrelation. The main tool for doing this is the so-called generalized Kojima-function. In addition to numerous examples we specialy emphasize the consideration of variational inequalities, complementarity problems and the standard problem of mathematical programming. Under natural assumptions on these problems we can model e.g. barrier-, penalty-, and SQP-Type-methods basing on Newton methods, and also methods using the so-called NCP-function exactly by means of special perturbations of the Kojima-function. Furthermore, by the explicit and natural choice of the perturbation parameters new methods of these kinds are introduced. The benefit of such a modelling is obvious, first of all due to the direct solution estimation (basing on stability properties of the Kojima-equation) and because the corresponding zeros can easily be interpreted as solutions of known subproblems. A further aspect considered in this paper is the detailed investigation of "nonsmooth cases". The theory of various generalized derivatives and resulting generalized Newton methods, which is introduced and investigated in the book "Nonsmooth Equations in Optimization" of B. Kummer and D. Klatte, is intensely used here. The crucial point is the applicability of the used generalized derivatives in practice, since they can be calculated exactly.

Diese Arbeit behandelt garantierte Fehlerkontrolle für elliptische partielle Differentialgleichungen anhand des Poisson-Modellproblems, des Stokes-Problems und des Hindernisproblems. Hierzu werden garantierte obere Schranken für den Energiefehler zwischen exakter Lösung und diskreten Finite-Elemente-Approximationen erster Ordnung entwickelt. Ein verallgemeinerter Ansatz drückt den Energiefehler durch Dualnormen eines oder mehrerer Residuen aus. Hinzu kommen berechenbare Zusatzterme, wie Oszillationen der gegebenen Daten, mit expliziten Konstanten. Für die Abschätzung der Dualnormen der Residuen existieren viele verschiedene Techniken. Diese Arbeit beschäftigt sich vorrangig mit Equilibrierungsschätzern, basierend auf Raviart-Thomas-Elementen, welche effiziente garantierte obere Schranken ermöglichen. Diese Schätzer werden mit einem Postprocessing-Verfahren kombiniert, das deren Effizienz mit geringem zusätzlichen Rechenaufwand deutlich verbessert. Nichtkonforme Finite-Elemente-Methoden erzeugen zusätzlich ein Inkonsistenzresiduum, dessen Dualnorm mit Hilfe diverser konformer Approximationen abgeschätzt wird. Ein Nebenaspekt der Arbeit betrifft den expliziten residuen-basierten Fehlerschätzer, der für gewöhnlich optimale und leicht zu berechnende Verfeinerungsindikatoren für das adaptive Netzdesign liefert, aber nur schlechte garantierte obere Schranken. Eine neue Variante, die auf den equilibrierten Flüssen des Luce-Wohlmuth-Fehlerschätzers basiert, führt zu stark verbesserten Zuverlässigkeitskonstanten. Eine Vielzahl numerischer Experimente vergleicht alle implementierten Fehlerschätzer und zeigt, dass effiziente und garantierte Fehlerkontrolle in allen vorliegenden Modellproblemen möglich ist. Insbesondere zeigt ein Modellproblem, wie die Fehlerschätzer erweitert werden können, um auch auf Gebieten mit gekrümmten Rändern garantierte obere Schranken zu liefern. / This thesis studies guaranteed error control for elliptic partial differential equations on the basis of the Poisson model problem, the Stokes equations and the obstacle problem. The error control derives guaranteed upper bounds for the energy error between the exact solution and different finite element discretisations, namely conforming and nonconforming first-order approximations. The unified approach expresses the energy error by dual norms of one or more residuals plus computable extra terms, such as oscillations of the given data, with explicit constants. There exist various techniques for the estimation of the dual norms of such residuals. This thesis focuses on equilibration error estimators based on Raviart-Thomas finite elements, which permit efficient guaranteed upper bounds. The proposed postprocessing in this thesis considerably increases their efficiency at almost no additional computational costs. Nonconforming finite element methods also give rise to a nonconsistency residual that permits alternative treatment by conforming interpolations. A side aspect concerns the explicit residual-based error estimator that usually yields cheap and optimal refinement indicators for adaptive mesh refinement but not very sharp guaranteed upper bounds. A novel variant of the residual-based error estimator, based on the Luce-Wohlmuth equilibration design, leads to highly improved reliability constants. A large number of numerical experiments compares all implemented error estimators and provides evidence that efficient and guaranteed error control in the energy norm is indeed possible in all model problems under consideration. Particularly, one model problem demonstrates how to extend the error estimators for guaranteed error control on domains with curved boundary.

Nous présentons une contribution autour de la modélisation des écoulements viscoplastiques. En vue d'applications réalistes telle que la simulation numérique des coulées de lave volcanique, le travail se concentre particulièrement sur les fluides complexes dont la rhéologie dépend fortement de grandeurs physiques telle que la température ou la concentration en particule. Nous développons un nouvel algorithme de résolution numérique des équations de Herschel-Bulkley combinant une méthode de Lagrangien augmenté à paramètre d'augmentation variable, une méthode des caractéristiques d'ordre 2 et une adaptation de maillage automatique. Sur des problèmes stationnaires ou en évolution tel que le problème test de la cavité entraînée, il apporte une solution efficace pour garantir à la fois une précision numérique élevée et un temps de calcul raisonnable. Cet algorithme est ensuite étendue et adapté au cas des rhéologies non-isothermes et aux suspensions. Concernant la simulation numérique des coulées de lave volcanique, nous détaillons une méthode de réduction par analyse asymptotique des équations de Herschel-Bulkley pour des écoulements de faible épaisseur sur une topographie arbitraire. Elle permet alors de décrire ces écoulements tridimensionnels de fluides viscoplastiques à surface libre par des équations bidimensionnelles surfaciques. Cette approche est ensuite étendue au cas non-isotherme en y ajoutant l'équation de la chaleur et des dépendances thermiques sur la rhéologie. Par intégration verticale de l'équation de la chaleur, on retrouve un modèle bidimensionnel. Le modèle non-isotherme est validé sur une expérience de dôme réalisée en laboratoire et une simulation numérique est réalisée autour d'une coulée qui a eu lieu sur le volcan du Piton de la Fournaise à la Réunion, en décembre 2010. La comparaison donne des résultats qui sont de notre point de vue satisfaisants et encourageants. / We present a contribution about modeling of viscoplastic flows. For realistic applications such as numerical simulation of volcanic lava flows, the work focuses particularly on complex fluids whose rheology strongly depends on physical quantities such as temperature or the particle concentration. We develop a new numerical resolution algorithm of Herschel-Bulkley's equations combining an augmented Lagrangian method with variable augmentation parameter, a second order characteristic method and an auto-adaptive mesh procedure. On stationary or evolving problems as the lid-driven cavity flow benchmark, it provides an effective solution to ensure both a high numerical accuracy within a reasonable computing time. This algorithm is then extended and adapted to the case of non-isothermal rheological and suspensions. On the numerical simulation of volcanic lava flows, we describe a method of reducing by asymptotic analysis of the Herschel-Bulkley's equations for thin flows on arbitrary topography. It allows to describe the three-dimensional flows of viscoplastic fluid with free surface by bidimensional surface equations. This approach is then extended to the non-isothermal case by adding the heat equation and thermal dependencies on rheology. By vertical integration of the heat equation, a two-dimensional model is maintained . The non-isothermal model is validated on a laboratory experiment of dome and a numerical simulation is performed on a December 2010 Piton de la Fournaise lava flow from La Réunion island. In our view, the comparison gives satisfactory and encouraging results.

Nous présentons une contribution autour de la modélisation des écoulements viscoplastiques. En vue d'applications réalistes telle que la simulation numérique des coulées de lave volcanique, le travail se concentre particulièrement sur les fluides complexes dont la rhéologie dépend fortement de grandeurs physiques telle que la température ou la concentration en particule. Nous développons un nouvel algorithme de résolution numérique des équations de Herschel-Bulkley combinant une méthode de Lagrangien augmenté à paramètre d'augmentation variable, une méthode des caractéristiques d'ordre 2 et une adaptation de maillage automatique. Sur des problèmes stationnaires ou en évolution tel que le problème test de la cavité entraînée, il apporte une solution efficace pour garantir à la fois une précision numérique élevée et un temps de calcul raisonnable. Cet algorithme est ensuite étendue et adapté au cas des rhéologies non-isothermes et aux suspensions. Concernant la simulation numérique des coulées de lave volcanique, nous détaillons une méthode de réduction par analyse asymptotique des équations de Herschel-Bulkley pour des écoulements de faible épaisseur sur une topographie arbitraire. Elle permet alors de décrire ces écoulements tridimensionnels de fluides viscoplastiques à surface libre par des équations bidimensionnelles surfaciques. Cette approche est ensuite étendue au cas non-isotherme en y ajoutant l'équation de la chaleur et des dépendances thermiques sur la rhéologie. Par intégration verticale de l'équation de la chaleur, on retrouve un modèle bidimensionnel. Le modèle non-isotherme est validé sur une expérience de dôme réalisée en laboratoire et une simulation numérique est réalisée autour d'une coulée qui a eu lieu sur le volcan du Piton de la Fournaise à la Réunion, en décembre 2010. La comparaison donne des résultats qui sont de notre point de vue satisfaisants et encourageants. / We present a contribution about modeling of viscoplastic flows. For realistic applications such as numerical simulation of volcanic lava flows, the work focuses particularly on complex fluids whose rheology strongly depends on physical quantities such as temperature or the particle concentration. We develop a new numerical resolution algorithm of Herschel-Bulkley's equations combining an augmented Lagrangian method with variable augmentation parameter, a second order characteristic method and an auto-adaptive mesh procedure. On stationary or evolving problems as the lid-driven cavity flow benchmark, it provides an effective solution to ensure both a high numerical accuracy within a reasonable computing time. This algorithm is then extended and adapted to the case of non-isothermal rheological and suspensions. On the numerical simulation of volcanic lava flows, we describe a method of reducing by asymptotic analysis of the Herschel-Bulkley's equations for thin flows on arbitrary topography. It allows to describe the three-dimensional flows of viscoplastic fluid with free surface by bidimensional surface equations. This approach is then extended to the non-isothermal case by adding the heat equation and thermal dependencies on rheology. By vertical integration of the heat equation, a two-dimensional model is maintained . The non-isothermal model is validated on a laboratory experiment of dome and a numerical simulation is performed on a December 2010 Piton de la Fournaise lava flow from La Réunion island. In our view, the comparison gives satisfactory and encouraging results.

Tikhonov-type regularization of linear and nonlinear ill-posed problems in abstract spaces under sparsity constraints gained relevant attention in the past years. Since under some weak assumptions all regularized solutions are sparse if the l1-norm is used as penalty term, the l1-regularization was studied by numerous authors although the non-reflexivity of the Banach space l1 and the fact that such penalty functional is not strictly convex lead to serious difficulties. We consider the case that the sparsity assumption is narrowly missed. This means that the solutions may have an infinite number of nonzero but fast decaying components. For that case we formulate and prove convergence rates results for the l1-regularization of nonlinear operator equations. In this context, we outline the situations of Hölder rates and of an exponential decay of the solution components.