Concatenating construction of the multisymplectic schemes for 2+1-dimensional sine-Gordon equation

In this paper, taking the 2+1-dimensional sine-Gordon equation as an example, we present the concatenating method to construct the multisymplectic schemes. The-method is to discretizee independently the PDEs in different directions with symplectic schemes, so that the multisymplectic schemes can be constructed by concatenating those symplectic schemes. By this method, we can construct multisymplectic schemes, including some widely used schemes with an accuracy of any order. The numerical simulation on the collisions of solitons are also proposed to illustrate the efficiency of the multisymplectic schemes.     

Multisymplectic Fourier pseudo-spectral integrators for Klein-Gordon-Schrödinger equations

A multisymplectic Fourier pseudo-spectral scheme,which exactly preserves the discrete multisymplectic conservation law,is presented to solve the Klein-Gordon-Schrdinger equations.The scheme is of spectral accuracy in space and of second order in time.The scheme preserves the discrete multisymplectic conservation law and the charge conservation law.Moreover,the residuals of some other conservation laws are derived for the geometric numerical integrator.Extensive numerical simulations illustrate the numerical behavior of the multisymplectic scheme,and demonstrate the correctness of the theoretical analysis.     

Multisymplectic Structure and Multisymplectic Scheme for the Nonlinear Ware Equation

Abstract The multisvmplectic structure of the nonlinear wave equation is derived directly from the variationalprinciple. In the numerical aspect,we present a multisymplectic nine points scheme which is equivalent to themultisymplectic Preissman scheme.A series of numerical results are reported to illustrate the effectiveness ofthe scheme.     

A multisymplectic explicit scheme for the modified regularized long-wave equation

In this paper, we derive a new 10-point multisymplectic scheme for the modified regularized long-wave equation. The new scheme is an explicit scheme in the sense that the third time level does not include nonlinear terms. Numerical results indicate that the new scheme not only provides satisfied numerical solutions, but also preserves three invariants of motion very well.     

Multisymplecticity of the centred box scheme for a class of hamiltonian PDEs and an application to quasi-periodically solitary waves

In this paper, the multisymplectic integrator for a class of Hamiltonian PDEs depending explicitly on time and spatial variables (nonautonomous Hamiltonian PDEs) is defined, and the multisymplecticity of the centred box scheme for this kind of Hamiltonian PDEs is proven. We give an application of the result to (periodic) quasi-periodic variable coefficient Korteweg-de Vries (qpKdV) equation, which is known to have a physical application in the propagation of surface waves in straits or channels with quasi-periodic varying depth and width in the time direction. We derive a multisymplectic scheme for a qpKdV equation in terms of the multisymplecticity of the centred box scheme, then make use of it to simulate numerically the (periodically) quasi-periodically solitary wave of the equation. Numerical experiments are presented in illustration of the multisymplectic scheme of qpKdV equation stemming the centred box discretization.     

MULTISYMPLECTIC FOURIER PSEUDOSPECTRAL METHOD FOR THE NONLINEAR SCHRODINGER EQUATIONS WITH WAVE OPERATOR

In this paper, the multisymplectic Fourier pseudospectral scheme for initial-boundary value problems of nonlinear SchrSdinger equations with wave operator is considered. We investigate the local and global conservation properties of the multisymplectic discretization based on Fourier pseudospectral approximations. The local and global spatial conservation of energy is proved. The error estimates of local energy conservation law are also derived. Numerical experiments are presented to verify the theoretical predications.     

A multisymplectic integrator for the periodic nonlinear Schrödinger equation

A multisymplectic integrator for the periodic nonlinear Schrödinger equation is presented in this paper. Its accuracy is proved. By introducing a norm, we investigate its nonlinear stability. We also discuss the relationship between this multisymplectic integrator and two variational integrators which are derived by using the discrete multisymplectic field theory and the finite element method.     

Acceleration of partitioned coupling schemes for problems of thermoelasticity

A partitioned coupling scheme for problems of thermo-elasticity at finite strains is presented. The coupling between the mechanical and thermal field is one of the most important multi-physics problem. Typically two different strategies are used to find an accurate solution for both fields: Partitioned or staggered coupling schemes, in which the mechanics and heat transfer is treated as a single field problem, or a monolithic solution of the full problem. Monolithic formulations have the drawback of a non-symmetric system which may lead to extremely large computational costs. Because partitioned schemes avoid this problem and allow for numerical formulations which are more flexible, we consider a staggered coupling algorithm which decouples the mechanical and the thermal field into partitioned symmetric sub-problems by means of an isothermal operator-split. In order to stabilize and to accelerate the convergence of the partitioned scheme, two different methods are employed: dynamic relaxation and a reduced order model quasi-Newton method. A numerical simulation of a quasi-static problem is presented investigating the performance of accelerated coupling schemes. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Error analysis of discrete conservation laws for Hamiltonian PDEs under the central box discretizations

