[논문 리뷰] Symmetric parareal algorithms for Hamiltonian systems

Xing Dai, Bris, C. Le외 2명|arXiv (Cornell University)|2010. 11. 29.

Numerical methods for differential equations참고 문헌 32인용 수 10

한 줄 요약

이 논문은 하미르토니안 시스템에 대해 대칭성과 기하적 투영을 조합하여 장시간 시뮬레이션에서 에너지와 궤적 정확도를 유지하는 대칭 파라레얼 알고리즘을 소개한다. 대칭 시간 스텝과 에너지 다양체 위에 대한 거칠은 수준의 투영을 통합함으로써, 1000개의 프로세서를 사용할 때 순차적 해법에 가까운 정확도를 확보하면서도 약 66배의 속도 향상을 달성하였으며, 표준 파라레얼 알고리즘에 비해 안정성과 기하학적 정확도가 크게 향상되었다.

ABSTRACT

The parareal in time algorithm allows to efficiently use parallel computing for the simulation of time-dependent problems. It is based on a decomposition of the time interval into subintervals, and on a predictor-corrector strategy, where the propagations over each subinterval for the corrector stage are concurrently performed on the processors. In this article, we are concerned with the long time integration of Hamiltonian systems. Geometric, structure-preserving integrators are preferably employed for such systems because they show interesting numerical properties, in particular excellent preservation of the total energy of the system. Using a symmetrization procedure and/or a (possibly also symmetric) projection step, we introduce here several variants of the original plain parareal in time algorithm [Lions, Maday and Turinici 2001, Baffico, Bernard, Maday, Turinici and Zerah 2002, Bal and Maday 2002] that are better adapted to the Hamiltonian context. These variants are compatible with the geometric structure of the exact dynamics, and are easy to implement. Numerical tests on several model systems illustrate the remarkable properties of the proposed parareal integrators over long integration times. Some formal elements of understanding are also provided.

연구 동기 및 목표

• 하미르토니안 시스템에 적용했을 때 표준 파라레얼 알고리즘이 보이는 열악한 장시간 안정성과 에너지 보존 문제를 해결하기 위해.
• 하미르토니안 역학의 기하학적 구조, 특히 심플렉틱성과 시간 역행성과 같은 성질을 존중하는 파라레얼 변종을 개발하기 위해.
• 대칭화와 대칭 투영 단계를 결합하여 장시간 시뮬레이션에서 수렴성과 불변성 보존을 향상시키기 위해.
• 제안된 알고리즘이 하미르토니안 문제에서 순차적 정밀 해법과 유사한 정확도를 유지하면서도 상당한 병렬 속도 향상을 달성할 수 있음을 보여주기 위해.

제안 방법

• 예측-보정 단계를 대칭화하여 시간 역행성과 기하학적 구조 보존을 향상시키는 방식으로 표준 파라레얼 알고리즘을 변형한다.
• 각 시간 하위구간의 끝에서 거칠은 수준에서 해를 일정 에너지 다양체 위로 투영하는 대칭 투영 단계를 도입한다.
• 좋은 기하학적 성질을 지닌 정밀 해법(예: 심플렉틱 적분기)과 예측 단계를 위한 더 거칠고 단순한 해법을 사용한다.
• 계산 비용을 줄이기 위해 투영을 거칠은 수준에서만 수행하며, 하위구간당 1~2회 반복으로 충분하다.
• 정밀 보정을 하위구간들 간에 동시에 수행하면서도 대칭성과 기하학적 일致성을 유지함으로써 알고리즘이 병렬성을 유지하도록 한다.
• 대칭화와 투영을 결합하여 표준 파라레얼이 하미르토니안 시스템에서 애로를 겪는 공진 현상을 효과적으로 해결하고, 에너지 외의 불변성까지 안정화한다.

실험 결과

연구 질문

• RQ1 파라레얼 알고리즘의 대칭화는 하미르토니안 시스템에서 기하학적 구조 보존에 기여하는가?
• RQ2 거칠은 수준에서 에너지 다양체 위에 대칭 투영을 적용하면 장시간 안정성과 정확도가 향상되는가?
• RQ3 대칭화와 투영의 조합은 장시간 시뮬레이션에서 수렴성과 불변성 보존에 어떤 영향을 미치는가?
• RQ4 제안된 대칭 파라레얼 알고리즘은 정밀 순차 해법과 유사한 정확도를 유지하면서 어떤 정도의 속도 향상을 달성할 수 있는가?

주요 결과

• 기하학적 투영을 통한 대칭 파라레얼 알고리즘은 장시간 통합 시간 동안 순차적 정밀 해법과 유사한 궤적 정확도를 달성한다.
• 1000개의 프로세서를 사용할 경우, 약 66배의 속도 향상을 달성하였으며, 15회의 파라레얼 반복 내에 수렴하였다.
• 비선형 대칭 투영에 필요한 평균 반복 수는 하위구간당 단지 1.12회로, 계산 오버헤드가 매우 낮음을 시사한다.
• 알고리즘은 총 에너지 뿐 아니라 다른 불변성까지 안정화하여 장기 시뮬레이션의 신뢰성을 향상시킨다.
• 대칭화와 투영의 조합은 표준 파라레얼이 하미르토니안 시스템에서 애로를 겪는 공진 현상을 효과적으로 해결한다.
• 모델 시스템에 대한 수치 실험은 이 방법의 강건성과 장시간 통합에서 뛰어난 기하학적 성질을 확인한다.