计算的稳定无旋表面孤独的重力波解欧拉方程(均匀,不可压缩和完美的液体)。定义的波弗劳德数Fr和结果是大约15数字准确。方法适用于所有但最高的波浪,即振幅/深度比小于0.796。
剧情简介:
SolitaryGravityWave (Fr, [], 1);%绘制结果只
[z, ws, fs, SWP] = SolitaryGravityWave (Fr);%输出结果在表面和参数
SWP [Z, ws, fs, W, F P) = SolitaryGravityWave (Fr, Z);%表面和体积输出
SWP [Z, ws, fs, W, F P) = SolitaryGravityWave (Fr, Z, 1);

输入:
Fr:弗劳德数(必须是一个标量)。
Z:复杂的横坐标字段所需的内部流体(默认Z = [])。
Z应该严格低于表面。1 < =图像放大(Z) <η(真实(Z))
y =埃塔(x)在自由表面的方程。
PF:情节国旗。如果绘制PF = 1最终结果,如果PF ~ = 1没有绘制(默认)。

输出(无量纲量):
z:表面复杂的横坐标,即。x + i *埃塔。
ws:复杂的表面速度,即。,u - v *。
fs:复杂的潜在的表面,即。,+我*ψφ。
SWP:孤波参数,即。
SWP(1) =波振幅,马克斯(eta)
SWP(2) =波质量
SWP(3) =循环
SWP(4) =冲动
SWP(5) =动能
SWP(6) =势能
W:复杂的速度在大部分横坐标Z。
F:在大部分复杂的潜在横坐标Z。
P:压力在大部分横坐标Z。
答:复杂的加速度在大部分横坐标Z (= dW / dt)。

例子:
z = SolitaryGravityWave (1.25);
情节(真实(z),图像放大(z))

编辑更多细节和m文件
看看转载:
http://hal.archives-ouvertes.fr/hal-00786077/