TY - JOUR

T1 - Stability of leapfrogging vortex pairs

T2 - A semi-analytic approach

AU - Behring, Brandon M.

AU - Goodman, Roy H.

N1 - Publisher Copyright:
© 2019 American Physical Society.
Copyright:
Copyright 2020 Elsevier B.V., All rights reserved.

PY - 2019/12/26

Y1 - 2019/12/26

N2 - We investigate the stability of a one-parameter family of periodic solutions of the four-vortex problem known as "leapfrogging" orbits. These solutions, which consist of two pairs of identical yet oppositely signed vortices, were known to W. Gröbli [Ph.D. thesis, Georg-August-Universität Göttingen, 1877] and A. E. H. Love [Proc. London Math. Soc. 1, 185 (1883)10.1112/plms/s1-25.1.185] and can be parameterized by a dimensionless parameter α related to the geometry of the initial configuration. Simulations by Acheson [Eur. J. Phys. 21, 269 (2000)10.1088/0143-0807/21/3/310] and numerical Floquet analysis by Tophøj and Aref [Phys. Fluids 25, 014107 (2013)10.1063/1.4774333] both indicate, to many digits, that the bifurcation occurs when 1/α=φ2, where φ is the golden ratio. This study aims to explain the origin of this remarkable value. Using a trick from the gravitational two-body problem, we change variables to render the Floquet problem in an explicit form that is more amenable to analysis. We then implement G. W. Hill's [Acta Math. 8, 1 (1886)10.1007/BF02417081] method of harmonic balance to high order using computer algebra to construct a rapidly converging sequence of asymptotic approximations to the bifurcation value, confirming the value found earlier.

AB - We investigate the stability of a one-parameter family of periodic solutions of the four-vortex problem known as "leapfrogging" orbits. These solutions, which consist of two pairs of identical yet oppositely signed vortices, were known to W. Gröbli [Ph.D. thesis, Georg-August-Universität Göttingen, 1877] and A. E. H. Love [Proc. London Math. Soc. 1, 185 (1883)10.1112/plms/s1-25.1.185] and can be parameterized by a dimensionless parameter α related to the geometry of the initial configuration. Simulations by Acheson [Eur. J. Phys. 21, 269 (2000)10.1088/0143-0807/21/3/310] and numerical Floquet analysis by Tophøj and Aref [Phys. Fluids 25, 014107 (2013)10.1063/1.4774333] both indicate, to many digits, that the bifurcation occurs when 1/α=φ2, where φ is the golden ratio. This study aims to explain the origin of this remarkable value. Using a trick from the gravitational two-body problem, we change variables to render the Floquet problem in an explicit form that is more amenable to analysis. We then implement G. W. Hill's [Acta Math. 8, 1 (1886)10.1007/BF02417081] method of harmonic balance to high order using computer algebra to construct a rapidly converging sequence of asymptotic approximations to the bifurcation value, confirming the value found earlier.

UR - http://www.scopus.com/inward/record.url?scp=85077808501&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85077808501&partnerID=8YFLogxK

U2 - 10.1103/PhysRevFluids.4.124703

DO - 10.1103/PhysRevFluids.4.124703

M3 - Article

AN - SCOPUS:85077808501

VL - 4

JO - Physical Review Fluids

JF - Physical Review Fluids

SN - 2469-990X

IS - 12

M1 - 124703

ER -