2022年FLUENT知识点解析 .pdf

《2022年FLUENT知识点解析 .pdf》由会员分享，可在线阅读，更多相关《2022年FLUENT知识点解析 .pdf（39页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1 一、基本设置1Double Precision的选择启动设置如图，这里着重说说Double Precision（双精度）复选框，对于大多数情况，单精度求解器已能很好的满足精度要求，且计算量小，这里我们选择单精度。然而对于以下一些特定的问题，使用双精度求解器可能更有利1。a. 几何特征包含某些极端的尺度（如非常长且窄的管道），单精度求解器可能不能足够精确地表达各尺度方向的节点信息。b. 如果几何模型包含多个通过小直径管道相互连接的体，而某一个区域的压力特别大 （因为用户只能设定一个总体的参考压力位置），此时，双精度求解器可能更能体现压差带来的流动（如渐缩渐扩管的无粘与可压缩流动模拟） 。c.

2、对于某些高导热系数比或高宽纵比的网格，使用单精度求解器可能会遇到收敛性不佳或精确度不足不足的问题，此时，使用双精度求解器可能会有所帮助。1 李鹏飞 ,徐敏义 ,王飞飞 .精通 CFD 工程仿真与案例实战： FLUENT GAMBIT ICEM CFD TecplotM. 北京,人民邮电出版社 ,2011:114-116 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 1 页，共 39 页 - - - - - - - - - 2 2网格光顺化用光滑和交换的方式改善网格

3、： 通过 Mesh下的 Smooth/Swap来实现，可用来提高网格质量，一般用于三角形或四边形网格，不过质量提高的效果一般般， 影响较小，网格质量的提高主要还是在网格生成软件里面实现， 所以这里不再用光滑和交换的方式改善网格，其原理可参考 FLUENT 全攻略 （已下载）。3Pressure-based 与 Density-based 求解器设置如图。下面说一说Pressure-based 和 Density-based的区别：Pressure-Based Solver是 Fluent 的优势，它是基于压力法的求解器，使用的是压力修正算法，求解的控制方程是标量形式的，擅长求解不可压缩流动，对

4、于可压流动也可以求解；Fluent 6.3以前的版本求解器，只有Segregated Solver和 Coupled Solver ，其实也是 Pressure-Based Solver的两种处理方法；名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 2 页，共 39 页 - - - - - - - - - 3 Density-Based Solver是 Fluent 6.3新发展出来的，它是基于密度法的求解器，求解的控制方程是矢量形式的，主要离散格式有Roe ，AU

5、SM+，该方法的初衷是让Fluent 具有比较好的求解可压缩流动能力，但目前格式没有添加任何限制器，因此还不太完善；它只有Coupled 的算法；对于低速问题，他们是使用Preconditioning方法来处理，使之也能够计算低速问题。Density-Based Solver下肯定是没有 SIMPLEC ，PISO这些选项的，因为这些都是压力修正算法，不会在这种类型的求解器中出现的；一般还是使用Pressure-Based Solver 解决问题。基于压力的求解器适用于求解不可压缩和中等程度的可压缩流体的流动问题。而基于密度的求解器最初用于高速可压缩流动问题的求解。虽然目前两种求解器都适用于各

6、类流动问题的求解（从不可压缩流动到高度可压缩流动） ，但对于高速可压缩流动而言，使用基于密度的求解器通常能获得比基于压力的求解器更为精确的结果。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 3 页，共 39 页 - - - - - - - - - 4 4axisymmetric和 axisymmetric swirl 从字面的意思很好理解axisymmetric和 axisymmetric swirl的差别：axisymmetric ：是轴对称的意思，也就是关于一

7、个坐标轴对称，2D的 axisymmetric问题仍为 2D问题。而 axisymmetric swirl：是轴对称旋转的意思，就是一个区域关于一条坐标轴回转所产生的区域，这产生的将是一个回转体， 是 3D的问题。在 Fluent 中使用这个，是将一个 3D的问题简化为 2D问题，减少计算量，需要注意的是，在Fluent 中，回转轴必须是x 轴。5操作工况参数（ Operating Conditions）操作压力的介绍名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第

