流体Demo92_岳杰顺_高超音速飞行器的气动加热.doc

Copyright Copyright ANSYS教学算例集 高超音速飞行器的气动加热仿真分析 撰写：曾元圆－北京航空航天大学 审核： 校对： TIME \@ yyyy年M月d日 2017年12月18日 目 录 TOC \o 1-3 \h \z \u 1. 摘要 1 2. 返回舱模型介绍 1 3. 前处理 1 4. 求解设置 3 4.1 启动Fluent 3 4.2 模型设置 4 4.3 添加材料 5 4.4 边界条件设置 5 4.5 求解方法设置 7 4.6 流场初始化 9 4.7 求解 9 5. 后处理 11 5.1 启动Results模块 11 5.2 创建参考面 12 5.3 云图显示 12 5.4 矢量图显示 13 5.5 热流计算 13 PAGE 15 PAGE 1 摘要 返回舱在返回大气的过程中与大气层发生剧烈的摩擦，产生大量的热，巨大的动能转化为气体的内能。这种气动热效应甚至高温下的化学反应，会严重影响返回舱内部设备的性能，防热设备占返回舱总重量的20%以上。因此，利用Fluent软件研究高超声速下返回舱的气动热对设计返回舱具有重要意义。 返回舱模型介绍 本案例所模拟的是类Apollo返回舱在高马赫数条件下的气动热。图1是模型示意图，球冠曲率半径为175.26mm，拐角曲率半径为7.30mm，倒锥角为32.50°，最大截面半径为146.05mm，返回舱模型体长为96.04mm，小头直径为44.74mm。 模型示意图（单位：inches[mm]） 前处理 由于三维外形的结构化网格生成比较费时，可采用非结构网格进行计算。在ICEM CFD中建立该三维模型，并划分网格。先生成壳网格；为了捕捉壁面处的气动热，再在壁面处采用棱柱网格进行法向加密；最后自动生成体网格。根据相关文献，为了满足气动热计算网格要求，壁面第一层网格雷诺数量级为10，从而第一层网格间距应在1.0e-5m量级（可根据需求在ICEM中进行调整）。 流场区域采用一个六面体区域，尺寸为长4000mm，宽3000mm，高3000mm。在这个尺寸下流场外边界已经恢复到环境条件，可直接采用压力远场边界条件。 将流场边界分别命名为进口inlet（图中开口处）、出口outlet（黄色部分）、远场farfield（粉红色部分）、舱体壁面（绿色部分）分为五个部分s1、s2、s3、s4、s5，如图2所示。由于s2、s4代表的拐角区域流场变化较为剧烈，故在这两部分对网格进行了加密。图3是返回舱的网格模型。 由于超声速流动在靠近返回舱壁面及尾迹区流场参数变化较大，对这些区域的网格进行加密后可得到较好的数值解，故在壁面及尾迹的一个六面体区域内进行网格加密，如图4所示。 将ICEM CFD划分完成的网格导出，存为capsule_volume.msh文件。 模型命名示意图 网格示意图 加密区示意图 求解设置 启动Fluent 双击Fluent，保持默认设置。将工作目录转换到网格文件所在的目录。 Fluent 启动界面 打开几何文件并选择求解器类型。 File→Read→Case，选择生成的网格capsule_volume.uns。 Grid→Check，检查网格质量，注意Minimum Volume应大于0. Grid→Scale，定义缩放因子。由于在CATIA中造型的时候使用的单位是mm，而尺寸是风洞模型尺寸。我们在fluent中进行全尺寸模拟，故需要对模型尺寸进行修改。在网格导入fluent之后，默认以m作单位。所以此时返回舱长96.04m，直径146.05m，在x、y、z三个方向输入缩放因子0.0267，把尺寸缩小到返回舱最大直径处为3.9m，这与真实的阿波罗返回舱的尺寸是接近的。 在Solver→Type中选择基于密度的（Density-Based）求解器。如图6所示。 网格缩放及选择求解器类型 模型设置 打开模型树中的【Model】→【Viscous】选择Laminar模型，如图7；双击【Energy】打开能量方程开关。有文献表明热流计算表现出对湍流模型不敏感，故这里可以只使用层流模型进行计算。 Model设置 添加材料 双击模型树中的【Materials】, 在Density处选择ideal-gas，在Viscosity处选择sutherland。 材料添加 边界条件设置 来流马赫数取为5，选取40km的大气环境，查标准大气表得在40km高空下的压强为287.13Pa，温度为250.35K。壁面温度边界条件取为300K等温壁，速度无滑移条件。 编辑【Boundary Condition】中的farfield、inlet、outlet，这三个边界采用相