The central box scheme has been the most successful of the multisymplectic integrators for Hamiltonian PDEs. In this paper, we investigate conservative properties of the central box scheme for Hamiltonian PDEs and derive the error formulas of discrete local and global conservation laws of energy and momentum. We apply these results to the nonlinear Schrödinger equation and Klein-Gordon equation. Numerical experiments are presented to verify the theoretical predications.     

Efficient local structure‐preserving schemes for the RLW‐type equation

The multisymplectic schemes have been used in numerical simulations for the RLW‐type equation successfully. They well preserve the local geometric property, but not other local conservation laws. In this article, we propose three novel efficient local structure‐preserving schemes for the RLW‐type equation, which preserve the local energy exactly on any time‐space region and can produce richer information of the original problem. The schemes will be mass‐ and energy‐preserving as the equation is imposed on appropriate boundary conditions. Numerical experiments are presented to verify the efficiency and invariant‐preserving property of the schemes. Comparisons with the existing nonconservative schemes are made to show the behavior of the energy affects the behavior of the solution.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1678–1691, 2017     

Two Novel Classes of Arbitrary High-Order Structure-Preserving Algorithms for Canonical Hamiltonian Systems

In this paper, we systematically construct two classes of structure-preserving schemes with arbitrary order of accuracy for canonical Hamiltonian systems. The one class is the symplectic scheme, which contains two new families of parameterized symplectic schemes that are derived by basing on the generating function method and the symmetric composition method, respectively. Each member in these schemes is symplectic for any fixed parameter. A more general form of generating functions is introduced, which generalizes the three classical generating functions that are widely used to construct symplectic algorithms. The other class is a novel family of energy and quadratic invariants preserving schemes, which is devised by adjusting the parameter in parameterized symplectic schemes to guarantee energy conservation at each time step. The existence of the solutions of these schemes is verified. Numerical experiments demonstrate the theoretical analysis and conservation of the proposed schemes.     

MKdV方程的多辛Fourier拟谱方法

基于谱微分矩阵方法,给出MKdV方程的多辛Fourier拟谱格式及其相应多辛离散守恒律,证明了它等价于通常的Fourier拟谱格式．数值结果表明,格式对于长时间计算具有稳定性与高精度．     

GSDBM方程的多辛几何和多辛Preissman格式研究（英文）

The multi-symplectic geometry for the GSDBM equation is presented in this paper. The multi-symplectic formulations for the GSDBM equation are presented and the local conservation laws are shown to correspond to certain well-known Hamiltonian functionals. The multi-symplectic discretization of each formulation is exemplified by the multisymplectic Preissmann scheme. The numerical experiments are given, and the results verify the efficiency of the Preissmann scheme.     

Multisymplecticity of the centred box discretization for hamiltonian PDEs with m ≥ 2 space dimensions

In this letter, it is shown that the centred box discretization for Hamiltonian PDEs with m ≥ 2 space dimensions is multisymplectic in the sense of Bridges and Reich in [1–6]. Multisymplectic discretizations for the generalized KP equation and the wave equation with 2 space dimensions, respectively, are given. A multisymplectically numerical scheme of the wave equation is derived.     

Synchronization of hyperchaotic oscillators via single unidirectional chaotic-coupling

In this paper, synchronization of two hyperchaotic oscillators via a single variable’s unidirectional coupling is studied. First, the synchronizability of the coupled hyperchaotic oscillators is proved mathematically. Then, the convergence speed of this synchronization scheme is analyzed. In order to speed up the response with a relatively large coupling strength, two kinds of chaotic coupling synchronization schemes are proposed. In terms of numerical simulations and the numerical calculation of the largest conditional Lyapunov exponent, it is shown that in a given range of coupling strengths, chaotic-coupling synchronization is quicker than the typical continuous-coupling synchronization. Furthermore, A circuit realization based on the chaotic synchronization scheme is designed and Pspice circuit simulation validates the simulated hyperchaos synchronization mechanism.     

Asymptotic‐preserving Godunov‐type numerical schemes for hyperbolic systems with stiff and nonstiff relaxation terms

We devise a new class of asymptotic‐preserving Godunov‐type numerical schemes for hyperbolic systems with stiff and nonstiff relaxation source terms governed by a relaxation time ε. As an alternative to classical operator‐splitting techniques, the objectives of these schemes are twofold: first, to give accurate numerical solutions for large, small, and in‐between values of ε and second, to make optional the choice of the numerical scheme in the asymptotic regime ε tends to zero. The latter property may be of particular interest to make easier and more efficient the coupling at a fixed spatial interface of two models involving very different values of ε. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

MULTISYMPLECTIC COMPOSITION INTEGRATORS OF HIGH ORDER

A composition method for constructing high order multisymplectic integrators is presented in this paper. The basic idea is to apply composition method to both the time and the space directions. We also obtain a general formula for composition method.     

二维抛物型方程的高精度多重网格解法

提出了数值求解二维抛物型方程的一种新的高精度加权平均紧隐格式，利用Fourier分析方法证明了该格式是无条件稳定的，为了克服传统迭代法在求解隐格式是收敛速度慢的缺陷，利用了多重网格加速技术，大大加快了迭代收敛速度，提高了求解效率，数值实验结果验证了方法的精确性和可靠性。