单元生死技术

[昝磊]

请问OpenFOAM中有类似ANSYS中的单元生死技术的模块吗？

[东岳]

有，但是用的很少，基本只是个demo，不能用来计算，除非你只需要个框架来进行代码开发。如果你急需类似的网格处理方式进行真实的模拟，不如用OpenFOAM v1712中的overset。

[昝磊]

@李东岳 谢谢您的回复! 请问您说的overset可以实现下面的过程吗?

图中有两套网格，红色网格是背景网格，黑色方框也是网格。通过将黑色方框网格运动，使背景网格在黑色方框外的部分不断增加，并作为求解域（即只有黑色方框外的红色网格才是求解域) . 并且希望新增红色网格的物理量(如密度)是事先指定的,而不是由上一时间步插值计算出来的.
请问这个过程overset可以实现吗? 如果不能, 请问有什么办法可以实现呢?
谢谢!

[东岳]

你这个确实用生死单元更好处理，类似下面的这种东西吧？要是几何就是这么简单的话可以搞。我建议你看看这部分的代码 https://github.com/OpenFOAM/OpenFOAM-4.x/tree/master/src/dynamicMesh/polyTopoChange

下面这个动画是我几年前搞的，现在有点忘了，但是不是很难。主要是用的人太少。你是我遇到的第二个。

[昝磊]

@李东岳 李老师您好！您给的代码链接就是生死单元的demo吗？
谢谢！

[东岳]

是的，我几年前做的。

[昝磊]

@李东岳 谢谢您，帮了我大忙了！😀

[我是河滩]

@昝磊 你好，生死单元你学会了吗？如何学习Openfoam中动网格的相关代码？感觉文件很多，不知从何看起？

[史浩]

fvMatrix::setValues 和 fvMatrix::setValuesFromList这两个函数将对应单元设置为固定值，和生死单元很像。动网格也可以，但是动网格会牵扯到通量守恒检查，并行的时候也容易出错，不过效果也不错。

[昝磊]

@史浩 嗯，现在磕磕绊绊是做完了，不过计算效率是真的低 我目前也没时间再细化了，以后再研究吧。

[昝磊]

@史浩 没有搞这个生死单元，还是用了动态铺层的网格。

[东岳]

@昝磊 大作发表了发上来让大家学习下 用OpenFOAM算的么

[昝磊]

@东岳 哈哈，是OpenFOAM.