Spectral interpolation in semi-implicit semi-Lagrangian methods for shallow water equations on the sphere

기술 요약: 구면에서의 얕은 물 방정식을 위한 준암시적 준라그랑주 방법에서의 스펙트럴 보간

문제 제기
준암시적 준라그랑주 (SISL) 방법은 코란트 - 프리드리히스 - 레위 (CFL) 조건이 허용하는 것보다 큰 시간 간격을 사용할 수 있기 때문에 운영 전 지구 기상 및 기후 모델의 표준으로 자리 잡았습니다. 그러나 현재 운영 체계에는 상당한 수치적 불일치가 존재합니다. 예측 변수의 도함수에 대한 공간 이산화는 일반적으로 스펙트럴 정확도를 가진 방법 (구면 조화 변환 [SHT] 또는 더블 푸리에 구 [DFS] 방법 등) 을 사용하는 반면, 격자 밖 이탈 지점에서 변수를 평가하는 데 필요한 보간 단계는 3 차 또는 5 차 텐서 곱 라그랑주 다항식과 같은 저차 기법에 의존합니다.

이러한 불일치는 수치 정확도에서 불일치를 초래합니다. 저차 보간은 종종 계산 비용 문제 (스펙트럴 보간의 $O(n^4)$에 비해 $O(n^2)$의 복잡도) 와 단순한 대류에 대한 이론적 오차 경계로 정당화되지만, 얕은 물 방정식 (SWE) 과 같은 복잡한 시스템에 대한 이러한 절충을 정당화하는 제한된 수치적 증거가 문헌에 존재합니다. furthermore, 저차 보간은 과도한 수치 확산과 위상 오차를 유발하는 것으로 알려져 있어, 장기간의 적분 시간 동안 정확도와 보존 특성을 저하시킬 수 있습니다.

방법론
본 연구는 구면에서의 SWE 를 위한 SISL 체계 내에서 스펙트럴 정확도를 가진 보간에 대한 최초의 수치적 조사를 제시합니다. 저자들은 안정된 푸리에 기저와 2 시간 단계 SETTLS 체계를 사용하는 요시무라 [56] 가 개발한 DFS 기반 SISL 모델을 기반으로 연구를 진행했습니다.

핵심 방법론적 혁신은 공간 도함수에 사용된 동일한 DFS 푸리에 기저에 기반한 스펙트럴 보간 체계로 표준 라그랑주 보간을 대체하는 것입니다. 계산 비용 장벽을 해결하기 위해 저자들은 비균일 고속 푸리에 변환 (NUFFT) 을 활용합니다.

• 보간 메커니즘: 예측 변수 (특히 준암시적 시스템의 $R_v$ 및 $R_h$ 항) 는 부분 규칙성 DFS 기저를 사용하여 전개됩니다. 이러한 전개는 비격자 이탈 지점에서 평가됩니다.
• 복잡도 감소: $O(n^2)$개의 이탈 지점에서 이변수 푸리에 급수를 직접 평가하는 것은 $O(n^4)$의 연산이 필요합니다. 급수를 복소 지수 형태로 재구성하고 NUFFT (특히 FINUFFT 라이브러리) 를 활용함으로써, 계산 복잡도는 $O(n^2 \log^2(1/\epsilon) + n^2 \log n)$로 감소됩니다. 여기서 $\epsilon$은 정확도 허용 오차입니다. 이로써 보간 단계는 나머지 스펙트럴 SISL 체계와 계산적으로 일관되게 됩니다.
• 구현: 저자들은 두 가지 맥락에서 이 스펙트럴 보간을 구현했습니다:
  1. DFS 기반 SISL 모델 (요시무라의 안정된 푸리에 기저 사용).
  2. SHT 기반 SISL 모델 (가우스 - 르장드르 격자 사용). 여기서 DFS 방법은 비등간격 위도 지점으로부터 이변수 푸리에 보간자를 생성하도록 적응되었습니다.

주요 기여 및 결과
저자들은 표준 SWE 테스트 케이스 (윌리엄슨 TC1, TC2, TC5, TC6 및 갈레프스키 테스트 케이스) 세트를 사용하여 새로운 스펙트럴 보간 체계를 표준 저차 라그랑주 보간 (LAG) 과 비교 평가했습니다. 결과는 몇 가지 뚜렷한 장점을 보여줍니다:

