FLUENT知识点解析良心出品必属精品.docx

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

2、缩渐扩管的 无粘与可压缩流动模拟）。c.对于某些高导热系数比或高宽纵比的网格， 使用单精度求解器 可能会遇到收敛性不佳或精确度不足不足的问题， 此时，使用双精度 求解器可能会有所帮助。1李鹏飞，徐敏义，王飞飞.精通CFD工程仿真与案例实战：FLUENT GAMBIT ICEM CFD TecplotM. 北京, 人民邮电出版社，2011:114-116FLUENT Launcher口 ！回 4W2.网格光顺化用光滑和交换的方式改善网格：通过Mesh下的Smooth/Swap来实 现，可用来提高网格质量，一般用于三角形或四边形网格，不过质量 提高的效果一般般，影响较小，网格质量的提高主要还是在网

3、格生成 软件里面实现，所以这里不再用光滑和交换的方式改善网格，其原理 可参考FLUEN性攻略（已下载）。3.Pressure-based 与 Density-based求解器设置如图。下面说一说 Pressure-based和Density-based 的区别：Pressure-Based Solver 是Fluent的优势，它是基于压力法的求 解器，使用的是压力修正算法，求解的控制方程是标量形式的，擅长 求解不可压缩流动，对于可压流动也可以求解； Fluent 6.3 以前的版本求解器，只有 Segregated Solver 和Coupled Solver ,其实也 是Pressure-B

4、ased Solver 的两种处理方法；Density-Based Solver 是Fluent 6.3 新发展出来的，它是基于 密度法的求解器，求解的控制方程是矢量形式的，主要离散格式有 Roe, AUSM +该方法的初衷是让Fluent具有比较好的求解可压缩流 动能力，但目前格式没有添加任何限制器，因此还不太完善；它只有 Coupled的算法；对于低速问题，他们是使用 Preconditioning 方法 来处理，使之也能够计算低速问题。Density-Based Solver 下肯定 是没有SIMPLEC PISO这些选项的，因为这些都是压力修正算法，不 会在这种类型的求解器中出现的；一

5、般还是使用 Pressure-BasedSolver解决问题。基于压力的求解器适用于求解不可压缩和中等程度的可压缩流体 的流动问题。而基于密度的求解器最初用于高速可压缩流动问题的求 解。虽然目前两种求解器都适用于各类流动问题的求解 （从不可压缩流动到高度可压缩流动），但对于高速可压缩流动而言，使用基于密 度的求解器通常能获得比基于压力的求解器更为精确的结果。GeneralMeshSc：dbe“， Check Report QuiitySolverType Pressure-eased Density iBasedVeJodty Formubtion q AbsoluteRelativeTime

6、2D Spacea Steady9 Planar, TransientAxisymmetricAxisYTTimetric Swirl口 Svity | units HeJp4.axisymmetric 和 axisymmetric swirl从字面的意思很好理解axisymmetric和axisymmetric swirl 的 差别：axisymmetric :是轴对称的意思，也就是关于一个坐标轴对称，2D的 axisymmetric 问题仍为 2D 问题。而 axisymmetric swirl :是 轴对称旋转的意思，就是一个区域关于一条坐标轴回转所产生的区域，这产生的将是一个回转体，是

7、3D的问题。在Fluent中使用这个, 是将一个3D的问题简化为2D问题，减少计算量，需要注意的是，在 Fluent中，回转轴必须是x轴。5.操作工况参数(Operating Conditions )操作压力的介绍关于参考压力的设定，首先需了解有关压力的一些定义。 ANSYSFLUEN并1有以下几个压力，即Static Pressure （静压）、DynamicPressure （动压）与 Total Pressure （总压）；Absolute Pressure（绝对压力）、Relative Pressure （参考压力）与 Operating Pressure （操作压力）。这些压力间的关

8、系为，Total Pressure （总压）=Static Pressure （静压）+Dynamic Pressure （动压）；Absolute Pressure （绝对压力）=Operating Pressure （操作压力）+Gauge Pressure （表压）。其中，静压、动压和总压是流体力学中关于压力的概念。静压是测量到的压力，动压是有关速度动能的压力，是流动速度能量的体现。而绝对压力、操作压力和表压是 FLUENTH入的压力参考量，在 ANSYS FLUENT,所有设定的压力都默认为表压。这是考虑到计算 精度的问题。操作压力的设定设定操作压力时需要注意的事项如下：对于不可压缩理

