2024/07/25 by Artur Gesla, Gesla, Artur, Yohann Duguet +5
Physics and Astronomy · Engineering · Mathematics · #Quantum chaos and dynamical systems #Spacecraft Dynamics and Control #Advanced Differential Equations and Dynamical Systems
paper · pdf · doi:10.48550/arxiv.2407.18230
We analyse and compare several algorithms to compute numerically periodic solutions of high-dimensional dynamical systems and investigate their Floquet stability without building the monodromy matrix. The solution and its perturbation are discretised in time either using finite differences, Fourier-Galerkin or Chebyshev expansions. The resulting nonlinear set of equations describing the periodic orbit is solved using a Newton-Raphson algorithm. The linearised equations determining the stability lead to a generalised eigenvalue problem. Unlike the Fourier-Galerkin method, the use of Chebyshev polynomials or finite differences has the advantage that the relevant Floquet exponents are directly given without the well known issue of having to sort out the eigenvalues. The speed of convergence of these three methods is illustrated with examples from the Lorenz system, the Langford system and a two-dimensional thermal convection flow inside a differentially heated cavity. This last example demonstrates the potential of the newly proposed Chebyshev expansion for large-scale problems arising from the discretisation of the incompressible Navier-Stokes equations.