[논문 리뷰] Compatible finite element methods for numerical weather prediction

Colin J. Cotter, Andrew T. T. McRae|arXiv (Cornell University)|2014. 01. 03.

Advanced Numerical Methods in Computational Mathematics참고 문헌 6인용 수 4

한 줄 요약

이 논문은 유한요소외삽계산 기반의 호환 가능한 유한요소 방법을 제안하여, 비정규 구형 메esh에서 회전하는 얕은 수심 방정식의 안정적이고 보존적이며 정확한 이산화를 가능하게 한다. 주요 기여는 정확한 지구적 모드와 함께 에너지 및 에너지분산량을 보존하는 증명 가능한 안정성 있는 2차 수렴 형식을 제공하는 것으로, 이는 이코사헤드럴 메쉬에서 윌리엄슨 테스트 케이스를 통해 입증되었다.

ABSTRACT

This article takes the form of a tutorial on the use of a particular class of mixed finite element methods, which can be thought of as the finite element extension of the C-grid staggered finite difference method. The class is often referred to as compatible finite elements, mimetic finite elements, discrete differential forms or finite element exterior calculus. We provide an elementary introduction in the case of the one-dimensional wave equation, before summarising recent results in applications to the rotating shallow water equations on the sphere, before taking an outlook towards applications in three-dimensional compressible dynamical cores.

연구 동기 및 목표

• 호환 가능한 유한요소를 사용하여 수치 날씨 예측을 위한 안정적이고 보존적인 유한요소 방법을 개발한다.
• 유한요소 이산화에서의 가짜 모드 분지 문제를 해결하기 위해 속도와 압력 자유도 간의 적절한 균형을 확보한다.
• 유한요소외삽계산을 활용하여 C-그리드 유한차분의 유사성 성질을 비정규 및 비수직 메쉬로 확장한다.
• 구상에서 얕은 수심 방정식에 대한 에너지 및 에너지분산량을 보존하는 형식을 가능하게 한다.
• 호환 가능한 유한요소 공간을 사용하여 3차원 압축성 역학 코어의 기초를 마련한다.

제안 방법

• 벡터 미적분학의 미분 구조를 유지하는 호환 가능한 유한요소 공간(예: BDFM1-DG1, RTk-DG1)을 사용하여, 컬 및 디브의 이산 대응식을 보장한다.
• 속도의 연속 법선 성분과 압력의 비연속 성분을 포함하는 혼합 유한요소 설정을 적용하여 안정성과 일致성을 확보한다.
• 변분 형식을 통해 약한 해를 유도하고, 부분적 통합과 잠재성 운동량의 일致성 투영을 통해 보존 법칙을 강제한다.
• 비선형 항에 대해 뉴턴 반복을 사용하는 반의존적 시간 이산화를 구현하여 안정성과 정확도를 확보한다.
• 수렴성과 해 품질 테스트를 위해 적응적 정밀도를 가진 이코사헤드럴 구형 메쉬를 사용한다.
• 에너지 보존을 유지하면서도 격자 척도에서 에너지분산량을 소산시키기 위해 잠재성 운동량의 이동을 수정한다. 이는 격자 척도 진동을 방지한다.

실험 결과

연구 질문

• RQ1 호환 가능한 유한요소 방법은 비정규 구형 메쉬에서 얕은 수심 방정식에 대해 안정적이고 보존적이며 정확한 해를 제공할 수 있는가?
• RQ2 속도와 압력 간의 전역 자유도 비율이 가짜 관성중력 또는 로스비 모드의 존재에 어떤 영향을 미치는가?
• RQ3 C-그리드 유한차분 방법과 유사하게 호환 가능한 유한요소를 사용하여 에너지 및 에너지분산량을 보존하는 형식을 구성할 수 있는가?
• RQ4 윌리엄슨 테스트 케이스 5와 같은 벤치마크 테스트 케이스에 적용했을 때 호환 가능한 유한요소 방법의 수렴 속도와 정확도는 어떠한가?
• RQ5 적절한 유한요소 공간(속도, 밀도, 잠재온도)을 사용하여 이 프레임워크를 3차원 압축성 유동으로 확장할 수 있는가?

주요 결과

• BDFM1-DG1 유한요소 쌍은 이코사헤드럴 메쉬의 정밀화에 따라 윌리엄슨 테스트 케이스 5에서 높이 필드에 대해 2차 수렴을 달성한다. 이는 정밀화된 메쉬에서의 수렴 플롯을 통해 확인되었다.
• 속도와 압력 간의 2:1 전역 자유도 비율 덕분에 지구적 모드를 정확히 유지한다. 이는 수치 날씨 예측에 있어 핵심적인 성질이다.
• 격자 척도에서 에너지분산량 소산을 포함한 에너지 및 에너지분산량 보존 형식은 빈도 합성 시뮬레이션에서 격자 척도 진동을 효과적으로 방지한다.
• 잠재성 운동량 보존 법칙은 이산 변분 형식에 일致성 있게 투영되어, 이산 시스템이 기본 물리 법칙을 존중함을 보장한다.
• 15일째의 속도 및 높이 필드는 산맥을 통과하면서 현실적인 진화를 보이며, 물리적으로 비합리적인 진동이 없고, 잠재성 운동량 필드는 50일 동안 잘 해상도를 유지한다.
• 이 프레임워크는 3차원으로 확장 가능하며, 미분 연산자를 유지하는 네 공간 구조(V0부터 V3까지)를 제안하여 안정적이고 보존적인 이산화를 가능하게 한다.