9、想气体的流动，操作压力的设定直接影响流体 密度的计算，因为对于理想气体而言，流动的密度由理想气体 方程获得，理想气体方程中的压力为操作压力。对于低马赫数的可压缩流动而言， 相比绝对静压，总压降是很 小的，因此其计算精度很容易受到数值截断误差的影响。 需要 采取措施来避免此误差的形成，ANSYS FLUENT过采用表压（由绝对压力减去操作压力）的形式来避免截断误差的形成，操作压力一般等于流场中的平均总压。对于高马赫数可压缩流动的求解而言，因为此时的压力比低马 赫可压缩流动的大得多，所以求解过程中的截断误差的影响不 大，可以不设定表压。由于 ANSYS FLUENT所有需输入的压 力都为表压，因此

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

11、 ANSYS FLUENT每次迭代后要调整表压值。这个过程通过使用参考压力位置处（或该 位置附近）节点的压力完成。因此，参考压力位置处的表压应一直为 0。如果使用了压力边界条件，则不会使用到上述关系，因此参考压 力位置不被使用。参考压力位置默认为等于或接近（0, 0, 0）的节点中心位置。实 际计算中可能需要设置参考压力位置到绝对静压已知的位置处。在 Operating Conditions 对话框中的 Reference Pressure Location 选项组中设置新的参考压力位置的 x, y, z的坐标即可。如果要考虑某一方向的加速度，如重力，可以勾选 Gravity复选 框。对于VO

12、F计算，应当选择 Specified Operating Density ,并且 在Operating Density 下为最轻相设置密度。这样做排除了水力静 压的积累，提高了 round-off精度为动量平衡。同样需要打开 Implicit Body Force ,部分平衡压力梯度和动量方程中体积力，提 高解的收敛性。Reference Pressure Location （参考压强位置）应是位于流体永 远是100%勺某一相（空气）的区域，光滑和快速收敛是其基本条件。二、求解模型的设定1.流动模型的设置无粘模型理想流体是一种设想的没有粘性的流体，在流动时各层之间没有 相互作用的切应力，即没有

13、内摩擦力。十分明显，理想流体对于切向 变形没有任何抗拒能力。应该强调指出，真正的理想流体在客观实际 中是不存在的，它只是实际流体在某些条件下的一种近似模型。在Inviscid 流动模型应用方面，无粘流动忽略了粘性对流动的影 响，这对高雷诺数的流动是合适的，因为高雷诺数流动惯性力的作用 远大于粘性力的作用，粘性力可以忽略，所以可以将其考虑成无粘流 动。无粘流动的求解更快，其激波在某些值上预测的偏高。无粘流动 能对流动状态和激波位置进行快速预测。马赫数与激波马赫数的定义是v M = a它表示流体的流动速度与当地声速之比，是一个无量纲的参量。对应于M 1这三种情况的流动分别称为亚声速流、声 速流和超

14、声速流。当马赫数很小时，速度的相对变化只能引起很小的 密度相对变化，但当马赫数很大时，则将引起较大的密度相对变化， 这也说明了马赫数是流体压缩性的一个表征。当飞机、炮弹和火箭以超音速飞行时，或者发生强爆炸、强爆震 时，气流受到急剧的压缩，压强和密度突然显著增加，这时所产生的 压强扰度将比声速大得多的速度传播，波阵面所到之处气流的各种参 数都将发生显著变化，参数突跃。这样一个强间断面叫做激波阵面。渐缩渐扩管的流动是计算流体力学模拟的经典问题之一。在这类流动中，激波的出现是流动中可压缩效应的体现。 精确的激波模拟是 CFDW究的热点之一。为了更好捕捉压力梯度，需要采用较细的网格 并结合合适的数值模

