付费下载
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、一、基本设置1. Double Precision 的选择启动设置如图,这里着重说说Double Precision (双精度)复选框,对于大多数情况,单精度求解器已能很好的满足精度要求,且计算量小,这里我们选择单精度。然而对于以下一些特定的问题,使用双精度求解器可能更有利1。a. 几何特征包含某些极端的尺度(如非常长且窄的管道),单精度求解器可能不能足够精确地表达各尺度方向的节点信息。b. 如果几何模型包含多个通过小直径管道相互连接的体,而某一个区域的压力特别大(因为用户只能设定一个总体的参考压力位置),此时,双精度求解器可能更能体现压差带来的流动(如 渐缩渐扩管的无粘与可压缩流动模拟)。c
2、. 对于某些高导热系数比或高宽纵比的网格,使用单精度求解器可能会遇到收敛性不佳或精确度不足不足的问题,此时,使用双精度求解器可能会有所帮助。FLUENT LauncherOptasDouble P redsianPracsssrg OipHiaris Serial -ParallelANSYSDimersioniod2D"r- II-"If I 3DDisplay Options叼 DispIdA MeshAfier Pleading 団 Embed Graphics Windows0 V/orkbench Color S chenrie圈 Show More Opt io
3、 ns0K1DeiaJtCancelf Help 生2. 网格光顺化用光滑和交换的方式改善网格:通过Mesh下的Smooth/Swap来实现,可用来提高网格质量,一般用于三角形或四边形网格,不过质量提高的效果一般般,影响较小,网格质量的 提高主要还是在网格生成软件里面实现,所以这里不再用光滑和交换的方式改善网格,其原理可参考FLUENT全攻略(已下载)。3. Pressure-based 与 Density-based求解器设置如图。下面说一说Pressure-based和Density-based的区别:P ressure-Based Solver是Flue nt的优势,它是基于压力法的求解
4、器,使用的是压力修正 算法,求解的控制方程是标量形式的,擅长求解不可压缩流动,对于可压流动也可以求解;Flue nt 6.3以前的版本求解器,只有Segregated Solver和 Cou pled Solver,其实也是P ressure-Based Solver的两种处理方法;Den sity-Based Solver是Flue nt 6.3新发展出来的,它是基于密度法的求解器,求解的控 制方程是矢量形式的,主要离散格式有Roe, AUSM+,该方法的初衷是让Flue nt具有比较1 李鹏飞,徐敏义,王飞飞精通CFD工程仿真与案例实战:FLUENT GAMBIT ICEM CFD Tec
5、plotM.北京,人民邮电出版社,2011:114-116好的求解可压缩流动能力,但目前格式没有添加任何限制器,因此还不太完善;它只有Coupled的算法;对于低速问题,他们是使用Preconditioning方法来处理,使之也能够计算低速问题。Density-Based Solver下肯定是没有 SIMPLEC,PISO这些选项的,因为这些都是 压力修正算法,不会在这种类型的求解器中出现的;一般还是使用P ressure-Based Solver解决问题。基于压力的求解器适用于求解不可压缩和中等程度的可压缩流体的流动问题。而基于密度的求解器最初用于高速可压缩流动问题的求解。虽然目前两种求解器
6、都适用于各类流动问题的求解(从不可压缩流动到高度可压缩流动),但对于高速可压缩流动而言,使用基于密度的求解器通常能获得比基于压力的求解器更为精确的结果。GeneralMesh Check Report Quality Display" SolverTypeI* Ftessure-eased 广 Density-Based呻odtv Formulation 唾 Absolute广 RelatjveTirne11 Steady 'C Transient四5p电e"4 PlrtarC Axisyrnmetrk匚 Axismmetric 5wirtunits. Gravit
7、y4. axisymmetric 禾R axisymmetric swirl从字面的意思很好理解axisymmetric和axisymmetric swirl的差别:axisymmetric :是轴对称的意思,也就是关于一个坐标轴对称,2D的axisymmetric问题仍为2D问题。而axisymmetric swirl :是轴对称旋转的意思,就是一个区域关于一条坐标轴 回转所产生的区域,这产生的将是一个回转体,是3D的问题。在Flue nt中使用这个,是将一个3D的问题简化为2D问题,减少计算量,需要注意的是,在Flue nt中,回转轴必须是x轴。5.操作工况参数(Op erati ng C
8、on ditio ns ) 操作压力的介绍关于参考压力的设定,首先需了解有关压力的一些定义。ANSYS FLUENT中有以下几个压力,即 Static Pressure(静压)、Dynamic Pressure (动压)与 Total Pressure(总压);Absolute P ressure (绝对压力)、Relative P ressure (参考压力)与 Op erat ing P ressure (操作压力)。这些压力间的关系为,Total Pressure (总压)=Static Pressure (静压)+Dynamic Pressure(动压);Absolute Pressu
9、re (绝对压力)=Operating Pressure (操作压力)+Gauge Pressure (表 压)。静压是测量到的压力,动压是其中,静压、动压和总压是流体力学中关于压力的概念。 有关速度动能的压力,是流动速度能量的体现。而绝对压力、操作压力和表压是FLUENT引入的压力参考量,在 ANSYS FLUENT 中,所有设定的压力都默认为表压。这是考虑到计算精度的问题。 操作压力的设定设定操作压力时需要注意的事项如下:对于不可压缩理想气体的流动,操作压力的设定直接影响流体密度的计算,因为对于理想气体而言,流动的密度由理想气体方程获得,理想气体方程中的压力为操作压力。对于低马赫数的可压缩
10、流动而言,相比绝对静压,总压降是很小的,因此其计算精度很容易受到数值截断误差的影响。需要采取措施来避免此误差的形成,ANSYSFLUENT通过采用表压(由绝对压力减去操作压力)的形式来避免截断误差的形 成,操作压力一般等于流场中的平均总压。对于高马赫数可压缩流动的求解而言,因为此时的压力比低马赫可压缩流动的大得多,所以求解过程中的截断误差的影响不大,可以不设定表压。由于ANSYSFLUENT中所有需输入的压力都为表压,因此此时可以将操作压力设定为0(这样可以最小化由于压力脉动而引起的误差),使表压与绝对压力相等。如果密度设定为常数或者其值由通过温度变化的函数获得,操作压力并没有在计算密度的过程
11、中被使用。默认的操作压力为101325Pa。操作压力的设定主要基于两点考虑,一是流动马赫数的大小,二是密度计算方法。表格1操作压力的推荐设置密度关系式马赫数操作压力理想气体定律大于0.10或约等于流场的平均压力理想气体定律小于0.1约等于流场的平均压力关于温度的函数不可压缩不使用常数不可压缩不使用不可压缩的理想气体不可压缩约等于流场的平均压力关于参考压力位置的设定对于不涉及任何压力边界条件的不可压缩流动,ANSYS FLUENT在每次迭代后要调整表压值。这个过程通过使用参考压力位置处(或该位置附近)节点的压力完成。因此,参考 压力位置处的表压应一直为0。如果使用了压力边界条件,则不会使用到上述
12、关系,因此参考压力位置不被使用。参考压力位置默认为等于或接近(0,0,0)的节点中心位置。实际计算中可能需要设置参考压力位置到绝对静压已知的位置处。在Operating Conditions对话框中的 RefereneePressure Location选项组中设置新的参考压力位置的x, y, z的坐标即可。如果要考虑某一方向的加速度,如重力,可以勾选Gravity复选框。对于 VOF 计算,应当选择 Sp eeified Op erat ing Den sity,并且在 Op erat ing Den sity 下为 最轻相设置密度。这样做排除了水力静压的积累,提高了round-off精度为
13、动量平衡。同样需要打开Imp licit Body Force,部分平衡压力梯度和动量方程中体积力,提高解的收敛性。(空Referenee Pressure Location (参考压强位置)应是位于流体永远是100%的某一相气)的区域,光滑和快速收敛是其基本条件。Op erating ConditionsGravityOperdtng Pressure (jDascal)1013250 GravityGravitational Accel era ton Crr/s2)-9,S|P|Reference Pressure Locatjori回PIV Cm /s2)BoLissiiesq Par
14、ametersOperdtrig Temperature Ic)Pressure2B3J6Variable-Dena ty Par ametersI Specified Operating DenatyOKCancel Help二、求解模型的设定1. 流动模型的设置无粘模型即真正所以理想流体是一种设想的没有粘性的流体, 在流动时各层之间没有相互作用的切应力, 没有内摩擦力。十分明显,理想流体对于切向变形没有任何抗拒能力。应该强调指出,的理想流体在客观实际中是不存在的,它只是实际流体在某些条件下的一种近似模型。在In viscid流动模型应用方面,无粘流动忽略了粘性对流动的影响,这对高雷诺数的流
15、 动是合适的,因为高雷诺数流动惯性力的作用远大于粘性力的作用,粘性力可以忽略,可以将其考虑成无粘流动。无粘流动的求解更快,其激波在某些值上预测的偏高。无粘流动能对流动状态和激波位置进行快速预测。马赫数与激波 马赫数的定义是它表示流体的流动速度与当地声速之比, 和M 1这三种情况的流动分别称为亚声速流、 的相对变化只能引起很小的密度相对变化, 化,这也说明了马赫数是流体压缩性的一个表征。a对应于M是一个无量纲的参量。对应于M 1,M 1 声速流和超声速流。 当马赫数很小时,速度 但当马赫数很大时,则将引起较大的密度相对变当飞机、炮弹和火箭以超音速飞行时,或者发生强爆炸、强爆震时,气流受到急剧的压
16、缩,压强和密度突然显著增加,这时所产生的压强扰度将比声速大得多的速度传播,波阵面所到之处气流的各种参数都将发生显著变化,参数突跃。这样一个强间断面叫做激波阵面。渐缩渐扩管的流动是计算流体力学模拟的经典问题之一。在这类流动中,激波的出现是流动中可压缩效应的体现。精确的激波模拟是CFD研究的热点之一。为了更好捕捉压力梯度,需要采用较细的网格并结合合适的数值模拟和格式。很多实际模拟中,局部网格的自适应会很有帮助。层流模型流动有层流和湍流之分,判断湍流的标准可以参考2,这里写出内流的判断标准:Re2300对于内流而言,一般大多数流动都是湍流, 一般不使用湍流模型。而对一些外流而言(如 外掠平板或是外掠
17、障碍物),则很有可能是层流运动。 湍流模型的评价与选择a. k 湍流模型这里我们使用的湍流模型是 Standard k模型,这种模型应用较多,计算量适中,有较多数据积累和比较高的精度,对于曲率较大和压力梯度较强等复杂流动模拟效果欠佳。一般工程计算都使用该模型,其收敛性和计算精度能满足一般的工程计算要求,但模拟旋流和绕流时有缺陷。壁面函数的选择对于有壁面的流动,当主流为充分发展湍流时,根据离壁面法线距离不同,可将流动划 分为壁面区(或称内区、近壁区)和核心区(或称外区)核心区是完全湍流区,为充分发展的湍流。在动量、热量和质量的交换过程中粘性力起主要作 因此流动几乎可以看成层流流动,且在平行于壁面
18、方向上的速在壁面区,由于有壁面的影响, 流动与核心区不同。 壁面区可分为3个子层:粘性底层、 过渡层和对数率层。粘性底层是一个紧贴壁面的极薄层, 用,而湍流切应力可以忽略, 度呈线性分布。在此层中,粘性力和湍流切应力的作用相当, 流动状况较为 实际工程计算中由于过渡层厚度极小, 可不考虑此层,直接过渡层处于粘性底层之外, 复杂,很难用公式或定律表述。 以对数率层的方法处理。对数率层处于近壁区的最外层,粘性力的影响不明显, 湍流切应力占主要地位, 流动处于充分发展的湍流状态,流速分布接近对数律。壁面区内不同子层的高度和速度可以沿壁面法向的无量纲高度和无量纲速度表达。2 李鹏飞,徐敏义,王飞飞.精
19、通CFD工程仿真与案例实战:FLUENT GAMBIT ICEM CFD TecplotM. 北京, 人民邮电出版社,2011:122其中,U是流体的时均速度,U是壁面摩擦速度,UW是壁面切应力,y是壁面的垂直距离。在605时,区域为粘性底层,此时速度沿壁面法线方向呈线性分布,即y 300时,流动处于对数率层,此时速度沿壁面法线方向呈对数率分布,即u 2.5lny 5.5。使用标准的或者非平衡的壁面函数(Re 106)。O标准壁面函对于大多数高雷诺数情况壁面函数法的本质是,对于湍流核心区的流动使用k模型求解,而在壁面区并不进行求解,直接使用半经验公式得出该区域的速度等物理量。FLUENT提供了
20、多种壁面函数处理方式,如标准壁面函数法、非平衡壁面函数法和增 强壁面处理。标准壁面函数法利用对数校正法提供了必需的壁面边界条件(对于平衡湍流边界层) 而非平衡壁面函数法用来改善高压力梯度、分离、再附和滞止等情况下的结果。数法和非平衡壁面函数法都允许在近避免区域上使用较粗的网格。优点缺点标准壁面函数 法应用较多,计算量小,有较 高的精度适合高雷诺数流动, 对低雷诺数流动问题, 有 压力梯度、高度蒸腾和大的体积力、低雷诺数 和高速三维流动问题不适合非平衡壁面函数法考虑了压力梯度,可以计算 分离,在附着以及撞击问题对低雷诺数流动问题,有较强压力梯度、强体 积力及强三维性问题不适合增强壁面处理不依赖壁
21、面法则,对于复杂 流动,特别是低雷诺数流动 很适合要求网格密,因而要求计算机处理时间长,内存大增强壁面处理选项把混合边界模型和两层边界模型结合起来, 近壁面现象很适合,湍流模型在内层上得到了修正。表格2几种壁面处理方法的比较对低雷诺数流动或者复杂Q VbcoiJ5 ModdMmJsIModel C«nstaritsI _ fl riVKod(j LaminarII Sp alart Allmar as. (1 eqr 'O«k*ep£i1on (2 eqn)匚«t-o*necia (3eqn)I _ * TransfLon k-kJ -omega
22、 (3 eqnj»_ * TransfLori SST 斗 eqn'JI Revriolds Stress (5 ew)S仙也pB低 Simd目ton (SAS)k<epsilon Model 佔"standard (_-I FtNQI * RealizableNeerWall Tretmert崎n Standard Wall FiincliOriS(_ “ ScaleWe Wall FunctionsIJ Non-Equilibrium Wai Fijnctioris 厂 * Enhanced Wall Treatjrwnt 厂 N User-Defined
23、 Wall FupictonsOptionsZU WiscoLB HeatingZlCurvature Correctionc0.09Ci-£pdon1.44CZ-EpsibnTKE FraudtJ NumberInoneTPtardtJ NumbersTXE Prandtl Nuntjernone>TDR. Prand-H Numbernor>e<&ierqy Prandtl NumbernoneUser 'Pefined FunctionsTurbulent ViscoartyHfclp 2. 多相流模型 VOF模型该模型通过求解单独的动量方程和
24、处理穿过区域的每一流体的容积比来模拟两种或三 种不能混合的流体。典型的应用包括流体喷射、流体中大泡运动、流体在大坝坝口的流动、 气液界面的稳态和瞬态处理等。一般而言VOF主要适用于非稳态的多相流模型,仅对某些特定问题的多相流模型的稳态问题能够适用。VOF方法适用于计算空气和水这样不能互相掺混的流体流动,对于分层流和活塞流,最方便的就是选择 VOF模型。需要注意的是,对于湍流模型的设置,VOF不能用于无粘流, 也不能用大涡模拟3。Geo-Rec on struct 格式Geo-Reconstruct格式(在Solution Methods中设置)是一种较为精确的追踪自由表面的 计算格式,广泛地应
25、用于瞬变流的VOF问题中,但必须注意的要使用该格式VOF模型必须使用显示离散格式(在 VOF模型设置选项设置)。Body Force Formulati on为提高解的收敛性,对于涉及到表面张力的计算,建议在Body Force Formulation中勾选Imp licit Body Force。这样做由于压力梯度和动量方程中表面张力的部分平衡,从而提高 解的收敛性。3李进良,李承曦,胡仁喜.精通FLUENT63流场分析M.北京,化学工业出版社,2009:231-236ModelIVIukiphagehkmbef of Etilefian PhasesO Off tTJ Volurrw of
26、 HurdO MixtureO Ejlerian'.VetStea;v Level SetSctierre.fi'iiEiphdtJlmpbot open diaonel Flow ll Open Channel Wave BC n Zonal DiscretijstQriVftlunw Fration Paramt也肱OptiopieVolume Fracton Cutoff |le-06匚ffljrant hhjmbefBody Force Formulatiori onpiicit Bod* ForceHelp Mixture模型这是一种简化的多相流模型,用于模拟各种有不
27、同速度的多相流,但是假定了在短空间尺度上局部的平衡。相之间的耦合应当是很强的。它也用于模拟有强烈耦合的各向同性多相 流和各向以相同速度运动的多相流。典型的应用包括沉降(sedimentation)、气旋分离器、低载荷作业下的多粒子流动、气相容积率很低的泡状流。Mixture P arameters一般需要勾选 Mixture Parameters中的Slip Velocity复选框,以此来求解滑移速度模型, 因为在多相流中各种组分的速度有很大不同。对于求解一个均匀的多相流问题可以选择不做滑移速度的计算,可以在 mixture parameters选项下将slip velocity关掉。Mfid
28、eJMulti ph 曰旨 e ModelNumber of Eulerian Phasesg OffQ Votume of Flud (i/MiitueCl Eulerian內讨4妒刁mMixture Psrarreters0 Slifj VeltHzitvBody Faroe Porrmjiation0 Impliot Body ForceOKCarioel Help Eulerian模型该模型可以模拟多相分离流及相互作用的相,相可以是液体、 气体、固体。与在离散相模型中Eulerian-Lagrangian方案只用于离散相不同,在多相流模型中Eulerian方案用于模型中的每一项。3.
29、固化与熔化模型FLUENT 采用“焓多孔度(enthalpy-porosity )”技术模拟流体的固化和熔化(Solidificatio n/Melt ing )过程。在流体的固化和熔化问题中,流场可以分成流体区域、固体 区域和两者之间的糊状区域。“焓一多孔度”技术采用的计算策略是将流体在网格单元内占有的体积百分比定义为多孔度(P orosity),并将流体和固体并存的糊状区域看作多孔介质区 进行处理。在流体的固化过程中,多孔度从1降低到0;反之,在熔化过程中,多孔度则从0升至1。“焓一多孔度”技术通过在动量方程中添加汇项(即负的源项)模拟因固体材料 存在而出现的压强降。“焓多孔度”技术可以模
30、拟的问题包括纯金属或二元合金中的固化、熔化问题、连续铸造加工过程等。计算中可以计算固体材料与壁面之间因空气的存在而产生的热阻,固化、 熔化过程中组元的输运等等。需要注意的是,在求解固化、熔化问题的过程中,只能采用分离算法,只能与 VOF模型配合使用,不能计算可压缩流,不能单独设定固体材料和流体材 料的性质,同时在模拟带反应的组元输运过程时,无法将反应区限制在流体区域,而是在全流场进行反应计算。 Parameters 定义在Parameters下面定义 Mushy Zone Constant (糊状区域常数)。这个常数的取值范围一般在104到107之间,取值越大沉降曲线就越陡峭,固化过程的计算速
31、度就越快,但是取 值过大容易引起计算振荡,因此需要在计算中通过试算获得最佳数值。 Materials设置在 Materials (材料)面板上,定义 Melting Heat (熔化热)、Solidus Temperature (固相 点温度)和Liquidus Temperature (液相点温度)。如果计算中涉及组元输运过程,则必须同 时定义溶剂的融解温度(Melting Temperature ),同时需要定义熔化物的液相线相对于浓度的斜率(Slope of Liquidus Line )、分配系数(Partition Coefficient )和固体中的扩散速率 (Diffusion
32、in Solid、等参数。设置边界条件除了常规的边界条件设置,对于固化和熔化问题还有一些特殊设置,其中包括:在计算壁面接触热阻时设置接触热阻(Con tact Resista nee)。这个参数在Wall (壁面、面板中的Thermal Conditions (热力学条件、下给定。Shear Condition (剪切条件)下选择如果需要定义壁面上表面张力对温度的梯度,则在 Marangoni Stress ( Marangoni 应力)选项。如果计算拉出速度,则在边界条件中的速度边界条件将被用于拉出速度的计算。Iidiflegation and MeltingMcxlel0 Eclid if
33、i cstidTeb ngP 申宙 mehes陽加心Sr”以r咖讪100000Bd(± DiffusionOlndtide Pull VeJodtiesOKCancel Help三、相设置相设置一般用于多相流的设置,对于相设置,这里主要讲一下In teraction的设置,如图:PbasesPhaseswater - P nrnay P hase air - Seccmdary Phasecdif-.rritftractan.14Interaction 设置Interaction设置用来定义两相的相互作用,其有多个选项卡,如图。Pfidie Interaction"bM I
34、 MML I吓h Sj岳vh那 I阀 Gjg 曲*茸 II jH I 冷出站也口v| Ti/tkETlD町t4r 1*jbU«rlrL«rKtDr.| C口附3让 2dJay"论LHso*Cnori Hff-pDrag选项卡针对每对物相,在下拉菜单中选择阻力函数。其中包括schiller-naumann模型、morsi-alexander模型、symmetric (对称)模型等用于流体与流体之间阻力计算的模型,也 包括wen-yu模型、gidaspow模型、syamlal-obrien模型等用于液体与固体之间阻力计算的 模型,还包括 syamlal-obrien-
35、symmetric模型用于固体与固体之间的阻力计算。除此之外, 还可以将阻力函数定义为constant (常数),或者选择user-defined (用户定义)由用户自己定义阻力函数。如果计算中不需要设定阻力,还可以选择none (不计阻力)选项。阻力设置的相关原理比较复杂,可参考帮助,一般保持默认的schiller-naumann设置不变。Surface Tension 选项卡Surface Tension选项卡用来定义表面张力,如果相包含壁面粘附, 可勾选“ Wall Adhesion ”复选框。四、Cell Zone ConditionFrame Motion 选项对于流体,可以通过 F
36、rame Motion选项确定坐标运动方式(如离心泵内部流体的旋转 使用运动参考系模型),如图:Flluiid7one NameHjid. 7Mstsrial NameairEdit 巨 rraiTkC Mo bon Source T cftra匚 Mesh Motion Fixed VaJurtI IPorajs ZoneLDFReference Fra |iMesh Motior | Paroue Zone | Embedded LES | Rtarton | SoifCft Terms | Fixed '/alues | Mulhiihast |Relative Speofiot
37、iorRelflthe To Cd Tore absokjtsZone Motior FunctionnoneRolatiorAxE 0广pinX fni任EonctintoonstantRolatiqral Vcloci:¥Sfcd (ra/s)orstanlCopy 10 Mesh Motwnirransktionfll VeodtvconstantconstantOK I I Lanoel HelpRotation-Axis Origin 指转轴,Rotational Velocity 值旋转速度。五、边界条件设置1.速度进口边界条件(Velocity Inlet )速度进口边
38、界条件用进口处流场速度及相关流动变量作为边界条件,在速度进口边界条件中,流场进口边界的驻点参数是不固定的。 为了满足进口处的速度条件, 驻点参数将在一 定范围内波动。需要注意的是,因为这种条件中允许驻点参数浮动,所以速度进口边界条件仅适用于不不要让速度进口可压流,如果用于可压流,则可能导致出现非物理解。同时还要注意的是, 条件过于靠近进口内侧的固体障碍物,这样会使驻点参数的不均匀程度大大增加。湍流参数的设置对于一般的流动边界条件,均涉及到湍流参数的定义,在Turbule nee Sp ecificatio n Method(湍流定义方法)下拉列表中,可以简单地用一个常数来定义湍流参数,即通过给
39、定湍流强度、湍流粘度比、水力直径或湍流特征长在边界上的值来定义流场边界上的湍流。a. 湍流强度(Turbulenee Intensity )湍流强度定义如下:|Ju2 V2 w2上式中U , v,w是速度脉动量,u是平均速度。比较常用的是Intensity and Hydraulic Diameter,湍流强度与水力直径的确定有相应的计 算方法,这里只是采用估算来加以确定。内流问题进口处的湍流强度取决于上游流动状态。动,则进口处可以使用低湍流强度。如果上游是充分发展的湍流,到几个百分点。如果管道中的流动是充分发展的湍流,则湍流强度 可以用下面公式计算得到,这个公式是从管流经验公式得到的:如果上
40、游是没有充分发展的未受扰流 则进口处湍流强度可以达I ( turbulenee intensity)I ± 0.16(ReDH) 1/8 u其中:ReH为按等效水力直径Dh计算得到的雷诺数。湍流强度小于1%时,可以认为湍流强度是比较低的,而在湍流强度大于10%时,则可以认为湍流强度是比较高的。比如,当雷诺数为50000时,代入上述公式可得到湍流强度为4%,默认的湍流强度设置为 5% (中等强度)。勺In let7ore Marnevelocity-inlet-5Velocity Specifiedtion MetfiodMagntude, Moriral to BoundaryRef
41、ierence FrarneVelodtv Magnitude On/s)0.2constantni tial Gauge Pressure pascal)lOconstantSxafiaton MethodSupersonic/lMoirertuni | iherTal Radiation | Speoes | DPM | P1ultiDhase| UDSTurbulenceIntiensitv jnd Hydraulic DianneberTurfciJent Ihhensity (鳴)Hydraulic Diameter &n)32OK Cmcgl Help Ib.湍流的长度尺度
42、与水力直径湍流的长度尺度(Turbulenee Length Scale)与水力直径(Hydraulic Diameter )是设置湍 流的重要参数。l则是与大涡结构相关的物理量。在 所以l是受到管道尺寸制约的几何。湍流能量主要集中在大涡结构中,而湍流长度尺度 充分发展的管流中,因为漩涡尺度不可能大于管道直径, 湍流长度尺度l与管道物理尺寸L关系可以表示为:l 0.07L式中的比例因子0.07是充分发展管流中混合长的最大值,而L则是管道直径。在管道截面不是圆形时,L可以取为管道的水力直径。水力直径是在管内流动(internal pipe flow )中引入的,其目的是为了给非圆管流动取一个合适
43、的特征长度来计算其雷诺数。四倍的湿横截面面积与湿圆周长度之商。管道直比在这湍流的特征长取决于对湍流发展具有决定性影响的几何尺度。在上面的讨论中,径是决定湍流发展过程的唯一长度量。如果在流动中还存在其他对流动影响更大的物体,如在管道中存在一个障碍物,而障碍物对湍流的发生和发展过程起着重要的干扰作用。种情况下,湍流特征长就应该取为障碍物的特征长度。从上面的分析可知,虽然上式对于大多数管道流动是适用的,但并不是普遍适用的,某些情况下可以进行调整。在FLUENT中选择特征长L或湍流长度尺度l的方法如下:1)对于充分发展的内流,可以用Inten sity and Hydraulic Diameter (
44、湍流强度与水力直径)方法定义湍流,其中湍流特征长度就是Hydraulic Diameter (水力直径)DH。2) 对于导向叶片或分流板下游的流场,可以用Intensity and Hydraulic Diameter (湍流强度与水力直径)定义湍流,并在Hydrauli Diameter (水力直径)中将导向叶片或分流板的开口部分的长度L定义为特征长度。Intensity and3)如果进口处的流动为受到壁面限制且带有湍流边界层的流动,可以在Length Scale面板中用边界层厚度99通过公式l 0.4 99计算得到湍流长度尺度l。最后在Turbulenee Length Scale (湍
45、流长度尺度)中输入 l的值。c.湍流动能和湍流耗散率在使用各种k模型对湍流进行计算时,需要给定进口边界上的湍流动能k( TurbulentKinetic Energy )和湍流耗散率(Turbulent Dissipation Rate )的估算值。目前没有理论上的精确计算这两个参数的公式,只能通过试验得到。但不可能对各种各样的流动都去做试验,因此,我们必须借助文献中已有的近似公式来估算(许多商用CFD软件也是这样处理的)。对于没有已知条件的情况, 可根据湍流强度I和特征长度L,由下式粗略估计k和 的分布:32k 2(Uavgl)I 3/203/4 l式中C为湍流模型中的一个经验常数,其值约等
46、于0.09。在没有直接输入湍流动能k和湍流耗散率的情况下,可以用Intensity and HydraulicDiameter (湍流强度与水力直径)或 Intensity and Length Scale (湍流强度与长度尺度)等办 法,利用上述公式确定湍流动能k和湍流耗散率。2.压强进口边界条件( P ressure In let )压强进口边界条件用于定义进口流体的压强,可用于不可压缩和可压缩流动。当进口压强已知,而流动速度或流量未知时,可使用压强进口边界条件。压强进口边界条件也可用于定义外部或非受限流动的“自由边界”。pres sure Inlet£ori& Name
47、inletM&mertufn | Tjia-mal Radiation | Spews | DPM | MultipHase | UDS |Specification, Method intensity 占皿 Hydraulic DiameterTurbuJent Inbensaty (%)Hydraulic Diametfir (m)OKCancelH«lp '定义总压与静压首先道总压(1PPo 72(动量)选项卡内,Refere nee Frame为参考值,有绝对值( Absolute )与Relative to Adjaeent Zone )两个选项供选择,一般
48、保持默认的绝对值,V2P )与静压(po )的关系如下(根据伯努力积分)在 Momentum相对于临近区域值(Gauge Total Pressure (表总压)文本框中输入总压的值。静压在FLUENT中被称为 Supersonic/lnitial Gauge Pressure (超音速/初始表压),如果进 口流动是超音速的或者是准备压强进口边界条件进行计算的初始化工作,则必须定义静压。在流场为亚音速时,FLUENT将忽略Supersonic/lnitial Gauge Pressure (超音速/初始表压)的输入数据,而用驻点参数求出静压。(例如,对于空气,当速度为 100m/S时,可先求出
49、2 1 1.225 1002 6125,静压等于总压减去动压得到。)21其动压为12对于高雷诺数的计算, 由于采用压强入口边界条件,湍流强度并不易求出, 这里给出推荐范围1%5% (当然也可以采用预估的方法估计一下)。3.压强出口边界条件( Pressure Outlet )压强出口边界条件在流场出口边界上定义静压,而静压的值仅在流场为亚音速时使用。如果在出口边界上流场达到超音速,则边界上的压强将从流场内部通过插值得到。其他流场变量均从流场内部通过插值获得。在压强出口边界上还需要定义“回流(backflow )”条件。回流条件是在压强出口边界上出现回流使用的边界条件,这样计算将更容易收敛。FL
50、UENT在压强出口边界条件上可以使用径向平衡条件,同时可以给定预期的流量。有关回流的湍流参数的设置如下图所示,需要说明的是这些湍流参数只有当回流发生时才被使用。但是即使在计算结果中没有回流出现,也应该将出口条件用真实流场的值设定, 这样可以在计算过程中出现回流时加速收敛。Pressure OutletZone NaneoutletMomentum | Ihermal Radiation | Speaes | DPM | MJliphase | UDSoorstantGauge FV ess LIT已(pas call) |qBackflow Direchon Specificabon Meth
51、od formal to BoirdaryI |Ra Jal EquiBarium Prsssure Di£tjibutnnQAverage Pressure SpeciationI llarael Flow FtaterurbulenceSpecifcation Method intensity and Hydraulic DiametervBackflow Turbulent Intensitv)Backflow Hydraulfc Dianne ter (m)4,BtJ|P1P0K Cancel Hep4出流边界条件(Outflow)出流边界条件用于模拟在求解前流速和压力未知的
52、出口边界。在该边界上,用户不需要定义任何内容(除非模拟辐射传热、粒子的离散相及多口出流)。该边界条件适用于出口处的流动是完全发展的情况。 所谓完全发展,意味着出流面上的流动情况由区域内部外推得到, 且对上游流动没有影响。出游边界条件不能用于可压流动,也不能与压力进口边界一起使用(压力进口边界条件可与压力出口边界条件一起使用)(占总出口边界条件的设置比较简单,只需给定所指定的出流边界上流体的流出量权重流出量的百分比)。如果系统只有一个出口,则直接输入“1”即可。注意:在使用出流边界条件时,如果在计算过程中,在出流边界上的任何一点有回流, 计算的收敛性都会受到影响,尤其在进行湍流计算时,这种现象比
53、较明显。这里,可尝试使用压力出口边界条件代替出流边界条件。5.内部界面(in terior )与交界面(in terface )这两种面用于两个区域的交界处,在此界面上不需要用户输入任何内容,只需要指定其位置。一般内部节点单元的表面都默认是内部界面。围的非旋转流体区域)的界面处,将两个区域“隔开” 内容,只需要指定其位置。我们注意到,相接壁面在导入 面生成另外一个相应的shadow wall,将相接壁面由wall 与 shadow wall 合并为 interior 类型。内部界面实际是两个区域公用一个界面(只有一个面)内部界面(interior)边界条件用在两个区域(如水泵中同叶轮一起旋转的
54、流体区域与周 。在该边界上,不需要用户输入任何wall。交界面(in terface)也是两个区需要在flue ntFLUENT 时,FLUENT 会为该壁 更改为interior时,FLUENT会将域的交界,不同的是有两个面成对出现。 中定义。交界面网格划分可以不一样, 这个功能使得划分网格变得容易许多,对于交界面,必须是重合或部分重合, flue nt会自动在重合的部分进行变量的插值和传递。 对于复杂的模型可以采取分块划分。如果可以熟练的将两个区域的网格划分出来,并使得交界处的节点统一(合并节点),那么可以不用采用in terface,直接采用in terior是最好的选择,因为采用这种方式即提高计算速度又提高计算精 度。相反,采用in terface由于插值过程的存在,必然会降低计算速度和精度。注意:in terior界面处的网格必须节点保持一致,in terface可以不用一致,但要保证in terface是成对出现的。6.壁面边界条件(wall )一般来说,壁面边界条件指定后不用修改,但有些情况会做一些相应的设置,下面是一些特殊例子。离心泵旋转的叶轮,因叶轮要随内部流体旋转,故其需设置moving wa
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 元旦节放假通知模板
- 辽宁省朝阳市凌源市2027届六年级数学第一学期期末达标测试试题含解析
- 2027届河北省保定市蠡县蠡吾镇北漳小学三年级数学第一学期期末联考模拟试题含解析
- 河北省张家口市崇礼县2027届六上数学期末统考模拟试题含解析
- 广西防城港市东兴市2027届三上数学期末学业水平测试试题含解析
- 导游个人上半年工作总结
- 张家口市崇礼县2027届三年级数学第一学期期末复习检测模拟试题含解析
- 2026年中国初冷器市场调查研究报告
- 港口项目AAA级安全文明标准化工地创建方案
- 抹灰厚度抽检记录表
- 2025年机场服务人员招聘面试参考题库及答案
- TZDTX 0002-2023 专用铁路企业安全生产标准化建设规范
- (人教A版)必修一高一数学上册第三章:函数的概念与性质重点题型复习(原卷版)
- 2025年广西壮族自治区纪委监委公开遴选公务员笔试试题及答案解析
- 磁共振压脂技术原理与应用
- 口岸建设资金管理办法
- DG-TJ08-2144-2025 公路养护工程质量检验评定标准
- 辅助生殖妇女妊娠管理
- 教师专业发展 课件 第5-9章 教师专业伦理-影响教师专业发展的外部因素
- 堤防工程施工规范
- 《PLC应用项目工单实践教程》课件 模块7 S7-1500系列PLC顺序控制设计法的应用
评论
0/150
提交评论