8、4 页，共 39 页 - - - - - - - - - 5 关于参考压力的设定，首先需了解有关压力的一些定义。ANSYS FLUENT 中有以下几个压力，即Static Pressure（静压） 、Dynamic Pressure （动压）与 Total Pressure（总压） ；Absolute Pressure（绝对压力）、 Relative Pressure（参考压力）与 Operating Pressure（操作压力）。这些压力间的关系为， Total Pressure（总压）=Static Pressure（静压） +Dynamic Pressure （动压）；Absolute

9、Pressure（绝对压力）=Operating Pressure（操作压力） +Gauge Pressure （表压） 。其中，静压、动压和总压是流体力学中关于压力的概念。静压是测量到的压力，动压是有关速度动能的压力， 是流动速度能量的体现。而绝对压力、操作压力和表压是FLUENT 引入的压力参考量，在ANSYS FLUENT中，所有设定的压力都默认为表压。这是考虑到计算精度的问题。操作压力的设定设定操作压力时需要注意的事项如下：对于不可压缩理想气体的流动， 操作压力的设定直接影响流体密度的计算， 因为对于理想气体而言， 流动的密度由理想气体方程获得，理想气体方程中的压力为操作压力。对于低马

10、赫数的可压缩流动而言，相比绝对静压， 总压降是很小的，因此其计算精度很容易受到数值截断误差的影响。需要采取措施来避免此误差的形成，ANSYS FLUENT通过采用表压名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 5 页，共 39 页 - - - - - - - - - 6 （由绝对压力减去操作压力）的形式来避免截断误差的形成，操作压力一般等于流场中的平均总压。对于高马赫数可压缩流动的求解而言，因为此时的压力比低马赫可压缩流动的大得多， 所以求解过程中的截断误差的影

11、响不大，可以不设定表压。由于ANSYS FLUENT中所有需输入的压力都为表压，因此此时可以将操作压力设定为0（这样可以最小化由于压力脉动而引起的误差） ，使表压与绝对压力相等。如果密度设定为常数或者其值由通过温度变化的函数获得，操作压力并没有在计算密度的过程中被使用。默认的操作压力为101325Pa 。操作压力的设定主要基于两点考虑，一是流动马赫数的大小，二是密度计算方法。表格 1 操作压力的推荐设置密度关系式马赫数操作压力理想气体定律大于 0.1 0 或约等于流场的平均压力理想气体定律小于 0.1 约等于流场的平均压力关于温度的函数不可压缩不使用常数不可压缩不使用不可压缩的理想气体不可压缩

12、约等于流场的平均压力名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 6 页，共 39 页 - - - - - - - - - 7 关于参考压力位置的设定对于不涉及任何压力边界条件的不可压缩流动，ANSYS FLUENT在每次迭代后要调整表压值。 这个过程通过使用参考压力位置处（或该位置附近）节点的压力完成。因此，参考压力位置处的表压应一直为0。如果使用了压力边界条件，则不会使用到上述关系，因此参考压力位置不被使用。参考压力位置默认为等于或接近（0，0，0）的节点中心

13、位置。实际计算中可能需要设置参考压力位置到绝对静压已知的位置处。在Operating Conditions对话框中的 Reference Pressure Location选项组中设置新的参考压力位置的x，y，z 的坐标即可。如果要考虑某一方向的加速度，如重力，可以勾选Gravity复选框。对于 VOF计算，应当选择 Specified Operating Density，并且在 Operating Density 下为最轻相设置密度。这样做排除了水力静压的积累，提高了round-off精度为动量平衡。同样需要打开Implicit Body Force，部分平衡压力梯度和动量方程中体积力，提高

14、解的收敛性。Reference Pressure Location （参考压强位置）应是位于流体永远是 100% 的某一相（空气）的区域，光滑和快速收敛是其基本条件。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 7 页，共 39 页 - - - - - - - - - 8 二、求解模型的设定1流动模型的设置无粘模型理想流体是一种设想的没有粘性的流体，在流动时各层之间没有相互作用的切应力，即没有内摩擦力。十分明显，理想流体对于切向变形没有任何抗拒能力。 应该强调指出

