Search for: All records

We present a rigorous analysis of the transient evolution of nearly circular bilayer interfaces evolving under the thin interface limit, ε ≪ 1, of the mass preserving L2-gradient flow of the strong scaling of the functionalized Cahn–Hilliard equation. For a domain Ω ⊂ R2 we construct a bilayer manifold with boundary comprised of quasi-equilibria of the flow and a projection onto the manifold that associates functions u in an H2 tubular neighborhood of the manifold with an interface Γ embedded in Ω. The linearization of the flow about the manifold does not present a clear spectral separation of modes normal and tangential to the manifold. The dimension of the parameterization of the interfaces and the bilayer manifold controls both the normal coercivity of the manifold and the coupling between normal and tangential modes, both of which increase with this dimension. The key step in the analysis is the identification of a range of dimensions in which coercivity dominates the coupling, permitting the closure of the nonlinear estimates that establish the asymptotic stability of the manifold. Orbits originating in a thin, forward invariant, tubular neighborhood ultimately converge to an equilibrium associated to a circular interface. Projections of these orbits yield interfacial evolution equivalent at leading order to the regularized curve-lengthening motion characterized by normal motion against mean curvature, regularized by a higher order Willmore expression. The curve lengthening is driven by absorption of excess mass from the regions of Ω away from the interface, leading to high dimensional dynamics that are ill-posed in the ε → 0+ limit. 

Abstract Adaptive time stepping methods for metastable dynamics of the Allen–Cahn and Cahn–Hilliard equations are investigated in the spatially continuous, semi-discrete setting. We analyse the performance of a number of first and second order methods, formally predicting step sizes required to satisfy specified local truncation error $$\sigma $$ σ in the limit of small length scale parameter $$\epsilon \rightarrow 0$$ ϵ → 0 during meta-stable dynamics. The formal predictions are made under stability assumptions that include the preservation of the asymptotic structure of the diffuse interface, a concept we call profile fidelity. In this setting, definite statements about the relative behaviour of time stepping methods can be made. Some methods, including all so-called energy stable methods but also some fully implicit methods, require asymptotically more time steps than others. The formal analysis is confirmed in computational studies. We observe that some provably energy stable methods popular in the literature perform worse than some more standard schemes. We show further that when Backward Euler is applied to meta-stable Allen–Cahn dynamics, the energy decay and profile fidelity properties for these discretizations are preserved for much larger time steps than previous analysis would suggest. The results are established asymptotically for general interfaces, with a rigorous proof for radial interfaces. It is shown analytically and computationally that for most reaction terms, Eyre type time stepping performs asymptotically worse due to loss of profile fidelity. 

The strong functionalized Cahn–Hilliard equation models self assembly of amphiphilic polymers in solvent. It supports codimension one and two structures that each admit two classes of bifurcations: pearling, a short-wavelength in-plane modulation of interfacial width, and meandering, a long-wavelength instability that induces a transition to curve-lengthening flow. These two potential instabilities afford distinctive routes to changes in codimension and creation of non-codimensional defects such as end caps and Y-junctions. Prior work has characterized the onset of pearling, showing that it couples strongly to the spatially constant, temporally dynamic, bulk value of the chemical potential. We present a multiscale analysis of the competitive evolution of codimension one and two structures of amphiphilic polymers within the H−1 gradient flow of the strong Functionalized Cahn–Hilliard equation. Specifically we show that structures of each codimension transition from a curve lengthening to a curve shortening flow as the chemical potential falls through a corresponding critical value. The differences in these critical values quantify the competition between the morphologies of differing codimension for the amphiphilic polymer mass. We present a bifurcation diagram for the morphological competition and compare our results quantitatively to simulations of the full system and qualitatively to simulations of self-consistent mean field models and laboratory experiments. In particular we propose that the experimentally observed onset of morphological complexity arises from a transient passage through pearling instability while the associated flow is in the curve lengthening regime. 

There is a large literature of numerical methods for phase field models from materials science. The prototype models are the Allen-Cahn and Cahn-Hilliard equa- tions. We present four benchmark problems for these equations, with numerical results validated using several computational methods with different spatial and temporal discretizations. Our goal is to provide the scientific community with a reliable reference point for assessing the accuracy and reliability of future software for this important class of problem.