We consider numerical approximations for a phase-field dendritic crystal growth model, which is a highly nonlinear system that couples the anisotropic Allen–Cahn type equation and the heat equation. By combining the stabilized-Invariant Energy Quadratization method with a novel decoupling technique, the scheme requires solving only a sequence of linear elliptic equations at each time step, making it the first, to the best of the author’s knowledge, totally decoupled, linear, unconditionally energy stable scheme for the model. We further prove the unconditional energy stability rigorously and present various numerical simulations to demonstrate the stability and accuracy.
more »
« less
On a novel full decoupling, linear, second‐order accurate, and unconditionally energy stable numerical scheme for the anisotropic phase‐field dendritic crystal growth model
The anisotropic phase‐field dendritic crystal growth model is a highly nonlinear system that couples the anisotropic Allen–Cahn equation and the thermal equation together. Due to the high anisotropy and nonlinear couplings in the system, how to develop an accurate and efficient, especially a fully decoupled scheme, has always been a challenging problem. To solve the challenge, in this article, we construct a novel fully decoupled numerical scheme which is also linear, energy stable, and second‐order time accurate. The key idea to realize the full decoupling structure is to introduce an ordinary differential equation to deal with the nonlinear coupling terms satisfying the so‐called “zero‐energy‐contribution” property. This scheme is very effective and easy to implement since only a few fully decoupled elliptic equations with constant coefficients need to be solved at each time step. We rigorously prove the solvability of each step and the unconditional energy stability, and perform a large number of numerical simulations in 2D and 3D to demonstrate its stability and accuracy numerically.
more » « less- PAR ID:
- 10449325
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Engineering
- Volume:
- 122
- Issue:
- 16
- ISSN:
- 0029-5981
- Page Range / eLocation ID:
- p. 4129-4153
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this paper, we consider numerical approximations for a dendritic solidification phase field model with melt convection in the liquid phase, which is a highly nonlinear system that couples the anisotropic Allen-Cahn type equation, the heat equation, and the weighted Navier-Stokes equations together. We first reformulate the model into a form which is suitable for numerical approximations and establish the energy dissipative law. Then, we develop a linear, decoupled, and unconditionally energy stable numerical scheme by combining the modified projection scheme for the Navier-Stokes equations, the Invariant Energy Quadratization approach for the nonlinear anisotropic potential, and some subtle explicit-implicit treatments for nonlinear coupling terms. Stability analysis and various numerical simulations are presented.more » « less
-
High order accurate and convergent numerical scheme for the strongly anisotropic Cahn–Hilliard modelmore » « less
Abstract We propose and analyze a second order accurate in time, energy stable numerical scheme for the strongly anisotropic Cahn–Hilliard system, in which a biharmonic regularization has to be introduced to make the equation well‐posed. A convexity analysis on the anisotropic interfacial energy is necessary to overcome an essential difficulty associated with its highly nonlinear and singular nature. The second order backward differentiation formula temporal approximation is applied, combined with Fourier pseudo‐spectral spatial discretization. The nonlinear surface energy part is updated by an explicit extrapolation formula. Meanwhile, the energy stability analysis is enforced by the fact that all the second order functional derivatives of the energy stay uniformly bounded by a global constant. A Douglas‐Dupont type regularization is added to stabilize the numerical scheme, and a careful estimate ensures a modified energy stability with a uniform constraint for the regularization parameter . In turn, the combination with an appropriate treatment for the nonlinear double well potential terms leads to a weakly nonlinear scheme. More importantly, such an energy stability is in terms of the interfacial energy with respect to the original phase variable, which enables us to derive an optimal rate convergence analysis.
-
In this paper, we consider the numerical approximation for a phase field model of the coupled two-phase free flow and two-phase porous media flow. This model consists of Cahn– Hilliard–Navier–Stokes equations in the free flow region and Cahn–Hilliard–Darcy equations in the porous media region that are coupled by seven interface conditions. The coupled system is decoupled based on the interface conditions and the solution values on the interface from the previous time step. A fully discretized scheme with finite elements for the spatial discretization is developed to solve the decoupled system. In order to deal with the difficulties arising from the interface conditions, the decoupled scheme needs to be constructed appropriately for the interface terms, and a modified discrete energy is introduced with an interface component. Furthermore, the scheme is linearized and energy stable. Hence, at each time step one need only solve a linear elliptic system for each of the two decoupled equations. Stability of the model and the proposed method is rigorously proved. Numerical experiments are presented to illustrate the features of the proposed numerical method and verify the theoretical conclusions.more » « less
-
more » « less
A second‐order accurate, linear numerical method is analyzed for the Landau–Lifshitz equation with large damping parameters. This equation describes the dynamics of magnetization, with a non‐convexity constraint of unit length of the magnetization. The numerical method is based on the second‐order backward differentiation formula in time, combined with an implicit treatment for the linear diffusion term from the harmonic mapping part and explicit extrapolation for the nonlinear terms. Afterward, a projection step is applied to normalize the numerical solution at a point‐wise level. This numerical scheme has shown extensive advantages in the practical computations for the physical model with large damping parameters, which comes from the fact that only a linear system with constant coefficients (independent of both time and the updated magnetization) needs to be solved at each time step, and has greatly improved the numerical efficiency. Meanwhile, a theoretical analysis for this linear numerical scheme has not been available. In this paper, we provide a rigorous error estimate of the numerical scheme, in the discrete norm, under suitable regularity assumptions and reasonable ratio between the time step size and the spatial mesh size. In particular, the projection operation is nonlinear, and a stability estimate for the projection step turns out to be highly challenging. Such a stability estimate is derived in details, which will play an essential role in the convergence analysis for the numerical scheme, if the damping parameter is greater than 3.
