# American Institute of Mathematical Sciences

July  2021, 41(7): 3319-3341. doi: 10.3934/dcds.2020407

## Asymptotic analysis of a structure-preserving integrator for damped Hamiltonian systems

 1 Babeş-Bolyai University Cluj-Napoca, Str. Kogalniceanu 3, 400084 Cluj-Napoca, Romania 2 Romanian Institute of Science and Technology, 400022 Cluj-Napoca, Romania and, Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy, 400320 Cluj-Napoca, Romania 3 Institute for Numerical and Applied Mathematics, University of Göttingen, Lotzestr. 16-18, 37083 Göttingen, Germany

* Corresponding author: A. Viorel

Received  May 2020 Revised  September 2020 Published  December 2020

The present work deals with the numerical long-time integration of damped Hamiltonian systems. The method that we analyze combines a specific Strang splitting, that separates linear dissipative effects from conservative ones, with an energy-preserving averaged vector field (AVF) integrator for the Hamiltonian subproblem. This construction faithfully reproduces the energy-dissipation structure of the continuous model, its equilibrium points and its natural Lyapunov function. As a consequence of these structural similarities, both the convergence to equilibrium and, more interestingly, the energy decay rate of the continuous dynamical system are recovered at a discrete level. The possibility of replacing the implicit AVF integrator by an explicit Störmer-Verlet one is also discussed, while numerical experiments illustrate and support the theoretical findings.

Citation: Adrian Viorel, Cristian D. Alecsa, Titus O. Pinţa. Asymptotic analysis of a structure-preserving integrator for damped Hamiltonian systems. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3319-3341. doi: 10.3934/dcds.2020407
##### References:

show all references

