[논문 리뷰] Uniformly accurate numerical schemes for highly oscillatory Klein-Gordon and nonlinear Schrödinger equations

Philippe Chartier, Nicolas Crouseilles외 2명|arXiv (Cornell University)|2013. 08. 02.

Electromagnetic Simulation and Numerical Methods참고 문헌 14인용 수 4

한 줄 요약

이 논문은 높은 주파수를 갖는 콘드라지와 비선형 슈뢰딩거 방정식을 위한 균일 정확도(Uniformly Accurate, UA) 수치적 방법을 제안한다. 이는 빠른 시간 동역학과 느린 시간 동역학을 분리하기 위해 이중 척도 형식을 도입함으로써 달성된다. 이 방법은 메쉬나 시간 간격 조정 없이도 모든 영역—비진동, 중간, 고진동—에서 균일 정확도를 유지하며, 점근 모델에 의존하지 않는다. 또한 $\varepsilon$의 크기와 관계없이 정확도의 순서가 유지된다. 주요 기여는 모든 $\varepsilon \in (0,\varepsilon_0]$에 대해 안정적이고 균일 수렴하는 수치적 방법을 확보한 것이다. 이는 이중 척도 프레임워크 내에서 챔프먼-엔스코그 유사 전개를 적용함으로써 달성되며, $\varepsilon$가 작지 않은 경우에도 안정성과 정확도를 보장한다.

ABSTRACT

This work is devoted to the numerical simulation of nonlinear Schrödinger and Klein-Gordon equations. We present a general strategy to construct numerical schemes which are uniformly accurate with respect to the oscillation frequency. This is a stronger feature than the usual so called "Asymptotic preserving" property, the last being also satisfied by our scheme in the highly oscillatory limit. Our strategy enables to simulate the oscillatory problem without using any mesh or time step refinement, and the orders of our schemes are preserved uniformly in all regimes. In other words, since our numerical method is not based on the derivation and the simulation of asymptotic models, it works in the regime where the solution does not oscillate rapidly, in the highly oscillatory limit regime, and in the intermediate regime with the same order of accuracy. In the same spirit as in [5], the method is based on two main ingredients. First, we embed our problem in a suitable "two-scal" reformulation with the introduction of an additional variable. Then a link is made with classical strategies based on Chapman-Enskog expansions in kinetic theory despite the dispersive context of the targeted equations, allowing to separate the fast time scale from the slow one. Uniformly accurate (UA) schemes are eventually derived from this new formulation and their properties and performances are assessed both theoretically and numerically.

연구 동기 및 목표

• 모든 진동 주파수 $\varepsilon$—비진동, 중간, 고진동 영역—에서 균일 정확도를 유지하는 수치적 방법을 개발한다.
• 표준 수치 방법이 정확도를 확보하기 위해 $\Delta t \lesssim \varepsilon^{q/p}$ 가 필요로 하는 한계를 극복한다. 이는 $\varepsilon \to 0$일 때 실용적이지 않게 된다.
• 점근 모델이나 그 절단 오차에 의존하지 않고도 균일 정확도(Uniformly Accurate, UA)를 확보하는 방법을 구축한다.
• 정확도의 순서가 $\varepsilon$에 관계없이 균일하게 유지되도록 보장한다. 특히 $\varepsilon$가 작지 않은 경우에도 마찬가지다.
• 점근 보존(Asymptotically Preserving, AP)이면서 동시에 모든 영역—비진동, 중간, 고진동—에서 균일 수렴성을 유지하는 프레임워크를 제공한다. 특히 AP 방법이 실패할 수 있는 중간 영역에서도 성립한다.

제안 방법

