This thesis comprises the study of two moduli spaces of piecewise J-holomorphic curves. The main scheme is to consider a subdivision of the 2-sphere into a collection of small domains and to study collections of J-holomorphic maps into a symplectic manifold. These maps are coupled by Lagrangian boundary conditions. The work can be seen as finding a 2-dimensional analogue of the finite-dimensional path space approximation by piecewise geodesics on a Riemannian manifold (Q,g). For a nice class of target manifolds we consider tangent bundles of Riemannian manifolds and symplectizations of unit tangent bundles. Via polarization they provide a rich set of Lagrangians which can be used to define appropriate boundary value problems for the J-holomorphic pieces. The work focuses on existence theory as a pre-stage to global questions such as combinatorial refinement and the quality of the approximation. The first moduli space of lifted type is defined on a triangulation of the 2-sphere and consists of disks in the tangent bundle whose boundary projects onto geodesic triangles. The second moduli space of punctured type is defined on a circle packing domain and consists of boundary punctured disks in the symplectization of the unit tangent bundle. Their boundary components map into single fibers and at punctures the disks converge to geodesics. The coupling boundary conditions are chosen such that the piecewise problem always is Fredholm of index zero and both moduli spaces only depend on discrete data. For both spaces existence results are established for the J-holomorphic pieces which hold true on a small scale. Each proof employs a version of the implicit function theorem in a different setting. Here the argument for the moduli space of punctured type is more subtle. It rests on a connection to tropical geometry discovered by T. Ekholm for 1-jet spaces. The boundary punctured disks are constructed in the vicinity of explicit Morse flow trees which correspond to the limiting objects under degeneration of the boundary condition.

In der Finanzmathematik hat der Besitzer einer amerikanische Option das Recht aber nicht die Pflicht, eine Aktie innerhalb eines bestimmten Zeitraums, für einen bestimmten Preis zu kaufen oder zu verkaufen. Die Bewertung einer amerikanische Option wird als so genanntes optimale stopping Problem formuliert. Erfolgt die Modellierung des Aktienkurses durch eine geometrische Brownsche Bewegung, wird der Wert einer amerikanischen Option durch ein deterministisches freies Randwertproblem (FRWP), oder einer äquivalenten Variationsungleichung (VU) auf ganz R in gewichteten Sobolev Räumen gegeben. Um Standardmethoden der Numerischen Mathematik anzuwenden, wird das unbeschränkte Gebiet zu einem beschränkten Gebiet abgeschnitten. Mit Hilfe der Fourier-Transformation wird eine Integraldarstellung der Lösung die den freien Rand explizit beinhaltet, hergeleitet. Mittels dieser Integraldarstellung werden Abschneidefehlerschranken bewiesen. Danach werden gewichtete Poincare Ungleichungen mit expliziten Konstanten bewiesen. Der Abschneidefehler und die gewichtete Poincare Ungleichung ermöglichen, einen zuverlässigen a posteriori Fehlerschätzer zwischen der exakten Lösung der VU und der semidiskreten Lösung des penalisierten Problems auf R herzuleiten. Eine hinreichend glatte Lösung der VU garantiert die Konvergenz der Lösung des penaltisierten Problems zur Lösung der VU. Ein a priori Fehlerschätzer für den Fehler zwischen der exakten Lösung der VU und der semidiskreten Lösung des penaltisierten Problems beendet die numerische Analysis. Die eingeführten aposteriori Fehlerschätzer motivieren einen Algorithmus für adaptive Netzverfeinerung. Numerische Experimente zeigen die verbesserte Konvergenz des adaptiven Verfahrens gegenüber der uniformen Verfeinerung. Der zuverlässige a posteriori Fehlerschätzer ermöglicht es, den Abschneidepunkt so zu wählen, dass der Gesamtfehler (Diskretisierungsfehler plus Abschneidefehler) kleiner als eine gegebenen Toleranz ist. / Among the central concerns in mathematical finance is the evaluation of American options. An American option gives the holder the right but not the obligation to buy or sell a certain financial asset within a certain time-frame, for a certain strike price. The valuation of American options is formulated as an optimal stopping problem. If the stock price is modelled by a geometric Brownian motion, the value of an American option is given by a deterministic parabolic free boundary value problem (FBVP) or equivalently a non-symmetric variational inequality (VI) on weighted Sobolev spaces on R. To apply standard numerical methods, the unbounded domain R is truncated to a bounded one. Applying the Fourier transform to the FBVP yields an integral representation of the solution including the free boundary explicitely. This integral representation allows to prove explicit truncation errors. Since the VI is formulated within the framework of weighted Sobolev spaces, we establish a weighted Poincare inequality with explicit determined constants. The truncation error estimate and the weighted Poncare inequality enable a reliable a posteriori error estimate between the exact solution of the VI and the semi-discrete solution of the penalised problem on R. A sufficient regular solution provides the convergence of the solution of the penalised problem to the solution of the VI. An a priori error estimate for the error between the exact solution of the VI and the semi-discrete solution of the penalised problem concludes the numerical analysis. The established a posteriori error estimates motivates an algorithm for adaptive mesh refinement. Numerical experiments show the improved convergence of the adaptive algorithm compared to uniform mesh refinement. The reliable a posteriori error estimate including explicit truncation errors allows to determine a truncation point such that the total error (discretisation and truncation error) is below a given error tolerance.

