计算钝体涡激振动时发散（pimpleDyMFoam）

[李东岳]

@wying 动网格这面最大的问题就是网格变形过大导致的库朗数特别大，需要大幅度降低时间步长保证收敛。你需要保证时间步长充分的小。

debug的时候要检查边界条件是否有问题可以先忽略网格变形，先直接算看能不能算出大体的结果。如果正确的话说明不是边界条件的问题。

[WYing]

@李东岳 感谢李老师的回复！我按照您的建议，先计算了钝体绕流（不加动网格）的情况，发现能够得到合理的流场和大致结果，这应该说明了边界条件设置没有问题。那么就是动网格的设置出错了？我计算的是圆柱+不同长度分流板的流致振动，奇怪的是，大部分算例可以正常计算，只有某些长度的分流板会导致结果发散。所有算例的网格划分方式和参数设置都一样的，真是百思不得其解。。。

[李东岳]

@wying 可能某些长度的板子网格变形太大了

[WYing]

@李东岳 我也想过这个可能性。如果是因为网格变形增大导致的发散，那么至少在开始一段时间内变形较小的时候计算能正常进行。但是这些发散的case从一开始的流场就有问题，而不是说随着网格变形增大，逐渐发散。。

[李东岳]

@wying 如果发散的算例网格比较少、也不涉密、可以发上来我给你看看

[WYing]

@李东岳 感谢您的回复！ crash.zip 这是我的case文件，由于大小限制，文件夹内没有包含最后一个时间步，但是大概在time=1.8所有停止计算。欢迎您提出的任何意见或建议！

[李东岳]

@wying 不好意思明天我就要去杭州了... 国庆后回来 只能看看有没有其他大佬能给你看看了，

[WYing]

@李东岳 谢谢李老师的回复！我尝试使用高版本的OF计算，发现相同设置下算例可以正常计算，不确定这是不是低版本的OF涉及动网格时存在bug。不过有个case在所有版本OF下都发散。。。能麻烦李老师有空的话帮忙看一下吗？算例文件case.zip （大约在time=1.43发散）十分感谢！！

[李东岳]

好的 我给你瞅瞅最新的这个case.zip

[cresendo]

@wying 试试看对p不要进行松弛

[WYing]

@cresendo 谢谢您的建议！去掉p松弛之后，算例确实可以计算下去，但是CFL保持在20左右，且速度场和压力场也明显有问题。所以还是不太对，不过还是感谢您的帮助！

[cresendo]

@wying 那应该就是动网格的问题了，试试在dynamicMeshDict里将松弛因子调成0.4试试

[WYing]

@cresendo 请问您说的dynamicMeshDict里的松弛因子是 accelerationRelaxation 0.9 和 accelerationDamping 0.95 吗？

[cresendo]

@wying 前者

[WYing]

@cresendo 我按照您的建议，调整了accelerationRelaxation 0.4，计算可以正常进行，CFL数小于1， 流场看上去也正常。

但有一个小疑问，通过查询dynamicMeshDict文件的官方描述，得知accelerationRelaxation直接减小刚体的加速度，比如在某时间步solver求出的加速度为10m/s^2，如果设置accelerationRelaxation=0.4，那么实际应用的加速为4m/s^2，这是一个很大的松弛。而官方描述中建议的accelerationRelaxation范围是0.9-1，并且提到 “Be careful with this accelerationRelaxation. Too low of a value will mean that the Body does not respond to the fluid forces correctly”。因此我担心设置较小的加速度松弛会不会对结果产生太大影响？因为没有实验结果作为参照，所以无法判断结果的准确性。不知道您一般是如何处理这种误差的？谢谢！

[cresendo]

@wying 你说的这种松弛方式是不正确的，应该是上一时间步的0.6加上当前时间步的0.4，这样松弛以后获得最终的加速度。至于wikidict上说的这个建议，松弛应该影响是的最终稳定的速度，但是可能的确会导致响应运动不是那么准确。因此，还是建议你寻找一个标模实验验证一下（这种实验很多）

[WYing]

@cresendo 我重新查看了source code，aRel的松弛方式确实如您所说，感谢指出错误！我对比了不同aRel的结果，发现Rel = 0.4-0.8都可以正常运行且结果几乎相同，但是aRel >0.8之后会发散。我不确定为什么wikidict上说太小的aRel会造成不准确的响应，至少我的结果并不受aRel取值的影响。相反，恰当的aRel可以使本来发散的case收敛。我之前没有指定过aRel，也就是取默认值1.0，因此比较好奇为什么您会设定aRel=0.4而不是其他值，是根据模拟的经验吗？谢谢！

[cresendo]

@wying 0.4是比较保守的设置，我也是根据自带的DTCHull算例参考的。但是对于涡激运动问题，以我个人2D的尝试，0.9是可以跑的。我本人是做涡激运动的，这是我的songtao.chen@sjtu.edu.cn，有问题可以一起交流~

[WYing]

@cresendo 在 计算钝体涡激振动时发散（pimpleDyMFoam） 中说：

songtao.chen@sjtu.edu.cn

谢谢！非常感谢您的建议！