[논문 리뷰] Hamiltonian interpolation of splitting approximations for nonlinear PDEs

Erwan Faou, Benoît Grébert|arXiv (Cornell University)|2009. 12. 15.

Numerical methods for differential equations참고 문헌 3인용 수 4

한 줄 요약

이 논문은 다항 비선형성을 가진 비선형 해밀턴 PDE에 적용된 분할 적분기법에 대해 해밀턴 보간을 수립하며, 수정된 해밀턴 함수의 수치적 흐름이 지수적으로 긴 시간 동안 정확한 흐름과 밀접하게 일치함을 보여준다. 스펙트럼 필터링을 고려한 후행 오차 분석을 통해 수정된 해밀턴 함수를 구성함으로써, 공진 조건이 없고 적분 가능성이 있는 조건 하에서 에너지와 정규성의 장기적 유지가 보장된다.

ABSTRACT

We consider a wide class of semi linear Hamiltonian partial differential equa- tions and their approximation by time splitting methods. We assume that the nonlinearity is polynomial, and that the numerical tra jectory remains at least uni- formly integrable with respect to an eigenbasis of the linear operator (typically the Fourier basis). We show the existence of a modified interpolated Hamiltonian equation whose exact solution coincides with the discrete flow at each time step over a long time depending on a non resonance condition satisfied by the stepsize. We introduce a class of modified splitting schemes fulfilling this condition at a high order and prove for them that the numerical flow and the continuous flow remain close over exponentially long time with respect to the step size. For stan- dard splitting or implicit-explicit scheme, such a backward error analysis result holds true on a time depending on a cut-off condition in the high frequencies (CFL condition). This analysis is valid in the case where the linear operator has a discrete (bounded domain) or continuous (the whole space) spectrum.

연구 동기 및 목표

• 유계가 아닌 선형 연산자를 가진 반선형 해밀턴 PDE로 후행 오차 분석을 확장하여 유한차원 추정에서의 상수의 폭주 문제를 해결한다.
• 수치적 흐름이 지수적으로 긴 시간 동안 각 시간 단계에서 수정된 해밀턴 함수의 정확한 흐름과 일치하는 수정된 해밀턴 함수의 존재를 확립한다.
• 스펙트럼 필터링을 적용한 분할 스킴(예: 리-트로터 스킴)을 비선형 샤크호프 및 웨이브 방정식에 적용하여 장기적 안정성을 확보한다.
• 푸리에 공간에서의 공진 조건이 없고 균일한 적분 가능성이 있는 조건 하에서 수치적 해가 연속 흐름에 가까워지는 것을 분석한다.
• 스텝 크기에 따라 지수적으로 긴 시간 동안 에너지 보존 및 정규성에 대한 정량적 경계를 유도한다.

제안 방법

• 형식적 후행 오차 분석을 통해 해밀턴 함수의 선형 및 비선형 부분 간의 반복된 푸아송 괄호를 사용하여 $ h $의 멱급수 형태로 수정된 해밀턴 함수 $ H_h $를 구성한다.
• 선형 연산자의 유계성을 확보하고 무한대가 되는 고유값을 피하기 위해 스펙트럼 필터 함수 $ \alpha_h(x) $ 를 도입한다.
• 이 방법은 토러스(이산 스펙트럼)와 전체 공간(연속 스펙트럼) 모두에 적용 가능하며, 푸리에 기반 표현 및 가짜 스펙트럼 방법을 사용한다.
• 핵심적인 추정은 고주파 성분의 성장을 제어하기 위해 푸리에 공간에서의 영동역 조건과 해의 $ L^1 $-유계성을 기반으로 한다.
• 수정된 해밀턴 함수 방정식의 재귀적 해법은 빈도 구성에 대한 다중색인 적분을 사용하는 리-푸아송 괄호 체계를 활용한다.
• 수렴성과 오차 항의 지수적 소실을 보장하기 위해 공진 조건 $ |\Lambda(\boldsymbol{\tau})| \geq 2\pi - \delta $ 이 도입된다.

실험 결과

연구 질문

• RQ1 정확한 흐름이 각 시간 단계에서 수정된 해밀턴 함수 $ H_h $의 수치적 흐름과 일치하는 수정된 해밀턴 함수를 구성할 수 있는가?
• RQ2 어떤 조건 하에서 분할 스킴의 수치적 해가 지수적으로 긴 시간 동안 연속 흐름에 가까워지는가?
• RQ3 스펙트럼 필터링을 어떻게 활용하여 선형 연산자가 무한대인 PDE의 분할 방법을 안정화할 수 있는가?
• RQ4 푸리에 공간에서의 균일한 적분 가능성은 장기적 정확성을 보장하는 데 어떤 역할을 하는가?
• RQ5 이 수정된 해밀턴 함수 프레임워크를 통해 지수적으로 긴 시간 동안 에너지와 정규성을 어느 정도까지 유지할 수 있는가?

주요 결과

• 분할 스킴의 수치적 흐름은 $ nh \leq \exp(c h^{-1/2}) $ 인 경우, 지수적으로 작은 오차 $ \mathcal{O}(\exp(-c h^{-1/2})) $ 까지 수정된 해밀턴 함수 $ H_h $ 의 정확한 흐름과 일치한다.
• 비선형성의 세제곱형 NLS 방정식에서 $ \alpha_h(x) = \sqrt{h} \arctan(x/\sqrt{h}) $ 를 사용할 경우, 수정된 해밀턴 함수 $ H_h^{(1)} $ 는 지수적으로 긴 시간 동안 초기 에너지를 $ \mathcal{O}(h) $ 수준으로 유지한다.
• 초기 자료가 $ H^1(\mathbb{R}^3) $ 에 속해 있을 경우, 고주파 성분이 필터를 통해 제어됨에 따라 $ H^1 $-노름에서 해가 지수적으로 긴 시간 동안 바운드되어 있다.
• 해가 $ B^s_M $ 에 속해 있고 $ s \geq 3 $ 라면, 원래 해밀턴 함수 $ H $ 는 지수적으로 긴 시간 동안 $ \mathcal{O}(h) $ 수준으로 유지된다.
• 공진 조건 $ |\Lambda(\boldsymbol{\tau})| \leq 2\pi - \delta $ 는 형식적 급수의 수렴성과 수정된 시스템의 지수적 안정성을 보장한다.
• 이 방법은 적절한 스펙트럼 필터링 하에 이산 및 연속 스펙트럼 모두에 대해 동일하게 적용 가능하며, 전체 공간 $ \mathbb{R}^d $ 에도 적용 가능하다.