Zhengguang Liu

Zhengguang Liu, Xiaoli Li

In this paper, we propose a novel Lagrange Multiplier approach, named zero-factor (ZF) approach to solve a series of gradient flow problems. The numerical schemes based on the new algorithm are unconditionally energy stable with the original energy and do not require any extra assumption conditions. We also prove that the ZF schemes with specific zero factors lead to the popular SAV-type method. To reduce the computation cost and improve the accuracy and consistency, we propose a zero-factor approach with relaxation, which we named the relaxed zero-factor (RZF) method, to design unconditional energy stable schemes for gradient flows. The RZF schemes can be proved to be unconditionally energy stable with respect to a modified energy that is closer to the original energy, and provide a very simple calculation process. The variation of the introduced zero factor is highly consistent with the nonlinear free energy which implies that the introduced ZF method is a very efficient way to capture the sharp dissipation of nonlinear free energy. Several numerical examples are provided to demonstrate the improved efficiency and accuracy of the proposed method.

Zhengguang Liu, Aijie Cheng, Xiaoli Li et al.

In this paper, we develop fast procedures for solving linear systems arising from discretization of ordinary and partial differential equations with Caputo fractional derivative w.r.t time variable. First, we consider a finite difference scheme to solve a two-sided fractional ordinary equation. Furthermore, we present a fast solution technique to accelerate Toeplitz matrix-vector multiplications arising from finite difference discretization. This fast solution technique is based on a fast Fourier transform and depends on the special structure of coefficient matrices, and it helps to reduce the computational work from $O(N^{3})$ required by traditional methods to $O(Nlog^{2}N)$ and the memory requirement from $O(N^{2})$ to $O(N)$ without using any lossy compression, where $N$ is the number of unknowns. Two finite difference schemes to solve time fractional hyperbolic equations with different fractional order $γ$ are considered. We present a fast solution technique depending on the special structure of coefficient matrices by rearranging the order of unknowns. It helps to reduce the computational work from $O(N^2M)$ required by traditional methods to $O(N$log$^{2}N)$ and the memory requirement from $O(NM)$ to $O(N)$ without using any lossy compression, where $N=τ^{-1}$ and $τ$ is the size of time step, $M=h^{-1}$ and $h$ is the size of space step. Importantly, a fast method is employed to solve the classical time fractional diffusion equation with a lower coast at $O(MN$log$^2N)$, where the direct method requires an overall computational complexity of $O(N^2M)$. Moreover, the applicability and accuracy of the scheme are demonstrated by numerical experiments to support our theoretical analysis.