Abstract

In this work, a novel class of high-order energy-preserving algorithms are developed for simulating the coupled Klein-Gordon-Schrödinger equations. We introduce a Lagrange multiplier approach and derive a new system which warrants exactly the original total energy, instead of the modified quadratic energy in the previous energy quadratization approach. The Lagrange system is then reformulated equivalently into an optimization problem subject to PDE constraints. We discretize the optimization system by applying the Gauss collocation method together with the prediction-correction technique in time and the sine pseudo-spectral method in space, which leads to a family of fully discrete high-order energy-preserving schemes. In the numerical experiments, we compare this method with other previous methods and verify the high accuracy as well as its ability to maintain the original energy at each time layer. The feasibility and validity of new method are also witnessed for the KGS model in the nonrelativistic limit regime. Additionally, enlightened by this seminal idea, another type of high-order numerical algorithms conserving two invariants of mass and energy are also presented and relevant experimental results are listed to verify the power of the suggested method.

Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.