• 보조 변수 $\tau = t/\varepsilon$ 를 도입하여 원래 방정식을 이중 척도 형식으로 통합함으로써 $U^\varepsilon(t,\tau)$ 를 정의하고, $u^\varepsilon(t) = U^\varepsilon(t, t/\varepsilon)$ 를 만족하도록 한다.
• 원래 방정식을 $(t,\tau)$ 에 대한 비선형 편미분방정식으로 재구성함: $\partial_t U^\varepsilon + \frac{1}{\varepsilon} \partial_\tau U^\varepsilon = \mathcal{F}(t,\tau, U^\varepsilon)$ 를 통해 빠른 동역학과 느린 동역학을 분리한다.
• 이중 척도 체계에 형식적인 챔프먼-엔스코그 전개를 적용하여 해를 느린 성분과 빠른 성분으로 분해하고, $\varepsilon$ 의 거듭제곱에 따른 점근 전개를 가능하게 한다.
• 이중 척도의 구조를 유지하면서 시간 분할 또는 룬게-쿠타 방법을 활용해 전개된 체계를 근사함으로써 균일 정확도(UA) 수치적 방법을 구성한다. 이는 $\varepsilon$ 에 대해 균일하게 정확도 순서를 유지한다.
• 특히 고주파 성분이 존재할 경우에도 전역 오차의 균일 수렴성을 확보하기 위해, 이중 척도 시스템의 초기 조건을 2차 또는 3차까지 보정한다.
• 비선형 콘드라지 방정식과 비선형 슈뢰딩거 방정식에 대해 수치적 검증을 수행하고, 다양한 $\varepsilon$ 와 $\Delta t$ 에서 $H^1$ 노름 기반 오차를 비교한다.

실험 결과

연구 질문

• RQ1 모든 진동 매개변수 $\varepsilon$ 값—비진동, 중간, 고진동 영역—에서 정확도가 균일하게 유지되는 수치적 방법을 설계할 수 있는가? 특히 $\varepsilon$ 가 작지 않은 경우에도 마찬가지인가?
• RQ2 점근 모델을 유도하지 않고도 이중 척도 형식을 활용해 고진동 편미분방정식의 빠른 동역학과 느린 동역학을 분리할 수 있는가?
• RQ3 이중 척도 시스템에서의 초기 조건 보정이 전역 오차의 균일성에 미치는 영향은 무엇인가?
• RQ4 제안된 방법이 중간 영역에서 $\varepsilon$ 가 크지도 작지도 않은 경우에도 정확도의 순서가 $\varepsilon$ 에 대해 균일하게 유지되는가?
• RQ5 운동량 보존 구조가 없는 분산 방정식에서, 이중 척도 프레임워크 내에서 챔프먼-엔스코그 유사 전개를 활용해 균일 정확도 수치적 방법을 유도할 수 있는가?

주요 결과

• 제안된 균일 정확도(UA) 수치적 방법은 모든 $\varepsilon \in (0,\varepsilon_0]$ 영역—비진동($\varepsilon \sim 1$), 중간, 고진동($\varepsilon \to 0$) 영역—에서 정확도의 순서가 균일하게 유지된다.
• 수치 실험 결과, 이중 척도 시스템의 초기 조건을 2차 또는 3차로 보정한 경우 $H^1$ 상대 오차가 $\varepsilon$ 에 대해 균일하게 유지되지만, 보정이 없을 경우 오차가 심하게 악화됨을 확인하였다.
• 2차 UA 수치적 방법에서 초기 조건을 3차로 보정한 경우, 오차는 $\varepsilon$ 에 관계없이 균일하게 유계이며, 모든 $\varepsilon$ 값에서 $\mathcal{O}(\Delta t^2)$ 의 스케일링을 보인다.
• 초기 조건을 보정하지 않은 경우 오차는 $\varepsilon$ 에 대해 균일하지 않게 되며, 특히 $\varepsilon$ 가 작은 경우에 심해짐을 확인하여, 초기 조건 보정이 균일 수렴성 확보에 필수적임을 입증하였다.
• 공간에 대한 푸리에 모드의 시간 진동을 분석한 결과, 고차수 초기 조건 보정이 $U(t,\tau=0,x)$ 에서 고주파 성분을 억제함을 확인하였으며, 이는 이론적 기대와 일치한다. 즉, $k \geq 1$ 인 경우 $\mathcal{O}(\varepsilon^k)$ 정도의 성분만 진동성을 유지한다.
• UA 수치적 방법으로 계산한 해는 스트랭 분할로 계산한 기준 해와 스토로보스코픽 시점 $t_k = 2\pi k \varepsilon$ 에서 일치함을 확인하여, 이중 척도 형식의 일관성을 입증하였다.