早上好

我正在尝试开发一个非结构化 CFD 代码来解决有限体积（以单元为中心）上下文中的欧拉方程（学习目的）。我能够从 cgns 文件构建代码的一些基本信息。我现在拥有的是：

• 节点索引及其 'x' 'y' 坐标。
• 元素索引和构建这些元素的节点（索引）。
• 边界条件。

由于代码是二维的，因此每个元素实际上都是一个单元格，每个面实际上是形成该单元格的一条线。

我现在构建的基本数据结构是： + 节点数组，我在其中放置节点索引并在 x 和 y 坐标之间进行选择。+ 元素/节点连接（逆时针方式）其中（对于四边形）：

connec(1,nelem) = node 1
connec(2,nelem) = node 2
connec(3,nelem) = node 3
connec(4,nelem) = node 4

这些数据结构构建起来非常简单，但下一个对我来说非常模糊。我需要一个基于人脸的数据结构，其中每个人脸都“知道”形成左元素索引和右元素索引的节点。

最终结果将是：

faces(1,face_index) = node 1 of the face.
faces(2,face_index) = node 2 of the face.
faces(3,face_index) = left element.
faces(4,face_index) = right element.

我进行了很多搜索，但找不到一个参考资料来解释我构建这个面部数据结构的算法。Lohner 关于有限元的书对基于节点的数据结构有很好的解释，但对基于面的数据结构却没有。

任何人都有很好的参考来解释这些算法甚至算法本身？

谢谢 ！干杯