[논문 리뷰] High-order difference and pseudospectral methods for discontinuous problems

C. Markakis, Leor Barack|arXiv (Cornell University)|2014. 06. 18.

Differential Equations and Numerical Methods인용 수 4

한 줄 요약

이 논문은 함수와 그 도함수의 점프를 사전에 알고 있는 경우에 이를 고려해 라그랑주 보간을 보정함으로써 불연속 해를 가진 PDE를 해결하기 위한 고차수 수치 방법을 제시한다. 표준 유한차분 또는 편미분 스펙트럴 행렬을 이러한 점프를 포함하도록 수정함으로써, 불연속성 존재하에서도 스펙트럼 정확도를 확보할 수 있으며, 이는 이동하는 파면이나 충격파를 가진 진화 PDE의 효율적이고 정확한 해를 가능하게 한다.

ABSTRACT

High order finite-difference or spectral methods are typically problematic in approximating a function with a jump discontinuity. Some common remedies come with a cost in accuracy near discontinuities, or in computational cost, or in complexity of implementation. However, for certain classes of problems involving piecewise analytic functions, the jump in the function and its derivatives are known or easy to compute. We show that high-order or spectral accuracy can then be recovered by simply adding to the Lagrange interpolation formula a linear combination of the jumps. Discretizations developed for smooth problems are thus easily extended to nonsmooth problems. Furthermore, in the context of one-dimensional finite-difference or pseudospectral discretizations, numerical integration and differentiation amount to matrix multiplication. We construct the matrices for such operations, in the presence of known discontinuities, by operating on the corrected Lagrange formula. In a method-of-lines framework, this provides a simple and efficient way to obtain solutions with moving discontinuities to evolution partial differential equations.

연구 동기 및 목표

• 점프 불연속성이 있는 문제를 해결할 때 유한차분법 및 편미분 스펙트럴 방법에서 고차수 정확도가 손실되는 문제를 해결하기 위해.
• 기존의 매끄러운 문제에 대한 이산화 방법을 알려진 점프를 가진 비매끄러운 해를 다룰 수 있도록 단순하고 효율적인 방식으로 확장하기 위해.
• 이동 불연속성을 가진 시간에 의존하는 PDE에서 고차수 정확도를 유지하기 위해 방법-라인 프레임워크를 활용하기 위해.
• 점프 보정을 포함한 미분 및 적분 행렬 연산자를 구성함으로써 스펙트럼 수렴성을 유지하기 위해.

제안 방법

• 표준 라그랑주 보간 공식을 함수와 그 도함수의 점프를 선형 조합으로 추가함으로써 수정한다.
• 점프 불연속성으로 인한 보정 항은 점프가 알려져 있거나 쉽게 계산 가능한 것으로 가정하여 해석적으로 유도된다.
• 유한차분 또는 편미분 스펙트럴 행렬은 수정된 보간 공식을 적용하여 미분 및 적분 연산자를 재정의한다.
• 결과적으로 생성된 행렬은 이동 불연속성을 가진 시간에 의존하는 PDE를 해결하기 위해 방법-라인 프레임워크에 사용된다.
• 표준 스펙트럴 방법의 구조를 유지하면서도 불연속점 근처에서 고차수 정확도를 보장한다.
• 점프 정보를 기저 함수에 직접 통합함으로써, 적응 메esh 조정, 도메인 분할, 복잡한 인터페이스 추적 기법이 필요 없도록 한다.

실험 결과

연구 질문

• RQ1 알려진 점프 불연속성을 가진 PDE를 해결할 때 스펙트럴 및 유한차분 방법에서 고차수 정확도를 유지할 수 있는가?
• RQ2 표준 스펙트럴 행렬을 효율성이나 정확도를 훼손하지 않고 불연속성 정보를 통합할 수 있도록 어떻게 수정할 수 있는가?
• RQ3 점프 보정은 편미분 스펙트럴 이산화의 조건수 및 안정성에 어떤 영향을 미치는가?
• RQ4 이 방법은 이동하는 파면을 가진 시간 진화 PDE에 대해 방법-라인 프레임워크에 원활하게 통합될 수 있는가?

주요 결과

• 알려진 점프를 보간 공식에 통합함으로써, 점프 불연속성을 가진 문제에서 스펙트럼 정확도를 복원한다.
• 보정 항은 해석적으로 유도되어 표준 라그랑주 보간에 선형 조합으로 추가되며, 고차수 수렴성을 유지한다.
• 수정된 보간을 사용해 미분 및 적분 행렬 연산자를 구성함으로써, 진화 PDE의 효율적 해법이 가능해진다.
• 이 방법은 이동 불연속성 조건 하에서도 공간과 시간 모두에서 고차수 정확도를 확보할 수 있으며, 복잡한 수치 처리 기법이 필요로 하지 않는다.
• 계산적으로 효율적이며, 표준 스펙트럴 방법의 단순성도 유지하므로, 고비용의 메쉬 재조정이나 적응 메쉬 조정이 필요로 하지 않는다.