15、， 真正的理想流体在客观实际中是不存在的，它只是实际流体在某些条件下的一种近似模型。在 Inviscid流动模型应用方面， 无粘流动忽略了粘性对流动的影响，这对高雷诺数的流动是合适的，因为高雷诺数流动惯性力的作用远大于粘性力的作用， 粘性力可以忽略， 所以可以将其考虑成无粘流名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 8 页，共 39 页 - - - - - - - - - 9 动。无粘流动的求解更快， 其激波在某些值上预测的偏高。无粘流动能对流动状态和激波位置

16、进行快速预测。马赫数与激波马赫数的定义是vMa它表示流体的流动速度与当地声速之比，是一个无量纲的参量。对应于1M，1M和1M这三种情况的流动分别称为亚声速流、声速流和超声速流。 当马赫数很小时， 速度的相对变化只能引起很小的密度相对变化，但当马赫数很大时，则将引起较大的密度相对变化，这也说明了马赫数是流体压缩性的一个表征。当飞机、炮弹和火箭以超音速飞行时，或者发生强爆炸、强爆震时，气流受到急剧的压缩， 压强和密度突然显著增加，这时所产生的压强扰度将比声速大得多的速度传播，波阵面所到之处气流的各种参数都将发生显著变化，参数突跃。这样一个强间断面叫做激波阵面。渐缩渐扩管的流动是计算流体力学模拟的经

17、典问题之一。在这类流动中，激波的出现是流动中可压缩效应的体现。精确的激波模拟是CFD研究的热点之一。为了更好捕捉压力梯度，需要采用较细的网格并结合合适的数值模拟和格式。很多实际模拟中， 局部网格的自适应会很有帮助。层流模型名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 9 页，共 39 页 - - - - - - - - - 10 流动有层流和湍流之分，判断湍流的标准可以参考2，这里写出内流的判断标准：Re2300UD对于内流而言， 一般大多数流动都是湍流， 一般不

18、使用湍流模型。而对一些外流而言（如外掠平板或是外掠障碍物），则很有可能是层流运动。湍流模型的评价与选择a.k湍流模型这里我们使用的湍流模型是Standard k模型，这种模型应用较多，计算量适中，有较多数据积累和比较高的精度，对于曲率较大和压力梯度较强等复杂流动模拟效果欠佳。一般工程计算都使用该模型，其收敛性和计算精度能满足一般的工程计算要求，但模拟旋流和绕流时有缺陷。壁面函数的选择对于有壁面的流动，当主流为充分发展湍流时，根据离壁面法线距离不同，可将流动划分为壁面区 （或称内区、近壁区）和核心区（或称外区） 。核心区是完全湍流区，为充分发展的湍流。在壁面区，由于有壁面的影响，流动与核心区不同

19、。壁面区可分为 3 个子层：粘性底层、过渡层和对数率层。2 李鹏飞 ,徐敏义 ,王飞飞 .精通 CFD 工程仿真与案例实战： FLUENT GAMBIT ICEM CFD TecplotM. 北京,人民邮电出版社 ,2011:122 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 10 页，共 39 页 - - - - - - - - - 11 粘性底层是一个紧贴壁面的极薄层，在动量、热量和质量的交换过程中粘性力起主要作用， 而湍流切应力可以忽略， 因此流动几乎可以

20、看成层流流动，且在平行于壁面方向上的速度呈线性分布。过渡层处于粘性底层之外，在此层中，粘性力和湍流切应力的作用相当，流动状况较为复杂，很难用公式或定律表述。实际工程计算中由于过渡层厚度极小， 可不考虑此层，直接以对数率层的方法处理。对数率层处于近壁区的最外层，粘性力的影响不明显，湍流切应力占主要地位，流动处于充分发展的湍流状态， 流速分布接近对数律。壁面区内不同子层的高度和速度可以沿壁面法向的无量纲高度和无量纲速度表达。uuUyUy其中，u是流体的时均速度，U是壁面摩擦速度，wU，w是壁面切应力，y是壁面的垂直距离。在5y时，区域为粘性底层， 此时速度沿壁面法线方向呈线性分布，即uy。在603

21、00y时，流动处于对数率层，此时速度沿壁面法线方向呈对数率分布，即2.5ln5.5uy。壁面函数法的本质是， 对于湍流核心区的流动使用k模型求解，而在壁面区并不进行求解， 直接使用半经验公式得出该区域的速度等物理量。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 11 页，共 39 页 - - - - - - - - - 12 FLUENT 提供了多种壁面函数处理方式，如标准壁面函数法、非平衡壁面函数法和增强壁面处理。标准壁面函数法利用对数校正法提供了必需的壁面边界

