[논문 리뷰] Spectral accuracy for the Hahn polynomials

René Goertz, Philipp Öffner|arXiv (Cornell University)|2016. 09. 23.

Mathematical functions and polynomials참고 문헌 17인용 수 4

한 줄 요약

이 논문은 초구형 보존법에 대한 스펙트럴 방법의 맥락에서 하안 다항식의 스펙트럴 정확도를 확립한다. 이는 이산 수직 다항식 투영을 사용함으로써, 부드러운 함수에 대해 하안 다항식 전개가 스펙트럴 수렴을 달성함을 보여주며, 보다 바람직한 노드 분포와 재구성 기법을 통해 룬게 현상 문제를 해결한다.

ABSTRACT

We consider in this paper the Hahn polynomials and their application in numerical methods. The Hahn polynomials are classical discrete orthogonal polynomials. We analyse the behaviour of these polynomials in the context of spectral approximation of partial differential equations. We study series expansions $u=\sum_{n=0}^\infty \hat{u}_n ϕ_n$, where the $ϕ_n$ are the Hahn polynomials. We examine the Hahn coefficients and proof spectral accuracy in some sense. We substantiate our results by numericals tests. Furthermore we discuss a problem which arise by using the Hahn polynomials in the approximation of a function $u$, which is linked to the Runge phenomenon. We suggest two approaches to avoid this problem. These will also be the motivation and the outlook of further research in the application of discrete orthogonal polynomials in a spectral method for the numerical solution of hyperbolic conservation laws.

연구 동기 및 목표

• 부드러운 함수의 근사에서 하안 다항식 전개의 스펙트럴 수렴을 확립하기 위해.
• 등간격 그리드에서의 하안 다항식 보간에 나타나는 룬게 현상 문제를 해결하기 위해.
• 근사 안정성을 향상시키기 위해 대체 노드 분포 및 재구성 전략을 제안하기 위해.
• 유동성 높은 방법(예: 플럭스 재구성(FR) 및 재구성을 통한 수정 절차(CPR))에 이산 수직 다항식을 사용하기 위한 기초를 마련하기 위해.

제안 방법

• 스펙트럴 전개 $ u = \sum_{n=0}^\infty \hat{u}_n \varphi_n $ 에서 하안 다항식을 기저 함수로 사용하며, 초함수 $ Q_n(x;\alpha,\beta,N) = {}_3F_2(-n, n+\alpha+\beta+1, -x; \alpha+1, -N; 1) $ 로 정의된다.
• 연속 수직 다항식 방법에서 흔한 적분 오차를 피하기 위해, 합산을 통한 정확한 계수 계산을 위한 이산 내적을 적용한다.
• 이산 스칼라 tích을 사용한 투영 기반 접근법을 통해 계수 $ \hat{u}_n(t) = \frac{1}{\|\varphi_n\|^2} \sum_{x \in \Omega} h(x) u(x,t) \varphi_n(x) $ 를 계산함으로써 계수 계산의 정확성을 확보한다.
• 바르멘드 시스템 $ u(x_j,t) = \mathcal{V}(\hat{u}_n(t)) $ 을 통한 보간을 조사하며, 조건수 및 노드 선택에 주의를 기울인다.
• 등간격 그리드에서의 룬게 현상을 완화하기 위해 비등간격 노드(예: 체비셰프 유형)를 사용하는 것을 제안한다.
• 지배형 진동을 줄이기 위해 게겐바우어 재구성 기법을 탐색하며, 재구성 순서 $ m $ 과 절단 순서 $ N $ 의 비율을 분석한다.

실험 결과

연구 질문

• RQ1 초구형 보존법의 맥락에서 하안 다항식이 부드러운 함수의 근사에서 스펙트럴 정확도를 달성할 수 있는가?
• RQ2 등간격 그리드에서 하안 다항식 근사에 룬게 현상이 영향을 미치는가? 그리고 이를 완화할 수 있는가?
• RQ3 함수 $ u $ 나 매개변수 $ \alpha, \beta $ 가 어떤 조건을 만족할 경우, 절단된 하안 급수의 안정성과 정확성이 유지되는가?
• RQ4 게겐바우어 유형 재구성에서 재구성 순서 $ m $ 과 절단 순서 $ N $ 의 비율이 총 근사 오차에 어떤 영향을 미치는가?
• RQ5 이산 수직 다항식은 비정렬 메쉬에까지 확장되어 비정렬 고차원 방법에 응용될 수 있는가?

주요 결과

• 해석 함수에 대해 하안 다항식 전개는 스펙트럴 수렴을 달성하며, 계수가 $ n $ 의 어떤 거듭제곱보다도 더 빠르게 감소함을 확인하여 스펙트럴 정확도를 입증한다.
• 수치 실험 결과, $ m=10 $ 인 절단된 하안 급수는 특히 $ \alpha=\beta=0 $ 인 경우에 뚜렷한 점별 오차와 진동을 보이며, 룬게 현상의 존재를 시사한다.
• 매개변수 $ \alpha, \beta $ 의 선택은 근사 품질에 큰 영향을 미치며, 시험 사례에서 $ (5,0) $ 이 $ (0,0) $ 또는 $ (0.5,0.5) $ 보다 더 우수한 수렴을 보였다.
• 등간격 그리드에서의 보간은 조건수가 나쁜 바르멘드 행렬을 초래하고 불안정성을 유도하며, 이는 룬게 현상의 영향을 확인한다.
• 체비셰프 유형 노드 또는 기타 대체 노드 분포를 사용할 경우, 등간격 그리드에 비해 안정성이 향상되고 점별 오차가 감소한다.
• 통제된 $ m/N $ 비율을 가진 게겐바우어 재구성 접근법은 진동을 줄이고 수렴을 향상시킬 수 있으나, $ m $ 과 $ N $ 이 너무 빠르게 증가하면 수렴에 실패할 수 있다.