[논문 리뷰] A spectral-infinite-element solution of Poisson's equation: an application to self gravity

Hom Nath Gharti, Jeroen Tromp|arXiv (Cornell University)|2017. 06. 02.

Geophysics and Gravity Measurements참고 문헌 22인용 수 13

한 줄 요약

이 논문은 지구의 내부에 대한 스펙트럴 요소와 외부 중력장에 대한 매핑된 무한 요소를 결합한 스펙트럴-무한 요소 방법을 제시한다. 이는 비유한 영역에서 푸아송 방정식을 정확하고 효율적으로 해결할 수 있도록 한다. 방법은 스펙트럴 요소에서는 가우스-레지안드르-로바토 적분을, 무한 방향에서는 가우스-라다우 적분을 사용하여 고정밀도를 달성하며, 균일한 구와 PREM 모델에 대해 분석적 해와의 비교를 통해 최소한의 오차로 검증된다.

ABSTRACT

We solve Poisson's equation by combining a spectral-element method with a mapped infinite-element method. We focus on problems in geostatics and geodynamics, where Earth's gravitational field is determined by Poisson's equation inside the Earth and Laplace's equation in the rest of space. Spectral elements are used to capture the internal field, and infinite elements are used to represent the external field. To solve the weak form of Poisson/Laplace equation, we use Gauss-Legendre-Lobatto quadrature in spectral elements inside the domain of interest. Outside the domain, we use Gauss-Radau quadrature in the infinite direction, and Gauss-Legendre-Lobatto quadrature in the other directions. We illustrate the efficiency and accuracy of the method by comparing the gravitational fields of a homogeneous sphere and the Preliminary Reference Earth Model (PREM) with (semi-)analytical solutions.

연구 동기 및 목표

• 지구역학 및 지구역학적 문제에 관련된 비유한 영역에서 푸아송 방정식을 해결하기 위한 고차수 수치 방법을 개발한다.
• 중력장 시뮬레이션에서 무한거리에 위치한 경계 조건을 수치적으로 구현하는 데 도전한다.
• 지구물리학적 문제인 빙하기 후의 빙하 복귀와 지진파 전파 등에서 자중중력과 연속체역학 간의 정확한 결합을 가능하게 한다.
• 균일한 구와 PREM 지구 모델에 대해 반분석적 해와의 비교를 통해 방법을 검증한다.
• 최소한의 계산 오버헤드로 병렬 구현을 통해 방법의 효율성과 안정성을 입증한다.

제안 방법

• 내부 영역에 대한 스펙트럴 요소와 외부 필드에 대한 매핑된 무한 요소를 결합하며, 좌표 상승 기법을 사용해 무한 영역을 유한한 계산 요소로 변환한다.
• 스펙트럴 요소에서는 고차수 정밀도를 확보하기 위해 가우스-레지안드르-로바토 적분을 적용한다.
• 무한 방향에는 가우스-라다우 적분을 사용하고, 다른 방향에는 일관성과 효율성을 확보하기 위해 가우스-레지안드르-로바토 적분을 사용한다.
• 노드 기반 스펙트럴 요소를 사용해 푸아송 및 라플라스 방정식의 약한 형태를 해석하며, 대각형 질량 행렬을 통해 시간 영역 시뮬레이션을 효율적으로 수행한다.
• 지구 모델과 무한 요소 영역 사이에 전이 층을 추가하여 표면 근처의 정확도를 향상시킨다.
• PETSc를 사용한 병렬 구현과 ParaView/VTK를 통한 시각화를 통해 확장성과 재현 가능성을 확보한다.

실험 결과

연구 질문

• RQ1 스펙트럴-무한 요소 방법은 3차원 비유한 영역에서 자중중력에 대한 푸아송 방정식을 정확하게 해결할 수 있는가?
• RQ2 균일한 구와 PREM 지구 모델에 대해 반분석적 해와 비교했을 때 방법의 성능은 어떠한가?
• RQ3 요소당 단지 세 개의 가우스-레지안드르-로바토 점만을 사용할 경우 해의 정확도와 안정성에 어떤 영향을 미치는가?
• RQ4 지구역학적 시뮬레이션에 있어 연속체역학 해석기와 효율적으로 결합될 수 있는가?
• RQ5 무한 요소 층을 추가함으로써 지구 표면 근처 및 외부 공간에서 중력장 계산의 정확도에 어떤 영향을 미치는가?

주요 결과

• 스펙트럴-무한 요소 방법은 균일한 구의 중력포텐셜에 대해 분석적 해와 뛰어난 일치를 보이며, 반경 방향 프로파일에서 오차가 0.1% 미만이다.
• PREM 모델에 대해 계산된 중력 가속도는 반분석적 해와 일치하며, 핵-마노 경계까지 선형 증가하고 지구 표면까지 거의 일정한 값을 나타낸다.
• 해결이 매끄럽고 고속 수렴하는 경우, 요소당 단지 세 개의 가우스-레지안드르-로바토 점만으로도 방법이 정확하게 유지된다.
• r = 4000 m로 연장된 모델 위에 무한 요소 층을 추가하면 분석적 해와 구분할 수 없을 정도로 일치하는 결과를 도출하며, 이는 외부 공간에서의 방법의 강건성을 확인한다.
• 방법은 안정적이고 효율적이며, 무한 요소 층에 대한 계산 비용이 극히 낮아 대규모 지구물리학적 시뮬레이션에 적합하다.
• 이 방법은 연속체역학 해석기와 직접 결합이 가능하여 반복 절차를 제거하고 장주기 지진 및 빙하 등화 조정 시뮬레이션을 지원한다.