Interpolatory model reduction of quadratic-bilinear dynamical systems with quadratic-bilinear outputs

This paper extends interpolatory model reduction to systems with (up to) quadratic-bilinear dynamics and quadratic-bilinear outputs. These systems are referred to as QB-QB systems and arise in a number of applications, including fluid dynamics, optimal control, and uncertainty quantification. In the interpolatory approach, the reduced order models (ROMs) are based on a Petrov-Galerkin projection, and the projection matrices are constructed so that transfer function components of the ROM interpolate the corresponding transfer function components of the original system. To extend the approach to systems with QB outputs, this paper derives system transfer functions and sufficient conditions on the projection matrices that guarantee the aforementioned interpolation properties. Alternatively, if the system has linear dynamics and quadratic outputs, one can introduce auxiliary state variables to transform it into a system with QB dynamics and linear outputs to which known interpolatory model reduction can be applied. This transformation approach is compared with the proposed extension that directly treats quadratic outputs. The comparison shows that transformation hides the problem structure. Numerical examples illustrate that keeping the original QB-QB structure leads to ROMs with better approximation properties.