15、拟和格式。很多实际模拟中，局部网格的自适应 会很有帮助。层流模型流动有层流和湍流之分，判断湍流的标准可以参考2,这里写出内流的判断标准：UD P Re- 2300对于内流而言，一般大多数流动都是湍流，一般不使用湍流模型。 而对一些外流而言（如外掠平板或是外掠障碍物），则很有可能是层 流运动。湍流模型的评价与选择a. k-8湍流模型这里我们使用的湍流模型是 Standard k-名模型，这种模型应用 较多，计算量适中，有较多数据积累和比较高的精度，对于曲率较大 和压力梯度较强等复杂流动模拟效果欠佳。一般工程计算都使用该模 型，其收敛性和计算精度能满足一般的工程计算要求， 但模拟旋流和绕流时有缺陷

16、。壁面函数的选择对于有壁面的流动，当主流为充分发展湍流时，根据离壁面法线 距离不同，可将流动划分为壁面区（或称内区、近壁区）和核心区（或 称外区）。核心区是完全湍流区，为充分发展的湍流。在壁面区，由于有壁面的影响，流动与核心区不同。壁面区可分 为3个子层：粘性底层、过渡层和对数率层。2李鹏飞，徐敏义，王飞飞.精通CFD工程仿真与案例实战：FLUENT GAMBIT ICEM CFD TecplotM. 北京, 人民邮电出版社，2011:122粘性底层是一个紧贴壁面的极薄层，在动量、热量和质量的交换 过程中粘性力起主要作用，而湍流切应力可以忽略，因此流动几乎可 以看成层流流动，且在平行于壁面方向

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

18、界模型和两层边界模型结合起来，对 低雷诺数流动或者复杂近壁面现象很适合， 湍流模型在内层上得到了 修正。表格2 几种壁面处理方法的比较优点缺点标准壁面函数法应用较多，计算量 小，后较局的精度适合高雷诺数流动，对低雷诺数流 动问题，有压力梯度、高度蒸腾和 大的体积力、低雷诺数和高速三维 流动问题不适合非平衡壁面函数法考虑了压力梯度，可 以计算分离，在附着 以及撞击问题对低雷诺数流动问题，有较强压力梯度、强体积力及强三维性问题不适合增强壁面处理不依赖壁面法则，对 于复杂流动，特别是要求网格密，因而要求计算机处理 时间长，内存大低雷诺数流动很适12! VJscous ModelNearWall Tr

19、eatmentq Standard Wail Functions.Scalable Wall FunctionsNon-Equilibrium Wall Functions Enhanced Wall TreatmentUser-De fined 7all FunctionsOptionsViscous Heatng Curvature CorrectionTurbUent Vhcosity none kandti NumbersTO Prandti Number none TDR Prandtl Number none -Energy P日 ndti Number none OK I Can

20、cel | Help2.多相流模型1VO葭型该模型通过求解单独的动量方程和处理穿过区域的每一流体的容 积比来模拟两种或三种不能混合的流体。典型的应用包括流体喷射、 流体中大泡运动、流体在大坝坝口的流动、气液界面的稳态和瞬态处理等。一般而言VOF主要适用于非稳态的多相流模型， 仅对某些特定问题的多相流模型的稳态问题能够适用。VO昉法适用于计算空气和水这样不能互相掺混的流体流动，对 于分层流和活塞流，最方便的就是选择 VO破型。需要注意的是，对 于湍流模型的设置，丫0环能用于无粘流，也不能用大涡模拟3。Geo-Reconstruct 格式Geo-Reconstruct 格式（在 Solution

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

22、Qoff Volume of FluidC Mixture Eulerian Wet Steamurrber of Euierian PhasesCoupled level Set + V口F I I Level SetVolume Fraction ParametersSdierne ExplicitO IrnplidtVolume Fraction Cutoff le-06Courant Number产 |DefaultBody Force Formulaton implicit Body ForceOptions匚 Open Channel Flow匚 Open Channel Wave

23、 BCI Zonal DiscretizatjoriOK Cancel Help2Mixture模型这是一种简化的多相流模型，用于模拟各种有不同速度的多相流， 但是假定了在短空间尺度上局部的平衡。相之间的耦合应当是很强 的。它也用于模拟有强烈耦合的各向同性多相流和各向以相同速度运 动的多相流。典型的应用包括沉降(sedimentation )、气旋分离器、 低载荷作业下的多粒子流动、气相容积率很低的泡状流。Mixture Parameters一般需要勾选 Mixture Parameters 中的 Slip Velocity 复选框, 以此来求解滑移速度模型，因为在多相流中各种组分的速度有很大

