The higher order wave equation of KdV type, which describes many important physical phenomena, has been investigated widely in last several decades. In this work, multisymplectic formulations for the higher order wave equation of KdV type are presented, and the local conservation laws are shown to correspond to certain well-known Hamiltonian functionals. The multi-symplectic discretization of e...

The comprehensive generalization of summation-by-parts of Del Rey Fernández et al. (J. Comput. Phys., 266, 2014) is extended to approximations of second derivatives with variable coefficients. This enables the construction of second-derivative operators with one or more of the following characteristics: i) non-repeating interior stencil, ii) nonuniform nodal distributions, and iii) exclusion of...

We extend the Bloch-decomposition based time-splitting spectral method introduced in an earlier paper [13] to the case of (non-)linear Klein-Gordon equations. This provides us with an unconditionally stable numerical method which achieves spectral convergence in space, even in the case where the periodic coefficients are highly oscillatory and/or discontinuous. A comparison to a traditional pse...

*Correspondence: [email protected] 2School of Science, Jiangnan University, Lihu Road, Wuxi, 214122, China Full list of author information is available at the end of the article Abstract In the study of pattern formation in bi-stable systems, the extended Fisher-Kolmogorov (EFK) equation plays an important role. In this paper, a Fourier pseudo-spectral method for solving the EFK equation in ...

This paper presents a new numerical method to solve time fractional Fokker-Planck equation. The space dimension is discretized to the Gauss-Lobatto points, then we apply pseudo-spectral successive integration matrix for this dimension. This approach shows that with less number of points, we can approximate the solution with more accuracy. The numerical results of the examples are displayed.

Pseudo-spectral approximations are constructed for the model equations describing the population kinetics of human tumor cells in vitro and their responses to radiotherapy or chemotherapy. These approximations are more efficient than finite-difference approximations. The spectral accuracy of the pseudo-spectral method allows us to resolve the model with a much smaller number of spatial grid-poi...

Abstract   The present paper is devoted to implementation of the immersed boundary technique into the Fourier pseudo-spectral solution of the vorticity-velocity formulation of the two-dimensional incompressible Navier-Stokes equations. The immersed boundary conditions are implemented via direct modification of the convection and diffusion terms, and therefore, in contrast to some other similar ...