[论文解读] A self-adaptive mesh method for the Camassa-Holm equation

Bao‐Feng Feng, Ken‐ichi Maruno+1 more|arXiv (Cornell University)|May 16, 2009

Nonlinear Waves and Solitons参考文献 48被引用 5

一句话总结

该论文提出了一种针对Camassa-Holm方程的自适应移动网格方法，能够动态地在解的梯度剧烈区域（如peakons和cuspons）集中网格点。该方法具有可积性，在连续极限下保持N-孤子解的性质，通过在演化非均匀网格上求解三对角系统，仅用约100个网格点即可实现高精度，显著优于标准方法在解析奇异解方面的表现。

ABSTRACT

A self-adaptive moving mesh method is proposed for the numerical simulations of the Camassa-Holm equation. It is an integrable scheme in the sense that it possesses the exact N-soliton solution. It is named a self-adaptive moving mesh method, because the non-uniform mesh is driven and adapted automatically by the solution. Once the non-uniform mesh is evolved, the solution is determined by solving a tridiagonal linear system. Due to these two superior features of the method, several test problems give very satisfactory results even if by using a small number of grid points.

研究动机与目标

• 开发一种能够准确捕捉Camassa-Holm方程奇异解（如peakons和cuspons）的数值方法，这些解因具有剧烈梯度而对标准格式构成挑战。
• 确保数值格式在连续极限下保持可积性，通过保留N-孤子解，与方程内在的可积结构保持一致。
• 通过动态调整网格，使网格点集中在解变化迅速的区域，而非使用均匀网格，从而降低计算成本并提高精度。
• 在具有挑战性的测试案例（包括soliton-cuspon和cuspon-cuspon碰撞）以及产生peakon列的非精确初值条件下，展示该方法的有效性。

提出的方法

• 该方法基于Camassa-Holm方程的可积半离散化，通过hodograph变换和因变量变换，从变形的修正KP族的双线性形式推导而来。
• 采用随时间演化的非均匀网格，网格点根据解的梯度自动调整，尤其在cuspon和peakon结构附近。
• 在每个时间步长，通过求解从半离散方程导出的三对角线性系统来计算解，确保计算效率。
• 采用两种时间积分格式：改进的前向欧拉法和经典的四阶龙格-库塔法，后者表现出更优的精度。
• 网格自适应机制确保在快速变化的局部区域（如类似冲击波或尖点特征）实现高分辨率，而无需预先知晓解的结构。
• 该方法通过构造保持可积性，因为半离散格式可容纳精确的N-孤子解，且当网格尺寸趋于零时，收敛于连续的N-孤子解。

实验结果

研究问题

• RQ1 自适应移动网格方法是否能在准确解析peakon和cuspon解的同时，保持Camassa-Holm方程的可积性？
• RQ2 该方法在使用少量网格点时，对soliton-cuspon和cuspon-cuspon碰撞等复杂相互作用的模拟表现如何？
• RQ3 该方法能否处理导致peakon列生成的非精确初值，特别是在无色散极限下？
• RQ4 改进的前向欧拉法与四阶龙格-库塔法在该系统中演化非均匀网格时的相对性能如何？

主要发现

• 该方法仅需约100个网格点即可实现高精度，远少于标准方法（如有限体积格式）为达到类似分辨率所需的4096个点。
• 在单个peakon及peakon-peakon相互作用情形下，该方法以极少的网格点准确捕捉了解的动力学行为，并保持了预期的解形态。
• 在cuspon-cuspon和soliton-cuspon碰撞情形下，该方法成功解析了复杂的相互作用，包括尖点和剧烈梯度的形成。
• 当以非精确初值和小色散参数（c = 90）初始化时，该方法在t = 200时生成了七至八个近似peakon结构的列，与无色散极限下的理论预期一致。
• 四阶龙格-库塔法在精度上优于改进的前向欧拉法，尽管两者对光滑解均有效，但四阶龙格-库塔法在cuspon解上因刚性而失效。
• 该方法在无色散区域能够稳健捕捉从光滑初值中出现的peakon列，表明其在具有奇异孤子的可积方程中具有更广泛的应用潜力。