[논문 리뷰] A consistent approach to unstructured mesh generation for geophysical models

Adam S. Candy|arXiv (Cornell University)|2017. 03. 24.

Hydrology and Watershed Management Studies참고 문헌 44인용 수 14

한 줄 요약

이 논문은 Shingle라는 자체 포함형 파이썬 라이브러리를 사용하여 지구물리학적 모델에서 일관되고 비의도적이며 반복 가능한 비정형 메시를 생성하는 완전 자동화된, 결정론적이고 반복 가능한 프레임워크를 제시한다. 수심 및 해안선과 같은 단일이고 일관된 지공간 데이터 소스에서 메시 생성을 유도함으로써, 재현 가능성, 데이터 기원 추적, 복잡한 다중 척도 영역 간의 원활한 통합을 보장하며, 인간의 실수를 크게 줄이고 엄격한 모델 상호 비교를 가능하게 한다.

ABSTRACT

Geophysical model domains typically contain irregular, complex fractal-like boundaries and physical processes that act over a wide range of scales. Constructing geographically constrained boundary-conforming spatial discretizations of these domains with flexible use of anisotropically, fully unstructured meshes is a challenge. The problem contains a wide range of scales and a relatively large, heterogeneous constraint parameter space. Approaches are commonly ad hoc, model or application specific and insufficiently described. Development of new spatial domains is frequently time-consuming, hard to repeat, error prone and difficult to ensure consistent due to the significant human input required. As a consequence, it is difficult to reproduce simulations, ensure a provenance in model data handling and initialization, and a challenge to conduct model intercomparisons rigorously. Moreover, for flexible unstructured meshes, there is additionally a greater potential for inconsistencies in model initialization and forcing parameters. This paper introduces a consistent approach to unstructured mesh generation for geophysical models, that is automated, quick-to-draft and repeat, and provides a rigorous and robust approach that is consistent to the source data throughout. The approach is enabling further new research in complex multi-scale domains, difficult or not possible to achieve with existing methods. Examples being actively pursued in a range of geophysical modeling efforts are presented alongside the approach, together with the implementation library Shingle and a selection of its verification test cases.

연구 동기 및 목표

• 지구물리학적 모델의 비정형 메시 생성에서 재현 불가능성과 일관성 부족 문제를 해결하기 위해, 종종 비공식적이고 수작업이 많으며 반복하기 어려운 방식이 지배하는 상황을 해결한다.
• 단일이고 자율적인 지공간 데이터 소스에서부터 경계를 따르는 비정형 메시를 생성하는 결정론적이고 자동화된 파이프라인을 개발한다.
• 모든 메시 생성 단계가 완전히 추적 가능하고 반복 가능한 방식으로 보장됨으로써, 엄격한 모델 상호 비교와 강력한 데이터 기원 추적을 가능하게 한다.
• 기존의 정형 격자로는 미세 구조를 포착하지 못하는 복잡한 다중 척도 지구물리학적 영역(예: 빙하 캐비티 및 해안 지역)을 지원한다.
• 수작업 메시 설계를 표준 지공간 도구 기반의 확장 가능한 소프트웨어 기반 프레임워크로 대체함으로써 계산 오버헤드와 인간 실수를 감소시킨다.

제안 방법

• 메시 생성의 기초로 레이터라이즈된 지오이드 또는 수심도(예: GEBCO)와 같은 일관된 데이터 소스를 사용한다.
• 도메인 경계를 원천 데이터로 정의된 바탕으로 삼아, 델로아네 삼각분할을 통해 2차원 비정형 메시(T_h)를 생성한다.
• 수평 메시의 구조를 유지하면서, 진행하는 표면 알고리즘을 사용해 수직으로 확장하여 3차원 체적 메시(T)를 생성한다.
• 전체 메시 생성 과정은 Shingle라는 파이썬 기반 라이브러리에 코딩되어 있으며, 결정론적 실행과 버전 제어된 처리를 통해 완전한 재현 가능성을 확보한다.
• 고정된 노드 집합을 통해 체적 보존을 유지하면서, 젖고 마르기 프로시저나 시간에 따라 변하는 표면 제약 조건을 통한 동적 경계 조정을 지원한다.
• 강력하고 확장 가능한 성능을 확보하기 위해 기존의 지공간 라이브러리(GDAL, PROJ.4, NetCDF)와 수치 도구(NumPy, SciPy)를 활용한다.

실험 결과

연구 질문

• RQ1 불규칙하고 프랙탈 같은 경계를 가진 복잡한 지구물리학적 영역에 대해 완전 자동화되고 결정론적이며 반복 가능한 메시 생성 파이프라인을 개발할 수 있는가?
• RQ2 비정형 메시 모델에서 모델 초기화 및 강제 조건 파rameter 설정 시 데이터 기원과 일관성을 어떻게 확보할 수 있는가?
• RQ3 단일이고 자율적인 데이터 소스가 비공식적이고 모델 전용의 메시 구축 방법을 얼마나 효과적으로 대체할 수 있는가?
• RQ4 제안된 방법은 다중 척도 지구물리학적 시뮬레이션에서 재현 가능성 향상과 엄격한 모델 상호 비교를 어떻게 향상시키는가?
• RQ5 고정 노드 메시 정밀화가 체적 보존과 장기적 시뮬레이션 안정성에 미치는 영향은 무엇인가?

주요 결과

• Shingle 라이브러리는 단일 지공간 데이터 소스에서부터 완전 자동화되고 결정론적이며 반복 가능한 비정형 메시 생성을 가능하게 하여 수작업 간섭을 제거하고 오류를 감소시킨다.
• 모든 메시 생성 단계를 버전 제어된 해석 가능한 파이썬 워크플로우에 코딩함으로써 데이터 기원과 일관성을 보장한다.
• 이 방법으로 생성된 메시는 원본 수심도 및 해안선 데이터와 일관되며, 델로아네 삼각분할과 수직 확장을 통해 경계 일치성이 유지된다.
• 정형 격자가 미세 구조를 포착하지 못하는 복잡한 다중 척도 영역(예: 빙하 캐비티 및 해안 지역)을 지원한다.
• 수평 메시의 노드를 고정하고 선형 보간을 사용한 정밀화를 통해, 동적 경계 조건이 존재하는 상황에서도 체적 보존을 유지한다.
• 모든 제약 조건과 처리 단계가 완전히 문서화되고 재현 가능하므로, 효율적인 모델 상호 비교와 데이터 공유를 가능하게 한다.