[논문 리뷰] A compact higher-order finite-difference scheme for the wave equation can be strongly non-dissipative on non-uniform meshes

Alexander Zlotnik, Raimondas Čiegis|arXiv (Cornell University)|2020. 01. 01.

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

한 줄 요약

이 논문은 1D 파동방정식에 대한 밀도 높은 고차수 유한차분 스킴이 비균일 메시에서 복소 고유값이 발생할 경우 강력한 비소산성 및 지수적 불안정성을 띠게 되며, 이는 해의 노름이 시간에 따라 지수적으로 증가하기 때문에, 균일 메시에서의 고차 정확도에도 불구하고 실용적으로는 사용할 수 없음을 보여준다.

ABSTRACT

We study necessary conditions for stability of a Numerov-type compact higher-order finite-difference scheme for the 1D homogeneous wave equation in the case of non-uniform spatial meshes. We first show that the uniform in time stability cannot be valid in any spatial norm provided that the complex eigenvalues appear in the associated mesh eigenvalue problem. Moreover, we prove that then the solution norm grows exponentially in time making the scheme strongly non-dissipative and therefore impractical. Numerical results confirm this conclusion. In addition, for some sequences of refining spatial meshes, an excessively strong condition between steps in time and space is necessary (even for the non-uniform in time stability) which is familiar for explicit schemes in the parabolic case.

연구 동기 및 목표

• 비균일 공간 메시에서 1D 파동방정식에 대한 Numerov 유형의 밀도 높은 고차수 유한차분 스킴의 안정성 분석.
• 스킴이 안정하거나 불안정해지는 조건을 규명하며, 특히 복소 고유값이 나타나는 조건에 초점을 맞춘다.
• 복소 고유값이 나타날 경우, 시간에 따라 일관된 안정성이 불가능하며, 해의 노름이 지수적으로 증가함을 입증한다.
• 비균일 메시에서의 안정성을 확보하기 위해 철저한 시간단계 제약, 예를 들어 $\tau \leqslant c_0 h_{\min}^2$ 가 필요한지 조사한다.
• 세밀화된 비균일 메시에서의 수치 실험을 통해 이론적 결과를 검증한다.

제안 방법

• 비균일 공간 메시를 가진 1D 파동방정식에 대해 Numerov 방법을 기반으로 한 삼단계 밀도 높은 유한차분 스킴을 수립한다.
• 메시 가중 내적을 사용하여 이산 에너지 노름을 정의하고, 차분 연산자의 관련 고유값 문제를 통해 안정성 분석을 수행한다.
• 이산 고유값 문제에서 복소 고유값이 나타나면, 시간에 따라 해의 노름이 지수적으로 증가함을 증명한다.
• 안정성에 필요한 조건을 유도하여, 복소 고유값이 존재할 경우 어떤 공간 노름에서도 시간에 따라 일관된 안정성이 성립하지 않음을 보여준다.
• 메시 세분화를 연구하기 위해 가족형 세밀화된 비균일 메시 $\overline{\omega}_{h}^{0,K}$ 를 사용하고, $\tau$ 와 $h_{\min}$ 을 포함한 점점 더 세밀해지는 안정성 조건을 유도한다.
• 공간적으로 국소화된 초기 자료를 가진 수치 실험을 수행하여 시간에 따른 오차 증가를 계산하고, 이론적 예측된 지수적 불안정성을 확인한다.

실험 결과

연구 질문

• RQ1 밀도 높은 고차수 유한차분 스킴이 비균일 메시에서 파동방정식에 대해 언제 불안정해지는가?
• RQ2 이산 고유값 문제에서 복소 고유값이 나타날 경우, 시간에 따라 일관된 안정성을 달성할 수 있는가?
• RQ3 비균일 메시에서 비균일 시간 안정성을 확보하기 위해 시간단계 $\tau$ 와 최소 공간 단계 $h_{\min}$ 사이에 필요한 관계는 무엇인가?
• RQ4 오차의 지수적 증가율은 메시 세분화와 시간 지속 기간에 따라 어떻게 달라지는가?
• RQ5 수치 결과가 이론적 불안정성 예측을 어느 정도 확인하는가?

주요 결과

• 이산 고유값 문제에서 복소 고유값이 존재하면 해의 노름이 지수적으로 증가하여, 스킴이 강력한 비소산성 상태가 되며, 결과적으로 실용적으로는 사용할 수 없게 된다.
• 스펙트럼 분석을 통해 복소 고유값이 존재할 경우 어떤 공간 노름에서도 시간에 따라 일관된 안정성이 불가능함을 입증하였다.
• 일부 세밀화된 비균일 메시의 수열에 대해서는, 비균일 시간 안정성조차도 $\tau \leqslant c_0 h_{\min}^2$ 의 안정성 조건이 필수적임을 보여주었으며, 이는 포물형 방정식의 명시적 스킴과 유사하다.
• 수치 실험 결과에 따르면 치명적인 오차 증가가 발생한다: $T=2$ 에서 $K=20$ 에서 오차는 $1.443 \times 10^{-4}$ 에서 $K=80$ 에서는 $1.082 \times 10^5$ 로 증가하여 지수적 불안정성을 확인한다.
• $K=20$ 에서 관측된 실용적 증가율 $\varkappa_{pr} \approx 4.817$ 이 이론적 추정치 $\varkappa_0 \approx 4.579$ 와 매우 유사하여 이론 모델의 타당성을 검증한다.
• 스킴 내의 $\sigma$ 값이 무엇이든 복소 고유값이 존재하는 한 스킴은 항상 불안정하며, 핵심 안정성 조건에서는 $\sigma$ 와 무관하다.