第章-fluｅnt简单算例
————————————————————————————————作者: ————————————————————————————————日期：
?
2004-06 FLUENT流体工程仿真计算实例与应用韩占忠王敬兰小平北京理工大学出版社 第一章流体力学基础与fluent简介 第二章二维流动与传热的数值计算 第一节冷、热水混合器内部二维流动 第二节喷管内二维非定常流动 第三节三角翼的可压缩外部绕流 第四节三角翼不可压缩的外部绕流(空化模型应用) 第五节vof模型的应用 第六节组分传输与气体燃烧 第三章三维流动与传热的数值计算 第一节冷、热水混合器内的三维流动与换热 第二节粘性流体通过圆管弯头段的三维流动 第三节三维稳态热传导问题 第四节动网格问题 第五节叶轮机械的mixing plane模型 2004-09 计算流体动力学分析CFD软件原理与应用王福军清华大学出版社（偏重理论） 第1章计算流动力学基础知识 第2章基于有限体积法的控制方程离散 第3章基于SIMPLE算法的流场数值计算
第4章三维流模型及其在CFD中的应用 第5章边界条件的应用 第6章网格的生成 第7章FLUENT软件的基本用法 第8章CFD综合应用实例 2007-02 FLUENT技术基础与应用实例王瑞金张凯王刚清华大学出版社 第1章Fluent概述 第2章流体力学基础知识 第3章流体力学数值模拟基础 第4章Fluent软件介绍 第5章速度场的计算 第6章温度场的计算 第7章多相流模型 第8章凝固和融化模型 第9章可动区域中流动问题的模拟 第10章动网格模型 第11章UDF和UDS 第12章Fluent并行计算 第13章Tecplot软件 2008-07 Fluent高级应用与实例分析江帆，黄鹏清华大学出版社第1章 CFD基础 第2章Fluent基本介绍
FLUENT 教程 赵玉新 I、目录 第一章、开始 第二章、操作界面 第三章、文件的读写 第四章、单位系统 第五章、读入和操作网格 第六章、边界条件 第七章、物理特性 第八章、基本物理模型 第九章、湍流模型 第十章、辐射模型 第十一章、化学输运与反应流 第十二章、污染形成模型 第十三章、相变模拟 第十四章、多相流模型 第十五章、动坐标系下的流动 第十六章、解算器的使用 第十七章、网格适应 第十八章、数据显示与报告界面的产生 第十九章、图形与可视化 第二十章、Alphanumeric Reporting 第二十一章、流场函数定义 第二十二章、并行处理 第二十三章、自定义函数 第二十四章、参考向导 第二十五章、索引（Bibliography） 第二十六章、命令索引 II、如何使用该教程 概述 本教程主要介绍了FLUENT 的使用，其中附带了相关的算例，从而能够使每一位使用 者在学习的同时积累相关的经验。本教程大致分以下四个部分：第一部分包括介绍信息、用户界面信息、文件输入输出、单位系统、网格、边界条件以及物理特性。第二和第三部分包含物理模型，解以及网格适应的信息。第四部分包括界面的生成、后处理、图形报告、并行处理、自定义函数以及FLUENT 所使用的流场函数与变量的定义。 下面是各章的简略概括 第一部分： z开始使用：本章描述了FLUENT 的计算能力以及它与其它程序的接口。介绍了如何对具体的应用选择适当的解形式，并且概述了问题解决的大致步骤。在本章中，我们给出
了一个可以在你自己计算机上运行的简单的算例。 z使用界面：本章描述了用户界面、文本界面以及在线帮助的使用方法。同时也提供了远程处理与批处理的一些方法。（请参考关于特定的文本界面命令的在线帮助） z读写文件：本章描述了FLUENT 可以读写的文件以及硬拷贝文件。 z单位系统：本章描述了如何使用FLUENT 所提供的标准与自定义单位系统。 z读和操纵网格：本章描述了各种各样的计算网格来源，并解释了如何获取关于网格的诊断信息，以及通过尺度化（scale）、分区（partition）等方法对网格的修改。本章还描述了非一致（nonconformal）网格的使用. z边界条件：本章描述了FLUENT 所提供的各种类型边界条件，如何使用它们，如何定义它们and how to define boundary profiles and volumetric sources. z物理特性：本章描述了如何定义流体的物理特性与方程。FLUENT 采用这些信息来处理你的输入信息。 第二部分： z基本物理模型：本章描述了FLUENT 计算流体流动和热传导所使用的物理模型（包括自然对流、周期流、热传导、swirling、旋转流、可压流、无粘流以及时间相关流）。以及在使用这些模型时你需要输入的数据，本章也包含了自定义标量的信息。 z湍流模型：本章描述了FLUENT 的湍流模型以及使用条件。 z辐射模型：本章描述了FLUENT 的热辐射模型以及使用条件。 z化学组分输运和反应流：本章描述了化学组分输运和反应流的模型及其使用方法。本章详细的叙述了prePDF 的使用方法。 z污染形成模型：本章描述了NOx 和烟尘的形成的模型，以及这些模型的使用方法。 第三部分： z相变模拟：本章描述了FLUENT 的相变模型及其使用方法。 z离散相变模型：本章描述了FLUENT 的离散相变模型及其使用方法。 z多相流模型：本章描述了FLUENT 的多相流模型及其使用方法。 z Flows in Moving Zones（移动坐标系下的流动）：本章描述了FLUENT 中单一旋转坐标系，多重移动坐标系，以及滑动网格的使用方法。 z Solver 的使用：本章描述了如何使用FLUENT 的解法器（solver）。 z网格适应：本章描述了explains the solution-adaptive mesh refinement feature in FLUENT and how to use it 第四部分： z显示和报告数据界面的创建：本章描述了explains how to create surfaces in the domain on which you can examine FLUENT solution data z图形和可视化：本章描述了检验FLUENT 解的图形工具 z Alphanumeric Reporting：本章描述了如何获取流动、力、表面积分以及其它解的数据。 z流场函数的定义：本章描述了如何定义FLUENT 面板内出现的变量选择下拉菜单中的流动变量，并且告诉我们如何创建自己的自定义流场函数。 z并行处理：本章描述了FLUENT 的并行处理特点以及使用方法 z自定义函数：本章描述了如何通过用户定义边界条件，物理性质函数来形成自己的FLUENT 软件。 如何使用该手册 z根据你对CFD 以及FLUENT 公司的熟悉，你可以通过各种途径使用该手册 对于初学者，建议如下：
1. FlUENT 1.1 求解方面 1.1.1 floating point error是什么意思？怎样避免它？ Floating point error已经提过很多次了并且也已经对它讨论了许多。下面是在Fluent论坛上的一些答案： 从数值计算方面看，计算机所执行的运算在计算机内是以浮点数（floating point number）来表示的。那些由于用户的非法数值计算或者所用计算机的限制所引起的错误称为floating point error。 1）非法运算：最简单的例子是使用Newton Raphson方法来求解f(x)=0的根时，如果执行第N次迭代时有，x=x(N)，f’(x(N))=0，那么根据公式x(N+1)=x(N)-f(x(N))/ f’(x(N))进行下一次迭代时就会出现被0除的错误。 2）上溢或下溢：这种错误是数据太大或太小造成的，数据太大称为上溢，太小称为下溢。这样的数据在计算机中不能被处理器的算术运算单元进行计算。 3）舍入错误：当对数据进行舍入时，一些重的数字会被丢失并且不可再恢复。例如，如果对0.1进行舍入取整，得到的值为0，如果再对它又进行计算就会导致错误。 避免方法 计算和迭代我认为设一个比较小的时间步长会比较好的。或者改成小的欠松驰因子也会比较好。从我的经验来看，我把欠松驰因子设为默认值的1/3；降低欠松驰因子或使用耦合隐式求解；改变欠松驰因子，如果是非稳态问题可能是时间步长太大；改善solver-control-limits 比例或许会有帮助；你需要降低Courant数；如果仍然有错误，不选择compute from初始化求解域，然后单击init。再选择你想从哪个面初始化并迭代，这样应该会起作用。另外一个原因可能是courant数太大，就样就是说两次迭代之间的时间步太大并且计算结果变化也较大（残差高）。 网格问题当我开始缩放网格时就会发生这个错误。在Gambit中，所有的尺寸都是以mm 为单位，在fluent按scale按钮把它转换成m，然后迭代几百次时就会发生这种错误。但是当我不把网格缩放到m时，让它和在Gambit中一样，迭代就会成功；我认为你应当检查网格，你的网格数太多了，使用较少的网格问题就会解决；网格太多，计算机资源不够用，使使比较粗的网格。 边界条件在我的分析中，我设了一个wall边界条件来代迭axis边界条件，结果fluent拒绝计算并告诉我floating point error。你的边界条件不能代表真实的物理现象；错误的边界条件定义可能会导致floating point error。例如把内边界设成interior；一次我使用对称边界条件模拟2D区间时也遇到这种问题，我把symmetry设为axe symmetric，就发生了floating point error；检查你设的湍流参数，减小湍流强度，先进行50次迭代。 多处理器问题我近来在进行多处理器模拟时也遇到相似的问题。问题的解决方法是在单个处理器上运行，这样就运算得很好。 错误迭代以错误的条件来初始化，在开始迭代时就会发生floating point error。 1.1.2 coupled和segregated求解有什么区别？ Coupled会同时求解所有的方程（质量守恒方程、动量守恒方程和能量守恒方程）而不是单个方程求解（方程互相分离）。当速度和压力高度耦合（高压和高速）时应该使用耦合求解，但这样会需要较长的计算时间。 在耦合求解中，能量方程中总是包含组分扩散（Species Diffusion Term）项。
Tutorial 7. Modeling Flow Through Porous Media Introduction Many industrial applications involve the modeling of ow through porous media, such as _lters, catalyst beds, and packing. This tutorial illustrates how to set up and solve a problem involving gas ow through porous media. The industrial problem solved here involves gas ow through a catalytic converter. Catalytic converters are commonly used to purify emissions from gasoline and diesel engines by converting environmentally hazardous exhaust emissions to acceptable substances. Examples of such emissions include carbon monoxide (CO), nitrogen oxides (NOx), and unburned hydrocarbon fuels. These exhaust gas emissions are forced through a substrate, which is a ceramic structure coated with a metal catalyst such as platinum or palladium. The nature of the exhaust gas ow is a very important factor in determining the performance of the catalytic converter. Of particular importance is the pressure gradient and velocity distribution through the substrate. Hence CFD analysis is used to designe_cient catalytic converters: by modeling the exhaust gas ow, the pressure drop and the uniformity of ow through the substrate can be determined. In this tutorial, FLUENT is used to model the ow of nitrogen gas through a catalytic converter geometry, so that the ow _eld structure may be analyzed. This tutorial demonstrates how to do the following: _ Set up a porous zone for the substrate with appropriate resistances. _ Calculate a solution for gas ow through the catalytic converter using the pressurebased solver. _ Plot pressure and velocity distribution on speci_ed planes of the geometry. _ Determine the pressure drop through the substrate and the degree of non-uniformity of ow through cross sections of the geometry using X-Y plots and numerical reports. 许多工业应用都涉及通过多孔介质（如过滤器，催化剂床和填料）的流动模型。本教程说明如何建立和解决涉及气体通过多孔介质的问题。 这里解决的工业问题涉及通过催化转换器的气体流量。催化转化器通常用于通过将对环境有害的废气排放物转化为可接受的物质来净化汽油和柴油发动机的排放物。 这种排放的例子包括一氧化碳（CO），氮氧化物（NOx）和未燃烧的碳氢化合物燃料。这些废气排放物被迫通过衬底，该衬底是涂覆有诸如铂或钯的金属催化剂的瓷结构。 排气流量的性质是决定催化转化器性能的一个非常重要的因素。特别重要的是通过基底的压力梯度和速度分布。因此，使用CFD分析来设计催化转换器：通过对排气流量进行建模，可以确定通过基板的流量的压降和流量的均匀性。在本教程中，FLUENT 用于模拟通过催化转化器几何形状的氮气流量，从而可以分析流量结构。 本教程演示了如何执行以下操作： _设置具有适当阻力的基材的多孔区域。 _使用基于压力的解算器计算通过催化转化器的气体流量的解决方案。 _绘制几何体特定平面上的压力和速度分布。 _确定通过基材的压降和不均匀的程度 通过使用X-Y图和数字报告的几何横截面的流量。 Prerequisites This tutorial assumes that you are familiar with the menu structure in FLUENT and that you have completed Tutorial 1. Some steps in the setup and solution procedure will not be shown explicitly. 本教程假设您熟悉FLUENT中的菜单结构您已完成教程1.设置和解决方案过程中的某些步骤不会明确显示。
FLUENT中文手册（简化版） 本手册介绍FLUENT的使用方法，并附带了相关的算例。下面是本教程各部分各章节的简略概括。 第一部分： ?开始使用：描述了FLUENT的计算能力以及它与其它程序的接口。介绍了如何对具体的应用选择适当的解形式，并且概述了问题解决的大致步骤。在本章中给出了一个简单的算例。 ?使用界面：描述用户界面、文本界面以及在线帮助的使用方法，还有远程处理与批处理的一些方法。?读写文件：描述了FLUENT可以读写的文件以及硬拷贝文件。 ?单位系统：描述了如何使用FLUENT所提供的标准与自定义单位系统。 ?使用网格：描述了各种计算网格来源，并解释了如何获取关于网格的诊断信息，以及通过尺度化（scale）、分区（partition）等方法对网格的修改。还描述了非一致（nonconformal）网格的使用. ?边界条件：描述了FLUENT所提供的各种类型边界条件和源项，如何使用它们，如何定义它们等 ?物理特性：描述了如何定义流体的物理特性与方程。FLUENT采用这些信息来处理你的输入信息。 第二部分： ?基本物理模型：描述了计算流动和传热所用的物理模型（包括自然对流、周期流、热传导、swirling、旋转流、可压流、无粘流以及时间相关流）及其使用方法，还有自定义标量的信息。 ?湍流模型：描述了FLUENT的湍流模型以及使用条件。 ?辐射模型：描述了FLUENT的热辐射模型以及使用条件。 ?化学组分输运和反应流：描述了化学组分输运和反应流的模型及其使用方法，并详细叙述了prePDF 的使用方法。 ?污染形成模型：描述了NOx和烟尘的形成的模型，以及这些模型的使用方法。 第三部分： ?相变模拟：描述了FLUENT的相变模型及其使用方法。 ?离散相变模型：描述了FLUENT的离散相变模型及其使用方法。 ?多相流模型：描述了FLUENT的多相流模型及其使用方法。 ?移动坐标系下的流动：描述单一旋转坐标系、多重移动坐标系、以及滑动网格的使用方法。 ?解法器（solver）的使用：描述了如何使用FLUENT的解法器。 ?网格适应：描述了如何优化网格以适应计算需求。 第四部分： ?显示和报告数据界面的创建：本章描述了explains how to create surfaces in the domain on which you can examine FLUENT solution data ?图形和可视化：本章描述了检验FLUENT解的图形工具 ?Alphanumeric Reporting：本章描述了如何获取流动、力、表面积分以及其它解的数据。 ?流场函数的定义：本章描述了如何定义FLUENT面板内出现的变量选择下拉菜单中的流动变量，并且告诉我们如何创建自己的自定义流场函数。 ?并行处理：本章描述了FLUENT的并行处理特点以及使用方法 ?自定义函数：本章描述了如何通过用户定义边界条件，物理性质函数来形成自己的FLUENT软件。 如何使用该手册 对于初学者，建议从阅读“开始”这一章起步。 对于有经验的使用者，有三种不同的方法供你使用该手册：按照特定程序的步骤从按程序顺序排列的目录列表和主题列表中查找相关资料；从命令索引查找特定的面板和文本命令的使用方法；从分类索引查找特定类别信息（在线帮助中没有此类索引，只能在印刷手册中找到它）。 什么时候使用Support Engineer：Support Engineer能帮你计划CFD模拟工程并解决在使用FLUENT 中所遇到的困难。在遇到困难时我们建议你使用Support Engineer。但是在使用之前有以下几个注意事项：●仔细阅读手册中关于你使用并产生问题的命令的信息 ●回忆导致你产生问题的每一步 ●如果可能的话，请记下所出现的错误信息 ●对于特别困难的问题，保存FLUENT出现问题时的日志以及手稿。在解决问题时，它是最好的资源。
Fluent性能分析 仅仅就我接触过得谈谈对fluent的认识，并说说哪些用户适合用，哪些不适合 fluent对我来说最麻烦的不在里面的设置，因为我本身解决的就是高速流动可压缩N-S方程，而且本人也是学力学的，诸如边界条件设置等概念还是非常清楚的同时我接触的流场模拟，都不会有很特别的介质，所以设置起来很简单 对我来说，颇费周折的是gambit做图和生成网格，并不是我不会，而是gambit对作图要求的条件很苛刻，也就是说，稍有不甚，就前功尽弃，当然对于计算流场很简单的用户，这不是问题。有时候好几天生成不了的图形，突然就搞定了，逐渐我也总结了一点经验，就是要注意一些小的拐角地方的图形，有时候做布尔运算在图形吻合的地方，容易产生一些小的面最终将导致无法在此生成网格， fluent里面的计算方法是有限体积法，而且我觉得它在计算过程中为了加快收敛速度，采取了交错网格，这样，计算精度就不会很高。同时由于非结构网格，肯定会导致计算精度的下降，所以我一贯来认为在fluent里面选取复杂的粘性模型和高精度的 格式没有任何意义，除非你的网格做的非常好。 而且fluent5.5以前的版本（包括5。5），其物理模型，（比如粘性流体的几个模型）都是预先设定的，所以，对于那些做探索性或者检验新方法而进行的模拟，就不适合 用。 同时gambit做网格，对于粘性流体，特别是计算湍流尺度，或者做热流计算来说其网格精度一般是不可能满足的，除非是很小的计算区域 所以，用fluent做的比较复杂一点的流场（除了经典的几个基本流场） 其计算所得热流，湍流，以及用雷诺应力模拟的粘性都不可能是准确的， 这在物理上和计算方法已经给fluent判了死刑，有时候看到很多这样讨论的文章，觉得 大家应该从物理和力学的本质上考虑问题。 但是，fluent往往能计算出量级差不多的结果，我曾经做了一个复杂的飞行器热流计算，高超音速流场，得到的壁面热流，居然在量级上是吻合的，但是，从计算热流需要的壁面网格精度来判断，gambit所做的网格比起壁面网格所满足的尺寸的要大了至少2个数量级， 我到现在还不明白fluent是怎么搞的。 综上，我觉得，如果对付老板的一些工程项目，可以用fluent对付过去但是如果真的做论文，或者需要发表文章，除非是做一些技术性工作，比如优化计算 一般用fluent是不适合的。 我感觉fluent做力的计算是很不错的，做流场结构的计算，即使得出一些涡也不是流场本身性质的反应，做低速流场计算，fluent的优势在于收敛 速度快，但是低速流场计算，其大多数的着眼点在于对流场结构的探索，所以计算得到的结果就要好好斟酌一下了，高速流场的模拟中，一般着眼点在于气动力的结果，压力分布
计算流体力学作业FLUENT 模拟燃烧 问题描述：长为2m、直径为0.45m的圆筒形燃烧器结构如图1所示，燃烧筒壁上嵌有三块厚为0.0005 m，高0.05 m的薄板，以利于甲烷与空气的混合。燃烧火焰为湍流扩散火焰。在燃烧器中心有一个直径为0.01 m、长为0.01 m、壁厚为0.002 m的小喷嘴，甲烷以60 m/s的速度从小喷嘴注入燃烧器。空气从喷嘴周围以0.5 m/s的速度进入燃烧器。总当量比大约是0.76(甲烷含量超过空气约28%)，甲烷气体在燃烧器中高速流动，并与低速流动的空气混合，基于甲烷喷嘴直径的雷诺数约为5.7×103。 假定燃料完全燃烧并转换为:CH4+2O2→CO2+2H2O 反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。利用FLUENT的finite-rate化学反应模型对一个圆筒形燃烧器内的甲烷和空气的混合物的流动和燃烧过程进行研究。 1、建立物理模型，选择材料属性，定义带化学组分混合与反应的湍流流动边界条件 2、使用非耦合求解器求解燃烧问题 3、对燃烧组分的比热分别为常量和变量的情况进行计算，并比较其结果 4、利用分布云图检查反应流的计算结果 5、预测热力型和快速型的NO X含量 6、使用场函数计算器进行NO含量计算 一、利用GAMBIT建立计算模型 第1步启动GAMBIT，建立基本结构 分析：圆筒燃烧器是一个轴对称的结构，可简化为二维流动，故只要建立轴对称面上的
二维结构就可以了，几何结构如图2所示。 (1)建立新文件夹 在F盘根目录下建立一个名为combustion的文件夹。 (2)启动GAMBIT (3)创建对称轴 ①创建两端点。A(0,0,0)，B(2,0,0) ②将两端点连成线 (4)创建小喷嘴及空气进口边界 ①创建C、D、E、F、G点
FLUENT多相流模型 分类 1、气液或液液流动 气泡流动：连续流体中存在离散的气泡或液泡 液滴流动：连续相为气相，其它相为液滴 栓塞（泡状）流动：在连续流体中存在尺寸较大的气泡 分层自由流动：由明显的分界面隔开的非混合流体流动。 2、气固两相流动 粒子负载流动：连续气体流动中有离散的固体粒子 气力输运：流动模式依赖，如固体载荷、雷诺数和例子属性等。最典型的模式有沙子的流动，泥浆流，填充床以及各相同性流 流化床：有一个盛有粒子的竖直圆筒构成，气体从一个分散器进入筒内，从床底不断冲入的气体使得颗粒得以悬浮。 3、液固两相流动 泥浆流：流体中的大量颗粒流动。颗粒的stokes数通常小于1。大于1是成为流化了的液固流动。 水力运输：在连续流体中密布着固体颗粒 沉降运动：在有一定高度的盛有液体的容器内，初始时刻均匀散布着颗粒物质，随后，流体会出现分层。 4、三相流 以上各种情况的组合 多相流动系统的实例 气泡流：抽吸、通风、空气泵、气穴、蒸发、浮选、洗刷。 液滴流：抽吸、喷雾、燃烧室、低温泵、干燥机、蒸发、气冷、洗刷。 栓塞流：管道或容器中有大尺度气泡的流动 分层流：分离器中的晃动、核反应装置沸腾和冷凝 粒子负载流：旋风分离器、空气分类器、洗尘器、环境尘埃流动 气力输运：水泥、谷粒和金属粉末的输运 流化床：流化床反应器、循环流化床 泥浆流：泥浆输运、矿物处理 水力输运：矿物处理、生物医学、物理化学中的流体系统 沉降流动：矿物处理。 多相流模型的选择原则 1、基本原则
1)对于体积分数小于10%的气泡、液滴和粒子负载流动，采用离散相 模型。 2)对于离散相混合物或者单独的离散相体积率超出10%的气泡、液滴 和粒子负载流动，采用混合模型或欧拉模型。 3)对于栓塞流、泡状流，采用VOF模型 4)对于分层/自由面流动，采用VOF模型 5)对于气动输运，均匀流动采用混合模型，粒子流采用欧拉模型。 6)对于流化床，采用欧拉模型 7)泥浆和水力输运，采用混合模型或欧拉模型。 8)沉降采用欧拉模型 9)对于更一般的，同时包含多种多相流模式的情况，应根据最感兴趣 的流动特种，选择合适的流动模型。此时由于模型只是对部分流动特 征采用了较好的模拟，其精度必然低于只包含单个模式的流动。 2、混合模型和欧拉模型的选择原则 VOF模型适合于分层的或自由表面流，而混合模型和欧拉模型适合于流动中有相混合或分离，或者分散相的体积分数超过10%的情况（小于10%可使用离散相模型）。 1)如果分散相有宽广的分布（如颗粒的尺寸分布很宽），最好采用混 合模型，反之使用欧拉模型。 2)如果相间曳力规律一直，欧拉模型通常比混合模型更精确；若相间 曳力规律不明确，最好选用混合模型。 3)如果希望减小计算了，最好选用混合模型，它比欧拉模型少解一部 分方程；如果要求精度而不在意计算量，欧拉模型可能是更好的选择。 但是要注意，复杂的欧拉模型比混合模型的稳定性差，可能会遇到收 敛困难。
ANSYS 15.0 系列测试报告 FLUENT 15.0 VOF模型 测试人：崔亮安世亚太公司 测试时间：2013.12.01
1、仿真平台 HP Z820工作站，Intel Xeon E5-2690 * 2，内存64GB，2TB SATA硬盘。安装ANSYS 15.0 Preview3版本。 2、仿真模型 对某车型上带有底部隔板的油箱，在车辆加速时油箱内燃油晃动的瞬态过程进行瞬态仿真分析，网格单元数约10万，使用FLUENT的VOF模型计算空气和燃油的两相交界面。重点考察FLUENT 15.0中VOF模型的计算效率和两相交界面捕捉精度的提升。 测试案例的几何形状 测试案例的网格模型 3、试用情况 1）.稳定性 在整个试用过程中，软件保持稳定，未出现任何不流畅、死机、系统崩溃等情况。2）.流畅度 模型拖动、旋转、缩放等操作十分流畅，模型设定及求解过程操作十分流畅。 3）.效率 该模型使用0.0005秒的时间步长进行瞬态计算，共计算了2000步，共计1.0秒时长。使用15.0 Preview3版本所用的计算时间为3693秒。之前使用13.0版本计算该模型所用计算时间为4381秒。新版本提速15.7%。 4）.硬件资源调用情况 由于该模型网格数量较少，仅使用单核进行求解计算。在整个计算过程中，单核占用率达到100%，内存占用峰值约为400 MB。之前使用13.0版本计算该模型的内存占用峰值约
为450兆。新版本对内存的峰值占用约为旧版本的90%左右。 5）.计算精度 VOF模型的计算精度体现在两相交界面捕捉的清晰程度，15.0版本的交界面捕捉清晰程度比旧版本略有提升，对于一些较小的气泡有着更好的捕捉能力。 t=0.45s时，15.0版本和13.0版本计算的两相交界面对比 t=0.45s时，15.0版本和13.0版本计算的两相交界面对比 4、总结 在ANSYS 15.0 Preview3版本的试用过程中，对FLUENT 15.0中VOF模型的计算效率提升感到满意，相比较于旧版本，约有15%的计算速度提升，这对缩短仿真分析的周期有极大帮助；还有约10%的内存峰值占用量下降，这对于合理利用现有硬件资源进行更大规模的模型计算有着重要意义。此外，新版本VOF模型的计算精度也有所提升，两相交界面捕捉更加锐利，对于一些较小的气泡，相对于旧版本有着更好的捕捉能力
1.多相流动模式 我们可以根据下面的原则对多相流分成四类： ?气-液或者液-液两相流： o 气泡流动：连续流体中的气泡或者液泡。 o 液滴流动：连续气体中的离散流体液滴。 o 活塞流动: 在连续流体中的大的气泡 o 分层自由面流动：由明显的分界面隔开的非混合流体流动。 ?气-固两相流： o 充满粒子的流动：连续气体流动中有离散的固体粒子。 o 气动输运：流动模式依赖诸如固体载荷、雷诺数和粒子属性等因素。最典型的模式有沙子的流动，泥浆流，填充床，以及各向同性流。 o 流化床：由一个盛有粒子的竖直圆筒构成，气体从一个分散器导入筒内。从 床底不断充入的气体使得颗粒得以悬浮。改变气体的流量，就会有气泡不断 的出现并穿过整个容器，从而使得颗粒在床内得到充分混合。 ?液-固两相流 o 泥浆流：流体中的颗粒输运。液-固两相流的基本特征不同于液体中固体颗 粒的流动。在泥浆流中，Stokes 数通常小于1。当Stokes数大于1 时，流动成为流化（fluidization）了的液-固流动。 o 水力运输: 在连续流体中密布着固体颗粒 o 沉降运动: 在有一定高度的成有液体的容器内，初始时刻均匀散布着颗粒物 质。随后，流体将会分层，在容器底部因为颗粒的不断沉降并堆积形成了淤 积层，在顶部出现了澄清层，里面没有颗粒物质，在中间则是沉降层，那里 的粒子仍然在沉降。在澄清层和沉降层中间，是一个清晰可辨的交界面。 ?三相流(上面各种情况的组合) 各流动模式对应的例子如下： ?气泡流例子：抽吸，通风，空气泵，气穴，蒸发，浮选，洗刷 ?液滴流例子：抽吸，喷雾，燃烧室，低温泵，干燥机，蒸发，气冷，刷洗?活塞流例子：管道或容器内有大尺度气泡的流动 ?分层自由面流动例子：分离器中的晃动，核反应装置中的沸腾和冷凝 ?粒子负载流动例子：旋风分离器，空气分类器，洗尘器，环境尘埃流动 ?风力输运例子：水泥、谷粒和金属粉末的输运
Introduction:This tutorial demonstrates how to model the2D turbu-lent?ow across a circular cylinder using LES(Large Eddy Simula-tion),and compute?ow-induced noise(aero-noise)using FLUENT’s acoustics model. In this tutorial you will learn how to: ?Perform2D Large Eddy Simulation(LES) ?Set parameters for an aero-noise calculation ?Save surface pressure data for an aero-noise calculation ?Calculate aero-noise quantities ?Postprocess an aero-noise solution Prerequisites:This tutorial assumes that you are familiar with the menu structure in FLUENT,and that you have solved or read Tu-torial1.Some steps in the setup and solution procedure will not be shown explicitly. Problem Description:The problem considers turbulent air?ow over a2D circular cylinder at a free stream velocity U of69.19m/s. The cylinder diameter D is1.9cm.The Reynolds number based on the?ow parameters is about90000.The computational do-main(Figure3.0.1)extends5D upstream and20D downstream of the cylinder,and5D on both sides of it.If the computational domain is not taken wide enough on the downstream side,so that no reversed?ow occurs,the accuracy of the aero-noise prediction may be a?ected.The rule of thumb is to take at least20D on the downstream side of the obstacle. c Fluent Inc.June20,20023-1
Tutorial:Dam-Break Simulation Using FLUENT’s Volume of Fluid Model Purpose This tutorial examines the dam-break problem using the Volume of Fluid(VOF)multiphase model. This tutorial demonstrates how to do the following: ?Set up a dam-break problem. ?Choose the time step by estimating the maximum possible velocity of the interface and the grid cell dimension. ?Solve the problem using the VOF model. ?Manipulate the solution parameters. Prerequisites This tutorial assumes that you are familiar with the FLUENT interface and that you have a good understanding of basic setup and solution procedures.In this tutorial,you will use VOF multiphase model,so you should have some experience with it.This tutorial will not cover the mechanics of using this model;instead,it will focus on the application of this model to solve a dam-break problem. If you have not used this model before,it would be helpful to?rst refer to the FLUENT6.3 User’s Guide and the FLUENT6.2Tutorial Guide. Problem Description The initial setup of the dam-break problem is shown in Figure1. In this problem,a rectangular column of water,in hydrostatic equilibrium,is con?ned between two walls.Gravity is acting downwards with a magnitude of-9.81m/s2.At the beginning of the calculation,the right wall is removed and the water is allowed to?ow out to the horizontal wall.
1什么叫松弛因子？松弛因子对计算结果有什么样的影响？它对计算的收敛情况又有什 么样的影响？ 1、亚松驰（Under Relaxation）：所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减，以避免由于差值过大而引起非线性迭代过程的发散。用通用变量来写 出时，为松驰因子（Relaxation Factors）。《数值传热学-214》 2、FLUENT中的亚松驰：由于FLUENT所解方程组的非线性，我们有必要控制的变化。一般用亚松驰方法来实现控制，该方法在每一部迭代中减少了的变化量。亚松驰最简 单的形式为：单元内变量等于原来的值加上亚松驰因子a与变化的积, 分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程，包 括耦合解算器所解的非耦合方程（湍流和其他标量）都会有一个相关的亚松驰因子。在FLUENT中，所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题，但是对于一些特殊的非线性问题（如：某些湍流或者高Rayleigh数自然对流问题），在计算开始时要慎重减小亚松驰因子。使用默认的亚松驰因子开始计算是很好的习惯。如 果经过4到5步的迭代残差仍然增长，你就需要减小亚松驰因子。有时候，如果发现残差 开始增加，你可以改变亚松驰因子重新计算。在亚松驰因子过大时通常会出现这种情况。 最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文件，并对解的算法做几 步迭代以调节到新的参数。最典型的情况是，亚松驰因子的增加会使残差有少量的增加， 但是随着解的进行残差的增加又消失了。如果残差变化有几个量级你就需要考虑停止计算 并回到最后保存的较好的数据文件。注意：粘性和密度的亚松驰是在每一次迭代之间的。 而且，如果直接解焓方程而不是温度方程（即：对PDF计算），基于焓的温度的更新是要进行亚松驰的。要查看默认的亚松弛因子的值，你可以在解控制面板点击默认按钮。对于 大多数流动，不需要修改默认亚松弛因子。但是，如果出现不稳定或者发散你就需要减小 默认的亚松弛因子了，其中压力、动量、k和e的亚松弛因子默认值分别为0.2，0.5，0.5和0.5。对于SIMPLEC格式一般不需要减小压力的亚松弛因子。在密度和温度强烈耦合 的问题中，如相当高的Rayleigh数的自然或混合对流流动，应该对温度和/或密度（所用 的亚松弛因子小于1.0）进行亚松弛。相反，当温度和动量方程没有耦合或者耦合较弱时，流动密度是常数，温度的亚松弛因子可以设为1.0。对于其它的标量方程，如漩涡，组分，PDF变量，对于某些问题默认的亚松弛可能过大，尤其是对于初始计算。你可以将松弛因子设为0.8以使得收敛更容易。 SIMPLE与SIMPLEC比较 在FLUENT中，可以使用标准SIMPLE算法和SIMPLEC（SIMPLE-Consistent）算法，默认是SIMPLE算法，但是对于许多问题如果使用SIMPLEC可能会得到更好的结果，尤其是可以应用增加的亚松驰迭代时，具体介绍如下: 对于相对简单的问题（如：没有附加模型激活的层流流动），其收敛性已经被压力速
基于FLUENT 的波浪管道热传递耦合模拟 CFD 可以对热传递耦合的流体流动进行模拟。CFD 模拟可以观察到管道内部的流动行为和热传递，这样可以改进波浪壁面复杂通道几何形状中的热传递。 目的： （1） 创建由足够数量的完整波浪组成的波浪管道，提供充分发展条件； （2） 应用周期性边界条件创建波浪通道的一部分； （3） 研究不同湍流模型以及壁面函数对求解的影响； （4） 采用固定表面温度以及固定表面热流量条件，确定雷诺数与热特性之间的 关系。 问题的描述： 通道由重复部分构成，每一部分由顶部的直面和底部的正弦曲面构成，如图。 图1 管道模型 空气的流动特性如下： 质量流量： m=0.816kg/s; 密度: ρ=1kg/m 3; 动力粘度：μ=0.0001kg/(m ·s); 流动温度： Tb=300K ； 流体其他热特性选择默认项。 流动初试条件： x 方向的速度=0.816m/s ； 湍动能=1m 2/s 2; 湍流耗散率=1×105m 2/s 3。 所有湍流模型中均采用增强壁面处理。 操作过程： 一、 完整波浪管道模型的数值模拟 （1） 计算 Re=uH/v=0.816×1/ (0.0001/1) =8160 Cf/2=0.0359Re -0.2=0.0359× (8160)-0.2=0.0059259 0628.00059259.0816.02 =?==f t C u u y +=u t y/v y=0.00159
（2）创建网格 本例为波浪形管道，管道壁面为我们所感兴趣的地方所以要局部细化。入口和出口处的边界网格设置如图。 图2 边网格 生成面网格 图3 管道网格 （3）运用Fluent进行计算 本例涉及热传递耦合，所以在fluent中启动能量方程，如图。 图4 能量方程