The incompressible two-dimensional Euler equations on a sphere constitute a fundamental model in hydrodynamics. The long-time behaviour of solutions is largely unknown; statistical mechanics predicts a steady vorticity configuration, but detailed numerical results in the literature contradict this theory, yielding instead persistent unsteadiness. Such numerical results were obtained using artificial hyperviscosity to account for the cascade of enstrophy into smaller scales. Hyperviscosity, however, destroys the underlying geometry of the phase flow (such as conservation of Casimir functions), and therefore might affect the qualitative long-time behaviour. Here, we develop an efficient numerical method for long-time simulations that preserve the geometric features of the exact flow, in particular conservation of Casimirs. Long-time simulations on a non-rotating sphere then reveal three possible outcomes for generic initial conditions: the formation of either 2, 3 or 4 coherent vortex structures. These numerical results contradict the statistical mechanics theory and show that previous numerical results, suggesting 4 coherent vortex structures as the generic behaviour, display only a special case. Through integrability theory for point vortex dynamics on the sphere we present a theoretical model which describes the mechanism by which the three observed regimes appear. We show that there is a correlation between a first integral γ (the ratio of total angular momentum and the square root of enstrophy) and the long-time behaviour: γ small, intermediate and large yields most likely 4, 3 or 2 coherent vortex formations. Our findings thus suggest that the likely long-time behaviour can be predicted from the first integral γ.

A Casimir preserving scheme for long-time simulation of spherical ideal hydrodynamics

Modin, Klas
;
Viviani, Milo
2020

Abstract

The incompressible two-dimensional Euler equations on a sphere constitute a fundamental model in hydrodynamics. The long-time behaviour of solutions is largely unknown; statistical mechanics predicts a steady vorticity configuration, but detailed numerical results in the literature contradict this theory, yielding instead persistent unsteadiness. Such numerical results were obtained using artificial hyperviscosity to account for the cascade of enstrophy into smaller scales. Hyperviscosity, however, destroys the underlying geometry of the phase flow (such as conservation of Casimir functions), and therefore might affect the qualitative long-time behaviour. Here, we develop an efficient numerical method for long-time simulations that preserve the geometric features of the exact flow, in particular conservation of Casimirs. Long-time simulations on a non-rotating sphere then reveal three possible outcomes for generic initial conditions: the formation of either 2, 3 or 4 coherent vortex structures. These numerical results contradict the statistical mechanics theory and show that previous numerical results, suggesting 4 coherent vortex structures as the generic behaviour, display only a special case. Through integrability theory for point vortex dynamics on the sphere we present a theoretical model which describes the mechanism by which the three observed regimes appear. We show that there is a correlation between a first integral γ (the ratio of total angular momentum and the square root of enstrophy) and the long-time behaviour: γ small, intermediate and large yields most likely 4, 3 or 2 coherent vortex formations. Our findings thus suggest that the likely long-time behaviour can be predicted from the first integral γ.
2020
Settore MATH-04/A - Fisica matematica
Settore MATH-05/A - Analisi numerica
Rotating flows; Computational methods; Hamiltonian theory; 2D turbulence; Euler equations; Matrix hydrodynamics; Conservative numerical simulations
File in questo prodotto:
File Dimensione Formato  
2. A Casimir preserving scheme for long-time.pdf

accesso aperto

Tipologia: Published version
Licenza: Creative Commons
Dimensione 2.81 MB
Formato Adobe PDF
2.81 MB Adobe PDF

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11384/148909
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 23
  • ???jsp.display-item.citation.isi??? 23
  • OpenAlex 31
social impact