OpenFOAM 8中的波浪边界是否支持纯波浪?
-
@卡洛 在 OpenFOAM 8中的波浪边界是否支持纯波浪? 中说:
我没用过of8波流耦合,我用过1912自带的造波,waves2Foam,olaFlow都测试过上面说的湍流模型没问题
非常感谢,我试了以下,发现这个模型不是很稳定。原来用kOmega SST可以算的,改用kOmegaSSTbuoyang后很快就发散了。我在paraFoam里看了一下,主要是波峰前方的波面上出现了过高的k。你碰到过这种情况吗?
-
放个阶段性不完全的解决方案,参考这里,我从头到尾就没想到ddt在波浪模拟上会有这么大的影响。
-
@tidedrinker 朋友 这里这个设置会有用吗 https://www.openfoam.com/news/main-news/openfoam-v20-12/solver-and-physics
- Incompressible non-uniform density turbulent model for VOF
写的是减少相界面处湍流粘度过度而导致的波阻尼
-
这个帖子其实包括两个问题,首先是OpenFOAM自带的的波浪边界是否支持纯波浪。目前我基本没有头绪,可以确定的是波可以造出来,但波形完全不对,Airy波,二阶和五阶Stokes都不行。
然后是波浪衰减,我已经基本搞明白了改哪些参数有用了。第一个是fvSchemes里的ddt,一般默认是Euler,波浪这里最好改成Crank-Nicolson,也就是算例waves里的设置。这样的话时间步长在1/250波浪周期就有不错的结果(看趋势或许更长也可以,但我没试),如果用Euler可能要1/1000以下才行。我上边贴了一篇论文,里边有详细的比较。Value的话0.9就行,再高提升有限而且不稳定。
第二个是我波浪参数有问题,要么是波高太大,要么是波陡太大,总之增大波长或者降低波高(交界处网格尺寸不变)都能有效改善衰减,可能是因为这个波浪参数非线性太强。同时,将Airy波(一阶波)改成Stokes波也会有明显改善,但二阶波和五阶波差别不大。
现在的话,水深和波长不变,波高降到2.4,改用Stokes二阶波和Crank-Nicolson,连波流耦合导致的波浪变形一起算上,波高衰减不超过3%,已经能用了。
而且,如果用outletPhaseMeanVelocity做出口边界,用不用消波好像没区别。
-
@tidedrinker 谢谢这些有用的信息 那目前为止 开紊流模型 还会有影响吗?还是说还是得直接层流模型来算?
-
@tidedrinker 好的谢谢 请问另外几个RANS模型下的话都有小问题吗?
-
@tidedrinker 好的谢谢 明白
-
@tidedrinker 请问纯波浪问题有解决吗?我在学习wave算例的时候遇到相同的问题,如果右侧边界条件不能设置成出口的话,水位会一直上升,想试一下纯波浪能不能解决这个问题,但是不知道要设置什么样的边界条件
-
@funky番 没有解决,至少在基金会的框架下到of10都没解决,即使of10宣称可以搞纯波浪,照样没戏。
替代方案是用ESI的版本,我现在做纯波浪的算例,比如说越浪,就用of2112,挺好用。但也有问题,就是越浪之后入口平均水位会倾向保持不变,也就是说没有水的出口的话,计算域内的水会越来越多。在海工上问题不大,但重现实验如果水槽不大越浪很强可能会有问题。
如果要在纯波浪和波流耦合之间切换的话,我现在是用基金会版本加载olaflow的库取代掉controlDict里的libwaves,然后用边界不用求解器,现在看上去还行。
-
@tidedrinker 好嘞好嘞谢谢了,我去安装个ESI版本试试
-
@tidedrinker 前辈你好,我在使用olaFlow做一个二维的狭缝共振,但是把两个浮箱加上之后造出的波浪不向前传递是什么原因呀,详细的我放在了这个问题里面问题地址