Year of Award
Doctor of Philosophy (PhD)
Department of Mathematics.
Mathematical analysis, Numerical analysis, Simulation methods
Research on interfacial phenomena has a long history, which has attracted tremendous interest in recent years. One of the most successful tools is the phase-ﬁeld approach. As phase-ﬁeld models usually involve very complex dynamics and it is nontrivial to obtain analytical solutions, numerical methods have played an important role in various simulations. This thesis is mainly devoted to developing accurate, eﬃcient and robust numerical methods and the related numerical analysis for three representative phase-ﬁeld models, namely the Allen-Cahn equation, the Cahn-Hilliard equation and the thin ﬁlm models. The ﬁrst part of this thesis is mainly concentrated on the stability analysis for these three models, with particular attention to the Allen-Chan equation. We have established three stability criterion, i.e., nonlinear energy stability, L∞-stability and L2-stability. As shared by most phase-ﬁeld models, one of the intrinsic properties of the Allen- Cahn and the Cahn-Hilliard equations is that they satisfy a nonlinear stability re- lationship, usually expressed as a time-decreasing free energy functional. We have studied several stabilized temporal discretization for both the Allen-Cahn and the Cahn-Hilliard equations so that the relevant nonlinear energy stability can be pre- served. The corresponding temporal discretization schemes are linear and are of second-order accuracy. We also apply multi-step implicit-explicit methods to ap- proximate the Allen-Cahn equation. We demonstrate that by suitably choosing the parameters in multi-step implicit-explicit methods the nonlinear energy stability can be preserved. Apart from studying the energy stability for the Allen-Cahn equation, we also establish the numerical maximum principle for some fully discretized schemes. We further extend our analysis technique to the fractional-in-space Allen-Cahn equation. A more general Allen-Cahn-type equation with a nonlinear degenerate mobility and a logarithmic free energy is also considered. The third stability under investigation is the L2-stability. We prove that the con- tinuum Allen-Cahn equation satisﬁes a uniform Lp-stability. Furthermore, we show that both semi-discretized Fourier Galerkin and Fourier collocation methods can in- herit this stability for p = 2, i.e., L2-stability. Based on the established L2-stability, we accomplish the spectral convergence estimate for the Fourier Galerkin methods. We adopt the second-order Strang splitting schemes in the temporal direction with Fourier collocation methods to demonstrate the uniform L2-stability in the fully dis- cretized scheme. Another contribution of this thesis is to propose a p-adaptive spectral deferred correction methods for the long time simulations for all three models. We develop a high-order accurate and energy stable scheme to simulate the phase-ﬁeld models by combining the semi-implicit spectral deferred correction method and ﬁrst-order stabilized semi-implicit schemes. It is found that the accuracy improvement may aﬀect the overall energy stability. To compromise the accuracy and stability, a local p- adaptive strategy is proposed to enhance the accuracy by sacriﬁcing some local energy stability in an acceptable level. Numerical results demonstrate the high eﬀectiveness of our proposed numerical strategy. Keywords: Phase-ﬁeld models, Allen-Cahn equations, Cahn-Hilliard equations, thin ﬁlm models, nonlinear energy stability, maximum principle, L2-stability, adaptive simulations, stabilized semi-implicit schemes, ﬁnite diﬀerence, Fourier spectral meth- ods, spectral deferred correction methods, convex splitting
Yang, Jiang, "Numerical analysis and simulations for phase-field equations" (2014). Open Access Theses and Dissertations. 29.
The author retains all rights to this work. The author has signed an agreement granting HKBU a non-exclusive license to archive and distribute their thesis.