Schwerpunkt dieser Arbeit ist die Untersuchung einer Klasse von Diskretisierungsmethoden für nichtlineare optimale Steuerungsprobleme mit gewöhnlichen Differentialgleichungen und Steuerungsbeschränkung sowie die Durchführung von numerischen Experimente. Die theoretischen Untersuchungen basieren aus einem gekoppeltes Parametrisierungs-Diskretisierungsschema aus stückweise polinomialen Ansatz für die Steuerung und einen Runge-Kutta-Verfahren zur Integration der Zustands-Differentialgleichung. Die Konvergenzordnung der Lösung wird unter Regularitätsbedingung und Optimalitätsbedingung 2.Ordnung ermittelt. Außerdem wird eine Möglichkeit zur numerischen Berechnung der Gradienten über internen numerischen Differentiation erläutert. Schließlich werden einige numerischen Resultate gegeben und die Abhängigkeiten zur den theoretischen Konvergenzresultate diskutiert. / The focal point of this work is the investigation of a class of discretization methods for nonlinear optimal control problems governed by ordinary differential equations with control restrictions, as well as the implementation of some numerical experiments. The theoretical investigations are based on a coupledparameterization-discretization pattern, a piecewise linear parameterization representation of the control, and the application of a Runge Kutta method for the integration of the differential state equation. The rate of convergence of the solution is obtained with the help of regularity conditions and the second order optimality conditions. Furthermore, we also present in this paper a possibility of the numerical computation of the gradients via numerical differentiation. Finally some numerical results are given and their relationship to the theoretical convergence results are discussed.

Die Dissertation behandelt methodische und algorithmische Aspekte der Spektralelementemethode zur räumlichen Diskretisierung partieller Differentialgleichungen. Die Weiterentwicklung einer symmetriebasierten Faktorisierung ermöglicht effiziente Operatoren für Tetraederelemente. Auf Grundlage einer umfassenden Leistungsanalyse werden Engpässe in der Implementierung der Operatoren identifiziert und durch algorithmische Modifikationen der Methode eliminiert.

Techn. University, Diss., 2004--Braunschweig.

Die Dissertation behandelt methodische und algorithmische Aspekte der Spektralelementemethode zur räumlichen Diskretisierung partieller Differentialgleichungen. Die Weiterentwicklung einer symmetriebasierten Faktorisierung ermöglicht effiziente Operatoren für Tetraederelemente. Auf Grundlage einer umfassenden Leistungsanalyse werden Engpässe in der Implementierung der Operatoren identifiziert und durch algorithmische Modifikationen der Methode eliminiert.

