怎么让交界面更加平滑

[asininno]

在二相流模拟中，捕捉界面后怎么让界面更平滑来避免计算的不稳定呢，openfoam中经常通过mules来求解相方程，它起到的什么作用？

[余正东]

@asininno[0_1511445949233_TutorialmultiphaseInterFoam.pdf](正在上传 100%) 可能有用

[asininno]

@余正东 你这发的看不见啊

[余正东]

@asininno 我的权限不够…………你加我吧（qq：768620698）

[yhdthu]

以我目前对基于vof界面捕捉方法的理解，mules采用的是pde based sharping interface method，这个方法就是在相方程中加入了一个守恒的压缩项，离散这个非线性项用的是fct通量限制器，代数重构，这方法捕捉的界面还是太粗糙，就说计算界面的法向量也就一阶，esi最新发布的isoAdvector方法碉堡了，跟youngs提出的PLIC还不是一个套路，做的几何重构，计算量还要比PLIC小，任意网格下都可以达到2Order的精度，目前正在拓展这个方法，很有趣

[wwzhao]

@yhdthu 动网格和库朗数限制是摆在isoAdvector前面比较大的两个问题。

[yhdthu]

@wwzhao 动网格还没考虑到，请问Co数是个什么说法？

[千里之行]

@wwzhao 动网格主要分两种，一种是整体做刚性运动，网格单元之间没有相对位移。另一种是单元有变形，有拓扑变化，或者有动态增减网格。 isoAdvector 应用于第二类动网格可能需要一些改动，那么第一类动网格呢？我目前正在准备研究这个，想请教一下您对此的看法。

[李东岳]

@yhdthu 在 怎么让交界面更加平滑 中说：

请问Co数是个什么说法？

因为是个显性的求解，方程2.几有说，刚在粗略看了一遍这个文章，才发现这个文章里面没有耦合起来啊！就是求个标量传输。文章结论还说可以扩展到其他CFD方面。那感情好了，有的搞。

动网格我也没关注。

[yhdthu]

@李东岳 前辈，我看了源码，我的理解哈，这个相方程和ns方程的耦合是通过dVf的计算，alpha是靠这东西更新的，rhoPhi也是由他计算来的

[yhdthu]

@李东岳 附上Roenby的另一篇文章，里面提及了如何向interFoam里耦合:laughing:
IsoAdvector: Geometric VOF on general meshes

[wwzhao]

@yhdthu 显式求解对Co数较为严格，基本要求所有网格的Co数小于1。

[wwzhao]

@千里之行 你说的第一类网格在Star-CCM+里面有的，但对于带自由面的问题，只限于小幅度纵摇（俯仰角）情况。纵摇较大的情况不好处理。

[litong189456]

@yhdthu 看了一下他的文章，发现就是基于不可压缩流体来推导的方程，您觉得能够推广到可压缩流体吗？

[yuan_neu]

@yhdthu 新版本的openfoam有加入这个模型吗

[李东岳]

esi最新发布的isoAdvector方法碉堡了

如果你说的是这个，那么在OpenFOAM-1712里面有。