Structure-preserving reduced-order modeling of Korteweg–de Vries equation

Date

2021-10-01

Author

Uzunca, Murat
Karasözen, Bülent
Yıldız, Süleyman

Metadata

Show full item record

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Item Usage Stats

500
views

0
downloads

Computationally efficient, structure-preserving reduced-order methods are developed for the Korteweg–de Vries (KdV) equations in Hamiltonian form. The semi-discretization in space by finite differences is based on the Hamiltonian structure. The resulting skew-gradient system of ordinary differential equations (ODEs) is integrated with the linearly implicit Kahan's method, which preserves the Hamiltonian approximately. We have shown, using proper orthogonal decomposition (POD), the Hamiltonian structure of the full-order model (FOM) is preserved by the reduced-order model (ROM). The reduced model has the same linear–quadratic structure as the FOM. The quadratic nonlinear terms of the KdV equations are evaluated efficiently by the use of tensorial framework, clearly separating the offline–online cost of the FOMs and ROMs. The accuracy of the reduced solutions, preservation of the conserved quantities, and computational speed-up gained by ROMs are demonstrated for the one-dimensional single and coupled KdV equations, and two-dimensional Zakharov–Kuznetsov equation with soliton solutions.

Subject Keywords

Hamiltonian systems, Solitary waves, Kahan’s method, Energy preservation, Model order reduction, Tensor algebra

URI

https://www.scopus.com/inward/record.uri?partnerID=HzOxMe3b&scp=85104447301&origin=inward
https://hdl.handle.net/11511/90331

Journal

Mathematics and Computers in Simulation

DOI

https://doi.org/10.1016/j.matcom.2021.03.042

Collections

Graduate School of Applied Mathematics, Article

Suggestions

OpenMETU
Core

Energy preserving methods for lattice equations Erdem, Özge; Karasözen, Bülent (2010-11-27) Integral preserving methods, like the averaged vector field, discrete gradient and trapezoidal methods are to Poisson systems. Numerical experiments on the Volterra equations and integrable discretization of the nonlinear Schrodinger equation are presented.
Symplectic and multisymplectic Lobatto methods for the "good" Boussinesq equation AYDIN, AYHAN; Karasözen, Bülent (2008-08-01) In this paper, we construct second order symplectic and multisymplectic integrators for the "good" Boussineq equation using the two-stage Lobatto IIIA-IIIB partitioned Runge-Kutta method, which yield an explicit scheme and is equivalent to the classical central difference approximation to the second order spatial derivative. Numerical dispersion properties and the stability of both integrators are investigated. Numerical results for different solitary wave solutions confirm the excellent long time behavior ...
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.
Multi-symplectic integration of coupled non-linear Schrodinger system with soliton solutions AYDIN, AYHAN; Karasözen, Bülent (2009-01-01) Systems of coupled non-linear Schrodinger equations with soliton solutions are integrated using the six-point scheme which is equivalent to the multi-symplectic Preissman scheme. The numerical dispersion relations are studied for the linearized equation. Numerical results for elastic and inelastic soliton collisions are presented. Numerical experiments confirm the excellent conservation of energy, momentum and norm in long-term computations and their relations to the qualitative behaviour of the soliton sol...
Series solution of the wave equation in optic fiber Çıldır, Sema; Çakır, Serhat; Department of Physics (2007) In this study, the mapped Galerkin method was applied to solve the vector wave equation based on H field and to obtain the propagation constant in x y space. The vector wave equation was solved by the transformation of the infinite x y plane onto a unit square. Two-dimensional Fourier series expansions were used in the solutions. Modal fields and propagation constants of dielectric waveguides were calculated. In the first part of the study, all of the calculations were made in step index fibers. Transvers...

Citation Formats

M. Uzunca, B. Karasözen, and S. Yıldız, “Structure-preserving reduced-order modeling of Korteweg–de Vries equation,” Mathematics and Computers in Simulation, pp. 193–211, 2021, Accessed: 00, 2021. [Online]. Available: https://www.scopus.com/inward/record.uri?partnerID=HzOxMe3b&scp=85104447301&origin=inward.