[논문 리뷰] Eliminating spurious eigenvalues in the analysis of incompressible fluids and other systems of differential-algebraic equations

M. Lisa Manning, Bassam Bamieh외 1명|arXiv (Cornell University)|2007. 05. 10.

Numerical methods for differential equations참고 문헌 30인용 수 2

한 줄 요약

이 논문은 비압축성 유체 및 기타 미분-대수방정식(DAEs)의 안정성 분석에서 임의의 고유값을 제거하기 위한 기호 표기 프레임워크를 제안한다. 분석적 감소 없이 미분 및 대수적 제약 조건을 동시에 풀음으로써, 전통적인 수치적 접근 방식에서 흔히 발생하는 물리적으로 비현실적인 불안정성을 피하여 유체역학적 안정성 문제에서 정확한 고유값 계산을 보장한다.

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 시스템의 안정성 분석을 위한 일반적이고 안정적이며 쉽게 일반화할 수 있는 프레임워크를 제공한다.
• 이 프레임워크는 제약 조건을 명시적이고 일관적으로 다룰 수 있는 능력 덕분에 수치적 안정성 분석에 매우 적합하다.