单元生死技术

[李东岳]

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

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

[昝磊]

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

[李东岳]

是的，我几年前做的。

[昝磊]

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

[我是河滩]

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

[史浩]

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

[昝磊]

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

[昝磊]

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

[李东岳]

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

[昝磊]

@东岳 哈哈，是OpenFOAM.