[논문 리뷰] On the Rate of Error Growth in Time for Numerical Solutions of Nonlinear Dispersive Wave Equations

Hendrik Ranocha, Manuel Quezada de Luna외 1명|arXiv (Cornell University)|2021. 02. 15.

Advanced Mathematical Physics Problems참고 문헌 83인용 수 13

한 줄 요약

이 논문은 비선형 분산파 방정식의 수치해에서 시간에 따라 변화하는 오차 성장을 조사하며, 광범위한 수치 실험을 통해 보존성 수치적 방법—특히 비선형 불변량을 유지하는 방법—이 선형 오차 성장을 보이며, 비보존성 방법은 이차 오차 성장을 보임을 입증한다. 본 연구는 이와 같은 행동이 비산란성 초순수 체계 및 2차원 얕은 수면 방정식을 포함한 광범위한 방정식 클래스로 확장됨을 보이며, 보존성 방법이 정확성과 장기적 안정성에서 비보존성 방법을 일관되게 능가함을 확인한다.

ABSTRACT

We study the numerical error in solitary wave solutions of nonlinear dispersive wave equations. A number of existing results for discretizations of solitary wave solutions of particular equations indicate that the error grows quadratically in time for numerical methods that do not conserve energy, but grows only linearly for conservative methods. We provide numerical experiments suggesting that this result extends to a very broad class of equations and numerical methods.

연구 동기 및 목표

• 고립파 해에 대해 보존성 수치적 방법에서 관찰된 유리한 선형 오차 성장이 기존의 사례를 넘어서 더 넓은 범위의 비선형 분산파 방정식으로까지 확장되는지 조사하기 위해.
• 보존성 방법이 비선형 또는 비제곱형 불변량을 유지하는 경우에도 여전히 선형 오차 성장을 유지하는지 검토하기 위해.
• 고전적 고립파 해가 이동 대칭성을 가지지 않는 비산란성 1차 히퍼볼릭 체계에서 이와 같은 행동이 지속되는지 탐구하기 위해.
• 변동 수심을 가진 얕은 수면 방정식과 같은 2차원 체계에서 보존성 방법의 성능을 평가하기 위해.
• 보존성 방법이 점 渐진적으로만 아니라 짧은 시간 동안에도 전역 오차를 상당히 줄임을 입증하기 위해, 이론적 프레임워크에서 에너지나 질량 보존이 엄격히 보장되지 않는 경우에도 적용 가능함을 보여주기 위해.

제안 방법

• 시간에 따라 안정성과 보존성을 보장하는 공간 이산화를 위해 합성-통합(SBP) 연산자와 동시에 근사 조건(SAT)을 활용하기 위해.
• 이산 불변량을 보존하는 동시에 명시적 또는 저비용 암시적 시간 스텝을 허용하는 리라크세이션 룬게-쿠타 시간 적분 방법을 사용하기 위해.
• 고차원 공간 정확도를 확보하기 위해 푸리에 스펙트럼 근사 기반의 스펙트럼 콜로케이션 방법을 적용하기 위해.
• Fornberg-Whitham, Camassa-Holm, Degasperis-Procesi, Holm-Hone, BBM-BBM 및 얕은 수면 방정식을 포함한 다양한 방정식에 대해 보존성 및 비보존성 수치적 방법의 변종을 구현하기 위해.
• 보존성 및 비보존성 방법 간의 시간 오차 성장률을 비교하기 위해 장기 시뮬레이션을 수행하기 위해.
• 시간에 따라 수치해와 기준해의 차이의 L2 노름을 계산하여 오차 성장을 평가하기 위해 수치 실험을 실시하기 위해.

실험 결과

연구 질문

• RQ1 기존의 분산 방정식(KdV, NLS 등)에서 보존성 방법에 대해 관찰된 시간에 따른 선형 오차 성장이 더 넓은 범위의 비선형 분산파 방정식으로까지 확장되는가?
• RQ2 선형 오차 성장은 특정 형태의 보존 불변량에 의존하는가, 아니면 비제곱형 불변량에 대해서도 성립하는가?
• RQ3 비산란성 1차 히퍼볼릭 체계에서는 고전적 고립파 해가 존재하지 않지만, 보존성 수치적 방법이 여전히 뛰어난 오차 행동을 달성할 수 있는가?
• RQ4 변동 수심을 가진 얕은 수면 방정식과 같은 2차원 체계에서도 보존성 방법의 이점이 유지되는가?
• RQ5 복잡한 실제 체계인 얕은 수면 방정식과 같은 체계에서 공간 해상도와 수치 오차가 관측된 오차 성장률에 어떤 영향을 미치는가?

주요 결과

• Fornberg-Whitham, Camassa-Holm, Degasperis-Procesi, Holm-Hone 및 BBM-BBM 방정식에 대해 보존성 수치적 방법은 시간에 따라 선형 오차 성장을 보이며, 비보존성 방법은 이차 오차 성장을 보인다.
• 오차 성장 향상 효과는 제곱형 불변량에 국한되지 않으며, 본 연구는 비선형 불변량이 존재하는 경우에도 선형 오차 성장이 유지됨을 확인한다.
• 표준 고립파 대칭성이 없는 1차원 변수 계수 p-체계에서는 보존성 방법이 여전히 상당히 감소된 오차 성장을 보이며, 이는 이 현상이 이동 대칭성 해를 넘어서도 적용된다는 것을 시사한다.
• 변동 수심을 가진 2차원 얕은 수면 시뮬레이션에서, 보존성 방법은 동일한 계산 비용에서 비보존성 방법보다 정확도가 약 한 계급 높은 결과를 도출한다. 이는 정확한 고립파 해가 존재하지 않음에도 불구하고 성립한다.
• 짧은 시간 동안에도 보존성 방법은 비보존성 방법보다 전역 오차가 상당히 작게 나타나, 시간 스케일 전반에서 일관된 우월성을 보여준다.
• 본 연구는 질량 보존이 관측된 선형 오차 성장을 위해 필수적인 조건임을 계산적 증거로 제시하며, 질량 보존을 위반한 경우 이차 오차 성장이 발생함을 보여준다.