[논문 리뷰] Initialization and the quasi-geostrophic slow manifold

Stephen M. Cox, A. J. Roberts|ArXiv.org|2003. 03. 07.

Quantum chaos and dynamical systems참고 문헌 31인용 수 4

한 줄 요약

이 논문은 로렌츠의 L86 모델에서 비정상성의 천천히 움직이는 다성분만의 분석을 통해 대기 모델의 초기화를 조사하며, 가장 느린 불변 다성분만(SIM)과 진정한 천천히 움직이는 다성분만이 지수적으로 작은 항을 제외하고 동일하다는 것을 보여준다. 이는 예측의 왜곡을 최소화하는 물리적으로 타당한 초기화 방법을 제안하며, 중력파의 분산은 피할 수 없고, 초기 중력파 활동의 제곱에 반비례하며, 정규형은 파동 분리를 드러내고 초기화 절차에서 내재된 예측의 분산을 보여준다.

ABSTRACT

Atmospheric dynamics span a range of time-scales. The projection of measured data to a slow manifold, ${\cal M}$, removes fast gravity waves from the initial state for numerical simulations of the atmosphere. We explore further the slow manifold for a simple atmospheric model introduced by Lorenz and anticipate that our results will relevant to the vastly more detailed dynamics of atmospheres and oceans. Within the dynamics of the Lorenz model, we make clear the relation between a slow manifold $\cal M$ and the ``slowest invariant manifold'' ({\SIM}), which was constructed by Lorenz in order to avoid the divergence of approximation schemes for $\cal M$. These manifolds are shown to be identical to within exponentially small terms, and so the {\SIM} in fact shares the asymptotic nature of $\cal M$. We also investigate the issue of balancing initial data in order to remove gravity waves. This is a question of how to compute an ``initialized'' point on $\cal M$ whose subsequent evolution matches that from the measured initial data that in general lie off $\cal M$. We propose a choice based on the intuitive idea that the initialization procedure should not significantly alter the forecast. Numerical results demonstrate the utility of our initialization scheme. The normal form for Lorenz' atmospheric model shows clearly how to separate the dynamics of the different atmospheric waves. However, its construction demonstrates that {\em any} initialization procedure must eventually alter the forecast--the time-scale of the divergence between the initialized and the uninitialized solutions is inevitable and is inversely proportional to the square of initial level of gravity-wave activity.

연구 동기 및 목표

• 낮은 차원의 대기 모델에서 천천히 움직이는 다성분만과 가장 느린 불변 다성분만(SIM) 간의 관계를 명확히 하기.
• 측정된 초기 조건을 천천히 움직이는 다성분만에 투영하여 비물리적인 중력파를 억제함으로써 대기 데이터를 초기화하는 데 도전하는 것.
• 초기화 과정에서 피할 수 없는 분산을 고려하면서도, 최소한의 예측 왜곡을 유도하는 물리적으로 타당한 초기화 방법을 개발하는 것.
• L86 모델의 정규형을 분석하여 로시비와 중력파 역학의 분리를 드러내고, 초기화의 본질적 한계를 밝히는 것.

제안 방법

• 비정상성 분석과 정규형 이론을 사용하여 L86 모델의 천천히 움직이는 다성분만을 구성함으로써 빠르고 천천히 움직이는 역학을 분리한다.
• 천천히 움직이는 다성분만과 SIM을 비교하여, 지수적으로 작은 항을 제외하고는 동일하다는 것을 보이며, 이는 이전의 근사 방법에서의 핵심적인 불확실성을 해결한다.
• 특히 자코비 타원 함수를 포함하는 정규형 변환에서 발생하는 적분을 평가하기 위해 윤곽선 적분과 잔여물 계산법을 적용한다.
• 복소 주파수에서의 특이점으로 인해 수렴 반경이 제한되는 발산 급수를 사용하여 천천히 움직이는 다성분만의 점근적 전개를 유도한다.
• 예측 편차를 최소화하는 원리를 기반으로 한 새로운 초기화 방법을 제안하며, 정규형을 사용하여 초기 자료를 천천히 움직이는 다성분만에 투영한다.
• 초기 중력파 진폭의 제곱에 반비례하여 초기화된 해와 초기화되지 않은 해 사이의 분산 시간 스케일이 결정됨을 보여준다.

실험 결과

연구 질문

• RQ1 L86 모델에서 가장 느린 불변 다성분만(SIM)과 진정한 천천히 움직이는 다성분만은 수학적으로 동일한가? 만약 그렇다면, 어떤 정밀도까지 동일한가?
• RQ2 초기 대기 자료를 어떻게 천천히 움직이는 다성분만에 투영할 수 있으며, 이때 예측 왜곡을 최소화할 수 있는가?
• RQ3 비선형 상호작용은 초기화된 해와 초기화되지 않은 해 사이의 분산을 유발하는 데 어떤 역할을 하는가? 그리고 이는 초깃값 중력파 에너지에 어떻게 의존하는가?
• RQ4 정규형과 점근적 전개가 L86 모델에서 로시비파와 중력파 역학의 분리를 어떻게 드러내는가?
• RQ5 초기화 정확도의 기본적인 한계는 무엇이며, 이는 초깃값 파동 활동의 측정으로 어떻게 기술되는가?

주요 결과

• 가장 느린 불변 다성분만(SIM)과 진정한 천천히 움직이는 다성분만은 지수적으로 작은 항을 제외하고 동일하며, 이는 SIM이 천천히 움직이는 다성분만의 타당한 근사임을 확인한다.
• 초기화된 해와 초기화되지 않은 해 사이의 분산 시간 스케일은 초기 중력파 진폭의 제곱에 반비례하며, 이는 초기화 정확도의 기본적인 한계를 나타낸다.
• L86 모델의 정규형은 로시비파와 중력파 역학의 명확한 분리를 드러내며, 중력파 진폭이 로시비파 진폭의 함수로 표현된다.
• 천천히 움직이는 다성분만을 위한 점근적 급수는 복소 주파수에서의 특이점으로 인해 발산하지만, 정규형을 통해 비록 발산하더라도 물리적으로 의미 있는 근사를 가능하게 한다.
• 제안된 초기화 방법은 예측 왜곡을 최소화하며, 수치적으로 효과적임을 보이며, 비선형 파동 상호작용으로 인해 분산이 피할 수 없이 발생한다.
• 정규형 적분의 잔여물 분석은 천천히 움직이는 다성분만에 대한 주요 기여가 복소 평면의 이중극에서 온다는 것을 보여주며, 잔여물은 $ \frac{1}{mW^{*}} e^{-K^\prime / W^{*}} $ 비례로 스케일링된다.