skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: High-order symplectic Lie group methods on $ SO(n) $ using the polar decomposition
A variational integrator of arbitrarily high-order on the special orthogonal group \begin{document}$ SO(n) $$\end{document} is constructed using the polar decomposition and the constrained Galerkin method. It has the advantage of avoiding the second order derivative of the exponential map that arises in traditional Lie group variational methods. In addition, a reduced Lie–Poisson integrator is constructed and the resulting algorithms can naturally be implemented by fixed-point iteration. The proposed methods are validated by numerical simulations on \begin{document}$$ SO(3) $$\end{document}$ which demonstrate that they are comparable to variational Runge–Kutta–Munthe-Kaas methods in terms of computational efficiency. However, the methods we have proposed preserve the Lie group structure much more accurately and and exhibit better near energy preservation.  more » « less
Award ID(s):
1813635
PAR ID:
10354975
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Computational Dynamics
Volume:
0
Issue:
0
ISSN:
2158-2491
Page Range / eLocation ID:
0
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We establish an instantaneous smoothing property for decaying solutions on the half-line \begin{document}$$ (0, +\infty) $$\end{document} of certain degenerate Hilbert space-valued evolution equations arising in kinetic theory, including in particular the steady Boltzmann equation. Our results answer the two main open problems posed by Pogan and Zumbrun in their treatment of \begin{document}$ H^1 $$\end{document} stable manifolds of such equations, showing that \begin{document}$$ L^2_{loc} $$\end{document} solutions that remain sufficiently small in \begin{document}$$ L^\infty $$\end{document} (i) decay exponentially, and (ii) are \begin{document}$$ C^\infty $$\end{document} for \begin{document}$$ t>0 $$\end{document}, hence lie eventually in the \begin{document}$$ H^1 $$\end{document}$ stable manifold constructed by Pogan and Zumbrun. 
    more » « less
  2. Let \begin{document}$$ f_c(z) = z^2+c $$\end{document} for \begin{document}$$ c \in {\mathbb C} $$\end{document}. We show there exists a uniform upper bound on the number of points in \begin{document}$$ {\mathbb P}^1( {\mathbb C}) $$\end{document} that can be preperiodic for both \begin{document}$$ f_{c_1} $$\end{document} and \begin{document}$$ f_{c_2} $$\end{document}, for any pair \begin{document}$$ c_1\not = c_2 $$\end{document} in \begin{document}$$ {\mathbb C} $$\end{document}. The proof combines arithmetic ingredients with complex-analytic: we estimate an adelic energy pairing when the parameters lie in \begin{document}$$ \overline{\mathbb{Q}} $$\end{document}, building on the quantitative arithmetic equidistribution theorem of Favre and Rivera-Letelier, and we use distortion theorems in complex analysis to control the size of the intersection of distinct Julia sets. The proofs are effective, and we provide explicit constants for each of the results. 
    more » « less
  3. We propose a systematic procedure called the Clebsch canonization for obtaining a canonical Hamiltonian system that is related to a given Lie–Poisson equation via a momentum map. We describe both coordinate and geometric versions of the procedure, the latter apparently for the first time. We also find another momentum map so that the pair of momentum maps constitute a dual pair under a certain condition. The dual pair gives a concrete realization of what is commonly referred to as collectivization of Lie–Poisson systems. It also implies that solving the canonized system by symplectic Runge–Kutta methods yields so-called collective Lie–Poisson integrators that preserve the coadjoint orbits and hence the Casimirs exactly. We give a couple of examples, including the Kida vortex and the heavy top on a movable base with controls, which are Lie–Poisson systems on \begin{document}$$ \mathfrak{so}(2,1)^{*} $$\end{document} and \begin{document}$$ (\mathfrak{se}(3) \ltimes \mathbb{R}^{3})^{*} $$\end{document}, respectively. 
    more » « less
  4. For any finite horizon Sinai billiard map \begin{document}$ T $$\end{document} on the two-torus, we find \begin{document}$$ t_*>1 $$\end{document} such that for each \begin{document}$$ t\in (0,t_*) $$\end{document} there exists a unique equilibrium state \begin{document}$$ \mu_t $$\end{document} for \begin{document}$$ - t\log J^uT $$\end{document}, and \begin{document}$$ \mu_t $$\end{document} is \begin{document}$$ T $$\end{document}-adapted. (In particular, the SRB measure is the unique equilibrium state for \begin{document}$$ - \log J^uT $$\end{document}.) We show that \begin{document}$$ \mu_t $$\end{document} is exponentially mixing for Hölder observables, and the pressure function \begin{document}$$ P(t) = \sup_\mu \{h_\mu -\int t\log J^uT d \mu\} $$\end{document} is analytic on \begin{document}$$ (0,t_*) $$\end{document}. In addition, \begin{document}$$ P(t) $$\end{document} is strictly convex if and only if \begin{document}$$ \log J^uT $$\end{document} is not \begin{document}$$ \mu_t $$\end{document}-a.e. cohomologous to a constant, while, if there exist \begin{document}$$ t_a\ne t_b $$\end{document} with \begin{document}$$ \mu_{t_a} = \mu_{t_b} $$\end{document}, then \begin{document}$$ P(t) $$\end{document} is affine on \begin{document}$$ (0,t_*) $$\end{document}. An additional sparse recurrence condition gives \begin{document}$$ \lim_{t\downarrow 0} P(t) = P(0) $$\end{document}$. 
    more » « less
  5. Any \begin{document}$ C^d $$\end{document} conservative map \begin{document}$$ f $$\end{document} of the \begin{document}$$ d $$\end{document}-dimensional unit ball \begin{document}$$ {\mathbb B}^d $$\end{document}, \begin{document}$$ d\geq 2 $$\end{document}, can be realized by renormalized iteration of a \begin{document}$$ C^d $$\end{document} perturbation of identity: there exists a conservative diffeomorphism of \begin{document}$$ {\mathbb B}^d $$\end{document}, arbitrarily close to identity in the \begin{document}$$ C^d $$\end{document} topology, that has a periodic disc on which the return dynamics after a \begin{document}$$ C^d $$\end{document} change of coordinates is exactly \begin{document}$$ f $$\end{document}$. 
    more » « less