1. 정확도 및 수렴:

   • TC1 (코사인 종의 대류): 스펙트럴 DFS 보간은 LAG 와 동일한 공간 수렴 속도를 유지하면서 오차를 약 두 자릿수 낮게 달성했습니다. 장기간 적분 시간 (10 회전) 에 걸쳐 DFS 방법은 LAG 대비 오차를 세 자릿수 감소시켰습니다.
   • TC2 (지오스트로픽 흐름): LAG 보간은 (우세한 공간 오차로 인해) 예상되는 2 차 시간 수렴을 보이지 못했지만, DFS 보간 체계는 이론적 2 차 정확도를 달성했습니다.
   • TC5 및 TC6 (산악 흐름 및 로스비 - 하우리츠 파): DFS 방법은 특히 더 높은 해상도와 더 작은 시간 간격에서 LAG 에 비해 높이 필드에서 일관되게 더 낮은 오차를 생성했습니다. TC6 에서 DFS 방법은 파동 골의 형태를 더 잘 보존하여 고차 오일러 참조 해와 더 밀접하게 일치했습니다.

2. 보존 특성:

   • 스펙트럴 보간은 질량 및 에너지 보존을 크게 개선했습니다. TC1 에서 120 일 후 절대 질량 오차는 DFS 의 경우 $1.21 \times 10^{-11}$인 반면 LAG 의 경우 $4.34 \times 10^{-4}$였습니다.
   • TC2 및 TC5 에서 DFS 방법은 시간에 따라 안정된 질량 및 에너지 보존 오차를 유지한 반면, LAG 오차는 증가하거나 더 크게 진동하는 경향이 있었습니다.

3. 수치 확산 및 소규모 역학:

   • 갈레프스키 테스트 케이스: 이 테스트는 비선형 제트의 난류로의 진화를 다룹니다. 운동 에너지 스펙트럼의 스펙트럴 분석 결과, LAG 보간은 고 파수 에너지의 급격한 소산을 유발한 것으로 나타났습니다. 반면, DFS 보간 체계는 에너지 캐스케이드를 보존하여 고 파수 스펙트럼이 난류 회전 흐름에서 관찰되는 이론적 $k^{-5/3}$ 스케일링에 접근하도록 했습니다. LAG 보간은 과도한 수치 감쇠를 나타내며 이 행동을 포착하지 못했습니다.

4. 계산 효율성:

   • 벽시계 시간 비교 결과, 낮은 정확도 요구 사항의 경우 LAG 보간이 더 빠릅니다. 그러나 요구되는 정확도가 증가함에 따라 (더 엄격한 오차 허용 오차), DFS 보간 방법이 더 효율적이 되어 LAG 접근법보다 짧은 시간에 목표 정확도에 도달합니다. 이는 NUFFT 가속 스펙트럴 보간이 고정밀 응용 분야에 실행 가능함을 확인시켜 줍니다.

의의
본 논문은 SISL 체계에 스펙트럴 정확도를 가진 보간을 통합함으로써 공간 이산화와 대류 보간 사이의 수치적 불일치를 해결한다고 주장합니다. 이 연구는 이 접근법이 단순한 이론적 개선이 아니라 다음과 같은 실질적인 이점을 제공함을 보여줍니다:

• 스펙트럴 도함수 계산과 보간 단계 사이의 "정확도 격차"를 제거합니다.
• 수치 확산을 줄여 장기간 적분 기간 동안 미세 구조 및 파동과 같은 특징을 보존합니다.
• 장기 기후 시뮬레이션에 필수적인 질량 및 에너지 보존을 향상시킵니다.
• NUFFT 가속 덕분에 과도한 계산 비용을 치르지 않고 이러한 개선을 달성합니다.

저자들은 저차 보간이 실용적인 절충안이 되어 왔지만, 새로운 스펙트럴 접근법이 고정밀 장기 대기 역학 시뮬레이션을 위한 더 우수한 대안을 제공하며, 이를 더 복잡한 3 차원 역학 코어로 확장하는 동기를 부여한다고 결론지었습니다.