Functionally Fitted Energy-Preserving Methods for Solving Oscillatory Nonlinear Hamiltonian Systems

In the last few decades, numerical simulation for nonlinear oscillators has received a great deal of attention, and manyresearchers have been concerned with the design and analysis ofnumerical methods for solving oscillatory problems. In this paper, from the perspective of the continuous finite element method, we propose and analyze new energy-preserving functionally fitted methods, in particular trigonometrically fitted methods of an arbitrarily high order for solving oscillatory nonlinear Hamiltonian systems with a fixed frequency. To implement these new methods in a widespread way, they are transformed into a class of continuous-stage Runge--Kutta methods. This paper is accompanied by numerical experiments on oscillatory Hamiltonian systems such as the FPU problem and nonlinear Schrodinger equation. The numerical results demonstrate the remarkable accuracy and efficiency of our new methods compared with the existing high-order energy-preserving methods in the literature.

[1]  L. Brugnano,et al.  A simple framework for the derivation and analysis of effective one-step methods for ODEs , 2010, Appl. Math. Comput..

[2]  B. Hulme.  One-step piecewise polynomial Galerkin methods for initial value problems , 1972 .

[3]  J. M. Franco.  Exponentially fitted symplectic integrators of RKN type for solving oscillatory problems , 2007, Comput. Phys. Commun..

[4]  Wensheng Tang,et al.  Time finite element methods: A unified framework for numerical discretizations of ODEs , 2012, Appl. Math. Comput..

[5]  H. De Meyer,et al.  Exponentially fitted Runge-Kutta methods , 2000 .

[6]  D. H. Peregrine,et al.  Water waves, nonlinear Schrödinger equations and their solutions , 1983, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[7]  G. Quispel,et al.  Geometric integration using discrete gradients , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[8]  G. Quispel,et al.  Energy-preserving Runge-Kutta methods , 2009 .

[9]  Yuto Miyatake,et al.  An energy-preserving exponentially-fitted continuous stage Runge–Kutta method for Hamiltonian systems , 2014 .

[10]  Manuel Calvo,et al.  Symmetric and symplectic exponentially fitted Runge-Kutta methods of high order , 2010, Comput. Phys. Commun..

[11]  P. Betsch,et al.  Inherently Energy Conserving Time Finite Elements for Classical Mechanics , 2000 .

[12]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[13]  Bin Wang,et al.  Arbitrary-Order Trigonometric Fourier Collocation Methods for Multi-Frequency Oscillatory Systems , 2016, Found. Comput. Math..

[14]  Hans Van de Vyver.  A fourth-order symplectic exponentially fitted integrator , 2006, Comput. Phys. Commun..

[15]  E. Hairer.  Energy-preserving variant of collocation methods 1 , 2010 .

[16]  Kazufumi Ozawa,et al.  A functional fitting Runge-Kutta method with variable coefficients , 2001 .

[17]  Jianlin Xia,et al.  Explicit symplectic multidimensional exponential fitting modified Runge-Kutta-Nyström methods , 2012 .

[18]  Liviu Gr. Ixaru,et al.  P-stability and exponential-fitting methods for y″″ = f(x, y) , 1996 .

[19]  Ernst Hairer,et al.  Long-Time Energy Conservation of Numerical Methods for Oscillatory Differential Equations , 2000, SIAM J. Numer. Anal..

[20]  Ernst Hairer,et al.  Variable time step integration with symplectic methods , 1997 .

[21]  D. G. Bettis.  Numerical integration of products of fourier and ordinary polynomials , 1970 .

[22]  Nguyen Huu Cong,et al.  On functionally-fitted Runge–Kutta methods , 2006 .

[23]  Roger Grimshaw,et al.  Water Waves , 2021, Mathematics of Wave Propagation.

[24]  T. E. Simos,et al.  An exponentially-fitted Runge-Kutta method for the numerical integration of initial-value problems with periodic or oscillating solutions , 1998 .

[25]  Linda R. Petzold,et al.  Numerical solution of highly oscillatory ordinary differential equations , 1997, Acta Numerica.

[26]  Manuel Calvo,et al.  Structure preservation of exponentially fitted Runge-Kutta methods , 2008 .

[27]  Elena Celledoni,et al.  Preserving energy resp. dissipation in numerical PDEs using the "Average Vector Field" method , 2012, J. Comput. Phys..

[28]  O. Gonzalez.  Time integration and discrete Hamiltonian systems , 1996 .

[29]  W. Gautschi.  Numerical integration of ordinary differential equations based on trigonometric polynomials , 1961 .

[30]  Robert D. Skeel,et al.  Does variable step size ruin a symplectic integrator , 1992 .

[31]  C. Bottasso.  A new look at finite elements in time: a variational interpretation of Runge-Kutta methods , 1997 .

[32]  Elena Celledoni,et al.  Energy-Preserving Integrators and the Structure of B-series , 2010, Found. Comput. Math..

[33]  L. Brugnano,et al.  Hamiltonian Boundary Value Methods ( Energy Preserving Discrete Line Integral Methods ) 1 2 , 2009 .

[34]  Donald A. French,et al.  Continuous finite element methods which preserve energy properties for nonlinear problems , 1990 .

[35]  Jian,et al.  MULTISYMPLECTIC FOURIER PSEUDOSPECTRAL METHOD FOR THE NONLINEAR SCHRODINGER EQUATIONS WITH WAVE OPERATOR , 2007 .

[36]  J. C. Simo,et al.  Assessment of Energy-momentum and Symplectic Schemes for Stiff Dynamical Systems , 2022 .

[37]  Donato Trigiante,et al.  Energy- and Quadratic Invariants-Preserving Integrators Based upon Gauss Collocation Formulae , 2012, SIAM J. Numer. Anal..

[38]  Yuto Miyatake.  A derivation of energy-preserving exponentially-fitted integrators for Poisson systems , 2015, Comput. Phys. Commun..

[39]  Manuel Calvo,et al.  On high order symmetric and symplectic trigonometrically fitted Runge-Kutta methods with an even number of stages , 2010 .

[40]  Xinyuan Wu,et al.  Extended RKN-type methods for numerical integration of perturbed oscillators , 2009, Comput. Phys. Commun..