We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-or...We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.展开更多
we investigate the structure-preserving numerical algorithm of the Degasperis Procesi equation which can be transformed into a bi-Hamiltonian form using the discrete variational derivative method. Based on two differe...we investigate the structure-preserving numerical algorithm of the Degasperis Procesi equation which can be transformed into a bi-Hamiltonian form using the discrete variational derivative method. Based on two different space discretization methods, the Fourier pseudospectral method and the wavelet collocation method, we develop two modified structure-preserving schemes under the periodic boundary condition. These proposed schemes are proved to preserve the Hamiltonian invariants theoretically and numerically. Meanwhile, the numerical results confirm that they can simulate the propagation of solitons effectively for a long time.展开更多
The two-component Camassa–Holm equation includes many intriguing phenomena. We propose a multi-symplectic compact method to solve the two-component Camassa–Holm equation. Based on its multi-symplectic formulation, t...The two-component Camassa–Holm equation includes many intriguing phenomena. We propose a multi-symplectic compact method to solve the two-component Camassa–Holm equation. Based on its multi-symplectic formulation, the proposed method is derived by the sixth-order compact finite difference method in spatial discretization and the symplectic implicit midpoint scheme in temporal discretization. Numerical experiments finely describe the velocity and density variables in the two-component integrable system and distinctly display the evolvement of the singular solutions. Moreover, the proposed method shows good conservative properties during long-time numerical simulation.展开更多
基金Supported by the National Natural Science Foundation of China under Grant Nos 11571366 and 11501570the Open Foundation of State Key Laboratory of High Performance Computing of China+1 种基金the Research Fund of National University of Defense Technology under Grant No JC15-02-02the Fund from HPCL
文摘We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11571366 and 11501570the Open Foundation of State Key Laboratory of High Performance Computing of China
文摘we investigate the structure-preserving numerical algorithm of the Degasperis Procesi equation which can be transformed into a bi-Hamiltonian form using the discrete variational derivative method. Based on two different space discretization methods, the Fourier pseudospectral method and the wavelet collocation method, we develop two modified structure-preserving schemes under the periodic boundary condition. These proposed schemes are proved to preserve the Hamiltonian invariants theoretically and numerically. Meanwhile, the numerical results confirm that they can simulate the propagation of solitons effectively for a long time.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11571366 and 11501570the Open Foundation of State Key Laboratory of High Performance Computing of China+1 种基金the Research Fund of the National University of Defense Technology under Grant No JC15-02-02the Fund from HPCL
文摘The two-component Camassa–Holm equation includes many intriguing phenomena. We propose a multi-symplectic compact method to solve the two-component Camassa–Holm equation. Based on its multi-symplectic formulation, the proposed method is derived by the sixth-order compact finite difference method in spatial discretization and the symplectic implicit midpoint scheme in temporal discretization. Numerical experiments finely describe the velocity and density variables in the two-component integrable system and distinctly display the evolvement of the singular solutions. Moreover, the proposed method shows good conservative properties during long-time numerical simulation.