网格变细后液面不正常

[Joann]

在算一个船在静水中的阻力系数，船不动，水流速1.18米每秒，按道理应该不会有很大的自由液面波动。
在ddtSchemes用的是LocalEular所以算是用的LTSInterFoam，两套网格，前一套网格量235200，后一套细一点是672000，前一套的自由液面挺正常，阻力系数也逐渐收敛。后一套网格算得的自由液面运动很大，阻力系数变化一直在大跳就没有收敛趋势。两个计算除了网格blockMesh里变了其他都不变。
如图，同一时间的两套网格的液面比较
以及较细网格的大跳阻力系数，令人迷茫？？？这不是没收敛的问题这就没打算收敛
为啥会这样？

[Joann]

@Joann 图上传不正常，我过会再上传试试

[Joann]

较粗的网格：


较细的网格：

[李东岳]

细网格看起来像是发散了... 原因不详

[Joann]

@李东岳 实在没找到原因，我把细网格的算例上传，老师如果有时间可以麻烦帮我看看吗？
链接：https://pan.baidu.com/s/1D-genYmkrjnflowqFpA6zQ
提取码：qhwt

[Joann]

@joann 问题应该就在网格上，因为只有网格变了，猜想是不是细网格的z方向上分得太细了和x，y方向上的长度相差较大，所以就发散了····

[一朵苍凉]

@joann
初始相分数场有没有检查一下，我之前也有遇到过这种问题，相分数场设置正确就正常了。

[Joann]

@一朵苍凉 是在setField之后看一下0秒时的alpha是否正常？

[一朵苍凉]

@joann
嗯，然后计算几秒后也可以再查看一下。

[Joann]

@一朵苍凉 好的我试试暂停几次看看alpha场，您之前是遇到初始alpha场不正常所以后面也发散了吗？
我看alpha场发现在跑了一会后出现船头轻微甲板上浪的情况，是否是因为出现过甲板上浪影响到后面的计算，于是就液面异常并且发散了。

[李东岳]

细网格算例网格数多少？

[Joann]

@李东岳 blockMesh里的ncell好像是67万左右，粗网格是23.5万左右。

[李东岳]

@joann 那个太难debug了 跑不动..

[Joann]

@李东岳 我再找找原因，不行就换网格

[Joann]

@李东岳 说到换网格，想问一下应该怎么选取正确的locationInMesh这个点。

经常卡在snappyHexMesh里的locationInMesh的这个点的选取上，这个点要求不能在面上要在单元体内，但是我选取了不在面上的点也会报错，有几次遇到的报错是说有多个锚点，可在sHMdict里确实就一个locationInMesh.

[一朵苍凉]

@joann 我之前遇到过相分数场不对导致受力振荡然后发散的情况，静水阻力一般甲板不会上浪吧，你还可以检查一下其他边界条件是不是有问题。

[Joann]

突然想到这个帖子，通过调整边界层厚度可以达到收敛了。