CFD中文网

最近在做风浪相互作用的模拟时遇到一个问题：在气液剪切过程中，湍流的研究是非常重要的，我在加入湍流模型之后，计算结果耗散特别严重，但是没有气体作用（nowind）时，SSTk-w湍流模型能够较准确的模拟，从下面的水位图中也可以看出来。并且湍动能k和湍流粘度nut在出口outlet逐渐向计算域内部扩散导致波高衰减（见下图），但是不知道具体原因。单纯的平底地形结果也没问题，但是遇到到这种斜坡地形就有衰减。
目前我已经尝试过of自带的几乎所有RANS类的湍流模型，计算结果均耗散严重（以SSTk-w为例），也尝试过对边界条件和壁面边界层的各种处理方式，还是效果不理想，想请教一下各位大佬，想解决这个问题，我还应该从哪些地方入手？
二维模型如图：
模型图.png
水位图：
水位图.png
不同时刻湍动能k分布：
k-12s.png
k-30s.png
k-60s.png
不同时刻湍流粘度nut分布：
nut-12s.png
nut-30s.png
nut-60s.png

https://www.dataphysics-instruments.com/us/knowledge/understanding-interfaces/surface-tension/
这里有一句

Thus, for liquids the surface tension and the surface energy are equal

至于动能，k就是湍动能啊，如果不只是要湍动能，那k怎么来的肯定要涉及动能啊。

@东岳 李博是不是在植入这个算法，我等你:ok:

@东岳 我现在接触的还很浅。。以后碰到这方面问题再来请教~~

主任招博士后不，我给您发个简历

或者你应该在intU[1]=...之前加个if(cellID > -1).

@haining-luo
我现在用ansa处理极复杂三维模型，保存.mesh格式面网格，用fluent meshing或star ccm生成多面体网格或切割体网格，挺方便。
如果模型太复杂，搞科研可以用icem，工程上不太建议用ICEM生成非结构网格。

@hotairballon1997

在您回帖的链接中好像没有看到关于这个问题的解答，

对呀，要分类讨论，就像你说的，是用 PBM，还是拉格朗日。
我发的那个链接里 MPPIC 和 DPMFoam 中用到的粒子都是粒子包（particle parcel），这种用粒子包去代替粒子的方法我没了解过，但肯定已经不是物理粒子了。所以要看哪个层面，哪种模型的碰撞聚合了。

@东岳 多谢老师！

@shrine 没关系，无关紧要这些文件

已经解决。
（在一个视频中提到这一点）网格统计问题有三种情况：
第一种，是显示问题的网格数目很少，就几个，占总体的0.0002%。对于这种情况直接更改限定值。
第二种，是显示问题的网格数目很多，占总体的超过10%，甚至更多，对于这种情况需要更改MeshData的参数设置，重新对网格进行绘制。
必须保证，不能出现负网格。

你这个是三维的LES单相流模拟，应该不难分析，常见的不收敛的原因有很多，上面大神们说的时间步长、网格质量、松弛因子都是很不错也很常见的调试首选项，还有就是边界类型也挺关键，再就是几何域合不合适，上下游延伸段够不够，比如存在很狭小的流动域之类的、动网格更新策略，都可以作为改进收敛性的切入点。

确实和概率论有关，谢谢，记录一下

\begin{equation}
\left( \begin{matrix}
k\\
x
\end{matrix}\right )=\frac{k!}{(k-x)!x!}
\end{equation}

@李东岳
谢谢老师，虽然只是一个小建议和思路，却给我节约了大量的时间，目前已经在OF7中调通了mutiphaseInterFoam求解器。:xiexie:

报的错都没贴出来，人眼debug啊:135:

从m_lg=0开始试，一项一项往里加东西，看加到哪一项出的问题

@下里巴人 看到这个帖子想起我之前做那个流化床模拟，加了跟管子就死活不收敛了，估计和你这个一样的问题。
话说有尺度差异的时候出来的这个网格好呆……明明圆形的楞给画成多边形了。
我想起Gmsh他们组最近有个博士答辩的东西就是做curve mesh的，不知道能不能用在计算上了。。

哦，明白了。

@cyberk 抱歉，这个问题我没有遇到过，但我看了你之前用ICEM画的网格，这个六叶Ruston桨的搅拌槽的网格不太复杂的呀，这个你可以试试只画1/6个区域，然后对称过去，画完整体之后再去openfoam中用topoSet划分搅拌的区域

Really? 那得谢谢王博士！

@coolhhh 好的，谢谢大佬~