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



  • 各位前辈好,
    本人最近在学习使用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步的时候已经振荡到无界然后发散了。而且浮力有越来越小的趋势。:143: 目前没有对矩阵求解和方程离散格式做任何改动,不太明白如此振荡的原因在哪里?请有经验的前辈指点一二,感激不尽。:xinxin3: mesh_setting2.png mesh_setting.png 1619345017(1).png 1e746636d8aed9b43b83169178d9625.jpg



  • 可以贴上alpha.water流场云图看看,而且最后一张图的曲线代表什么意思?(看不太清楚)



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


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1