22、条件（对于平衡湍流边界层）。而非平衡壁面函数法用来改善高压力梯度、分离、再附和滞止等情况下的结果。标准壁面函数法和非平衡壁面函数法都允许在近避免区域上使用较粗的网格。对于大多数高雷诺数情况使用标准的或者非平衡的壁面函数（6Re10） 。增强壁面处理选项把混合边界模型和两层边界模型结合起来，对低雷诺数流动或者复杂近壁面现象很适合，湍流模型在内层上得到了修正。表格 2 几种壁面处理方法的比较优点缺点标准壁面函数法应用较多，计算量小，有较高的精度适合高雷诺数流动， 对低雷诺数流动问题，有压力梯度、高度蒸腾和大的体积力、 低雷诺数和高速三维流动问题不适合非平衡壁面函数法考虑了压力梯度， 可以计算分离，

23、 在附着以及撞击问题对低雷诺数流动问题， 有较强压力梯度、强体积力及强三维性问题不适合增强壁面处理不依赖壁面法则， 对于复杂流动， 特别是要求网格密， 因而要求计算机处理时间长，内存大名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 12 页，共 39 页 - - - - - - - - - 13 低雷诺数流动很适合2多相流模型VOF模型该模型通过求解单独的动量方程和处理穿过区域的每一流体的容积比来模拟两种或三种不能混合的流体。典型的应用包括流体喷射、流体中大泡运动、

24、 流体在大坝坝口的流动、 气液界面的稳态和瞬态处名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 13 页，共 39 页 - - - - - - - - - 14 理等。一般而言 VOF主要适用于非稳态的多相流模型，仅对某些特定问题的多相流模型的稳态问题能够适用。VOF方法适用于计算空气和水这样不能互相掺混的流体流动，对于分层流和活塞流，最方便的就是选择VOF模型。需要注意的是，对于湍流模型的设置， VOF不能用于无粘流，也不能用大涡模拟3。Geo-Reconstru

25、ct 格式Geo-Reconstruct 格式（在 Solution Methods中设置）是一种较为精确的追踪自由表面的计算格式，广泛地应用于瞬变流的VOF问题中，但必须注意的要使用该格式VOF模型必须使用显示离散格式 （在VOF模型设置选项设置）。Body Force Formulation 为提高解的收敛性，对于涉及到表面张力的计算，建议在Body Force Formulation 中勾选 Implicit Body Force。这样做由于压力梯度和动量方程中表面张力的部分平衡，从而提高解的收敛性。3 李进良 , 李承曦 , 胡仁喜 . 精通 FLUENT.6.3 流场分析 M. 北京

26、 , 化学工业出版社 , 2009:231-236 名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 14 页，共 39 页 - - - - - - - - - 15 Mixture 模型这是一种简化的多相流模型， 用于模拟各种有不同速度的多相流，但是假定了在短空间尺度上局部的平衡。相之间的耦合应当是很强的。它也用于模拟有强烈耦合的各向同性多相流和各向以相同速度运动的多相流。典型的应用包括沉降（sedimentation ） 、气旋分离器、低载荷作业下的多粒子流动、气

27、相容积率很低的泡状流。Mixture Parameters 一般需要勾选 Mixture Parameters中的 Slip Velocity复选框，以此来求解滑移速度模型， 因为在多相流中各种组分的速度有很大不名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 15 页，共 39 页 - - - - - - - - - 16 同。对于求解一个均匀的多相流问题可以选择不做滑移速度的计算，可以在 mixture parameters选项下将 slip velocity关掉

28、。Eulerian模型该模型可以模拟多相分离流及相互作用的相，相可以是液体、气体、固体。与在离散相模型中Eulerian-Lagrangian方案只用于离散相不同，在多相流模型中Eulerian方案用于模型中的每一项。3固化与熔化模型FLUENT 采用“焓多孔度（ enthalpy-porosity） ”技术模拟流体的固化和熔化（ Solidification/Melting）过程。在流体的固化和熔化问题中，流场可以分成流体区域、 固体区域和两者之间的糊状区域。“焓多孔度”技术采用的计算策略是将流体在网格单元内占有的体名师归纳总结 精品学习资料 - - - - - - - - - - - -

