[논문 리뷰] Notes on error estimates for Galerkin approximations of the 'classical' Boussinesq system and related hyperbolic problems

Dimitrios Antonopoulos, Vassilios A. Dougalis|arXiv (Cornell University)|2010. 08. 25.

Advanced Mathematical Physics Problems참고 문헌 21인용 수 10

한 줄 요약

이 논문은 고전적인 일차원 Boussinesq 시스템과 분산성 또는 점성 정규화를 갖는 관련 초월적 문제에 대한 Galerkin-유한요소 근사에 대해 엄밀한 오차 추정을 수립한다. 시간 적분에 명시적 Runge-Kutta 방법을 사용한 준연속 및 완전 연속 스킴을 통해 저자들은 최적 수렴 속도를 증명한다: 균일 메쉬에서는 해 성분에 대해 $L^2$에서 $O(h^2)$이며, 준균일 메쉬에서는 $\eta$에 대해 $O(h^{3/2})$, $u$에 대해 $O(h^2)$이다. 이는 일致한 수렴 차수를 보이는 수치 실험을 통해 검증되었다.

ABSTRACT

We consider the `classical' Boussinesq system in one space dimension and its symmetric analog. These systems model two-way propagation of nonlinear, dispersive long waves of small amplitude on the surface of an ideal fluid in a uniform horizontal channel. We discretize an initial-boundary-value problem for these systems in space using Galerkin-finite element methods and prove error estimates for the resulting semidiscrete problems and also for their fully discrete analogs effected by explicit Runge-Kutta time-stepping procedures. The theoretical orders of convergence obtained are consistent with the results of numerical experiments that are also presented.

연구 동기 및 목표

• 속도에 대한 동차 딜리클레 경계 조건을 갖는 고전적 일차원 Boussinesq 시스템에 대한 준연속 및 완전 연속 Galerkin 근사에 대한 사전 오차 추정을 수립하기 위해.
• 분산성 또는 점성 정규화를 갖는 Boussinesq 시스템 및 관련 초월적 시스템에 대해 Galerkin-유한요소 방법과 명시적 Runge-Kutta 시간 스텝을 조합한 수렴 행동을 분석하기 위해.
• 균일 및 준균일 공간 메쉬 상에서의 수치 실험을 통해 이론적 수렴 속도를 검증하기 위해.
• 기존의 Boussinesq 시스템이 네 개의 경계 조건이 필요로 하는 데 반해, 이 시스템은 오직 두 개의 경계 조건(속도에 대한)만 필요로 하여 잘 정의됨을 보여주기 위해.

제안 방법

• 구간 $[0,1]$에서 연속(piecewise) 선형 함수를 사용한 Galerkin-유한요소 방법을 통한 공간 이산화.
• 안정성 및 수렴 분석을 위해 시간 스텝 크기 $k = h$ 및 $k = Δ x$를 갖는 명시적 Runge-Kutta 스킴을 사용한 시간 이산화.
• 에너지 추정 및 안정성 경계를 도출하기 위해 $L^2$ 기반의 소볼레프 노름과 약한 공식을 사용한 준연속 문제에 대한 분석.
• 공간 근사 오차를 제한하기 위해 $L^2$ 투영 및 보간 오차 추정을 적용.
• 초기-경계값 문제에 대한 약한 해의 존재성과 유일성을 확립하기 위해 포물선 정규화 기법을 적용.
• 제조된 정확한 해를 갖는 완전 연속 스킴을 통한 수치 검증을 통해 이론적 수렴 차수를 확인.

실험 결과

연구 질문

• RQ1 속도에 대한 동질적 딜리클레 경계 조건을 갖는 고전적 Boussinesq 시스템에 대한 Galerkin-유한요소 근사의 최적 수렴 속도는 무엇인가?
• RQ2 분산성 또는 점성 정규화가 존재할 때, 균일 메쉬와 준균일 메쉬 간의 수렴 속도는 어떻게 다를까?
• RQ3 유사한 오차 추정 프레임워크를 $u_{xxt}$ 또는 $u_{xx}$ 항을 갖는 관련 초월적 시스템으로 확장할 수 있는가?
• RQ4 명시적 Runge-Kutta 스킴에 의한 시간 이산화가 완전 연속 스킴의 전체 수렴 차수에 미치는 영향은 무엇인가?
• RQ5 이론적 오차 경계는 실질적 시뮬레이션에서 수치 결과와 어떻게 비교되는가?

주요 결과

• 표면 고도 $\eta$의 $L^2$ 오차는 균일 메쉬에서는 $O(h^2)$로 수렴하고, 준균일 메쉬에서는 $O(h^{3/2})$로 수렴한다.
• 속도 $u$는 균일 및 준균일 메쉬 모두에서 $O(h^2)$로 수렴하며, 이는 이론적 예측과 일致한다.
• 수치 실험은 이론적 수렴 차수를 확인한다: 균일 메쉬에서는 $\eta$에 대해 $O(h^2)$, $u$에 대해 $O(h^2)$이며, 준균일 메쉬에서는 $\eta$에 대해 $O(h^{3/2})$이다.
• 점성 초월 시스템 $\eta_t + u_x = 0$, $u_t + \eta_x - u_{xx} = 0$에 대해, $\eta$의 $L^2$ 오차는 준균일 메쉬에서는 $O(\Delta x)$, 균일 메쉬에서는 $O(h^{3/2})$이며, $u$는 $O(h^2)$로 수렴한다.
• 점성 시스템의 수렴 속도는 선형 분산 케이스와 일致하여 분석 프레임워크의 강건성을 확인한다.
• 수치 결과는 이론적 오차 추정을 강력히 지지하며, 모든 테스트 케이스에서 관측된 수렴 차수와 예측된 속도가 매우 유사하다.