OpenFOAM怎么保证Sf是由owner指向neighbour的？

[Voynich]

想用python编一个简单的非结构网格simple算法加深理解，发现openfoam网格的Sf计算出来后不需要判断方向自然就是由owner指向neighbour的，这是怎么保证的？其他网格软件也有类似的feature吗？

[Cp_Zhao]

盲猜： 构建面是按照单元号从小到大遍历的，面方向为指向单元外侧。故面肯定是由小编号单元先建立的，指向小编号单元外侧（即从owner向neighbour）

[李东岳]

想用python编一个简单的非结构网格simple算法加深理解

大佬这玩意整的高端，相当屌。你这个搞出来，一个FVM求解器的雏形就出来了。搞出来之后开源么？这个相当流弊了。在这个基础上做速度压力耦合啥的，一个CFD软件后端雏形就出来了。

发现openfoam网格的Sf计算出来后不需要判断方向自然就是由owner指向neighbour的，这是怎么保证的？

对于非三角形，openfoam里面面的顶点编号首先要保证按照顺序来，不能从对角线来。在保证顺序转圈的基础上，这些点的计算面$\bfS_f$的公式默认调用的是右手定则。因此在给定你每个面的按照顺序转圈的顶点编号后，其计算出来的面矢量，就是右手定则的指向。因为你的$\bfS_f$已经算出来了，那么这个指向就出来了。每个面都是从小网格编号指向大网格编号。小网格就是owner，大编号的网格就是nei。

其他软件不知道。我也好奇。

[Voynich]

@Cp_Zhao 感谢回复

[Voynich]

@李东岳 感谢回复。我是想先写个简单的特定算例学习一下算法，格式边界条件啥的都直接写死，然后再慢慢做的通用一点，能搞出来肯定开源。