Higher Order Phase Averaging for Highly Oscillatory Systems
Werner Bauer,Colin Cotter,Beth Wingate
DOI: https://doi.org/10.1137/21M1430546
2022-09-02
Multiscale Modeling and Simulation
Abstract:Multiscale Modeling &Simulation, Volume 20, Issue 3, Page 936-956, September 2022. We introduce a higher order phase averaging method for nonlinear oscillatory systems. Phase averaging is a technique to filter fast motions from the dynamics while still accounting for their effect on the slow dynamics. Phase averaging is useful for deriving reduced models that can be solved numerically with more efficiency, since larger timesteps can be taken. Recently, Haut and Wingate [SIAM J. Sci. Comput., 36 (2014) pp. A693–A713] introduced the idea of computing finite window numerical phase averages in parallel as the basis for a coarse propagator for a parallel-in-time algorithm. In this contribution, we provide a framework for higher order phase averages that aims to better approximate the unaveraged system while still filtering fast motions. While the basic phase average assumes that the solution is independent of changes of phase, the higher order method expands the phase dependency in a basis which the equations are projected onto. In this new framework, the original numerical phase averaging formulation arises as the lowest order version of this expansion in which the nonlinearity is projected onto the space of functions that are independent of the phase. Our new projection onto functions that are [math]th degree polynomials in the phase gives rise to higher order corrections to the phase averaging formulation. We illustrate the properties of this method on an ODE that describes the dynamics of a swinging spring due to Lynch (2002). Although idealized, this model shows an interesting analogy to geophysical flows as it exhibits a slow dynamics that arises through the resonance between fast oscillations. On this example, we show convergence to the nonaveraged (exact) solution with increasing approximation order also for finite averaging windows. At zeroth order, our method coincides with a standard phase average, but at higher order it is more accurate in the sense that solutions of the phase averaged model track the solutions of the unaveraged equations more accurately.
mathematics, interdisciplinary applications,physics, mathematical