This thesis comprises the study of two moduli spaces of piecewise J-holomorphic curves. The main scheme is to consider a subdivision of the 2-sphere into a collection of small domains and to study collections of J-holomorphic maps into a symplectic manifold. These maps are coupled by Lagrangian boundary conditions. The work can be seen as finding a 2-dimensional analogue of the finite-dimensional path space approximation by piecewise geodesics on a Riemannian manifold (Q,g). For a nice class of target manifolds we consider tangent bundles of Riemannian manifolds and symplectizations of unit tangent bundles. Via polarization they provide a rich set of Lagrangians which can be used to define appropriate boundary value problems for the J-holomorphic pieces. The work focuses on existence theory as a pre-stage to global questions such as combinatorial refinement and the quality of the approximation. The first moduli space of lifted type is defined on a triangulation of the 2-sphere and consists of disks in the tangent bundle whose boundary projects onto geodesic triangles. The second moduli space of punctured type is defined on a circle packing domain and consists of boundary punctured disks in the symplectization of the unit tangent bundle. Their boundary components map into single fibers and at punctures the disks converge to geodesics. The coupling boundary conditions are chosen such that the piecewise problem always is Fredholm of index zero and both moduli spaces only depend on discrete data. For both spaces existence results are established for the J-holomorphic pieces which hold true on a small scale. Each proof employs a version of the implicit function theorem in a different setting. Here the argument for the moduli space of punctured type is more subtle. It rests on a connection to tropical geometry discovered by T. Ekholm for 1-jet spaces. The boundary punctured disks are constructed in the vicinity of explicit Morse flow trees which correspond to the limiting objects under degeneration of the boundary condition.

Diese Arbeit beschäftigt sich mit verschiedenen mathematischen Fragestellungen hinsichtlich der Modellierung, Analysis und numerischen Simulation von Gasnetzen. Hierbei liegt der Fokus auf der mathematischen Handhabung von partiellen differential-algebraischen Gleichungen, die mit algebraischen Gleichungen gekoppelt sind. Diese bieten einen einfachen Zugang hinsichtlich der Modellierung von dynamischen Strukturen auf Netzen Somit sind sie insbesondere für Gasnetze geeignet, denen im Zuge der steigenden Bedeutung von erneuerbaren Energien ein gestiegenes Interesse seitens der Öffentlichkeit, Politik und Wissenschaft entgegen gebracht wird. Wir führen zunächst die gängigsten Elemente, die in Gasnetzen benötigt werden ein und formulieren zwei PDAE-Klassen für solche Netze: Eine für reine Rohrnetze, und eine, die zusätzliche Elemente wie Verdichter und Widerstände beinhaltet. Des Weiteren untersuchen wir die Sensitivität der Lösung der Rohrnetz-PDAE hinsichtlich Störungen. Dabei berücksichtigen wir Störungen, die nicht nur den dynamischen Teil der PDAE beeinflussen, sondern auch Störungen in den algebraischen Gleichungen und weisen Stabilitätseigenschaften für die Lösung der PDAE nach. Darüber hinaus beschäftigen wir uns mit einer neu entwickelten, an die Netztopologie angepassten Ortsdiskretisierung, welche die Stabilitätseigenschaften der PDAE auf DAE Systeme überträgt. Des Weiteren zeigen wir, wie sich die Gasnetz-DAE zu einer gewöhnlichen Differentialgleichung, welche die inhärente Dynamik der DAE widerspiegelt entkoppeln lässt. Dieses entkoppelte System kann darüber hinaus direkt aus den Topologie- und Elementinformationen des Netzes aufgestellt werden. Abschließend demonstrieren wir die Ergebnisse an Benchmark-Gasnetzen. Dabei vergleichen wir sowohl die entkoppelte Differentialgleichung mit dem ursprünglichen DAE System, zeigen aber auch, welche Vorteile die an die Netztopologie angepasste Ortsdiskretisierung gegenüber existierenden Verfahren besitzt. / This thesis addresses several aspects regarding modelling, analysis and numerical simulation of gas networks. Hereby, our focus lies on (partial) differential-algebraic equations, thus systems of partial and ordinary differential equations which are coupled by algebraic equations. These coupled systems allow an easy approach towards the modelling of dynamic structures on networks. Therefore, they are well suited for gas networks, which have gained a rise of attention in society, politics and science due to the focus towards renewable energies. We give an introduction towards gas network modelling that includes the most common elements that also appear in real gas networks and present two PDAE systems: One for pipe networks and one that includes additional elements like resistors and compressors. Furthermore, we investigate the impact of perturbations onto the pipe network PDAE, where we explicitly allow perturbations to affect the system in the differential as well as in the algebraic components. We conclude that the solution of the PDAE possesses stability properties. In addition, this thesis introduces a new spatial discretisation that is adapted to the net- work topology. This topology-adapted semi-discretisation results in a DAE which possesses the same perturbation behaviour as the space continuous PDAE. Furthermore, we present a topology based decoupling procedure that allows to reformulate the DAE as an ordinary differential equation (ODE), which represents the inherent dynamics of the DAE system. This ODE, together with a decoupled set of algebraic equations, can be derived from the topology and element information directly. We conclude by demonstrating the established results for several benchmark networks. This includes a comparison of numerical solutions for the decoupled ODE and the DAE system. In addition we present the advantages of the topology-adapted spatial discretisation over existing well established methods.

