毕业论文.2.doc

沈阳化工大学学士学位论文论文题目：张小岑热能090109110112内容摘要第一章 文献综述1.1课题背景螺旋管因其具有传热效率高、加工制作方便等优点，在高效蒸汽发生器和冷却器、核反应堆、电站锅炉、船舶动力、石油化工、航天航空、微电子器件冷却、先进燃料电池系统冷却、食品制药、以及制冷与低温技术等领域得到了非常广泛的应用螺旋流道作为一种强化换热通道，由于其结构紧凑，广泛应用于动力、石油和化工等换热设备中fl.流体在螺旋管内流动时，由于离心力的作用，会在垂直于主流方向的横截面上产生二次流，二次流相对于轴向主流速度在数量级上虽然较小，但其与轴向主流复合成螺旋式的前进运动，破坏了边界层，强化了换热，因此研究螺旋流道内流场结构及换热特性对工程实际应用具有重要的意义。1.2螺旋管的研究现状 到目前为止，国内外已对管内两相流动与传热进行了长期的研究探索，但大部分研究是以水一空气或水一水蒸汽为工质对直管进行的，许多关于直管内的两相流流型、空隙率、流动阻力和传热特性的数据和计算关联式己被业内人士广泛认可和应用。我国西安交通大学在这方面进行了长期系统的研究，取得了显著的成果。而对矩形螺旋管内的两相流动与传热的研究则非常少，现有的一些文献主要是关于矩形螺旋管内的凝结传热特性的研究，少数的关于矩形螺旋管内流动沸腾特性的研究，而关于矩形卧式螺旋管内两相流极为少见。 1974年，美国加利福尼亚大学的Rowland教授和他的博士后Molina13在“自然”杂志上发表文章，指出正在广泛使用的CFC类制冷工质会消耗地球周围的臭氧层，对地球生物和人类带来危害。之后，对螺旋管和其它环保制冷工质的两相流动与传热特性的研究，以及以螺旋管为实验工质进行的模化研究得到了国内外广泛的重视，也已取得了不少的成果，大多数也是针对直管进行的。从上世纪初以来，很多学者对螺旋管内的单相流动与传热进行了实验研究，也有一些关于螺旋管和其它曲线管内单相流动与传热的数值模拟研究的文献报道。但对于螺旋管内的两相流动与传热特性的研究相对较少，并且其中多数是针对空气一水或水一水蒸汽进行的。圆形螺旋风管系采用宽为1 Ocm的被锌俐带绕制成，有单层与两层风管中夹保温材料构成的双层保温风管两种型式.圆形螺旋风管的摩擦阻力的大小决定于咬口接缝制作模具.实验曲线表明:由于各种口径风管的制作模具不同，所以摩擦阻力曲线排列无理论规律，必须对其各别测定.另外又rIE明，Ii 115螺淀风管覃漆阻灯虽咯大千普通镀锌ill板摩阻，但很接近.目前此风管月多于船用空调系统。众所周知，在通风空调工程中应用的通风管道内，其空气流动状态，绝大多数是在过渡区域内，风速为5 -40m/s范围内.因此当雷诺数Re与风管直径D为定值时，入值主要取决于风管表面的粗糙度K值.由于圆形螺旋风管是采用宽为l0cm的镀锌钢带制成，整个风管具有螺距为l0cm的咬口接缝，所以粗糙度K值不能简单采用材料表面的绝对粗糙度，应考虑咬口接缝的影响因素，而应采用当量水力粗糙度表示。因此咬口接缝加工的质量，直接影响到当量水力粗糙度数值的大小，亦即影晌螺旋风管的摩擦阻力系数.由于各生产厂家的加工水平与选用材料不同，所以要对各厂生产的螺旋风管进行摩擦阻力的实验研究，取得其摩擦阻力系数，制作计算图表，以供用户使用。由于流场特性主要取决于流道结构和操作参数，因此，由于截面形状不同，三角形流道内流体的流动特性肯定有别于其他截面形状的螺旋流道.同时在实际的工程应用中，也存在三角形螺旋流道的典型应用，如角钢夹套，它是将角钢螺旋缠绕焊接在筒体状容器的外侧，角钢与筒体外壁构成一个三角形截面螺旋流道，流道内通入加热或冷却介质，可加热或冷却容器内的物料。 与半圆形螺旋流道流体流动的阻力相比，当Re数相同且半圆形螺旋流道螺距、曲率半径及当量直径按模型1取值时，三角形螺旋流道的阻力较小，流体流动性能好，且随Re增大，三角形螺旋流道的优势愈加明显.在所研究的Re范围内，三角形螺旋流道的阻力系数厂为半圆管螺旋流道阻力系数的84.1%-99.5%.螺旋管由于其结构紧凑和高换热性能等特点，在发电机组、核电设备、制冷和食品加工等工业场所被作为热交换器的基本结构而广泛使用。与直管相比，螺旋管具有较高的对流换热效率。人们对圆形横截面螺旋管已经进行了一系列的理论、数值及实验研究。采用数值分析方法，研究了扭矩对矩形截面螺旋管内充分发展态层流流动的影响。结果表明，螺旋管的曲率及扭矩是影响层流换热性能的主要参数。在航空航天领域，如航空发动机、火箭发力机以及冲压发动机的冷却系统中，矩形横截面内螺旋管冷却结构由于加工技术相对简单，已被广泛使用。目前，对于矩形横截面螺旋管的研究仗部分局限于层流流动，对于湍流的研究非常有限。1.3本课题的研究方向学会使用Cambit和Fluent软件，使用第二章Gambit和Fluent软件的介绍由于我需要学习FLUENT来做毕业设计，王老师给了我一本书，王瑞金的FLUENT技术基础与应用实例，当然，学这本书之前必须要有两个条件：第一，具有流体力学的基础，第二，有FLUENT安装软件可以应用。然后就照着书上二维的计算例子，一个例子，一个步骤地去学习，然后学习三维，再针对具体你所遇到的项目进行针对性的计算。不能急于求成，从前理器GAMBIT，到通过FLUENT进行仿真，再到后处理，如TECPLOT，进行循序渐进的学习，坚持，效果是非常显著的。如果身边有懂得FLUENT的老师，那么遇到问题向老师请教是最有效的方法，碰到不懂的问题也可以上网或者查找相关书籍来得到答案。2.1Gambit软件的介绍 GAMBIT软件是面向CFD的专业前处理器软件，它包含全面的几何建模能力,既可以在GAMBIT内直接建立点、线、面、体几何，也可以从主流的CAD/CAE系统如PRO/E、UGII、IDEAS、CATIA、SOLIDWORKS、ANSYS、PATRAN导入几何和网格，GAMBIT强大的布尔运算能力为建立复杂的几何模型提供的极大的方便。GAMBIT具有灵活方便的几何修正功能，当从接口中导入几何时会自动的合并重合的点、线、面；GAMBIT在保证原始几何精度的基础上通过虚拟几何自动的缝合小缝隙，这样既可以保证几何精度，又可以满足网格划分的需要。GAMBIT功能强大的网格划分工具，可以划分出包含边界层等CFD特殊要求的高质量的网格。GAMBIT中专有的网格划分算法可以保证在较为复杂的几何区域可以直接划分出高质量的六面体网格。GAMBIT中的TGRID方法可以在极其复杂的几何区域中可以划分出与相邻区域网格连续的完全非结构化的网格，GAMBIT网格划分方法的选择完全是智能化的，当你选择一个几何区域后GAMBIT会自动选择最合适的网格划分算法，是网格划分过程变的极为容易。GAMBIT可以生成FLUENT5、FLUENT4.5、FIDAP、POLYFLOW、NEKTON、ANSYS等求解器所需要的网格。2.1.1Gambit软件的特点 ACIS内核基础上的全面三维几何建模能力，通过多种方式直接建立点、线、面、体，而且具有强大的布尔运算能力，ACIS内核已提高为ACIS R12。该功能大大领先于其它CAE软件的前处理器； 可对自动生成的Journal文件进行编辑，以自动控制修改或生成新几何与网格； 可以导入PRO/E、UG、CATIA、SOLIDWORKS、ANSYS、PATRAN等大多数CAD/CAE软件所建立的几何和网格。导入过程新增自动公差修补几何功能，以保证GAMBIT与CAD软件接口的稳定性和保真性，使得几何质量高，并大大减轻工程师的工作量； 新增PRO/E、CATIA等直接接口， 使得导入过程更加直接和方便； 强大的几何修正功能，在导入几何时会自动合并重合的点、线、面；新增几何修正工具条，在消除短边、缝合缺口、修补尖角、去除小面、去除单独辅助线和修补倒角时更加快速、自动、灵活，而且准确保证几何体的精度； G/TURBO模块可以准确而高效的生成旋转机械中的各种风扇以及转子、定子等的几何模型和计算网格； 强大的网格划分能力，可以划分包括边界层等CFD特殊要求的高质量网格。GAMBIT中专用的网格划分算法可以保证在复杂的几何区域内直接划分出高质量的四面体、六面体网格或混合网格； 先进的六面体核心(HEXCORE)技术是GAMBIT所独有的，集成了笛卡尔网格和非结构网格的优点，使用该技术划分网格时更加容易，而且大大节省网格数量、提高网格质量； 居于行业领先地位的尺寸函数（Size function）功能可使用户能自主控制网格的生成过程以及在空间上的分布规律，使得网格的过渡与分布更加合理，最大限度地满足CFD分析的需要； GAMBIT可高度智能化地选择网格划分方法，可对极其复杂的几何区域划分出与相邻区域网格连续的完全非结构化的混合网格； 新版本中增加了新的附面层网生成器，可以方便地生成高质量的附面层网格； 可为FLUENT、POLYFLOW、 FIDAP、ANSYS等解算器生成和导出所需要的网格和格式。2.1.2 Gambit的操作步骤步骤1：启动Gambit软件步骤2：文件的创建及求解器的选择步骤3：创建控制点步骤4：创建边步骤5：创建面步骤6：体的创建步骤7：实体网格的划分步骤8：边界条件类型的指定步骤9：Mesh文件的输出2.2 Fluent软件的介绍Fluent是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%，凡是和流体、热传递和化学反应等有关的工业均可使用。它具有丰富的物理模型、先进的数值方法和强大的前后处理功能。CFD商业软件FLUENT，是通用CFD软件包，用来模拟从不可压缩到高度可压缩范围内的复杂流动。由于采用了多种求解方法和多重网格加速收敛技术，因而FLUENT能达到最佳的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网格技术及成熟的物理模型，使FLUENT在转换与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面有广泛应用。目前与FLUENT配合最好的标准网格软件是ICEM，而不是早已过时的GAMBIT。2.2.1 Fluent FLUENT软件特点(1) Fluent 软件采用基于完全非结构化网格的有限体积法，而且具有基于网格节点和网格单元的梯度算法；(2)定常/非定常流动模拟，而且新增快速非定常模拟功能； (3) Fluent T软件中的动/变形网格技术主要解决边界运动的问题，用户只需指定初始网格和运动壁面的边界条件，余下的网格变化完全由解算器自动生成。网格变形方式有三种：弹簧压缩式、动态铺层式以及局部网格重生式。其局部网格重生式是Fluent 所独有的，而且用途广泛，可用于非结构网格、变形较大问题以及物体运动规律事先不知道而完全由流动所产生的力所决定的问题；(4) Fluent 软件具有强大的网格支持能力，支持界面不连续的网格、混合网格、动/变形网格以及滑动网格等。值得强调的是，FLUENT软件还拥有多种基于解的网格的自适应、动态自适应技术以及动网格与网格动态自适应相结合的技术；(5) Fluent 软件包含三种算法：非耦合隐式算法、耦合显式算法、耦合隐式算法，是商用软件中最多的；(6) Fluent 软件包含丰富而先进的物理模型，使得用户能够精确地模拟无粘流、层流、湍流。湍流模型包含Spalart-Allmaras模型、k-模型组、k-模型组、雷诺应力模型(RSM)组、大涡模拟模型(LES)组以及最新的分离涡模拟(DES)和V2F模型等。另外用户还可以定制或添加自己的湍流模型；(7)适用于牛顿流体、非牛顿流体；含有强制/自然/混合对流的热传导，固体/流体的热传导、辐射；化学组份的混合/反应；(8)自由表面流模型，欧拉多相流模型，混合多相流模型，颗粒相模型，空穴两相流模型，湿蒸汽模型；融化溶化/凝固；蒸发/冷凝相变模型；风扇，散热器，以热交换器为对象的集中参数模型(9)离散相的拉格朗日跟踪计算；非均质渗透性、惯性阻抗、固体热传导，多孔介质模型（考虑多孔介质压力突变）；惯性或非惯性坐标系，复数基准坐标系及滑移网格；动静翼相互作用模型化后的接续界面；基于精细流场解算的预测流体噪声的声学模型；质量、动量、热、化学组份的体积源项；丰富的物性参数的数据库；(10)磁流体模块主要模拟电磁场和导电流体之间的相互作用问题；连续纤维模块主要模拟纤维和气体流动之间的动量、质量以及热的交换问题；(11)高效率的并行计算功能，提供多种自动/手动分区算法；内置MPI并行机制大幅度提高并行效率。另外，FLUENT特有动态负载平衡功能，确保全局高效并行计算；(12) Fluent 软件提供了友好的用户界面，并为用户提供了二次开发接口（UDF）；(13) Fluent FLUENT软件采用C/C+语言编写，从而大大提高了对计算机内存的利用率。(14)在CFD软件中，Fluent软件是目前国内外使用最多、最流行的商业软件之一。Fluent的软件设计基于CFD计算机软件群的概念,针对每一种流动的物理问题的特点，采用适合于它的数值解法在计算速度、稳定性和精度等各方面达到最佳。由于囊括了Fluent Dynamical International比利时PolyFlow和Fluent Dynamical International(FDI)的全部技术力量(前者是公认的在黏弹性和聚合物流动模拟方面占领先地位的公司,后者是基于有限元方法CFD软件方面领先的公司 ),因此Fluent具有以上软件的特点。2.2.2 Fluent的求解器求解步骤1：Fluent求解器的选择步骤2：网格的相关操作步骤3：选择计算模型步骤4：定义流体的物理性质步骤5：设置边界条件步骤6：求解方法的设置及其控制步骤7：Tecplot后处理第三章 Fluent实例简介一、模拟实验目的和内容本模拟实验的目的主要有3个：（1）学生初步了解并掌握Fluent求解问题的一般过程，主要包括前处理、计算、后处理三个部分。（2）理解计算机求解问题的原理，即通过对系统进行离散化，从而求解代数方程组，求得整个系统区域的场分布。（3）模拟系统总的传热量并与傅立叶导热定律的求解结果相比较，验证数值模拟的可靠性。实验内容主要包括：（1）模拟一维稳态导热平板内的温度分布。（2）模拟一维稳态导热总的传热量。二、实例简介图1-1 导热计算区域示意图如图1-1所示，平板的长宽度远远大于它的厚度，平板的上部保持高温，平板的下部保持低温。平板的长高比为30，可作为一维问题进行处理。需要求解平板内的温度分布以及整个稳态传热过程的传热量。三、实例操作步骤1. 利用Gambit对计算区域离散化和指定边界条件类型步骤1：启动Gambit软件并建立新文件在路径C:Fluent.Incntbinntx86下打开gambit文件（双击后稍等片刻），其窗口布局如然后是建立新文件，操作为选择FileNew打开对话框。在ID文本框中输入onedim作为文件名，然后单击Accept按纽，在随后显示的对话框中单击Yes按纽保存。步骤2：创建几何图形选择OperationGeometryFace ，打开对话框。 在Width内输入30，在Height中输入1，在Direction下选择+X+Y坐标系，然后单击Apply，并在Global Control下点击 ，则出现几何图形。步骤3：网格划分（1）边的网格划分 当几何区域确定之后，接下来就需要对几何区域进行离散化，即进行网格划分。选择OperationMeshEdge，打开对话框。在Edges后面的黄色对话框中选中edge.1和edge.3。也可以采用Shift鼠标左键的方法选中edge.1和edge.3。然后在Spacing中选择Interval count，在其左边的对话框中输入100，即将这两个边各划分成100个等份。最后点击Apply确认。则出网格划分。 采用同样的方法对面的其它边进行网格划分，设定edge.2和edge.4的Spacing对应的数值为10，注意Spacing的类型仍然为Interval count，可以得到所示面上各边的网格划分。（2）面的网格划分对边进行网格划分实际上是对计算区域的边界进行离散化，计算区域的内部同样需要进行离散化，需要对计算区域进行面网格划分。选择OperationMeshFace ，打开对话框。在Faces后面的黄色框中选中face.1，选中之后，可以看到面上的边均变成红色，表示选择成功。对话框中的其它选项均保持默认值，此时Spacing的类型为Interval size，它左边的默认值为1。点击Apply确认。步骤4：边界条件类型的指定在指定边界条件之前，需要选定一个求解器，因为不同求解器的边界类型不一样。这里选择SolveFluent5/6，选择之后Gambit布局窗口标题栏中的Solve:Generic将变成Solve: Fluent5/6。选择OperationZone，打开对话框,指定边界条件的类型。首先指定面的上边为热源。具体操作为在Name右边的白色框中输入heat，选择Entity下面的类型为Edges，然后在Edges右边的黄色对话框中选择热源对应的边edge.3，点击Apply之后就将edge3定义成了热源。用同样的方法可以将下边定义成冷源cold。左右两条边可以不需要定义，保持Gambit默认即可。都定义完之后，可以得到边界名称和边界类型。步骤5：指定计算区域的类型Gambit默认的计算区域的类型为流体，而这里墙体内部的材料为固体，因此需要设置。设置方法为：选择OperationZone，打开窗口，选择Type为Solid，选择Entity为Faces，并在Faces右边的黄色对话框中选择面face.1，然后点击应用，即将计算区域的类型指定为固体区域。步骤6：网格文件的输出选择FileExportMesh打开输出文件的对话框。注意只有选择了Export 2-D(X-Y) Mesh选项之后才能输出为.msh文件。点击Accept之后，窗口下面的Transcript内出现Mesh was successfully written to onedim.msh，表示网格文件输出成功。2. 利用Fluent求解器进行求解利用Gambit软件绘制出几何图形、划分网格、指定边界类型以及输出Mesh文件，然后用Fluent将网格文件导入，便可以对其进行数值求解。步骤1：网格文件的读入、检查及显示启动Fluent的2D求解器之后，首先需要对网格文件进行读入并检查。启动Fluent后，在Versions中选择2d，点击Run按钮即可。（1）网格文件的读入选择FileReadCase在C:Fluent.Incntbinntx86下找到onedim.msh文件并将其读入（2）检查网格文件选择GridCheck对网格文件进行检查，这里要注意最小的网格体积（minimum volume）值一定要大于0。（3）显示网格选择DisplayGrid，出现网格显示对话框。网格文件的各个部分的显示可以通过Surfaces下面列表框中某个部分是否选中来控制。Surfaces下面列表框中的都被选中，此时单击Display，就会看到网格形状。步骤2：选择计算模型一维导热模型的控制方程只有能量方程，只需要选择DefineModelsEnergy，然后在出现的对话框中选中Energy Equation,单击OK即完成了方程的选择。步骤3：定义固体的物理性质选择DefineMaterials, 打窗口，在Material Type选项中选择solid，Fluent默认的固体材料为铝aluminum，我们假定平板的材料为铝，材料的属性取默认值，点击Change/Create按钮，再点击Close即可。步骤4：设置边界条件选择DefineBoundary Conditions，对计算区域的边界条件进行具体设置。对热源heat的边界类型wall点击set，出现的对话框，将默认的Thermal Condition下的heat Flux改为第一类边界条件Temperature，在Temperature右边的白色文本框内输入310。用同样的方法对冷源进行设置，其温度为300。即热源和冷源的温度差为10K。步骤5：求解设置（1）初始化选择SolveInitializeInitialize，打开对话框。依次点击Init、Apply和Close按钮。（2）残差设置选择SolveMonitorsResidual，打开对话框。选择Options下面的Plot复选项，则可在计算时动态地显示计算残差。并将energy右边的残差设定为1e-08，然后点击OK按钮（3）迭代计算选择SolveIterate，打开对话框。设置Number of Iterations 为200。然后单击Iterate按钮，就会显示计算过程。步骤6：保存结果选择FileWriteCase & Data，保存所有的设置和所有的数据。四、模拟实验结果经过上面的迭代计算，就可以查看模拟计算的结果。模拟结果的主要包括三个方面：（1）平板内部的温度分布；（2）平板内部的温度梯度；（3）平板总的传热量。（1）平板的温度分布选择DisplayContours，出现对话框，在Contours of 下选择Temperature 和Static Temperature，单击Display出现一个窗口，按住鼠标中间向右拖动将等温度图适当放大（图形的缩放、移动可以通过Display - Mouse Button来打开Mouse Buttons（鼠标按键）面板进行设定。），即可得到温度分布。在Contours窗口中选中Options中的Filled，可以得到的温度分布云图。（3）平板的总传热量选择ReportFluxes，打开对话框，在Options下选择Total Heat Transfer Rate， Boundaries下选择heat，然后单击Compute即可得到平板的总热流量为60726.6W。根据傅立叶导热定律计算的理论结果为60720W，相对误差为0.01%，表明结果正确。Fluent保存和编辑图形的方法：左键（或右键）点击显示窗口左上角的图标，点开后最下面有三个选项：Page Setup、Print及Copy to Clipboard，选择Page Setup，出现如下图所示窗口。按照上面窗口的设置完成后点击OK，再选择Copy to Clipboard，再到WORD中粘贴，即可得到彩色白底图形。在Page Setup中，Picture Format（图形格式）一栏中可以将图形格式设为Vector（矢量）或Raster（光栅）。其中Vector（矢量）格式清晰度高，但操作速度较慢，Raster（光栅）格式清晰度稍差，但操作速度较快，可以根据自己的需要决定图形格式。第四章数值模拟及计算4.1模型的建立需要以下建立模型（1） 管径长轴a=10，短轴b=10（2） 螺旋管半径Rc=50，螺距p=20实际操作步骤（1）利用Gambit对计算区域离散化和指定边界条件类型步骤1：启动Gambit软件并建立新文件在路径C:Fluent.Incntbinntx86下打开gambit文件然后是建立新文件，操作为选择FileNew打开入图4-1所示的对话框。图4-1建立新文件在ID文本框中输入onedim作为文件名，然后单击Accept按纽，在随后显示的图4-2对话框中单击Yes按纽保存。图4-2 确认保存对话框步骤2：创建几何图形选择OperationGeometryVertex打开图4-3所示的对话框。图4-3 创建点的对话框在Global内输入x，y，z各控制点坐标，单击apply按钮，改点就会在控制框内显示出来，重复这一操作，可以得出四个点。选择OperationGeometryedge打开图4-4所示的对话框。图4-4创建线的对话框 在Verticies内选中点，然后单击apply，就会出现一跳线 重复这个操作就会得到以下图形选择OperationGeometryedge 打开如图4-5所示的对话框图4-5 创建螺旋线参照图中的数据选项，单击apply按钮生成4-6所示的图形。图4-6几何图形的显示选择OperationGeometryFace ，打开图4-7所示的对话框。图4-7 面的创建在edges内选中4条边，然后单击apply，就会出现一个面 。选择选择OperationGeometryVolume，打开图4-8所示的对话框。图4-8体的创建步骤3：网格划分（1）面的网格划分对边进行网格划分实际上是对计算区域的边界进行离散化，计算区域的内部同样需要进行离散化，需要对计算区域进行面网格划分。选择OperationMeshFace ，打开图4-9所示的对话框。图4-9 面网格划分对话框在Faces后面的黄色框中选中face.1，选中之后，可以看到面上的边均变成红色，表示选择成功。对话框中的其它选项均保持默认值，此时Spacing的类型为Interval size，它左边的默认值为1。点击Apply确认。（2）体的网格划分选择OperationMeshvolmues ，打开图4-10所示的对话框图4-10体网格划分对话框在Volumes后面的黄色框中选中Volumes.1，选中之后，可以看到面上的边均变成红色，表示选择成功。对话框中的其它选项均保持默认值，此时Spacing的类型为Interval size，它左边的默认值为1。点击Apply确认。就得到如图4-11所示的几何图形。图4-11 网格的划分步骤4：边界条件类型的指定在指定边界条件之前，需要选定一个求解器，因为不同求解器的边界类型不一样。这里选择SolveFluent5/6，选择之后Gambit布局窗口标题栏中的Solve:Generic将变成Solve: Fluent5/6。选择OperationZone，打开对话框,指定边界条件的类型。首先指定面的上边为热源。具体操作为在Name右边的白色框中输入inlet，选择Entity下面的类型为faces，然后在faces右边的黄色对话框中选择热源对应的边faces.1，在type内选择类型为VELOCITY_INLET，点击Apply之后就将inlet定义成了入口面。用同样的方法可以将下边定义出口面outlet和内边加热面nei，保持Gambit默认即可。都定义完之后，可以得到图4-12的边界名称和边界类型。图4-12边界条件的指定步骤5：网格文件的输出选择FileExportMesh打开输出文件的对话框，如图1-15所示。图4-13输出文件对话框注意只有选择了Export 2-D(X-Y) Mesh选项之后才能输出为.msh文件。点击Accept之后，窗口下面的Transcript内出现Mesh was successfully written to onedim.msh，表示网格文件输出成功。2. 利用Fluent求解器进行求解利用Gambit软件绘制出几何图形、划分网格、指定边界类型以及输出Mesh文件，然后用Fluent将网格文件导入，便可以对其进行数值求解。步骤1：网格文件的读入、检查及显示启动Fluent的3D求解器之后，首先需要对网格文件进行读入并检查。启动Fluent后出现下面的窗口，在Versions中选择3d，点击Run按钮即可。（1）网格文件的读入选择FileReadCase在C:Fluent.Incntbinntx86下找到onedim.msh文件并将其读入，如图1-16所示。（2）检查网格文件选择GridCheck对网格文件进行检查，这里要注意最小的网格体积（minimum volume）值一定要大于0。（3）显示网格选择DisplayGrid，出现网格显示对话框，如图4-15所示。图4-15 网格显示对话框网格文件的各个部分的显示可以通过Surfaces下面列表框中某个部分是否选中来控制。如图1-17所示的Surfaces下面列表框中的都被选中，此时单击Display，就会看到的网格形状。步骤2：选择计算模型选择DefineModelsEnergy，然后在出现的如图4-16所示的对话框中选中Energy Equation,单击OK即完成了方程的选择。图4-16能量方程的选择对话框步骤3：定义水的物理性质选择DefineMaterials, 打开如图4-17所示窗口，在Material Type选项中选择fluid，Fluent默认的流体为空气air，我们假定流体为水water，材料的属性取默认值，点击Change/Create按钮，再点击Close即可。图4-17水的属性步骤4：设置边界条件选择DefineBoundary Conditions，对计算区域的边界条件进行具体设置。对流体的类型选择水，inlet的边界类型VELOCITY_INLET点击set，出现图4-18所示的对话框，将默认的速度改为0.0896。nei边的边界类型为wall，Thermal Condition下Temperature右边的白色文本框内输入350K如图4-19所示。图4-18 速度边界条件的设定图4-19热源边界条件的设定步骤5：求解设置（1）初始化选择SolveInitializeInitialize，打开如图4-20所示的对话框。Compute内选择inlet，依次点击Init、Apply和Close按钮。图4-20 初始化对话框（2）残差设置选择SolveMonitorsResidual，打开如图4-21所示的对话框。选择Options下面的Plot复选项，则可在计算时动态地显示计算残差。将x，y，z的残差设定为1e-06,并将energy右边的残差设定为1e-06，然后点击OK按钮。图4-21 残差设置对话框（3）迭代计算选择SolveIterate，打开如图4-22所示的对话框。设置Number of Iterations 为200。然后单击Iterate按钮，就会显示图4-22所示的计算过程。图4-22 迭代设置对话框 图1-25 迭代求解过程（3） SurfaceIso_Surface和Iso_clip创建面，在Line/Rake中创建线。（4） ReportSurface Integrals计算平均温度，阻力系数和Nu。步骤6：保存结果选择FileWriteCase & Data，保存所有的设置和所有的数据。Fluent保存和编辑图形的方法：左键（或右键）点击显示窗口左上角的图标，点开后最下面有三个选项：Page Setup、Print及Copy to Clipboard，选择Page Setup，出现如下图所示窗口。按照上面窗口的设置完成后点击OK，再选择Copy to Clipboard，再到WORD中粘贴，即可得到彩色白底图形。在Page Setup中，Picture Format（图形格式）一栏中可以将图形格式设为Vector（矢量）或Raster（光栅）。其中Vector（矢量）格式清晰度高，但操作速度较慢，Raster（光栅）格式清晰度稍差，但操作速度较快，可以根据自己的需要决定图形格式。第五章 计算结果及分析5.1.雷诺数对矩形螺旋管流动特性和传热特性的影响建立以下模型，其中截面边长a=10，边长b=10，螺旋半径Rc=50，螺距h=20，截取z=0 x0 y30的速度图，分别取雷诺数Re=500，1000，1500，2000，2500对比其流动特性和传热特性。5.1.1雷诺数对矩形截面的螺旋管流动特性的影响。（1）雷诺数对速度场的影响X方向 Re=500 Re=1000 Re=1500 Re=2000 Re=2500Y方向 Re=500 Re=1000 Re=1500 Re=2000 Re=2500Z方向 Re=500 Re=1000 Re=1500 Re=2000Re=2500从图5.1中可以看出，螺旋管的截面在X向速度等值线主要分布在两个区域中，分别位于长轴两端及中心线区域，长轴两端区域速度等值线大致相同且速度方向相反。在Y向速度等值线则沿长轴和短轴之间的一条斜线上分成形状大致相同且速度相反的两个区域。而在Z轴上。最