计算流体力学报告.docx

计算流体力学大作业 系 （院）： 宇航学院专 业： 航空宇航推进理论与工程 2016年6月目录1.题目要求12.核心参数选取23.网格划分34.参数设置35.计算结果41.工况点1：粗网格、基于密度、AB入口、总压0.2MPa42.工况点2：粗网格、基于密度、AB入口、总压6MPa53.工况点3：粗网格、基于密度、BC入口、总压0.2MPa74.工况点4：粗网格、基于密度、BC入口、总压6MPa85.工况点5：粗网格、基于压力、AB入口、总压0.2MPa106.工况点6：粗网格、基于压力、AB入口、总压6MPa117.工况点7：粗网格、基于压力、BC入口、总压0.2MPa138.工况点8：粗网格、基于压力、BC入口、总压6MPa149.工况点9：细网格、基于密度、AB入口、总压0.2MPa1610.工况点10：细网格、基于密度、AB入口、总压6MPa1711.工况点11：细网格、基于密度、BC入口、总压0.2MPa1912.工况点12：细网格、基于密度、BC入口、总压6MPa2013.工况点13：细网格、基于压力、AB入口、总压0.2MPa2214.工况点14：细网格、基于压力、AB入口、总压6MPa2315.工况点15：细网格、基于压力、BC入口、总压0.2MPa2516.工况点16：细网格、基于压力、BC入口、总压6MPa2617.计算结果分析2718.工况点18：粗网格、基于压力、AB入口、总压6MPa、二阶精度2819.推力计算301.工况点1：细网格、基于压力、AB入口、总压0.2MPa302.工况点2：细网格、基于压力、AB入口、总压6MPa301.题目要求大作业2(3专业) 几何参数：关键参数为喷管喉部半径Rt(单位为MM)。该参数由每个同学的学号最后3位数字确定。设最后3位数字为abc，则：R = 10*max(a，b，c) + min(a，b，c)喷管入口半径Rc = 3*Rt，喷管出口半径Re = 3*Rt，燃烧室长度L= 10*Rt，轴对称坐标系。流体属性：空气(理想气体、可压缩)，只计算层流。设置工作压力为0以免混淆绝对压力与相对压力。1.AHGF为对称轴，EF为压力出口边界，背压 ，BCDE为绝热壁面，AB为压力入口边界，来流总温3500K，总压分别取0.2MPa和6MPa；初始计算时可设置压力入口边界静压为其总压的90%；2.变更AB为壁面边界，BC为压力入口边界，来流总温3500K，总压分别取0.2MPa和6MPa。3.分别采用粗网格和细网格进行计算，建议粗、细网格下轴向网格点数分别为100、200，径向网格点数分别为30、60。4.分别采用基于密度的算法和基于压力的算法。输出结果：1.要求在Tecplot中显示温度场、压力场、马赫数分布、流线图；2.比较粗网格和细网格下流场的差别。3.比较采用基于密度和基于压力的算法的收敛情况。4.上述16种工况只要求一阶精度，但选取AB为压力入口边界，总压6MPa，粗网格条件下做基于压力算法的二阶精度计算，并将其与细网格结果比较；5.选取AB为压力入口边界、细网格、基于压力算法的2种工况计算该发动机的推力，推力可由ABCDEF边界上的压力对轴向面积进行积分减去EF边界的轴向面积与背压 的乘积而得到。6.当入口边界总压为0.2MPa时，出口EF边界有可能同时出现入流和出流，导致计算无法完全收敛，无需强求收敛。2.核心参数选取由于本人学号是SY1515333，学号后三位是333，所以核心参数：Rt= 10*max(a，b，c) + min(a，b，c)=33mm喷管入口半径Rc = 3*Rt=99mm喷管出口半径Re = 3*Rt=99mm燃烧室长度L= 10*Rt=330mm3.网格划分题目要求网格分为粗网格和细网格，所以本文给出了两个网格（将入口AB和BC合为一个网格），见图1、2。图 1粗网格图 2细网格其中应节点设置要求，粗网格网格数为3073，细网格网格数为12143，相较下变化明显。4.参数设置由于题目要求分别求解基于密度和基于压力两种情况下的解，所以算法会有不同，见表1。表格 1参数设置计算模型层流模型计算模型选择能量方程和Laminar模型求解方法选择Implicit（基于密度）SIMPLE（基于压力）精度一阶精度（一例采用二阶进度）入口边界条件压力入口6MPa或0.2MPa、温度（3500K）出口边界条件压力出口0.1MPa壁面边界条件无滑移壁面流体介质空气（理想气体）5.计算结果1.工况点1：粗网格、基于密度、AB入口、总压0.2MPa图 3压力云图图 4温度云图图 5 马赫数云图图 6流线图2.工况点2：粗网格、基于密度、AB入口、总压6MPa图 7压力云图图 8温度云图图 9马赫数云图图 10流线图3.工况点3：粗网格、基于密度、BC