469 research outputs found

    Weighted Birkhoff Averages and the Parameterization Method

    Full text link
    This work provides a systematic recipe for computing accurate high order Fourier expansions of quasiperiodic invariant circles in area preserving maps. The recipe requires only a finite data set sampled from the quasiperiodic circle. Our approach, being based on the parameterization method, uses a Newton scheme to iteratively solve a conjugacy equation describing the invariant circle. A critical step in properly formulating the conjugacy equation is to determine the rotation number of the quasiperiodic subsystem. For this we exploit a the weighted Birkhoff averaging method. This approach facilities accurate computation of the rotation number given nothing but the already mentioned orbit data. The weighted Birkhoff averages also facilitate the computation of other integral observables like Fourier coefficients of the parameterization of the invariant circle. Since the parameterization method is based on a Newton scheme, we only need to approximate a small number of Fourier coefficients with low accuracy to find a good enough initial approximation so that Newton converges. Moreover, the Fourier coefficients may be computed independently, so we can sample the higher modes to guess the decay rate of the Fourier coefficients. This allows us to choose, a-priori, an appropriate number of modes in the truncation. We illustrate the utility of the approach for explicit example systems including the area preserving Henon map and the standard map. We present example computations for invariant circles with period as low as 1 and up to more than 100. We also employ a numerical continuation scheme to compute large numbers of quasiperiodic circles in these systems. During the continuation we monitor the Sobolev norm of the Parameterization to automatically detect the breakdown of the family.Comment: 38 pages, 15 figure

    Critical homoclinics in a restricted four body problem: numerical continuation and center manifold computations

    Full text link
    The present work studies the robustness of certain basic homoclinic motions in an equilateral restricted four body problem. The problem can be viewed as a two parameter family of conservative autonomous vector fields. The main tools are numerical continuation techniques for homoclinic and periodic orbits, as well as formal series methods for computing normal forms and center stable/unstable manifold parameterizations. After careful numerical study of a number of special cases we formulate several conjectures about the global bifurcations of the homoclinic families.Comment: 38 pages, 21 figures, fixed several typos, expanded the introduction and added a new appendix about numerical continuation of orbit
    corecore