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.

A computationally efficient autoregressive method for generating phase screens with frozen flow and turbulence in optical simulations

Fractal phase screen generation algorithm for atmospheric turbulence

Production of Phase Screens for Simulation of Atmospheric Turbulence

Modeling astronomical adaptive optics performance with temporally filtered Wiener reconstruction of slope data

Data driven synthetic wavefront generation for boundary layer data

Generation of a large scale atmospheric turbulence based on an improved fractal algorithm

Modelling astronomical adaptive optics performance with temporally-filtered Wiener reconstruction of slope data

Gaussian phase autocorrelation as an accurate compensator for FFT-based atmospheric phase screen simulations

Atmospheric Phase Screen Simulation Using Zernike Polynomial

Image based automatic wind velocity profiler for adaptive optics

Advection of the image point in probabilistically-reconstructed phase spaces

Prediction of Optical Wave Phase Fluctuations in a Turbulent Atmosphere Based On Current Database

Real-Time Screen Space Rendering Method for Particle-Based Multiphase Fluid Simulation

Accurate Zernike-Corrected Phase Screens for Arbitrary Power Spectra

Atmospheric Turbulence with Kolmogorov Spectra: Software Simulation, Real-Time Reconstruction and Compensation by Means of Adaptive Optical System with Bimorph and Stacked-Actuator Deformable Mirrors

Zonal Representation-Based Residual-Wavefront Reconstruction Using Multiframe Shack-Hartmann Measurements

A Framework for Iterative Phase Retrieval Technique Integration into Atmospheric Adaptive Optics—Part I: Wavefront Sensing in Strong Scintillations

Single-Layer Phase Screen With Pointing Errors for Free Space Optical Communication

Kinetic-Based Multiphase Flow Simulation

A generalized approach to compensate for low and high frequency errors in FFT based phase screen simulations

Higher Order Phase Averaging for Highly Oscillatory Systems