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 introduce a multiple interval Chebyshev-Gauss-Lobatto spectral collocation method for the initial value problems of the nonlinear ordinary differential equations (ODES). This method is easy to implement and possesses the high order accuracy. In addition, it is very stable and suitable for long time calculations. We also obtain the hp-version bound on the numerical error of the multiple inter...

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 ...

the purpose of this study is to present an approximate numerical method for solving high order linear fredholm-volterra integro-differential equations in terms of rational chebyshev functions under the mixed conditions. the method is based on the approximation by the truncated rational chebyshev series. finally, the effectiveness of the method is illustrated in several numerical examples. the p...

We introduce a method for constructing an element-by-element sparse approximate inverse (SAI) preconditioner designed to be effective in a massively-parallel spectral element modeling environment involving nonsymmetric systems. This new preconditioning approach is based on a spectral optimization of a low-resolution preconditioned system matrix (PSM). We show that the local preconditioning matr...

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.

A short discussion of optimal control methods is presented including indirect, direct shooting, and direct transcription methods. Next the basics of multiple-interval pseudospectral methods are given independent of the numerical scheme to highlight the fundamentals. The two numerical schemes discussed are the Legendre pseudospectral method with LGL nodes and the Chebyshev pseudospectral method ...

In the first part of this paper, we survey results that are associated with three types of Laplacian matrices:difference, normalized, and signless. We derive eigenvalue and eigenvector formulaes for paths and cycles using circulant matrices and present an alternative proof for finding eigenvalues of the adjacency matrix of paths and cycles using Chebyshev polynomials. Even though each results i...