pde2path - A Matlab package for continuation and bifurcation in 2D elliptic systems

pde2path is a free and easy to use Matlab continuation/bifurcation package for elliptic systems of PDEs with arbitrary many components, on general two dimensional domains, and with rather general boundary conditions. The package is based on the FEM of the Matlab pdetoolbox, and is explained by a number of examples, including Bratu’s problem, the Schnakenberg model, Rayleigh-Benard convection, and von Karman plate equations. These serve as templates to study new problems, for which the user has to provide, via Matlab function files, a description of the geometry, the boundary conditions, the coefficients of the PDE, and a rough initial guess of a solution. The basic algorithm is a one parameter arclength-continuation with optional bifurcation detection and branch-switching. Stability calculations, error control and mesh-handling, and some elementary time-integration for the associated parabolic problem are also supported. The continuation, branch-switching, plotting etc are performed via Matlab command-line function calls guided by the AUTO style. The software can be downloaded from www.staff.uni-oldenburg.de/hannes.uecker/pde2path , where also an online documentation of the software is provided such that in this paper we focus more on the mathematics and the example systems.

[1]  H. D. Mittelmann.  Multilevel continuation techniques for nonlinear boundary value problems with parameter dependence , 1986 .

[2]  Y. Kuznetsov.  Elements of Applied Bifurcation Theory , 2023, Applied Mathematical Sciences.

[3]  Hannes Uecker,et al.  Coupled Mode Equations and Gap Solitons for the 2D Gross-Pitaevskii equation with a non-separable periodic potential , 2008, 0810.4499.

[4]  W. Beyn,et al.  Continuation of Low-Dimensional Invariant Subspaces in Dynamical Systems of Large Dimension , 2001 .

[5]  V. M. Lashkin.  Two-dimensional multisolitons and azimuthons in Bose-Einstein condensates , 2007, 0710.2745.

[6]  James Demmel,et al.  Continuation of Invariant Subspaces in Large Bifurcation Problems , 2008, SIAM J. Sci. Comput..

[7]  Michael J. Ward,et al.  A Metastable Spike Solution for a Nonlocal Reaction-Diffusion Model , 2000, SIAM J. Appl. Math..

[8]  Eusebius J. Doedel,et al.  Lecture Notes on Numerical Analysis of Nonlinear Equations , 2007 .

[9]  P K Maini,et al.  Bifurcating spatially heterogeneous solutions in a chemotaxis model for biological pattern generation. , 1991, Bulletin of mathematical biology.

[10]  Hannes Uecker,et al.  Numerical Results for Snaking of Patterns over Patterns in Some 2D Selkov-Schnakenberg Reaction-Diffusion Systems , 2013, SIAM J. Appl. Dyn. Syst..

[11]  N. Rashevsky,et al.  Mathematical biology , 1961, Connecticut medicine.

[12]  K. Georg.  MATRIX-FREE NUMERICAL CONTINUATION AND BIFURCATION* , 2001 .

[13]  Y. Kivshar,et al.  Vector azimuthons in two-component Bose-Einstein condensates , 2009 .

[14]  Martin Golubitsky,et al.  Boundary conditions and mode jumping in the buckling of a rectangular plate , 1979 .

[15]  R. Seydel.  Practical Bifurcation and Stability Analysis , 1994 .

[16]  P. Schütz,et al.  Oscillations of fronts and front pairs in two- and three-component reaction-diffusion systems , 1996 .

[17]  S. Gurevich,et al.  Transition from bright to dark dissipative solitons in dielectric barrier gas-discharge , 2007 .

[18]  H. Uecker,et al.  Localized patterns, stationary fronts, and snaking in bistable ranges of spots and stripes , 2013 .

[19]  Willy Govaerts,et al.  Numerical methods for bifurcations of dynamical equilibria , 1987 .

[20]  G. Bratu.  Sur les équations intégrales non linéaires , 1913 .

[21]  Roger Pierre,et al.  Finite element analysis of the buckling and mode jumping of a rectangular plate , 1997 .

[22]  J. Schnakenberg,et al.  Simple chemical reaction systems with limit cycle behaviour. , 1979, Journal of theoretical biology.

[23]  Frank Uhlig,et al.  Numerical Algorithms with C , 1996 .

[24]  C.-S. Chien,et al.  Mode Jumping In The Von Kármán Equations , 2000, SIAM J. Sci. Comput..

[25]  E. Knobloch,et al.  Mode Interactions in Large Aspect Ratio Convection , 1997 .

[26]  Mikhailov,et al.  Cellular structures in catalytic reactions with global coupling. , 1996, Physical review letters.

[27]  E. Allgower,et al.  Numerical Continuation Methods , 1990 .