IDEAS home Printed from https://ideas.repec.org/a/eee/apmaco/v266y2015icp656-674.html
   My bibliography  Save this article

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

Author

Listed:
  • Kuehn, Christian

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 ach 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 modeling. 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 exploiting 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. We conclude with a section on future software development issues that would be helpful to be addressed to simplify interfaces to allow for more efficient, time-saving, gluing computation for dynamical systems analysis of PDEs in the near future.

Suggested Citation

  • Kuehn, Christian, 2015. "Efficient gluing of numerical continuation and a multiple solution method for elliptic PDEs," Applied Mathematics and Computation, Elsevier, vol. 266(C), pages 656-674.
  • Handle: RePEc:eee:apmaco:v:266:y:2015:i:c:p:656-674
    DOI: 10.1016/j.amc.2015.05.120
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0096300315007511
    Download Restriction: Full text for ScienceDirect subscribers only

    File URL: https://libkey.io/10.1016/j.amc.2015.05.120?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    As the access to this document is restricted, you may want to search for a different version of it.

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Christian Kuehn & Cinzia Soresina, 2020. "Numerical continuation for a fast-reaction system and its cross-diffusion limit," Partial Differential Equations and Applications, Springer, vol. 1(2), pages 1-26, April.

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:apmaco:v:266:y:2015:i:c:p:656-674. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    We have no bibliographic references for this item. You can help adding them by using this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Catherine Liu (email available below). General contact details of provider: https://www.journals.elsevier.com/applied-mathematics-and-computation .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.