Energy preserving model order reduction of the nonlinear Schrodinger equation

Download
2018-12-01
An energy preserving reduced order model is developed for two dimensional nonlinear Schrodinger equation (NLSE) with plane wave solutions and with an external potential. The NLSE is discretized in space by the symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting system of Hamiltonian ordinary differential equations are integrated in time by the energy preserving average vector field (AVF) method. The mass and energy preserving reduced order model (ROM) is constructed by proper orthogonal decomposition (POD) Galerkin projection. The nonlinearities are computed for the ROM efficiently by discrete empirical interpolation method (DEIM) and dynamic mode decomposition (DMD). Preservation of the semi-discrete energy and mass are shown for the full order model (FOM) and for the ROM which ensures the long term stability of the solutions. Numerical simulations illustrate the preservation of the energy and mass in the reduced order model for the two dimensional NLSE with and without the external potential. The POD-DMD makes a remarkable improvement in computational speed-up over the POD-DEIM. Both methods approximate accurately the FOM, whereas POD-DEIM is more accurate than the POD-DMD.
ADVANCES IN COMPUTATIONAL MATHEMATICS

Suggestions

Average Vector Field Splitting Method for Nonlinear Schrodinger Equation
Akkoyunlu, Canan; Karasözen, Bülent (2012-05-02)
The energy preserving average vector field integrator is applied to one and two dimensional Schrodinger equations with symmetric split-step method. The numerical results confirm the long-term preservation of the Hamiltonians, which is essential in simulating periodic waves.
Time-Space Adaptive Method of Time Layers for the Advective Allen-Cahn Equation
UZUNCA, MURAT; Karasözen, Bülent; Sariaydin-Filibelioglu, Ayse (2015-09-18)
We develop an adaptive method of time layers with a linearly implicit Rosenbrock method as time integrator and symmetric interior penalty Galerkin method for space discretization for the advective Allen-Cahn equation with nondivergence-free velocity fields. Numerical simulations for convection dominated problems demonstrate the accuracy and efficiency of the adaptive algorithm for resolving the sharp layers occurring in interface problems with small surface tension.
Structure preserving reduced order modeling for gradient systems
Akman Yildiz, Tugba; UZUNCA, MURAT; Karasözen, Bülent (2019-04-15)
Minimization of energy in gradient systems leads to formation of oscillatory and Turing patterns in reaction-diffusion systems. These patterns should be accurately computed using fine space and time meshes over long time horizons to reach the spatially inhomogeneous steady state. In this paper, a reduced order model (ROM) is developed which preserves the gradient dissipative structure. The coupled system of reaction-diffusion equations are discretized in space by the symmetric interior penalty discontinuous...
Energy Stable Discontinuous Galerkin Finite Element Method for the Allen-Cahn Equation
Karasözen, Bülent; Sariaydin-Filibelioglu, Ayse; Yücel, Hamdullah (2018-05-01)
In this paper, we investigate numerical solution of Allen-Cahn equation with constant and degenerate mobility, and with polynomial and logarithmic energy functionals. We discretize the model equation by symmetric interior penalty Galerkin (SIPG) method in space, and by average vector field (AVF) method in time. We show that the energy stable AVF method as the time integrator for gradient systems like the Allen-Cahn equation satisfies the energy decreasing property for fully discrete scheme. Numerical result...
Exact solution of Schrodinger equation for Pseudoharmonic potential
Sever, Ramazan; TEZCAN, CEVDET; Aktas, Metin; Yesiltas, Oezlem (2008-02-01)
Exact solution of Schrodinger equation for the pseudoharmonic potential is obtained for an arbitrary angular momentum. The energy eigenvalues and corresponding eigenfunctions are calculated by Nikiforov-Uvarov method. Wavefunctions are expressed in terms of Jacobi polynomials. The energy eigenvalues are calculated numerically for some values of l and n with n <= 5 for some diatomic molecules.
Citation Formats
B. Karasözen, “Energy preserving model order reduction of the nonlinear Schrodinger equation,” ADVANCES IN COMPUTATIONAL MATHEMATICS, pp. 1769–1796, 2018, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/32097.