EGU21-7543, updated on 04 Mar 2021
EGU General Assembly 2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.

Higher order phase averaging for big timesteps

Werner Bauer and Colin Cotter
Werner Bauer and Colin Cotter
  • Imperial College London, Department of Mathematics, London, United Kingdom of Great Britain – England, Scotland, Wales

We introduce a higher order phase averaging method for nonlinear PDEs. Our method is suitable for highly oscillatory systems of nonlinear PDEs that generate slow motion through resonance between fast frequencies, such as is the case for rotating fluids with small but finite Rossby number. Phase averaging is a technique to filter fast motions from the dynamics whilst still accounting for their effect on the slow dynamics. In the small Rossby number limit of the phase averaged rotating shallow water equations, one recovers the quasi-geostrophic equations (as shown by Schochet, Majda and others). Peddle et al. 2017, Haut and Wingate 2014, have shown that phase averaging at finite Rossby number allows to take larger timesteps than would otherwise be possible. This was used as a coarse propagator (large timesteps at lower accuracy) for a Parareal method where corrections were made using a standard timestepping method with small timesteps.

In this contribution, we introduce an additional phase variable in the exponential time integrator that allows us to derive arbitrary order averaging methods that can be used as more accurate corrections to the basic phase averaged model, without needing small timesteps. We envisage their use as part of a time-parallel algorithm based on deferred corrections to the basic average. We illustrate the properties of this method on an ODE that describes the dynamics of a swinging spring, a model due to Peter Lynch. Although idealized, this model shows an interesting analogy to geophysical flows as it exhibits a high sensitivity of small scale oscillation on the large scale dynamics. On this example, we show convergence to the non-averaged (exact) solution with increasing approximation order also for finite averaging windows. At zeroth order, our method coincides with that in Peddle et al. 2017, Haut and Wingate 2014, but at higher order it is more accurate in the sense that it better approximates the faster oscillations around the slow manifold.

How to cite: Bauer, W. and Cotter, C.: Higher order phase averaging for big timesteps, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-7543,, 2021.

Display materials

Display file

Comments on the display material

to access the discussion