This thesis deals with regularization parameter selection methods in the context of Tikhonov-type regularization with Poisson distributed data, in particular the reconstruction of images, as well as with the identification of the volatility surface from observed option prices. In Part I we examine the choice of the regularization parameter when reconstructing an image, which is disturbed by Poisson noise, with Tikhonov-type regularization. This type of regularization is a generalization of the classical Tikhonov regularization in the Banach space setting and often called variational regularization. After a general consideration of Tikhonov-type regularization for data corrupted by Poisson noise, we examine the methods for choosing the regularization parameter numerically on the basis of two test images and real PET data. In Part II we consider the estimation of the volatility function from observed call option prices with the explicit formula which has been derived by Dupire using the Black-Scholes partial differential equation. The option prices are only available as discrete noisy observations so that the main difficulty is the ill-posedness of the numerical differentiation. Finite difference schemes, as regularization by discretization of the inverse and ill-posed problem, do not overcome these difficulties when they are used to evaluate the partial derivatives. Therefore we construct an alternative algorithm based on the weak formulation of the dual Black-Scholes partial differential equation and evaluate the performance of the finite difference schemes and the new algorithm for synthetic and real option prices.

This thesis deals with regularization parameter selection methods in the context of Tikhonov-type regularization with Poisson distributed data, in particular the reconstruction of images, as well as with the identification of the volatility surface from observed option prices. In Part I we examine the choice of the regularization parameter when reconstructing an image, which is disturbed by Poisson noise, with Tikhonov-type regularization. This type of regularization is a generalization of the classical Tikhonov regularization in the Banach space setting and often called variational regularization. After a general consideration of Tikhonov-type regularization for data corrupted by Poisson noise, we examine the methods for choosing the regularization parameter numerically on the basis of two test images and real PET data. In Part II we consider the estimation of the volatility function from observed call option prices with the explicit formula which has been derived by Dupire using the Black-Scholes partial differential equation. The option prices are only available as discrete noisy observations so that the main difficulty is the ill-posedness of the numerical differentiation. Finite difference schemes, as regularization by discretization of the inverse and ill-posed problem, do not overcome these difficulties when they are used to evaluate the partial derivatives. Therefore we construct an alternative algorithm based on the weak formulation of the dual Black-Scholes partial differential equation and evaluate the performance of the finite difference schemes and the new algorithm for synthetic and real option prices.