[논문 리뷰] Numerical solutions of Korteweg de Vries and Korteweg de Vries-Burger's equations using computer programming

Mehri Sajjadian|arXiv (Cornell University)|2012. 09. 09.

Nonlinear Waves and Solitons참고 문헌 11인용 수 10

한 줄 요약

이 논문은 Korteweg-de Vries(KdV) 및 Korteweg-de Vries-Burgers(KdVB) 방정식을 수치적으로 해결하기 위해 sinc-보간법과 θ-가중 유한차분 방법을 조합한 방법을 제안한다. 이 방법은 결과로 생기는 행렬 체계를 Mathematica를 사용해 풀어 최대 절대 오차가 1.85×10⁻¹² 수준까지 낮춘 높은 정확도를 달성하며, KdV 방정식에 대해 세 가지 보존법칙을 검증한다.

ABSTRACT

In this paper, numerical and solitonic solutions of Korteweg de Vries(KdV) and Korteweg de Vries-Burger's (KdVB) equations with initial and boundary conditions are calculated by sinc-collocation method. The basis of method is sinc functions. First, discretizing time derivative of KdV and KdVB's equations using a classic finite difference formula and space derivatives by θ- weighted scheme between successive two time lev- els is applied, then Sinc functions are used to solve these two equations. Mathematica programming is used to solve matrix representation of these equations. KdV equation describes behavior of traveling waves which is a third order non-linear partial differential equation (PDE). Maximum absolute errors are given in Tables. The figures show approximate solutions of these two equations. Three conservation laws for KdV's equation are obtained.

연구 동기 및 목표

• 초기 조건과 경계 조건을 갖는 Korteweg-de Vries(KdV) 및 Korteweg-de Vries-Burgers(KdVB) 방정식을 해결하기 위한 강력한 수치적 스킴을 개발하기 위해.
• 이러한 비선형 편미분방정식에 대해 고정밀 수치적 해를 도출하기 위해 sinc-보간법을 적용하기 위해.
• 제안된 수치적 방법을 사용하여 KdV 방정식의 세 보존법칙이 어떻게 유지되는지 검증하기 위해.
• 다양한 매개변수를 활용한 수치 실험을 통해 스킴의 효율성과 정확도를 평가하기 위해.
• Mathematica 기반 행렬 해법과 오차 분석을 통해 이 방법의 계산적 매력도를 입증하기 위해.

제안 방법

• KdV 및 KdVB 방정식의 시간 도함수를 고전적인 유한차분 공식을 사용해 이산화하기 위해.
• 안정성과 정확도를 확보하기 위해 연속적인 시간 단계 사이의 공간 도함수에 대해 θ-가중 스킴을 적용하기 위해.
• 보간을 위해 사인 함수를 기저 함수로 사용하여 공간 근사에서 스펙트럴 정확도를 달성하기 위해.
• 반연속 시스템의 행렬 표현을 구성하고 Mathematica를 사용해 수치적으로 해를 구하기 위해.
• 다양한 매개변수 설정과 도메인 구성 조건에서 KdV 및 KdVB 방정식에 대해 이 방법을 구현하기 위해.
• KdV 방정식의 보존법칙을 계산하고 분석하기 위해 이 방법을 활용하여 질량, 운동량, 에너지의 불변성을 검증하기 위해.

실험 결과

연구 질문

• RQ1 sinc-보간법과 θ-가중 시간 이산화를 조합한 방법이 KdV 및 KdVB 방정식을 해석할 때 높은 정확도를 달성할 수 있는가?
• RQ2 제안된 스킴에서 최대 절대 오차는 시간 간격과 공간 그리드 크기의 변화에 따라 어떻게 변화하는가?
• RQ3 KdV 방정식의 세 보존법칙이 수치적 스킴 하에서 얼마나 잘 유지되는가?
• RQ4 분산(μ), 점성(ν), 비선형성(ε) 등의 매개변수가 다양할 경우 이 방법의 성능은 어떻게 되는가?
• RQ5 기타 수치적 방법과 비교했을 때 sinc-보간법 접근의 계산 효율성과 수렴 행동은 어떠한가?

주요 결과

• KdVB 방정식에 대해 높은 공간 및 시간 해상도(n=16, δt=0.00001) 조건에서 최대 절대 오차가 1.85×10⁻¹² 수준까지 낮아졌다.
• KdV 방정식의 경우 질량, 운동량, 에너지 보존법칙이 모두 유지되어 스킴의 물리적 일관성을 확인하였다.
• 수치적 결과는 L∞ 및 L2 오차가 시간이 지남에 따라 점차 증가함을 보여주었으며, 한 테스트 케이스에서 T=0.0009일 때 L∞ 오차가 1.66×10⁻¹¹에 도달하였다.
• μ=0.001, ν=0.001, ε=1, θ=0.5 등의 다양한 매개변수 조합에서도 높은 정확도와 안정성을 보였다.
• 그림들은 KdV 및 KdVB 방정식에 대해 부드럽고 정확한 근사 해를 보여주며, 솔리톤 유사 행동이 유지됨을 확인하였다.
• sinc-보간법에서 발생하는 행렬 체계를 Mathematica로 효율적으로 처리함으로써 실용적 구현 가능성을 확인하였다.