Linear model reduction using SPOD modes

Peter Frame,Cong Lin,Oliver Schmidt,Aaron Towne

2024-05-10

Abstract:The majority of model reduction approaches use an efficient representation of the state and then derive equations to temporally evolve the coefficients that encode the state in the representation. In this paper, we instead employ an efficient representation of the entire trajectory of the state over some time interval and solve for the coefficients that define the trajectory on the interval. We use spectral proper orthogonal decomposition (SPOD) modes, in particular, which possess properties that make them suitable for model reduction and are known to provide an accurate representation of trajectories. In fact, with the same number of total coefficients, the SPOD representation is substantially more accurate than any representation formed by specifying the coefficients in a spatial (e.g., POD) basis for the many time steps that make up the interval. We develop a method to solve for the SPOD coefficients that encode the trajectories in forced linear dynamical systems given the forcing and initial condition, thereby obtaining the accurate representation of the trajectory. We apply the method to two examples, a linearized Ginzburg-Landau problem and an advection-diffusion problem. In both, the error of the proposed method is orders of magnitude lower than both POD-Galerkin and balanced truncation applied to the same problem, as well as the most accurate solution within the span of the POD modes. The method is also fast, with CPU time comparable to or lower than both benchmarks in the examples we present.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to represent state trajectories efficiently and accurately in linear dynamical systems. Traditional methods usually first find an effective compressed representation of the state of a system at specific time points, and then derive equations for the evolution of these representation coefficients over time. However, this method is less efficient for representing trajectories over the entire time interval because it is necessary to specify mode coefficients for each time step, and these coefficients are highly correlated between adjacent time steps.

To overcome this limitation, this paper proposes a new method, that is, using Spectral Proper Orthogonal Decomposition (SPOD) modes to represent trajectories over the entire time interval. SPOD modes have the following advantages:

1. **High efficiency**: Compared with the traditional spatial POD mode, using the same total number of coefficients, the SPOD representation can represent trajectories more accurately.
2. **Low error**: In two example problems (the linearized Ginzburg - Landau problem and the convection - diffusion problem), the error of this method is several orders of magnitude lower than that of the POD - Galerkin and balanced truncation methods.
3. **Fast calculation**: The CPU time of this method is comparable to or lower than that of the benchmark method.

Specifically, the paper develops an algorithm to quickly solve for SPOD coefficients given initial conditions and external forces, thereby obtaining an accurate representation of the trajectory. Through this method, the accuracy of the model can be significantly improved while maintaining computational efficiency.