计算流体力学试编.docx

计算流体力学大作业系 （院）： 宇航学院 专 业： 航空宇航推进理论与工程 2016年6月1.题目要求大作业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入口、总压0.2MPa图11压力云图图12温度云图图13马赫数云图图14流线图4.工况点4：粗网格、基于密度、BC入口、总压6MPa图15压力云图图16温度云图图17马赫数云图图18流线图5.工况点5：粗网格、基于压力、AB入口、总压0.2MPa图19压力云图图20温度云图图21马赫数云图图22流线图6.工况点6：粗网格、基于压力、AB入口、总压6MPa图23压力云图图24温度云图图25马赫数云图图26流线图7.工况点7：粗网格、基于压力、BC入口、总压0.2MPa图27压力云图图28温度云图图29马赫数云图图30流线图8.工况点8：粗网格、基于压力、BC入口、总压6MPa图31压力云图图32温度云图图33马赫数云图图34流线图9.工况点9：细网格、基于密度、AB入口、总压0.2MPa图35压力云图图36温度云图图37马赫数云图图38流线图10.工况点10：细网格、基于密度、AB入口、总压6MPa图39压力云图图40温度云图图41马赫数云图图42流线图11.工况点11：细网格、基于密度、BC入口、总压0.2MPa图43压力云图图44温度云图图45马赫数云图图46流线图12.工况点12：细网格、基于密度、BC入口、总压6MPa图47压力云图图48温度云图图49马赫数云图图50流线图13.工况点13：细网格、基于压力、AB入口、总压0.2MPa图51压力云图图52温度云图图53马赫数云图图54流线图14.工况点14：细网格、基于压力、AB入口、总压6MPa图55压力云图图56