Efficient gluing of numerical continuation and a multiple solution method for elliptic PDEs

DOI10.1016/J.AMC.2015.05.120MaRDI QIDQ669408zbMATHOpenAlexFDO

Authors Christian Kuehn

Publication date 15 March 2019

Published in Applied Mathematics and Computation (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1406.6900

zbMATH Keywords

minimaxbifurcationnumerical continuationsymmetry-breakinggluing computationpde2path

Mathematics Subject Classification ID

Semilinear elliptic equations (35J61) Bifurcations in context of PDEs (35B32) Boundary value problems for second-order elliptic equations (35J25) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Numerical bifurcation problems (65P30)

Abstract: Numerical continuation calculations for ordinary differential equations (ODEs) are, by now, an established tool for bifurcation analysis in dynamical systems theory as well as across almost all natural and engineering sciences. Although several excellent standard software packages are available for ODEs, there are - for good reasons - no standard numerical continuation toolboxes available for partial differential equations (PDEs), which cover a broad range of different classes of PDEs automatically. A natural approach to this problem is to look for efficient gluing computation approaches, with independent components developed by researchers in numerical analysis, dynamical systems, scientific computing and mathematical modelling. In this paper, we shall study several elliptic PDEs (Lane-Emden-Fowler, Lane-Emden-Fowler with microscopic force, Caginalp) via the numerical continuation software pde2path and develop a gluing component to determine a set of starting solutions for the continuation by exploting the variational structures of the PDEs. In particular, we solve the initialization problem of numerical continuation for PDEs via a minimax algorithm to find multiple unstable solution. Furthermore, for the Caginalp system, we illustrate the efficient gluing link of pde2path to the underlying mesh generation and the FEM MatLab pdetoolbox. Even though the approach works efficiently due to the high-level programming language and without developing any new algorithms, we still obtain interesting bifurcation diagrams and directly applicable conclusions about the three elliptic PDEs we study, in particular with respect to symmetry-breaking. In particular, we show for a modified Lane-Emden-Fowler equation with an asymmetric microscopic force, how a fully connected bifurcation diagram splits up into C-shaped isolas on which localized pattern deformation appears towards two different regimes.

Recommendations

• pde2path -- a Matlab package for continuation and bifurcation in 2D elliptic systems
• Continuation and bifurcation in nonlinear PDEs - algorithms, applications, and experiments
• A parallel and adaptive continuation method for semilinear bifurcation problems
• Application of numerical continuation to compute all solutions of semilinear elliptic equations
• On homotopy continuation method for computing multiple solutions to the Henon equation

Cites work

• scientific article; zbMATH DE number 3130242 (Why is no real title available?)
• scientific article; zbMATH DE number 3930251 (Why is no real title available?)
• scientific article; zbMATH DE number 3720301 (Why is no real title available?)
• scientific article; zbMATH DE number 3788145 (Why is no real title available?)
• scientific article; zbMATH DE number 47535 (Why is no real title available?)
• scientific article; zbMATH DE number 3640828 (Why is no real title available?)
• scientific article; zbMATH DE number 1215244 (Why is no real title available?)
• scientific article; zbMATH DE number 1234388 (Why is no real title available?)
• scientific article; zbMATH DE number 2000348 (Why is no real title available?)
• scientific article; zbMATH DE number 1522917 (Why is no real title available?)
• scientific article; zbMATH DE number 2163592 (Why is no real title available?)
• scientific article; zbMATH DE number 3806869 (Why is no real title available?)
• scientific article; zbMATH DE number 1569072 (Why is no real title available?)
• scientific article; zbMATH DE number 815352 (Why is no real title available?)
• A Local Minimax-Newton Method for Finding Multiple Saddle Points with Symmetries
• A NUMERICAL TOOLBOX FOR HOMOCLINIC BIFURCATION ANALYSIS
• A SURVEY OF METHODS FOR COMPUTING (UN)STABLE MANIFOLDS OF VECTOR FIELDS
• A high-linking algorithm for sign-changing solutions of semilinear elliptic equations
• A local min-max-orthogonal method for finding multiple solutions to noncooperative elliptic systems
• A minimax method for finding multiple critical points and its applications to semilinear PDEs
• A mountain pass method for the numerical solution of semilinear elliptic problems
• ALGORITHMS AND VISUALIZATION FOR SOLUTIONS OF NONLINEAR ELLIPTIC EQUATIONS
• Adaptive Finite Element Methods for Parabolic Problems I: A Linear Model Problem
• An analysis of a phase field model of a free boundary
• An optimal bound on the number of interior spike solutions for the Lin-Ni-Takagi problem
• An overview of the Trilinos project
• BIFURCATION TRACKING ALGORITHMS AND SOFTWARE FOR LARGE SCALE APPLICATIONS
• Bifurcation: analysis, algorithms, applications. Proceedings of the Conference at the University of Dortmund, August 18-22, 1986
• Computing Geodesic Level Sets on Global (Un)stable Manifolds of Vector Fields
• Computing absolute and essential spectra using continuation
• Continuation of Bifurcations in Periodic Delay‐Differential Equations Using Characteristic Matrices
• Continuation of Quasi-periodic Invariant Tori
• Convergence Results of a Local Minimax Method for Finding Multiple Critical Points
• Deterministic continuation of stochastic metastable equilibria via Lyapunov equations and ellipsoids
• Elements of applied bifurcation theory
• Elliptic partial differential equations of second order
• Finite Volume Methods for Hyperbolic Problems
• From first Lyapunov coefficients to maximal canards
• Generalized Ginzburg-Landau and Cahn-Hilliard equations based on a microforce balance
• Introduction to Numerical Continuation Methods
• Isogeometric analysis. Toward integration of CAD and FEA
• Isolas of 2-pulse solutions in homoclinic snaking scenarios
• Large amplitude stationary solutions to a chemotaxis system
• Localized Hexagon Patterns of the Planar Swift–Hohenberg Equation
• Long time behavior of solutions to the Caginalp system with singular potential
• MATCONT
• MULTIPLE PARAMETER CONTINUATION: COMPUTING IMPLICITLY DEFINED k-MANIFOLDS
• Multi-bubble nodal solutions for slightly subcritical elliptic problems in domains with symmetries
• NUMERICAL ANALYSIS AND CONTROL OF BIFURCATION PROBLEMS (I): BIFURCATION IN FINITE DIMENSIONS
• NUMERICAL ANALYSIS AND CONTROL OF BIFURCATION PROBLEMS (II): BIFURCATION IN INFINITE DIMENSIONS
• Non-existence of positive solutions of Lane-Emden systems
• Numerical Bifurcation Methods and their Application to Fluid Dynamics: Analysis beyond Simulation
• Numerical bifurcation study of superconducting patterns on a square
• Numerical computation of bifurcations in large equilibrium systems in \textsc{Matlab}
• Numerical continuation methods for dynamical systems. Path following and boundary value problems.
• Numerical continuation of canard orbits in slow-fast dynamical systems
• Numerical results for snaking of patterns over patterns in some 2D Selkov-Schnakenberg reaction-diffusion systems
• ON THE NUMERICAL CONTINUATION OF ISOLAS OF EQUILIBRIA
• On finding multiple solutions to a singularly perturbed Neumann problem
• On homotopy continuation method for computing multiple solutions to the Henon equation
• On the Generalized Emden–Fowler Equation
• On the number of interior peak solutions for a singularly perturbed Neumann problem
• Practical bifurcation and stability analysis: from equilibrium to chaos.
• Recipes for continuation
• Riemannian geometry and geometric analysis
• Snakes, Ladders, and Isolas of Localized Patterns
• Snaking and isolas of localised states in bistable discrete lattices
• Snaking of radial solutions of the multi-dimensional Swift-Hohenberg equation: a numerical study
• Some generalizations of the Caginalp phase-field system
• The Bordering Algorithm and Path Following Near Singular Points of Higher Nullity
• The Caginalp phase-field system with coupled dynamic boundary conditions and singular potentials
• The Mathematical Theory of Finite Element Methods
• The Moment Map: Nonlinear Dynamics of Density Evolution via a Few Moments
• The Numerical Computation of Connecting Orbits in Dynamical Systems
• The geometry of mixed-mode oscillations in the Olsen model for the peroxidase-oxidase reaction
• The immersed boundary method
• The world of the complex Ginzburg-Landau equation
• Time integration and steady-state continuation for 2D lubrication equations
• To Snake or Not to Snake in the Planar Swift–Hohenberg Equation
• Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
• pde2path -- a Matlab package for continuation and bifurcation in 2D elliptic systems

Cited in

(15)

• Hopf Bifurcation and Time Periodic Orbits with pde2path – Algorithms and Applications
• On the influence of cross-diffusion in pattern formation
• A meeting point of entropy and bifurcations in cross-diffusion herding
• A \(J\)-\textit{integral}-based arc-length solver for brittle and ductile crack propagation in finite deformation-finite strain hyperelastic solids with an application to graphene kirigami
• pde2path -- a Matlab package for continuation and bifurcation in 2D elliptic systems
• Numerical continuation for a fast-reaction system and its cross-diffusion limit
• Pseudospectral methods for computing the multiple solutions of the Schrödinger equation
• On numerical solutions for a class of relativistic quasilinear Schrödinger equations
• Random attractors for stochastic partly dissipative systems
• Numerical challenges for the understanding of localised solutions with different symmetries in non-local hyperbolic systems
• Mixed Fourier-Legendre spectral methods for the multiple solutions of the Schrödinger equation on the unit disk
• Numerical continuation and SPDE stability for the 2D cubic-quintic Allen-Cahn equation
• Methods of continuation and their implementation in the COCO software platform with application to delay differential equations
• Numerical continuation and bifurcation for differential geometric PDEs
• Continuation and bifurcation in nonlinear PDEs - algorithms, applications, and experiments

Describes a project that uses

Uses Software

• MATCONT
• Matlab
• AUTO
• DDE-BIFTOOL
• Trilinos
• FEniCS
• AUTO-07P
• XPPAUT
• COCO
• CL_MATCONTL
• PyDSTool
• CL_MATCONT
• pde2path
• COLSYS
• FiPy
• HomCont
• Torcont
• SfePy
• Python

This page was built for publication: Efficient gluing of numerical continuation and a multiple solution method for elliptic PDEs

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q669408)