[논문 리뷰] High-order generalized-$\alpha$ methods

Quanling Deng, Pouria Behnoudfar외 1명|arXiv (Cornell University)|2019. 02. 14.

Numerical methods for differential equations참고 문헌 8인용 수 4

한 줄 요약

이 논문은 제3차 정확도를 갖는 일반화된-$\alpha$ 시간 적분법을 제안하며, 무조건적 안정성과 사용자 제어가 가능한 고주파수 감쇠를 유지한다. 이는 이중차 정확도 일반화된-$\alpha$ 방법을 이산 형식에서 고차수 테일러 전개를 통해 확장함으로써 달성된다. 이 방법은 원래 방법의 무조건적 안정성과 단일 스텝 구조를 유지하면서 제3차 정확도를 확보한다. 수치적 감쇠는 파rameter $\alpha_m$ 및 $\alpha_f$를 통해 제어된다.

ABSTRACT

The generalized-$\\alpha$ method encompasses a wide range of time integrators. The method possesses high-frequency dissipation while minimizing unwanted low-frequency dissipation and the numerical dissipation can be controlled by the user. The method is unconditionally stable and is of second-order accuracy in time. We extend the second-order generalized-$\\alpha$ method to third-order in time while the numerical dissipation can be controlled in a similar fashion. We establish that the third-order method is unconditionally stable. We discuss a possible path to the generalization to higher order schemes. All these high-order schemes can be easily implemented into programs that already contain the second-order generalized-$\\alpha$ method.

연구 동기 및 목표

• 이중차 정확도 일반화된-$\alpha$ 방법을 시간에 대해 제3차 정확도로 확장하면서도 그 핵심 특성을 유지하고자 한다.
• 새로운 고차수 스킴에서 무조건적 안정성을 유지하고자 한다.
• 일반화된-$\alpha$ 방법의 상징적인 특징인 고주파수 수치적 감쇠에 대한 사용자 제어를 유지하고자 한다.
• 기존의 이중차 정확도 일반화된-$\alpha$ 코드베이스에 원활하게 통합될 수 있도록 하고자 한다.
• 더 높은 차수의 스킴로 일반화하기 위한 프레임워크를 제공하고자 한다.

제안 방법

• 이산 형식에서 $U_{n}^{\alpha_f}$ 및 $V_{n}^{\alpha_m}$를 선형 보조단계가 아니라 고차수 테일러 전개로 간주함으로써 일반화된-$\alpha$ 방법을 재구성한다.
• 제3차 정확도를 달성하기 위해 $U_{n+\alpha_f}$ 및 $V_{n+\alpha_m}$의 테일러 전개에서 고차수 항을 나타내는 추가적인 파rameter $\alpha_f$ 및 $\alpha_m$을 도입한다.
• 자유 파rameter $\alpha_m$ 및 $\alpha_f$를 결정하기 위해 테일러 급수 분석을 사용하여 시간 스텝 스킴을 유도하고 최적의 정확도를 확보한다.
• 증폭 행렬을 구성하고 고유값을 분석하여 무조건적 안정성 영역을 확립한다.
• 고주파수 고유값 $\rho_\infty$를 $[0,1]$ 범위의 사용자 정의 값으로 설정함으로써 수치적 감쇠를 제어한다.
• 기호적 및 수치적 분석을 통해 무조건적 안정성을 검증하고 안정성 및 원하는 감쇠를 위한 파rameter 영역을 규명한다.

실험 결과

연구 질문

• RQ1 일반화된-$\alpha$ 방법은 무조건적 안정성을 잃지 않고 시간에 대해 제3차 정확도로 확장될 수 있는가?
• RQ2 고차수 시간 적분법에서 사용자 제어가 가능한 고주파수 감쇠는 어떻게 유지할 수 있는가?
• RQ3 어떤 파rameter 선택 $\alpha_m$ 및 $\alpha_f$가 제3차 정확도와 무조건적 안정성을 동시에 확보하는가?
• RQ4 결과로 도출된 스킴은 기존의 이중차 정확도 일반화된-$\alpha$ 코드에 최소한의 변경으로 구현될 수 있는가?
• RQ5 제3차 일반화된-$\alpha$ 방법의 안정성 영역는 무엇이며, 이는 이중차 버전과 비교해 어떻게 다른가?

주요 결과

• 제안된 제3차 일반화된-$\alpha$ 방법은 해와 속도 항의 고차수 테일러 전개를 통해 시간에 대해 공식적으로 제3차 정확도를 달성한다.
• 모든 $\alpha_m \geq \frac{7}{12}$ 및 $\frac{1}{2} \leq \alpha_f \leq \alpha_m - \frac{1}{12}$ 에 대해 무조건적 안정성이 보장되며, 안정성 영역는 기호적 및 수치적으로 확인되었다.
• 고주파수에서의 수치적 감쇠는 파rameter $\rho_\infty$에 의해 제어되며, $0 \leq \rho_\infty \leq 1$ 에 대해 $\alpha_m = \frac{13 + 20\rho_\infty - 5\rho_\infty^2}{12(\rho_\infty + 1)^2}$ 및 $\alpha_f = \frac{1 + 3\rho_\infty}{2(\rho_\infty + 1)^2}$ 로 설정된다.
• 증폭 행렬 고유값은 무조건적 안정성 영역 전역에서 크기가 1 이내로 제한되며, 이는 안정성을 확인한다.
• 수치적 플롯을 통해 제3차 스킴이 이중차 스킴보다 더 큰 무조건적 안정성 영역를 갖는다.
• 이 방법은 여전히 단일 스텝, 암시적 스킴이며, 기존의 이중차 정확도 일반화된-$\alpha$ 솔버에 최소한의 코드 변경으로 쉽게 통합될 수 있다.