29、- - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 16 页，共 39 页 - - - - - - - - - 17 积百分比定义为多孔度（porosity ） ，并将流体和固体并存的糊状区域看作多孔介质区进行处理。在流体的固化过程中，多孔度从1 降低到 0；反之，在熔化过程中，多孔度则从0 升至 1。 “焓多孔度”技术通过在动量方程中添加汇项（即负的源项） 模拟因固体材料存在而出现的压强降。“焓多孔度”技术可以模拟的问题包括纯金属或二元合金中的固化、熔化问题、 连续铸造加工过程等。计算中可以计算固体材料与壁面之间因空气的存在而产生的热阻，固化

30、、熔化过程中组元的输运等等。需要注意的是，在求解固化、熔化问题的过程中，只能采用分离算法，只能与 VOF模型配合使用，不能计算可压缩流，不能单独设定固体材料和流体材料的性质，同时在模拟带反应的组元输运过程时，无法将反应区限制在流体区域，而是在全流场进行反应计算。Parameters 定义在 Parameters 下面定义 Mushy Zone Constant （糊状区域常数）。这个常数的取值范围一般在104到 107之间，取值越大沉降曲线就越陡峭， 固化过程的计算速度就越快， 但是取值过大容易引起计算振荡，因此需要在计算中通过试算获得最佳数值。Materials设置在 Materials （

31、材料） 面板上， 定义 Melting Heat （熔化热）、 Solidus Temperature（固相点温度） 和 Liquidus Temperature（液相点温度）。如果计算中涉及组元输运过程，则必须同时定义溶剂的融解温度（Melting Temperature） ，同时需要定义熔化物的液相线相对于浓度名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 17 页，共 39 页 - - - - - - - - - 18 的斜率（ Slope of Liquid

32、us Line） 、分配系数（ Partition Coefficient）和固体中的扩散速率（Diffusion in Solid ）等参数。设置边界条件除了常规的边界条件设置，对于固化和熔化问题还有一些特殊设置，其中包括：在计算壁面接触热阻时设置接触热阻（Contact Resistance ） 。 这个参数在 Wall （壁面）面板中的 Thermal Conditions（热力学条件）下给定。 Shear Condition （剪切条件）下选择Marangoni Stress （Marangoni 应力）选项。如果计算拉出速度，则在边界条件中的速度边界条件将被用于拉出速度的计算。名师归

33、纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 18 页，共 39 页 - - - - - - - - - 19 三、相设置相设置一般用于多相流的设置，对于相设置，这里主要讲一下Interaction的设置，如图：Interaction设置Interaction设置用来定义两相的相互作用，其有多个选项卡，如图。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - -

34、 - - 第 19 页，共 39 页 - - - - - - - - - 20 Drag 选项卡针对每对物相，在下拉菜单中选择阻力函数。其中包括schiller-naumann 模型、morsi-alexander 模型、symmetric（对称）模型等用于流体与流体之间阻力计算的模型，也包括wen-yu 模型、gidaspow 模型、syamlal-obrien 模型等用于液体与固体之间阻力计算的模型，还包括 syamlal-obrien-symmetric 模型用于固体与固体之间的阻力计算。 除此之外，还可以将阻力函数定义为constant （常数） ，或者选择 user-defined

35、（用户定义） 由用户自己定义阻力函数。如果计算中不需要设定阻力，还可以选择none（不计阻力）选项。阻力设置的相关原理比较复杂，可参考帮助，一般保持默认的schiller-naumann设置不变。Surface Tension选项卡Surface Tension选项卡用来定义表面张力，如果相包含壁面粘附，可勾选“ Wall Adhesion ”复选框。四、Cell Zone Condition Frame Motion 选项对于流体，可以通过Frame Motion 选项确定坐标运动方式（如离心泵内部流体的旋转使用运动参考系模型），如图：名师归纳总结 精品学习资料 - - - - - - -

36、- - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 20 页，共 39 页 - - - - - - - - - 21 Rotation-Axis Origin指转轴， Rotational Velocity值旋转速度。五、边界条件设置1速度进口边界条件（ Velocity Inlet）名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 21 页，共 39 页 - - - - - - - - - 22

