Search

We propose an energy-optimized invariant energy quadratization method to solve the gradient flow models in this paper, which requires only one linear energy-optimized step to correct the auxiliary variables on each time step. In addition to inheriting the benefits of the baseline and relaxed invariant energy quadratization method, our approach has several other advantages. Firstly, in the process of correcting auxiliary variables, we can directly solve linear programming problem by the energy-optimized technique, which greatly simplifies the nonlinear optimization problem in the previous relaxed invariant energy quadratization method. Secondly, we construct new linear unconditionally energy stable schemes by applying backward Euler formulas and Crank-Nicolson formula, so that the accuracy in time can reach the first- and second-order. Thirdly, comparing with relaxation technique, the modified energy obtained by energy-optimized technique is closer to the original energy, meanwhile the accuracy and consistency of the numerical solutions can be improved. Ample numerical examples have been presented to demonstrate the accuracy, efficiency and energy stability of the proposed schemes., Comment: 35 pages, 39 figures

The nonlocal Allen-Cahn equation with nonlocal diffusion operator is a generalization of the classical Allen-Cahn equation. It satisfies the energy dissipation law and maximum bound principle (MBP), and is important for simulating a series of physical and biological phenomena involving long-distance interactions in space. In this paper, we construct first- and second-order (in time) accurate, unconditionally energy stable and MBP-preserving schemes for the nonlocal Allen-Cahn type model based on the stabilized exponential scalar auxiliary variable (sESAV) approach. On the one hand, we have proved the MBP and unconditional energy stability carefully and rigorously in the fully discrete levels. On the other hand, we adopt an efficient FFT-based fast solver to compute the nearly full coefficient matrix generated from the spatial discretization, which improves the computational efficiency. Finally, typical numerical experiments are presented to demonstrate the performance of our proposed schemes.

The nonlocal Cahn-Hilliard (NCH) equation with nonlocal diffusion operator is more suitable for the simulation of microstructure phase transition than the local Cahn-Hilliard (LCH) equation. In this paper, based on the exponential semi-implicit scalar auxiliary variable (ESI-SAV) method, the highly effcient and accurate schemes in time with unconditional energy stability for solving the NCH equation are proposed. On the one hand, we have demostrated the unconditional energy stability for the NCH equation with its high-order semi-discrete schemes carefully and rigorously. On the other hand, in order to reduce the calculation and storage cost in numerical simulation, we use the fast solver based on FFT and FCG for spatial discretization. Some numerical simulations involving the Gaussian kernel are presented and show the stability, accuracy, efficiency and unconditional energy stability of the proposed schemes.