gambit里面网格质量检查方法及fluent计算udf的导入步骤.doc

第一部分：gambit中网格质量的检查方法如何检查网格质量，用什么指标来说明网格好不好呢？怎么控制？一般是什么原因造成的? 一般也就是，网格的角度，网格变形的梯度等等吧 判断网格质量的方面有很多，不知你用的是什么软件，下面总结的是针对Gambit帮助文件的简单归纳，不同的软件有不同的评价单元质量的指标，使用时最好仔细阅读帮助文件。Area单元面积，适用于2D单元，较为基本的单元质量特征。 Aspect Ratio长宽比，不同的网格单元有不同的计算方法，等于1是最好的单元，如正三角形，正四边形，正四面体，正六面体等；一般情况下不要超过5：1. Diagonal Ratio对角线之比，仅适用于四边形和六面体单元，默认是大于或等于1的，该值越高，说明单元越不规则，最好等于1，也就是正四边形或正六面体。 Edge Ratio长边与最短边长度之比，大于或等于1，最好等于1，解释同上。 EquiAngle Skew通过单元夹角计算的歪斜度，在0到1之间，0为质量最好，1为质量最差。最好是要控制在0到0.4之间。 EquiSize Skew通过单元大小计算的歪斜度，在0到1之间，0为质量最好，1为质量最差。2D质量好的单元该值最好在0.1以内，3D单元在0.4以内。 MidAngle Skew通过单元边中点连线夹角计算的歪斜度，仅适用于四边形和六面体单元，在0到1之间，0为质量最好，1为质量最差。 Size Change相邻单元大小之比，仅适用于3D单元，最好控制在2以内。 Stretch伸展度。通过单元的对角线长度与边长计算出来的，仅适用于四边形和六面体单元，在0到1之间，0为质量最好，1为质量最差。 Taper锥度。仅适用于四边形和六面体单元，在0到1之间，0为质量最好，1为质量最差。 Volume单元体积，仅适用于3D单元，划分网格时应避免出现负体积。 Warpage翘曲。仅适用于四边形和六面体单元，在0到1之间，0为质量最好，1为质量最差。 另外，在Fluent中的窗口键入：grid quality 然后回车，Fluent能检查网格的质量，主要有以下三个指标：1.Maxium cell squish: 如果该值等于1，表示得到了很坏的单元；2.Maxium cell skewness: 该值在0到1之间，0表示最好，1表示最坏；3.Maxium aspect-ratio: 1表示最好。第二部分：fluent中导入UDF进行计算的基本步骤：1. 导入mesh，Filereadcase2. Grid-Scale 调整x， y的系数，达到需要尺寸3. 因为是可压缩流动，DefineModel-Solver选density based，其他保持不变4. 当考虑温度时，需要选中energy equation DefineModelEnergy5. Materials选择需要的气体 DefineMaterial，比如这里用氮气，点击copy6. 选择好氮气之后，在密度一栏设置成理想气体。然后点击Changeclose7. Define-Operating conditions里Operating pressure设为零8. 导入UDF：Defineuser Defined-functions- compiled-，在跳出来的对话框里，点击Add，选择自己编写的udf，如下图。然后点击Build，主窗口会出现以下信息：再在compiled- UDF 页面点击Load，主窗口有以下信息：表明编译通过，已经可以使用了9. DefineUser Defined-function hooks界面，点击Edit，选择自己编的UDF里的adjust程序，如下图：10. DefineUser DefinedMemory界面，设定UDMI个数，根据自己的程序。11. 设置滑移边界条件：DefineBoundary Conditions 选择壁面，然后点击Set，进行设置： 壁面类型：Momentum 页面，设置滑移速度，从下拉菜单中找到对应的UDF，如下图Thermal 页面，设置温度跃变，从下拉菜单中找到对应的UDF，如下图：12. 设置源项Boundary：DefineBoundary Conditions 选择fluids 然后点击Set进入先选择气体种类，然后在source Terms前打勾，然后点击X Momentum右边的Edit，出现如下界面：把0变为1，然后从下拉菜单中找到合适的UDF另外两个一样处理13. 设置入口条件：DefineBoundary Conditions 选择inlet 然后点击Set进入设置总压和表压，如果要设置温度，点击Thermal页面，进行设置。14. 设置出口条件：DefineBoundary Conditions 选择outlet 然后点击Set进入根据出口压强进行设置15. 初始化流场: SolveInitializeinitionalize:根据经验，稍微调整，速度不要太大，温度不要太低。16. 设置残差收敛准则： SolveMonitorsResidual可以在左上角的Plot前打勾，这样可以实时看到残差曲线，同时设置Residual，一般设成10E-617. 设置计算时的Courant 数SolveControlsSolutionCourant数一开始小一点，有利于收敛，然后