[논문 리뷰] Effects of Horizontal Discretization on Triangular and Hexagonal Grids on Linear Baroclinic and Symmetric Instabilities

Steffen Maaß, Sergey Danilov|arXiv (Cornell University)|2026. 01. 07.

Oceanographic and Atmospheric Processes인용 수 0

한 줄 요약

논문은 삼각 및 육각 격자(B 및 C 스태거)에서의 수평 이산화가 선형 바로클리닉 및 대칭 불안정성에 미치는 영향을 분석하여 격자 기하학과 관련된 잘못된 모드가 나타나고 이를 완화하기 위해 점도/확산을 조정해야 함을 보여준다.

ABSTRACT

As global ocean general circulation models are run at eddy-permitting resolutions, reproducing accurate growth rates of baroclinic instabilities is a major concern when choosing a discretization of the equations of motion. From this viewpoint, we analyze discretizations on triangular and hexagonal grids with different types of variable staggering used in several ocean circulation models. By extending the linear baroclinic instability analysis in the Eady configuration to discretizations on more complex grids, several numerical subtleties are revealed. In comparison to discretizations on quadrilateral grids, the analyzed discretizations are less robust against unstable spurious modes, partly created by the mesh geometry. Some of the subtleties arise because spurious modes on staggered triangular and hexagonal grids do not adhere to Galilean invariance. As a consequence, their growth rates demonstrate a dependence on the alignment between the background flow and the grid, as well as the strength of a uniform background flow. The interactions with spurious modes become more significant on the axis of symmetric instabilities where the physical and spurious branches of instability are more difficult to separate in wavenumber space. Our analysis shows that in most cases moderate biharmonic viscosity and diffusion suppress spurious branches. However, one needs to carefully calibrate the viscosity and diffusivity parameters for each of the considered discretizations in order to achieve this.

연구 동기 및 목표

• 삼각 및 육각 격자 이산화(B 및 C 스태거)가 Eady 구성에서 선형 바로클리닉 및 대칭 불안정성을 재현하는 방법을 평가한다.
• 격자 기하학 및 스태거 차이로 인해 발생하는 수치적 미묘함과 잘못된 모드를 식별한다.
• 다양한 이산화에서 잘못된 불안정성 가지를 억제하는 데 있어 점도와 확산의 효과를 평가한다.

제안 방법

• 삼각 및 육각 격자에 대한 선형 바로클리닉 불안정성 분석(Eady 문제)을 확장한다.
• 비격자(unstructured) 메쉬에서 B 및 C 격자에 대한 이산적인 수평 연산자와 푸리에 기호를 개발한다.
• 다양한 전달(advective) 스킴(FDV, FDCRE, AVI, ASC)의 비교와 그것이 안정성 및 잘못된 모드에 미치는 영향을 분석한다.
• 경사 격자 고려 및 재구성 스킴에 따른 부유도(부유) 전파를 분석한다.
• 완전한 분석을 위한 hexagonal C 격자와 triangular A 격자에 대한 푸리에 기호 분석 보조 노트북을 제공한다.

실험 결과

연구 질문

• RQ1 삼각형 B- 및 C-격자 이산화가 Eady 구성에서 바로클리닉 및 대칭 불안정성의 성장률을 어떻게 재현하는가?
• RQ2 격자 기하학과 스태거로 인해 삼각형/육각형 격자에서 어떤 잘못된 모드가 발생하며, 이는 배경 흐름 정렬에 따라 어떻게 달라지는가?
• RQ3 다양한 이산화에서 점도와 확산이 잘못된 불안정성 가지를 어느 정도까지 억제하는가?
• RQ4 다양한 수평 운동량 및 부유도 이산화가 잘못된 모드의 존재 및 거동에 어떤 영향을 미치는가?

주요 결과

• 스태거된 삼각 및 육각 격자에서 속도 및 스칼라 자유도 간의 불균형과 격자 기하학으로 인해 잘못된 모드가 발생한다.
• 잘못된 모드는 배경 흐름과 격자의 정렬 및 균일한 배경 흐름의 세기에 의존할 수 있으며, 이는 일부 스킴에서 갈리안 불변성의 손실을 시사한다.
• 적당한 바이하모닉 점도와 확산은 일반적으로 이산화 전반에 걸쳐 잘못된 불안정성 가지를 억제하지만, 격자별로 보정이 필요하다.
• 일부 advective 스킴(FDCRE, ASC)은 물리적 가지의 정확성을 높이지만 반대 방향의 잘못된 기하 모드를 유발할 수 있다.
• 잘못된 모드와의 상호 작용은 대칭 불안정성 축 근처에서 물리적 가지와 잘못된 가지가 같은 파수 공간에서 겹치는 영역에서 더 크게 나타난다.
• A-그리드(삼각형) 배열은 잘못된 불안정에 덜 취약하며, 육각형 및 삼각형 B/C 격자는 신중한 처리가 필요하다.