Scholarlink

<p>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.</p>

错误：搜索内容不能为空，请输入英文关键词
错误：关键词超出字数限制，请精简
高级检索

Computation and stability analysis of periodic orbits using finite differences, Fourier or Chebyshev spectral expansions in time

  • 期刊论文

  • 发表日期：

    2025年9月26日

查看全文
  • Artur Gesla,
  • Yohann Duguet,
  • Patrick Le Quéré,
  • Laurent Martin Witkowski

摘要

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.