[논문 리뷰] A sparse spectral method for Volterra integral equations using orthogonal polynomials on the triangle

Timon S. Gutleb, Sheehan Olver|arXiv (Cornell University)|2019. 06. 10.

Fractional Differential Equations Solutions참고 문헌 29인용 수 9

한 줄 요약

이 논문은 삼각형 영역 위에서 이변량 수직 다항식을 사용하여 제1종 및 제2종 볼테라 적분방정식을 해결하기 위한 희소 스펙트럼 방법을 제시한다. 가중 카비에 다항식 기저에서 볼테라 연산자의 희소성과 다변량 클렌샤우 알고리즘을 활용함으로써, 비결합 핵을 포함한 경우에도 지수 수렴과 높은 효율성을 달성하며, 토피트스 연산자와의 연결을 통한 엄밀한 수렴 증명이 이루어진다.

ABSTRACT

We introduce and analyse a sparse spectral method for the solution of Volterra integral equations using bivariate orthogonal polynomials on a triangle domain. The sparsity of the Volterra operator on a weighted Jacobi basis is used to achieve high efficiency and exponential convergence. The discussion is followed by a demonstration of the method on example Volterra integral equations of the first and second kind with known analytic solutions as well as an application-oriented numerical experiment. We prove convergence for both first and second kind problems, where the former builds on connections with Toeplitz operators.

연구 동기 및 목표

• 제1종 및 제2종 볼테라 적분방정식을 해결하기 위한 효율적이고 고차수의 수치적 방법을 개발하기 위해.
• 기존의 해법이 결합형 핵에 국한되어 있는 한계를 극복하기 위해.
• 삼각형 영역에서 이변량 자비에 다항식 기저에서 볼테라 연산자의 희소성을 활용하여 계산 효율성을 높이기 위해.
• 연산자 이론적 도구를 사용하여 제1종 및 제2종 문제에 대한 엄밀한 수렴 이론을 수립하기 위해.
• 핵의 다항식 근사 및 변화 이론을 통해 비다항형 핵에의 적용 가능성을 입증하기 위해.

제안 방법

• 이 방법은 삼각형 영역에서 이변량 수직 다항식을 사용하여 핵과 해 함수를 스펙트럼 기저로 표현한다.
• 가중 카비에 다항식 기저에서 볼테라 적분 연산자의 희소성을 활용하여 높은 계산 효율성을 달성한다.
• 볼테라 적분과 유도되는 선형 시스템을 효율적으로 계산하기 위해 클렌샤우 알고리즘의 다변량 변형을 사용한다.
• 이 방법은 적분방정식을 밴드 구조를 가진 선형 시스템으로 변환하여 희소 행렬 기법을 통해 빠른 해를 구현한다.
• 수렴 분석을 위해, 볼테라 연산자를 토피트스 연산자의 컴act한 변화에 연결함으로써, 특히 제1종 방정식에 대해 분석한다.
• 핵의 다항식 근사를 통해 비다항형 핵에 대한 확장이 가능하며, 변화 이론을 활용한다.

실험 결과

연구 질문

• RQ1 삼각형 위에서 이변량 수직 다항식을 기반으로 한 스펙트럼 방법은 제1종 및 제2종 볼테라 적분방정식에 대해 고정밀도와 효율성을 달성할 수 있는가?
• RQ2 카비에 다항식 기저에서 볼테라 연산자의 희소성을 어떻게 활용하여 계산 비용을 줄일 수 있는가?
• RQ3 기존의 많은 스펙트럼 해법이 결합형 핵에 국한되어 있는 것과 달리, 이 방법은 비결합형 핵에 적용 가능한가?
• RQ4 특히 제1종 방정식에 대해 이 방법의 수렴 이론적 기반은 무엇인가?
• RQ5 수렴성을 유지하면서 이 방법을 비다항형 핵으로 확장할 수 있는가?

주요 결과

• 이 방법은 제1종 및 제2종 볼테라 적분방정식 모두에서 지수 수렴을 달성하며, 토피트스 연산자와의 연결 및 유한 섹션 방법을 통한 수렴 증명이 이루어졌다.
• 볼테라 연산자는 삼각형에서 카비에 다항식 기저에서 밴드 구조를 보이며, 이로 인해 유도된 선형 시스템을 효율적으로 해결할 수 있다.
• 기존의 결합형 핵에 국한된 해법과는 달리, 기저와 정의역의 선택 덕분에 더 넓은 핵의 클래스에 적용 가능하다.
• 제1종 방정식의 경우, 수렴은 연산자의 프레드홀름 성질에 의존하며, 이는 K(x,x) ≠ 0 이고 핵이 x=y에서 특이하지 않을 때 성립한다.
• 핵이 점차 증가하는 차수의 다항식으로 근사될 때, 해 시퀀스 uM은 M→∞일 때 진짜 해 u로 수렴하며, 변화 이론에서 유도된 오차 한계가 존재한다.
• 수치 실험을 통해 간단한 문제와 적용 중심의 예제, 포함된 특이 방정식까지 고정밀도와 높은 효율성을 확인하였다.