skip to main content


Title: A novel Decoupled and stable scheme for an anisotropic phase-field dendritic crystal growth model
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
Award ID(s):
1720212
NSF-PAR ID:
10100274
Author(s) / Creator(s):
Date Published:
Journal Name:
Applied mathematics letters
ISSN:
1873-5452
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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
  2. 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
  3. 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.

     
    more » « less
  4. Purpose This paper aims to present an unconditionally energy-stable scheme for approximating the convective heat transfer equation. Design/methodology/approach The scheme stems from the generalized positive auxiliary variable (gPAV) idea and exploits a special treatment for the convection term. The original convection term is replaced by its linear approximation plus a correction term, which is under the control of an auxiliary variable. The scheme entails the computation of two temperature fields within each time step, and the linear algebraic system resulting from the discretization involves a coefficient matrix that is updated periodically. This auxiliary variable is given by a well-defined explicit formula that guarantees the positivity of its computed value. Findings Compared with the semi-implicit scheme and the gPAV-based scheme without the treatment on the convection term, the current scheme can provide an expanded accuracy range and achieve more accurate simulations at large (or fairly large) time step sizes. Extensive numerical experiments have been presented to demonstrate the accuracy and stability performance of the scheme developed herein. Originality/value This study shows the unconditional discrete energy stability property of the current scheme, irrespective of the time step sizes. 
    more » « less
  5. null (Ed.)
    In this paper we propose and analyze a finite difference numerical scheme for the Poisson-Nernst-Planck equation (PNP) system. To understand the energy structure of the PNP model, we make use of the Energetic Variational Approach (EnVarA), so that the PNP system could be reformulated as a non-constant mobility H − 1 H^{-1} gradient flow, with singular logarithmic energy potentials involved. To ensure the unique solvability and energy stability, the mobility function is explicitly treated, while both the logarithmic and the electric potential diffusion terms are treated implicitly, due to the convex nature of these two energy functional parts. The positivity-preserving property for both concentrations, n n and p p , is established at a theoretical level. This is based on the subtle fact that the singular nature of the logarithmic term around the value of 0 0 prevents the numerical solution reaching the singular value, so that the numerical scheme is always well-defined. In addition, an optimal rate convergence analysis is provided in this work, in which many highly non-standard estimates have to be involved, due to the nonlinear parabolic coefficients. The higher order asymptotic expansion (up to third order temporal accuracy and fourth order spatial accuracy), the rough error estimate (to establish the ℓ ∞ \ell ^\infty bound for n n and p p ), and the refined error estimate have to be carried out to accomplish such a convergence result. In our knowledge, this work will be the first to combine the following three theoretical properties for a numerical scheme for the PNP system: (i) unique solvability and positivity, (ii) energy stability, and (iii) optimal rate convergence. A few numerical results are also presented in this article, which demonstrates the robustness of the proposed numerical scheme. 
    more » « less