[논문 리뷰] Local and Parallel Finite Element Algorithm Based On Multilevel Discretization for Eigenvalue Problem

Yu Li, Xiaole Han외 2명|arXiv (Cornell University)|2014. 01. 20.

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

한 줄 요약

이 논문은 다중 격자 해석과 국소 경계값 문제의 병렬 해법을 기반으로, 미세 체계 고유값 문제를 일련의 거친 체계 고유값 문제와 국소 경계값 문제로 감소시켜 고유값 문제를 해결하기 위한 다수준 국소 및 병렬 유한요소 알고리즘을 제안한다. 이 방법은 프로세서당 최적의 계산 작업을 유지하면서도 최적 수렴 속도를 달성하여 기존의 유한요소 접근 방식에 비해 효율성이 크게 향상된다.

ABSTRACT

A local and parallel algorithm based on the multilevel discretization is proposed in this paper to solve the eigenvalue problem by the finite element method. With this new scheme, solving the eigenvalue problem in the finest grid is transferred to solutions of the eigenvalue problems on the coarsest mesh and a series of solutions of boundary value problems by using the local and parallel algorithm. The computational work in each processor can reach the optimal order. Therefore, this type of multilevel local and parallel method improves the overall efficiency of solving the eigenvalue problem. Some numerical experiments are presented to validate the efficiency of the new method.

연구 동기 및 목표

• 과학 계산에서 발생하는 대규모 고유값 문제를 효율적으로 해결하기 위한 수치적 방법을 개발하기 위해.
• 미세 메esh에서의 유한요소 고유값 문제를 해결하는 데 드는 계산 비용을 줄이기 위해.
• 국소적이고 병렬적인 계산을 통해 프로세서당 계산 작업을 최소화하면서도 최적의 수렴 속도를 달성하기 위해.
• 다수준 보정 및 국소 유한요소 기법을 고유값 문제에 확장하여 확장성 향상을 도모하기 위해.

제안 방법

• 이 방법은 미세 체계 고유값 문제를 일련의 거친 체계 고유값 문제와 국소 경계값 문제로 분해하기 위해 다수준 이산화 프레임워크를 사용한다.
• 각 수준에서 국소적이고 병렬적인 유한요소 기법을 적용하여 경계값 문제를 독립적이고 동시적으로 해결한다.
• 최적의 수렴 순서를 유지하기 위해 다수준 보정 체계를 사용하여 해를 반복적으로 보정한다.
• 계산 작업이 프로세서당 최적임을 보장하여 문제 크기에 따라 효율적으로 스케일링된다.
• 선형 및 2차 유한요소를 모두 사용하여 구현하였으며, 정렬된 메쉬와 비정렬 메쉬에서 수치적 검증을 수행하였다.
• 다중 고유값에 대한 확장과 적응형 메쉬 정밀화 기법과의 결합을 지원한다.

실험 결과

연구 질문

• RQ1 다수준 국소 및 병렬 알고리즘이 유한요소 고유값 문제에 대해 최적 수렴성을 달성할 수 있는가?
• RQ2 제안된 방법이 기존의 유한요소 해법에 비해 계산 비용을 줄일 수 있는가?
• RQ3 해의 정확도를 유지하면서도 프로세서당 계산 작업을 최적화할 수 있는가?
• RQ4 가장 거친 메쉬 크기 $H$와 국소 패치 크기 $\delta$의 선택이 수렴성과 효율성에 어떤 영향을 미치는가?
• RQ5 이 방법은 다중 고유값과 비선형 고유값 문제로 확장될 수 있는가?

주요 결과

• 제안된 다수준 국소 및 병렬 알고리즘은 선형 유한요소의 경우 수렴 순서 2, 2차 유한요소의 경우 수렴 순서 4를 달성하여 이론적 기대와 일치한다.
• 수치 실험 결과, 첫 다섯 고유값의 오차 감소가 약 2.0의 최적 순서를 따르며, 메쉬 크기가 감소할수록 수렴 속도가 점차 2.0에 수렴하는 것으로 확인되었다.
• 첫 번째 고유함수의 경우 $H^1$-오차 감소는 약 0.9996의 최적 수렴 순서를 보이며, 거의 최적 수렴임을 나타낸다.
• 특히 2차 유한요소에서는 $\delta$가 작아져도 수렴 순서가 유지되며, 4차 수렴이 보존됨을 확인하였다.
• 프로세서당 계산 작업이 최적이며, 국소적이고 병렬적인 계산 특성 덕분에 복수의 프로세서에서 효율적으로 스케일링된다.
• 변동하는 확산 및 위치 에너지 항을 포함한 수치 예제를 통해 다중 고유값과 복잡한 계수 행렬에 대해서도 알고리즘이 효과적이고 정확함을 입증하였다.