A Kernel Based Unconditionally Stable Scheme for Nonlinear Parabolic Partial Differential Equations

In this paper,a class of high order numerical schemes is proposed to solve the nonlinear parabolic equations with variable coefficients.This method is based on our previous work[11]for convection-diffusion equations,which relies on a special kernel-based formulation of the solutions and successive convolution.However,disadvantages appear when we extend the previous method to our equations,such as inefficient choice of parameters and unprovable stability for high-dimensional problems.To overcome these difficulties,a new kernel-based formulation is designed to approach the spatial derivatives.It maintains the good properties of the original one,including the high order accuracy and unconditionally stable for one-dimensional problems,hence allowing much larger time step evolution compared with other explicit schemes.In additional,without extra computational cost,the proposed scheme can enlarge the available interval of the special parameter in the formulation,leading to less errors and higher efficiency.Moreover,theoretical investigations indicate that it is unconditionally stable for multi-dimensional problems as well.We present numerical tests for one-and two-dimensional scalar and system,demonstrating the designed high order accuracy and unconditionally stable property of the scheme.

Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method （FEM） based on a proper orthogo- nal decomposition （POD） theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save mem- ory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element （FE） scheme is shown to be un- conditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations （FDEs）.

On precise time integration method for non-classically damped MDOF systems

In the complex mode superposition method, the equations of motion for non-classically damped multiple-degree-of-freedom （MDOF） discrete systems can be transferred into a combination of some generalized SDOF complex oscillators. Based on the state space theory, a precise recurrence relationship for these complex oscillators is set up; then a delicate general solution of non-classically damped MDOF systems, completely in real value form, is presented in this paper. In the proposed method, no calculation of the matrix exponential function is needed and the algorithm is unconditionally stable. A numerical example is given to demonstrate the validity and efficiency of the proposed method.

Alternating segment explicit-implicit scheme for nonlinear third-order KdV equation

A group of asymmetric difference schemes to approach the Korteweg-de Vries （KdV） equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.

ENERGY STABLE NUMERICAL METHOD FOR THE TDGL EQUATION WITH THE RETICULAR FREE ENERGY IN HYDROGEL

Here we focus on the numerical simulation of the phase separation about macromolecule microsphere composite （MMC） hydrogel. The model is based on time-dependent Ginzburg- Landau （TDGL） equation with the reticular free energy. An unconditionally energy stable difference scheme is proposed based on the convex splitting of the corresponding energy functional. In the numerical experiments, we observe that simulating the whole process of the phase separation requires a considerably long time. We also notice that the total free energy changes significantly in initial stage and varies slightly in the following time. Based on these properties, we apply the adaptive time stepping strategy to improve the computational efficiency. It is found that the application of time step adaptivity can not only resolve the dynamical changes of the solution accurately but also significantly save CPU time for the long time simulation.

A High Order Adaptive Time-Stepping Strategy and Local Discontinuous Galerkin Method for the Modified Phase Field Crystal Equation

In this paper,we will develop a first order and a second order convex splitting,and a first order linear energy stable fully discrete local discontinuous Galerkin(LDG)methods for the modified phase field crystal(MPFC)equation.In which,the first order linear scheme is based on the invariant energy quadratization approach.The MPFC equation is a damped wave equation,and to preserve an energy stability,it is necessary to introduce a pseudo energy,which all increase the difficulty of constructing numerical methods comparing with the phase field crystal(PFC)equation.Due to the severe time step restriction of explicit timemarchingmethods,we introduce the first order and second order semi-implicit schemes,which are proved to be unconditionally energy stable.In order to improve the temporal accuracy,the semi-implicit spectral deferred correction(SDC)method combining with the first order convex splitting scheme is employed.Numerical simulations of the MPFC equation always need long time to reach steady state,and then adaptive time-stepping method is necessary and of paramount importance.The schemes at the implicit time level are linear or nonlinear and we solve them by multigrid solver.Numerical experiments of the accuracy and long time simulations are presented demonstrating the capability and efficiency of the proposed methods,and the effectiveness of the adaptive time-stepping strategy.

An Adaptive Time-Stepping Strategy for the Cahn-Hilliard Equation

This paper studies the numerical simulations for the Cahn-Hilliard equation which describes a phase separation phenomenon.The numerical simulation of the Cahn-Hilliardmodel needs very long time to reach the steady state,and therefore large time-stepping methods become useful.The main objective of this work is to construct the unconditionally energy stable finite difference scheme so that the large time steps can be used in the numerical simulations.The equation is discretized by the central difference scheme in space and fully implicit second-order scheme in time.The proposed scheme is proved to be unconditionally energy stable and mass-conservative.An error estimate for the numerical solution is also obtained with second order in both space and time.By using this energy stable scheme,an adaptive time-stepping strategy is proposed,which selects time steps adaptively based on the variation of the free energy against time.The numerical experiments are presented to demonstrate the effectiveness of the adaptive time-stepping approach.