OpenFOAM根据自带求解器计算自己算例-生成网格-边界条件总结.pdfVIP

采用OpenFOAM 自带求解器计算算例 ——网格生成边界条件 自己的公众号：一介UDFer，不定期更新个人Fluent+UDF 和 OpenFOAM 学习经验，交流学习。 说明：此处的应用只限于使用OpenFOAM 自带的求解器计算自己想计算的 case，对求解的离散格式，求解方法，湍流模型，燃烧模型等均未做修改。 1 OpenFOAM 的case 结构 OpenFOAM （OF ）和fluent 有很大的区别，OF 完全是靠代码设置，fluent 有 很好的用户界面，通过下拉菜单选择。OF 自带了很多求解器，这些求解器适用 于固定的工况，比如说物理模型的几何结构及边界条件，方程的离散方式及求解 方式等。以XiFoam （OF 中的预混燃烧求解器）为例，OpenFOAM 计算的case 构架如下： 1.1 System 文件夹： controlDict 文件：对时间步、输入输出时间、场数据读取和写入进行控制 fvSchemes 文件：指定有限体积法的离散格式 fvSolution 文件：指定方程组矩阵求解器、残差以及其他算法控制 1.2 constant 文件夹： 存在网格文件（ployMesh 文件夹下） 湍流模型选择字典文件、流动特性参数、流体热物理性参数，化学反应相关 参数，重力等 1.3 0 文件夹： 存放初始条件字典文件（速度、压力、温度，其他与流动或者燃烧相关的参 数初始值） 2 生成网格 如果使用OF 自带的求解器计算自己想要计算的case，第一步是生成相应的 网格文件。 采用代码画网格，费时间，并且对于几何结构比较复杂的几何结构也不适合， 所以可以采用专业的网格生成软件画网格，然后转化为OpenFOAM 网格文件。 将自己画好的网格（此处指 *.mesh 文件），采用命令fluentMeshToFoam 转化成 OF 网格。命令详细介绍参考/index.php/FluentMeshToFoam ； 这里我常用的是 fluentMeshToFoam + 空格+ 网格文件名字+ 空格+-scale+ 空格 +0.001 。如下：fluentMeshToFoam 1mm.mesh –scale 0.001 ；意思是将 1mm.mesh 这个网格文件转化成 OF 网格，并将网格尺寸改为mm （我画网格时 通常采用m 为单位）。 3 设置边界条件 3.1 检查边界类型 检查转化后的边界条件是否是自己所需要的边界条件（在 case 文件下下的 constant/ polymesh/ boundary 文件）。 OpenFOAM 中基本的边界类型如下： 这几类是定义在polymesh 文件夹下boundary 文件中的，比如我们定义的入 口和出口等，type 类型为patch ，这一部分通常是由fluentMeshToFoam 自动生成 的，我们做计算之前需要做一下检查。声明：需要注意的是，这些种类只是声明 边界类型的，和边界条件还有很大的差异，边界条件是要定义具体的值（固定的 值、固定的梯度或者其他）。 常用的有： Patch: 出口和入口类型采用该类型。 Wall ：各种壁面壁面。 Empty:用在2D 计算，OpenFOAM 没有2D 和3D 之分，计算2D 是将第三 维度设成empty 。 实际代码中如下： 代码中outlet，inlet，lowerWall 是画网格是给边界起的名字，type 定义的是 边界的类型，type 的下两行是该边界的区域。 3.2 修改0 文件夹下的初始各个文件下的边界条件 说明：此处的边界条件和boundary 文件下的边界不一样，boundary 文件下 的边界类型只是定义了边界的类型，并没有对边界处的参数值进行定义，边界处 的参数值定义均在0 文件夹下各个相应的文件中定义。 将 0 文件夹下的所有自带的初始条件文件中的边界类型调整成和 boundary 中一样的类型，并且赋予其相应的边界值。 3 问题 这里边最大的问题是OpenFOAM 的边界条件都是怎么定义的，相应的代码 是什么？需要赋予什么样的值？比如说，速度入口处的速