[논문 리뷰] An Efficient Multigrid Method for Ground State Solution of Bose-Einstein Condensates

Hehu Xie, Fei Xu외 1명|arXiv (Cornell University)|2017. 11. 04.

Advanced Numerical Methods in Computational Mathematics참고 문헌 20인용 수 5

한 줄 요약

이 논문은 유한요소를 통해 고르스-피타예프스키 방정식을 해결함으로써 보즈-아인슈타인 응축물의 기본 상태를 계산하기 위한 효율적인 다계층 방법을 제안한다. 비선형 반복 구현 방식이 새로 도입되어 비선형성에 관계없이 渐近적으로 최적의 계산 작업을 달성하며, 이는 비선형성에 의존하지 않는 최적 수렴성과 거의 선형 복잡도를 유지한다. 이는 다계층 방법으로 선형 문제를 해결할 때와 동일한 효율성을 보인다.

ABSTRACT

An efficient multigrid method is proposed to compute the ground state solution of Bose-Einstein condensations by the finite element method based on the combination of the multigrid method for nonlinear eigenvalue problem and an efficient implementation for the nonlinear iteration. The proposed numerical method not only has the optimal convergence rate, but also has the asymptotically optimal computational work which is independent from the nonlinearity of the problem. The independence from the nonlinearity means that the asymptotic estimate of the computational work can reach almost the same as that of solving the corresponding linear boundary value problem by the multigrid method. Some numerical experiments are provided to validate the efficiency of the proposed method.

연구 동기 및 목표

• 보즈-아인슈타인 응축물에서 고르스-피타예프스키 방정식(GPE)의 기본 상태 해를 계산하기 위한 더 효율적인 수치적 방법을 개발하기 위해.
• 기존의 GPE 고유값 문제에 대한 다계층 방법에서 비선형성에 의존하는 계산 비용 문제를 해결하기 위해.
• 비선형성 매개수 ζ에 거의 의존하지 않는 渐近적으로 최적의 계산 복잡도를 달성하기 위해.
• 다계층 프레임워크에 효율적인 비선형 반복 기법을 통합하여 확장성과 성능을 향상시키기 위해.

제안 방법

• 비차원화된 GPE 고유값 문제에 대해 유한요소 이산화를 적용하며, 이는 이차형식 a(⋅,⋅)와 b(⋅,⋅)를 포함한 약한 형태를 사용한다.
• 일치보정단계(one-correction-step) 다계층 알고리즘을 적용하며, 보조적인 내재된 유한요소 공간 계층에서 보정을 계산한다.
• 각 수준에서 비선형 시스템을 분리하고 가속화하기 위해 텐서 기반 연산을 사용한 비선형 반복의 효율적 구현을 도입한다.
• 이전 연구에서 제안된 다수준 보정 기법을 활용하지만, 반복적인 고비용 비선형 해법을 피하기 위해 특수화된 비선형 해법을 통합하여 향상시킨다.
• 특성가소성 메쉬 세분화를 위해 잔차 기반 후행 오차 추정기를 사용하여 비凸 영역 내의 특이점을 다룬다.
• 계산 작업을 분석하여 거의 최적의 스케일링을 보이며, 다계층 방법으로 선형 문제를 해결할 때와 유사한 복잡도에 도달한다.

실험 결과

연구 질문

• RQ1 다계층 방법을 사용하여 비선형 GPE 고유값 문제를 해결할 때, 비선형성 강도 ζ에 거의 의존하지 않는 계산 비용을 달성할 수 있는가?
• RQ2 GPE에 대한 다계층 방법에서 비선형 반복을 효율적으로 구현하면 고비용의 계산 오버헤드를 피할 수 있는가?
• RQ3 제안된 방법은 균일 메쉬와 적응형 메쉬 환경에서 모두 최적 수렴률을 유지하면서 渐近적으로 최적의 작업 복잡도를 달성하는가?
• RQ4 이 방법은 나비에-스토크스 방정식이나 단계장 모델과 같은 다항식 비선형성을 갖는 다른 문제로 확장 가능한가?

주요 결과

• 제안된 다계층 방법은 비선형성 매개수 ζ에 거의 의존하지 않는 渐近적으로 최적의 계산 작업을 달성하며, 해당 선형 문제를 해결할 때의 복잡도에 가까워진다.
• 1D, 2D, 3D에서의 수치 실험 결과, 새로운 방법의 CPU 시간은 메쉬 크기와 거의 선형적으로 증가하지만, 큰 ζ 값에서도 원래 방법([26] 참조)이 강한 ζ 의존성을 보이는 것과 대비된다.
• 비볼록 영역인 L자형 도메인에서 적응형 메쉬 세분화를 사용할 경우, 최적 수렴률과 정확도를 유지하며 특이성에 대한 강건성을 확인한다.
• 오차 추정치와 수렴 순서는 이론적 예측과 일치하며, 후행 오차 추정기는 적응형 메쉬 세분화를 효과적으로 이끌어낸다.
• 새로운 방법의 계산 작업은 특히 고비선형성에서 기존 다계층 방법([26])보다 크게 낮아 효율성 향상을 확인한다.
• 이 방법은 나비에-스토크스 방정식 및 단계장 모델과 같은 다항식 비선형성을 갖는 다른 비선형 문제로도 확장 가능하며, 일반 비선형 고유값 문제의 전처리기로 활용될 수 있다.