Abstract:We discuss model reduction for a particular class of quadratic-bilinear (QB) descriptor systems. The main goal of this article is to extend the recently studied interpolation-based optimal model reduction framework for QBODEs [Benner et al. '16] to a class of descriptor systems in an efficient and reliable way. Recently, it has been shown in the case of linear or bilinear systems that a direct extension of interpolation-based model reduction techniques to descriptor systems, without any modifications, may lead to poor reduced-order systems. Therefore, for the analysis, we aim at transforming the considered QB descriptor system into an equivalent QBODE system by means of projectors for which standard model reduction techniques for QBODEs can be employed, including aforementioned interpolation scheme. Subsequently, we discuss related computational issues, thus resulting in a modified algorithm that allows us to construct \emph{near}--optimal reduced-order systems without explicitly computing the projectors used in the analysis. The efficiency of the proposed algorithm is illustrated by means of a numerical example, obtained via semi-discretization of the Navier-Stokes equations.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an effective model order reduction method for a specific type of Quadratic - Bilinear (QB) descriptor systems. Specifically, the author hopes to extend the recently studied interpolation - based optimal model order reduction framework from QBODEs (Quadratic - Bilinear Ordinary Differential Equations) to descriptor systems to ensure high efficiency and reliability when dealing with such complex systems.
### Specific Background and Challenges of the Problem
1. **Computational Cost of High - Fidelity Models**:
Numerical simulations of high - fidelity dynamic systems are usually very expensive and inefficient, which makes it impractical or even impossible to directly use these high - fidelity models for engineering design research.
2. **Importance of Model Order Reduction**:
Model Order Reduction (MOR) aims to construct surrogate models (i.e., reduced - order models) that, although of lower complexity, can reproduce the dynamic behavior of high - fidelity systems. For linear systems, MOR techniques are very mature and widely used in many fields.
3. **Challenges of Nonlinear Systems**:
For nonlinear systems (such as quadratic - bilinear systems), although some MOR techniques have been successfully applied, there are still many challenges, especially for descriptor systems (DAEs) containing algebraic constraints. Directly applying interpolation - based model order reduction techniques to descriptor systems without modification may lead to poorly performing reduced - order models.
4. **Application of Navier - Stokes Equations**:
The paper pays special attention to how to apply these techniques to the semi - discretized model of Navier - Stokes equations. The Navier - Stokes equations are used to describe fluid dynamics, and their numerical solution usually involves complex quadratic - bilinear descriptor systems.
### Overview of the Solution
To address the above challenges, the author proposes an iterative model order reduction scheme. The main steps are as follows:
1. **Conversion to an Equivalent QBODE System**:
Transform the original quadratic - bilinear descriptor system into an equivalent quadratic - bilinear ordinary differential equation system (QBODE) through a projection transformation, so that standard model order reduction techniques can be applied.
2. **Improved Algorithm**:
An improved algorithm is proposed, which can construct an approximately optimal reduced - order system without explicitly computing the projection matrix. This greatly improves computational efficiency and avoids the computational burden caused by explicitly computing the projection matrix.
3. **Numerical Verification**:
Numerical experiments are carried out through a lid - driven cavity model based on Navier - Stokes equations to verify the effectiveness of the proposed method.
### Summary of Mathematical Formulas
- System Form:
\[
E_{11} \dot{v}(t) = A_{11} v(t) + A_{12} p(t) + H v(t) \otimes v(t) + \sum_{k = 1}^m N_k v(t) u_k(t) + B_1 u(t),
\]
\[
0 = A_{21} v(t) + B_2 u(t), \quad v(0) = v_0,
\]
\[
y(t) = C_1 v(t) + C_2 p(t).
\]
- System After Projection Transformation:
\[
\Pi_l E_{11} \Pi_r \dot{\tilde{v}}(t) = \Pi_l A_{11} \Pi_r \tilde{v}(t) + \Pi_l H (\Pi_r \tilde{v}(t) \otimes \Pi_r \tilde{v}(t)) + \sum_{k = 1}^m \Pi_l N_k \Pi_r \tilde{v}(t) u_k(t) + \Pi_l B_1 u(t),
\]
\[
y(t) = C v(t) + C_H \Pi_r \tilde{v}(t) \o