37、速度进口边界条件用进口处流场速度及相关流动变量作为边界条件，在速度进口边界条件中，流场进口边界的驻点参数是不固定的。为了满足进口处的速度条件，驻点参数将在一定范围内波动。需要注意的是，因为这种条件中允许驻点参数浮动，所以速度进口边界条件仅适用于不可压流，如果用于可压流， 则可能导致出现非物理解。同时还要注意的是， 不要让速度进口条件过于靠近进口内侧的固体障碍物，这样会使驻点参数的不均匀程度大大增加。湍流参数的设置对于一般的流动边界条件，均涉及到湍流参数的定义，在Turbulence Specification Method （湍流定义方法）下拉列表中，可以简单地用一个常数来定义湍流参数，即通过

38、给定湍流强度、 湍流粘度比、水力直径或湍流特征长在边界上的值来定义流场边界上的湍流。a. 湍流强度（ Turbulence Intensity）湍流强度定义如下：222uvwIu上式中u，v，w是速度脉动量，u是平均速度。比较常用的是 Intensity and Hydraulic Diameter，湍流强度与水力直径的确定有相应的计算方法，这里只是采用估算来加以确定。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 22 页，共 39 页 - - - - - - -

39、 - - 23 内流问题进口处的湍流强度取决于上游流动状态。如果上游是没有充分发展的未受扰流动， 则进口处可以使用低湍流强度。如果上游是充分发展的湍流， 则进口处湍流强度可以达到几个百分点。如果管道中的流动是充分发展的湍流，则湍流强度I（turbulence intensity）可以用下面公式计算得到，这个公式是从管流经验公式得到的：1/80.16(Re)HDuIu其中：ReHD为按等效水力直径HD计算得到的雷诺数。湍流强度小于 1时，可以认为湍流强度是比较低的，而在湍流强度大于 10时，则可以认为湍流强度是比较高的。比如，当雷诺数为 50000时，代入上述公式可得到湍流强度为4% ，默认的湍

40、流强度设置为 5% （中等强度）。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 23 页，共 39 页 - - - - - - - - - 24 b. 湍流的长度尺度与水力直径湍流的长度尺度（ Turbulence Length Scale）与水力直径（Hydraulic Diameter）是设置湍流的重要参数。湍流能量主要集中在大涡结构中， 而湍流长度尺度l则是与大涡结构相关的物理量。 在充分发展的管流中， 因为漩涡尺度不可能大于管道直径，所以l是受到管道尺寸制

41、约的几何。 湍流长度尺度l与管道物理尺寸L关系可以表示为：0.07lL式中的比例因子 0.07 是充分发展管流中混合长的最大值，而L则是管道直径。在管道截面不是圆形时，L可以取为管道的水力直径。水力直径是在管内流动（ internal pipe flow ）中引入的，其目的是为了给非圆管流动取一个合适的特征长度来计算其雷诺数。四倍的湿横截面面积与湿圆周长度之商。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 24 页，共 39 页 - - - - - - - - -

42、 25 湍流的特征长取决于对湍流发展具有决定性影响的几何尺度。在上面的讨论中， 管道直径是决定湍流发展过程的唯一长度量。如果在流动中还存在其他对流动影响更大的物体，比如在管道中存在一个障碍物，而障碍物对湍流的发生和发展过程起着重要的干扰作用。在这种情况下，湍流特征长就应该取为障碍物的特征长度。从上面的分析可知，虽然上式对于大多数管道流动是适用的，但并不是普遍适用的，在某些情况下可以进行调整。在 FLUENT 中选择特征长L或湍流长度尺度l的方法如下：1）对于充分发展的内流，可以用Intensity and Hydraulic Diameter （湍流强度与水力直径）方法定义湍流，其中湍流特征长

43、度就是Hydraulic Diameter（水力直径）HD。2）对于导向叶片或分流板下游的流场，可以用Intensity and Hydraulic Diameter （湍流强度与水力直径） 定义湍流，并在 Hydrauli Diameter（水力直径）中将导向叶片或分流板的开口部分的长度L定义为特征长度。3）如果进口处的流动为受到壁面限制且带有湍流边界层的流动，可以在 Intensity and Length Scale 面板中用边界层厚度99通过公式990.4l计算得到湍流长度尺度l。 最后在 Turbulence Length Scale（湍流长度尺度）中输入l的值。c.湍流动能和湍流耗