24、不同。对于求解一个均匀的多相流问题可以选择不做滑移速度的计算,可以在 mixture parameters 选项下将 slip velocity 关掉。Eulerian 模型该模型可以模拟多相分离流及相互作用的相，相可以是液体、气体、固体。与在离散相模型中Eulerian-Lagrangian 方案只用于离散 相不同，在多相流模型中Eulerian方案用于模型中的每一项。3.固化与熔化模型FLUEN邛用“始多孔度(enthalpy-porosity )”技术模拟流体 的固化和熔化(Solidification/Melting )过程。在流体的固化和熔 化问题中，流场可以分成流体区域、固体区域和

25、两者之间的糊状区域。 “始多孔度”技术采用的计算策略是将流体在网格单元内占有的体积百分比定义为多孔度（porosity ）,并将流体和固体并存的糊状区 域看作多孔介质区进行处理。在流体的固化过程中，多孔度从 1降低到0;反之，在熔化过程中，多孔度则从 0升至1。“始多孔度” 技术通过在动量方程中添加汇项（即负的源项）模拟因固体材料存在 而出现的压强降。“始-多孔度”技术可以模拟的问题包括纯金属或二元合金中的 固化、熔化问题、连续铸造加工过程等。计算中可以计算固体材料与 壁面之间因空气的存在而产生的热阻， 固化、熔化过程中组元的输运等等。需要注意的是，在求解固化、熔化问题的过程中，只能采用分 离

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

27、ture （固相点温度） 和 Liquidus Temperature （液相点温度）。 如果计算中涉及组元输运过程，则必须同时定义溶剂的融解温度（Melting Temperature）,同时需要定义熔化物的液相线相对于浓度的斜率(Slope of Liquidus Line )、分配系数(PartitionCoefficient )和固体中的扩散速率(Diffusion in Solid )等参数。设置边界条件除了常规的边界条件设置，对于固化和熔化问题还有一些特殊设置，其中包括：在计算壁面接触热阻时设置接触热阻( ContactResistance )。这个参数在 Wall (壁面)面板中

28、的 Thermal Conditions(热力学条件)下给定。ShearCondition （剪切条件）下选择 Marangoni Stress （Marangoni 应力） 选项。如果计算拉出速度，则在边界条件中的速度边界条件将被用于拉出速度的计算。口 Solidification and Melting三、相设置相设置一般用于多相流的设置，对于相设置，这里主要讲一下Interaction 的设置，如图：PhasesPhaseswater - Primary Phaseair Secondary PhaseFit. Interacton,1-HelpInteraction 设置Interac

29、tion 设置用来定义两相的相互作用，具有多个选项卡,如图。Drag选项卡针对每对物相，在下拉菜单中选择阻力函数。其中包括schiller-naumann 模型、morsi-alexander 模型、symmetric （对称） 模型等用于流体与流体之间阻力计算的模型，也包括 wen-yu模型、gidaspow模型、syamlal-obrien 模型等用于液体与固体之间阻力计算的模型，还包括syamlal-obrien-symmetric 模型用于固体与固体之间的阻力计算。除此之外，还可以将阻力函数定义为constant （常 数），或者选择user-defined （用户定义）由用户自己定义

30、阻力函数。如果计算中不需要设定阻力，还可以选择 none （不计阻力）选项。阻力设置的相关原理比较复杂，可参考帮助，一般保持默认的 schiller-naumann 设置不变。Surface Tension 选项卡Surface Tension选项卡用来定义表面张力，如果相包含壁面粘附，可勾选“ Wall Adhesion ”复选框。 四、Cell Zone ConditionFrame Motion 选项对于流体，可以通过FrameMotion选项确定坐标运动方式（如离 心泵内部流体的旋转使用运动参考系模型），如图：Rotation-Axis Origin 指转轴，Rotational Velocity 值旋转速度。五、边界条件设置1.速度进口边界条件(Velocity Inlet )速度进口边界条件用进口