[논문 리뷰] Error estimates of finite difference schemes for the Korteweg-de Vries equation

Clémentine Courtès, Fré́dé́ric Lagoutière외 1명|arXiv (Cornell University)|2017. 12. 06.

Advanced Mathematical Physics Problems참고 문헌 44인용 수 24

한 줄 요약

이 논문은 Korteweg-de Vries (KdV) 방정식을 푸는 유한차분 스킴에 대한 수렴 속도를 확립한다. 비선형 플럭스에는 명시적 Rusanov 스킴을, 분산 항에는 4점 θ-스킴을 사용한다. θ ≥ 1/2일 때는 초순수 조건 하에서 H^s(ℝ)에서의 일阶 수렴을 증명하며, s ≥ 6 이고, 초순수 조건 하에서 성립한다. θ < 1/2일 때는 '에어리' CFL 조건 하에서 수렴이 보장되며, 더러운 초기 자료에 대해서는 수렴 차수가 낮아진다. 수치 결과는 s ≥ 3일 때 최적의 수렴을 보여준다.

ABSTRACT

This article deals with the numerical analysis of the Cauchy problem for the Korteweg-de Vries equation with a finite difference scheme. We consider the Rusanov scheme for the hyperbolic flux term and a 4-points $\ heta$-scheme for the dispersive term. We prove the convergence under a hyperbolic Courant-Friedrichs-Lewy condition when $\ heta\\geq \\frac{1}{2}$ and under an "Airy" Courant-Friedrichs-Lewy condition when $\ heta&lt;\\frac{1}{2}$. More precisely, we get the first order convergence rate for strong solutions in the Sobolev space $H^s(\\mathbb{R})$, $s \\geq 6$ and extend this result to the non-smooth case for initial data in $H^s(\\mathbb{R})$, with $s\\geq \\frac{3}{4}$ , to the price of a loss in the convergence order. Numerical simulations indicate that the orders of convergence may be optimal when $s\\geq3$.

연구 동기 및 목표

• KdV 방정식에 적용된 유한차분 스킴에 대한 엄밀한 수렴 속도 분석이 부족한 점을 다스린다. 특히 저규칙성 초기 자료에 대해.
• 초기 자료의 규칙성에 따라 수렴 속도를 정의함으로써, 전통적인 부드러운 해를 넘어서는 분석을 수행한다.
• θ-스킴의 시간 스텝 파rameter θ가 안정성과 수렴에 미치는 영향을 분석하며, θ ≥ 1/2와 θ < 1/2를 구분한다.
• 분산 효과를 활용하여 수치 해석과 분산 PDE 이론을 연결함으로써 오차 추정치에서 분산 보정 효과를 활용한다.
• 세차 분산 시스템, 예를 들어 abcd-시스템으로의 확장을 위한 기초를 마련한다.

제안 방법

• 비선형 플럭스 항 ∂x(u²/2)를 명시적 Rusanov 스킴으로 이산화하여 초순수 조건 하에서 안정성을 확보한다.
• 세차 분산 항 ∂x³u를 4점 θ-스킴으로 근사하며, 안정성과 정확도는 θ의 값에 따라 달라진다.
• 합성-부분합 기법과 이산 적분 by parts를 사용하여 이산 에너지 추정치를 유도함으로써 비선형 및 분산 항을 제어한다.
• 고차 이산 도함수를 저차 도함수로 제한하기 위해 이산 Gagliardo-Nirenberg 유형 부등식을 활용한다.
• 안정성 추정치를 연속적인 잘 정의된 문제 결과와 KdV 방정식의 분산 보정 효과를 결합하여 전역 오차를 제어한다.
• 이산 L²-보존 법칙과 이산 ℓ²-노름 추정치를 사용하여 수치 해의 진화를 제어한다.

실험 결과

연구 질문

• RQ1 초기 자료가 H^s(ℝ)에서 제한된 규칙성을 가지며 s ≥ 3/4일 때, KdV 방정식에 대한 유한차분 스킴의 수렴 속도는 무엇인가?
• RQ2 분산 항에 대한 θ-스킴에서 θ의 선택이 스킴의 안정성과 수렴에 미치는 영향은 무엇인가?
• RQ3 θ < 1/2일 때 수렴을 위해 필요한 CFL 조건은 무엇이며, 기존의 초순수 조건과 어떻게 다를까?
• RQ4 저규칙성 초기 자료에 대해 분산 보정 효과를 활용함으로써 수렴 속도를 향상하거나 유지할 수 있는가?
• RQ5 이론적 수렴 속도가 수치 시뮬레이션에서 관찰되는가, 특히 s ≥ 3일 때?

주요 결과

• θ ≥ 1/2일 때는 초순수 조건 하에서 s ≥ 6인 강한 해에 대해 H^s(ℝ) 노름에서 일계 수렴을 달성한다.
• s ≥ 3/4인 H^s(ℝ)의 초기 자료에 대해, 규칙성 손실로 인해 수렴 차수가 낮아지지만 여전히 수렴이 보장된다.
• θ < 1/2일 때는 '에어리' CFL 조건 하에서 수렴이 증명되며, 이는 초순수 조건보다 덜 엄격하지만 여전히 안정성을 보장한다.
• 수치 시뮬레이션 결과는 s ≥ 3일 때 수렴 차수가 최적임을 시사하며, 이는 이론적 한계가 이 영역에서 정확함을 시사한다.
• 저규칙성 초기 자료에 대해 오차를 제어하기 위해 분산 보정 효과가 핵심적으로 기여한다. 이는 제한된 부드러움에도 불구하고 수렴을 가능하게 한다.
• 이 분석은 abcd-시스템과 같은 다른 세차 분산 시스템으로 확장 가능하며, 이는 연속 수준에서 더 강력한 잘 정의된 문제 결과가 필요하다.