商用软件GAMBIT网格数据解析和应用.pdfVIP

商用软件GAMBIT 网格数据解析和应用 1, 2 1 2 刘晶 ，张敏 ，John C. Chai 1 南京理工大学动力工程学院，南京(210094) 2 南洋理工大学机械与航天学院，新加坡(639798) 摘 要：利用商业软件GAMBIT 非结构化网格生成技术，解决非FLUENT 用户使用不同解 算器时，网格生成耗费大量机时的问题。通过转换GAMBIT 生成的网格数据结构，连接编 制求解导热方程数值计算程序，并且进行数值计算和与精确解比较，以此证明本方法的合理 性。 关键词: GAMBIT ，非结构网格，热传导 引言 网格生成是数值计算的第一步。在工程数值模拟中，数值计算结果的最终精度及计算过 程效率，主要取决于所生成的网格与所采用的算法。一般而论，所有的各种生成网格方法在 一定的条件下都有其优越性及弱点，各种求解算法也各有适应范围。一个成功而高效的数值 [1] 计算，只有在网格的生成及求解算法这两者之间有良好的匹配时才能实现 。 数值计算中采用的网格可以大致分为结构化网格与非结构化网格两类。结构化网格的特 点是每个节点与其邻点之间的连接关系固定不变，且隐含在所生成的网格中，不必专门设置 数据去确认节点与邻点之间的这种关系。而非结构化网格单元与节点的编号就无固定规则可 遵循，它消除了结构网格中节点的结构性限制，节点和单元的分布可控性好，因而能较好地 处理边界，适用于模拟真实复杂外型。但它的缺点是网格生成过程比较费时，离散时所形成 代数方程求解的收敛特性差。 为了克服非结构化网格生成复杂的缺点，本文对商用软件GAMBIT 的非结构化网格数 据进行解析应用，在非结构化网格上用单元中心有限体积法[2,3]进行热传导数值计算。全文 共分三个部分，第一部分是用C/C++编程语言读取商用软件GAMBIT 网格生成的数据文件， 提取有效信息进行数据的转换和重组；第二部分用 FORTRAN 编写解算器，来验证所得到 网格的正确性；最后给出两个热传导算例，进行数值解与精确解的比较，并对计算结果进行 分析讨论。 1. GAMBIT 数据结构的应用 一些商用化的网格生成软件，如ICEM_CFD 、TGrid、GAMBIT 等是面向CFD 专业的 前处理器软件[4]，它们具有强大的网格划分功能，可以划分出包含边界层的 CFD 特殊要求 的高质量网格。GAMBIT 网格划分方法的选择完全是智能化，当你选择一个几何区域后， GAMBIT 会自动选择最合适的网格划分算法，使网格划分过程变的极为容易。GAMBIT 中 的网格划分算法可以保证在较为复杂的几何区域内直接划分出高质量的六面体网格。 GAMBIT 中可使用的求解器（SOLVER）有FLUENT5/6,FLUENT4,GENERTIC 等，它 们输出的文件类型有：*.FDNEUT，*.MSH，*.NEU，*.CDB，*.GRD。在这些文件中存储的 是所生成非结构化网格的点、面、体信息，不同类型的文件给出的信息各异。本研究利用 *.NEU 和*.MSH 这两种类型的文件来再生成网格。*.NEU 中给出了点（NUMNP ）和控制单 元（NELEM ）总数，以及网格点和单元的信息；*.MSH 采用的是十六进制，记录的是网格 点、面和单元关联信息，包括边界条件。 -1- 用C/C++编程语言进行信息的转换和重组时，要重新生成点（X （I ），Y （I ），Z （I ））、 面（NFTYPE （I ）、NFNODE(J,I) ）和体（NCTYPE （I ）、NCNODE （J,I ））的信息。对于 GAMBIT，计算区域边界上，面的四个点的序号服从左手螺旋定则