“…Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a coupled set of steady equations in integer harmonics that can be tackled with the acceleration techniques a↵orded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6], and Time-Spectral [7,8,9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques [10,11] can be applied to provide the ability to perform design optimization without resorting to costly unsteady adjoint methods.…”