网格变细后液面不正常

[李东岳]

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

[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]

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