A variable stepsize exponential multistep integrator, with contour integral approximation of the operator-valued exponential functions, is proposed for solving semilinear parabolic equations with nonsmooth initial data. By this approach, the exponential 𝑘-step method would have 𝑘 th -order convergence in approximating a mild solution, possibly nonsmooth at the initial time. In consistency with the theoretical analysis, a numerical example shows that the method can achieve high-order convergence in the maximum norm for semilinear parabolic equations with discontinuous initial data.

In this paper, we propose a family of time-stepping schemes for approximating general nonlinear Schrödinger equations. The proposed schemes all satisfy both mass and energy conservation (in a modified form for the latter). Truncation and dispersion error analyses are provided for four proposed schemes. Efficient fixed-point iterative solvers are also constructed to solve the resulting nonlinear discrete problems. As a byproduct, an efficient one-step implementation of the BDF schemes is obtained as well. Extensive numerical experiments are presented to demonstrate the convergence and the capability of capturing the blow-up time of the proposed schemes.

No abstract

A new type of low-regularity integrator is proposed for Navier-Stokes equations, coupled with a stabilized finite element method in space. Unlike the other low-regularity integrators for nonlinear dispersive equations, which are all fully explicit in time, the proposed method is semi-implicit in time in order to preserve the energy-decay structure of NS equations. First-order convergence of the proposed method is established independent of the viscosity coefficient µ, under weaker regularity conditions than other existing numerical methods, including the semi-implicit Euler method and classical exponential integrators. Numerical results show that the proposed method is more accurate than the semi-implicit Euler method in the viscous case µ = O(1), and more accurate than the classical exponential integrator in the inviscid case µ → 0.

No abstract