[논문 리뷰] An improvement of the product integration method for a weakly singular Hammerstein equation

Laurence Grammont, Hanane Kaboul|arXiv (Cornell University)|2016. 04. 04.

Numerical methods for differential equations참고 문헌 13인용 수 4

한 줄 요약

이 논문은 약한 특이성을 가진 헤머스타인 적분방정식을 해결하기 위한 새로운 선형화-이산화 접근법을 제안한다. 여기서 뉴턴 방법을 먼저 적용하여 비선형 연산자 방정식을 선형화한 후, 공간 이산화에 제품 통합을 사용한다. 주요 기여는 이산화 파arameter에 관계없이 유도된 반복 근사가 정확한 해로 수렴함을 증명한 것이다. 이는 고전적 이산화-선형화 접근법에서 내재된 정확도 한계를 극복한다.

ABSTRACT

We present a new method to solve nonlinear Hammerstein equations with weakly singular kernels. The process to approximate the solution, followed usually, consists in adapting the discretization scheme from the linear case in order to obtain a nonlinear system in a finite dimensional space and solve it by any linearization method. In this paper, we propose to first linearize, via Newton method, the nonlinear operator equation and only then to discretize the obtained linear equations by the product integration method. We prove that the iterates, issued from our method, tends to the exact solution of the nonlinear Hammerstein equation when the number of Newton iterations tends to infinity, whatever the discretization parameter can be. This is not the case when the discretization is done first: in this case, the accuracy of the approximation is limited by the mesh size discretization. A Numerical example is given to confirm the theorical result.

연구 동기 및 목표

• 이산화를 선형화 이전에 수행할 경우 메쉬 크기로 인해 정확도가 제한되는 고전적 제품 통합 방법의 한계를 해결하기 위해.
• 약한 특이성을 가진 커널을 가진 비선형 헤머스타인 방정식의 수치적 해 정확도를 향상시키기 위해.
• 거의 근사해가 메쉬 크기와 무관하게 정확한 해로 수렴함을 보장하기 위해.
• 선형화 후 이산화를 수행하는 새로운 방법의 이론적 기반을 마련하기 위해.
• 수치 실험을 통해 새로운 방법이 고전적 접근법보다 더 나은 수렴성을 보임을 보여주기 위해.

제안 방법

• 비선형 연산자 방정식 $\varphi - K(\varphi) = y$ 에 직접 뉴턴 방법을 적용하여 선형화된 방정식의 수열을 생성한다.
• 각 선형화된 방정식을 조각별 선형 보간을 사용한 제품 통합 방법으로 이산화한다.
• 특히 약한 특이성을 가진 커널에 적합한 제품 트라프레즈로이드 법칙을 사용하여 수치적 적분을 수행한다.
• 각 뉴턴 반복에서 유도된 선형 시스템을 표준 반복 해법기로 풀어낸다.
• 이산화를 선형화 이후에 수행함으로써 메쉬 크기와 무관하게 수렴성을 높일 수 있도록 보장한다.
• 로그 특이성을 포함한 적분을 정확하게 평가하기 위해 특이성 제거 기법을 활용한다.

실험 결과

연구 질문

• RQ1 선형화를 이산화 이전에 수행할 경우 고전적 접근법보다 더 나은 수렴 성질을 보이는가?
• RQ2 이산화 파arameter가 고정되고 거친 상태일 때도 뉴턴 반복 근사가 정확한 해로 수렴하는가?
• RQ3 약한 특이성을 가진 커널에 대해 제안된 선형화-이산화 방법의 이론적 수렴 성질은 어떠한가?
• RQ4 새로운 방법에서 해의 정확도는 뉴턴 반복 수와 메쉬 크기 사이에 어떻게 영향을 받는가?
• RQ5 반복 수가 증가함에 따라 새로운 방법이 고전적 이산화-선형화 방법보다 오차 감소 측면에서 더 우수한가?

주요 결과

• 제안된 선형화-이산화 방법은 반복 수가 무한히 증가함에 따라 이산화 파arameter에 관계없이 뉴턴 반복 근사가 정확한 해로 수렴함을 보장한다.
• 반대로 고전적 이산화-선형화 방법은 오직 근사해 $\psi_n$ 으로 수렴하며, 이 오차는 메쉬 크기 $n$ 으로 제한되며 진짜 해 $\varphi$ 로 수렴하지 않는다.
• 수치 실험 결과, 새로운 방법은 조밀한 메쉬가 아닐지라도 고전적 방법보다 훨씬 낮은 오차를 달성함을 확인하였다 (예: $n=50$).
• 이산화 파arameter $n$ 이 증가함에 따라 주어진 정밀도에 도달하기 위해 필요한 뉴턴 반복 수가 감소하여 계산 비용을 절감한다.
• 낮은 순서의 제품 통합(예: 트라프레즈로이드 법칙)에서도 효과적으로 작동하므로, 향후 고차수 적분 법칙을 적용해 효율성을 더욱 향상시킬 수 있음을 시사한다.
• 이론적 수렴성은 연산자의 프레셰 미분 가능성과 반복 연산자의 스펙트럼 반경 분석을 통해 확립되었다.