High order tracer variance stable transport with low order energy conserving dynamics for the thermal shallow water equations

기술 요약: 저차 에너지 보존 역학과 결합된 고차 트레이서 분산 안정 수송

문제 제기
지구물리 유체 역학의 수치 모델링은 종종 계산 효율성과 물리적 충실도 사이의 절충을 요구합니다. 저차 수치 기법 (유한 차분, 유한 체적, 유한 요소) 은 조건수가 낮아 빠른 수렴과 행렬 조립을 용이하게 하므로 암시적 대기 솔버에 빈번히 활용됩니다. 그러나 이러한 기법들은 고주파수 아티팩트를 억제하기 위해 업윈딩 (upwinding) 에 의존하는 경향이 있어, 과도한 수치 감쇠를 초래하며 난류 특성의 해상도를 저하시킵니다. 반면, 불연속 갈레르킨 (Discontinuous Galerkin) 과 같은 고차 방법은 관성 운동을 해결하고 정확한 수치 플럭스를 구성하는 데 효과적이지만, 계산 비용이 많이 들며 종종 저차 암시적 솔버의 견고한 안정성 특성을 결여하고 있습니다. 또한, 표준 고차 수송 기법들은 불안정성을 유발하지 않고 트레이서 분산 (시스템의 수치 엔트로피) 을 보존하는 데 어려움을 겪는 반면, 저차 기법들은 비정준 해밀토니안 시스템에서 장기 안정성에 필요한 에너지 보존 구조를 유지하지 못하는 경우가 많습니다.

방법론
저자들은 열 얕은 물 방정식의 맥락에서, 역학을 위한 저차 혼합 유한 요소 솔버와 열역학적 트레이서의 물질 수송을 위한 고차 불연속 갈레르킨 (DG) 기법을 결합한 하이브리드 수치 프레임워크를 제안합니다.

• 저차 역학: 역학 솔버는 방정식의 에너지 보존 구조를 유지하는 혼합 유한 요소 형식 (최저 차수 다항식) 을 활용합니다. 이 솔버는 조밀하지 않은 메쉬에서 작동하며 속도 ($u$), 깊이 ($h$), 그리고 깊이 가중 부력 ($S$) 에 대한 예측 변수를 계산합니다.
• 고차 수송: 부력 수송은 중첩되고 정제된 메쉬에서 작동하는 고차 DG 기법에 의해 처리됩니다. 부력은 플럭스 변수 ($S=sh$) 가 아닌 물질 변수 ($s$) 로 취급됩니다. 이 기법은 이산 형식이 트레이서 분산을 명백히 보존할 수 있도록 연속 수준에서 대칭적 확장 (skew-symmetric expansion) 을 이류 연산자에 적용합니다.
• 결합 전략: 두 솔버는 중첩된 메쉬 계층을 통해 결합됩니다. 저차 메쉬는 고차 수송 메쉬 내에 내장되어, 저차 요소의 셀 중심이 고차 요소의 가우스 - 르장드르 구적점과 일치하도록 합니다. 이는 자유도 간 1:1 매핑을 생성합니다.
  • 고차 부력장은 각 시간 단계 시작 시 저차 예측 변수로부터 진단됩니다.
  • 고차 수송 방정식은 시간 단계 동안 평균화된 저차 질량 플럭스와 깊이 장을 사용하여 진행됩니다.
  • 결과적으로 얻어진 고차 부력은 역학 솔버를 업데이트하기 위해 저차 공간으로 투영됩니다.
• 보존 특성:
  • 에너지: 결합은 플럭스 형 부력 수송과 압력 구배의 반대칭성을 유지하여 저차 역학에 대한 이산 에너지 보존을 보장합니다.
  • 트레이서 분산: 고차 수송 기법은 공간적으로 트레이서 분산을 보존하도록 설계되었습니다. 안정성을 보장하고 비선형 에일리어싱으로 인한 그리드 규모 노이즈를 방지하기 위해 업윈딩 항이 포함되며, 이는 트레이서 분산의 무제한 증가가 아닌 제어된 감쇠 (댐핑) 를 허용합니다.
