船舶阻力求解发散求助(静网格)

[一朵苍凉]

各位前辈好，
本人最近在学习使用OpenFOAM进行船舶阻力模拟，参照的是DTCHull算例的设置，对网格进行了重新划分，湍流变量的初始值进行了重新计算。参照的是https://scholar.google.com/scholar?q=Uncertainty+analysis+in+ship+resistance+prediction+using+OpenFOAM&hl=zh-CN&as_sdt=0&as_vis=1&oi=scholart
该论文中提供的公式，对KOmegaSST模型的k、omega和nut值进行了重新计算。船的航速为1.057m/s，计算得到k的初始值为0.00067，omega的初始值为59.07，nut的初始值为0.00002。DTCHull这个算例在限制船舶所有自由度的情况下，差不多求解50步就基本收敛了，而且几乎没有振荡。而我的算例阻力振荡非常严重，求解到20000步的时候已经振荡到无界然后发散了。而且浮力有越来越小的趋势。 目前没有对矩阵求解和方程离散格式做任何改动，不太明白如此振荡的原因在哪里？请有经验的前辈指点一二，感激不尽。

[cresendo]

可以贴上alpha.water流场云图看看，而且最后一张图的曲线代表什么意思？（看不太清楚）

[一朵苍凉]

最后一张图是X方向的受力，也就是船的阻力。
谢谢回复 ，如您所言，setFields确实有问题，我设置的值正好等于流体域边界值，当时没有检查其实边界处的相分数没有完全赋值。将setFields中的box范围改成稍大于流体域就可以正常求解了。我本来还以为是数值振荡这么严重呢，都想离散格式精度调低看看了。现在看来，还是物理模型导致的振荡才更严重。