We investigate the multisymplectic Euler box scheme for the Korteweg-de Vries (KdV) equation. A new completely explicit six-point scheme is derived. Numerical experiments of the new scheme with comparisons to the Za...We investigate the multisymplectic Euler box scheme for the Korteweg-de Vries (KdV) equation. A new completely explicit six-point scheme is derived. Numerical experiments of the new scheme with comparisons to the Zabusky-Kruskal scheme, the multisymplectic 12-point scheme, the narrow box scheme and the spectral method are made to show nice numerical stability and ability to preserve the integral invariant for long-time integration.展开更多
The Hamiltonian formulations of the linear 'good' Boussinesq (L.G.B.) equation and the multi-symplectic formulation of the nonlinear 'good' Boussinesq (N.G.B.) equation are considered. For the multi-sy...The Hamiltonian formulations of the linear 'good' Boussinesq (L.G.B.) equation and the multi-symplectic formulation of the nonlinear 'good' Boussinesq (N.G.B.) equation are considered. For the multi-symplectic formulation, a new fifteen-point difference scheme which is equivalent to the multi-symplectic Preissmann integrator is derived. We also present numerical experiments, which show that the symplectic and multi-symplectic schemes have excellent long-time numerical behavior.展开更多
This paper considers the multi-symplectic formulations of the generalized fifth-order KdV equation in Hamiltonian space. Recurring to the midpoint rule, it presents an implicit multi-symplectic scheme with discrete mu...This paper considers the multi-symplectic formulations of the generalized fifth-order KdV equation in Hamiltonian space. Recurring to the midpoint rule, it presents an implicit multi-symplectic scheme with discrete multi-symplectic conservation law to solve the partial differential equations which are derived from the generalized fifth-order KdV equation numerically. The results of the numerical experiments show that this multi-symplectic algorithm is good in accuracy and its long-time numerical behaviour is also perfect.展开更多
The dynamic analysis on the ultra-large spatial structure can be simplified drastically by ignoring the flexibility and damping of the structure.However,these simplifications will result in the erroneous estimate on t...The dynamic analysis on the ultra-large spatial structure can be simplified drastically by ignoring the flexibility and damping of the structure.However,these simplifications will result in the erroneous estimate on the dynamic behaviors of the ultra-large spatial structure.Taking the spatial beam as an example,the minimum control energy defined by the difference between the initial total energy and the final total energy in the assumed stable attitude state of the beam is investigated by the structure-preserving method proposed in our previous studies in two cases:the spatial beam considering the flexibility as well as the damping effect,and the spatial beam ignoring both the flexibility and the damping effect.In the numerical experiments,the assumed simulation interval of three months is evaluated on whether or not it is long enough for the spatial flexible damping beam to arrive at the assumed stable attitude state.And then,taking the initial attitude angle and the initial attitude angle velocity as the independent variables,respectively,the minimum control energies of the mentioned two cases are investigated in detail.From the numerical results,the following conclusions can be obtained.With the fixed initial attitude angle velocity,the minimum control energy of the spatial flexible damping beam is higher than that of the spatial rigid beam when the initial attitude angle is close to or far away from the stable attitude state.With the fixed initial attitude angle,ignoring the flexibility and the damping effect will underestimate the minimum control energy of the spatial beam.展开更多
Local structure-preserving algorithms including multi-symplectic, local energy- and momentum-preserving schemes are proposed for the generalized Rosenau-RLW-KdV equation based on the multi-symplectic Hamiltonian formu...Local structure-preserving algorithms including multi-symplectic, local energy- and momentum-preserving schemes are proposed for the generalized Rosenau-RLW-KdV equation based on the multi-symplectic Hamiltonian formula of the equation. Each of the present algorithms holds a discrete conservation law in any time-space region. For the original problem subjected to appropriate boundary conditions, these algorithms will be globally conservative. Discrete fast Fourier transform makes a significant improvement to the computational efficiency of schemes. Numerical results show that the proposed algorithms have satisfactory performance in providing an accurate solution and preserving the discrete invariants.展开更多
The mixed state of two-band II-superconductor is analyzed by the multi-symplectic method. As to the Ginzburg-Landau equation depending on time that describes the mixed state of two-band II-superconductor, the multi-sy...The mixed state of two-band II-superconductor is analyzed by the multi-symplectic method. As to the Ginzburg-Landau equation depending on time that describes the mixed state of two-band II-superconductor, the multi-symplectic formulations with several conservation laws: a multi-symplectic conservation law, an energy con- servation law, as well as a momentum conservation law, are presented firstly; then an eighteen points scheme is constructed to simulate the multi-symplectic partial differential equations (PDEs) that are derived from the Ginzburg-Landau equation; finally, based on the simulation results, the volt-ampere characteristic curves are obtained, as well as the relationships between the temperature and resistivity of a suppositional two-band II-superconductor model under different magnetic intensi- ties. From the results of the numerical experiments, it is concluded that the notable property of the mixed state of the two-band II-superconductor is that: The trans- formation temperature decreases and the resistivity ρ increases rapidly with the increase of the magnetic intensity B. In addition, the simulation results show that the multi-symplectic method has two remarkable advantages: high accuracy and excellent long-time numerical behavior.展开更多
Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and ...Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and the corresponding multi-symplectic conservation property is proved. The backward error analysis shows that the explicit multi-symplectic scheme has good accuracy. The sine-Gordon equation and the Klein-Gordon equation are simulated by an explicit multi-symplectic scheme. The numerical results show that the new explicit multi-symplectic scheme can well simulate the solitary wave behaviors of the nonlinear wave equation and approximately preserve the relative energy error of the equation.展开更多
In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-...In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws-a multi-symplectic conservation law (CLS), a local energy conservation law (ECL) as well as a local momentum conservation law (MCL) --is constructed to discrete the PDEs that are derived from the membrane free vibration equation. The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior,展开更多
We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler ...We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler method for temporal discretization. It is verified that the proposed method has corresponding discrete multi-symplectic conservation laws. Numerical simulations indicate that the proposed scheme is characterized by excellent conservation.展开更多
In this paper we propose stochastic multi-symplectic conservation law for stochastic Hamiltonian partial differential equations,and develop a stochastic multisymplectic method for numerically solving a kind of stochas...In this paper we propose stochastic multi-symplectic conservation law for stochastic Hamiltonian partial differential equations,and develop a stochastic multisymplectic method for numerically solving a kind of stochastic nonlinear Schrodinger equations.It is shown that the stochasticmulti-symplecticmethod preserves themultisymplectic structure,the discrete charge conservation law,and deduces the recurrence relation of the discrete energy.Numerical experiments are performed to verify the good behaviors of the stochastic multi-symplectic method in cases of both solitary wave and collision.展开更多
We propose a multi-symplectic wavelet splitting equations. Based on its mu]ti-symplectic formulation, method to solve the strongly coupled nonlinear SchrSdinger the strongly coupled nonlinear SchrSdinger equations can...We propose a multi-symplectic wavelet splitting equations. Based on its mu]ti-symplectic formulation, method to solve the strongly coupled nonlinear SchrSdinger the strongly coupled nonlinear SchrSdinger equations can be split into one linear multi-symplectic subsystem and one nonlinear infinite-dimensional Hamiltonian subsystem. For the linear subsystem, the multi-symplectic wavelet collocation method and the symplectic Euler method are employed in spatial and temporal discretization, respectively. For the nonlinear subsystem, the mid-point symplectic scheme is used. Numerical simulations show the effectiveness of the proposed method during long-time numerical calculation.展开更多
In this paper, using the concatenating method, a series of local structure-preserving algorithms are obtained for the Klein-Gordon-Zakharov equation, including four multisymplectic algorithms, four local energy-preser...In this paper, using the concatenating method, a series of local structure-preserving algorithms are obtained for the Klein-Gordon-Zakharov equation, including four multisymplectic algorithms, four local energy-preserving algorithms, four local momentumpreserving algorithms;of these, local energy-preserving and momentum-preserving algorithms have not been studied before. The local structure-preserving algorithms mentioned above are more widely used than the global structure-preserving algorithms, since local preservation algorithms can be preserved in any time and space domains, which overcomes the defect that global preservation algorithms are limited to boundary conditions. In particular, under appropriate boundary conditions, local preservation laws are global preservation laws.Numerical experiments conducted can support the theoretical analysis well.展开更多
In this paper,the complex multi-symplectic method and the implementation of the generalized sinhGordon equation are investigated in detail.The multi-symplectic formulations of the generalized sinh-Gordon equation in H...In this paper,the complex multi-symplectic method and the implementation of the generalized sinhGordon equation are investigated in detail.The multi-symplectic formulations of the generalized sinh-Gordon equation in Hamiltonian space are presented firstly.The complex method is introduced and a complex semi-implicit scheme with several discrete conservation laws(including a multi-symplectic conservation law(CLS),a local energy conservation law(ECL) as well as a local momentum conservation law(MCL)) is constructed to solve the partial differential equations(PDEs) that are derived from the generalized sinh-Gordon equation numerically.The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior and high accuracy.展开更多
In this paper, based on the concatenating method, we present a unified framework to construct a series of local structure-preserving algorithms for the Korteweg-de Vries (KdV) equation, including eight multi-symplec...In this paper, based on the concatenating method, we present a unified framework to construct a series of local structure-preserving algorithms for the Korteweg-de Vries (KdV) equation, including eight multi-symplectic algorithms, eight local energy-conserving algo- rithms and eight local momentum-conserving algorithms. Among these algorithms, some have been discussed and widely used while the most are new. The outstanding advantage of these proposed algorithms is that they conserve the local structures in any time-space re- gion exactly. Therefore, the local structure-preserving algorithms overcome the restriction of global structure-preserving algorithms on the boundary conditions. Numerical experiments are conducted to show the performance of the proposed methods. Moreover, the unified framework can be easily applied to many other equations.展开更多
Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic ...Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations (PDEs) derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.展开更多
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, multi- symplectic formulations for the higher order wa...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, multi- symplectic 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 each formulation is calculated by the multi-symplectic Fourier pseudospectral scheme. Numerical experiments are carried out, which verify the efficiency of the Fourier pseudospectral method.展开更多
基金Supported by the National Baslc Research Programme under Grant No 2005CB321703, and the National Natural Science Foundation of China under Grant Nos 40221503, 10471067 and 40405019.
文摘We investigate the multisymplectic Euler box scheme for the Korteweg-de Vries (KdV) equation. A new completely explicit six-point scheme is derived. Numerical experiments of the new scheme with comparisons to the Zabusky-Kruskal scheme, the multisymplectic 12-point scheme, the narrow box scheme and the spectral method are made to show nice numerical stability and ability to preserve the integral invariant for long-time integration.
基金Supported by State Key Laboratory of Scientific/Engineering Computing Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences and Foundation of office of overseas Chinese affair under the state council (02
文摘The Hamiltonian formulations of the linear 'good' Boussinesq (L.G.B.) equation and the multi-symplectic formulation of the nonlinear 'good' Boussinesq (N.G.B.) equation are considered. For the multi-symplectic formulation, a new fifteen-point difference scheme which is equivalent to the multi-symplectic Preissmann integrator is derived. We also present numerical experiments, which show that the symplectic and multi-symplectic schemes have excellent long-time numerical behavior.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10572119, 10772147 and 10632030)the Doctoral Program Foundation of Education Ministry of China (Grant No 20070699028)+1 种基金the National Natural Science Foundation of Shaanxi Province of China (Grant No 2006A07)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment
文摘This paper considers the multi-symplectic formulations of the generalized fifth-order KdV equation in Hamiltonian space. Recurring to the midpoint rule, it presents an implicit multi-symplectic scheme with discrete multi-symplectic conservation law to solve the partial differential equations which are derived from the generalized fifth-order KdV equation numerically. The results of the numerical experiments show that this multi-symplectic algorithm is good in accuracy and its long-time numerical behaviour is also perfect.
基金The research is supported by the National Natural Science Foundation of China(11672241,11972284,11432010)Fund for Distinguished Young Scholars of Shaanxi Province(2019JC-29)Fund of the Youth Innovation Team of Shaanxi Universities,the Seed Foundation of Qian Xuesen Laboratory of Space Technology,and the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment(GZ1605).
文摘The dynamic analysis on the ultra-large spatial structure can be simplified drastically by ignoring the flexibility and damping of the structure.However,these simplifications will result in the erroneous estimate on the dynamic behaviors of the ultra-large spatial structure.Taking the spatial beam as an example,the minimum control energy defined by the difference between the initial total energy and the final total energy in the assumed stable attitude state of the beam is investigated by the structure-preserving method proposed in our previous studies in two cases:the spatial beam considering the flexibility as well as the damping effect,and the spatial beam ignoring both the flexibility and the damping effect.In the numerical experiments,the assumed simulation interval of three months is evaluated on whether or not it is long enough for the spatial flexible damping beam to arrive at the assumed stable attitude state.And then,taking the initial attitude angle and the initial attitude angle velocity as the independent variables,respectively,the minimum control energies of the mentioned two cases are investigated in detail.From the numerical results,the following conclusions can be obtained.With the fixed initial attitude angle velocity,the minimum control energy of the spatial flexible damping beam is higher than that of the spatial rigid beam when the initial attitude angle is close to or far away from the stable attitude state.With the fixed initial attitude angle,ignoring the flexibility and the damping effect will underestimate the minimum control energy of the spatial beam.
基金supported by the National Natural Science Foundation of China(Grant Nos.11201169 and 61672013)the Foundation of Jiangsu Key Laboratory for Numerical Simulation of Large Scale Complex Systems(Grant No.201606)
文摘Local structure-preserving algorithms including multi-symplectic, local energy- and momentum-preserving schemes are proposed for the generalized Rosenau-RLW-KdV equation based on the multi-symplectic Hamiltonian formula of the equation. Each of the present algorithms holds a discrete conservation law in any time-space region. For the original problem subjected to appropriate boundary conditions, these algorithms will be globally conservative. Discrete fast Fourier transform makes a significant improvement to the computational efficiency of schemes. Numerical results show that the proposed algorithms have satisfactory performance in providing an accurate solution and preserving the discrete invariants.
基金the National Natural Science Foundation of China (Grant Nos. 10572119, 10772147 and 10632030)the Doctoral Program Founda-tion of Education Ministry of China (Grant No. 20070699028)+1 种基金the Natural Science Foundation of Shaanxi Province of China (Grant No. 2006A07)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment
文摘The mixed state of two-band II-superconductor is analyzed by the multi-symplectic method. As to the Ginzburg-Landau equation depending on time that describes the mixed state of two-band II-superconductor, the multi-symplectic formulations with several conservation laws: a multi-symplectic conservation law, an energy con- servation law, as well as a momentum conservation law, are presented firstly; then an eighteen points scheme is constructed to simulate the multi-symplectic partial differential equations (PDEs) that are derived from the Ginzburg-Landau equation; finally, based on the simulation results, the volt-ampere characteristic curves are obtained, as well as the relationships between the temperature and resistivity of a suppositional two-band II-superconductor model under different magnetic intensi- ties. From the results of the numerical experiments, it is concluded that the notable property of the mixed state of the two-band II-superconductor is that: The trans- formation temperature decreases and the resistivity ρ increases rapidly with the increase of the magnetic intensity B. In addition, the simulation results show that the multi-symplectic method has two remarkable advantages: high accuracy and excellent long-time numerical behavior.
基金Project supported by the National Natural Science Foundation of China(Nos.11161017,11071251,and 10871099)the National Basic Research Program of China(973 Program)(No.2007CB209603)+1 种基金the Natural Science Foundation of Hainan Province(No.110002)the Scientific Research Foun-dation of Hainan University(No.kyqd1053)
文摘Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and the corresponding multi-symplectic conservation property is proved. The backward error analysis shows that the explicit multi-symplectic scheme has good accuracy. The sine-Gordon equation and the Klein-Gordon equation are simulated by an explicit multi-symplectic scheme. The numerical results show that the new explicit multi-symplectic scheme can well simulate the solitary wave behaviors of the nonlinear wave equation and approximately preserve the relative energy error of the equation.
基金the National Natural Science Foundation of China(Nos.10632030 and 10572119)Program for New Century Excellent Talents of Ministry of Education of China(No.NCET-04-0958)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment
文摘In this paper, the multi-symplectic formulations of the membrane free vibration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws-a multi-symplectic conservation law (CLS), a local energy conservation law (ECL) as well as a local momentum conservation law (MCL) --is constructed to discrete the PDEs that are derived from the membrane free vibration equation. The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior,
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 10971226, 91130013, and 11001270)the National Basic Research Program of China (Grant No. 2009CB723802)
文摘We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler method for temporal discretization. It is verified that the proposed method has corresponding discrete multi-symplectic conservation laws. Numerical simulations indicate that the proposed scheme is characterized by excellent conservation.
基金supported by the NNSFC(No.11001009)supported by the Director Foundation of GUCAS,the NNSFC(No.11071251)supported by the Foundation of CAS and the NNSFC(No.11021101,No.91130003).
文摘In this paper we propose stochastic multi-symplectic conservation law for stochastic Hamiltonian partial differential equations,and develop a stochastic multisymplectic method for numerically solving a kind of stochastic nonlinear Schrodinger equations.It is shown that the stochasticmulti-symplecticmethod preserves themultisymplectic structure,the discrete charge conservation law,and deduces the recurrence relation of the discrete energy.Numerical experiments are performed to verify the good behaviors of the stochastic multi-symplectic method in cases of both solitary wave and collision.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10971226,91130013,and 11001270)the National Basic Research Program of China(Grant No.2009CB723802)+1 种基金the Research Innovation Fund of Hunan Province,China (Grant No.CX2011B011)the Innovation Fund of National University of Defense Technology,China(Grant No.B120205)
文摘We propose a multi-symplectic wavelet splitting equations. Based on its mu]ti-symplectic formulation, method to solve the strongly coupled nonlinear SchrSdinger the strongly coupled nonlinear SchrSdinger equations can be split into one linear multi-symplectic subsystem and one nonlinear infinite-dimensional Hamiltonian subsystem. For the linear subsystem, the multi-symplectic wavelet collocation method and the symplectic Euler method are employed in spatial and temporal discretization, respectively. For the nonlinear subsystem, the mid-point symplectic scheme is used. Numerical simulations show the effectiveness of the proposed method during long-time numerical calculation.
基金supported by the National Natural Science Foundation of China(11801277,11771213,12171245)。
文摘In this paper, using the concatenating method, a series of local structure-preserving algorithms are obtained for the Klein-Gordon-Zakharov equation, including four multisymplectic algorithms, four local energy-preserving algorithms, four local momentumpreserving algorithms;of these, local energy-preserving and momentum-preserving algorithms have not been studied before. The local structure-preserving algorithms mentioned above are more widely used than the global structure-preserving algorithms, since local preservation algorithms can be preserved in any time and space domains, which overcomes the defect that global preservation algorithms are limited to boundary conditions. In particular, under appropriate boundary conditions, local preservation laws are global preservation laws.Numerical experiments conducted can support the theoretical analysis well.
基金Supported by the National Natural Science Foundation of China (Grant Nos.10572119,10772147,and 10632030)the Doctoral Program Foundation of Education Ministry of China (Grant No.20070699028)+2 种基金the China Postdoctoral Science Foundation (Grant No. 20090450170)the National Natural Science Foundation of Shaanxi Province of China (Grant No.2006A07)the NPU Foundation for Fundamental Research and the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment (Grant No.GZ0802)
文摘In this paper,the complex multi-symplectic method and the implementation of the generalized sinhGordon equation are investigated in detail.The multi-symplectic formulations of the generalized sinh-Gordon equation in Hamiltonian space are presented firstly.The complex method is introduced and a complex semi-implicit scheme with several discrete conservation laws(including a multi-symplectic conservation law(CLS),a local energy conservation law(ECL) as well as a local momentum conservation law(MCL)) is constructed to solve the partial differential equations(PDEs) that are derived from the generalized sinh-Gordon equation numerically.The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior and high accuracy.
文摘In this paper, based on the concatenating method, we present a unified framework to construct a series of local structure-preserving algorithms for the Korteweg-de Vries (KdV) equation, including eight multi-symplectic algorithms, eight local energy-conserving algo- rithms and eight local momentum-conserving algorithms. Among these algorithms, some have been discussed and widely used while the most are new. The outstanding advantage of these proposed algorithms is that they conserve the local structures in any time-space re- gion exactly. Therefore, the local structure-preserving algorithms overcome the restriction of global structure-preserving algorithms on the boundary conditions. Numerical experiments are conducted to show the performance of the proposed methods. Moreover, the unified framework can be easily applied to many other equations.
基金supported by the National Natural Science Foundation of China (Nos. 10772147 and10632030)the Ph. D. Program Foundation of Ministry of Education of China (No. 20070699028)+2 种基金the Natural Science Foundation of Shaanxi Province of China (No. 2006A07)the Open Foundationof State Key Laboratory of Structural Analysis of Industrial Equipment (No. GZ0802)the Foundation for Fundamental Research of Northwestern Polytechnical University
文摘Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations (PDEs) derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.
文摘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, multi- symplectic 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 each formulation is calculated by the multi-symplectic Fourier pseudospectral scheme. Numerical experiments are carried out, which verify the efficiency of the Fourier pseudospectral method.