We introduce two exponential-type integrators for the “good” Boussinesq equation. They are of orders one and two, respectively, and they require lower spatial regularity of the solution compared to classical exponential integrators. For the first integrator, we prove first-order convergence in H^r for solutions in H^{r+1} with r>1/2. This new integrator even converges (with lower order) in H^r for solutions in H^r, i.e., without any additional smoothness assumptions. For the second integrator, we prove second-order convergence in H^r for solutions in H^{r+3} with r>1/2 and convergence in L^2 for solutions in H^3. Numerical results are reported to illustrate the established error estimates. The experiments clearly demonstrate that the new exponential-type integrators are favorable over classical exponential integrators for initial data with low regularity.

We present and analyze two numerical methods for the logarithmic Schrödinger equation (LogSE) consisting of a regularized splitting method and a regularized conservative Crank–Nicolson finite difference method (CNFD). In order to avoid numerical blow-up and/or to suppress round-off error due to the logarithmic nonlinearity in the LogSE, a regularized logarithmic Schrödinger equation (RLogSE) with a small regularized parameter 0<ε≪1 is adopted to approximate the LogSE with linear convergence rate O(ε). Then we use the Lie–Trotter splitting integrator to solve the RLogSE and establish its error bound O(τ^{1/2}ln(ε^{−1})) with τ>0 the time step, which implies an error bound at O(ε+τ^{1/2}ln(ε^{−1})) for the LogSE by the Lie–Trotter splitting method. In addition, the CNFD is also applied to discretize the RLogSE, which conserves the mass and energy in the discretized level. Numerical results are reported to confirm our error bounds and to demonstrate rich and complicated dynamics of the LogSE.

We present a regularized finite difference method for the logarithmic Schr ̈odingerequation (LogSE) and establish its error bound. Due to the blowup of the logarithmic nonlinearity, i.e., ln ρ→ −\infty when ρ→0+v with ρ = |u|^2 being the density and u being the complex-valuedwave function or order parameter, there are significant difficulties in designing numerical methodsand establishing their error bounds for the LogSE. In order to suppress the roundoff error and toavoid blowup, a regularized LogSE (RLogSE) is proposed with a small regularization parameter 0 < ε << 1 and linear convergence is established between the solutions of RLogSE and LogSE interm of ε. Then a semi-implicit finite difference method is presented for discretizing the RLogSEand error estimates are established in terms of the mesh sizehand time stepτas well as the smallregularization parameterε. Finally numerical results are reported to illustrate our error bounds.

We present two uniformly accurate numerical methods for discretizing the Zakharovsystem (ZS) with a dimensionless parameter 0< ε ≤ 1, which is inversely proportional to theacoustic speed. In the subsonic limit regime, i.e., 0< ε << 1, the solution of ZS propagates waves with O(ε)- andO(1)-wavelengths in time and space, respectively, and/or rapid outgoing initial layerswith speed O(1/ε) in space due to the singular perturbation of the wave operator in ZS and/or theincompatibility of the initial data. By adopting an asymptotic consistent formulation of ZS, wepresent a time-splitting exponential wave integrator (TS-EWI) method via applying a time-splittingtechnique and an exponential wave integrator for temporal derivatives in the nonlinear Schr ̈odingerequation and wave-type equation, respectively. By introducing a multiscale decomposition of ZS, wepropose a time-splitting multiscale time integrator (TS-MTI) method. Both methods are explicitand convergent exponentially in space for all kinds of initial data, which is uniformly for ε ∈ (0,1].The TS-EWI method is simpler to be implemented and it is only uniformly and optimally accuratein time for well-prepared initial data, while the TS-MTI method is uniformly and optimally accuratein time for any kind of initial data. Extensive numerical results are reported to show their efficiencyand accuracy, especially in the subsonic limit regime. Finally, the TS-MTI method is applied tostudy numerically convergence rates of ZS to its limiting models when ε→0+.

We propose an explicit numerical method for the periodic Korteweg–de Vries equation. Our method is based on a Lawson-type exponential integrator for time integration and the Rusanov scheme for Burgers’ nonlinearity. We prove first-order convergence in both space and time under a mild Courant–Friedrichs–Lewy condition τ=O(h)⁠, where τ and h represent the time step and mesh size for solutions in the Sobolev space H^3((−π,π))⁠, respectively. Numerical examples illustrating our convergence result are given.