44、散率名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 25 页，共 39 页 - - - - - - - - - 26 在使用各种k模型对湍流进行计算时，需要给定进口边界上的湍流动能k（Turbulent Kinetic Energy） 和湍流耗散率（Turbulent Dissipation Rate）的估算值。目前没有理论上的精确计算这两个参数的公式，只能通过试验得到。 但不可能对各种各样的流动都去做试验，因此，我们必须借助文献中已有的近似公式来估算（许多商用CF

45、D软件也是这样处理的） 。对于没有已知条件的情况，可根据湍流强度I和特征长度L，由下式粗略估计k和的分布：23()2avgkuI3/23/4kCl式中C为湍流模型中的一个经验常数，其值约等于0.09。在没有直接输入湍流动能k和湍流耗散率的情况下，可以用Intensity and Hydraulic Diameter（湍流强度与水力直径）或Intensity and Length Scale（湍流强度与长度尺度）等办法，利用上述公式确定湍流动能k和湍流耗散率。2压强进口边界条件（ Pressure Inlet）压强进口边界条件用于定义进口流体的压强，可用于不可压缩和可压缩流动。当进口压强已知，而

46、流动速度或流量未知时，可使用压名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 26 页，共 39 页 - - - - - - - - - 27 强进口边界条件。 压强进口边界条件也可用于定义外部或非受限流动的“自由边界”。定义总压与静压首先道总压（p）与静压（0p）的关系如下（根据伯努力积分） ：2012ppv在 Momentum （动量）选项卡内， Reference Frame 为参考值，有绝对值（ Absolute ）与相对于临近区域值（ Relative t

47、o Adjacent Zone）两个选项供选择，一般保持默认的绝对值，Gauge Total Pressure （表总压）文本框中输入总压的值。静压在 FLUENT 中被称为 Supersonic/Initial Gauge Pressure（超音速/ 初始表压）， 如果进口流动是超音速的或者是准备压强进口边界条件进行计算的初始化工作，则必须定义静压。在流场为亚音速时，名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - - - - - - - - - - - - - - - 第 27 页，共 39 页 - - - - - - -

48、 - - 28 FLUENT 将忽略 Supersonic/Initial Gauge Pressure（超音速 /初始表压）的输入数据，而用驻点参数求出静压。（例如，对于空气，当速度为100/m s时，可先求出其动压为22111.225 100612522v，静压等于总压减去动压得到。 ）对于高雷诺数的计算，由于采用压强入口边界条件，湍流强度并不易求出，这里给出推荐范围1%5%（当然也可以采用预估的方法估计一下） 。3压强出口边界条件（ Pressure Outlet）压强出口边界条件在流场出口边界上定义静压，而静压的值仅在流场为亚音速时使用。 如果在出口边界上流场达到超音速，则边界上的压强

49、将从流场内部通过插值得到。其他流场变量均从流场内部通过插值获得。在压强出口边界上还需要定义“回流（backflow ） ”条件。回流条件是在压强出口边界上出现回流使用的边界条件，这样计算将更容易收敛。 FLUENT 在压强出口边界条件上可以使用径向平衡条件，同时可以给定预期的流量。有关回流的湍流参数的设置如下图所示，需要说明的是这些湍流参数只有当回流发生时才被使用。 但是即使在计算结果中没有回流出现，也应该将出口条件用真实流场的值设定，这样可以在计算过程中出现回流时加速收敛。名师归纳总结 精品学习资料 - - - - - - - - - - - - - - -精心整理归纳 精选学习资料 - -

50、 - - - - - - - - - - - - - 第 28 页，共 39 页 - - - - - - - - - 29 4出流边界条件（ Outflow ）出流边界条件用于模拟在求解前流速和压力未知的出口边界。在该边界上，用户不需要定义任何内容（除非模拟辐射传热、粒子的离散相及多口出流）。该边界条件适用于出口处的流动是完全发展的情况。 所谓完全发展，意味着出流面上的流动情况由区域内部外推得到，且对上游流动没有影响。 出游边界条件不能用于可压流动，也不能与压力进口边界一起使用 （压力进口边界条件可与压力出口边界条件一起使用） 。出口边界条件的设置比较简单，只需给定所指定的出流边界上流体的流出