[论文解读] The Trigonometric Cubic B-spline Algorithm for Burgers' Equation

İdris Dağ, Özlem Ersoy Hepson+1 more|arXiv (Cornell University)|Jul 21, 2014

Advanced Numerical Analysis Techniques参考文献 15被引用 6

一句话总结

本文提出一种结合三角函数三次B样条配置法与Crank-Nicholson格式的数值方法，用于求解Burgers方程，在捕捉陡峭前缘和激波解方面表现出高精度与稳定性。数值实验表明，该方法在低黏性情况下性能优越，其结果与精确解高度吻合，且在误差范数方面优于现有的基于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

研究动机与目标

• 开发一种鲁棒的数值格式以求解Burgers方程，尤其适用于产生陡峭前缘和激波波的低黏性情况。
• 研究三角函数三次B样条（CTB）在非线性PDE配置方法中的有效性。
• 从精度与收敛性角度，比较所提出的CTB配置法与现有B样条法及Galerkin法的性能。
• 展示该方法在时间演化过程中保持解光滑性与稳定性的能力。
• 提供一种基于三角函数基函数的可靠替代传统B样条配置法的方法。

提出的方法

• 空间域采用三角函数三次B样条（CTB）基函数进行离散化，其具有二阶连续可微性，可提供高阶精度。
• 时间方向采用Crank-Nicholson格式进行离散，确保二阶精度与无条件稳定性。
• 采用配置法，其中解被近似为CTB函数的线性组合，并在选定的配置点上强制残差为零。
• 得到的非线性方程组通过迭代方法求解，且在每个时间步长施加边界与初始条件。
• 该方法利用CTB配置矩阵的对角结构，使得每个时间步长的线性系统可高效求解。
• CTB基函数通过引入形状参数的三角函数定义，确保光滑性，并提升对振荡或陡峭解的逼近性能。

实验结果

研究问题

• RQ1 与标准三次B样条相比，三角函数三次B样条是否能在存在激波波的情况下为Burgers方程提供更精确、更稳定的解？
• RQ2 当黏性参数λ较小时，CTB配置法在捕捉陡峭前缘解方面表现如何？
• RQ3 CTB配置法在不同空间与时间离散化下的收敛行为及误差分布如何？
• RQ4 与现有B样条配置法及Galerkin法相比，该方法在精度与计算效率方面表现如何？
• RQ5 使用三角函数基函数是否能改善非线性PDE中行波解的逼近效果？

主要发现

• 所提出的CTB配置法具有高精度，如表5所示，在λ = 0.01且h = 1/36时，最大绝对误差处于10^-3量级。
• 在t = 0.5时，该方法的数值解与精确解高度一致，绝大多数空间点的相对误差低于0.002。
• 当λ = 0.001时，该方法成功捕捉到向右传播的尖锐激波前缘，与理论分析一致。
• L∞误差范数保持有界，并随空间分辨率提高而减小，表明方法具有收敛性与稳定性。
• 在陡峭梯度区域，该方法在精度上优于或至少不逊于现有的三次B样条配置法与Galerkin法。
• 图形结果（图5–8）证实了该方法在时间演化过程中保持解光滑性，并能准确解析激波轮廓。