TY - GEN A1 - Dellnitz, Michael A1 - Hohmann, Andreas T1 - A Subdivision Algorithm for the Computation of Unstable Manifolds and Global Attractors N2 - Each invariant set of a given dynamical system is part of the global attractor. Therefore the global attractor contains all the potentially interesting dynamics, and, in particular, it contains every (global) unstable manifold. For this reason it is of interest to have an algorithm which allows to approximate the global attractor numerically. In this article we develop such an algorithm using a subdivision technique. We prove convergence of this method in a very general setting, and, moreover, we describe the qualitative convergence behavior in the presence of a hyperbolic structure. The algorithm can successfully be applied to dynamical systems of moderate dimension, and we illustrate this fact by several numerical examples. T3 - ZIB-Report - SC-95-11 Y1 - 1995 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1770 ER - TY - GEN A1 - Dellnitz, Michael A1 - Golubitsky, Martin A1 - Hohmann, Andreas A1 - Stewart, Ian T1 - Spirals in Scalar Reaction-Diffusion Equations N2 - Spiral-like patterns are an often observed phenomenon in chemical experiments such as the Belousov-Zhabotinskii reaction. The talk is concerned with a new PDE model whose solutions have the form of rotating spirals. In contrast to previous approaches it is based on a {\em scalar\/} reaction diffusion equation defined on a disk. A particular choice of boundary conditions leads to a non-selfadjoint operator which permits non-trivial dynamics. We study this equation using a combination of equivariant bifurcation theory and numerical simulations. The latter involves the direct simulation of the time dependent system as well as the computation of rotating waves and their stability. T3 - ZIB-Report - SC-95-06 Y1 - 1995 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1720 ER - TY - GEN A1 - Hohmann, Andreas T1 - Multilevel Newton h-p Collocation N2 - We derive a simple accuracy matching strategy for inexact Gauss Newton methods and apply it to the numerical solution of boundary value problems of ordinary differential equations by collocation. The matching strategy is based on an affine contravariant convergence theorem, i. e. , the characteristic constants are invariant under affine transformations of the domain. The inexact Gauss Newton method is applied to an integral formulation of the BVP. As discretization for the arising linear subproblems we employ adaptive collocation at Gaussian nodes with varying local orders and stepsizes. The grids are chosen via adaptive refinement and order selection. T3 - ZIB-Report - SC-94-25 Y1 - 1994 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1530 ER - TY - GEN A1 - Wulff, Claudia A1 - Hohmann, Andreas A1 - Deuflhard, Peter T1 - Numerical Continuation of Periodic Orbits with Symmetry. N2 - We consider periodic orbits of autonomous parameter dependent ODE's. Starting from a shooting algorithm for the numerical computation of periodic orbits via an adaptive Poincar\'e-section we develop a pathfollowing algorithm for periodic solutions based on a tangential continuation method with implicit reparametrization. For ODE's equivariant w.r.t. a finite group we show that spatial as well as spatio-temporal symmetries of periodic orbits can be exploited within the (multiple) shooting context. We describe how turning points, period doubling bifurcations and Hopf points along the branch of periodic solutions can be handled. Furthermore equivariant Hopf points and generic secondary bifurcations of periodic orbits with $ Z_m$-symmetry are treated. We tested the code with standard examples, e.g., the period doubling cascade in the Lorenz equations. To show the efficiency of the described methods we also used the program for an application from electronics, a ring oscillator with $n $ inverters. In this example the exploitation of symmetry reduces the amount of work for the continuation of periodic orbits from ${\cal O}(n^2)$ to ${\cal O}(n)$ T3 - ZIB-Report - SC-94-12 Y1 - 1994 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1421 ER - TY - GEN A1 - Hohmann, Andreas T1 - Object Oriented Design of Multilevel Newton and Continuation Methods. N2 - We present a set of C++ classes that realize an abstract inexact Gauss Newton method in combination with a continuation process for the solution of parameter dependent nonlinear problems. The object oriented approach allows the continuation of different types of solutions within the same framework. \\{\bf Abstract: }We present a set of C++ classes that realize an abstract inexact Gauss Newton method in combination with a continuation process for the solution of parameter dependent nonlinear problems. The object oriented approach allows the continuation of different types of solutions within the same framework. T3 - ZIB-Report - SC-94-04 Y1 - 1994 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1342 ER - TY - GEN A1 - Schütte, Christof A1 - Hohmann, Andreas A1 - Dinand, Manfred T1 - Numerical Simulation of Relaxation Oscillations of Waveguide-Lasers. N2 - An analysis of relaxation oscillations in local Er-doped optically pumped lasers is reported. It is based on a time dependent rate equation model for a quasi-two-level-system with wavelength dependent emission- and absorption cross-sections. For the first time a numerically reliable simulation of the characteristic laser behaviour was possible: the onset and decay of the oscillations, the time-dependent repetition period and the steady state signal output power. The characteristic waveguide parameters, as the erbium-concentration profile, the polarization dependent pump- and signal mode intensity profiles, the scattering losses, the cavity length and the front and rear reflectivities were all taken into account. The basic formulas are general and can also be used for Er-doped fiber lasers. Mathematically the problem can be characterized as a large boundary value problem, which can approximately be replaced by a stiff initial value problem of ordinary differential equations. The used algorithmic replacement procedure is motivated and discussed in detail. Here, pump- and signal evolution versus time are presented for an planar Er-diffused $\rm Ti$:$\rm LiNbO_{3}$ waveguide laser. The numerically obtained results show a nearly quantitative agreement with experimental investigations. Simultanously they supply knowledge about non-measureable (space-dependent population dynamic of the Er-atoms) and till today not measured data (dynamical response of the laser by a sharp peak in the external pump). T3 - ZIB-Report - SC-93-14 Y1 - 1993 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1106 ER - TY - GEN A1 - Hohmann, Andreas T1 - An Adaptive Continuation Method for Implicitly Defined Surfaces. N2 - A new method for the numerical aproximation of an implicitly defined surface is presented. It is a generalization of the Euler- Gauss-Newton method for implicitly defined (one- parameter) curves to the case of (two-parameter) surfaces. The basic task in the more general case is an efficient combination of modern CAGD techniques (such as triangular Bernstein-Bzier patches and the nine parameter Hermite interpolant) and the rank deficient Gauss-Newton method. T3 - ZIB-Report - SC-91-20 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-705 ER - TY - GEN A1 - Gatermann, Karin A1 - Hohmann, Andreas T1 - Hexagonal Lattice Dome - Illustration of a Nontrivial Bifurcation Problem. N2 - The deformation of a hexagonal lattice dome under an external load is an example of a parameter dependent system which is equivariant under the symmetry group of a regular hexagon. In this paper the mixed symbolic-numerical algorithm SYMCON is applied to analyze its steady state solutions automatically showing their different symmetry and stability properties. T3 - ZIB-Report - SC-91-08 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-587 ER - TY - GEN A1 - Gatermann, Karin A1 - Hohmann, Andreas T1 - Symbolic Exploitation of Symmetry in Numerical Path-following. N2 - Parameter-dependent systems of nonlinear equations with symmetry are treated by a combination of symbolic and numerical computations. In the symbolic part of the algorithm the complete analysis of the symmetry occurs, and it is here where symmetrical normal forms, symmetry reduced systems, and block diagonal Jacobians are computed. Given a particular problem, the symbolic algorithm can create and compute through the list of possible bifurcations thereby forming a so-called tree of decisions correlated to the different types of symmetry breaking bifurcation points. The remaining part of the algorithm deals with the numerical pathfollowing based on the implicit reparametrisation as suggested and worked out by Deuflhard/Fiedler/Kunkel. The symmetry preserving bifurcation points are computed using recently developed augmented systems incorporating the use of symmetry. {\bf Keywords:} pathfollowing, mixed symbolic-numeric algorithm, parameter-dependent, nonlinear systems, linear representations. T3 - ZIB-Report - SC-90-11 KW - pathfollowing KW - mixed symbolic-numeric algorithm KW - parameter dependent KW - nonlinear systems KW - linear representations Y1 - 1990 N1 - No preprint available; the article appeared in: IMPACT Comput. Sci. Eng. 3 (1991) pp. 330-365 ER - TY - THES A1 - Hohmann, Andreas T1 - Inexact Gauss Newton Methods for Parameter Dependent Nonlinear Problems. N2 - A new approach to inexact Gauss Newton methods for the solution of underdetermined nonlinear problems is presented. It is based on a convergence theorem being invariant under affine transformations of the domain and results in an easily implementable accuracy matching strategy for the arising linear subproblems which guarantees the quadratic convergence. Thanks to the weak assumptions on the given nonlinear problem, the results provide a general framework for multilevel Newton and continuation methods. As an example, a new multilevel Newton h-p collocation method for boundary value problems of ordinary differential equations is developed. It combines the inexact Newton method with a linear collocation solver using adaptive refinement and variable orders. The performance of the resulting C++ class library {\sc Cocon} is demonstrated by some numerical examples including singular perturbed problems. In addition, the new method is applied to a realistic railway bogie model in which a branch of periodic solutions emanates from a branch of fixed points at a Hopf bifurcation. T3 - ZIB-Report - TR-93-13 Y1 - 1993 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5049 ER -