[논문 리뷰] Solution Techniques for the Stokes System: A priori and a posteriori modifications, resilient algorithms

Markus Huber, Lorenz John외 4명|arXiv (Cornell University)|2015. 11. 18.

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

한 줄 요약

이 논문은 스토크스 시스템을 해결하는 저순서 유한요소법에 대해 사전 에너지 보정 및 사후 유량 보정 기법을 제안하며, 근사 공간을 확장하지 않고도 정확도와 질량 보존성을 향상시킨다. 이 방법들은 전역 오차 오염을 억제하고 국소 보존성을 보장하며, 내재된 다중격자 해법은 초대규모 계산 환경에서의 고장 내성 보장하여 복잡한 지구물리학적 맨틀 대류 시뮬레이션을 포함한 대규모 문제에서 최적 수렴성과 병렬 성능을 달성한다.

ABSTRACT

This article proposes modifications to standard low order finite element approximations of the Stokes system with the goal of improving both the approximation quality and the parallel algebraic solution process. Different from standard finite element techniques, we do not modify or enrich the approximation spaces but modify the operator itself to ensure fundamental physical properties such as mass and energy conservation. Special local a~priori correction techniques at re-entrant corners lead to an improved representation of the energy in the discrete system and can suppress the global pollution effect. Local mass conservation can be achieved by an a~posteriori correction to the finite element flux. This avoids artifacts in coupled multi-physics transport problems. Finally, hardware failures in large supercomputers may lead to a loss of data in solution subdomains. Within parallel multigrid, this can be compensated by the accelerated solution of local subproblems. These resilient algorithms will gain importance on future extreme scale computing systems.

연구 동기 및 목표

• 재입구 모서리 특이성으로 인한 스토크스 유동 시뮬레이션에서의 전역 오차 오염 문제를 유한요소 공간을 확장하지 않고 해결한다.
• 해법 복잡도를 증가시키지 않고 동일 순서 유한요소 방법에서 국소 질량 보존성을 달성한다.
• 미래 엑사스케일 시스템에서의 하드웨어 장애 가능성을 고려한 내재된 고장 내성 해법을 개발한다.
• 복잡한 지구물리학적 유동 문제에서 다중격자 해법의 최적 수렴성과 병렬 확장성을 보장한다.
• 운반과 스토크스 유동을 포함하는 다물리학 시뮬레이션에서 물리적으로 일관되고 보존적인 결합을 가능하게 한다.

제안 방법

• 재입구 모서리에서의 특이성을 고려해 강성 행렬에 국소 사전 수정을 적용하여 에너지 내적곱을 보정함으로써 이산 에너지 표현을 향상시킨다.
• 등순서 유한요소 근사에 대해 사후 유량 보정을 사용하여 국소 질량 보존성을 강제함으로써 안정화 또는 확장된 공간이 필요 없도록 한다.
• 보정된 방법을 Krylov 가속을 갖춘 기하학적 다중격자 프레임워크에 통합하여 최적의 대수적 수렴 속도를 확보한다.
• 장애가 발생한 부분 도메인에서 하위문제를 가속하여 해를 복구할 수 있도록 내재된 다중격자 알고리즘을 설계한다.
• 보정된 유한요소 방법을 명시적 시간 적분과 결합하여 유동-운반 문제에서의 확산-대류 방정식에 대한 유량 보정을 수행한다.
• 대규모 시스템에서 병렬 다중격자와 대수적 다중격자 예비조건화를 사용하여 해법을 구현함으로써 높은 병렬 효율성과 고장 내성을 확보한다.

실험 결과

연구 질문

• RQ1 재입구 모서리 특이성으로 인한 스토크스 유동에서의 전역 오차 오염은 유한요소 공간을 확장하지 않고 제거될 수 있는가?
• RQ2 등순서 유한요소 방법에서 해법 복잡도를 증가시키지 않고 국소 질량 보존성을 달성할 수 있는가?
• RQ3 초대규모 계산 환경에서 하드웨어 장애에 대한 다중격자 해법은 어떻게 내재화할 수 있는가?
• RQ4 사전 에너지 보정은 이산 시스템에서 최적의 수렴 속도와 개선된 에너지 표현을 이끌 수 있는가?
• RQ5 스토크스 방정식과 운반 방정식의 보존적 결합은 장기적인 지구물리학적 시뮬레이션의 물리적 일관성에 어떤 영향을 미치는가?

주요 결과

• 사전 에너지 보정은 재입구 모서리가 있는 영역에서 스토크스 시스템의 전역 오차 오염을 제거하며, 메esh 확장 없이도 최적 수렴 속도를 달성한다.
• 사후 유량 보정은 등순서 유한요소에서 국소 질량 보존성을 보장하며, 비보존적 방법에서 관찰된 최대 100%에 이르는 물리적으로 비합리적인 해의 과도한 상승 현상을 억제한다.
• 보정된 방법은 최적의 대수적 수렴성을 유지하며, 대규모 시스템에서도 뛰어난 병렬 확장성을 보인다.
• 내재된 다중격자 알고리즘은 하위문제 해를 가속함으로써 국소 데이터 손실 상황에서 복구할 수 있으며, 향후 엑사스케일 아키텍처에서 고장 내성 운영이 가능하다.
• 비동일 점성도 또는 온도 의존 점성도를 가진 지구물리학적 시뮬레이션에서는 보존적 결합이 비물리적 강제항을 방지하고 속도 및 온도장의 물리적 일관성을 유지한다.
• 1억 개의 테트라헤드론 메쉬에서, 복잡한 점성도 변화와 열대류를 포함한 맨틀 대류 시뮬레이션을 11,500단계에 걸쳐 안정적이고 유한한 해로 성공적으로 수행하였다.