建议大家一起讨论一下湍流边界条件该如何设置
本人也是刚刚接触Fluent系列软件不久，在应用它来求解CFD问题时遇到了不少问题，也得到了很多宝贵经验，其中湍流边界条件的设置就是一个很棘手的问题。
最近对该问题总结经验如下：
在入口、出口或远场边界流入流域的流动，FLUENT需要指定输运标量的值。本节描述了对于特定模型需要哪些量，并且该如何指定它们。也为确定流入边界值最为合适的方法提供了指导方针。
使用轮廓指定湍流参量
在入口处要准确的描述边界层和完全发展的湍流流动，你应该通过实验数据和经验公式创建边界轮廓文件来完美的设定湍流量。如果你有轮廓的分析描述而不是数据点，你也可以用这个分析描述来创建边界轮廓文件，或者创建用户自定义函数来提供入口边界的信息。一旦你创建了轮廓函数，你就可以使用如下的方法：
λSpalart-Allmaras模型：在湍流指定方法下拉菜单中指定湍流粘性比，并在在湍流粘性比之后的下拉菜单中选择适当的轮廓名。通过将m_t/m和密度与分子粘性的适当结合，F LUENT为修改后的湍流粘性计算边界值。
λk-e模型：在湍流指定方法下拉菜单中选择K和Epsilon并在湍动能（Turb. Kinetic E nergy）和湍流扩散速度（Turb. Dissipation Rate）之后的下拉菜单中选择适当的轮廓名。λ雷诺应力模型：在湍流指定方法下拉菜单中选择K和Epsilon并在湍动能（Turb. Kin etic Energy）和湍流扩散速度（Turb. Dissipation Rate）之后的下拉菜单中选择适当的轮廓名。在湍流指定方法下拉菜单中选择雷诺应力部分，并在每一个单独的雷诺应力部分之后的下拉菜单中选择适当的轮廓名。
湍流量的统一说明
在某些情况下流动流入开始时，将边界处的所有湍流量指定为统一值是适当的。比如说，在进入管道的流体，远场边界，甚至完全发展的管流中，湍流量的精确轮廓是未知的。
在大多数湍流流动中，湍流的更高层次产生于边界层而不是流动边界进入流域的地方，因此这就导致了计算结果对流入边界值相对来说不敏感。然而必须注意的是要保证边界值不是非物理边界。非物理边界会导致你的解不准确或者不收敛。对于外部流来说这一特点尤其突出，如果自由流的有效粘性系数具有非物理性的大值，边界层就会找不到了。
你可以在使用轮廓指定湍流量一节中描述的湍流指定方法，来输入同一数值取代轮廓。你也可以选择用更为方便的量来指定湍流量，如湍流强度，湍流粘性比，水力直径以及湍流特征尺度，下面将会对这些内容作一详细叙述。
湍流强度I定义为相对于平均速度u_avg的脉动速度u^'的均方根。
小于或等于1%的湍流强度通常被认为低强度湍流，大于10%被认为是高强度湍流。从外界，测量数据的入口边界，你可以很好的估计湍流强度。例如：如果你模拟风洞试验，自由流的湍流强度通常可以从风洞指标中得到。在现代低湍流风洞中自由流湍流强度通常低到0.0
5%。.
对于内部流动，入口的湍流强度完全依赖于上游流动的历史，如果上游流动没有完全发展或者没有被扰动，你就可以使用低湍流强度。如果流动完全发展，湍流强度可能就达到了百分之几。完全发展的管流的核心的湍流强度可以用下面的经验公式计算：
例如，在雷诺数为50000是湍流强度为4%
湍流尺度l是和携带湍流能量的大涡的尺度有关的物理量。在完全发展的管流中，l被管道的尺寸所限制，因为大涡不能大于管道的尺寸。L和管的物理尺寸之间的计算关系如下：
其中L为管道的相关尺寸。因子0.07是基于完全发展湍流流动混合长度的最大值的，对于非圆形截面的管道，你可以用水力学直径取代L。
如果湍流的产生是由于管道中的障碍物等特征，你最好用该特征长度作为湍流长度L而不是用管道尺寸。
注意：公式并不是适用于所有的情况。它只是在大多数情况下得很好的近似。对于特定流动，选择L和l的原则如下：
λ对于完全发展的内部流动，选择强度和水力学直径指定方法，并在水力学直径流场中指定L=D_H。
λ对于旋转叶片的下游流动，穿孔圆盘等，选择强度和水力学直径指定方法，并在水力学直径流场中指定流动的特征长度为L
λ对于壁面限制的流动，入口流动包含了湍流边界层。选择湍流强度和长度尺度方法并使用边界层厚度d_99来计算湍流长度尺度l，在湍流长度尺度流场中输入l=0.4 d_99这个值
湍流粘性比m_t/m直接与湍流雷诺数成比例(Re_t ?k^2/(e n))。Re_t在高湍流数的边界层，剪切层和完全发展的管流中是较大的(100到1000)。然而，在大多数外流的自由流边界层中m_t/m相当的小。湍流参数的典型设定为1 < m_t/m <10。
要根据湍流粘性比来指定量，你可以选择湍流粘性比（对于Spalart-Allmaras模型）或者强度和粘性比（对于k-e模型或者RSM）。
推导湍流量的关系式
要获得更方便的湍流量的输运值，如：I, L,或者m_t/m，你必须求助于经验公式，下面是F LUENT中常用的几个有用的关系式。要获得修改的湍流粘性，它和湍流强度I长度尺度l
有如下关系：
在Spalart-Allmaras模型中，如果你要选择湍流强度和水力学直径来计算l可以从前面的公式中获得。
湍动能k和湍流强度I之间的关系为：
其中u_avg为平均流动速度
除了为k和e指定具体的值之外，无论你是使用湍流强度和水力学直径，强度和长度尺度或者强度粘性比方法，你都要使用上述公式。
如果你知道湍流长度尺度l你可以使用下面的关系式：
其中是湍流模型中指定的经验常数（近似为0.09），l的公式在前面已经讨论了。
除了为k和e制定具体的值之外，无论你是使用湍流强度和水力学直径还是强度和长度尺度，你都要使用上述公式。
E的值也可以用下式计算，它与湍流粘性比m_t/m以及k有关：
其中是湍流模型中指定的经验常数（近似为0.09）。
除了为k和e制定具体的值之外，无论你是使用湍流强度和水力学直径还是强度和长度尺度，你都要使用上述公式。
如果你是在模拟风洞条件，在风洞中模型被安装在网格和/或金属网格屏下游的测试段，你可以用下面的公式：
其中，是你希望的在穿过流场之后k的衰减(比方说k入口值的10%), 自由流的速度是流域内自由流的流向长度Equation 9是在高雷诺数各向同性湍流中观察到的幂率衰减的线性近似。它是基于衰减湍流中k的精确方程U ?k/?x = - e.
如果你用这种方法估计e，你也要用方程7检查结果的湍流粘性比m_t/m，以保证它不是太大。
虽然这不是FLUENT内部使用的方法，但是你可以用它来推导e的常数自由流值，然后你可以用湍流指定方法下拉菜单中选择K和Epsilon直接指定。在这种情况下，你需要使用方程3从I来计算k。
当使用RSM时，如果你不在雷诺应力指定方法的下拉列表中使用雷诺应力选项，明显的制定入口处的雷诺应力值，它们就会近似的由k的指定值来决定。湍流假定为各向同性，保证
以及
(下标a不求和).
如果你在雷诺应力指定方法下拉列表中选择K或者湍流强度，FLUENT就会使用这种方法。
对大涡模拟（LES）指定入口湍流
大涡模拟模型一节中所描述的LES速度入口中指定的的湍流强度值，被用于随机扰动入口处速度场的瞬时速度。它并不指定被模拟的湍流量。正如大涡模拟模型中介绍的边界条件中所描述的，通过叠加每个速度分量的随机扰动来计算流动入口边界处的随机成分.
压力入口边界条件
压力入口边界条件用于定义流动入口的压力以及其它标量属性。它即可以适用于可压流，也可以用于不可压流。压力入口边界条件可用于压力已知但是流动速度和/或速率未知的情况。这一情况可用于很多实际问题，比如浮力驱动的流动。压力入口边界条件也可用来定义外部或无约束流的自由边界。对于流动边界条件的概述，请参阅流动入口和出口一节。
压力入口边界条件的输入
综述
对于压力入口边界条件你需要输入如下信息
驻点总压
λ驻点总温
λ流动方向
λ静压
λ湍流参数（对于湍流计算）
λ辐射参数(对于使用P-1模型、DTRM模型或者DO模型的计算)
λ化学组分质量百分比(对于组分计算)
λ混合分数和变化(对于PDF燃烧计算)
λ程序变量(对于预混和燃烧计算)
λ离散相边界条件(对于离散相的计算)
λ次要相的体积分数(对于多相计算)
所有的值都在压力入口面板中输入(Figure 1)，该面板是从边界条件打开的。
Figure 1: 压力入口面板
压力输入和静压头
压力场(p_s^')和压力输入(p_s^' or p_0^')包括静压头r_0 g x。也就是FLUENT 以下式定义的压力：
或者
这一定义允许静压头放进体积力项(r - r_0)g中考虑，而且当密度一致时，从压力计算中排除了。因此你的压力输入不因该考虑静压的微分，压力(p^'_s)的报告也不会显示静压的任何影响。有关浮力驱动流动的内容请参阅浮力驱动流动和自然对流的信息
定义总压和总温
在压力入口面板中的Gauge Total Pressure field输入总压值。总温会在Total Temperat ure field中设定。记住，总压值是在操作条件面板中定义的与操作压力有关的的总压值。不可压流体的总压定义为：
对于可压流体为：
其中：p_0 =总压
p_s = 静压
M = 马赫数
c = 比热比(c_p/c_v)
如果模拟轴对称涡流，方程1中的v包括了旋转分量。如果相邻区域是移动的（即：如果使用旋转参考坐标系，多重参考坐标系，混合平面或者滑移网格），而且你是使用分离解算器。那么方程1中的速度（或者方程3中的马赫数）将是绝对的，或者相对与网格速度。这依赖于解算器面板中绝对速度公式是否激活。对于耦合解算器，方程1中的速度（或者方程3中的马赫数）通常是在绝对坐标系下的速度。
看了你的问题我也碰到过，你的初始是如何设置的？
FLUENT的初始化面板中有一项是设置从哪个地方开始计算(compute from), 选择从不同的边界开始计算有很大的区别吗?
一般是选取ALL ZONE，即所有区域的平均处理，通常也可选择有代表性的进口
(如多个进口时)进行初始化。对于一般流动问题，初始值的设定并不重要，
因为计算容易收敛。但当几何条件复杂，而且流动速度高变化快(如音速流动)，
初始条件要仔细选择。如果不收敛，还应试验不同的初始条件，甚至逐次改变边界
条件最后达到所要求的条件。
2 要判断自己模拟的结果是否是正确的，似乎解的收敛性要比那些初始条件和
边界条件更重要，可以这样理解吗？也就是说，对于一个具体的问题，初始条件
和边界条件的设定并不是唯一的，为了使解收敛，需要不断调整初始条件和
边界条件直到解收敛为止，是吗？如果解收敛了，是不是就可以基本确定模拟
的结果是正确的呢？
对于一个具体的问题，边界条件的设定当然是唯一的，只不过初始化时可以选择
不同的初始条件（指定常流），为了使解的收敛比较好，我一般是逐渐的调节边界
条件到额定值（"额定值"是指你题目中要求的入口或出口条件，例如计算一个管
内流动，要求入口压力和温度为10MPa和3000K，那么我开始叠代时选择入口压力和
温度为1MPa和500K（假设，这看你自己问题了），等流场计算的初具规模、收敛的
较好了，再逐渐调高压力和温度，经过好几次调节后最终到达额定值10MPa和3000K，这样比一开始就设为10MPa和3000K收敛的要好些）这样每次叠代可以比较容易收敛，每次调节后不用再初始化即自动调用上次的解为这次的初始解，然后继续叠代。
即使解收敛了，这并不意味着就可以基本确定模拟的结果是正确的，还需要和实验
的结果以及理论分析结果进行对比分析。
怎样判断计算结果是否收敛!
1)观察点处的值不再随计算步骤的增加而变化；
2)各个参数的残差随计算步数的增加而降低，最后趋于平缓；
3)要满足质量守恒（计算中不牵涉到能量）或者是质量与能量守恒（计算中牵涉到能量）。
特别要指出的是，即使前两个判据都已经满足了，也并不表示已经得到合理的收敛解了，因为，如果松弛因子设置得太紧，各参数在每步计算的变化都不是太大，也会使前两个
判据得到满足。此时就要再看第三个判据了。还需要说明的就是,一般我们都希望在收敛的情况下，残差越小越好，但是残差曲线是全场求平均的结果，有时其大小并不一定代表计算结果的好坏，有时即使计算的残差很大，但结果也许是好的，关键是要看计算结果是否符合物理事实，即残差的大小与
模拟的物理现象本身的复杂性有关，必须从实际物理现象上看计算结果。比如说本斑
最近在算的一个全机模型,在大攻角情况下,解震荡得非常厉害，而且残差的量级也总
下不去,但这仍然是正确的,为什么呢,因为大攻角下实际流动情形就是这样的,不断有
涡的周期性脱落,流场本身就是非定常的,所以解也是波动的,处理的时候取平均就可以呢:)
残差
是cell各个Face的通量之和，当收敛后，理论上当单元体内没有源相时各个面流入的通
量也就是对物理量的输运之和应该为0。最大残差或者RSM残差反映流场与所要模拟流场（指收敛后
应该得到的流场，当然收敛后得到的流场与真实流场之间还是存在一定的差距）的差距，残差越小
越好，由于存在数值精度问题，不可能得到0残差，对于单精度计算一般应该低于初始残差1e-03
以下为好，但还要看具体问题。
一般在Fluent里可以添加进出口流量监控（print or plot），当残差收敛到一定程度后，
还要看进出口流量是否达到稳定平衡，才可以确认收敛与否。
残差在较高位震荡，需要检查边界条件是否合理，其次检查初始条件是否合适，比如在有激波的
流场，初始条件不合适，会带来流场的震荡。有时流场可能有分离或者回流，这本身是非定常现象，
计算时残差会在一定程度上发生震荡，这时如果进出口流量是否达到稳定平衡，也可以认为流场收敛了
（前提是要消除其他不合理因数）。另外Fluent缺损地采用多重网格，在计算后期，将多重网格设置
为零可以避免一些波长的残差在细网格上发生震荡。
在进行稳态计算时候，开始残差线是一直下降的，可是到后来各种残差线都显示为波形波动，是不是不收敛阿？
有些复杂或流动环境恶劣情形下确实很难收敛。计算的精度（2阶），网格太疏，网格质量太差，等都会使残差波动。经常遇到，一开始下降，然后出现波动，可以降低松弛系数，我的问题就能收敛，但如果网格质量不好，是很难的。通常，计算非结构网格，如果问题比较复杂，会出现这种情况，建议作网格时多下些功夫。理论上说，残差的震荡是数值迭代在计算域内传递遭遇障碍物反射形成周期震荡导致的结果，与网格亚尺度雷诺数有关。例如，通常压力边界是主要的反射源，换成OUTFLOW边界会好些。这主要根据经验判断。所以我说网格和边界条件是主要因素。
不知道你进行计算时，残差是如何变化的，我个人认为总是出现
turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 2931 cells
reversed flow in 130 faces on outflow 8.
reversed flow in 267 faces on outflow 9.
turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 2838 cells
reversed flow in 115 faces on outflow 8.
reversed flow in 267 faces on outflow 9.
turbulent viscosity limited to viscosity ratio of 1.000000e+005 in 2743 cells
这样的提示结果如果叠代到左后不出现上面提示，结果汉可以接受
reversed flow in 115 faces on outflow 8.
reversed flow in 267 faces on outflow 9.
其实计算有回流现象倒没什么，只是湍流粘度比总是超限结果不可信吧！
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）项。
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多相流模型 分类 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 沉降运动: 在有一定高度的成有液体的容器内，初始时刻均匀散布着颗粒物 质。随后，流体将会分层，在容器底部因为颗粒的不断沉降并堆积形成了淤 积层，在顶部出现了澄清层，里面没有颗粒物质，在中间则是沉降层，那里 的粒子仍然在沉降。在澄清层和沉降层中间，是一个清晰可辨的交界面。 ?三相流(上面各种情况的组合) 各流动模式对应的例子如下： ?气泡流例子：抽吸，通风，空气泵，气穴，蒸发，浮选，洗刷 ?液滴流例子：抽吸，喷雾，燃烧室，低温泵，干燥机，蒸发，气冷，刷洗?活塞流例子：管道或容器内有大尺度气泡的流动 ?分层自由面流动例子：分离器中的晃动，核反应装置中的沸腾和冷凝 ?粒子负载流动例子：旋风分离器，空气分类器，洗尘器，环境尘埃流动 ?风力输运例子：水泥、谷粒和金属粉末的输运
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可能会得到更好的结果，尤其是可以应用增加的亚松驰迭代时，具体介绍如下: 对于相对简单的问题（如：没有附加模型激活的层流流动），其收敛性已经被压力速
两相流：通常把含有大量固体或液体颗粒的气体或液体流动称为两相流；其中含有多种尺寸组颗粒群为一个“相”，气体或液体为另一“相”，由此就有气—液，气—固，液—固等两相流之分。 两相流的研究：对两相流的研究有两种不同的观点：一是把流体作为连续介质，而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。 引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日Langrangian 坐标或物质坐标；以变形后瞬时坐标为自变量称为欧拉Eulerian 坐标或空间坐标。 一.离散相模型 FLUENT在求解连续相的输运方程的同时，在拉格朗日坐标下模拟流场中离散相的第二相；← 离散相模型解决的问题：煤粉燃烧、颗粒分离、喷雾干燥、液体燃料的燃烧等；← 应用范←围：FLUENT中的离散相模型假定第二相体积分数一般说来要小于10-12%（但颗粒质量承载率可以大于10-12%，即可模拟离散相质量流率等/大于连续相的流动）；不适用于模拟在连续相中无限期悬浮的颗粒流问题，包括：搅拌釜、流化床等； 颗粒-颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑；← 湍流中颗粒处理的两种模型：Stochastic Tracking，应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响；Cloud Tracking，运用统计方法来跟踪颗粒围绕某一平均轨道的湍流扩散。通过计算颗粒的系统平均运动方程得到颗粒的某个“平均轨道”← 二.多相流模型 FLUENT中提供的模型： VOF模型(Volume of Fluid Model)← 混合模型(Mixture Model)← 欧拉模型(Eulerian Model)← 1.VOF模型(Volume of Fluid Model) VOF模型用来处理没有相互穿插的多相流问题，在处理两相流中，假设计算的每个控制容积中第一相的体积含量为α1，如果α1=0，表示该控制容积中不含第一相，如果α1=1，则表示该控制容积中只含有第一相，如果0←<α1<1，表示该控制容积中有两相交界面； VOF方法是用体积率函数表示流体自由面的位置和流体所占的体积，其方法占内存小，是一种简单而有效的方法。← 2.混合模型(Mixture Model) 用混合特性参数描述的两相流场的场方程组称为混合模型；← 考虑了界面传递特性以及两相间的扩散作用和脉动作用；使用了滑移速度的概念，允许相以不同的速度运动；← 用于模拟各相有不同速度的多相流；也用于模拟有强烈耦合的各向同性多相流和各相以相同速度运动的多相流；← 缺点：界面特性包括不全，扩散和脉动特性难于处理。← 3.欧拉模型(Eulerian Model) 欧拉模型指的是欧拉—欧拉模型；← 把颗粒和气体看成两种流体，空间各点都有这两种流体各自不同的速度、温度和密度，这些流体其存在在同一空间并相互渗透，但各有不同的体积分数，相互间有滑移；←
3.1计算流体力学基础与FLUENT软件介绍 3.1.1计算流体力学基础 计算流体力学(Computational Fluid Dynamics，简称CFD)是利用数值方法通过计算机求解描述流体运动的数学方程，揭示流体运动的物理规律，研究定常流体运动的空间物理特性和非定常流体运动的时空物理特征的学科[}ss}。其基本思想可以归纳为:把原来在时间域和空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关十这些离散点上场变量之间的关系的代数方程组，然后求解代数方程组获得场变量的近似值[f=}}l 计算流体力学可以看作是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒方程)控制下对流动的数值仿真。通过这种数值仿真，可以得到流场内各个位置上的基本物理量(如速度、压力、温度和浓度等)的分布以及这些物理量随时间的变化规律。 还可计算出相关的其它物理量，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合还可进行结构优化设计等。 过去，流体力学的研究主要有实验研究和理论分析两种方法。实验研究主要以实验为研究手段，得到的结果真实可信，是理论分析和数值计算的基础，其重要性不容低估。然}fu实验往往受到模型尺寸、流场扰动和测量精度等的限制，有时可能难以通过实验的方法得到理想的结果。此外，实验往往经费投入较大、人力和物力耗费较大及周期较长;理论分析方法通常是利用简化的流动模型假设，给出所研究问题的解析解或简化方程。然}fu随着时代的发展，这些方法已不能很好地满足复杂非线性流体运动规律的研究。理论分析方法的优点是所得结果具有普遍适用性，各种影响因素清晰可见，是指导试验研究和验证新的数值计算方法的理论基础。但是，它往往要求对计算对象进行抽象和简化，才有可能得出理论解。}fU对十非线性情况，只有少数流动才能得到解析结果。 计算流体力学方法很好地克服了前面两种方法的弱点，与传统的理论分析方法、实验研究方法一同组成了研究流体流动问题的完整体系。计算流体力学的发展，先后经历 2 FLUENT软件介绍 FLUENT软件是由美国FLUENT公司开发的著名的CFD计算分析软件，在航空、航天、透平机械、汽车、船舶、机械、化工、石化、计算机、半导体、能源、医学等领域得到了广泛的应用。能够解决流动、传热、化学反应、燃烧、多相流、旋涡流动等问题。 FLUENT软件研究的流动模型包括了定常和非定常流动，层流(包括各种非牛顿流模型)，紊流(包括最先进的紊流模型)，不可压缩和可压缩流动，传热和化学反应等。FLUENT软件设计基于“CFD计算机软件群的概念”，针对每一种流动的物理问题的特点，采用适合于它的数值解法在计算速度、稳定性和精度等各方面达到最佳。不同领域的计算软件组合起来，成为CFD软件群，从而高效率地解决各个领域的复杂流动的计算问题，在各软件之间可以方便地进行数值交换，采用统一的前后处理工具，省去了科研工作者在计算方法、编程、前后处理等方面投入的重复、低效的劳动，而可以将主要精力用十物理问题本身的探索上。 流体有限体积法(Finite V olume Method，简称FVM)是目前计算流体动力学领域内应用最普遍的一种对偏微分方程组的离散方法。FLUENT软件就是采用C语言编写的基于非结构化网格和有限体积法的通用CFD求解器，它推出了多种优化的物理模型，如定常和非定常流动;层流(包括各种非牛顿流模型);紊流(包括最先进的紊流模型);不可压缩和可压缩流动;传热;化学反应等。对每一种物理问题的流动特点，有适合它的数值解法，用户可对显式或隐式差分格式进行选择，以期在计算速度、稳定性和精度等方面达到最佳。 在FLUENT 5.0之后的版本中，都采用GAMBIT的专用前处理软件。GAMBIT软件是面向CFD的专业前处理器软件，它包含全面的几何建模能力，也可以从主流的CAD/CAE软件导入几何体和网格，GAMBIT强大的布尔运算能力为建立复杂的几何模型提供的极大的方便。GAMBIT功能强大的网格划分工具，可以划分出包含边界层等CFD特殊要求的高质量的网格。GAMBIT中专有的网格划分算法可以保证在较为复杂的几何区域直接划分出高质量的六面体网格。GAMBIT中的TGRID方法可以在极其复杂的几何区域中划分出与相邻区域网格连续的完全非结构化的网格，GAMBIT网格划分方法的选择完全是智能化的，在选择一个几何区域后GAMBIT会自动选择最合适的网格划分算法，使网格划分过程变得极为容易。 通用CFD软件包，用来模拟从不可压缩到高度可压缩范围内的复杂流动。由于采用了多种求解方法和多重网格加速收敛技术，因而FLUENT能达到最佳的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网格技术及成熟的物理模型，使FLUENT在转捩与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面有广泛应用。
FLUENT软件学习报告 一、软件简介 CFD商业软件FLUENT，是通用CFD软件包，用来模拟从不可压缩到高度可压 FLUENT能达到最佳的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网格技术及成熟的物理模型，使FLUENT在转换与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面有广泛应用。 从本质上讲，FLUENT只是一个求解器。FLUENT本身提供的主要功能包括导入网格模型、提供计算的物理模型、施加边界条件和材料特性、求解和后处理。FLUENT支持的网格生成软件包括GAMBIT、TGRid、prePDF、GeoMesh和其他CAD/CAE软件包。 二、软件使用方法 本学习报告将以一简单算例—台阶运动演示FLUENT软件与GAMBIT及CAD 的结合使用。 2.1 物理模型 二维后台阶运动的计算区域如图2-1所示。计算区域为0.4m×1.2m,台阶长度为0.2m，高度为0.1m。 2.2在CAD中生成几何模型 在CAD中按下列步骤生成如图2-1几何模型：
（1）绘制求解区域形状。 （2）调用PEDIT命令，将构成台阶及边界的线生成多段线。 （3）调用REGION命令，将多段线形成的封闭区间生成区域。 （4）调用EXPORT命令，将绘图结果导出为ASCI格式文件命名为台阶，以便在GAMBIT中进行后续处理。 图2-2是在AUTOCAD中绘制的后台阶绕流的几何模型，该结果包含一个局域。 2.3在GAMBIT中划分网格 在AUTOCAD中生成了一个二维台阶的几何模型，该模型包含一个区域，现在转入到GAMBIT中进行网格划分。 按照导入几何模型、生成流体区域、划分网格、定义边界类型和区域类型的步骤完成GAMBIT划分网格的工作。网格划分完成后输出保存为MSH格式的网格文件。绘制结果如图2-3. 图2-3 网格
建筑群风场fluent计算分析大作业 一、建筑群风场分析目的及计算模型的选取 1、建筑群风场分析的目的 随着城市人口的集中和建筑技术的发展，建筑物之间的间距也变得越来越小，这些建筑物对周围环境风场的影响较大，风力载荷正成为建筑群设计中必须考虑的重要因素。风对建筑物以及建筑物周围环境的影响具体表现为以下几点： （1）在建筑物比较密集的地方，建筑物改变了原来的风场，在相同条件下，建筑物周围的局部风速增大。 （2）风力载荷是一种随机载荷, 受建筑物高度、风向、风的强度以及持续时间的影响很大。高层建筑物周围的局部负压过大, 使得建筑物顶局部掀起或装饰玻璃破碎、脱 落。 （3）建筑物的外轮廓形状一般都是非流线形的, 因而流场不可避免地伴随有分离流动、涡的脱落和振荡.这些现象会在建筑物的居室产生严重的噪音, 更严重时还会引起结 构和流体的耦合震荡。 因此，研究建筑群风场的速度分布、静压分布十分有必要。 2、计算模型选择 当建筑物是钝体，空气绕过钝体时的风场和绕过流线体时存在着分离流和剪切层的非定常振动，钝体周围流场十分复杂，是由撞击、分离、再附、环绕和旋涡等确定的。另外，建筑物通常建造在大气边界层。在大气边界层中，气流质点运动杂乱无章，气流流动表现为湍流状态。湍流是由大小不同尺度的涡体叠合而成，对时间和空间都是非线性的随机运动的，因此使用湍流模型解决此问题。在湍流模型中，基于Reynolds时均的Realizable K-ε模型能
在整体上很好地反映出建筑物表面风压的变化趋势，模拟结果与试验值相差较小且计算效率高，所以选用该模型。 二、计算模型设置 1、边界条件选取 计算流域入流处采用 FLUENT 中的速度进口边界条件（velocity-inlet ）。边界条件用于定义在流动进口处的流动速度及相关其它标量型流动变量。该边界条件适用于不可压缩流动，对于可压缩流动问题时会使得入口处的总温度或总压有一定的波动，导致非物理结果，所以可压缩流问题不适合采用速度进口边界条件。本文为不可压缩流，可采用以 velocity-inlet 边界。需对流动速度 v 、k 和ε定义。设建筑物所在地形为B 类地形，其风场为B 类风场，10m 高度处、10min 平均的基本风压为 w0=0.35kPa ，相应的标准高度处平均风速为u0＝23.7m/s ；则人口处的湍流强度I 、湍流动能K 和湍流耗散率 ε的具体表达式如下 21.5()k u I =?； 3342 0.09k l ε=； 0.250.3150.1()5450450z I z z -≤??=?<≤?? 其中z ，u 分别是流域中任意高度和对应的平均风速，z 由模型底部开始算起。l 为湍流积分尺度，采用经验公式0.5100(l z =。平均风速剖面、湍动能，c 和耗散率 值采用Fluent 提供的UDF 编程与Fluent 作接口实现。 出口采用完全发展出流边界条件（outflow ）。Outflow 边界条件用于出流边界上的压力或速度都未知的情况，适用于出口处流动是完全发展的情况。 计算流域顶部和两侧采用对称边界条件（symmetry ），适用于流动及传热场是对称的
第四章Fluent后处理 利用FLUENT 提供的图形工具可以很方便的观察CFD 求解结果，并得到满意的数据和图形，用来定性或者定量研究整个计算。本章将重点介绍如何使用这些工具来观察您的计算结果。 1 生成基本图形 在FLUENT中能够方便的生成网格图、等值线图、剖面图,速度矢量图和迹线图等图形来观察计算结果。下面将介绍如何产生这些图形。 一、生成网格图 生成网格或轮廓线视图的步骤 (1)打开网格显示面板 菜单：Display –〉Grid... 图4-1 网格显示对话框 (2)在表面列表中选取表面。点击表面列表下的Outline 按钮来选择所有“外”表面。如果所有的外表面都已经处于选中状态，单击该按钮将使所有外表面处于未选中的状态。点击表面列表下的Interior 按钮来选择所有“内”表面。同样，如果所有的内表面都已经处于选中状态，单击该按钮将使所有内表面处于未选中的状态。 (3)根据需要显示的内容，可以选择进行下列步骤： 1)显示所选表面的轮廓线，在图4-1所示的对话框中进行如下设置：在Options 项选择Edges，在Edge Type 中选择Outline。 2)显示网格线,在Options 选择Edges，在Edge Type 中选择ALL。 3)绘制一个网格填充图形,在Options 选择Faces。显示选中面的网格节点,在Options 选择Nodes。
(4)设置网格和轮廓线显示中的其它选项。 (5)单击Display 按钮，就可以在激活的图形窗口中绘制选定的网格和轮廓线。 二、绘制等值线和轮廓图 生成等值线和轮廓的步骤： 通过图4-2 所示的等值线对话框来生成等值线和轮廓。 菜单：Display –〉Contours... 图4-2 等值线对话框 生成等值线或轮廓的基本步骤如下： (1) 在Contours Of 下拉列表框中选择一个变量或函数作为绘制的对象。首先在上面的列表中选择相关分类；然后在下面的列表中选择相关变量。 (2) 在Surfaces 列表中选择待绘制等值线或轮廓的平面。对于2D情况，如果没有选取任何面，则会在整个求解对象上绘制等值线或轮廓。对于3D情况，至少需要选择一个表面。 (3) 在Levels 编辑框中指定轮廓或等值线的数目。最大数为100。 (4) 如果需要生成一个轮廓视图，请在Option 中选中Draw Profiles 选项。在轮廓选项对话框中(如图4-3)，可以如下定义轮廓：