[论文解读] Solitary wave solutions of FKPP equation using Homogeneous balance method(HB method)

Yirui Yang, Wei Kou+2 more|arXiv (Cornell University)|Sep 24, 2020

Fluid Dynamics Simulations and Interactions被引用 4

一句话总结

本文应用齐次平衡（HB）方法推导费舍尔-科尔莫戈罗夫-彼得罗夫斯基-皮斯丘诺夫（FKPP）方程的精确孤波解，该方程为非线性反应-扩散方程。通过基于HB原理构造函数变换并求解所得常微分方程，作者获得了以双曲正割和正切函数表示的显式解，并通过图形分析和参数约束进行了验证。

ABSTRACT

In this paper, we use the homogeneous balance(HB) method is used to construct function transformation to solve the nonlinear development equation||Fisher-Kolomogror-Pertrovskii-Piskmov equation (FKPP equation), then the exact solution of FKPP equation is obtained, and the solutionis converted into the form of isolated wave solution. Finally, we have shown the solution of FKPPequation in some picture forms, and the rationality of the solution in this paper can be verified byusing the pictures.

研究动机与目标

• 开发一种系统方法，利用齐次平衡（HB）方法求解非线性FKPP方程。
• 通过基于HB原理的函数变换，将FKPP方程转化为可解的常微分方程（ODE）。
• 推导以双曲函数（sech和tanh）表示的精确孤波解。
• 通过图形表示和与原方程的一致性检验验证解析解。
• 确定参数约束（k, c, a, b），以获得物理解释明确且数学上一致的解。

提出的方法

• 将HB方法应用于假设一个涉及函数φ(x,t)高阶偏导数的解形式，其中最高阶导数的阶数N = 2。
• 假设解的形式为u(x,t) = f''(φ)φ_x² + a f'(φ)φ_x + b，其中f(φ)通过非线性项与最高阶导数项之间的平衡关系确定。
• 求得函数f(φ)为f(φ) = -6 ln φ，该结果在HB原理下可导出可解的常微分方程组。
• 通过引入φ(x,t) = 1 + e^(kx + ct)，将解简化为包含sech²和tanh函数的形式，从而实现孤波的显式表示。
• 通过平衡最高阶项并确保与原FKPP方程的一致性，推导出参数k和c。
• 最终解表示为u(x,t) = -3/2 k² sech²[½(kx + ct)] + 3ak tanh[½(kx + ct)] - 3ak + b，其中特定参数值通过代入FKPP方程确定。

实验结果

研究问题

• RQ1 齐次平衡方法能否有效应用于推导FKPP方程的精确孤波解？
• RQ2 f(φ)和φ(x,t)的何种函数形式能满足FKPP方程中非线性项与最高阶导数项之间的平衡条件？
• RQ3 孤波解的显式解析形式在双曲函数（sech和tanh）下如何表示？
• RQ4 参数k, c, a, b如何影响孤波解的形状与动力学特性？
• RQ5 所推导的解能否通过代入和图形分析验证为精确解？

主要发现

• 齐次平衡方法成功地将FKPP方程约化为可解的常微分方程组，得到以双曲函数表示的精确解。
• 解的形式为u(x,t) = -3/2 k² sech²[½(kx + ct)] + 3ak tanh[½(kx + ct)] - 3ak + b，其中k和c由参数a和b确定。
• 当b = 0，c = 5/6，且k = -1/(6a)时，解变为u(x,t) = -1/(24a²) sech²[½(-x/(6a) + 5t/6)] + (1/2) tanh[½(-x/(6a) + 5t/6)] + 1/2。
• 当b = 1，c = -5/6，且k = 1/(6a)时，解为u(x,t) = -1/(24a²) sech²[½(x/(6a) - 5t/6)] - (1/2) tanh[½(x/(6a) - 5t/6)] + 1/2。
• 参数a被约束为±1/√6，以满足原FKPP方程，确保所推导解的一致性。
• 图形表示（图1）证实了解的孤波特性，显示了正负a值下局部化且稳定的波形。