[论文解读] Eliminating spurious eigenvalues in the analysis of incompressible fluids and other systems of differential-algebraic equations

M. Lisa Manning, Bassam Bamieh+1 more|arXiv (Cornell University)|May 10, 2007

Numerical methods for differential equations参考文献 30被引用 2

一句话总结

本文提出一种描述符记法框架，用于消除不可压缩流体及其他微分代数方程（DAE）稳定性分析中的虚假特征值。通过在不进行解析约化的前提下同时求解微分方程和代数约束，该方法避免了传统数值方法中常见的非物理解稳定性，确保了流体动力学稳定性问题中特征值计算的准确性。

ABSTRACT

We describe a general framework for avoiding spurious eigenvalues -- unphysical unstable eigenvalues that often occur in hydrodynamic stability problems. In two example problems, we show that when system stability is analyzed numerically using {\em descriptor} notation, spurious eigenvalues are eliminated. Descriptor notation is a generalized eigenvalue formulation for differential-algebraic equations that explicitly retains algebraic constraints. We propose that spurious eigenvalues are likely to occur when algebraic constraints are used to analytically reduce the number of independent variables in a differential-algebraic system of equations before the system is approximated numerically. In contrast, the simple and easily generalizable descriptor framework simultaneously solves the differential equations and algebraic constraints and is well-suited to stability analysis in these systems.

研究动机与目标

• 解决不可压缩流体流动数值稳定性分析中虚假特征值这一长期存在的问题。
• 识别虚假特征值的根源，即在数值近似前对 DAE 中的变量进行解析约化所导致。
• 提出一种描述符记法框架，在数值求解过程中保留代数约束，从而提高特征值计算的准确性。
• 展示描述符方法在消除流体动力学稳定性问题中非物理不稳定性方面的有效性。
• 提供一种通用且易于扩展的方法，适用于广泛类别的微分代数系统。

提出的方法

• 本文采用描述符记法，一种广义特征值公式化方法，可显式保留 DAE 中的代数约束。
• 避免对独立变量数量进行解析约化，从而在数值近似过程中保持系统完整结构。
• 该方法在统一的矩阵框架内同时求解微分方程和代数约束。
• 通过从描述符公式化中导出的广义特征值问题来分析稳定性。
• 将该方法应用于两个流体动力学稳定性问题，以验证其有效性。
• 该框架设计为通用且可轻松扩展至流体动力学以外的其他 DAE 系统。

实验结果

研究问题

• RQ1 为何在不可压缩流体系统的数值稳定性分析中会出现虚假特征值？
• RQ2 在 DAE 中对变量进行解析约化如何导致非物理特征值的出现？
• RQ3 基于描述符的公式化能否在保持稳定性分析精度的同时消除虚假特征值？
• RQ4 该描述符框架是否可推广至流体动力学以外的其他微分代数系统？
• RQ5 在数值近似过程中保留代数约束对特征值计算有何影响？

主要发现

• 当在数值近似前解析地消除代数约束时，虚假特征值会随之出现，从而引入非物理不稳定性。
• 通过在数值求解过程中保持系统完整结构，描述符框架成功消除了虚假特征值。
• 该方法通过在不进行变量约化的情况下同时求解微分方程和代数方程，维持了稳定性分析的准确性。
• 该方法在两个涉及不可压缩流体流动的测试案例中表现有效，展现出鲁棒性和可靠性。
• 描述符公式化为 DAE 系统的稳定性分析提供了一种通用、稳定且易于推广的框架。
• 由于能够显式且一致地处理约束，该框架非常适用于数值稳定性分析。