[논문 리뷰] Galerkin Least-Squares Stabilization in Ice Sheet Modeling - Accuracy, Robustness, and Comparison to other Techniques

Christian Helanow, Josefin Ahlkrona|arXiv (Cornell University)|2017. 02. 27.

Cryospheric studies and observations참고 문헌 44인용 수 6

한 줄 요약

이 연구는 빙하류 유한요소 모델링에서 갈레르킨 최소제곱(GLS) 안정화 방법을 평가하며, 내부 페널티 및 압력 페널티와 같은 다른 방법들과 비교한다. GLS는 안정화 파arameter 선택에 민감하며, 특히 수직 속도에 대해 빙하 경계 근처에서 진동을 유발함을 발견하였다. 내부 페널티 방법은 복잡한 기하구조에서 더 뛰어난 안정성과 정확도를 보이며, 특히 우수한 성능을 발휘한다.

ABSTRACT

We investigate the accuracy and robustness of one of the most common methods used in glaciology for the discretization of the $\mathfrak{p}$-Stokes equations: equal order finite elements with Galerkin Least-Squares (GLS) stabilization. Furthermore we compare the results to other stabilized methods. We find that the vertical velocity component is more sensitive to the choice of GLS stabilization parameter than horizontal velocity. Additionally, the accuracy of the vertical velocity component is especially important since errors in this component can cause ice surface instabilities and propagate into future ice volume predictions. If the element cell size is set to the minimum edge length and the stabilization parameter is allowed to vary non-linearly with viscosity, the GLS stabilization parameter found in literature is a good choice on simple domains. However, near ice margins the standard parameter choice may result in significant oscillations in the vertical component of the surface velocity. For these cases, other stabilization techniques, such as the interior penalty method, result in better accuracy and are less sensitive to the choice of the stabilization parameter. During this work we also discovered that the manufactured solutions often used to evaluate errors in glaciology are not reliable due to high artificial surface forces at singularities. We perform our numerical experiments in both FEniCS and Elmer/Ice.

연구 동기 및 목표

• 등차수 유한요소 방법에 Galerkin Least-Squares(GLS) 안정화를 적용한 빙하 시뮬레이션의 정확도와 안정성을 평가하는 것.
• 특히 수직 속도에서 발생하는 수치 오차가 빙하 표면 위치와 부피 예측으로 어떻게 전파되는지 조사하는 것.
• 기본 벤치마크 및 실제 빙하 문제에서 GLS 안정화를 다른 안정화 방법—압력 페널티, 내부 페널티, 압력 글로벌 프로젝션, 국소 프로젝션—과 비교하는 것.
• 빙하학적 유한요소 모델에서 오차 평가를 위한 제작된 해법(manufactured solutions)의 신뢰성 평가.
• 실제 빙하학적 환경에서 표준 GLS 형식의 최적 안정화 파arameter 도출 및 제한 사항 규명

제안 방법

• FEniCS와 Elmer/Ice 두 유한요소 코드를 사용하여 수치 실험을 수행하여 결과의 상호 검증을 가능하게 한다.
• GLS 안정화는 등차수 선형 유한요소를 사용한 p-Stokes 방정식에 적용되며, 요소 크기와 점성도에 따라 변하는 안정화 파arameter는 Franca와 Frey(1992)의 결과에서 유도된다.
• 두 가지 오차 측정 전략을 사용한다: 수렴 테스트를 위한 제작된 해법과 2차원에서의 검증을 위한 고해상도 기준 해.
• 성능 평가를 위해 ISMIP-HOM 벤치마크와 실제 출구 빙하 설정을 사용하여 메쉬 정밀도 및 점성도 변화 조건에서의 성능을 시험한다.
• 동일한 조건에서 내부 페널티, 압력 페널티, 압력 글로벌 프로젝션, 국소 프로젝션 등의 대안 안정화 기법을 구현하고 비교한다.
• 메쉬 정밀도 연구는 수직 및 수평 수렴을 중점으로 하며, 특히 특이점과 빙하 경계 근처에서 수직 속도 성분의 거동을 분석한다.

실험 결과

연구 질문

• RQ1 GLS 안정화 파arameter 선택이 빙하 모델의 수직 속도 및 표면 위치 정확도에 어떤 영향을 미치는가?
• RQ2 제작된 해법이 빙하학적 유한요소 시뮬레이션에서 신뢰할 수 있는 오차 추정치를 제공하는가?
• RQ3 GLS, 내부 페널티, 압력 페널티 등 다양한 안정화 유한요소 방법 간 정확도, 안정성, 안정화 파arameter 민감도 측면에서의 비교는 어떻게 되는가?
• RQ4 수직 속도 성분의 수치 오차가 장기적인 빙하 표면 변화 및 부피 예측에 어떤 영향을 미치는가?
• RQ5 GLS 안정화가 비물리적 진동을 유도하는 상황는 언제이며, 내부 페널티와 같은 대안 방법이 이러한 문제를 완화할 수 있는가?

주요 결과

• 수직 속도 성분은 수평 성분보다 GLS 안정화 파arameter 선택에 훨씬 더 민감하며, 파arameter가 잘못 선택되면 빙하 경계 근처에서 진동이 발생한다.
• Franca와 Frey(1992)에서 제안한 표준 GLS 안정화 파aram터는 세포 크기를 최소 모서리 길이로 정의하고 점성도가 비선형적으로 변할 경우 단순 도메인에서는 잘 작동하지만, 실제 출구 빙하 시뮬레이션에서는 실패한다.
• 제작된 해법은 특이점에서 인위적인 표면 힘을 유도하여 오차 추정치가 과도하게 높아지고 수직 방향으로 초수렴 현상이 발생함으로써, 모델 검증에 대한 신뢰성을 저하시킨다.
• 내부 페널티 방법은 GLS에 비해 더 뛰어난 안정성과 파arameter 선택에 대한 낮은 민감도를 보이며, 특히 높은 점성도 대비가 있는 복잡한 기하구조에서 유리하다.
• 압력 글로벌 프로젝션 방법은 정확도는 높지만, 대규모 빙하 시뮬레이션에는 계산적으로 비현실적이다.
• 국소 프로젝션 안정화 방법은 GLS와 유사한 성능을 보이나 구현이 더 복잡하며, 비선형 변형형은 향상된 성능을 제공할 수 있으나 추가 연구가 필요하다.