[논문 리뷰] A Multilevel Correction Scheme for Nonsymmetric Eigenvalue Problems by Finite Element Methods

Hehu Xie, Zhimin Zhang|arXiv (Cornell University)|2015. 05. 23.

Advanced Numerical Methods in Computational Mathematics참고 문헌 22인용 수 8

한 줄 요약

이 논문은 비대칭 고유값 문제에 대해 유한요소법을 사용한 다수준 보정 체계를 제안하며, 고유값 상승도 1이 되어야 한다는 조건 없이 효율적이고 정확한 근사가 가능하도록 한다. 가장 굵은 메esh에서만 작은 고유값 문제를 풀고, 더 세밀한 공간에서의 소스 문제를 풀음으로써 최적 수렴 속도를 달성하고 계산 효율을 크게 향상시킨다. 특히 다중격자 또는 적응형 메쉬 개선과 결합할 경우 더욱 효과적이다.

ABSTRACT

A multilevel correction scheme is proposed to solve defective and nodefective of nonsymmetric partial differential operators by the finite element method. The method includes multi correction steps in a sequence of finite element spaces. In each correction step, we only need to solve two source problems on a finer finite element space and two eigenvalue problems on the coarsest finite element space. The accuracy of the eigenpair approximation is improved after each correction step. This correction scheme improves overall efficiency of the finite element method in solving nonsymmetric eigenvalue problems.

연구 동기 및 목표

• 유체역학 및 안정성 분석에서 발생하는 비대칭 고유값 문제를 효율적으로 해결하기 위한 유한요소법을 개발한다.
• 이전에 대칭 문제에 사용된 다수준 보정 기법을 고유값 상승도가 1이 되어야 한다는 조건 없이 비대칭 연산자로 확장한다.
• 세밀한 유한요소 공간에서의 대규모 고유값 문제를 풀이하는 데 드는 비용을 줄임으로써 계산 효율을 향상시킨다.
• 다중격자 및 적응형 유한요소 방법과의 원활한 통합을 가능하게 하여 최적 수렴 속도와 병렬 처리를 달성한다.

제안 방법

• 점점 더 세밀해지는 포함 관계를 가진 유한요소 공간의 수열을 사용한다.
• 각 보정 단계에서 이전 고유함수를 하중 벡터로 사용하여 확장된 유한요소 공간에서 소스 문제를 푼다.
• 보정된 고유함수를 포함한 공간에서 새로운 고유값 문제를 풀며, 고유값 계산에는 오직 가장 굵은 공간만 사용한다.
• 반복적으로 보정 단계를 거치며 고유값과 고유함수 성분을 동시에 개선함으로써 고유쌍 근사치를 향상시킨다.
• 세밀한 격자에서의 전체 고유값 문제를 대체로, 더 세밀한 공간에서의 저비용 소스 문제와 가장 굵은 공간에서의 작은 고유값 문제로 대체함으로써 전체 고유값 문제를 풀지 않는다.
• 경계값 문제에 대한 다중격자 솔버와 ZZ 복구 기반 후행 오차 추정기를 사용한 적응형 메쉬 개선과 호환된다.

실험 결과

연구 질문

• RQ1 고유값 상승도가 1이 되어야 한다는 가정 없이 비대칭 고유값 문제에 대해 다수준 보정 체계를 설계할 수 있는가?
• RQ2 유한요소법을 사용하여 대규모 비대칭 고유값 문제를 풀이하는 데 드는 계산 비용을 어떻게 줄일 수 있는가?
• RQ3 특이 고유함수를 가진 문제(예: L자형 도메인)에 적용했을 때 다수준 보정 체계가 최적 수렴 속도를 유지할 수 있는가?
• RQ4 미분 연산자에 복소수 벡터장이 포함되었을 경우 이 방법의 성능은 어떠한가?
• RQ5 다중격자 및 적응형 유한요소 기법과 효과적으로 통합하여 효율을 향상시킬 수 있는가?

주요 결과

• 특이 고유함수를 가진 문제(예: 비볼록 도메인인 L자형 도메인)에 대해서도 다수준 보정 체계는 고유값 및 고유함수 근사에 대해 최적 수렴 속도를 달성한다.
• 수치적 결과는 일정한 메쉬에서 첫 번째 고유값에 대해 수렴 속도가 $h^{4/3}$의 순서로 유지됨을 확인하였으며, 이는 이론적 기대와 일치한다.
• 적응형 메쉬에서는 오차가 기대되는 비율로 감소함을 확인하였고, 초기 메쉬의 영향은 최종 정확도에 거의 영향을 주지 않았다.
• 복소수 벡터장이 포함된 미분 연산자에 대해서도 방법이 잘 작동하며, 첫 여섯 고유값에 대해 수렴 속도가 매우 빠르게 나타났다.
• 다중공간 접근법은 첫 여섯 고유값의 합을 계산할 때 다중격자 방법보다 더 높은 정확도를 제공함으로써 고정밀도 환경에서의 우수성을 입증하였다.
• 이 체계는 매우 효율적이다: 가장 굵은 공간에서의 작은 고유값 문제와 더 세밀한 공간에서의 소스 문제만 풀어도 가장 세밀한 공간에서 전체 고유값 문제를 풀었을 때와 동일한 정확도를 달성한다.