[논문 리뷰] On exact solutions and numerics for cold, shallow, and thermocoupled ice sheets

Ed Bueler, Jed Brown|arXiv (Cornell University)|2006. 10. 16.

Cryospheric studies and observations참고 문헌 10인용 수 4

한 줄 요약

이 논문은 냉각된 얕은 열적결합 얼음빙하 모델에 대한 정확한 해석적 해를 제시하여 수치적 방법의 엄밀한 검증을 가능하게 한다. 주어진 두께 및 온도 분포가 결합된 질량 및 에너지 수지 방정식을 만족하도록 하기 위해 필요한 보상적 강수량과 열원 함수를 유도하며, EISMINT 시뮬레이션에서 발생하는 ' spokes'(스포크)의 기원을 온도 이송 기법의 불안정성으로 설명하는 C 언어 구현과 오차 분석을 제공한다.

ABSTRACT

This three section report can be regarded as an extended appendix to (Bueler, Brown, and Lingle 2006). First we give the detailed construction of an exact solution to a standard continuum model of a cold, shallow, and thermocoupled ice sheet. The construction is by calculation of compensatory accumulation and heat source functions which make a chosen pair of functions for thickness and temperature into exact solutions of the coupled system. The solution we construct here is ``TestG'' in (Bueler and others, 2006) and the steady state solution ``Test F'' is a special case. In the second section we give a reference C implementation of these exact solutions. In the last section we give an error analysis of a finite difference scheme for the temperature equation in the thermocoupled model. The error analysis gives three results, first the correct form of the Courant-Friedrichs-Lewy (CFL) condition for stability of the advection scheme, second an equation for error growth which contributes to understanding the famous ``spokes'' of (Payne and others, 2000), and third a convergence theorem under stringent fixed geometry and smoothness assumptions.

연구 동기 및 목표

• 열적결합 얼음빙하 모델에 대해 냉각된 얕은 얼음 근사 조건 하에서 정확한 해석적 해를 구성하기.
• 주어진 두께 및 온도 분포를 가진 알려진 해를 제공함으로써 얼음빙하 모델의 수치적 검증을 엄밀히 가능하게 하기.
• 특히 CFL 조건과 오차 증가에 초점을 맞춰 온도 방정식에 대한 유한차분 스킴의 안정성과 수렴성 분석하기.
• EISMINT II 시뮬레이션에서 나타나는 '스포크'(편향)의 기원을 점탄성 열 발생에 대한 온도에 대한 국소 리프시츠 상수와 연결하여 설명하기.
• 모델 검증 및 테스트에 사용할 수 있도록 정확한 해의 참조 C 구현 제공하기.

제안 방법

• 두께 $ H(t,r) = H_s(r) + \phi(r)\gamma(t) $ 를 시간과 반경에 따라 변하는 형태로 지정하며, $ \phi(r) $ 는余弦형 편미분을 갖는다.
• 온도 프로파일 $ T(t,r,z) = T_s(r) \frac{\nu(t,r) + H(t,r)}{\nu(t,r) + z} $ 을 가정하며, $ \nu $ 는 $ H $ 와 $ T_s $ 로부터 유도되어 경계 조건을 만족하도록 한다.
• 속도장 $ \mathbf{U} $ 는 정수 $ n $ 을 가진 글렌 유동법칙을 사용하여 해석적으로 계산되며, 적분 항등식 $ \int e^x x^n dx = p_n(x) e^x $ 를 활용한다. 여기서 $ p_n $ 은 재귀적으로 정의된 다항식이다.
• 지배 미분방정식에 정확한 필드를 대입하여 필요한 강제항인 보상적 강수량과 열원 $ \Sigma_c $ 를 유도한다.
• 온도 방정식에 대한 유한차분 스킴을 분석하여 올바른 CFL 조건과 국소 절삭 오차를 포함한 오차 증가 부등식을 유도한다.
• C 구현은 격점에서 모든 필드와 강제항의 평가를 포함하여 정확한 해의 참조 코드를 제공한다.

실험 결과

연구 질문

• RQ1 주어진 두께 및 온도 분포가 열적결합 얼음빙하 방정식의 정확한 해가 되도록 하기 위해 필요한 보상적 강수량과 열원 함수의 명시적 형태는 무엇인가?
• RQ2 온도에 대한 점탄성 열 발생의 국소 리프시츠 상수가 수치 스킴에서의 오차 증가에 어떤 영향을 미치는가?
• RQ3 왜 EISMINT II 시뮬레이션에서 '스포크'(편향)가 따뜻한 얼음 영역에서 나타나며, 이는 흐름 법칙의 온도 의존성과 연결될 수 있는가?
• RQ4 이 모델에서 온도 방정식의 이송 스킴에 대한 안정성을 확보하기 위한 올바른 CFL 조건은 무엇인가?
• RQ5 유한차분 스킴의 오차는 시간이 지남에 따라 어떻게 증가하며, 수렴이 보장되는 조건은 무엇인가?

주요 결과

• 정확한 해 'Test G' 는 주어진 $ H(t,r) $ 와 $ T(t,r,z) $ 가 결합된 질량 및 에너지 수지 방정식을 만족하도록 하기 위해 필요한 강제항 $ \Sigma_c $ 와 $ M $ 을 유도함으로써 구성된다.
• 이송 스킴에 대한 올바른 CFL 조건은 $ \Delta t \leq \min_{ijkl} \left\{ \frac{\Delta z}{|w_{ijkl}|} \right\}^{-1} $ 로 유도되며, 이는 수직 속도 크기에 따라 달라진다.
• 스킴의 오차 증가는 시간 간격, 국소 점탄성 열 발생, 절삭 오차를 포함한 국소 부등식에 의해 지배되며, 이 중 첫 번째 항이 오차 증폭에 가장 큰 영향을 미친다.
• 국소 리프시츠 상수 $ L_{ijkl} = \Sigma_{ijkl} \frac{Q}{R T_{ijkl}^2} $ 는 EISMINT II 시뮬레이션에서 '스포크'가 나타나는 이유를 설명한다. 이는 냉각된 얼음에서 따뜻한 얼음으로 전이될 때 급격히 증가하기 때문이다.
• 수렴 정리는 엄격한 가정 하에서 성립한다: 고정된 기하구조와 정확한 해의 매끄러움 조건이 충족될 경우, $ \Delta x \to 0 $ 일 때 유한차분 스킴이 수렴함을 보장한다.
• C 구현은 정확한 해를 위한 완전한 기능성 참조 코드를 제공하여 수치 모델이 알려진 해석 기준에 대해 검증 가능하게 한다.