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.

In [H. Guo, Q. Zhang, J. Wang, Applied Mathematics and Computation, 259 (2015), 88-105], a nonconservative local discontinuous Galerkin (LDG) method for both flow and transport equations was introduced for the one-dimensional coupled system of compressible miscible displacement problem. In this paper, we will continue our effort and develop a conservative LDG method for the problem in two space dimensions. Optimal error estimates in $L^{\infty}(0, T; L^{2})$ norm for not only the solution itself but also the auxiliary variables will be derived. The main difficulty is how to treat the inter-element discontinuities of two independent solution variables (one from the flow equation and the other from the transport equation) at cell interfaces. Numerical experiments will be given to confirm the accuracy and efficiency of the scheme.

No abstract uploaded!

No abstract uploaded!

In this article the authors present a new solver for the Poisson problem on spherical shells. The method is derived by application of the Fast Fourier Transform (FFT) in all three spherical variables. Through a particular change of variables valid on spherical shells, this new approach avoids a variable coefficient on the radial differential operator, resulting in a constant-coefficient formulation that can be fast-diagonalized via FFT. The authors describe the Fourier expansion in the angular variables and the fast diagonalization of the radial derivative using the FFT, including second-, fourth- and sixth-order finite difference approximations. An algorithmic description of the sixth-order solver is described, and numerical tests that demonstrate the theoretical rate of convergence are presented.