In this thesis, we investigate the numerical solution of time-dependent nonlinear Schrödinger equations (more specifically, the Gross-Pitaevskii equation) that appear in the modeling of Bose-Einstein condensates. Since the model is known to conserve important physical invariants, such as mass and energy of the condensate, our goal is to study the importance of reproducing the conservation on the discrete level. The reliability of conservative, compared to non-conservative, methods shall be studied through high order collocation methods for the time discretization and finite element-based space discretizations. In particular, this includes symplectic discontinuous Galerkin time-stepping methods, as well as Continuous Petrov-Galerkin methods. The methods shall be tested for a problem with a known analytical solution, namely two interacting solitons in 1D. This problem is a suitable choice due to its high sensitivity to oscillations of the energy and difficulty to approximate for large time scales.
| Identifer | oai:union.ndltd.org:UPSALLA1/oai:DiVA.org:su-194703 |
| Date | January 2021 |
| Creators | Riera, Pau |
| Publisher | Stockholms universitet, Fysikum, KTH Royal Institute of Technology |
| Source Sets | DiVA Archive at Upsalla University |
| Language | English |
| Detected Language | English |
| Type | Student thesis, info:eu-repo/semantics/bachelorThesis, text |
| Format | application/pdf |
| Rights | info:eu-repo/semantics/openAccess |
Page generated in 0.0017 seconds