搜档网
当前位置:搜档网 › Fluent经典问题答疑

Fluent经典问题答疑

Fluent经典问题答疑
Fluent经典问题答疑

Fluent经典问题及答疑

1 对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢?(#61)

2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。(13楼)

3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不同?(#80)

4 常见离散格式的性能的对比(稳定性、精度和经济性)(#62)

5 在利用有限体积法建立离散方程时,必须遵守哪几个基本原则?(#81)

6 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么?(#130)

7 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难?(#55)

8 什么叫边界条件?有何物理意义?它与初始条件有什么关系?(#56)

9 在一个物理问题的多个边界上,如何协调各边界上的不同边界条件?在边界条件的组合问题上,有什么原则?

10 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别?(#143)

11 在网格生成技术中,什么叫贴体坐标系?什么叫网格独立解?(#35)

12 在GAMBIT的foreground和background中,真实体和虚实体、实操作和虚操作四个之间是什么关系?

13 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量?及其在做网格时大致注意到哪些细节?(#38)

14 画网格时,网格类型和网格方法如何配合使用?各种方法有什么样的应用范围及做网格时需注意的问题?(#169)

15 对于自己的模型,大多数人有这样的想法:我的模型如何来画网格?用什么样的方法最简单?这样做网格到底对不对?(#154)

16 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?(#40)

17 依据实体在GAMBIT建模之前简化时,必须遵循哪几个原则?(#170)

18 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理?b、计算域内的内部边界如何处理(2D)?(#128)

19 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些?(#127)

20 何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT是怎样使用区域的?(#41)

21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?(9楼)

22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?(7楼)

23 在FLUENT运行过程中,经常会出现“turbulence viscous rate”超过了极限值,此时如

何解决?而这里的极限值指的是什么值?修正后它对计算结果有何影响?(#28)

24 在FLUENT运行计算时,为什么有时候总是出现“reversed flow”?其具体意义是什么?有没有办法避免?如果一直这样显示,它对最终的计算结果有什么样的影响?(#29)

25 燃烧过程中经常遇到一个“头疼”问题是计算后温度场没什么变化?即点火问题,解决计算过程中点火的方法有哪些?什么原因引起点火困难的问题? (#183)

26 什么叫问题的初始化?在FLUENT中初始化的方法对计算结果有什么样的影响?初始化中的“patch”怎么理解?(12楼)

27 什么叫PDF方法?FLUENT中模拟煤粉燃烧的方法有哪些?(#197)

28 在利用prePDF计算时出现不稳定性如何解决?即平衡计算失败。

29 在prePDF运行时,报告中会出现提示信息和错误有哪些?并分析错误的原因,提出解决的方法?

30 FLUENT运行过程中,出现残差曲线震荡是怎么回事?如何解决残差震荡的问题?残差震荡对计算收敛性和计算结果有什么影响?(#30)

31 数值模拟过程中,什么情况下出现伪扩散的情况?以及对于伪扩散在数值模拟过程中如何避免?(#31)

32 FLUENT轮廓(contour)显示过程中,有时候标准轮廓线显示通常不能精确地显示其细节,特别是对于封闭的3D物体(如柱体),其原因是什么?如何解决?(#51)

33 如果采用非稳态计算完毕后,如何才能更形象地显示出动态的效果图?(#52)

34 在FLUENT的学习过程中,通常会涉及几个压力的概念,比如压力是相对值还是绝对值?参考压力有何作用?如何设置和利用它?(#32)

35 在FLUENT结果的后处理过程中,如何将美观漂亮的定性分析的效果图和定量分析示意图插入到论文中来说明问题?(#53)

36 在DPM模型中,粒子轨迹能表示粒子在计算域内的行程,如何显示单一粒径粒子的轨道(如20微米的粒子)?(#159)

37 在FLUENT定义速度入口时,速度入口的适用范围是什么?湍流参数的定义方法有哪些?各自有什么不同?(#57)

38 在计算完成后,如何显示某一断面上的温度值?如何得到速度矢量图?如何得到流线?(#95)

39 分离式求解器和耦合式求解器的适用场合是什么?分析两种求解器在计算效率与精度方面的区别。(#58)

40 在处理高速空气动力学问题时,采用哪种耦合求解器效果更好?为什么?(#68)

41 近20多年来,用于超音速流动的湍流模型主要有哪些?各之间模型有什么不同?

42 超音速燃烧反应的模型有哪些?它们有什么特点?

43 FLUENT中常用的文件格式类型:dbs,msh,cas,dat,trn,jou,profile等有什么用处?(#22)

44 在计算区域内的某一个面(2D)或一个体(3D)内定义体积热源或组分质量源。如何把这个zone定义出来?而且这个zone仍然是流体流动的。(#96)(#147)

45 FLUENT进行化学反应计算时模型的选择、求解器的选择以及相关参数的设置需要哪些问题?用FLUENT如何进行化学反应的计算?

46 如何选择单、双精度解算器的选择?(#50)

47 求解器为flunet5/6在设置边界条件时,specify boundary types下的types中有三项关于interior,interface,internal设置,在什么情况下设置相应的条件?它们之间的区别是什么?interior好像是把边界设置为内容默认的一部分;interface是两个不同区域的边界区,比如说离心泵的叶轮旋转区和叶轮出口的交界面;internal;请问以上三种每个的功能?最好能举

一两个例子说明一下,因为这三个都是内部条件吧,好像用的很多。(#48)

48 FLUENT并行计算中Flexlm如何对多个License的管理?(#46)

49 在“solver”中2D 、axisymmetric和axisymmetric swirl如何区别?对于2D和3D各有什么适用范围?(#37)

50 在设置速度边界条件时,提到了“Velocity formulation(Absolute和Relative)”都是指的动量方程的相对速度表示和绝对速度表示,这两个速度如何理解?(#59)

51 对于出口有回流的问题,在出口应该选用什么样的边界条件(压力出口边界条件、质量出口边界条件等)计算效果会更好?(#42)

52 对于不同求解器,离散格式的选择应注意哪些细节?实际计算中一阶迎风差分与二阶迎风差分有什么异同?(#69)

53 对于FLUENT的耦合解算器,对时间步进格式的主要控制是Courant数(CFL),那么Courant数对计算结果有何影响?(#43)

54 在分离求解器中,FLUENT提供了压力速度耦和的三种方法:SIMPLE,SIMPLEC及PISO,它们的应用有什么不同?(#44)

55 对于大多数情况,在选择选择压力插值格式时,标准格式已经足够了,但是对于特定的某些模型使用其它格式有什么特别的要求?

(#60)

56 计算流体力学中在设定初始条件和边界条件的时候总是要先选择一组湍流参数,并给出其初值。如何选择并给出这些初值呢?有什么经验公式或者别的好的办法吗?(#73)57 讨论在数值模拟过程中采用四面体网格计算效果好,还是采用六面体网格更妙呢?(#70)

58 如何将自己用C语言编辑的程序导入到FLUENT中?在利用UDF编写程序时需注意哪些问题?(#157)

59 在UDF中compiled型的执行方式和interpreted型的执行方式有什么不同?(#72)

60 在用gambit的时候,导入pro/e的stp文件后,在消去最短边的时候,有些最短边不能消去,其是空间线段,用面merge的方法和连接点的方法都不行,请问该怎么消去这类短边?(#144)

61 FLUENT help和GAMBIT help能教会我们(特别是刚入门的新手)学习什么基本知识?(#126)

62 FLUENT如何做汽车外流场计算的模拟?并且怎么可以得到汽车的阻力系数和升力系数?(#170)

63 FLUENT模拟飞行器外部流场,最高MA多少时就不准确了?MA达到一定的程度做模拟需注意哪些问题?(#125)

64 在用gambit建模,保存成*.msh文件时总是出现No entity的错误:Continuum Entity fluid does not contain any valid entity and is not written! Boundary Entity wall does not contain any validentity and is not written! 不知道是什么问题?产生的原因是什么?如何解决?(#150) 65 在做燃烧模拟的时候,入口燃料温度定义为蒸发/离解开始时的温度(也就是,为离散相材料指定的蒸发温度“Vaporization Temperature”),这是指水分蒸发温度吗?一般是多少?(#196)

66 在计算煤粉燃烧时遇到这样的问题:

Warning: volatile + combustible fraction for lignite is greater than 1.0shell conduction zones 如何解决?

67 FLUENT控制方程是无因次的还是有因次的?如果是无因次的,怎么无因次的?

68 做飞机设计时,经常计算一些翼型,可是经常出现计算出来的阻力是负值,出现负值究

竟是什么原因,是网格的问题还是计算参数设置的问题?(#71)

69 FLUENT中的Turbulent intensify是如何定义的?该值应该是小于等于100%,可是我的计算中该值达到400%,不知为何?

70 边界条件中湍流强度怎么设置:入口边界条件中的湍流强度和出口边界条件中的回流湍流强度怎么设置?是取默认值10%吗?(#135)

71 关于Injection中的Total Flow rate:injection 选surface,此时选了好几个面(面积不一定完全相同,但颗粒的入口速度相同),那Total Flow Rate 是指几个面的总流量还是某一个面的啊?只能处理完全相同的面吗?(#160)

72 FLUENT中能不能做插值:在ansys中的模型节点坐标和FLUENT中模型的节点坐标不一致,能不能在FLUENT中对所需要的ansys的节点进行插值?

73 在模拟气固两相流时,斯托克斯假设是什么?dpm模型中的粒子流是指什么?

74 大概需要划分100万个左右的单元,且只计算稳态流动,请问这样的问题PC机上算的了吗?如果能算至少需要怎样的计算机配置呢?(#76)

75 在FLUENT中粘压力应该属于压力还是粘性力?

76 GAMBIT划分三维网格后,怎样知道结点数?如何知道总生成多少网格(整个模型)?(#78) (#153)

77 在FLUENT的后处理中可以显示一个管道的。某个标量的。圆截面平均值沿管道轴线(中心线)的变化曲线吗?何显示空间某一点的数值呀(比如某一点温度)?(#138)

78 如何求解密闭容器中温度升高引起的压力升高(氦在一段直径为20mm,长度为100mm 的圆管内。初始温度为5.5K;然后设置一侧管壁温度为定温70K加热。5.5-10K之间,氦的膨胀系数>0.17。)?79 VOF计算时,出现错误:

Error: Too many (4258) VOF sub-timesteps. The velocity field is probably diverging. Pleasecheck the solution, and reduce the time-step if necessary.

Error Object: ()

该怎么改?采用的时间步长并不大,改小也还是这个错误.

80 如何在gambit中输入cad和Pro/e的图形?如何将FLUNET的结果EXPORT成ANSYS 的文件?(#123)

81 入口和出口处的k和epsilon值怎么设置?(#161)

82 是否能用UDF自定义NOx的生成率:燃料NOx的生成,FLUENT中自带的模型外,怎么用过自己UDF定义的模型呢?主要就是添加NOx输运方程中的源项。UDF的帮助文件中有个例子,就是计算热力NOx的,请问怎么样用它来计算燃料NOx,文献中的NOx生成率大多是用质量分数定义的,而UDF中的例子用的却是摩尔浓度啊,请问如何转换啊?

83 如何获得混合流体中每个cell上的各组分的摩尔分数?一混合气体在气道中流动,要获得每个cell中的每个物质的摩尔分数,应当用什么宏啊?比如要获得混合气体每cell的压力,可用C_P(c,t),可每个物质的摩尔分数该如何获得呢?

84 用FLUENT打开以前做好的case,也能显示流场,但是用鼠标中键点击放大的时候出现错误,不停地显示下面信息

Error: Floating point error: invalid number

Error Object: ()

如果试着关闭流场窗口,就不停地显示下面的信息:

Error message from graphics function Show_Selection_Source:

Can't 'Show' - the 'locater' has been deleted

这个问题是什么原因引起的?如何解决?

85 组分定义:请问气道中流体为氢气和氮气,该流体定义时,该定义成混和气体呢,还是

多相气体?(#162)

86 Interface后不能初始化:两网格之间采用INTERFACE,为何在FLUNET初始化时提示:Couldn’tallocate fine level coefficient matrix!是什么原因?(#170)

87 courant数:在模拟高压的流场的时候,迭代的时候总是自动减小其数值,这是什么原因造成的,为什么?怎么修改?(#75)

88 在模拟计算过程中出现定义的步长不同,结果就不同。这个现象与fluent模拟时出现的“网格数不同,结果也不同”现象相似,请问这种一般问题出在哪里?是累计误差的原因吗?(#170)

89 如何在FLUENT6.1.22 在win32 rshd 集群上搭建并行平台?

90 在FLUENT中tmerge与fuse的区别是什么?

91 计算中出现:

Error:>(great-than):invalid argument[2]:wrong type[not a number]

Error Object:nan

是什么意思?还有计算中出现:

turbulent viscosity limited to viscosity ratio of 1.000000e+05 in 1069 cells

有人说只要Solve-control-limits增大可消除,但是这样做不影响计算结果吗?因为帮助上说是不合理才作限制的?

92 关于wall边界条件的设定:在设定wall类型的边界面的时候,有些时候可供选择是:(1)设定热流密度;(2)设定温度;(3)Couple;有时候却是有五个的选择的:(1)热流密度;(2)设定温度;(3)对流换热(convection),(4)radiation;(5)mixed。选择couple的,可是在第二种情况的时候,却没有这个选项,因为是对流换热的问题,如果选择(convection),又要设置什么(free stream)的温度,这些怎么设置呢?

93 粒子轨迹数与粒子流率有关系吗?如果粒子按直径分级,而每个入射点有几组粒子入射(为了模拟粒子的湍流扩散),每组粒子、不同直径的粒子与质量流率有什么关系呢?离散相的质量流率是如何分配到入射的粒子中去的?比如有N个粒子击中某壁面,如何由这N 个粒子的特征(如粒子直径,质量流率,粒子速度等)得出与总的粒子质量流率之间的关系呢?(比如这几个粒子的质量占总的离散项流率的多少)

94 把带网格的几个volume,copy到另一处,但原来split的界面,现在都变成了wall,怎么才能把wall变成内部流体呢?(19楼)

95 可以用左键转动云图,但想用中键拖动其位置时,Fluent显示如下错误信息:

Error message from

graphics function Show_Selection_Source:

Can't 'Show' - the 'locater' has been deleted

这样有什么问题呢? (#122)

96 在利用dpm模型时,射流源可以设成多组分的混合物吗?模拟喷嘴射出炭粉,无论速度多大,射出的距离都是那么远一点?为什么?

97 在udf中,U,V,W代表的速度,分别代表什么方向的,直角坐标还是柱坐标?(#141)

98 Gambit的网格相连问题:如果物体是由两个相连的模型所结合,一个的网格划分比较密、另一个比较稀疏,用Gambit有办法将两个网格密度不同的物体,相连在一起吗?(#49)99利用FLUENT可以模拟出氯化钠溶液析出氯化钠晶体的过程吗?如何可以采用什么样的模型?

100 在FLUENT里定义流体的密度时,定义为不可压理想流体是用在什么地方呀,讲义上说是用于可变密度的不可压流动,不知如何理解?(#121)

101 已经建好的模型,想修改一些尺寸,但不知道顶点的座标,请问如何在gambit中显示

点的座标?(#79)

102 在FLUENT模拟以后用display下的操作都无法显示,不过刚开始用的是好的,然后就不行了,为什么?(#17楼)

103 能否同时设置进口和出口都为压力的边界条件?在这样的边界条件设置情况下发现没有收敛,研究的物理模型只是知道进口和出口的压力,不知道怎么修改才能使其收敛?(#77)

104 在FLUENT计算时,有时候计算时间会特别长,为了避免断电或其它情况影响计算,应设置自动保存功能,如何设置自动保存功能?在非定常计算中读入自动保存文件时如下出现问题:

Writing "F:\propane\16\160575.cas"...

Error: sopenoutputfile&: unable to open file for output

Error Object: "F:\propane\16\160575.cas"

Error: Error writing "F:\propane\16\160575.cas".

Error Object: #f

非定常的,算了一段之后停下来,改天继续算的时候,自动保存的时候出现问题,请问如何解决?(#113)

105 gambit划分时运动部分与静止部分交接面:一个系统的两块,运动部分与静止部分交接部分近似认为没有空隙(无限小,虽然实际上是不可能的),假设考虑做成一个实体,那么似乎要一起运动或静止;假设分开做成两个实体,那么交接处的两个不完全重合的面要设为W ALL还是什么呢,设成W ALL不就不能过流了吗?(#54)

106 在计算模拟中,continuity总不收敛,除了加密网格,还有别的办法吗?别的条件都已经收敛了,就差它自己了,还有收敛的标准是什么?是不是到了一定的尺度就能收敛了,比如10-e5具体的数量级就收敛了。(#134)

107 离散相后处理:如何将浓度单位由(kg/m3)改为(个/m3)?fluent的宏中如何改变浓度单位?

108 想把gambit的图形保存成图片,可是底色总是黑色,怎么改为白色呀。用windows中画图板的反色,好像失真很多。如何处理?(#66)

109 FLUENT中对于稳态流动问题,得到稳定的流场后,在入口处加入示踪剂,计算示踪剂浓度在出口处的变化,得出RTD曲线。加入示踪剂后,求解时,是只解传质方程和能量方程还是所有方程(包括动量方程、连续性方程)一起求解?

如何在FLUENT中加入示踪剂?

110 在分析一个转轮时,想求得转轮的转矩,不知道fluent中有什么方法可以提供该数据。本来想到用叶片上面的压力乘半径,然后做积分运算,但是由于叶片正反壁面统一定义的,即全部定义为wall-rn1,所以分不出方向来了。(#120)

111 如何在gambit中实现坐标轴的变换:有一个三维的网格,想在柱坐标中实现,可是gambit 中一直显示直角坐标?(#97)

112 FLUENT能计算压力脉动吗?有人说:DYNAMIC PRESSURE=0.5*密度*V*V;可是,在我的算例中,进出口分别设为PRESSURE INLET 及PRESURE OUTLET,进口稳态压力24MPa.进出口流速完全相等。计算结果:进口的DYNAMIC PRESSURE=0.5*密度*V*V=5+5E帕,而出口的DYNAMIC PRESSURE=2+4E帕这就很疑惑:到底怎样计算压力脉动,压力脉动=DYNAMIC PRESSURE?(#152)

113 利用vof非稳态求解,结果明显没有收敛的情况下,为什么就开始提示收敛,虽然可以不管它,继续算下去达到收敛。但是求解怎么会提前收敛?(#114)

114 非稳态过程的收敛问题:用FLUENT计算非稳态问题,是不是在计算时必须保证在每

个时间步timestep里都要收敛才行,否则计算结果就不对呢?也就是说,在iteration选项里,max iteration pertime step设为一个值,比如500,就是如果500次迭代后仍未收敛,进入下一步迭代,那对结果会有什么影响。(#172)

115 如何确定燃烧终止或火焰熄灭?1.根据化学反应速率判断,化学反应速率小于某个值,可以认为化学反应终止。但是这个标准怎么取?2.根据温度判断,认为火焰温度温度小于某个值,火焰熄灭了。同样这个“熄灭温度”怎么取?3.根据温度梯度判断,认为温度骤降时,燃烧就终止;同样存在一个标准问题。

116 在Gambit中如何将两个dbs文件到入:把炉膛分成了三个dbs文件,现在想导入两个dbs文件,在Gambit中进行操作,但好象使用open命令就只能open一个dbs文件,请问这要怎么处理?(#118)

117 颗粒轨迹显示问题:用dpm模型进行计算,计算完成后,在particle tracks 中显示颗粒的轨迹,若是style选为line就无法显示,若是选为point就可以显示,这是为什么?

118 使用revolve命令时,给出点、旋转轴以及角度、高度等参数,按理论其轨线应该是光滑的曲线,但放大后发现显示的轨线并不光滑,这是为什么?任何曲线都是用多个线段逼近的,在autocad中有专门的系统变量如isolines等来控制分段数或线框密度,那么在gambit 中是否也有这样的操作呢?

119 用GAMBIT生成网格时要是出现负值怎么办啊?有什么办法可以改正吗,只能将网格重新画吗?

120 scale是把你所画模型中的单位转化为Fluent默认的m,而unite是根据你自己的需要转化单位,也就是把Fluent中默认的m转画为其他的单位,两中方法对计算没有什么影响吗?(#116)

121 GAMBIT处理技巧:两个圆内切产生的尖角那个面如何生成网格质量才比较好?(#115)

122 关于wall zone上的力和力矩的计算,菜单Report/Forces可以计算某一墙区域上的力和力矩。按照帮助中的说明28.3.1:Alphanumeric Reporting-->Forces on Boundaries-->Computing Forces andMoments,力分为压力和粘性力并给出了压力的计算公式,怎么算出粘性力?123 出口的NOx值(ppm)改怎么平均?计算锅炉出口的NOx的值,用面积加权平均还是该用质量加权平均啊?

124 FLUENT中,对计算结果进行面积分比较容易实现,对整个流场的体积分也很简单,但是,如果想对流场中某一部分区域进行体积分,请问如何实现?

125 在定义PDF的时候,如果FLUENT自带的热动力学数据库理没有自己想要的燃料组分,如何对其进行修改呢?

126 能否用fluent能否模拟催化反应?怎样定义边界才能使反应只发生在定义的催化剂表面附近,在气流到达催化剂所在位置之前不让反应发生。

127 用size function作网格时很难控制网格的数目。比如size function的什么条件除了growth rate都没改变,但是原来设置的growth rate为4.5,可是将这个数变大为6或减少为3,网格的数目都比4.5时要多,这是什么原因?(#156)

128 在gambit中对一体积成功的进行了体网格,网格进行了examine mesh,也没有什么问题,可当要进行边界类型(boundary type)的设定时,却发现type 只有node,element_side 两项,没有什么wall,pressure_outlet等。为何无法定义边界?(#16)

129 流固耦合的边界条件:在FLUENT里面,wall与腔体内的流体换热产生换热耦合问题时,wall怎么定义啊?wall在FLUENT的boundary conditions里只有heat flux, temperature, convection, radiation和mixed 5种选项没有coupling?在materials里面改个什么地方就可以了? (#171)

130 用FLUENT来模拟旋流燃烧器射流情况:燃烧器出口应用柱坐标,而整个炉膛使用直角坐标,如何把这两种坐标整合在一起?

131 pdf设置煤粉的低位发热量时,那个值可以自己任意指定吗?设置时按照例子设位3.5e+7时ok,可设置低于3.0e+7的时候,计算pdf时出现

(*)- ERROR RETURNING FROM EQUILIBRIUM SOLVER

(*)- ERROR IN STOICHIOMETRIC MIXTURE FRACTION CALCULATION的情况,这是怎么回事?

132 网格数量和内存之间的关系是什么?(#33)

133 如何在FLUENT中进行密度的选择?(#142)

134-141 (#27)

134 Tecplot数据文件的结构是什么样的?

135 Tecplot怎样创建等高图?

136 怎样用Tecplot创建切片图?

137 怎样用Tecplot创建有影等高图?

138 怎样用Tecplot创建流线图?

139 怎样用Tecplot创建向量图?

140 怎样用Tecplot创建动画?

141 怎样用Tecplot创建X-Y图?

142 什么是多孔介质;在那些方面应用?(#136)

143 有没有介绍多孔介质的专业书籍? (#137)

144 有关多孔介质的边界设定在哪有介绍?

145 有没有计算多孔介质的算例或文章?

146 多孔介质中ppcm是指的什么具体意思?

147 用fluent能否模拟催化反应?怎样定义边界才能使反应只发生在定义的催化剂表面附近,在气流到达催化剂所在位置之前不让反应发生。

148 如何判断那个时刻化学反应已经终止?

149 fluent能够模拟离子反应吗? 比如钡离子与硫酸根离子反应生成硫酸钡之类的。

150 在气固两相流模型中,网格尺寸应该比颗粒画得大,还是小?这两种哪个算得结果更准确些?

151 网格右手规则是什么?

152 用非预混模型计算LPG(液化气,主要是丙烷和丁烷)燃烧,用什么辐射模型合适,P1还是DO或者是其他的。

153 用非预混模型计算LPG燃烧,能量总是不能收敛,怎样改进?

154 y+是结果算出来之后查看的,具体在生成网格时,如何来控制网格满足一定的y+要求呢?(#36)

155 如何区分层流和紊流?以什么为标准来区分呢?从层流过渡到紊流的标准是什么?(#67)

156 解释一下普朗特混合长度和湍流特征长度的区别啊?(#163)

157 gambit不响应键盘的原因。(#82)

158 gambit不能正常启动的原因有哪些?(#45)

159 在fluent中如何设置工作目录?在Gambit中如何设置工作目录?(#39)

160 在计算过程中其它指数都收敛了,就continuity不收敛是怎么回事?在初始化设置中,那些项影响continuity的收敛?(#74)

1.在gambit中对一体积成功的进行了体网格,网格进行了examine mesh,也没有什么问题,可当要进行边界类型(boundary type)的设定时,却发现type 只有node,element_side两项,没有什么wall,pressure_outlet等。为何无法定义边界?

===> 因为没有选择求解器为fluent 5/6

2.在FLUENT模拟以后用display下的操作都无法显示,不过刚开始用的是好的,然后就不行了,为什么?

===> DirectX 控制面板中的“加速”功能禁用即可

3.把带网格的几个volume,copy到另一处,但原来split的界面,现在都变成了wall,怎么才能把wall变成内部流体呢?

===〉直接边界面定义为interior即可

第3题:在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不同?

注:我将原题目的提问顺序进行了修改调整,这样更利于回答。

4.FLUENT中常用的文件格式类型:dbs,msh,cas,dat,trn,jou,profile等有什么用处?

在Gambit目录中,有三个文件,分别是default_id.dbs,jou,trn文件,对Gambit运行save,将会在工作目录下保存这三个文件:default_id.dbs,default_id.jou,default_id.trn。jou文件是gambit命令记录文件,可以通过运行jou文件来批处理gambit命令;

dbs文件是gambit默认的储存几何体和网格数据的文件;

trn文件是记录gambit命令显示窗(transcript)信息的文件;

msh文件可以在gambit划分网格和设置好边界条件之后export中选择msh文件输出格式,该文件可以被fluent求解器读取。

Case文件包括网格,边界条件,解的参数,用户界面和图形环境。

Data文件包含每个网格单元的流动值以及收敛的历史纪录(残差值)。Fluent自动保存文件类型,默认为date和case文件

Profile文件边界轮廓用于指定求解域的边界区域的流动条件。例如,它们可以用于指定入口平面的速度场。

读入轮廓文件,点击菜单File/Read/Profile...弹出选择文件对话框,你就可以读入边界轮廓文件了。

写入轮廓文件,你也可以在指定边界或者表面的条件上创建轮廓文件。例如:你可以在一个算例的出口条件中创建一个轮廓文件,然后在其它算例中读入该轮廓文件,并使用出口轮廓作为新算例的入口轮廓。要写一个轮廓文件,你需要使用Write Profile面板(Figure 1),菜单:File/Write/Profile...

是关于Tecplot软件使用的;在这里给Tecplot新手推荐一个学习的方法:

如果你对Tecplot一点都不熟悉的话,别紧张,没关系的,你可以直接看Tecplot的动画Demo 很快就能入门。

Demo网页:在“开始”->"所有程序“->"Tecplot 10"->"Tutorials"->"Tutorial"

打开网页,里面有Tecplot使用的大部分功能的动画教程,花上一二十分钟,就能对Tecplot 有所了解,并可以作初步的使用了。

这八道问题的答案基本都在里面了。

以下是Demo的内容:

Let's take care of the warning "turbulent viscosity limited to viscosity ratio****" which is not physical. This problem is mainly due to one of the following:

1)Poor mesh quality(i.e.,skewness > 0.85 for Quad/Hex, or skewness > 0.9 for Tri/Tetra elements). {what values do you have?}

2)Use of improper turbulent boudary conditions.

3)Not supplying good initial values for turbulent quantities.

出现这个警告,一般来讲,最可能的就是网格质量的问题,尤其是Y+值的问题;在划分网格的时候要注意,第一层网格高度非常重要,可以使用NASA的Viscous Grid Space Calculator 来计算第一层网格高度;如果这方面已经注意了,那就可能是边界条件中有关湍流量的设置问题,关于这个,本版中已经有专门的帖子进行了讨论,Fluent培训的教程中也有讲到,请大家参考。

回答问题:4 常见离散格式的性能的对比(稳定性、精度和经济性)。

请参考王福军的书《计算流体动力学分析—CFD理论与应用》

5.在利用有限体积法建立离散方程时,必须遵守哪几个基本原则?

具体参考王福军的书《计算流体动力学—CFD软件原理与应用》的第52-54页,这里只作简短介绍。

在利用有限体积法建立离散方程时,必须遵守如下四条基本原则:

1.控制体积界面上的连续性原则;

2.正系数原则;

3.源项的负斜率线性化原则;

4.主系数等于相邻节点系数之和原则。

有限体积法的四条基本原则.pdf

6 .流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么?

答:这个问题的范畴好大啊。简要的说一下个人的理解吧:流场数值求解的目的就是为了得到某个流动状态下的相关参数,这样可以节省实验经费,节约实验时间,并且可以模拟一些不可能做实验的流动状态。主要方法有有限差分,有限元和有限体积法,好像最近还有无网格法和波尔兹曼法(格子法)。基本思路都是将复杂的非线性差分/积分方程简化成简单的代数方程。相对来说,有限差分法对网格的要求较高,而其他的方法就要灵活的多

第7题:可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难?

注:这个问题不是一句两句话就能说清楚的,大家还是看下面的两篇小文章吧,摘自《计算流体力学应用》,读完之后自有体会。

8 什么叫边界条件?有何物理意义?它与初始条件有什么关系?

边界条件与初始条件是控制方程有确定解的前提。

边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。对于任何问题,都需要给定边界条件。

初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况,对于瞬态问题,必须给定初始条件,稳态问题,则不用给定。对于边界条件与初始条件的处理,直接影响计算结果的精度。

在瞬态问题中,给定初始条件时要注意的是:要针对所有计算变量,给定整个计算域内各单元的初始条件;初始条件一定是物理上合理的,要靠经验或实测结果。

10 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别?PS:这个问题想来应该是比较基础的问题,既然没人回答,我就插几句吧;嘿嘿。

我们知道很多描述物理问题的控制方程最终就可以归结为偏微分方程,描述流动的控制方程也不例外。

从数学角度,一般将偏微分方程分为椭圆型(影响域是椭圆的,与时间无关,且是空间内的闭区域,故又称为边值问题),双曲型(步进问题,但依赖域仅在两条特征区域之间),抛物型(影响域以特征线为分界线,与主流方向垂直;具体来说,解的分布与瞬时以前的情况和边界条件相关,下游的变化仅与上游的变化相关;也称为初边值问题);

从物理角度,一般将方程分为平衡问题(或稳态问题),时间步进问题。

两种角度,有这样的关系:椭圆型方程描述的一般是平衡问题(或稳态问题),双曲型和抛物型方程描述的一般是步进问题。

至于具体的分类方法,大家可以参考一般的偏微分方程专著,里面都有介绍。关于各种不同近似水平的流体控制方程的分类,大家可以参考张涵信院士编写《计算流体力学—差分方法的原理与应用》里面讲的相当详细。

三种类型偏微分方程的基本差别如下:

1)三种类型偏微分方程解的适定性(即解存在且唯一,并且解稳定)要求对定解条件有不同的提法;

2)三种类型偏微分方程解的光滑性不同,对定解条件的光滑性要求也不同;

椭圆型和抛物型方程的解是充分光滑的,因此对定解条件的光滑性要求不高。而双曲型方程允许有所谓的弱解存在(如流场中的激波),即解的一阶导数可以不连续,所以对定解条件的光滑性要求很高,这也正是采用有限元法求解双曲型方程困难较多的原因之一。

3)三种类型偏微分方程的影响区域和依赖区域不一样。

在双曲型和抛物型方程所控制的流场中,某一点的影响区域是有界的,可采用步进求解。如对双曲型方程求解时,为了与影响区域的特征一致,采用上风格式比较适宜。而椭圆型方程的影响范围遍及全场,必须全场求解,所采用的差分格式也要采用相应的中心格式。

以上只是一些较为肤浅的概念,如想深入,可参考相关的偏微分方程及数值计算等书籍。

可压缩Euler及Navier-Stokes方程数值解

描述无粘流动的基本方程组是Euler方程组,描述粘性流动的基本方程组是Navier-Stokes 方程组。用数值方法通过求解Euler方程和Navier-Stokes方程模拟流场是计算流体动力学的重要内容之一。由于飞行器设计实际问题中的绝大多数流态都具有较高的雷诺数,这些流动粘性区域很小,由对流作用主控,因此针对Euler方程发展的计算方法,在大多数情况下对Navier-Stokes方程也是有效的,只需针对粘性项用中心差分离散。

用数值方法求解无粘Euler方程组的历史可追溯到20世纪50年代,具有代表性的方法是1952年Courant等人以及1954年Lax和Friedrichs提出的一阶方法。从那时开始,人们发展了大量的差分格式。Lax和Wendroff的开创性工作是非定常Euler(可压缩Navier-Stokes)方程组数值求解方法发展的里程碑。二阶精度Lax-Wendroff格式应用于非线性方程组派生出了一类格式,其共同特点是格式空间对称,即在空间上对一维问题是三点中心格式,在时间上是显式格式,并且该类格式是从时间空间混合离散中导出的。该类格式中最流行的是MacCormack格式。

采用时空混合离散方法,其数值解趋近于定常时依赖于计算中采用的时间步长。尽管由

时间步长项引起的误差与截断误差在数量级上相同,但这却体现了一个概念上的缺陷,因为在计算得到的定常解中引进了一个数值参数。将时间积分从空间离散中分离出来就避免了上述缺陷。常用的时空分别离散格式有中心型格式和迎风型格式。空间二阶精度的中心型格式(一维问题是三点格式)就属于上述范畴。该类格式最具代表性的是Beam-Warming隐式格式和Jameson等人采用的Runge-Kutta时间积分方法发展的显式格式。迎风型差分格式共同特点是所建立起的特征传播特性与差分空间离散方向选择的关系是与无粘流动的物理特性一致的。第一个显式迎风差分格式是由Courant等人构造的,并推广为二阶精度和隐式时间积分方法。基于通量方向性离散的Steger-Warming和Van Leer矢通量分裂方法可以认为是这类格式的一种。该类格式的第二个分支是Godunov方法,该方法在每个网格步求解描述相邻间断(Riemann问题)的当地一维Euler方程。根据这一方法Engquist、Osher和Roe等人构造了一系列引入近似Riemann算子的格式,这就是著名的通量差分方法。

对于没有大梯度的定常光滑流动,所有求解Euler方程格式的计算结果都是令人满意的,但当出现诸如激波这样的间断时,其表现确有很大差异。绝大多数最初发展起来的格式,如Lax-Wendroff格式中心型格式,在激波附近会产生波动。人们通过引入人工粘性构造了各种方法来控制和限制这些波动。在一个时期里,这类格式在复杂流场计算中得到了应用。然而,由于格式中含有自由参数,对不同问题要进行调整,不仅给使用上带来了诸多不便,而且格式对激波分辨率受到影响,因而其在复杂流动计算中的应用受到了一定限制。

另外一种方法是力图阻止数值波动的产生,而不是在其产生后再进行抑制。这种方法是建立在非线性限制器的概念上,这一概念最初由Boris和Book及Van Leer提出,并且通过Harten发展的总变差减小(TVD, Total Variation Diminishing)的重要概念得以实现。通过这一途径,数值解的变化以非线性的方式得以控制。这一类格式的研究和应用,在20世纪80年代形成了一股发展浪潮。1988年,张涵信和庄逢甘利用热力学熵增原理,通过对差分格式修正方程式的分析,构造了满足熵增条件能够捕捉激波的无波动、无自由参数的耗散格式(NND格式)。该类格式在航空航天飞行器气动数值模拟方面得到了广泛应用。

1987年,Harten和Osher指出,TVD格式最多能达到二阶精度。为了突破这一精度上的限制引入了实质上无波动(ENO)格式的概念。该类格式“几乎是TVD”的,Harten因此推断这些格式产生的数值解是一致有界的。继Harten和Osher之后,Shu和Osher将ENO格式从一维推广到多维。J.Y.Yang在三阶精度ENO差分格式上也做了不少工作。1992年,张涵信另辟蹊径,在NND格式的基础上,发展了一种能捕捉激波的实质上无波动、无自由参数的三阶精度差分格式(简称ENN格式)。1994年,Liu、Osher和Chan发展了WENO(Weighted Essentially Non-Oscillatory)格式。WENO格式是基于ENO格式构造的高阶混合格式,它在保持了ENO格式优点的同时,计算流场中虚假波动明显减少。此后,Jiang提出了一种新的网格模板光滑程度的度量方法。目前高阶精度格式的研究与应用是计算流体力学的热点问题之一。

不可压缩Navier-Stokes方程求解

不可压缩流体力学数值解法有非常广泛的需求。从求解低速空气动力学问题,推进器内部流动,到水动力相关的液体流动以及生物流体力学等。满足这么广泛问题的研究,要求有与之相应的较好的物理问题的数学模型以及鲁棒的数值算法。

相对于可压缩流动,不可压缩流动的数值求解困难在于,不可压缩流体介质的密度保持常数,而状态方程不再成立,连续方程退化为速度的散度为零的方程。由此,在可压缩流动

的计算中可用于求解密度和压力的连续方程在不可压缩流动求解中仅是动量方程的一个约束条件,由此求解不可压缩流动的压力称为一个困难。求解不可压缩流动的各种方法主要在于求解不同的压力过程。

目前,主要有两类求解不可压缩流体力学的方法,原始变量方法和非原始变量方法。求解不可压缩流动的原始变量方法是将Navier-Stokes方程写成压力和速度的形式,进行直接求解,这种形式已被广为应用。非原始变量方法主要有Fasel提出的流函数-涡函数法、Aziz 和Hellums提出的势函数-涡函数方法。在求解三维流动问题时,上述每一个方法都需要反复求解三个Possion方程,非常耗时。原始变量方法可以分为三类:第一种方法是Harlow 和Welch首先提出的压力Possion方程方法。该方法首先将动量方程推进求得速度场,然后利用Possion方程求解压力,这一种方法由于每一时间步上需要求解Possion方程,求解非常耗时。第二种方法是Patanker和Spalding的SIMPLE(Semi-Implicit Method for Pressure-Linked Equation)法,它是通过动量方程求得压力修正项对速度的影响,使其满足速度散度等于零的条件作为压力控制方程。第三种方法是虚拟压缩方法,这一方法是Chorin 于1967年提出的。该方法的核心就是通过在连续方程中引入一个虚拟压缩因子,再附加一项压力的虚拟时间导数,使压力显式地与速度联系起来,同时方程也变成了双曲型方程。这样,方程的形式就与求解可压缩流动的方程相似,因此,许多求解可压缩流动的成熟方法都可用于不可压缩流动的求解。

目前,由于基于求解压力Possion方程的方法非常复杂及耗时,难于求解具体的工程实际问题,因此用此方法解决工程问题的工作越来越少。工程上常用的主要是SIMPLE方法和虚拟压缩方法。

10. 椭圆形方程求解区中各点上的值是相互影响的,因而各节点上的代数方程必须联立求解,而不能先解得区域中某一部分上的值后再去确定其余地区上的值

抛物型方程描写物理上的一类步进问题,这类问题中因变量与时间有关或问题中有类似时间的变量,因而又称初值问题。其求解区域是个开区间,计算时从已知的初值出发,逐步向前推进,依次获得适合于给定边界条件的解.

双曲型方程数值求解也是一个步进过程。物理学中的波动过程,空气动力学中的无粘流体稳态超音速流动及无粘流体的非稳态流动,都是双曲型问题。

第11题:关于贴体网格的定义,随便找本CFD的书上面基本上都有介绍,在此就不多言。在此谈谈网格独立性的问题。

可以参考帖子:https://www.sodocs.net/doc/3417026085.html,/dvbbs/dispbbs.asp?boardID=61&ID=1208&page=5

数值计算的与实验值之间的误差来源只要有这几个:物理模型近似误差(无粘或有粘,定常与非定常,二维或三维等等),差分方程的截断误差及求解区域的离散误差(这两种误差通常统称为离散误差),迭代误差(离散后的代数方程组的求解方法以及迭代次数所产生的误差),舍入误差(计算机只能用有限位存储计算物理量所产生的误差)等等。在通常的计算中,离散误差随网格变细而减小,但由于网格变细时,离散点数增多,舍入误差也随之加大。

第13题:具体内容可以参考Gambit Documentation中的Quality Type Definitions章节。

判断网格质量的方面有:

Area单元面积,适用于2D单元,较为基本的单元质量特征。

Aspect Ratio长宽比,不同的网格单元有不同的计算方法,等于1是最好的单元,如正三角

形,正四边形,正四面体,正六面体等;一般情况下不要超过5:1.

Diagonal Ratio对角线之比,仅适用于四边形和六面体单元,默认是大于或等于1的,该值越高,说明单元越不规则,最好等于1,也就是正四边形或正六面体。

Edge Ratio长边与最短边长度之比,大于或等于1,最好等于1,解释同上。

EquiAngle Skew通过单元夹角计算的歪斜度,在0到1之间,0为质量最好,1为质量最差。最好是要控制在0到0.4之间。

EquiSize Skew通过单元大小计算的歪斜度,在0到1之间,0为质量最好,1为质量最差。2D质量好的单元该值最好在0.1以内,3D单元在0.4以内。

MidAngle Skew通过单元边中点连线夹角计算的歪斜度,仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。

Size Change相邻单元大小之比,仅适用于3D单元,最好控制在2以内。

Stretch伸展度。通过单元的对角线长度与边长计算出来的,仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。

Taper锥度。仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。V olume单元体积,仅适用于3D单元,划分网格时应避免出现负体积。

Warpage翘曲。仅适用于四边形和六面体单元,在0到1之间,0为质量最好,1为质量最差。

以上只是针对Gambit帮助文件的简单归纳,不同的软件有不同的评价单元质量的指标,使用时最好仔细阅读帮助文件。

另外,在Fluent中的窗口键入:grid quality 然后回车,Fluent能检查网格的质量,主要有以下三个指标:

1.Maxium cell squish: 如果该值等于1,表示得到了很坏的单元;

2.Maxium cell skewness: 该值在0到1之间,0表示最好,1表示最坏;

3.Maxium 'aspect-ratio': 1表示最好。

以上的一些只是简略提要,具体的请参考相关资料。

15 对于自己的模型,大多数人有这样的想法:我的模型如何来画网格?用什么样的方法最简单?这样做网格到底对不对?

我一般首先关注的是我的模型中那一部分的流态情况最复杂,需要加密网格;那一部分得流动情况是我最关注的,需要加密网格。

之后根据上述情况开始考虑怎样对模型进行分块:分块的原则,一方面要保证流态复杂得部位要独立进行网格加密,流动平稳的区域要相应加大网格以节省资源,另一方面,要进行保证分块后的各个子模块结构比较规整,尽可能采用结构网格,提高计算精度。

最后,进行网格划分,开始关注网格划分得顺序,要先对那些线和面进行网格设定约束,保证各个区域满足不同得网格疏密,之后再对各个不同的体积分块进行划分,我一般采用从内到外,从最关注得体到非关注得体的顺序进行划分

16在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?

这个问题就是非连续性网格的设置,一般来说就是把两个交接面设置为一对interface。另外,作此操作可能出现的问题及可供参考的解决方法为:

问题:把两个面(其中一个实际是由若干小面组成,将若干小面定义为了group了)拼接在一起,也就是说两者之间有流体通过,两个面个属不同的体,网格导入到fluent时,使用interface 时出现网格check的错误,将interface的边界条件删除,就不会发生网格检查的错误,如何

将两个面的网格相连?

原因:interface后的两个体的交接面,fluent以将其作为内部流体处理(非重叠部分默认为wall,合并后网格会在某些地方发生畸变,导致合并失败,也可能准备合并的两个面几何位置有误差,应该准确的在同一几何位置(合并的面大小相等时),在合并之前要合理分块。

解决方法:为了避免网格发生畸变(可能一个面上的网格跑到另外的面上了),可以一面网格粗,一面网格细避免;再者就是通过将一个面的网格直接映射到另一面上的,两个面默认为interior.也可以将网格拼接一起.

17 依据实体在GAMBIT建模之前简化时,必须遵循哪几个原则?

答:最根本的原则就是简化后对实际流动影响不大

我觉得在建模前首先要考虑你模型的结构,物理意义上模型是否为轴对称结构或对称结构,如果是的话看能否简化为二维问题,因为二维问题不管从建模上还是求解上都远远方便与三维模型,而且也能达到相应的精度。

其次,在有些梯度比较大的地方这些问题不能简化,像有很多拐角的地方往往存在一些集中,这些不能忽略。

18 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理?b、计算域内的内部边界如何处理(2D)?

答:gambit默认为wall,一般情况下可以到fluent再修改边界类型。内部边界如果是split 产生的,那么就不需再设定了,如果不是,那么就需要设定为interface或者是internal

19 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些?

答:要得到一个问题的定解就需要有定解条件,而边界条件就属于定解条件。也就是说边界条件确定了结果。不同的流体介质有不同的物理属性,也就会得到不同的结果,所以必须指定区域类型。对于gambit来说,默认的区域类型是fluid,所以一般情况下不需要再指定了。

第20题:何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT是怎样使用区域的?

Fluid Zone是一个单元组,是求解域内所有流体单元的综合。所激活的方程都要在这些单元上进行求解。向流体区域输入的信息只是流体介质(材料)的类型。对于当前材料列表中没有的材料,需要用户自行定义。注意,多孔介质也当作流体区域对待。

Solid Zone也是一个单元组,只不过这组单元仅用来进行传热计算,不进行任何的流动计算。作为固体处理的材料可能事实上是流体,但是假定其中没有对流发生,固体区域仅需要输入材料类型。

Fluent中使用Zone的概念,主要是为了区分分块网格生成,边界条件的定义等等;

21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?

可以采用残差控制面板来显示;或者采用通过某面的流量控制;如监控出口上流量的变化;采用某点或者面上受力的监视;涡街中计算达到收敛时,绕流体的面上受的升力为周期交变,而阻力为平缓的直线。

怎样判断计算结果是否收敛?

1、观察点处的值不再随计算步骤的增加而变化;

2、各个参数的残差随计算步数的增加而降低,最后趋于平缓;

3、要满足质量守恒(计算中不牵涉到能量)或者是质量与能量守恒(计算中牵涉到能量)。特别要指出的是,即使前两个判据都已经满足了,也并不表示已经得到合理的收敛解了,因为,如果松弛因子设置得太紧,各参数在每步计算的变化都不是太大,也会使前两个判据得到满足。此时就要再看第三个判据了。

还需要说明的就是,一般我们都希望在收敛的情况下,残差越小越好,但是残差曲线是全场求平均的结果,有时其大小并不一定代表计算结果的好坏,有时即使计算的残差很大,但结果也许是好的,关键是要看计算结果是否符合物理事实,即残差的大小与模拟的物理现象本身的复杂性有关,必须从实际物理现象上看计算结果。比如说一个全机模型,在大攻角情况下,解震荡得非常厉害,而且残差的量级也总下不去,但这仍然是正确的,为什么呢,因为大攻角下实际流动情形就是这样的,不断有涡的周期性脱落,流场本身就是非定常的,所以解也是波动的,处理的时候取平均就可以呢:)

22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?

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可能会得到更好的结果,尤其是可以应用增加的亚松驰迭代时,具体介绍如下:

对于相对简单的问题(如:没有附加模型激活的层流流动),其收敛性已经被压力速度耦合所限制,你通常可以用SIMPLEC算法很快得到收敛解。在SIMPLEC中,压力校正亚松驰因子通常设为1.0,它有助于收敛。但是,在有些问题中,将压力校正松弛因子增加到 1.0可能会导致不稳定。对于所有的过渡流动计算,强烈推荐使用PISO算法邻近校正。它允许你使用大的时间步,而且对于动量和压力都可以使用亚松驰因子1.0。对于定常状态问题,具有邻近校正的PISO并不会比具有较好的亚松驰因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。当你使用PISO邻近校正时,对所有方程都推荐使用亚松驰因子为1.0或者接近1.0。如果你只对高度扭曲的网格使用PISO倾斜校正,请设定动量和压力的亚松驰因子之和为1.0比如:压力亚松驰因子0.3,动量亚松驰因子0.7)。如果你同时使用PISO的两种校正方法,推荐参阅PISO邻近校正中所用的方法。

24 这个问题的意思是出现了回流,这个问题相对于湍流粘性比的警告要宽松一些,有些case 可能只在计算的开始阶段出现这个警告,随着迭代的计算,可能会消失,如果计算一段时间之后,警告消失了,那么对计算结果是没有什么影响的,如果这个警告一直存在,可能需要作以下处理:

1.如果是模拟外部绕流,出现这个警告的原因可能是边界条件取得距离物体不够远,如果边界条件取的足够远,该处可能在计算的过程中的确存在回流现象;对于可压缩流动,边界最好取在10倍的物体特征长度之处;对于不可压缩流动,边界最好取在4倍的物体特征长度之处。

2.如果出现了这个警告,不论对于外部绕流还是内部流动,可以使用pressure-outlet边界条件代替outflow边界条件改善这个问题。

25 燃烧过程中经常遇到一个“头疼”问题是计算后温度场没什么变化?即点火问题,解决计算过程中点火的方法有哪些?什么原因引起点火困难的问题?

这个问题确实比较常见,大概说一下自己的解决方法。

1)点火问题。一般来说,对简单的问题我采取初始化时给一个较高温度的办法,当然,PATCH也是很好的办法。对那些希望局部燃烧的问题,我是画网格的时候予以区分,然后关闭反应项。

2)不着火的问题。原因比较复杂,可能是点火不成功,也可能是其他原因。对一个甲烷燃烧问题,或许可以修改活化能来解决。另外就是对具体问题需要选合适的燃烧模型。

计算后温度场没有变化,说明化学反应根本没有发生!

fluent解决计算过程中的点火问题,是采用"patch"的方法在已经初始化的流场中设定一个高温的点火区,该选项中设定的温度值一定要高于燃料的最低点火温度(例如甲烷最低点火温度为800K,这里就可以将其设定为2000K或3000K),其余参数都可以和初始化时一样。不要担心点火区的高温会对整个计算的影响!

在燃料点火和燃烧中,最关键的是燃烧反应机理,fluent中可以采用一步或两步总体反应(fluent自带)或通过import导入正确的详细反应机理。在具体设置过程中,一定要注意:

在species Model中,Options选项中共包含以下五个选项:

inlet Diffusion;

Diffusion Energy Source,;

Full Multicomponent Diffusion;

Thermal Diffusion;

KINetics from reaction Design。

对于没有专门购买KINetics软件包的朋友来说,只能选前面四项,即通过fluent自带的化学动力学程序来计算燃烧反应机理。如果选了最后一项,而你又没有专门购买Reaction Design的化学动力学软件包,那在计算过程中燃烧反应一定不会发生,也就是不能点火、燃烧,计算后温度场自然就不会变化!(我最初也郁闷怎么就不能点火、燃烧呢!后来我的一个朋友提醒我,才搞定)

26 什么叫问题的初始化?在FLUENT中初始化的方法对计算结果有什么样的影响?初始化中的“patch”怎么理解?

问题的初始化就是在做计算时,给流场一个初始值,包括压力、速度、温度和湍流系数等。理论上,给的初始场对最终结果不会产生影响,因为随着跌倒步数的增加,计算得到的流场会向真实的流场无限逼近,但是,由于Fluent等计算软件存在像离散格式精度(会产生离散误差)和截断误差等问题的限制,如果初始场给的过于偏离实际物理场,就会出现计算很难收敛,甚至是刚开始计算就发散的问题。因此,在初始化时,初值还是应该给的尽量符合实际物理现象。这就要求我们对要计算的物理场,有一个比较清楚的理解。

初始化中的patch就是对初始化的一种补充,比如当遇到多相流问题时,需要对各相的参数进行更细的限制,以最大限度接近现实物理场。这些就可以通过patch来实现,patch 可以对流场分区进行初始化,还可以通过编写简单的函数来对特定区域初始化

27 什么叫PDF方法?FLUENT中模拟煤粉燃烧的方法有哪些?

PDF方法为概率密度函数法

FLUENT中模拟煤粉燃烧的方法有非预混燃烧、有限速率等

概率密度函数输运输运方程方法(PDF方法)是近年来逐步建立起来的描述湍流两相流动的新模型方法。所谓的概率密度函数(Probability Density Function,简称PDF)方法是基于湍流场随机性和概率统计描述,将流场的速度、温度和组分浓度等特征量作为随机变量,研究其概率密度函数在相空间的传递行为的研究方法。PDF模型介于统观模拟和细观模拟之间,是从随机运动的分子动力论和两相湍流的基本守恒定律出发,探讨两相湍流的规律,因此可作为发展双流体模型框架内两相湍流模型的理论基础。它实质上是沟通E-L模型和E-E模型的桥梁,可以用颗粒运动的拉氏分析通过统计理论,即PDF方程的积分建立封闭的E-E两相湍流模型。

非预混湍流燃烧过程的正确模拟要求同时模拟混合和化学反应过程。FLUENT 提供了四种反应模拟方法:即有限率反应法、混合分数PDF 法、不平衡(火焰微元)法和预混燃烧法。火焰微元法是混合分数PDF 方法的一种特例。该方法是基于不平衡反应的,混合分数PDF 法不能模拟的不平衡现象如火焰的悬举和熄灭,NOx 的形成等都可用该方法模拟。但由于该方法还未完善,在FLUENT 只能适用于绝热模型。

对许多燃烧系统,辐射式主要的能量传输方式,因此在模拟燃烧系统时,对辐射能量的传输的模拟也是非常重要的。在FLUENT 中,对于模拟该过程的模型也是非常全面的。包括DTRM、P-1、Rosseland、DO 辐射模型,还有用WSGG 模型来模拟吸收系数。

在查资料时看到的,可能回答不详细!

30

请参考本版帖子:https://www.sodocs.net/doc/3417026085.html,/dvbbs/dispbbs.asp?boardID=61&ID=1246&page=2 这个问题我也一直在想,看到simwe上有人这么回答的:(顺便把相类似的问题的解答也放在一起,方便大家一起解决这类的问题。)

一. 残差波动的主要原因:1、高精度格式;2、网格太粗;3、网格质量差;4、流场本身边界复杂,流动复杂;5、模型的不恰当使用。

二. 问:在进行稳态计算时候,开始残差线是一直下降的,可是到后来各种残差线都显示为波形波动,是不是不收敛阿?

答:有些复杂或流动环境恶劣情形下确实很难收敛。计算的精度(2 阶),网格太疏,网格质量太差,等都会使残差波动。经常遇到,一开始下降,然后出现波动,可以降低松弛系数,我的问题就能收敛,但如果网格质量不好,是很难的。通常,计算非结构网格,如果问题比较复杂,会出现这种情况,建议作网格时多下些功夫。理论上说,残差的震荡是数值迭代在计算域内传递遭遇障碍物反射形成周期震荡导致的结果,与网格亚尺度雷诺数有关。例如,通常压力边界是主要的反射源,换成OUTFLOW 边界会好些。这主要根据经验判断。所以我说网格和边界条件是主要因素。

三. 1、网格问题:比如流场内部存在尖点等突变,导致网格在局部质量存在问题,影响收敛。

2、可以调整一下courant number,courant number实际上是指时间步长和空间步长的相对关系,系统自动减小courant数,这种情况一般出现在存在尖锐外形的计算域,当局部的流速过大或者压差过大时出错,把局部的网格加密再试一下。

在fluent中,用courant number来调节计算的稳定性与收敛性。一般来说,随着courant number的从小到大的变化,收敛速度逐渐加快,但是稳定性逐渐降低。所以具体的问题,在计算的过程中,最好是把courant number从小开始设置,看看迭代残差的收敛情况,如果收敛速度较慢而且比较稳定的话,可以适当的增加courant number的大小,根据自己具体的问题,找出一个比较合适的courant number,让收敛速度能够足够的快,而且能够保持它的稳定性。

四. 另外,Fluent官方论坛上有一个管理员这么解释的:https://www.sodocs.net/doc/3417026085.html,/forum/viewtopic.php?t=1887&sid=5228995e41f4763b95368f660d1a6 71b

1. Residual pattern in an unsteady simulation will be oscillatory. At each time-step, Fluent will try to converge residual values below the default limits set in the Residual Monitors panel.

2. In a steady simulation, ideally the residuals should not oscillate.

Oscillating residuals may be due to wrong case setup. Following comments should help in reducing convergence trouble.

a. Make sure that grid is of good quality. Please refer the section 6.2.2 of Fluent 6.3 Users Guide for more details on the grid quality. https://www.sodocs.net/doc/3417026085.html,/fluent/doc/ori/html/ug/node155.htm

b. Recheck all the boundary condition.

c. If you are using turbulence model, please check the turbulence parameters on all inlet / outlet

fluent问地训练题目回答

1.pressure based 和density based Coupled会同时求解所有的方程(质量守恒方程、动量守恒方程和能量守恒方程)而不是单个方程求解(方程互相分离)。当速度和压力高度耦合(高压和高速)时应该使用耦合求解,但这样会需要较长的计算时间。 在耦合求解中,能量方程中总是包含组分扩散(Species Diffusion Term)项。 当使用segregated求解时,fluent允许指定固体材料的各项异性传导性。 求解方法主要根据要求解的模型来选择。Segregated方法是基于压力,而coupled求解是基于密度的。这样就使得segregated求解低速流动较好而coupled求解音速/超音速问题较好。我不推荐使用coupled求解所有低于马赫数4的流动(直到基于压力的coupled求解方法出现在下一个fluent版本中)。我曾经用segregated方法求解直到1.5马赫的问题,并且结果很好。但是速度越高,需要的网格就越多(因为segregated趋向于“平滑”波动),所以必须多加注意划分网格。 Coupled方法使用默认设置时往往是比较稳定的。Segregated方法常常对容许极根很敏感。当使用segregated方法求解时,不要提高turbulent viscosity ration limit(除非你根据过去的经验或者你的物理模型有很好的理由超过这个极限,但我从没有听说这样是比较理想的)。不要给压力和温度极限限定的合理的范围(例如Plimits=Pstatic+/-(2*dynamic pressure))来计算适当的温度。 1,pressure based 求解方法在求解不可压流体时,如果我们联立求解从动量方程和连续性方程离散得到的代数方程组,可以直接得到各速度分量及相应的压力值,但是要占用大量的计算内存,这一方法已可以在Fluent6.3中实现,所需内存为分离算法的1.5-2倍,同时Fluent6.3中的压力基耦合求解器也很适合求解带有激波的高速空气动力问题(可压流体),这是一个新变化。本人也在尝试用这个模型模拟一些噪声问题。 2,density based求解方法是针对可压流体设计的,因而更适合于可压流场的计算。以速度分量、密度(密度基)作为基本变量,压力则由状态方程求解。 Pressure-Based Solver是Fluent的优势,它是基于压力法的求解器,使用的是压力修正算法,求解的控制方程是标量形式的,擅长求解不可压缩流动,对于可压流动也可以求解;Fluent 6.3以前的版本求解器,只有Segregated Solver和Coupled Solver,其实也就是Pressure-Based Solver 的两种处理方法; Density-Based Solver是Fluent 6.3新发展出来的,它是基于密度法的求解器,求解的控制方程是矢量形式的,主要离散格式有Roe,AUSM+,该方法的初衷是让Fluent具有比较好的求解可压缩流动能力,但目前格式没有添加任何限制器,因此还不太完善;它只有Coupled的算法;对于低速问题,他们是使用Preconditioning方法来处理,使之也能够计算低速问题。 Density-Based Solver下肯定是没有SIMPLEC,PISO这些选项的,因为这些都是压力修正算法,不会在这种类型的求解器中出现的;一般还是使用Pressure-Based Solver解决问题。 2.连续性方程不收敛是怎么回事? 在计算过程中其它指数都收敛了,就continuity不收敛是怎么回事 这和Fluent程序的求解方法SIMPLE有关。SIMPLE根据连续方程推导出压力修正方法求解压力。由于连续方程中流场耦合项被过渡简化,使得压力修正方程不能准确反映流场的变化,从而导致该方程收敛缓慢。 你可以试验SIMPLEC方法,应该会收敛快些。 3.边界条件对应的一般设定方法 边界条件对应的一般设定方法:

Fluent经典问题及解答

Fluent经典问题及解答 1 对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢?(#61) 2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。(13楼) 3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不同?(#80) 4 常见离散格式的性能的对比(稳定性、精度和经济性)(#62) 5 在利用有限体积法建立离散方程时,必须遵守哪几个基本原则?(#81) 6 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么?(#130) 7 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难?(#55) 8 什么叫边界条件?有何物理意义?它与初始条件有什么关系?(#56) 9 在一个物理问题的多个边界上,如何协调各边界上的不同边界条件?在边界条件的组合问题上,有什么原则? 10 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别?(#143) 11 在网格生成技术中,什么叫贴体坐标系?什么叫网格独立解?(#35) 12 在GAMBIT的foreground和background中,真实体和虚实体、实操作和虚操作四个之间是什么关系? 13 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量?及其在做网格时大致注意到哪些细节?(#38) 14 画网格时,网格类型和网格方法如何配合使用?各种方法有什么样的应用范围及做网格时需注意的问题?(#169) 15 对于自己的模型,大多数人有这样的想法:我的模型如何来画网格?用什么样的方法最简单?这样做网格到底对不对?(#154) 16 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?(#40) 17 依据实体在GAMBIT建模之前简化时,必须遵循哪几个原则?(#170) 18 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理?b、计算域内的内部边界如何处理(2D)?(#128) 19 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些?(#127) 20 何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT是怎样使用区域的?(#41) 21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?(9楼) 22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?(7楼)

Fluent中常见问题

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可能会得到更好的结果,尤其是可以应用增加的亚松驰迭代时,具体介绍如下: 对于相对简单的问题(如:没有附加模型激活的层流流动),其收敛性已经被压力速度耦合所限制,你通常可以用SIMPLEC算法很快得到收敛解。在SIMPLEC中,压力校正亚松驰因子通常设为1.0,它有助于收敛。但是,在有些问题中,将压力校正松弛因子增加到1.0可能会导致不稳定。对于所有的过渡流动计算,强烈推荐使用PISO算法邻近校正。它允许你使用大的时间步,而且对于动量和压力都可以使用亚松驰因子1.0。对于定常状态问题,具有邻近校正的PISO并不会比具有较好的亚松驰因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。当你使用PISO邻近校正时,对所有方程都推荐使用亚松驰因子为1.0或者接近1.0。如果你只对高度扭曲的网格使用PISO 倾斜校正,请设定动量和压力的亚松驰因子之和为1.0比如:压力亚松驰因子0.3,动量亚

Fluent 经典问题

QUICK格式可能产生比二阶精度更好的结果。但是,一般情况下,用二阶精度就已足够,即使使用QUICK格式,结果也不一定好。乘方格式(Power-law Scheme)一般产生与一阶精度格式相同精度的结果。中心差分格式一般只用于大涡模拟,而且要求网格很细的情况。 53 对于FLUENT的耦合解算器,对时间步进格式的主要控制是Courant数(CFL),那么Courant 数对计算结果有何影响? courant number实际上是指时间步长和空间步长的相对关系,系统自动减小courant数,这种情况一般出现在存在尖锐外形的计算域,当局部的流速过大或者压差过大时出错,把局部的网格加密再试一下。 在Fluent中,用courant number来调节计算的稳定性与收敛性。一般来说,随着courant number的从小到大的变化,收敛速度逐渐加快,但是稳定性逐渐降低。所以具体的问题,在计算的过程中,最好是把courant number从小开始设置,看看迭代残差的收敛情况,如果收敛速度较慢而且比较稳定的话,可以适当的增加courant number的大小,根据自己具体的问题,找出一个比较合适的courant number,让收敛速度能够足够的快,而且能够保持它的稳定性。 54 在分离求解器中,FLUENT提供了压力速度耦和的三种方法:SIMPLE,SIMPLEC及PISO,它们的应用有什么不同? 在FLUENT中,可以使用标准SIMPLE算法和SIMPLEC(SIMPLE-Consistent)算法,默认是SIMPLE算法,但是对于许多问题如果使用SIMPLEC可能会得到更好的结果,尤其是可以应用增加的亚松驰迭代时,具体介绍如下: 对于相对简单的问题(如:没有附加模型激活的层流流动),其收敛性已经被压力速度耦合所限制,你通常可以用SIMPLEC算法很快得到收敛解。在SIMPLEC中,压力校正亚松驰因子通常设为1.0,它有助于收敛。但是,在有些问题中,将压力校正松弛因子增加到1.0可能会导致不稳定。 对于所有的过渡流动计算,强烈推荐使用PISO算法邻近校正。它允许你使用大的时间步,而且对于动量和压力都可以使用亚松驰因子1.0。对于定常状态问题,具有邻近校正的PISO 并不会比具有较好的亚松驰因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。 当你使用PISO邻近校正时,对所有方程都推荐使用亚松驰因子为1.0或者接近1.0。如果你只对高度扭曲的网格使用PISO倾斜校正,请设定动量和压力的亚松驰因子之和为1.0比如:压力亚松驰因子0.3,动量亚松驰因子0.7)。如果你同时使用PISO的两种校正方法,推荐参阅PISO邻近校正中所用的方法 55 对于大多数情况,在选择选择压力插值格式时,标准格式已经足够了,但是对于特定的某些模型使用其它格式有什么特别的要求? 压力插值方式的列表只在使用Pressure-based求解器中出现。一般情况下可选择Standard;对于含有高回旋数的流动,高 Rayleigh数的自然对流,高速旋转流动,多孔介质流动,高曲率计算区域等流动情况,选择PRESTO格式;对于可压缩流动,选择Second Order;当然也可以选择Second Order以提高精度;对于含有大体力的流动,选择Body Force Weighted。 注意:Second Order格式不可以用于多孔介质;在使用VOF和Mixture多相流模型时,只能

fluent经典问题整理

网格质量与那些因素有关? 网格质量本身与具体问题的具体几何特性、流动特性及流场求解算法有关。因此,网格质量最终要由计算结果来评判,但是误差分析以及经验表明,CFD计算对计算网格有一些一般性的要求,例如光滑性、正交性、网格单元的正则性以及在流动变化剧烈的区域分布足够多的网格点等。对于复杂几何外形的网格生成,这些要求往往并不可能同时完全满足。例如,给定边界网格点分布,采用Laplace 方程生成的网格是最光滑的,但是最光滑的网格不一定满足物面边界正交性条件,其网格点分布也很有可能不能捕捉流动特征,因此,最光滑的网格不一定是最好的网格。对计算网格的一个最基本的要求当然是所有网格点的Jacobian必须为正值,即网格体积必须为正,其他一些最常用的网格质量度量参数包括扭角(skew angle)、纵横比(aspect ratio、Laplacian)、以及弧长(arc length)等。通过计算、检查这些参数,可以定性的甚至从某种程度上定量的对网格质量进行评判。Parmley等给出了更多的基于网格元素和网格节点的网格质量度量参数。有限元素法关于插值逼近误差估计的理论,实际上也对网格单元的品质给出了基本的规定:即每个单元的内切球半径与外切球半径之,应该是一个适当的,与网格疏密无关的常数。 实体与虚体的区别 在建模中,经常会遇到实体、实面与虚体、虚面,虚体的计算域也可以进行计算并得到所需的结果。那么它们的区别是什么呢? 对于求解是没有任何区别的,只要你能在虚体或者实体上划分你需要的网格。关键是看你网格生成的质量如何,与实体虚体无关。 gambit的实体和虚体在生成网格和计算的时候对于结果没有任何影响,实体和虚体的主要区别有以下几点: 1.实体可以进行布尔运算但是虚体不能,虽然不能进行布尔运算,但是虚体存在merge,split 等功能。 2.实体运算在很多cad软件里面都有,但是虚体是gambit的一大特色,有了虚体以后,gambit 的建模和网格生成的灵活性增加了很多。 3.在网格生成的过程中,如果有几个相对比较平坦的面,你可以把它们通过merge合成一个,这样,作网格的时候,可以节省步骤,对于曲率比较大的面,可能生成的网格质量不好,这时候,你可以采取用split的方式把它划分成几个小面以提高网格质量。 在Fluent中进行非稳态(unsteady)计算时如何设置步长?

FLUENT论坛精华常见问题[1]

湍流与黏性有什么关系? 湍流和粘性都是客观存在的流动性质。 湍流的形成需要一定的条件,粘性是一切流动都具有的。 流体流动方程本身就是具非线性的。 NS方程中的粘性项就是非线性项,当然无粘的欧拉方程也是非线性的。 粘性是分子无规则运动引起的,湍流相对于层流的特性是由涡体混掺运动引起的。 湍流粘性是基于湍流体的parcel湍流混掺是类比于层流体中的分子无规则运动,只是分子无规则运动遥远弱些吧了。不过,这只是类比于,要注意他们可是具有不同的属性。 粘性是耗散的根源,实际流体总是有耗散的。 而粘性是制约湍流的。 LANDAU说,粘性的存在制约了湍流的自由度。 湍流粘性系数和层流的是不一样的,层流的粘性系数基本可认为是常数,可湍流中层流底层中粘性系数很小,远小于层流时的粘性系数;而在过渡区,与之相当,在一个数量级;在充分发展的湍流区,又远大于层流时的粘性系数.这是鮑辛内斯克1987年提出的。 1 FLUENT的初始化面板中有一项是设置从哪个地方开始计算(compute from),选择从不同的边界开始计算有很大的区别吗?该怎样根据具体问题选择从哪里计算呢?比如有两个速度入口A和B,还有压力出口等等,是选速度入口还是压力出口?如果选速度入口,有两个,该选哪个呀?有没有什么原则标准之类的东西? 一般是选取ALL ZONE,即所有区域的平均处理,通常也可选择有代表性的进口(如多个进口时)进行初始化。对于一般流动问题,初始值的设定并不重要,因为计算容易收敛。但当几何条件复杂,而且流动速度高变化快(如音速流动),初始条件要仔细选择。如果不收敛,还应试验不同的初始条件,甚至逐次改变边界条件最后达到所要求的条件。 2 要判断自己模拟的结果是否是正确的,似乎解的收敛性要比那些初始条件和边界条件更重要,可以这样理解吗?也就是说,对于一个具体的问题,初始条件和边界条件的设定并不是唯一的,为了使解收敛,需要不断调整初始条件和边界条件直到解收敛为止,是吗?如果解收敛了,是不是就可以基本确定模拟的结果是正确的呢? 对于一个具体的问题,边界条件的设定当然是唯一的,只不过初始化时可以选择不同的初始条件(指定常流),为了使解的收敛比较好,我一般是逐渐的调节边界条件到额定值("额定值"是指你题目中要求的入口或出口条件,例如计算一个管内流动,要求入口压力和温度为10MPa和3000K,那么我开始叠代时选择入口压力和温度为1MPa和500K(假设,这看你自己问题了),等流场计算的初具规模、收敛的较好了,再逐渐调高压力和温度,经过好几次调节后最终到达额定值10MPa和3000K,这样比一开始就设为10MPa 和3000K收敛的要好些)这样每次叠代可以比较容易收敛,每次调节后不用再初始化即自动调用上次的解为这次的初始解,然后继续叠代。即使解收敛了,这并不意味着就可以基本确定模拟的结果是正确的,还需要和实验的结果以及理论分析结果进行对比分析。 连续性方程不收敛是怎么回事? 在计算过程中其它指数都收敛了,就continuity不收敛是怎么回事 这和Fluent程序的求解方法SIMPLE有关。SIMPLE根据连续方程推导出压力修正方法求解压力。由于连续方程中流场耦合项被过渡简化,使得压力修正方程不能准确反映流场的变化,从而导致该方程收敛缓慢。你可以试验SIMPLEC方法,应该会收敛快些。 边界条件对应的一般设定方法 边界条件对应的一般设定方法: *Genaeral--- pressure inlet;pressure outlet *Compressible flows---mass flow inlet;pressure far-field *Incompressible ---velocity inlet;outflow

Fluent动网格----layering个一个简单实例(作者Snow)

Fluent动网格----layering个一个简单实例我这几天看了点动网格技术方面的东西,在学习过程中发现这方面的例子很少,自己也走了一些弯路。现在还好,弄明白了一些,能够应付现在我的工作。为了让更多学习者快速了解动网格,我打算尽量把我学习心得在这里和大家分享,这里给出一个layering的一个简单例子。 1.Gambit画网格 本例很简单,在Gambit里画一个10*10的矩形,网格间隔为1,也就是有100个网格,具体见下图。都学动网格的人了,不至于这个不会做! 这里需要注意一个问题:设置边界条件的时候,一定要把要移动的边单独设定,本例中一右边界作为移动的边,设成wall就可以,这里再后面需要制定。 2.编写UDF #include "udf.h" #include "unsteady.h" #include "stdio.h" #include "stdlib.h" /************************************************************/ real current_time = 0.0 ; Domain * domain ; Thread * thread ; real NV_VEC( origin ),NV_VEC( force ),NV_VEC( moment ) ; /************************************************************/ DEFINE_CG_MOTION(throttle,dt,vel,omega,time,dtime) { current_time = CURRENT_TIME ; vel[0] = 30; Message("time=%f omega=%f\n",current_time) ; }

FLUENT常见问题

如何区分层流和紊流?以什么为标准来区分呢?从层流过渡到紊流的标准是什么? 答:自然界中的流体流动状态主要有两种形式,即层流laminar和湍流(就是问题中所说的紊流)turbulence.层流是指流体在流动过程中两层之间没有相互混渗,而湍流是指流体不是处于分层流动状态。 对于圆管内流动,雷诺数小于等于2300,管流一定为层流,雷诺数大于等于8000到12000之间,管流一定为湍流,雷诺数大于2300而小于8 000时,流动处于层流与湍流的过渡区。 对于一般流动,在计算雷诺数时,可以用水力半径代替管径。 第40题:在处理高速空气动力学问题时,采用哪种耦合求解器效果更好?为什么? 高速空气动力学问题也属于可压缩流动的范围,在Fluent中原则上,使用Pressure-ba sed和Density-based求解器都可以。从历史根源上讲,基于压力的求解器以前主要用于不可压缩流动和微可压缩流动,而基于密度的求解器用于高速可压缩流动。现在,两种求解器都适用于从不可压到高速可压的很大范围流动,但总的来讲,当计算高速可压缩流动时,基于密度的求解器还是璧基于压力的求解器更有优势,因此,在使用Fluent计算高速可压缩流动时,推荐使用Density-based求解器。 也许有很多人对于Pressure-based和Density-based求解器的原理的认识还不够深,在此稍微介绍一下: 求解Navier-Stokes方程的计算方法根据连续方程的处理方式,可以分为密度法和压力法。不论是密度法还是压力法,速度场都是由动量方程所控制,差别在压力场的确定方法上,密度法是通过连续方程确定密度,再由状态方程换算压力,这一方法多用于可压缩流动,作一定修正后,也可用于低马赫数流动,而这一流动已被看做不可压缩流,但此时精度及鲁棒性都有所降低,对于湍流甚至会失去有效性。密度法的弱点正好是压力法的长处,压力法是通过压力方程或压力修正方程来获得压力场,由于其鲁棒性及有效性,得以广泛使用。该方法原是作为求解不可压缩流动发展起来的,但也可以推广到可压缩流的计算上。这两种方法在求解思路上也有所不同,密度法多用同步求解各变量,而压力法则常为顺序求解各变量。显然顺序求解的一个优势是便于补充方程而无需修改算法程序。 Fluent用户手册上,对于可压缩流动有以下需要注意的策略,在此就不再翻译了,以免曲解原意。 Solution Strategies for Compressible Flows The difficulties associated with solving compressible flows are a result of the hi gh degree of coupling between the flow velocity, density, pressure, and energy. Thi

[整理]fluent经典问题请问双CPU并行计算的效率问题.

fluent 经典问题请问双CPU并行计算的效率问题.txt27信念的力量在于即使身处逆境,亦能帮助你鼓起前进的船帆;信念的魅力在于即使遇到险运,亦能召唤你鼓起生活的勇气;信念的伟大在于即使遭遇不幸,亦能促使你保持崇高的心灵。发信人: rao (绕绕), 信区: NumComp 标题: [合集] 请问双CPU并行计算的效率问题 发信站: BBS 水木清华站 (Mon Jul 7 03:32:43 2003), 站内 ☆─────────────────────────────────────☆ xuzheng (天使暂时离开@_@反方向的钟) 于 (Fri Jul 4 11:03:44 2003) 提到: 大致上只有一个CPU在工作,或者两CPU占有率相当于一个CPU mpich1.2.5+fortran 怎么配置可以使两个CPU同时工作?? BOW ☆─────────────────────────────────────☆ luxz (panda--在热死和冻死边缘挣扎) 于 (Fri Jul 4 11:04:57 2003) 提到: mpirun -np 2 *.exe 【在 xuzheng (天使暂时离开@_@反方向的钟) 的大作中提到: 】 : 大致上只有一个CPU在工作,或者两CPU占有率相当于一个CPU : mpich1.2.5+fortran : 怎么配置可以使两个CPU同时工作?? : BOW ☆─────────────────────────────────────☆ xuzheng (天使暂时离开@_@反方向的钟) 于 (Fri Jul 4 11:06:27 2003) 提到: 不是,你误解了我的意思 再具体点说就是16个节点双CPU的集群,并行计算过程中 每个节点的CPU效率大概只有50%

Fluent经典问题及答疑1

Fluent经典问题及答疑1 1 对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢?(#61) 2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。(13楼) 3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不同?(#80) 4 常见离散格式的性能的对比(稳定性、精度和经济性)(#62) 5 在利用有限体积法建立离散方程时,必须遵守哪几个基本原则?(#81) 6 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么?(#130) 7 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难?(#55) 8 什么叫边界条件?有何物理意义?它与初始条件有什么关系?(#56) 9 在一个物理问题的多个边界上,如何协调各边界上的不同边界条件?在边界条件的组合问题上,有什么原则? 10 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别?(#143) 11 在网格生成技术中,什么叫贴体坐标系?什么叫网格独立解?(#35) 12 在GAMBIT的foreground和background中,真实体和虚实体、实操作和虚操作四个之间是什么关系? 13 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量?及其在做网格时大致注意到哪些细节?(#38) 14 画网格时,网格类型和网格方法如何配合使用?各种方法有什么样的应用范围及做网格时需注意的问题?(#169) 15 对于自己的模型,大多数人有这样的想法:我的模型如何来画网格?用什么样的方法最简单?这样做网格到底对不对?(#154) 16 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?(#40) 17 依据实体在GAMBIT建模之前简化时,必须遵循哪几个原则?(#170) 18 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理?b、计算域内的内部边界如何处理(2D)?(#128) 19 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些?(#127) 20 何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT是怎样使用区域的?(#41) 21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?(9楼) 22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?(7楼) 23 在FLUENT运行过程中,经常会出现“turbulence viscous rate”超过了极限值,此时如何解决?而这里的极限值指的是什么值?修正后它对计算结果有何影响?(#28) 24 在FLUENT运行计算时,为什么有时候总是出现“reversed flow”?其具体意义是什么?有没有

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全攻略(探索阶段)

GAMBIT使用说明 GAMBIT是使用FLUENT进行计算的第一个步骤。在GAMBIT 中我们将完成对计算模型的基本定义和初始化,并输出初始化结果供FLUENT的计算需要。以下是使用GAMBIT的基本步骤。 1.1定义模型的基本几何形状 如左图所示的按钮就是用于构造模型的基本几何形状的。当按下这个按钮时,将出现 如下5个按钮,它们分别是用以定义点、线、面、体的几何形状的。 值得注意的是我们定义这些基本的几何元素的一般是依照以下的顺序: 点——线(两点确定一线)——面(3线以上确定一面)——体(3面以上确定体)对各种几何元素的操作基本方式是:首先选中所要进行的操作,再定义完成操作所要的其他元素,作后点“APPL Y”按钮完成操作。以下不一一重复。 下面我们分别介绍各个几何元素的确定方法: 1.1.1点的操作 对点的操作在按下点操作按钮后进行(其他几何元素的操作也是这样)。点有以下几种主要操作 定义点的位置按钮,按下后出现下面对话框 Coordinate Sys.:用以选择已有坐标系中进行当前操 作的坐标系 Type:可以选择3种相对坐标系为当前坐标系:笛卡 儿坐标、柱坐标、球坐标。 以下通过在Global 中直接输入点的x、y、z值定义点, 注意这里的坐标值是绝对坐标值,而Local中输入的是相 对坐标值,一般我们使用绝对坐标值。 Label:为所定义的点命名。 在完成以上定义后就可以通过进行这个点 的定义,同时屏幕左半部的绘图区中将出现被定义的点。 用关闭此对话框。 查看所有点的几何参数按钮(在以后的操作中也可以查看其他元素的几何参数) 在Vertices栏中选择被查询的点,有两种选择方式(其他几 何元素的选择与此类似): ①按住shift键的同时用鼠标左键取点

fluent常见错误汇总 (1)

Fluent 计算错误汇总 1. .fluent 不能显示图像 在运行fluent 时,导入case 后,检查完grid ,在显示grid 时,总是出现这样的错误 Error message from graphics function Update_Display: Unable to Set OpenGL Rendering Context Error: FLUENT received a fatal signal (SEGMENTATION VIOLATION). Error Object: () 解决办法解决办法:: 右键单击快捷方式,把目标由x:fluent.incntbinntx86fluent.exe 改成: x:fluent.incntbinntx86fluent.exe 2d -driver msw 如果还有三维的,可以再建立一个快捷方式改成: x:fluent.incntbinntx86fluent.exe 3d -driver msw 这就可以直接调用了。如果不是以上原因引起的话,也有可能是和别的软件冲突,如MATLAB 等,这也会使fluent 无法显示图像。 Q1:GAMBIT 安装后无法运行,出错信息是“unable find Exceed X Server” A. GAMBIT 需要装EXCEED 才能用。 gambit 的运行:先运行命令提示符,输入gambit,回车 fluent 的运行:直接在开始-程序-Fluent Inc 里面 Q2:Fluent 安装后无法运行,出错信息是“unable find/open license.dat" A. FLUENT 和GAMBIT 需要把相应license.dat 文件拷贝到FLUENT.INC/license 目录下 Q3:出错信息:运行gambit 时提示找不到gambit 文件? A. FLUENT 和GAMBIT 推荐使用默认安装设置, 安装完GAMBIT 请设置环境变量, 设置办法“开始-程序-FLUENT INC-Set Environment" 另外设置完环境变量需要重启一下,否则仍会提示找不到环境变量。 Q4:使用Fluent 和Gambit 需要注意什么问题? A. 安装好FLUENT 和GAMBIT 最好设置一下用户默认路径 推荐设置办法,在非系统分区建一个目录,如d:\users a) win2k 用户在控制面板-用户和密码-高级-高级,在使用fluent 用户的配置文件 修改本地路径为d:\users,重起到该用户运行命令提示符,检查用户路径是否修改 b) xp 用户,把命令提示符发送到桌面快捷方式,右键单击命令提示符快捷方式 在快捷方式-起始位置加入D:\users,重起检查

fluent过来人经验谈之continuity不收敛的问题

continuity不收敛的问题 (1)连续性方程不收敛是怎么回事? 在计算过程中其它指数都收敛了,就continuity不收敛是怎么回事。 这和fluent程序的求解方法SIMPLE有关。SIMPLE根据连续方程推导出压力修正方法求解压力。由于连续方程中流场耦合项被过渡简化,使得压力修正方程不能准确反映流场的变化,从而导致该方程收敛缓慢。 你可以试验SIMPLEC方法,应该会收敛快些。 在计算模拟中,continuity总不收敛,除了加密网格,还有别的办法吗?别的条件都已经收敛了,就差它自己了,还有收敛的标准是什么?是不是到了一定的尺度就能收敛了,比如10-e5具体的数量级就收敛了 continuity是质量残差,具体是表示本次计算结果与上次计算结果的差别,如果别的条件收敛了,就差它。可以点report,打开里面FLUX选项,算出进口与出口的质量流量差,看它是否小于0.5%.如果小于,可以判断它收敛. (2) fluent残差曲线图中continuity是什么含义? 是质量守恒方程的反映,也就是连续性的残差。这个收敛的快并不能说明你的计算就一定正确,还要看动量方程的迭代计算。表示某次迭代与上一次迭代在所有cells积分的差值,continuty表示连续性方程的残差 (3) 正在学习Fluent,模拟圆管内的流动,速度入口,出口outflow 运行后xy的速度很快就到1e-06了,但是continuity老是降不下去,维持在1e-00和1e-03之间,减小松弛因子好像也没什么变化大家有什么建议吗? 你查看了流量是否平衡吗?在report->flux里面操作,mass flow rate,把所有进出口都选上,compute一下,看看nut flux是什么水平,如果它的值小于总

Fluent经典问题及答疑2

Fluent经典问题及答疑2 51 对于出口有回流的问题,在出口应该选用什么样的边界条件(压力出口边界条件、质量出口边界条件等)计算效果会更好?(#42) 52 对于不同求解器,离散格式的选择应注意哪些细节?实际计算中一阶迎风差分与二阶迎风差分有什么异同?(#69) 53 对于FLUENT的耦合解算器,对时间步进格式的主要控制是Courant数(CFL),那么Courant 数对计算结果有何影响?(#43) 54 在分离求解器中,FLUENT提供了压力速度耦和的三种方法:SIMPLE,SIMPLEC及PISO,它们的应用有什么不同?(#44) 55 对于大多数情况,在选择选择压力插值格式时,标准格式已经足够了,但是对于特定的某些模型使用其它格式有什么特别的要求? (#60) 56 计算流体力学中在设定初始条件和边界条件的时候总是要先选择一组湍流参数,并给出其初值。如何选择并给出这些初值呢?有什么经验公式或者别的好的办法吗?(#73) 57 讨论在数值模拟过程中采用四面体网格计算效果好,还是采用六面体网格更妙呢?(#70) 58 如何将自己用C语言编辑的程序导入到FLUENT中?在利用UDF编写程序时需注意哪些问题?(#157) 59 在UDF中compiled型的执行方式和interpreted型的执行方式有什么不同?(#72) 60 在用gambit的时候,导入pro/e的stp文件后,在消去最短边的时候,有些最短边不能消去,其是空间线段,用面merge的方法和连接点的方法都不行,请问该怎么消去这类短边?(#144) 61 FLUENT help和GAMBIT help能教会我们(特别是刚入门的新手)学习什么基本知识?(#126) 62 FLUENT如何做汽车外流场计算的模拟?并且怎么可以得到汽车的阻力系数和升力系数?(#170) 63 FLUENT模拟飞行器外部流场,最高MA多少时就不准确了?MA达到一定的程度做模拟需注意哪些问题?(#125) 64 在用gambit建模,保存成*.msh文件时总是出现No entity的错误:Continuum Entity fluid does not contain any valid entity and is not written! Boundary Entity wall does not contain any validentity and is not written! 不知道是什么问题?产生的原因是什么?如何解决?(#150) 65 在做燃烧模拟的时候,入口燃料温度定义为蒸发/离解开始时的温度(也就是,为离散相材料指定的蒸发温度“Vaporization Temperature”),这是指水分蒸发温度吗?一般是多少?(#196) 66 在计算煤粉燃烧时遇到这样的问题: Warning: volatile + combustible fraction for lignite is greater than 1.0shell conduction zones 如何解决? 67 FLUENT控制方程是无因次的还是有因次的?如果是无因次的,怎么无因次的? 68 做飞机设计时,经常计算一些翼型,可是经常出现计算出来的阻力是负值,出现负值究竟是什么原因,是网格的问题还是计算参数设置的问题?(#71) 69 FLUENT中的Turbulent intensify是如何定义的?该值应该是小于等于100%,可是我的计算中该值达到400%,不知为何? 70 边界条件中湍流强度怎么设置:入口边界条件中的湍流强度和出口边界条件中的回流湍流强度怎么设置?是取默认值10%吗?(#135) 71 关于Injection中的Total Flow rate:injection 选surface,此时选了好几个面(面积不一定完全相同,但颗粒的入口速度相同),那Total Flow Rate 是指几个面的总流量还是某一个面的啊?只能处理完全相同的面吗?(#160) 72 FLUENT中能不能做插值:在ansys中的模型节点坐标和FLUENT中模型的节点坐标不一致,能

相关主题