fluent算法的一些说明

FLUENT-manual中解算方法的一些说明，摘录翻译了其中比较重要的细节，希望对初学FLUENT的朋友在选择设置上提供一些帮助，不致走过多的弯路。)|3H+P0_*S#H5I9P4A/l!Q5w:b5u*c2t离散I7J'h)V2W#_1、QUICK格式仅仅应用在结构化网格上，具有比second-orderupwind更高的精度，当然，FLUENT也允许在非结构网格或者混合网格模型中使用QUICK格式，在这种情况下，非结构网格单元仍然使用second-orderupwind格式计算。2、MUSCL格式可以应用在任何网格和复杂的3维流计算，相比second-orderupwind，third-orderMUSCL可以通过减少数值耗散而提高空间精度，并且对所有的传输方程都适用。third-orderMUSCL目前在FLUENT中没有流态限制，可以计算诸如冲击波类的非连续流场。(c9j&df){#U3、有界中心差分格式boundedcentraldifferencing是LES默认的对流格式，当选择LES后，所有传输方程自动转换为boundedcentraldifferencing。4、lowdiffusiondiscretization只能用在亚音速流计算，并且只适用于implicit-time，对高Mach流，或者在explicittime公式下运行LES，必须使用second-orderupwind。5、改进的HRIC格式相比QUICK与secondorder为VOF计算提供了更高的精度，相比Geo-Reconstruct格式减少更多的计算花费。6、explicittimestepping的计算要求苛刻，主要用在捕捉波的瞬态行为，相比implicittimestepping精度更高，花费更少。但是下列情况不能使用explicittimestepping：/K#g7Y:v0W5L%E3Q&I-t（1）分离计算或者耦合隐式计算。explicittimestepping只能用于耦合显式计算。（2）不可压缩流计算。Explicittimestepping不能用于计算时间精度不可压缩流（如除了理想气体的气体定律）。不可压缩流计算必须在每个时间步迭代至收敛。（3）收敛加速。FASmultigrid与residualsmoothing在explicittimestepping条件下破坏时间精度。7、node-based平均格式比默认的cell-based格式在非结构网格特别是三角形和四面体网格的计算上更精确。)z,n;b6s/s8`分离解算器（solutioncontrols）1、当standardpressure插值格式无效的时候，可以考虑：*H%L)Z'p2Y5J${（1）linear格式，相邻单元的压力平均作为计算面压力。.b.^:_5s0K-s5~!?（2）second-order格式，通过2阶精度对流项重构面压力改进standard与linear格式，但是如果网格质量很差的话，计算会有问题。并且，second-order不适合于多孔介质引起的非连续压力梯度流以及VOF与mixture多相流计算。/F'D,b/\&w9F+L8o&z（3）body-force-weighted通过假设压力和体积力之间差异的标准梯度是常数来计算面压力。如果体积力在动量方程中优先知道的话，如浮力，轴对称旋转流计算，可以获得较好的效果。2、当模型中包含多孔介质，body-force-weighted格式只计算无孔面，并且考虑外体积力（gravity,swirl,Coriolis）以及由于密度的迅速改变而导致的压力梯度（naturalconvection,VOF）的非连续性。所有内部和外部的多孔面按照特定的格式处理，保证法向速度通过单元面的连续性而不管阻力是否连续。2r&H7\;W/v&l3、PRESTO!适用于所有类型的网格，但是对三角形和四面体网格，并不能提供比其他算法更高的精度。4、second-orderupwind与QUICK格式不适用于可压缩多相流中密度的定义。first-orderupwind用于可压缩相的计算，算术平均法用于不可压缩相的计算。由于计算稳定性的原因，推荐在计算可压缩流时，先使用first-order格式，然后转向高精度格式。5、PISO算法的目的是减少SIMPLE与SIMPLEC在求解压力修正方程过程中的反复计算，在每次迭代中需要占用更多的CPU时间，但是可以显著地的减少收敛所需要的迭代步数，特别是针对瞬时问题。#WB&P8G'@1W1q+f-l8f选择解算模式1、segregatedsolver常规上用做不可压缩流和轻微可压缩流计算，coupledsolver最初用做高速可压流的计算。当前这两种解算方式都适用于大范围的流体计算（从不可压到高度可压），但是coupledsolver在高速可压流计算中具备一定的优势。2、默认条件下，FLUENT使用segregatedsolver，但是高速可压流，具有强烈体积力（浮力，旋转力）的耦合流以及非常精细网格的流动，可以选择coupledimplicitsolver，计算中运动方程与能量方程耦合，收敛速度更快但是比segregatedsolver消耗更多的内存。如果电脑内存不足，可以使用segregatedsolver或者coupledexplicitsolver，但是coupledexplicitsolver需要更长的时间达到收敛。,W#M3\%@5M*[7a1L3、下列segregatedsolver物理模型不适用于coupledsolvers：（1）空化模型，VOF模型，混合多相流模型，Eulerian多相流模型*W2q@+G8F%MB&V（2）多孔介质（3）非预混燃烧模型，预混燃烧模型，不完全预混燃烧模型（4）PDF传输模型;v8H'|3j2];C（5）Soot与NOx模型7m;J5o:|,E&m;x.r,U（11）辐射模型（12）熔化/凝固模型（13）壳传导模型$E7|7R6S&U.k7Y（14）操作压力变化5Z,Y*r,G'Z&E（15）周期性流动4、下列不能使用segregatedsolver，必须使用coupledsolvers：$?;x'O.e9b9H2c（1）真实气体模型（非理想气体）8H3w9|*Q!M#A1C（2）自定义的真实气体模型4R!V(L%x3v*Z:C0[+]M（3）非反射边界条件(4)laminarflamesIH7y!F8{/t(x&{(F9Z&A%{'{9d*O$l*Aa离散格式)j#r-^1s(bu6VB,Pfirst-order格式具有较好的收敛性但是精度低，多数情况下计算开始应使用second-order。某些情况下可以先使用first-order然后转为second-order。如果second-order收敛困难，用first-order。在模型简单的结构网格计算上，first-order解算精度与second-order区别不大。)o2{)c!I~+T%A6QQUICK在结构网格条件下计算旋转流场与涡流比second-order提供更高的精度，通常情况下，second-order是足够有效的，使用QUICK不会提高精度。powerlaw也可能用到，但是只有first-order精度。中心差分仅仅用在湍流模型，并且网格必须足够精细而且局部Peclet数要小于1。/[/}/k5B0a%x压力插值格式选择(n.@,w:c)_;}1、如果问题包含大的体积力，推荐使用body-force-weighted。2、高涡流数，高Rayleigh数自然对流，高速旋流，多孔介质以及强烈的曲体流，推荐PRESTO!。&`0w$l4F$st9N3、可压缩流推荐second-order。4、当其他格式无效时使用second-order提高精度。5、second-order不能用在多孔介质以及多相流计算。密度插值格式选择:x%`;^9]2m!om0{1R*V6R多相可压缩流只能使用first-order。+q(p8x9B'j3c:^8m'V计算产生冲击波的可压缩流，强烈推荐在四边形、六面体或者混合网格条件下所有变量包括密度都使用QUICK。压力速度耦合方式选择SIMPLEC相比SIMPLE通常收敛更快，并且可以使用较大的松弛因子。但是比较复杂的湍流或者附加物理方程的计算，SIMPLEC与SIMPLE基本没有区别。强烈推荐在网格扭曲变形较大的情况下使用PISO！。!}$@%i/g3|7vaA(c-i*Z'x*jSIMPLEC与SIMPLE在稳态计算中比PISO！有效，PISO！在非稳态计算中更好。#l,?1l3r'V'm(Q7I7o+t8{7C$m,B网格除low-Reynolds-numberturbulence与LES模型外，湍流计算模式中需要使用wallfunctions，边界层使用粗网格，贴近壁面的网格不能过分精细，保证：30yplus150Fluent经典问题1对于刚接触到FLUENT新手来说，面对铺天盖地的学习资料和令人难读的FLUENThelp，如何学习才能在最短的时间内入门并掌握基本学习方法呢？学习任何一个软件，对于每一个人来说，都存在入门的时期。认真勤学是必须的，什么是最好的学习方法，我也不能妄加定论，在此，我愿意将我三年前入门FLUENT心得介绍一下，希望能给学习FLUENT的新手一点帮助。由于当时我需要学习FLUENT来做毕业设计，老师给了我一本书，韩占忠的《FLUENT流体工程仿真计算实例与应用》，当然，学这本书之前必须要有两个条件，第一，具有流体力学的基础，第二，有FLUENT安装软件可以应用。然后就照着书上二维的计算例子，一个例子，一个步骤地去学习，然后学习三维，再针对具体你所遇到的项目进行针对性的计算。不能急于求成，从前处理器GAMBIT，到通过FLUENT进行仿真，再到后处理，如TECPLOT，进行循序渐进的学习，坚持，效果是非常显著的。如果身边有懂得FLUENT的老师，那么遇到问题向老师请教是最有效的方法，碰到不懂的问题也可以上网或者查找相关书籍来得到答案。另外我还有本《计算流体动力学分析》王福军的，两者结合起来学习效果更好。2CFD计算中涉及到的流体及流动的基本概念和术语：理想流体和粘性流体；牛顿流体和非牛顿流体；可压缩流体和不可压缩流体；层流和湍流；定常流动和非定常流动；亚音速与超音速流动；热传导和扩散等。=61&ID=1411A.理想流体（IdealFluid）和粘性流体（ViscousFluid）：流体在静止时虽不能承受切应力，但在运动时，对相邻的两层流体间的相对运动，即相对滑动速度却是有抵抗的，这种抵抗力称为粘性应力。流体所具备的这种抵抗两层流体相对滑动速度，或普遍说来抵抗变形的性质称为粘性。粘性的大小依赖于流体的性质，并显著地随温度变化。实验表明，粘性应力的大小与粘性及相对速度成正比。当流体的粘性较小（实际上最重要的流体如空气、水等的粘性都是很小的），运动的相对速度也不大时，所产生的粘性应力比起其他类型的力如惯性力可忽略不计。此时我们可以近似地把流体看成无粘性的，这样的流体称为理想流体。十分明显，理想流体对于切向变形没有任何抗拒能力。这样对于粘性而言，我们可以将流体分为理想流体和粘性流体两大类。应该强调指出，真正的理想流体在客观实际中是不存在的，它只是实际流体在某些条件下的一种近似模型。B.牛顿流体（NewtonianFluid）和非牛顿流体（non-NewtonianFluid）：日常生活和工程实践中最常遇到的流体其切应力与剪切变形速率符合下式的线性关系，称为牛顿流体。而切应力与变形速率不成线性关系者称为非牛顿流体。图2-1（a）中绘出了切应力与变形速率的关系曲线。其中符合上式的线性关系者为牛顿流体。其他为非牛顿流体，非牛顿流体中又因其切应力与变形速率关系特点分为膨胀性流体（Dilalant），