[논문 리뷰] On superconvergence of Runge-Kutta convolution quadrature for the wave equation

Jens Markus Melenk, Alexander Rieder|arXiv (Cornell University)|2019. 04. 01.

Electromagnetic Scattering and Analysis참고 문헌 45인용 수 6

한 줄 요약

이 논문은 파동 방정식에 대한 룬게-쿠타 컨volution 퀼드레이처(RK-CQ)에서 초수렴 현상을 설명하기 위해 시간 도함수 형태의 공식을 도입하여 수렴 속도를 향상시킨다. 이를 통해 도입된 새로운 Dirichlet-to-Impedance 맵에 대한 추정식을 활용하여 주파수 의존도를 |s|의 거듭제곱으로 한 단계 감소시키며, 다각형 영역의 경우 로그 항을 포함한다. 이에 따라 고전적 경계를 초월한 향상된 수렴 속도를 달성한다.

ABSTRACT

The semidiscretization of a sound soft scattering problem modelled by the wave equation is analyzed. The spatial treatment is done by integral equation methods. Two temporal discretizations based on Runge-Kutta convolution quadrature are compared: one relies on the incoming wave as input data and the other one is based on its temporal derivative. The convergence rate of the latter is shown to be higher than previously established in the literature. Numerical results indicate sharpness of the analysis. The proof hinges on a novel estimate on the Dirichlet-to-Impedance map for certain Helmholtz problems. Namely, the frequency dependence can be lowered by one power of $\abs{s}$(up to a logarithmic term for polygonal domains) compared to the Dirichlet-to-Neumann map.

연구 동기 및 목표

• 파동 방정식에 대한 RK-CQ에서 관측된 초수렴 현상이 고전적 예측을 초월하는 이유를 설명하는 것.
• 입력 데이터로 들어오는 파동과 그 시간 도함수를 사용하는 두 가지 시간 이산화 방법을 분석하는 것.
• Dirichlet-to-Impedance 맵에 대한 더 날카운 경계를 유도하여 개선된 수렴의 이론적 기반을 마련하는 것.
• 구형 또는 반공간 기하구조에서의 알려진 경계를 일반적인 매끄럽거나 다각형 영역으로 확장하며 볼록성 가정 없이 적용하는 것.
• 고차수 Radau IIA 스킴을 사용한 비볼록(L-shaped) 영역에서의 수치 실험을 통해 이론적 결과를 검증하는 것.

제안 방법

• 경계 적분 방정식을 사용하여 파동 방정식 산란 문제를 공식화하고, 시간 반디스크리티제이션에 컨volution 퀄드레이처를 적용한다.
• 입력 데이터로 들어오는 파동의 시간 도함수를 사용하는 수정된 공식을 도입하여, 수정된 컨볼루션 기호 s⁻¹K(s)를 도출한다.
• Dirichlet-to-Impedance 맵에 대한 새로운 추정식을 유도하여 주파수 의존도를 |s|의 거듭제곱으로 한 단계 감소시키며, 다각형 영역의 경우 로그 인자까지 포함한다.
• 향상된 경계를 활용하여 수정된 방법의 수렴 속도가 기존 공식보다 한 단계 높아짐을 증명한다.
• 반디스크리트 설정에 이론을 적용하여 공간 이산화를 경계 요소법으로 가정하고 시간 이산화에 집중한다.
• 공간 이산화에 갈라킨 경계 요소 방법을 사용하고, 노름 등가성에 대해 연산자 V(1)을 사용하여 H⁻¹/²(Γ) 노름으로 오차를 계산한다.

실험 결과

연구 질문

• RQ1 왜 파동 방정식에 대한 RK-CQ의 수렴 속도가 특정 공식에서는 고전적 예측을 초월하는가?
• RQ2 개선된 경계 적분 연산자 추정을 통해 초수렴 현상이 엄밀하게 설명될 수 있는가?
• RQ3 비볼록 또는 다각형 기하구조에서도 개선된 수렴이 유지되는가? 기하구조는 수렴 속도에 어떤 영향을 미치는가?
• RQ4 들어오는 파동의 시간 도함수를 입력 데이터로 사용할 경우, 기존 공식 대비 수렴 차수는 어떻게 달라지는가?
• RQ5 이론적 추정이 정확한가? 고차수 시간 스텝 스킴에 대해 수치적으로 검증 가능한가?

주요 결과

• 들어오는 파동의 시간 도함수를 사용하는 RK-CQ 방법의 수렴 속도는 기존 공식보다 한 단계 높아지며, 기호 경계에서 순서 q+1−μ+1을 달성한다.
• Dirichlet-to-Impedance 맵의 경우, Dirichlet-to-Neumann 맵 대비 주파수 의존도가 |s|의 거듭제곱으로 한 단계 감소하며, 다각형 영역의 경우 로그 인자를 포함한다.
• 수치 실험을 통해 예측된 수렴 속도가 확인되었다: 3단계 Radau IIA 방법은 순서 3, 5단계 방법은 순서 5이며, 도함수를 사용한 방법은 로그 항을 제외한 고전적 순서에 도달한다.
• 비볼록 L자형 영역에서도 수렴 속도 향상 현상이 관측되어, 이론적 분석이 볼록 기하구조를 초월한 강건성을 보여준다.
• 이론적 분석이 정확함을 입증하며, 수치 결과가 예측된 수렴 차수와 일치함으로써 Dirichlet-to-Impedance 맵에 대한 개선된 경계가 검증된다.
• 초수렴 효과는 s⁻¹-가중 Dirichlet-to-Neumann 맵을 항등원과 Dirichlet-to-Impedance 맵의 합으로 분해함으로써 설명되며, 이는 더 나은 주파수 경계를 허용한다.