• 시간 적분: 시스템은 수송 단계에 대해 강하게 안정된 3 차 룽게 - 쿠타 (SSP-RK3) 적분기를 사용하고, 결합된 역학에 대해서는 (GMRES 를 통해 해결되는) 준-뉴턴 접근법을 사용합니다.

주요 기여

1. 하이브리드 형식: 저자들은 저차 에너지 보존 역학과 고차 트레이서 분산 보존 수송을 조화시키는 일관된 형식을 도입했습니다.
2. 물질 형 수송: 이전의 플럭스 형에 의존하던 트레이서 분산 보존 형식과 달리, 본 연구는 이산 이류 연산자를 물질 형으로 표현하여 저차 플럭스 형 역학과 효과적으로 결합합니다.
3. 메쉬 계층: 저차 셀 중심이 고차 가우스 - 르장드르 점과 정렬되는 특정 메쉬 계층의 구현은 복잡한 보간 오차 없이 직접적이고 효율적인 결합을 가능하게 합니다.
4. 안정성 메커니즘: 이 방법은 업윈딩을 사용한 고차 수송이 순수 저차 방법보다 더 풍부한 난류 구조를 보존하면서도 저차 업윈딩 특유의 과도한 감쇠 없이 안정적인 해를 제공함을 보여줍니다.

결과
이 방법은 세 가지 표준 테스트 사례를 통해 검증되었습니다:

1. 고체 회전: 고차 DG 기법은 업윈딩된 플럭스에 대해 4 차 수렴을, 중심 플럭스에 대해 3 차 수렴을 보여주었습니다. 질량 보존은 정확했으며, 트레이서 분산 보존 오차는 해상도에 따라 감소하여 기법의 안정성을 확인했습니다.
2. 열지오스트로픽 평형: 정상 상태 테스트에서 전체 결합 시스템은 저차 역학에 의해 제한된 모든 변수에 대해 2 차 수렴을 보였습니다. 이 기법은 질량, 깊이 가중 부력, 그리고 에너지에 대해 기계 정밀도 수준의 보존을 유지했으며, 상대 와도에서 장기적인 드리프트가 없었습니다.
3. 열적 불안정성 (비선형 난류): 매우 비선형적인 영역에서 업윈딩을 포함한 고차 수송 기법은 미세 구조를 보존하는 안정적인 해를 생성했습니다.
   • 비교: 순수 저차 중심 플럭스는 제어되지 않은 그리드 규모 에일리어싱으로 인해 불안정해졌습니다. 순수 저차 업윈딩 플럭스는 과도하게 확산되어 난류 구조를 흐리게 만들었습니다.
   • 하이브리드 성능: 결합된 고차/저차 기법은 안정성을 유지하면서 (트레이서 분산이 적당히 감소) 난류 해의 풍부함을 보존하여 두 가지 순수 저차 대안보다 우수한 성능을 발휘했습니다.

의의 및 주장
저자들은 이 접근법이 저차 암시적 솔버의 계산 효율성과 고차 수송 기법의 물리적 충실도 사이의 간극을 성공적으로 연결한다고 주장합니다. 시스템의 전체 정확도는 저차 역학에 의해 공식적으로 2 차로 제한되지만, 고차 트레이서 분산 보존 수송의 사용은 모델의 안정성과 난류 해의 보존을 크게 향상시킵니다.

저자들은 현재 구현이 기하학적 멀티그리드 방법에 비해 매우 대규모 병렬 도메인에서 효율적으로 확장되지 않을 수 있는 GMRES 솔버를 사용한다는 점을 겸손하게 지적하며, 이를 향후 확장의 잠재적 영역으로 제안합니다. 저자들은 결합 전략이 건조 및 습윤 대기를 위한 압축성 오일러 방정식과 같은 열역학적 트레이서를 포함하는 다른 비정준 해밀토니안 시스템에도 일반화될 수 있다고 주장합니다. 주요 성과는 에너지 보존과 트레이서 분산 안정성이 전통적인 저차 대기 모델에서 흔히 발견되는 과도한 감쇠나 불안정성 없이 결합 시스템에서 동시에 유지될 수 있음을 입증한 것입니다.