动网格和overset物体速度快速震荡导致崩溃的问题

[tower]

@xiezhuoyu 我也询问过别人，给的建议也是加密网格，但是我试着加密了下，还是相同的错误。 其实我不是很理解的是，从我po在上面的log文件来看压力虽然迭代的次数多，但是最终残差不是已经满足条件了吗，为何压力积分后的力还会不准确，出现这种震荡问题呢

[李东岳]

最近在做一个波浪中物体二自由度的运动案例，最开始使用普通的动网格来计算的，因为计算幅度可以接受，但是结果却发现在开始计算0.2s左右，计算就发散了，发散的原因是物体Z方向上的速度急剧增大，

这个算例要是小的话，可以发上来给你看看

[tower]

多谢@东岳老师，我把动网格和overset的案例文件都上传了，您看一下。
动网格我是用的waves2foam+v1712的组合 overset用的1906的版本
overset： oversetMesh.zip
动网格： dynamicMesh.zip

[tower]

各位，
我尝试增大PIMPLE里面循环次数来增加p_rgh的迭代次数，但是发现每次迭代的次数还是很多 不是一个递减的趋势
如下：

DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.72710724614e-07, Final residual = 8.72514217379e-09, No Iterations 382
time step continuity errors : sum local = 7.03336340027e-09, global = 1.9150124018e-09, cumulative = 3.04929773866e-07
DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.69104939702e-07, Final residual = 1.44876414351e-09, No Iterations 335
time step continuity errors : sum local = 1.11124908301e-09, global = 1.95050292461e-10, cumulative = 3.05124824159e-07
DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.2782727071e-07, Final residual = 1.06023228234e-09, No Iterations 607
time step continuity errors : sum local = 7.68585702259e-10, global = -3.36686161497e-10, cumulative = 3.04788137997e-07
DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.47947887396e-07, Final residual = 9.72671115058e-10, No Iterations 464
time step continuity errors : sum local = 6.73366904995e-10, global = 7.62059709157e-11, cumulative = 3.04864343968e-07
smoothSolver:  Solving for epsilon, Initial residual = 0.000146170038677, Final residual = 5.6866402182e-09, No Iterations 8
smoothSolver:  Solving for k, Initial residual = 0.00238502525139, Final residual = 9.25193944196e-09, No Iterations 20
ExecutionTime = 689.06 s  ClockTime = 695 s

尝试修改为层流之后，连续迭代，次数基本成一个下降趋势，但是受力情况还是不变，错误依旧

如下：

Updating StokesII wave model for patch inlet
DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.33810842593e-07, Final residual = 1.05803239609e-09, No Iterations 95
time step continuity errors : sum local = 2.20723157353e-08, global = 6.73895228141e-09, cumulative = 7.30414848751e-06
DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.58642796578e-09, Final residual = 9.05674029848e-10, No Iterations 23
time step continuity errors : sum local = 1.87310827339e-08, global = 5.04129159013e-09, cumulative = 7.3091897791e-06
DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.88155187954e-09, Final residual = 9.84985998753e-10, No Iterations 1
time step continuity errors : sum local = 1.99973561049e-08, global = 4.66563789626e-09, cumulative = 7.313855417e-06
DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.11812493747e-09, Final residual = 9.39334225107e-10, No Iterations 1
time step continuity errors : sum local = 1.91664654171e-08, global = 4.75521938246e-09, cumulative = 7.31861063638e-06
ExecutionTime = 8152.77 s  ClockTime = 8159 s

[李东岳]

我这没装wave2foam，所以.. 看看其他大佬有什么意见

[xiezhuoyu]

虽然我算例的浮体和你不一样，但是问题很类似：在波浪还没有传播到浮体附近之前，浮体就受力振荡，p_rgh的迭代次数不是递减的趋势，感觉我就是在踩你踩过的坑 @tower
overset： oversetMesh.zip 这个算例用的是OF1906自带的造波消波模块吧，不需要waves2Foam造波消波，@东岳 老师能不能给我们一点建议o_0，看看应该从什么角度着手尝试解决这个问题？谢谢！
也希望有其他前辈能够就这个问题指点一下，缩短小白踩坑的时间，感谢！！

[CFDngu]

我用overset计算过浮体大幅度运动的问题，没遇到你们说的问题。
第一种可能是设置问题，可以仔细检查一下算例的设置，包括overset网格的设置；
第二中可能就是waves2Foam与overset网格结合的问题。这个我也不确定，因为我之前算的时候采用了自己写的造波程序。

[xiezhuoyu]

@CFDngu 谢谢建议！！

[卡洛]

老哥，请问最终问题怎么解决的，我最近也遇到这个问题了

[OItoCFD]

老哥，我来歪个楼，问你一下，wave2foam作用是啥？原版OpenFOAM不就有各种波浪么，或者自己自建个波浪边界条件也行吧，为啥还要额外一个wave2foam呢？

[jlx]

@tower 老哥，请问下问题解决了么

[hachikoi]

老哥们，你们是如何解决这个问题
我在使用两套overset做浮体自衰减的时候也遇到类似的问题

[wangwang]

你好，请问这个问题解决了吗？

[wzlr2023]

大家好，请问这个问题解决了吗？我用overinerDyMFoam算一个二维方块浮体也遇到了这个问题，更改了很多参数都没解决。。。