Abstract:The numerical computation of heteroclinic connections in partial differential equations (PDEs) with a gradient structure, such as those arising in the modeling of phase transitions, is considered. Initially, a scalar reaction diffusion equation is studied; structural assumptions are made on the problem to ensure the existence of absorbing sets and, consequently, a global attractor. As a result of the gradient structure, it is known that, if all equilibria are hyperbolic, the global attractor comprises the set of equilibria and heteroclinic orbits connecting equilibria to one another. Thus it is natural to consider direct approximation of the set of equilibria and the connecting orbits.Results are proved about the Fourier spanning basis for branches of equilibria and also for certain heteroclinic connections; these results exploit the oddness of the nonlinearity. The reaction-diffusion equation is then approximated by a Galerkin spectral discretization to produce a system of ordinary differential equations (ODEs). Analogous results to those holding for the PDE are proved for the ODEs in particular, the existence and structure of the global attractor and appropriate spanning bases for the equilibria and certain heteroclinic connections, are studied. Heteroclinic connections in the system of ODEs are then computed using a generalization of known methods to cope with the gradient structure. Suitable parameterizations of the attractor are introduced and numerical continuation used to find families of connections on the attractor. Special connections, which are stable in certain Fourier spanning bases, are used as starting points for the computations.The methods used allow the calculation of connecting orbits that are unstable as solutions of the initial value problem, and thus provide a computational tool for understanding the dynamics of dissipative problems in a manner that could not be achieved by use of standard initial value methods. Numerical results are given for the Chafee-Infante problem and for the Cahn-Hilliard equation. A one-parameter family of PDEs connecting these two problems is introduced, and it is demonstrated numerically that the global attractor for the Chafee-Infante problem can be continuously deformed into that for the Cahn-Hilliard equation.

Basins of attraction and critical curves for Newton-type methods in a phase equilibrium problem

The Convergence and Fractals of Attractive Basins

Robust Approximation Algorithms for the Detection of Attraction Basins in Dynamical Systems

Multiple nontrivial solutions for a double phase system with concave-convex nonlinearities in subcritical and critical cases

Basins of convergence of equilibrium points in the generalized Hill problem

The Numerical Computation of Heteroclinic Connections in Systems of Gradient Partial Differential Equations

Critical Curve of the Non-Newtonian Polytropic Filtration Equations Coupled via Nonlinear Boundary Flux

Phase-space iterative solvers

Critical curves for degenerate parabolic equations coupled via non-linear boundary flux

On Local Behavior of Newton-Type Methods Near Critical Solutions of Constrained Equations

First and Second Order Unconditionally Energy Stable Schemes for Topology Optimization Based on Phase Field Method

Asymptotic Behavior and Numerical Simulations of a Conservative Phase-Field Model with Two Temperatures

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Global Well-Posedness

Convex Envelope Method for determining liquid multi-phase equilibria in systems with arbitrary number of components

Application of tetragonal curves to coupled Boussinesq equations

Multiple positive solutions for a double phase system with singular nonlinearity

From bubbles to clusters: Multiple solutions to the Allen--Cahn system

Anderson acceleration. Convergence analysis and applications to equilibrium chemistry

Mathematical and Numerical Aspects of a Phase-field Approach to Critical Nuclei Morphology in Solids

Sharp critical mass criteria for weak solutions to a degenerate cross-attraction system

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Local Well-Posedness