[论文解读] Numerical solutions of Korteweg de Vries and Korteweg de Vries-Burger's equations using computer programming

Mehri Sajjadian|arXiv (Cornell University)|Sep 9, 2012

Nonlinear Waves and Solitons参考文献 11被引用 10

一句话总结

本文提出一种结合sinc配置法与θ加权有限差分格式的数值方法，用于求解Korteweg-de Vries（KdV）方程和Korteweg-de Vries-Burgers（KdVB）方程。该方法利用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配置法，为这些非线性PDE提供高精度的数值解。
• 通过所提出的数值方法验证KdV方程中三个守恒律的保持性。
• 通过不同参数设置的数值实验，评估该格式的效率与精度。
• 通过基于Mathematica的矩阵求解与误差分析，展示该方法在计算上的优越性。

提出的方法

• 使用经典的有限差分公式对KdV与KdVB方程中的时间导数进行离散化。
• 在相邻时间层之间对空间导数应用θ加权格式，以确保格式的稳定性和精度。
• 使用sinc函数作为基函数进行配点，实现在空间逼近中的谱精度。
• 构建半离散化系统的矩阵表示，并利用Mathematica进行数值求解。
• 在不同参数设置与区域配置下，将该方法应用于KdV与KdVB方程。
• 利用该方法计算并分析KdV方程的守恒律，验证质量、动量与能量不变量的保持性。

实验结果

研究问题

• RQ1 sinc配置法结合θ加权时间离散化是否能在求解KdV与KdVB方程时实现高精度？
• RQ2 在所提出的格式中，最大绝对误差如何随时间步长与空间网格尺寸的变化而变化？
• RQ3 所提出的数值格式在多大程度上保持了KdV方程的三个守恒律？
• RQ4 当色散系数（μ）、黏性系数（ν）与非线性系数（ε）变化时，该方法在不同参数区域中的表现如何？
• RQ5 与其它数值方法相比，sinc配置法的计算效率与收敛特性如何？

主要发现

• 在空间与时间分辨率较高的情况下（n=16，δt=0.00001），所提方法求解KdVB方程的最大绝对误差低至1.85×10⁻¹²。
• 对于KdV方程，该方法成功保持了质量、动量与能量三个守恒律，证实了格式的物理一致性。
• 数值结果表明，L∞与L2误差随时间逐渐增长，在某一测试案例中，T=0.0009时L∞误差达到1.66×10⁻¹¹。
• 该方法在多种参数组合下（如μ=0.001，ν=0.001，ε=1，θ=0.5）均表现出高精度与稳定性。
• 图示展示了KdV与KdVB方程的光滑且精确的近似解，且孤子样行为得以保持。
• 在Mathematica中的计算实现高效处理了配置法产生的矩阵系统，证实了该方法在实际应用中的可行性。