[논문 리뷰] The application of the exact operational matrices for solving the Emden-Fowler equations, arising in astrophysics

Kourosh Parand, Hossayni, Sayyed A.외 1명|arXiv (Cornell University)|2013. 10. 07.

Fractional Differential Equations Solutions참고 문헌 72인용 수 8

한 줄 요약

이 논문은 천체물리학에서 발생하는 매우 비선형적인 Emden-Fowler 방정식을 해결하기 위해 베르슈타인 다항식 기반의 정확한 운영 행렬(EOMs)을 도입한다. 기존의 일반 운영 행렬(OOMs)에 비해 훨씬 높은 정확도를 보이며, 정확한 계수 근사 덕분에 거의 정확한 잔차 계산이 가능해져 빠른 오차 수렴을 이룬다. 다양한 비선형성의 테스트 케이스에서 잔차 및 오차 노름 모두에서 OOMs를 능가한다.

ABSTRACT

The objective of this paper is to apply the well-known exact operational matrices (EOMs) idea for solving the Emden-Fowler equations, illustrating the superiority of EOMs versus ordinary operational matrices (OOMs). Up to now, a few studies have been conducted on EOMs and the differential equations solved by them do not have high-degree nonlinearity and the reported results are not regarded as appropriate criteria for the excellence of the new method. So, we chose Emden-Fowler type differential equations and solved them by this method. To confirm the accuracy of the new method and to show the preeminence of EOMs versus OOMs, the norm1 of the residual and error function of both methods are evaluated for multiple $m$ values, where $m$ is the degree of the Bernstein polynomials. We reported the results in form of plots to illustrate the error convergence of both methods to zero and also to show the primacy of the new method versus OOMs. The obtained results have demonstrated the increased accuracy of the new method.

연구 동기 및 목표

• 고차수 비선형 미분방정식, 특히 복잡한 비선형성을 지닌 경우에 기존 일반 운영 행렬(OOMs)의 한계를 해결하기 위해.
• 최근에 소개된 정확한 운영 행렬(EOMs)을 비선형성이 강하고 천체물리학적 의미를 지닌 Emden-Fowler 방정식에 적용하기 위해.
• 이전에 OOMs가 정확한 결과를 도출하지 못한 방정식에 대해 EOMs의 포괄적인 벤치마크를 제공하기 위해.
• 고정밀 수치해(메이플의 dverk78를 활용)를 기반으로 진짜 오차를 계산하고, 다양한 다항식 차수에서 오차 및 잔차 노름을 비교하여 방법을 검증하기 위해.
• 다양한 비선형 항에 대해 노름1 오차 및 잔차 수렴을 정량적으로 비교함으로써 EOMs의 우수성을 입증하기 위해.

제안 방법

• 비선형 항 $ g(y(x)) $ 를 시리즈 운영 행렬을 통해 표현하기 위해 차수 $ m \cdot N $ 의 베르슈타인 다항식을 사용하며, 정확한 계수 근사를 가능하게 한다.
• 베르슈타인 기저 함수에 대한 적분, 미분 및 곱셈 연산에 대해 정확한 운영 행렬(EOMs)을 유도하여 행렬 연산의 수학적 정확성을 보장한다.
• 잔차가 기저 함수와 수직이 되도록 조건을 부과함으로써, 갈레르킨 방법을 적용하여 미분방정식을 대수방정식의 연립으로 변환한다.
• 무한 구간 $[0, \infty)$ 의 Emden-Fowler 방정식을 $[0,1]$ 에 정의된 베르슈타인 다항식과 호환 가능한 $[0, M]$ 으로 변환하기 위해 도메인 절단 및 변수 변환 기법을 사용한다.
• 수렴 특성을 평가하기 위해 증가하는 $ m $ 값에 따라 $ L^1 $-노름을 사용하여 잔차 및 오차 노름을 계산하고 플롯한다.
• 검증을 위해 고차수 연속 룬게-쿠타 방법(메이플의 dverk78)을 기준 해로 사용하여 진짜 오차를 계산한다.

실험 결과

연구 질문

• RQ1 정확한 운영 행렬(EOMs)이 일반 운영 행렬(OOMs)에 비해 매우 비선형적인 Emden-Fowler 방정식을 풀 때 정확도를 크게 향상시킬 수 있는가?
• RQ2 EOMs와 OOMs를 사용할 때 다항식 차수 $ m $ 이 증가함에 따라 오차 및 잔차 노름은 어떻게 수렴하는가?
• RQ3 EOMs는 Emden-Fowler 방정식에서 $ y^{3/2} $, $ y^{5/2} $, $ \sinh(y) $, $ \sin(y) $, $ e^y $ 와 같은 강한 비선형성에 얼마나 잘 대응할 수 있는가?
• RQ4 EOMs를 사용할 경우 잔차 함수가 도메인 전역에서 거의 0에 가까워지며, 이는 높은 해의 정밀도를 의미하는가?
• RQ5 해가 알려져 있지 않은 방정식에 대해 EOMs 접근법이 수치적 기준에 의존하여 신뢰성 있게 적용될 수 있는가?

주요 결과

• 모든 테스트 비선형성에서 EOM 기반 방법은 OOMs보다 훨씬 낮은 $ L^1 $-노름 오차 및 잔차를 달성하였으며, 오차 수렴 속도도 더 빠르게 나타났다.
• 사례 $ f(x) = 1 $, $ g(x) = y^{3/2}(x) $, $ N = 12 $ 에서, 모든 $ m $ 값에서 EOM 방법은 OOMs에 비해 현저히 감소된 오차 및 잔차 노름을 보였다.
• $ f(x) = 1 $, $ g(x) = \sin(y(x)) $, $ N = 12 $ 의 경우, 노름1 오차 및 잔차 플롯은 EOMs가 잔차를 0에 더 가깝게 만들고 오차가 더 빠르게 수렴함을 보여주었다.
• 이 방법은 OOMs가 행렬 연산에서 발생하는 근사 오류로 인해 달성할 수 없었던 거의 정확한 잔차 함수를 달성하였다.
• 고차수 룬게-쿠타 방법을 기준 해로 사용한 결과, EOM 결과는 항상 OOM 결과보다 기준해에 더 가까웠다.
• 시리즈 운영 행렬을 통해 $ g(y(x)) $ 가 고차수의 베르슈타인 다항식의 선형 조합으로 정확하게 표현될 수 있었으며, 이는 강한 비선형성을 다루는 데 핵심적인 역할을 하였다.