Sie haben Javascript deaktiviert!
Sie haben versucht eine Funktion zu nutzen, die nur mit Javascript möglich ist. Um sämtliche Funktionalitäten unserer Internetseite zu nutzen, aktivieren Sie bitte Javascript in Ihrem Browser.

Schnee auf dem Campus. Bildinformationen anzeigen

Schnee auf dem Campus.

Foto: Universität Paderborn, Johannes Pauly

Christian Offen, PhD, Massey University, New Zealand

 Christian Offen, PhD, Massey University, New Zealand

Angewandte Mathematik - Numerik und Steuerung


Technologiepark 21
33100 Paderborn
 Christian Offen, PhD, Massey University, New Zealand
08/2020 - 08/2023

Wissenschaftlicher Mitarbeiter

Universität Paderborn, Institut für Mathematik. Arbeitsgruppe Angewandte Mathematik, Numerik und Optimierung.

11/2016 - 09/2020


Massey University, Neuseeland. Dissertation: Analysis of Hamiltonian boundary value problems and symplectic integration, Betreuer: Dist. Prof. Robert McLachlan. (2019 - Hatherton Award,   2020 - Dean's List of Exceptional Theses)

10/2014 - 09/2016

Master of Science Mathematik

Universität Hamburg, Masterstudium. Masterarbeit: Seiberg-Witten Curves and Parabolic Affine Spheres (Differential Geometry), Betreuer: Prof. Dr. Vicente Cortés

10/2011 - 09/2014

Bachelor of Science Mathematik

Universität Hamburg, Bachelor Studium. Bachelorarbeit: Die Menge der Lösungen der Monge-Ampère-Gleichung in der endlich punktierten reellen Ebene (Differentialgeometrie), Betreuer: Prof. Dr. Vicente Cortés

Liste im Research Information System öffnen


Backward error analysis for variational discretisations of partial differential equations

R.I. McLachlan, C. Offen, in: arXiv:2006.14172, 2020

In backward error analysis, an approximate solution to an equation is compared to the exact solution to a nearby "modified" equation. In numerical ordinary differential equations, the two agree up to any power of the step size. If the differential equation has a geometric property then the modified equation may share it. In this way, known properties of differential equations can be applied to the approximation. But for partial differential equations, the known modified equations are of higher order, limiting applicability of the theory. Therefore, we study symmetric solutions of discretized partial differential equations that arise from a discrete variational principle. These symmetric solutions obey infinite-dimensional functional equations. We show that these equations admit second-order modified equations which are Hamiltonian and also possess first-order Lagrangians in modified coordinates. The modified equation and its associated structures are computed explicitly for the case of rotating travelling waves in the nonlinear wave equation.

Detection of high codimensional bifurcations in variational PDEs

L.M. Kreusser, R.I. McLachlan, C. Offen, Nonlinearity (2020), 33(5), pp. 2335-2363

Analysis of Hamiltonian boundary value problems and symplectic integration

C. Offen, Massey University, 2020

Ordinary differential equations (ODEs) and partial differential equations (PDEs) arise in most scientific disciplines that make use of mathematical techniques. As exact solutions are in general not computable, numerical methods are used to obtain approximate solutions. In order to draw valid conclusions from numerical computations, it is crucial to understand which qualitative aspects numerical solutions have in common with the exact solution. Symplecticity is a subtle notion that is related to a rich family of geometric properties of Hamiltonian systems. While the effects of preserving symplecticity under discretisation on long-term behaviour of motions is classically well known, in this thesis (a) the role of symplecticity for the bifurcation behaviour of solutions to Hamiltonian boundary value problems is explained. In parameter dependent systems at a bifurcation point the solution set to a boundary value problem changes qualitatively. Bifurcation problems are systematically translated into the framework of classical catastrophe theory. It is proved that existing classification results in catastrophe theory apply to persistent bifurcations of Hamiltonian boundary value problems. Further results for symmetric settings are derived. (b) It is proved that to preserve generic bifurcations under discretisation it is necessary and sufficient to preserve the symplectic structure of the problem. (c) The catastrophe theory framework for Hamiltonian ODEs is extended to PDEs with variational structure. Recognition equations for A-series singularities for functionals on Banach spaces are derived and used in a numerical example to locate high-codimensional bifurcations. (d) The potential of symplectic integration for infinite-dimensional Lie-Poisson systems (Burgers’ equation, KdV, fluid equations, . . . ) using Clebsch variables is analysed. It is shown that the advantages of symplectic integration can outweigh the disadvantages of integrating over a larger phase space introduced by a Clebsch representation. (e) Finally, the preservation of variational structure of symmetric solutions in multisymplectic PDEs by multisymplectic integrators on the example of (phase-rotating) travelling waves in the nonlinear wave equation is discussed.

Preservation of Bifurcations of Hamiltonian Boundary Value Problems Under Discretisation

R.I. McLachlan, C. Offen, Foundations of Computational Mathematics (2020), 20(6), pp. 1363-1400

We show that symplectic integrators preserve bifurcations of Hamiltonian boundary value problems and that nonsymplectic integrators do not. We provide a universal description of the breaking of umbilic bifurcations by nonysmplectic integrators. We discover extra structure induced from certain types of boundary value problems, including classical Dirichlet problems, that is useful to locate bifurcations. Geodesics connecting two points are an example of a Hamiltonian boundary value problem, and we introduce the jet-RATTLE method, a symplectic integrator that easily computes geodesics and their bifurcations. Finally, we study the periodic pitchfork bifurcation, a codimension-1 bifurcation arising in integrable Hamiltonian systems. It is not preserved by either symplectic on nonsymplectic integrators, but in some circumstances symplecticity greatly reduces the error.


    Symplectic integration of PDEs using Clebsch variables

    R.I. McLachlan, C. Offen, B.K. Tapley, Journal of Computational Dynamics (2019), 6(1), pp. 111-130

    Many PDEs (Burgers' equation, KdV, Camassa-Holm, Euler's fluid equations, …) can be formulated as infinite-dimensional Lie-Poisson systems. These are Hamiltonian systems on manifolds equipped with Poisson brackets. The Poisson structure is connected to conservation properties and other geometric features of solutions to the PDE and, therefore, of great interest for numerical integration. For the example of Burgers' equations and related PDEs we use Clebsch variables to lift the original system to a collective Hamiltonian system on a symplectic manifold whose structure is related to the original Lie-Poisson structure. On the collective Hamiltonian system a symplectic integrator can be applied. Our numerical examples show excellent conservation properties and indicate that the disadvantage of an increased phase-space dimension can be outweighed by the advantage of symplectic integration.


    Bifurcation of solutions to Hamiltonian boundary value problems

    R.I. McLachlan, C. Offen, Nonlinearity (2018), pp. 2895-2927

    A bifurcation is a qualitative change in a family of solutions to an equation produced by varying parameters. In contrast to the local bifurcations of dynamical systems that are often related to a change in the number or stability of equilibria, bifurcations of boundary value problems are global in nature and may not be related to any obvious change in dynamical behaviour. Catastrophe theory is a well-developed framework which studies the bifurcations of critical points of functions. In this paper we study the bifurcations of solutions of boundary-value problems for symplectic maps, using the language of (finite-dimensional) singularity theory. We associate certain such problems with a geometric picture involving the intersection of Lagrangian submanifolds, and hence with the critical points of a suitable generating function. Within this framework, we then study the effect of three special cases: (i) some common boundary conditions, such as Dirichlet boundary conditions for second-order systems, restrict the possible types of bifurcations (for example, in generic planar systems only the A-series beginning with folds and cusps can occur); (ii) integrable systems, such as planar Hamiltonian systems, can exhibit a novel periodic pitchfork bifurcation; and (iii) systems with Hamiltonian symmetries or reversing symmetries can exhibit restricted bifurcations associated with the symmetry. This approach offers an alternative to the analysis of critical points in function spaces, typically used in the study of bifurcation of variational problems, and opens the way to the detection of more exotic bifurcations than the simple folds and cusps that are often found in examples.

      Hamiltonian boundary value problems, conformal symplectic symmetries, and conjugate loci

      R.I. McLachlan, C. Offen, New Zealand Journal of Mathematics (2018), 48, pp. 83-99

      In this paper we continue our study of bifurcations of solutions of boundary-value problems for symplectic maps arising as Hamiltonian diffeomorphisms. These have been shown to be connected to catastrophe theory via generating functions and ordinary and reversal phase space symmetries have been considered. Here we present a convenient, coordinate free framework to analyse separated Lagrangian boundary value problems which include classical Dirichlet, Neumann and Robin boundary value problems. The framework is then used to prove the existence of obstructions arising from conformal symplectic symmetries on the bifurcation behaviour of solutions to Hamiltonian boundary value problems. Under non-degeneracy conditions, a group action by conformal symplectic symmetries has the effect that the flow map cannot degenerate in a direction which is tangential to the action. This imposes restrictions on which singularities can occur in boundary value problems. Our results generalise classical results about conjugate loci on Riemannian manifolds to a large class of Hamiltonian boundary value problems with, for example, scaling symmetries.

      Symplectic integration of boundary value problems

      R.I. McLachlan, C. Offen, Numerical Algorithms (2018), pp. 1219-1233

      Symplectic integrators can be excellent for Hamiltonian initial value problems. Reasons for this include their preservation of invariant sets like tori, good energy behaviour, nonexistence of attractors, and good behaviour of statistical properties. These all refer to {\em long-time} behaviour. They are directly connected to the dynamical behaviour of symplectic maps φ:M→M' on the phase space under iteration. Boundary value problems, in contrast, are posed for fixed (and often quite short) times. Symplecticity manifests as a symplectic map φ:M→M' which is not iterated. Is there any point, therefore, for a symplectic integrator to be used on a Hamiltonian boundary value problem? In this paper we announce results that symplectic integrators preserve bifurcations of Hamiltonian boundary value problems and that nonsymplectic integrators do not.

        Local intersections of Lagrangian manifolds correspond to catastrophe theory

        C. Offen, in: arXiv:1811.10165, 2018

        Two smooth map germs are right-equivalent if and only if they generate two Lagrangian submanifolds in a cotangent bundle which have the same contact with the zero-section. In this paper we provide a reverse direction to this classical result of Golubitsky and Guillemin. Two Lagrangian submanifolds of a symplectic manifold have the same contact with a third Lagrangian submanifold if and only if the intersection problems correspond to stably right equivalent map germs. We, therefore, obtain a correspondence between local Lagrangian intersection problems and catastrophe theory while the classical version only captures tangential intersections. The correspondence is defined independently of any Lagrangian fibration of the ambient symplectic manifold, in contrast to other classical results. Moreover, we provide an extension of the correspondence to families of local Lagrangian intersection problems. This gives rise to a framework which allows a natural transportation of the notions of catastrophe theory such as stability, unfolding and (uni-)versality to the geometric setting such that we obtain a classification of families of local Lagrangian intersection problems. An application is the classification of Lagrangian boundary value problems for symplectic maps.

        Liste im Research Information System öffnen

        Die Universität der Informationsgesellschaft