用Fluent对六个栅栏不同出口尺寸减压装置的数值模拟毕业设计.doc

用Fluent对六个栅栏不同出口尺寸减压装置的数值模拟毕业设计目录第一章 绪论11.1、减压节水装置研究背景及意义11.2、减压节水装置的设计原理11.3、减压节水装置设计21.4、应用前景和展望3第二章 CFD（计算流体动力学）42.1、软件总体介绍42.1.1、CFD的应用领域42.1.2、CFD的研究意义42.1.3、CFD的处理方法52.2、流体力学理论分析52.2.1、建立“力学模型”52.2.2、建立控制方程62.2.3、求解方程组62.2.4、对解进行分析解释62.3、数值计算62.4、发展前景7第三章 GAMBIT模型的创建和处理83.1、GAMBIT软件的简介83.1.1、功能83.1.2、特点83.2、GAMBIT几何建模93.3、GAMBIT网格的生成313.4、边界类型的设定32第四章 利用FLUENT 3D求解计算器进行求解374.1、FLUENT软件的简介374.1.1、FLUENT软件的应用374.1.2、基本特点374.1.3、程序的结构384.1.4、用FLUENT程序求解问题的步骤384.1.5、关于FLUENT求解器的说明394.2、利用FLUENT 3D求解计算器进行求解394.3、FLUENT数据处理。514.3.1、读入之前保存的Case和Data文件。514.3.2、显示网格。514.3.3、创建等值面。524.3.4、绘制压强图。54第五章 计算结果分析555.1、压强图555.2、绘制压强折线图。585.2.1、出口70mm各横截面压强595.2.2、出口90mm各横截面压强59第六章 结论61参考文献62致谢63外文原文64中文译文83第1章 绪论概述：本章主要介绍减压装置的研究意义，以及设计的原理、本次设计的的原理图。1.1、减压节水装置研究背景及意义地球表面的72%被水覆盖，储水量是很丰富的，共有14.5亿立方千米之多。尽管数量巨大，但是能直接被人们生产和生活利用的，却少得可怜。淡水资源仅占所有水资源的0.75%，有近70%的淡水固定在南极和格陵兰的冰层中，其余多为土壤水分或深层地下水，不能被人类利用。江河、湖泊、水库及浅层地下水等来源的水较易于开采供人类直接使用，但其数量不足世界淡水的1%，约占地球上全部水的0.007%。我国是一个缺水严重的国家。虽然我国的淡水资源总量为28000亿立方米，占全球水资源的6%，仅次于巴西、俄罗斯和加拿大，名列世界第四位。但是，我国的人均水资源量却只有2300立方米，仅为世界平均水平的1/4，在世界排第110位，是全球人均水资源最贫乏的国家之一。严峻的水资源问题一直是我国可持续发展中一个不能回避的难点，水资源短缺已经严重影响了我国的经济发展和社会和谐，成了我国无法回避的世纪挑战。在高楼大厦林立的今天，通过加压把水输送到各个楼层，然而出口的压强又不能太大，压强过大可能会导致水管破裂，这样既浪费宝贵水源又破坏管道。 所以研究减压装置具有节约水资源的保护管道的双重意义。 1.2、减压节水装置的设计原理根据有压管道消能原理。当进水断面突然缩小时，会产生局部的阻力损失，有水力学知识可知，该局部阻力损失可表示为： （1.1）其中，h1局部水头损失的阻力 局部水头损失的阻力系数 V水流速度 g重力加速度根据以上的原理可以看出，只需通过调节过水断面突然缩小的面积，就可以调节装置的效能效果了。本装置通过齿轮传动原理，带动内套筒的转动，从而实现了过水面积的调节。当主管道中的水流过此装置时，内外套筒装置上底面的圆形板状体使过流面积突然减小，起到了消能作用；在水流通过内外套筒栅格状侧边壁时，突然扩大形状的内外套筒侧边壁又使过流面积增大，起到了减压节水的作用。根据内外筒套重合程度的不同，消能效率将随之发生改变，达到不同的减压节水效果。1.3、减压节水装置设计 如图1-1所示，为减压节水装置本装置由底面直径为100mm，长350mm的圆柱体及内部减压结构套筒组成，内部装置上部结构是一直径为50mm，长3mm的小型圆柱体，下部结构是一外径为100mm，内径为70mm，宽度为3mm的圆环体，上下两部分由六条均匀分布的栅栏结构相连，栅栏结构宽度为13mm，厚度为3m（内部结构壁厚均为3mm）。100 35050 70/90 100 200 图 1-1 减压节水装置1.4、应用前景和展望在民用方面，本装置适合用于居住于高层建筑的低层用户使用，市场需求大。同时本装置具有可调节行，制作简单，成本低，安装操作方便。可以满足高层建筑中不同楼层用户的减压节水需要，具有很强的通用性。在水利工程方面，本装置作为一种消能装置，可以广泛应用于解决水库长距离输水线路的安全问题，农业灌溉系统中，及水电站、泵站的水击消除上，具有广阔的市场前景。若能成批量生产，前景广阔。 这篇论文的写作以及系统开发的过程，也是我越来越认识到自己知识与经验缺乏的过程。虽然，我尽可能地收集材料，竭尽所能运用自己所学的知识进行论文写作和系统开发，但论文还是存在许多不足之处,系统功能并不完备。第2章 CFD（计算流体动力学）概述：本章主要介绍CFD，软件(Computational Fluid Dynamics），即计算流体动力学, 简称CFD。CFD是近代流体力学，数值数学和计算机科学结合的产物，是一门具有强大生命力的边缘科学。它以电子计算机为工具，应用各种离散化的数学方法，对流体力学的各类问题进行数值实验、计算机模拟和分析研究，以解决各种实际问题。2.1、软件总体介绍2.1.1、CFD的应用领域 计算流体力学和相关的计算传热学，计算燃烧学的原理是用数值方法求解非线性联立的质量、能量、组分、动量和自定义的标量的微分方程组，求解结果能预报流动、传热、传质、燃烧等过程的细节，并成为过程装置优化和放大定量设计的有力工具。计算流体力学的基本特征是数值模拟和计算机实验，它从基本物理定理出发，在很大程度上替代了耗资巨大的流体动力学实验设备，在科学研究和工程技术中产生巨大的影响。是目前国际上一个强有力的研究领域, 是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术, 广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计、HAVC&R 等诸多工程领域,板翅式热换器设计是CFD 技术应用的重要领域之一。2.1.2、CFD的研究意义CFD 在最近20 年中得到飞速的发展, 除了计算机硬件工业的发展给它提供了坚实的物质基础外, 还主要因为无论分析的方法或实验的方法都有较大的限制, 例如由于问题的复杂性, 既无法作分析解, 也因费用昂贵而无力进行实验确定, 而CFD 的方法正具有成本低和能模拟较复杂或较理想的过程等优点。经过一定考核的CFD软件可以拓宽实验研究的范围, 减少成本昂贵的实验工作量。在给定的参数下用计算机对现象进行一次数值模拟相当于进行一次数值实验, 历史上也曾有过首先由CFD 数值模拟发现新现象而后由实验予以证实的例子。CFD软件一般都能推出多种优化的物理模型，如定常和非定常流动、层流、紊流、不可压缩和可压缩流动、传热、化学反应等等。对每一种物理问题的流动特点, 都有适合它的数值解法, 用户可对显式或隐式差分格式进行选择, 以期在计算速度、稳定性和精度等方面达到最佳。CFD软件之间可以方便地进行数值交换, 并采用统一的前、后处理工具, 这就省却了科研工作者在计算机方法、编程、前后处理等方面投入的重复、低效的劳动, 而可以将主要精力和智慧用于物理问题本身的探索上。2.1.3、CFD的处理方法CFD软件的一般结构由前处理、求解器、后处理三部分组成。前处理、求解器及后处理三大模块, 各有其独特的作用, 分别表示如下:前处理求解器后处理作用a. 几何模型b. 划分网格a. 确定CFD 方法的控制方程b. 选择离散方法进行离散c. 选用数值计算方法d. 输入相关参数速度场、温度场、压力场及其它参数的计算机可视化及动画处理目前比较好的CFD软件有：CFX、Fluent、CF design、除了Fluent、CF design是美国公司的软件外，其它三个都是英国公司的产品 。2.2、流体力学理论分析根据流体运动的普遍规律如质量守恒、动量守恒、能量守恒等，利用数学分析的手段，研究流体的运动，解释已知的现象，预测可能发生的结果。理论分析的步骤大致如下：2.2.1、建立“力学模型”一般做法是：针对实际流体的力学问题，分析其中的各种矛盾并抓住主要方面，对问题进行简化而建立反映问题本质的“力学模型”。流体力学中最常用的基本模型有：连续介质（见连续介质假设）、牛顿流体、不可压缩流体、理想流体（见粘性流体）、平面流动等。2.2.2、建立控制方程针对流体运动的特点，用数学语言将质量守恒、动量守恒、能量守恒等定律表达出来，从而得到连续性方程、动量方程和能量方程。此外，还要加上某些联系流动参量的关系式（例如状态方程），或者其他方程。这些方程合在一起称为流体力学基本方程组。流体运动在空间和时间上常有一定的限制，因此，应给出边界条件和初始条件。整个流动问题的数学模式就是建立起封闭的、流动参量必须满足的方程组，并给出恰当的边界条件和初始条件。2.2.3、求解方程组在给定的边界条件和初始条件下，利用数学方法，求方程组的解。由于这方程组是非线性的偏微分方程组，难以求得解析解，必须加以简化，这就是前面所说的建立力学模型的原因之一。力学家经过多年努力，创造出许多数学方法或技巧来解这些方程组（主要是简化了的方程组），得到一些解析解。2.2.4、对解进行分析解释求出方程组的解后，结合具体流动，解释这些解的物理含义和流动机理。通常还要将这些理论结果同实验结果进行比较，以确定所得解的准确程度和力学模型的适用范围。2.3、数值计算前面提到的采用简化模型后的方程组或封闭的流体力学基本方程组用数值方法求解。电子计算机的出现和发展，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性。数值方法可以部分或完全代替某些实验，节省实验费用。数值计算方法最近发展很快，其重要性与日俱增。2.4、发展前景从阿基米德到现在的二千多年，特别是从20世纪以来，流体力学已发展成为基础科学体系的一部分，同时又在工业、农业、交通运输、天文学、地学、生物学、医学等方面得到广泛应用。今后，人们一方面将根据工程技术方面的需要进行流体力学应用性的研究，另一方面将更深入地开展基础研究以探求流体的复杂流动规律和机理。后一方面主要包括：通过湍流的理论和实验研究，了解其结构并建立计算模式；多相流动；流体和结构物的相互作用；边界层流动和分离；生物地学和环境流体流动等问题；有关各种实验设备和仪器等第3章 GAMBIT模型的创建和处理概述：本章主要介绍GAMBIT软件的功能和特点以及软件的使用。并利用软件进行建模的过程。3.1、GAMBIT软件的简介3.1.1、功能GAMBIT是为了帮助分析者和设计者建立并网格化计算流体力学（CFD）模型和其它科学应用而设计的一个软件包。GAMBIT通过它的用户界面（GUI）来接受用户的输入。GAMBIT GUI简单而又直接的做出建立模型、网格化模型、指定模型区域大小等基本步骤，然而这对很多的模型应用已是足够了。面向CFD分析的高质量的前处理器，其主要功能包括几何建模和网格生成。由于GAMBIT本身所具有的强大功能，以及快速的更新，在目前所有的CFD前处理软件中，GAMBIT稳居上游。3.1.2、特点1.ACIS内核基础上的全面三维几何建模能力，通过多种方式直接建立点、线、面、体，而且具有强大的布尔运算能力，ACIS内核已提高为ACIS R12。该功能大大领先于其它CAE软件的前处理器；2. 可对自动生成的Journal文件进行编辑，以自动控制修改或生成新几何与网格；3. 可以导入PRO/E、UG、CATIA、SOLIDWORKS、ANSYS、PATRAN等大多数CAD/CAE软件所建立的几何和网格。导入过程新增自动公差修补几何功能，以保证GAMBIT与CAD软件接口的稳定性和保真性，使得几何质量高，并大大减轻工程师的工作量；4.新增PRO/E、CATIA等直接接口， 使得导入过程更加直接和方便；5.强大的几何修正功能，在导入几何时会自动合并重合的点、线、面；新增几何修正工具条，在消除短边、缝合缺口、修补尖角、去除小面、去除单独辅助线和修补倒角时更加快速、自动、灵活，而且准确保证几何体的精度；6.G/TURBO模块可以准确而高效的生成旋转机械中的各种风扇以及转子、定子等的几何模型和计算网格；7.强大的网格划分能力，可以划分包括边界层等CFD特殊要求的高质量网格。GAMBIT中专用的网格划分算法可以保证在复杂的几何区域内直接划分出高质量的四面体、六面体网格或混合网格；8.先进的六面体核心(HEXCORE)技术是GAMBIT所独有的，集成了笛卡尔网格和非结构网格的优点，使用该技术划分网格时更加容易，而且大大节省网格数量、提高网格质量；3.2、GAMBIT几何建模首先，启动GAMBIT软件，点击GAMBIT软件界面如图图2.1GAMBIT软件主界面新建一个文件夹后，开始进行GAMBIT建模.(以下所有单位均为m)。步骤如下：3.2.1、用鼠标左击操作工具栏（GAMBIT建模界面右上角）GEOMETRY；3.2.2、在GEOMETRY子操作板上用鼠标右击VOLUME命令按钮，在弹出下拉菜单中，左击按钮，打开Create Real Brick对话框图2.2Create Real Brick 对话框3.2.3、在Create Real Brick对话框中的Height中输0.35Rdius1中输入0.025adius2可保留空GAMBIT默认它跟Radius1有相同的值，其余保留默认值，设置如图图2.3圆柱参数设置图点击Apply，这时在左下角将显示信息说您已经创建一个物体“volume.1”,鼠标左击右下角的，此时在图形窗口上显示如图图2.4圆柱图 3.2.4按照第三步的方法，建立两个高度为3mnm,半径分别为50mm和35m的圆柱（或半径分别为50mm和45mm的圆柱），如图所示。图2.5新建两个高度为3mm半径50m和35mm的小圆柱图2.6新建两个高度为3mm半径50m和45mm的小圆柱3.2.5、将新建的两个小圆柱进行布尔运算，合成一个圆柱，具体操作为，鼠标右击Volume中的，选中，打开Subtract Real Volumes 对话框，如图2.7图2.7 Subtract Real Volumes 对话框在Volume中选择Volume2,即第4步建立的半径为50mm的圆柱，在Subtract Volumes中选择Volume3,即半径为35mm（或半径为45mm）的圆柱，如图2.8图2.8布尔运算设置点击Apply。3.2.6、移动圆柱。鼠标左击，打开Move/Copy Volumes 对话框，在Pick中选中Volume.2,Z轴移动0.2m。设置如图2.9图2.9 Move/Copy Volumes 对话框设置点击Apply，显示如图2.10图2.10圆环平移3.2.7、按第三步的方法新建一个高度为3mm，半径为25mm的圆盘Volume.3，如图2.11。图2.11半径25mm的圆盘Volume.33.2.8、按第六步移动圆柱Volume.3。 将圆盘Volume.3沿Z轴方向平移100m，设置如图2.12.图2.12平移设定点击Appiy.平移后，模型如图2.13图2.13圆盘平移3.2.9、在圆环和小圆盘上建矩形栅栏切面 步骤为：3.2.9.1在Z=0的圆面上新建两个长13mm，宽3mm的矩形。（face12） 鼠标左击，在Face中左击，打开Create Real Rectangular Face对话框，Widh中输入0.003，Height中输入0.013.设置如图2.14图2.14矩形设置点击Apply.模型如图2.15图2.15-1出口为70mm栅栏切面在地面的投影图2.15-2出口为90mm栅栏切面在地面的投影3.2.9.2、平移矩形face12.按第6步把矩形face12沿X轴方向平移0.02m,Z轴方向平移0.1m.到圆盘切面。设置如图2.16图2.16矩形平移设置点击Apply,模型如图2.17，2.18图2.17-1 出口为70mm,矩形栅栏切面平移图2.17-2 出口为90mm,矩形栅栏切面平移图2.18-1 出口为70mm, 矩形栅栏在圆盘上的投影图2.18-2出口为90mm,矩形栅栏在圆盘上的投影3.2.9.3、按1步在建一个尺寸相同的矩形face13,如图2.19图2.19-1出口为70mm,矩形栅栏切面平移图2.19-2出口为90mm,矩形栅栏切面平移按2步平移face13到圆柱切面上。把face13沿X轴方向平移0.047m,沿Z轴方向平移0.2m.设置如图2.20图2.20 矩形face13平移设置点击Apply,模型如图2.21，2.22图2.21-1出口为70mm,矩形栅栏切面平移图2.21-2出口为90mm,矩形栅栏切面平移图2.22-1出口为70mm,矩形face13在圆柱上的投影图2.22-2出口为90mm,矩形face13在圆柱上的投影3.2.10、点成线。把矩形face12与face13上的四个点相对应连接成线。在Geometry中左击（点成线界面），选择，打开Create Straight Edge对话框，如图2.23.图2.23 Create Straight Edge 点成线设置在Vertices中选择相对应的点。（11、15），（12、16），（13、17），（14、18）点击Apply,模型如图2.24图2.24点成线3.2.11、线成面。将形成栅栏的六个面中，每个面由四条线形成。在Geometry中选择，在弹出菜单中选择，打开Create Face from Wireframe 对话框，如图2.25.图2.25 Create Face from Wireframe 线成面界面在Edges中选择对应的四条线组成一个面。Edge(13、17、21、22)，（12、16、20、21），(11、15、79、20)，（14、18、19、22）。点击Apply，模型如图2.26图2.26 线成面3.2.12、面成体。将形成栅栏的六个面组成一个整体。在Geometry中选择，在弹出菜单中选择，打开Stitch Faces 界面，如图2.27.图2.27 Stitch Faces 线成面界面在Faces中选择face12.13.14.15.16.17六个面。点击Alppy,组成一个面Volume4.如图2.28图2.28 一条栅栏连接3.2.13、复制剩余的5个栅栏。 在Geometry中选择，选择，打开Move/Copy Volumes 对话框，设置如图2.29，（因为是六个栅栏所以每个旋转角度为60度，首先选择Volume4,点击Alppy,再选择Volume5进行下一次复制，依次复制出5个栅栏）。图2.29 栅栏旋转复制界面设置点击Alppy,模型如图2.30，2.31图2.30-1 出口70mm,六栅栏初步模型图2.30-2 出口90mm,六栅栏初步模型图2.31 出口70mm的模型切面投影图图2.32出口90mm的模型切面投影图3.2.14、利用布尔运算，将整个装置除内部结构外的部分（即实际运行中会被流体占据的部分）合成一个整体。鼠标右击Volume中的，选择，弹出Unite Real Volumes对话框，在Volumes中选择volume1.如图2.33.图2.33布尔运算界面点击Apply,如图2.34所示。图2.34 完成的模型这样，就完成了研究装置的GAMBIT模型的建立。3.3、GAMBIT网格的生成3.3.1、单机Operation ,选择Mesh中的，选择，打开Mesh Volumes对话框，各参数设置如图2.35图2.35体网格参数点击Apply，生成网格，如图2.36图2.36生成的体网络3.3.2、在Global Control 01栏中选择 EXAMINE MESF 命令按钮，打开Examine Mesh对话框如图2.37所示。33图2.37 Examine Mesh对话框（出口70mm）31图2.38 Examine Mesh对话框（出口90mm）在Display Type 中选择Range,在 Quality Type菜单中选择Equia Angle Skew,从绿色垂直条形码中点击一条来观察各个元素在确定的质量范围内，移动Upper 和 Lower滑块来重新定义质量网格范围。生成网格数目出口半径（mm）7090网格数目8747309596983.3.3保存。3.4、边界类型的设定GAMBIT生成网格后，需要定义边界条件，步骤如下。3.4.1、打开生成的网格，单机Operation中的命令按钮，选择Zones中的，在弹出Specify Boundary Types对话框，如图2.38所示。图2.38 Specify Boundary Types对话框3.4.2、设置入口边界类型。在Name中填入in，Type选择VELOCITY_INLET,Face中选中face.1,即入口断面，具体设置如图2.39图2.39入口断面边界设定点击Apply.3.4.3、设置出口边界类型。在Name中填入out,Type选择OUTFLOW,Faces中选中face.3,即出口断面，具体设置如图2.40图2.40出口边界条件的设定点击Apply.3.4.4、设置边缘边界类型。在Name中填入BOUNDARY, Type选择WALL, Faces中选中face.3和face.1以外的所有面，具体设置如图2.41，2.42.图2.41面的选择图2.42边界设定点击Apply.3.4.5、保存文件，并到处网格。图2.43网格文件导出第4章 利用FLUENT 3D求解计算器进行求解概述：本章主要介绍了FLUENT软件，并利用FLUENT软件进行模型的数值模拟计算。4.1、FLUENT软件的简介 Fluent是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%，凡是和流体、热传递和化学反应等有关的工业均可使用。它具有丰富的物理模型、先进的数值方法和强大的前后处理功能。4.1.1、FLUENT软件的应用CFD商业软件FLUENT，是通用CFD软件包，用来模拟从不可压缩到高度可压缩范围内的复杂流动。由于采用了多种求解方法和多重网格加速收敛技术，因而FLUENT能达到最佳的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网络技术及成熟的物理模型，使FLUENT在转换与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面有广泛应用。4.1.2、基本特点FLUENT软件具有以下特点：1.FLUENT软件采用基于完全非结构化网格的有限体积法，而且具有基于网格节点和网格单元的梯度算法；定常/非定常流动模拟，而且新增快速非定常模拟功能；2.FLUENT软件中的动/变形网格技术主要解决边界运动的问题，用户只需指定初始网格和运动壁面的边界条件，余下的网格变化完全由解算器自动生成。网格变形方式有三种：弹簧压缩式、动态铺层式以及局部网格重生式。其局部网格重生式是FLUENT所独有的，而且用途广泛，可用于非结构网格、变形较大问题以及物体运动规律事先不知道而完全由流动所产生的力所决定的问题；3.FLUENT软件具有强大的 网格支持能力，支持界面不连续的网格、混合网格、动/变形网格以及滑动网格等。值得强调的是，FLUENT软件还拥有多种基于解的网格的自适应、动态自适应技术以及动网格与网格动态自适应相结合的技术；4.FLUENT软件包含三种算法：非耦合隐式算法、耦合显式算法、耦合隐式算法，是商用软件中最多的；5.FLUENT软件包含丰富而先进的物理模型，使得用户能够精确地模拟无粘流、层流、湍流。湍流模型包含Spalart-Allmaras模型、k-模型组、k-模型组、雷诺应力模型(RSM)组、大涡模拟模型(LES)组以及最新的分离涡模拟(DES)和V2F模型等。另外用户还可以定制或添加自己的湍流模型；4.1.3、程序的结构FLUENT程序软件包由以下几个部分组成：1.GAMBIT用于建立几何结构和网格的生成。2.FLUENT用于进行流动模拟计算的求解器。3.prePDF用于模拟PDF燃烧过程。4.TGrid用于从现有的边界网格生成体网格。5.Filters(Translators)转换其他程序生成的网格，用于FLUENT计算。可以接口的程序包括：ANSYS，I-DEAS，NASTRAN，PATRAN等。4.1.4、用FLUENT程序求解问题的步骤利用FLUENT软件进行求解的步骤如下：1.确定几何形状，生成计算网格(用GAMBIT，也可以读入其他指定程序生成的网格)。2.输入并检查网格。3.选择求解器(2D或3D等)。4.选择求解的方程：层流或湍流(或无粘流)，化学组分或化学反应，传热模型等。确定其他需要的模型，如：风扇、热交换器、多孔介质等模型。5.确定流体的材料物性。6.确定边界类型及其边界条件。7.条件计算控制参数。8.流场初始化。9.求解计算。10.保存结果，进行后处理等。4.1.5、关于FLUENT求解器的说明1.FLUENT2D二维单精度求解器；2.FLUENT3D三维单精度求解器：3.FLUENT2ddp二维双精度求解器；4.FLUENT 3ddp三维双精度求解器。4.2、利用FLUENT 3D求解计算器进行求解打开FLUENT软件，在Versions中选择3d ,其他保留默认设置，如图4.1所示。图4.1 求解器选择单机Run，软件界面如图4.2 所示。图4.2 FLUENT 软件主界面4.2.1、依次点击FileReadCase,选择之前导出的网络文件ChuKou70.msh或（ChuKou90.msh）。如图4.3所示。图4.3 导出网络文件4.2.2、检查网络。依次点击GridCheck, FLUENT将会对网格进行检查，并将结果在信息反馈窗口中显示出来，如图4.4所示，其中注意最小体积（minimum volume）一项，要确保为正直，否则无法计算。图4.4检查网络4.2.3、依次点击 GridSmooth/Swap，打开Smooth/Swap Gird对话框，如图4.5所示。图4.5 Smooth/Swap Gird对话框点击Swap.信息反馈窗口中显示，再点击一次Swap，使得Number face swapped 变为0，如图4.6所示图 4.64.2.4、建立求解模型。1、依次点击DefineModelsSolve,打开Solve对话框，如图4.7所示。图4.7 Solve对话框保留默认设置，点击 OK .2、选择湍流模型。依次点击DefineModelsViscous，打开Viscous Models对话框，如图4.8所示。图4.8Viscous Models对话框Model选择k-e模型，各项设置保留默认，点击 OK.如图4.9所示。图4.9 求解器模型的选择4.2.5设置流体属性。依次点击DefineMaterials，打开Materials对话框，如图4.10所示。图4.10 Materials对话框1.点击Fluent Database，打开Fluent Database Materials对话框，如图4.11所示。图4.11Fluent Database Materials对话框2.在Fluent Database Materials对话框中选择water-liquid，点击Coppy。3.在Materials对话框中单机Chane/Create,关闭对话框。4.2.6设置边界条件。依次点击Define- Boundary Conditions，打开Boundary Conditions对话框，如图4.12所示。图4.12 Boundary Conditions对话框1.设置管道的入口边界条件。在Zone下拉列表中选取in,如图4.13所示。在Type下拉列表中选取Velocity-inlet，点击Set，打开边界条件设置对话框如图4.14所示。图4.13 Boundary Conditions对话框点击Set，打开边界条件设置对话框如图4.14所示。图4.14 Velocity-inlet对话框在Velocity Magnitude（m/s）中输入0.3，其余设置保留默认。点击 OK .2.设置管道出口边界条件。在Zone下拉列表中选取out,在Type下拉列表中选取outflow,如图4.15所示。图4.15 Boundary Conditions对话框点击Set，打开边界条件设置对话框如图4.16所示。图4.16 出口边界条件点击 OK 。3.设置液体边界条件。在Zone下拉列表中选取fluid,在Type下拉列表中选取fluid,如图4.17所示。图4.17 流体类型的选择在Material Name下拉列表中选择water-liquid,其余设置保留默认，点击 OK。4.2.7、求解初始化。依次点击SloveInitializeInitialize，打开Solution Initialization对话框，如图4.18所示。图4.18 Solution Initialization对话框在Computer From中选择in，点击 Init。4.2.8、设置求解器参数。依次点击SloveMonitorsResidual。打开Residual Monitors对话框，如图4.19所示。图4.19 Residual Monitors对话框在Options中勾选Print和Plot，其余保留默认设置，如图4.19,点击 OK .4.2.9、依次点击FileWriteCase,保存数据。4.2.10、迭代计算。依次点击SloveIterate，打开Iterate对话框，如图4.20所示。图4.20 Iterate对话框在Number of Iterations中填入300，即进行300次迭代计算。点击Iterate，开始迭代，出现残差图。经过一段事件后，在迭代102次后收敛，如图4.21（4.22）所示。图4.21 速度为0.3m/s，出口为70mm时的残差图图4.22 速度为0.3m/s，出口为90mm时的残差图依次点击FileWriteDate，保存数据。4.3、FLUENT数据处理。用FLUENT初步计算后，还需要对所得数据进行进一步处理，才能得到我们所需要的数据。4.3.1、读入之前保存的Case和Data文件。依次点击FileRead- Case&Data，读取文件。4.3.2、显示网格。点击DisplayGrid，弹出Grid Display 对话框，如图4.23所示。图4.23 Grid Display 对话框点击Display，窗口显示如图4.24所示。图4.24 显示网格4.3.3、创建等值面。点击Surface -Iso-Surface，打开Iso-Surface对话框，如图4.25所示。图4.25 Iso-Surface对话框取Z=0.05的横截面为等值面，参数设置如图4.26所示。图4.26 Z=0.05等值面设定 点击Greate,创建等值面。依次创建Z=0.1,Z=0.15,Z=0.2,Z=0.3等值面。4.3.4、绘制压强图。取等值面Z=0.05的压强，点击DisplaContours,打开Contours对话框，参数设置如图4.27所示。图4.27 Contours对话框设定点击Display。压强图如下。图4.28 k-e模型出口70mm,速度0.3m/s,z=0.05截面的压强图第五章 计算结果分析概述：对计算结果进一步分析，在流速为0.3m/s时，对出口为70mm和90mm的装置的效能进行分析。取z=0.05,z=0.1,z=0.15,z=0.2,z=0.3的等值面。5.1、压强图取z=0.05,z=0.1,z=0.15,z=0.2,z=0.3的等值面图5.1 k-e模型出口70mm,速度0.3m/s,z=0.05截面的压强图图5.2 k-e模型出口70mm,速度0.3m/s,z=0.1截面的压强图图5.3 k-e模型出口70mm,速度0.3m/s,z=0.15截面的压强图图5.4 k-e模型出口70mm,速度0.3m/s,z=0.2截面的压强图图5.5 k-e模型出口70mm,速度0.3m/s, z=0.2截面的压强图图5.6 k-e模型出口70mm,速度0.3m/s,z=0.3截面的压强图5.2、绘制压强折线图。依次点击ReportSurface Intergrals打开Surface Intergrals对话框，如图5.7所示图5.7 Surface Intergrals对话框在Report Type中选择Area-Weighted Average；在Filed Variable中选择Pressure（压强），Total Pressure;在Surfaces中选择in、out、z=0.05、z=0.1、z=0.15、z=0.2、z=0.3。点击Compute。得到数据如下表。5.2.1、出口70mm各横截面压强表1 、出口70mm各横截面压强出口70mm压强图位置（m)00.050.10.150.20.250.30.35压强（Pa）46.984188-3.213038929.59367-435.40082-670.01044-2256.8242-2275.0291-2295.8图5.12 出口70mm各横截面压强5.2.2、出口90mm各横截面压强表2、出口90mm各横截面压强出口90mm压强图位置（m)00.050.10.150.20.250.30.35压强（Pa）48.146492-3.77182840.921375-481.02695-635.92529-970.93683-997.06195-1019.48图5.13 出口90mm各横截面压强小结：两条压强图总体都是呈下降趋势；都是在0.1处压强开始迅速下降，在出口处降到最低；出口70mm的减压装置比90mm的减压装置减压效果更好表3、出口不同模型的减压效果出口不同模型的减压效果 In（Pa）Out(Pa)压强差(Pa)70mm46.984188-2295.82342.78418890mm48.146492-1019.481067.626492由表5可以得出出水口70mm的装置比出水口90mm的减压装置减压效果更好。第六章 结论通过对减压装置进行物理建模，然后对模型精心数值模拟计算，绘制减压装置各断面的压力和速度图，在进水速度为0.3m/s,减压装置其余属性保持相同时得出以下结论：1.水流经过减压装置后，压力明显降低；2.压强在管内的变化情况是在0.1处开始呈下降趋势；在减压装置处，压力下降幅度更大；3.两种出口尺寸不同的减压装置都可以起到减压的效果，都是在减压装置所在区域压强开始下降；但是出口70mm的装置比出口90mm的装置减压效果更好；参考文献1 张红霞. 可调式减压消能装置 D.太原：太原理工大学，20092 白丹，宋立勋，郭霖. 滴灌减压消能件D.西安：西安理工大学，20123 徐志通，余雪松. 高层建筑压力流雨水排水系统立管减压消能装置R.北京：北京泰宁科创雨水利用技术股份有限公司4 曲景学,杨永全,张建民,许唯临.消能孔板空化特性的数值模拟J.四川大学学报(工程科学版),2001,(第3期). 5 余挺,田忠,王韦,许唯临.收缩式洞塞泄洪洞的消能和空化特性J.水利学报,2011,(第2期).6 付宁宁. 热管在油浸式变压器中温度场分布的Fluent数值模拟D.河北：河北工业大学7 王立新. 回路型重力热管蒸发腔沸腾传热的FLUENT数值分析.D.浙江：浙江大学8 胡俊明. 基于Fluent的波浪辐射与绕射问题数值模拟研究.D.哈尔滨：哈尔滨工程大学9 张萌萌. 基于FLUENT的北极海冰三维温度场数值仿真和融化模拟.D.武汉：武汉理工大学10 宫汝志. 汽轮发电机通风冷却风扇CFD设计与优化.D.哈尔滨：哈尔滨工业大学11 张振江. 船舶CFD网格自动生成技术的开发及其应用.D.上海：上海交通大学12 王新龙. 轻烃外输泵叶轮内部流场的CFD模拟研究.D.东北：东北石油大学13 曹应佳. 基于FLUENT的波节管换热器性能数值仿真及其结构优化.D.武汉：武汉理工大学14 郭春雨. 螺旋桨与舵附推力鳍相互干扰水动力性能数值计算.D.哈尔滨：哈尔滨工业大学15 高宏适. 数值模拟技术在金属塑性加工中的应用.N. 世界金属导报, 2012-11-2016 刘凡. 油藏数值模拟数模研究.R. 中国石油报, 2010-08-0417 Arvay; Adam. Proton exchange membrane fuel cell modeling and simulation using Ansys FluentD. Arizona State University.bChemistry.18 Sweeten; Benjamin. CFD analysis of UAVs using VORSTAB, FLUENT, and advanced aircraft analysis software.D. The University of Kansas19 Dakin; Glenn H. IT fluent information literacy skill: A Delphi study of an essential course offeringD. University of Phoenix致谢毕业论文暂告收尾，这也意味着我在太原理工大学现代科技学院的四年的学习生活既将结束。回首既往，自己一生最宝贵的时光能于这样的校园之中，能在众多学富五车、才华横溢的老师们的熏陶下度过，实是荣幸之极。在这四年的时间里，我在学习上和思想上都受益非浅。这除了自身努力外，与各位老师、同学和朋友的关心、支持和鼓励是分不开的论文的写作是枯燥艰辛而又富有挑战的。老师的谆谆诱导、同学的出谋划策及家长的支持鼓励，是我坚持完成论文的动力源泉。在此，我特别要感谢我的导师张老师、孙老师以及秦师兄和李师兄。从论文的选题、文献的采集、框架的设计、结构的布局到最终的论文定稿，从内容到格式，从标题到标点，她都费尽心血。没有张老师、孙老师的辛勤栽培、孜孜教诲，以及两位师兄的指导和改正，就没有我论文的顺利完成。感谢我们小组的各位同学，与他们的交流使我受益颇多。最后要感谢我的家人以及我的朋友们对我的理解、支持、鼓励和帮助，正是因为有了他们，我所做的一切才更有意义；也正是因为有了他们，我才有了追求进步的勇气和信心。时间的仓促及自身专业水平的不足，整篇论文肯定存在尚未发现的缺点和错误。恳请阅读此篇论文的老师、同学，多予指正，不胜感激！外文原文SAND EROSION MODEL IMPROVEMENT FOR