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

From MaRDI portal
Publication:669408

DOI10.1016/J.AMC.2015.05.120zbMATH Open1410.65480arXiv1406.6900OpenAlexW1856394478MaRDI QIDQ669408FDOQ669408


Authors: Christian Kuehn Edit this on Wikidata


Publication date: 15 March 2019

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

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.


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




Recommendations




Cites Work


Cited In (12)

Uses Software





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)