[논문 리뷰] Galerkin Method for the numerical solution of the RLW equation by using exponential B-splines

Melis Zorşahin Görgülü, İdris Dağ외 1명|arXiv (Cornell University)|2015. 04. 21.

Fractional Differential Equations Solutions참고 문헌 33인용 수 4

한 줄 요약

이 논문은 얕은 수면파 및 플라즈마 파동을 모델링하는 비선형 편미분방정식인 정규화된 장파(RLW) 방정식을 해결하기 위해 지수 B-스플라인을 사용하는 갈레르킨 유한요소법을 제안한다. 이 방법은 높은 정확도와 불변량의 보존성을 달성하여, 콜로케이션 방법보다 우수하고 기존 갈레르킨 접근법과 비교해도 성능이 유사하거나 향상됨을 입증하였다. 단일 파동 전파, 상호작용, 그리고 경계력에 의한 파동 생성과 같은 벤치마크 테스트에서 성능을 평가하였다.

ABSTRACT

In this paper, the exponential B-spline functions are used for the numerical solution of the RLW equation. Three numerical examples related to propagation of single solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.

연구 동기 및 목표

• 비선형 분산파를 모델링하는 RLW 방정식을 해결하기 위한 안정적인 수치적 방법을 지수 B-스플라인 기반으로 개발한다.
• 기존의 스플라인 기반 방법을 향상시키기 위해 다항식 스플라인보다 급격한 기울기와 진동 행동을 더 효과적으로 다룰 수 있는 지수 B-스플라인의 능력을 활용한다.
• 단일 파동 전파, 두 파동의 상호작용, 경계력에 의한 파동 생성이라는 세 가지 벤치마크 문제를 통해 방법의 정확도와 보존 성질을 검증한다.
• 기존의 수치 기법과의 비교를 통해 오차 노름과 불변량의 보존 측면에서 우수하거나 동등한 성능을 보임을 입증한다.

제안 방법

• 지수 B-스플라인 함수를 갈레르킨 유한요소법의 시험 함수로 사용하여 RLW 방정식의 해를 근사한다.
• RLW 방정식의 약한 형태를 적용하여 비선형 대수방정식계를 도출하고, 반복적으로 이를 해석한다.
• 형태 매개변수 p를 사용하여 hyperbolic 함수로 정의된 지수 B-스플라인은 급격한 변화나 감쇠를 보이는 해의 보다 우수한 근사가 가능하다.
• 시간 적분은 크랭크-니콜슨 스킴을 사용하여 시간에 대해 두 번째 차수 정확도와 무조건적 안정성을 확보한다.
• 영역 외부의 연장된 키노트를 사용하여 경계 조건을 구현하고, 경계에서 동차 뉴먼 조건을 적용한다.
• 장기적인 안정성과 정확도 평가를 위해 시뮬레이션 전반에 걸쳐 $ C_1 $, $ C_2 $, $ C_3 $ 등의 보존량을 모니터링한다.

실험 결과

연구 질문

• RQ1 기존의 다항식 B-스플라인 또는 콜로케이션 방법과 비교해 지수 B-스플라인이 RLW 방정식의 더 정확하고 안정적인 수치적 해를 제공할 수 있는가?
• RQ2 장시간 시뮬레이션 동안 지수 B-스플라인 갈레르킨 방법이 $ C_1 $, $ C_2 $, $ C_3 $ 불변량을 얼마나 잘 보존하는가?
• RQ3 이 방법은 단일 파동 전파, 두 파동 상호작용, 경계력에 의한 파동 생성 시뮬레이션에서 정확도가 어느 정도인가?
• RQ4 특히 파동 프로파일과 진폭을 정확히 포착할 수 있도록, 굵은 공간 및 시간 메esh에서도 높은 정확도를 유지하는가?
• RQ5 기존 갈레르킨 및 콜로케이션 방법과 비교해 이 방법의 성능과 오차 노름 측면에서 어떤가?

주요 결과

• 단일 파동 전파 테스트에서 지수 B-스플라인 갈레르킨 방법(EBSGM)은 $ L_\infty $ 오차 노름이 $ 1.90 \times 10^{-3}$ 이하로 유지되어 높은 정확도를 보였다.
• 두 파동 상호작용 테스트에서는 $ C_1 $, $ C_2 $, $ C_3 $ 불변량이 시간이 지남에 따라 상대 오차 0.1~0.3퍼센트 이내로 유지되어 뛰어난 보존 성질을 입증하였다.
• 파동 생성 테스트에서는 진폭 $ U_0 = 2 $ 인 단일 파동의 형성이 성공적으로 시뮬레이션되었고, $ t = 100 $ 에서 진폭이 약 3.76로 안정화되었다. 이는 이전 연구와 일치하는 결과이다.
• 계산된 $ t = 100 $ 시점의 파동 진폭(3.76)은 이전 연구들 [10, 25, 34]의 결과와 매우 유사하여 이 방법의 신뢰성을 확인하였다.
• 특히 파동 전파 중 진폭 오차를 최소화하는 데서 콜로케이션 기반 지수 B-스플라인 방법보다 뛰어난 정확도를 보였다.
• $ C_1 $, $ C_2 $, $ C_3 $ 의 보존은 장시간 시뮬레이션 동안 안정적으로 유지되었으며, $ C_1 $ 이 8.94 에서 78.26 으로 약간 증가하는 것으로 나타나 방법의 강건성을 보여주었다.