[논문 리뷰] The Trigonometric Cubic B-spline Algorithm for Burgers' Equation

İdris Dağ, Özlem Ersoy Hepson외 1명|arXiv (Cornell University)|2014. 07. 21.

Advanced Numerical Analysis Techniques참고 문헌 15인용 수 6

한 줄 요약

이 논문은 버거스 방정식을 해결하기 위해 삼각법적 입자 큐빅 B-스플라인 보간법과 크랭크-니콜슨 스킴을 조합한 방법을 제안한다. 이는 급격한 전 fronts 및 충격파 해를 정확하고 안정적으로 포착하는 데 뛰어난 성능을 발휘한다. 수치 실험에서 특히 낮은 점성도의 경우에서 뛰어난 성능을 보이며, 정확한 해와 매우 유사한 결과를 도출하고, 기존의 B-스플라인 기반 방법보다 오차 노름에서 뛰어난 성능을 보였다.

ABSTRACT

The cubic Trigonometric B-spline(CTB) functions are used to set up the collocation method for finding solutions of the Burgers' equation. The effect of the CTB in the collocation method is sought by studying two text problems. The Burgers' equation is fully-discretized using the Crank-Nicholson method for the time discretizion and CTB function for discretizion of spatial variable. Numerical examples are performed to show the convenience of the method for solutions of Burgers equation

연구 동기 및 목표

• 낮은 점성도로 인해 급격한 전 fronts 및 충격파를 형성하는 경우를 포함한 버거스 방정식을 해결하기 위한 강력한 수치적 스킴을 개발하는 것.
• 비선형 편미분방정식에 대한 보간법에서 삼각법적 입자 큐빅 B-스플라인(CTB)의 효과성을 조사하는 것.
• 기존의 B-스플라인 및 갈레르킨 방법과 비교하여 제안된 CTB 보간법의 정확도 및 수렴 특성에 대해 평가하는 것.
• 시간에 따라 해의 부드러움과 안정성을 유지할 수 있는 능력을 입증하는 것.
• 삼각함수 기반 기저 함수를 사용하여 전통적인 B-스플라인 보간법의 신뢰할 수 있는 대안을 제공하는 것.

제안 방법

• 공간 영역은 두 번 연속 미분 가능한 삼각법적 입자 큐빅 B-스플라인(CTB) 기저 함수를 사용하여 이산화하며, 이는 고차 정확도를 제공한다.
• 시간 이산화에는 크랭크-니콜슨 방법을 적용하여 2차 정밀도와 무조건적 안정성을 확보한다.
• 보간법에서는 해를 CTB 함수의 선형 조합으로 근사하고, 잔차를 선택된 보간점에서 0이 되도록 강제한다.
• 유도된 비선형 방정식 시스템은 반복적으로 해결되며, 각 시간 단계에서 경계 조건과 초기 조건이 적용된다.
• CTB 보간 행렬의 대각선 구조를 활용하여 각 시간 단계에서 선형 시스템을 효율적으로 해결한다.
• CTB 기저는 형상 파라미터를 가진 삼각함수로 정의되며, 진동하거나 급격한 해를 갖는 경우의 부드러움과 개선된 근사 성능을 보장한다.

실험 결과

연구 질문

• RQ1 충격파가 존재하는 상황에서 표준 큐빅 B-스플라인에 비해 삼각법적 입자 큐빅 B-스플라인(CTB)이 버거스 방정식의 해를 더 정확하고 안정적으로 제공할 수 있는가?
• RQ2 점성도 파라미터 λ가 작을 때 CTB 보간법이 급격한 전 fronts 해를 어떻게 포착하는가?
• RQ3 다양한 공간 및 시간 이산화에 대해 CTB 보간법의 수렴 특성과 오차 분포는 어떠한가?
• RQ4 기존의 B-스플라인 보간법 및 갈레르킨 방법과 비교하여 제안된 방법의 정확도 및 계산 효율성은 어떠한가?
• RQ5 삼각함수 기반 기저 함수의 사용이 비선형 편미분방정식에서 이동파 해의 근사에 기여하는가?

주요 결과

• 표본 5에 나타낸 바와 같이, λ = 0.01 및 h = 1/36일 때 최대 절대 오차가 약 10^-3 수준으로 매우 높은 정확도를 확보하였다.
• t = 0.5일 때, 수치적 해는 정확한 해와 매우 유사하며, 대부분의 공간 점에서 상대 오차가 0.002 이하로 유지되었다.
• λ = 0.001일 경우, 방법은 오른쪽으로 진행되는 날카로운 충격파 전 fronts 형성과 전파를 성공적으로 포착하였으며, 해석적 기대와 일치하였다.
• L∞ 오차 노름은 유한하고 공간 해상도가 높아질수록 감소하여 수렴성과 안정성을 보였다.
• 특히 급격한 기울기를 가지는 영역에서 기존의 큐빅 B-스플라인 보간법 및 갈레르킨 방법보다 정확도에서 뛰어나거나 동등한 성능을 보였다.
• 그림 5~8의 시각적 결과는 방법이 시간에 따라 해의 부드러움을 유지하고 충격 프로파일을 정확히 해석할 수 있음을 확인하였다.