Finite Difference Models of Dynamical Systems with Quadratic Right-Hand Side
-
Published:2024-01-04
Issue:1
Volume:12
Page:167
-
ISSN:2227-7390
-
Container-title:Mathematics
-
language:en
-
Short-container-title:Mathematics
Author:
Malykh Mikhail12ORCID, Gambaryan Mark1ORCID, Kroytor Oleg1ORCID, Zorin Alexander1ORCID
Affiliation:
1. Department of Computational Mathematics and Artificial Intelligence, RUDN University, 117198 Moscow, Russia 2. Meshcheryakov Laboratory of Information Technologies, Joint Institute for Nuclear Research, 141980 Dubna, Russia
Abstract
Difference schemes that approximate dynamic systems are considered discrete models of the same phenomena that are described by continuous dynamic systems. Difference schemes with t-symmetry and midpoint and trapezoid schemes are considered. It is shown that these schemes are dual to each other, and, from this fact, we derive theorems on the inheritance of quadratic integrals by these schemes (Cooper’s theorem and its dual theorem on the trapezoidal scheme). Using examples of nonlinear oscillators, it is shown that these schemes poses challenges for theoretical research and practical application due to the problem of extra roots: these schemes do not allow one to unambiguously determine the final values from the initial values and vice versa. Therefore, we consider difference schemes in which the transitions from layer to layer in time are carried out using birational transformations (Cremona transformations). Such schemes are called reversible. It is shown that reversible schemes with t-symmetry can be easily constructed for any dynamical system with a quadratic right-hand side. As an example of such a dynamic system, a top fixed at its center of gravity is considered in detail. In this case, the discrete theory repeats the continuous theory completely: (1) the points of the approximate solution lie on some elliptic curve, which at Δt→0 turns into an integral curve; (2) the difference scheme can be represented using quadrature; and (3) the approximate solution can be represented using an elliptic function of a discrete argument. The last section considers the general case. The integral curves are replaced with closures of the orbits of the corresponding Cremona transformation as sets in the projective space over R. The problem of the dimension of this set is discussed.
Funder
Russian Science Foundation RUDN University Strategic Academic Leadership Program
Subject
General Mathematics,Engineering (miscellaneous),Computer Science (miscellaneous)
Reference40 articles.
1. Hairer, E., Wanner, G., and Nørsett, S.P. (2008). Solving Ordinary Differential Equations I, Springer. [3rd ed.]. 2. Hairer, E., Wanner, G., and Lubich, C. (2000). Geometric Numerical Integration. Structure-Preserving Algorithms for Ordinary Differential Equations, Springer. 3. Brugnano, L., and Iavernaro, F. (2016). Line Integral Methods for Conservative Problems, CRC Press. 4. da Veiga, L.B., Lipnikov, K., and Manzini, G. (2014). The Mimetic Finite Difference Method for Elliptic Problems, Springer. 5. Castillo, J.E., and Miranda, G.F. (2013). Mimetic Discretization Methods, CRC.
|
|