##### References:
The (eventual) exponential decay of the potential energy (a) as well as state-spaces trajectories (b) for decreasing step sizes $h = 1 \text{ and } 0.1$, both with identical initial conditions $u_0 = 0.01, v_0 = 0$. The double-well potential is depicted in c)
Capturing an energy plateau: a) the AVF Splitting algorithm (black) compared to the trapezoidal rule (blue) for $h = 1$ (and benchmark Runge-Kutta (red)); b) AVF Splitting algorithm (black) compared to the conformal symplectic algorithm (28) (green) for $h = 0.1$. The contour lines of the nonconvex potential (33) are depicted in c)
Erroneous total energy oscillations: a) AVF Splitting algorithm (black) compared to the trapezoidal rule (blue) for $h = 0.1$ (and benchmark Runge-Kutta (red)); b) AVF Splitting algorithm (black) compared to the conformal symplectic algorithm (28) (green) for $h = 0.01$. The contour lines of the Rosenbrock potential are depicted in c)
 [1] Makoto Okumura, Daisuke Furihata. A structure-preserving scheme for the Allen–Cahn equation with a dynamic boundary condition. Discrete & Continuous Dynamical Systems, 2020, 40 (8) : 4927-4960. doi: 10.3934/dcds.2020206 [2] Alain Haraux. Some applications of the Łojasiewicz gradient inequality. Communications on Pure & Applied Analysis, 2012, 11 (6) : 2417-2427. doi: 10.3934/cpaa.2012.11.2417 [3] Qi Hong, Jialing Wang, Yuezheng Gong. Second-order linear structure-preserving modified finite volume schemes for the regularized long wave equation. Discrete & Continuous Dynamical Systems - B, 2019, 24 (12) : 6445-6464. doi: 10.3934/dcdsb.2019146 [4] Lin Lu, Qi Wang, Yongzhong Song, Yushun Wang. Local structure-preserving algorithms for the molecular beam epitaxy model with slope selection. Discrete & Continuous Dynamical Systems - B, 2021, 26 (9) : 4745-4765. doi: 10.3934/dcdsb.2020311 [5] Yuto Miyatake, Tai Nakagawa, Tomohiro Sogabe, Shao-Liang Zhang. A structure-preserving Fourier pseudo-spectral linearly implicit scheme for the space-fractional nonlinear Schrödinger equation. Journal of Computational Dynamics, 2019, 6 (2) : 361-383. doi: 10.3934/jcd.2019018 [6] Takeshi Fukao, Shuji Yoshikawa, Saori Wada. Structure-preserving finite difference schemes for the Cahn-Hilliard equation with dynamic boundary conditions in the one-dimensional case. Communications on Pure & Applied Analysis, 2017, 16 (5) : 1915-1938. doi: 10.3934/cpaa.2017093 [7] Zhuchun Li, Yi Liu, Xiaoping Xue. Convergence and stability of generalized gradient systems by Łojasiewicz inequality with application in continuum Kuramoto model. Discrete & Continuous Dynamical Systems, 2019, 39 (1) : 345-367. doi: 10.3934/dcds.2019014 [8] Jann-Long Chern, Zhi-You Chen, Yong-Li Tang. Structure of solutions to a singular Liouville system arising from modeling dissipative stationary plasmas. Discrete & Continuous Dynamical Systems, 2013, 33 (6) : 2299-2318. doi: 10.3934/dcds.2013.33.2299 [9] Vincent Giovangigli, Lionel Matuszewski. Structure of entropies in dissipative multicomponent fluids. Kinetic & Related Models, 2013, 6 (2) : 373-406. doi: 10.3934/krm.2013.6.373 [10] Canghua Jiang, Zhiqiang Guo, Xin Li, Hai Wang, Ming Yu. An efficient adjoint computational method based on lifted IRK integrator and exact penalty function for optimal control problems involving continuous inequality constraints. Discrete & Continuous Dynamical Systems - S, 2020, 13 (6) : 1845-1865. doi: 10.3934/dcdss.2020109 [11] Toshiko Ogiwara, Danielle Hilhorst, Hiroshi Matano. Convergence and structure theorems for order-preserving dynamical systems with mass conservation. Discrete & Continuous Dynamical Systems, 2020, 40 (6) : 3883-3907. doi: 10.3934/dcds.2020129 [12] Hideo Kubo. Asymptotic behavior of solutions to semilinear wave equations with dissipative structure. Conference Publications, 2007, 2007 (Special) : 602-613. doi: 10.3934/proc.2007.2007.602 [13] Jinying Ma, Honglei Xu. Empirical analysis and optimization of capital structure adjustment. Journal of Industrial & Management Optimization, 2020, 16 (3) : 1037-1047. doi: 10.3934/jimo.2018191 [14] Hong Zhang, Fei Yang. Optimization of capital structure in real estate enterprises. Journal of Industrial & Management Optimization, 2015, 11 (3) : 969-983. doi: 10.3934/jimo.2015.11.969 [15] Ahmet Sahiner, Gulden Kapusuz, Nurullah Yilmaz. A new smoothing approach to exact penalty functions for inequality constrained optimization problems. Numerical Algebra, Control & Optimization, 2016, 6 (2) : 161-173. doi: 10.3934/naco.2016006 [16] Changjun Yu, Kok Lay Teo, Liansheng Zhang, Yanqin Bai. A new exact penalty function method for continuous inequality constrained optimization problems. Journal of Industrial & Management Optimization, 2010, 6 (4) : 895-910. doi: 10.3934/jimo.2010.6.895 [17] Tetsuya Ishiwata, Kota Kumazaki. Structure preserving finite difference scheme for the Landau-Lifshitz equation with applied magnetic field. Conference Publications, 2015, 2015 (special) : 644-651. doi: 10.3934/proc.2015.0644 [18] Marius Durea, Elena-Andreea Florea, Radu Strugariu. Henig proper efficiency in vector optimization with variable ordering structure. Journal of Industrial & Management Optimization, 2019, 15 (2) : 791-815. doi: 10.3934/jimo.2018071 [19] Guangmei Shao, Wei Xue, Gaohang Yu, Xiao Zheng. Improved SVRG for finite sum structure optimization with application to binary classification. Journal of Industrial & Management Optimization, 2020, 16 (5) : 2253-2266. doi: 10.3934/jimo.2019052 [20] Jing-An Cui, Xinyu Song. Permanence of predator-prey system with stage structure. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 547-554. doi: 10.3934/dcdsb.2004.4.547

2019 Impact Factor: 1.338