An official website of the United States government
In contrast to the classical Allen-Cahn equation, the conservative Allen-Cahn equation with a nonlocal Lagrange multiplier not only satisfies the maximum bound principle (MBP) and energy dissipation law but also ensures mass conservation. Many existing schemes often fail to preserve all these properties at the discrete level or require high regularity in time on the exact solution for convergence analysis. In this paper, we construct a new class of low regularity integrators (LRIs) for time discretization of the conservative Allen-Cahn equation by repeatedly using Duhamel's formula. The proposed first- and second-order LRI schemes are shown to conserve mass unconditionally and satisfy the MBP under some time step size constraints. Temporal error estimates for these schemes are derived under a low regularity assumption that the exact solution is only Lipschitz continuous in time, followed by a rigorous proof for energy stability of the corresponding time-discrete solutions. Various numerical experiments and comparisons in two and three dimensions are presented to verify the theoretical results and illustrate the performance of the LRI schemes, especially when the interfacial parameter approaches zero.
more »
« less
Second-order stabilized semi-implicit energy stable schemes for bubble assemblies in binary and ternary systems
In this paper, we propose some second-order stabilized semi-implicit methods for solving the Allen-Cahn-Ohta-Kawasaki and the Allen-Cahn-Ohta-Nakazawa equations. In the numerical methods, some nonlocal linear stabilizing terms are introduced and treated implicitly with other linear terms, while other nonlinear and nonlocal terms are treated explicitly. We consider two different forms of such stabilizers and compare the difference regarding the energy stability. The spatial discretization is performed by the Fourier collocation method with FFT-based fast implementations. Numerically, we verify the second order temporal convergence rate of the proposed schemes. In both binary and ternary systems, the coarsening dynamics is visualized as bubble assemblies in hexagonal or square patterns.
more »
« less
- Award ID(s):
- 2142500
- PAR ID:
- 10325843
- Date Published:
- Journal Name:
- Discrete & Continuous Dynamical Systems - B
- Volume:
- 0
- Issue:
- 0
- ISSN:
- 1531-3492
- Page Range / eLocation ID:
- 0
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The physics-informed neural networks (PINNs) has been widely utilized to numerically approximate PDE problems. While PINNs has achieved good results in producing solutions for many partial differential equations, studies have shown that it does not perform well on phase field models. In this paper, we partially address this issue by introducing a modified physics-informed neural networks. In particular, they are used to numerically approximate Allen- Cahn-Ohta-Kawasaki (ACOK) equation with a volume constraint. Technically, the inverse of Laplacian in the ACOK model presents many challenges to the baseline PINNs. To take the zero- mean condition of the inverse of Laplacian, as well as the volume constraint, into consideration, we also introduce a separate neural network, which takes the second set of sampling points in the approximation process. Numerical results are shown to demonstrate the effectiveness of the modified PINNs. An additional benefit of this research is that the modified PINNs can also be applied to learn more general nonlocal phase-field models, even with an unknown nonlocal kernel.more » « less
-
TheIsing modelof statistical physics has served as a keystone example of phase transitions, thermodynamic limits, scaling laws, and many other phenomena and mathematical methods. We introduce and explore anIsing game, a variant of the Ising model that features competing agents influencing the behavior of the spins. With long-range interactions, we consider a mean-field limit resulting in a nonlocal potential game at the mesoscopic scale. This game exhibits a phase transition and multiple constant Nash-equilibria in the supercritical regime. Our analysis focuses on a sharp interface limit for which potential minimizing solutions to the Ising game concentrate on two of the constant Nash-equilibria. We show that the mesoscopic problem can be recast as a mixed local/nonlocal space-time Allen-Cahn type minimization problem. We prove, using a Γ-convergence argument, that the limiting interface minimizes a space-time anisotropic perimeter type energy functional. This macroscopic scale problem could also be viewed as a problem of optimal control of interface motion. Sharp interface limits of Allen-Cahn type functionals have been well studied. We build on that literature with new techniques to handle a mixture of local derivative terms and nonlocal interactions. The boundary conditions imposed by the game theoretic considerations also appear as novel terms and require special treatment.more » « less
-
null (Ed.)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.more » « less
-
In this paper we study the global well-posedness of the Allen–Cahn–Ohta–Nakazawa model with two fixed nonlinear volume constraints. Utilizing the gradient flow structure of its free energy, we prove the existence and uniqueness of the solution by following De Giorgi’s minimizing movement scheme in a novel way.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3934/dcdsb.2021246
