Fluent求解参数设置知识分享

Fluent求解参数设置知识分享
Fluent求解参数设置知识分享

F l u e n t求解参数设置

求解参数设置(Solution Methods/Solution Controls):

在设置完计算模型和边界条件后,即可开始求解计算了,因为常会出现求解不收敛或者收敛速度很慢的情况,所以就要根据具体的模型制定具体的求解策略,主要通过修改求解参数来完成。在求解参数中主要设置求解的控制方程、选择压力速度耦合方法、松弛因子、离散格式等。

在VOF模型中,PISO比较适合于不复杂的流体,SIMPLE和SIMPLEC适合于可压缩的流体或者处于封闭域中的流体。

? 求解的控制方程:

在求解参数设置中,可以选择所需要

求解的控制方程。可选择的方程包括

Flow(流动方程)、Turbulence(湍流方

程)、Energy(能量方程)、Volume

Fraction(体积分数方程)等。在求解过程

中,有时为了得到收敛的解,先关闭一

些方程,等一些简单的方程收敛后,再

开启复杂的方程一起计算。

? 选择压力速度耦合方法:

在基于压力求解器中,FLUENT提供了压力速度耦合的4种方法,即SIMPLE、SIMPLEC(SIMPLE.Consistent)、PISO以及Coupled。定常状态计算一般使用SIMPLE或者SIMPLEC方法,对于过渡计算推荐使用PISO方法。PISO方法还可以用于高度倾斜网格的定常状态计算和过渡计算。需要注意的是压力速度耦合只用于分离求解器,在耦合求解器中不可以使用。

在FLUENT中,可以使用标准SIMPLE算法和SIMPLEC算法,默认是SIMPLE算法,但对于许多问题如果使用SIMPLEC可能会得到更好的结果,尤其是可以应用增加的亚松弛迭代时。

对于相对简单的问题(如没有附加模型激活的层流流动),其收敛性可以被压力速度耦合所限制,用户通常可以使用SIMPLEC算法很快得到收敛解。在SIMPLEC算法中,压力校正亚松弛因子通常设为1.0,它有助于收敛,但是,在有些问题中,将压力校正松弛因子增加到1.0可能会导致流动不稳定,对于这种情况,则需要使用更为保守的亚松弛或者使用SIMPLE算法。对于包含湍流或附加物理模型的复杂流动,只要用压力速度耦合做限制,SIMPLEC就会提高收敛性,它通常是一种限制收敛性的附加模拟参数,在这种情况下,SIMPLE 和SIMPLEC会给出相似的收敛速度。

对于所有的过渡流动计算,推荐使用PISO算法邻近校正。它允许用户使用大的时间步,而且对于动量和压力都可以使用亚松弛因子1.0。对于定常状态问题,具有邻近校正的PISO并不会比具有较好的亚松弛因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO 倾斜校正。

当使用PISO邻近校正时,对所有方程都推荐使用亚松弛因子为1.0或者接近1.0。如果只对高度扭曲的网格使用PISO倾斜校正,则要设定动量和压力的亚松弛因子之和为1.0(例如,压力亚松弛因子0.3,动量亚松弛因子0.7)。

? 松弛因子:

在求解过程中,控制变量的变化是很必要的,这就是通过松弛因子来实现的。它控制变量在每次迭代中的变化,也就是说,变量的新值为原值加上变化量乘以松弛因子。

松弛因子可控制收敛的速度和改善收敛

的状况:

?为1,相当于不用松弛因子

?大于1,为超松弛因子人,加快收敛速

?小于1,欠松弛因子,改善收敛的条件

一般来讲,大家都是在收敛不好的时

候,采用一个较小的欠松弛因子。

Fluent里面用的是欠松弛,主要防止两次

迭代值相差太大引起发散。

松弛因子的值在0~1之间,越小表示两次迭代值之间变化越小,也就越稳定,但收敛也就越慢。

使用默认的亚松弛因子开始计算是很好的习惯,对于大多数流动,不需要修改默认亚松弛因子,如果经过4~5步的迭代,残差仍然增长,就需要减小亚松弛因子。压力、动量、k和ε的亚松弛因子默认值分别为0.3、0.7、0.8 和 0.8。对于SIMPLEC格式一般不需要减小压力的亚松弛因子。在密度和温度强烈耦合

的问题中(如相当高的Rayleigh曲数的自然或混合对流流动),应该对温度或密度(所用的亚松弛因子小于1.0 )的亚松弛因子进行设置。当温度和动量方程没有耦合或者耦合较弱时,流动密度是常数,温度的亚松弛因子可以设1.0。对于其他的标量方程,如漩涡、组分、PDF变量,对于某些问题默认的亚松弛因子可能过大,尤其是对于初始计算,可以将松弛因子设为0.8以使收敛更容易。? 离散格式:

当流动与网格对齐时,如使用四边形/六面体网格模拟层流流动,使用一阶精度离散格式是可以接受的。但当流动斜穿网格线时,一阶精度格式将产生明显的离散误差(数值扩散)。因此,对于2D三角形及3D四面体网格,注意要使用二阶精度格式,对复杂流动更是如此。一般来讲,在一阶精度格式下容易收敛,但精度较差,有时,为了加快计算速度,可先在一阶精度格式下计算,然后再转到二阶精度格式下计算。如果使用二阶精度格式遇到难于收敛的情况,则考虑改换一阶精度格式来计算。

对于转动及有旋流的计算,在使用四边形/六面体网格时,具有三阶精度的QUICK格式可能产生比二阶精度更好的结果。但是,一般情况下,用二阶精度就已足够,即使使用QUICK格式,结果也不一定好。乘方格式一般产生与一阶精度格式相同的精度结果。中心差分格式一般只用于大涡模拟模型,而且要求网格很细的情况。

Fluent后处理(DOC)

第四章Fluent后处理 利用FLUENT 提供的图形工具可以很方便的观察CFD 求解结果,并得到满意的数据和图形,用来定性或者定量研究整个计算。本章将重点介绍如何使用这些工具来观察您的计算结果。 1 生成基本图形 在FLUENT中能够方便的生成网格图、等值线图、剖面图,速度矢量图和迹线图等图形来观察计算结果。下面将介绍如何产生这些图形。 一、生成网格图 生成网格或轮廓线视图的步骤 (1)打开网格显示面板 菜单:Display –〉Grid... 图4-1 网格显示对话框 (2)在表面列表中选取表面。点击表面列表下的Outline 按钮来选择所有“外”表面。如果所有的外表面都已经处于选中状态,单击该按钮将使所有外表面处于未选中的状态。点击表面列表下的Interior 按钮来选择所有“内”表面。同样,如果所有的内表面都已经处于选中状态,单击该按钮将使所有内表面处于未选中的状态。 (3)根据需要显示的内容,可以选择进行下列步骤: 1)显示所选表面的轮廓线,在图4-1所示的对话框中进行如下设置:在Options 项选择Edges,在Edge Type 中选择Outline。 2)显示网格线,在Options 选择Edges,在Edge Type 中选择ALL。 3)绘制一个网格填充图形,在Options 选择Faces。显示选中面的网格节点,在Options 选择Nodes。

(4)设置网格和轮廓线显示中的其它选项。 (5)单击Display 按钮,就可以在激活的图形窗口中绘制选定的网格和轮廓线。 二、绘制等值线和轮廓图 生成等值线和轮廓的步骤: 通过图4-2 所示的等值线对话框来生成等值线和轮廓。 菜单:Display –〉Contours... 图4-2 等值线对话框 生成等值线或轮廓的基本步骤如下: (1) 在Contours Of 下拉列表框中选择一个变量或函数作为绘制的对象。首先在上面的列表中选择相关分类;然后在下面的列表中选择相关变量。 (2) 在Surfaces 列表中选择待绘制等值线或轮廓的平面。对于2D情况,如果没有选取任何面,则会在整个求解对象上绘制等值线或轮廓。对于3D情况,至少需要选择一个表面。 (3) 在Levels 编辑框中指定轮廓或等值线的数目。最大数为100。 (4) 如果需要生成一个轮廓视图,请在Option 中选中Draw Profiles 选项。在轮廓选项对话框中(如图4-3),可以如下定义轮廓:

学习fluent(流体常识及软件计算参数设置)

luent 中一些问题 ( 目录 ) 离散化的目的 计算区域的离散及通常使用的网格 控制方程的离散及其方法 各种离散化方法的区别 8 9 10在GAMBIT 中显示的“check 主要通过哪几种来判断其网格的质量?及其在做网格时大 致注意到哪些细节? 11 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克 服这种情况呢? 12在设置GAMBIT 边界层类型时需要注意的几个问题: a 、没有定义的边界线如何处理? b 、计算域内的内部边界如何处理( 2D )? 13 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪 些? 14 20 何为流体区域( fluid zone )和固体区域( solid zone )?为什么要使用区域的概念? FLUENT 是怎样使用区域的? 15 21 如何监视 FLUENT 的计算结果?如何判断计算是否收敛?在 FLUENT 中收敛准则是 如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些 参数?解决不收1 如何入门 2 CFD 2.1 2.2 2. 3 2.4 2.5 2.6 计算中涉及到的流体及流动的基本概念和术语 理想流体( Ideal Fluid )和粘性流体( Viscous Fluid ) 牛顿流体( Newtonian Fluid )和非牛顿流体( non-Newtonian Fluid ) 可压缩流体 ( Compressible Fluid )和不可压缩流体( Incompressible Fluid ) 层流( Laminar Flow )和湍流( Turbulent Flow ) 定常流动( Steady Flow )和非定常流动( Unsteady Flow ) 亚音 速流动 (Subsonic) 与超音速流动( Supersonic ) 热传导( Heat Transfer )及扩散 ( Diffusion ) 2.7 3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常 使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有 什么不 同? 3.1 3.2 3.3 3.4 4 常见离散格式的性能的对比(稳定性、精度和经济性) 5 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是 什 么? 6 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反 而比 可压缩流动有更多的困难? 6.1 可压缩 Euler 及 Navier-Stokes 方程数值解 6.2 不可压缩 Navier-Stokes 方程求解 什么叫边界条件?有何物理意义?它与初始条件有什么关系? 在数值计算中,偏微分方程的 双曲型方程、椭圆型方程、抛物型方程有什么区别? 在网格生成技术中,什么叫贴体坐标 系?什么叫网格独立解?

学习fluent (流体常识及软件计算参数设置)

luent中一些问题----(目录) 1 如何入门 2 CFD计算中涉及到的流体及流动的基本概念和术语 2.1 理想流体(Ideal Fluid)和粘性流体(Viscous Fluid) 2.2 牛顿流体(Newtonian Fluid)和非牛顿流体(non-Newtonian Fluid) 2.3 可压缩流体(Compressible Fluid)和不可压缩流体(Incompressible Fluid) 2.4 层流(Laminar Flow)和湍流(Turbulent Flow) 2.5 定常流动(Steady Flow)和非定常流动(Unsteady Flow) 2.6 亚音速流动(Subsonic)与超音速流动(Supersonic) 2.7 热传导(Heat Transfer)及扩散(Diffusion) 3 在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不 同? 3.1 离散化的目的 3.2 计算区域的离散及通常使用的网格 3.3 控制方程的离散及其方法 3.4 各种离散化方法的区别 4 常见离散格式的性能的对比(稳定性、精度和经济性) 5 流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么? 6 可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难? 6.1 可压缩Euler及Navier-Stokes方程数值解 6.2 不可压缩Navier-Stokes方程求解 7 什么叫边界条件?有何物理意义?它与初始条件有什么关系? 8 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别? 9 在网格生成技术中,什么叫贴体坐标系?什么叫网格独立解? 10 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量?及其在做网格时大致注意到哪些细节? 11 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢? 12 在设置GAMBIT边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理? b、计算域内的内部边界如何处理(2D)? 13 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些? 14 20 何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT是怎样使用区域的? 15 21 如何监视FLUENT的计算结果?如何判断计算是否收敛?在FLUENT中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收

FLUENT中的求解器、算法和离散方法

FLUENT中的求解器、算法和离散方法 作为一个非科班出身的CFD工程师,一开始常常被CFD软件里各种概念搞的晕头转向。最近终于静下心来看了看CFD理论的书,理清了一些概念。就此写一遍博文,顺便整理一下所学内容。 I 求解器: FLUENT中求解器的选择在如下图所示界面中设置: FLUENT中的求解器主要是按照是否联立求解各控制方程来区分的,详见下图:

II 算法: 算法是求解时的策略,即按照什么样的方式和步骤进行求解。FLUENT中算法的选择在如下图所示的界面中设置:

这里简单介绍一下SIMPLE、SIMPLEC、PISO等算法的基本思想和适用范围。 SIMPLE算法:基本思想如前面讲求解器的那张图中解释分离式求解器的例子所示的一样,这里再贴一遍: 1.假设初始压力场分布。 2.利用压力场求解动量方程,得到速度场。 3.利用速度场求解连续性方程,使压力场得到修正。 4.根据需要,求解湍流方程及其他方程 5.判断但前计算是否收敛。若不收敛,返回第二步。 简单说来,SIMPLE算法就是分两步走:第一步预测,第二步修正,即预测-修正。 SIMPLC算法:是对SIMPLE算法的一种改进,其计算步骤与SIMPLE算法相同,只是压力修正项中的一些系数不同,可以加快迭代过程的收敛。 PISO算法:比SIMPLE算法增加了一个修正步,即分三步:第一步预测,第二步修正得到一个修正的场分布,第三步在第二步基础上在进行一侧修正。即预测-修正-修正。PISO算法在求解瞬态问题时有明显优势。对于稳态问题可能SIMPLE 或SIMPLEC更合适。 如果你实在不知道该如何选择,就保持FLUENT的默认选项好了。因为默认选项可以很好解决70%以上的问题,而且对于大部分出了问题的计算来说,也很少是因为算法选择不恰当所致。 III 离散方法: 离散方法是指按照什么样的方式将控制方程在网格节点离散,即将偏微分格式的控制方程转化为各节点上的代数方程组。FLUENT中离散方法的选择在如下图所示的界面中设置:

fluent图形后处理技巧

在图的图的标题栏上右键,先在page setup中选择color,然后选copy to clipboard 就可以了,不用截图。 你可以这样子,没必要colormap一定非得在左边,是吧?如果你的模型是扁长型的话,你可以这样子:在fluent中display>options ,在option panel中的右下角,在colormap alignment 中选bottom。然后在显示的图形界面中将图放大,并将其拖到靠近colormap的地方,再继续我之前帖子中的操作就可以了。 数据可以在显示图形时调整好,然后不要关闭调整好的窗口,连续导入不同的数据进行显示就可以了..或者可以采用tecplot来进行后处理,图片会漂亮些.... File-hardcopy-调整一下即可 不用改,复制到word里背景直接就变成白色了 生成图片使用file下的hardcopy命令,有一个选项是背景色翻转,你虽然看到的是黑色,输出图片背景是白色 的。还有一种方式就是显示也希望是白色背景,使用命令display>set>colors>background 把gambit的背景变成白色 在edit的default的graphic的windows-background-color中把black修改成white,然后modify f luent中默认的图形背景颜色为黑色,这对于要发表的图形很不利,因此很多人希望背景为白色,那么可以使用如下命令:Lf ile-》hardcopy设置格式选择为jpg,color选项之后save那么图形就是希望的白色背景。我发现似乎转化成jpg之后没有运行时候显示的清晰,略微模糊一些,大家可以实验其他设置选择,以求得最好的效果zV>3}D另外可以在控制台命令行输入display/set/color回车之后就显示哪些可以设置的选择,敲进比如background之后就可以改变了,提醒一下单纯改变背景为黑色会使得legnd变成一个梯子,其数字会消失。you should change foreground from white to black .this can be done at he same dislay/set/colors> as the background.p<> 好怎么去掉FLUENT图形显示的黑色背景,一般都建议用抓图后反色背景。另外还有数据显示范围比较小,数据显示相同,色轴没有差别的情况。 本人通过摸索,发现这两个问题可以直接在FLUENT里设置。

辐射和对流模型Fluent参数设置

辐射和对流模型Fluent参数设置 1.读入***.mesh文件,并对网格文件进行进行检查,Grid→cheek,主要看最小体积和最小面积不能为负,之后进行刻度转换,Grid→scale,在Gmbit 里面建模默认尺寸为米,与实际尺寸之间要进行转化,如下图: 2.选择求解器,Define→Models→sover……根据情况选择,如上图:接着选择辐射模型,Define→Models→Radiation,如下图,当Radiation Model面板上 点击ok时,会出现一个信息提示框,告诉你新 的材料物性被添加了,你将在后面设置物性参 数,因此现在只需单击ok确认这个信息即可, 如下图: 注意:当你激活辐射模型后,Fluent会自动打开能量求解器,如下图: 不用再Define→Models→Energy……

3.设置流体粘性,由于模型中空气流速比较大,设成双方程模型:如下图: 4.设置操作条件,此模型此有流体,属有重力情况,Define→Operating Conditions,选中 Gravity.Y方向加速度设置为-9.8 2 m,击OK确定。 /s 设置工作温度,在后面要激活的Boussinesq model要用到,(Boussinesq model:

考虑温度变化而忽略压强变化引起的密度变化叫做Boussinesq 假设) 5. 定义材料并设置其物理属性 Define →Material …… 先定义空气物性,要定义成有浮力的,取Boussinesq 选项。 Density=1.1653/m kg ,()k kg j C p ?=/1005 Thermal Conductivity=0.0267()k m w ?/,Material Type :fluid ; Thermal Expansion Coefficient =0.0033()k /1。 通过滚动条使先前面板中不可见的物性显示出来。在Scattering Coefficient 和Scattering Phase Function 中保持默认值,在要解决的问题中不涉及到散射问题;设定热扩散系数(用boussinesq 模型时)为1e-5K -1。单击Change/Create ,关闭Materials 面板。 6.设置边界条件Define → Boundary Conditions ……

Fluent求解参数设置

求解参数设置(Solution Methods/Solution Controls): 在设置完计算模型和边界条件后,即可开始求解计算了,因为常会出现求解不收敛或者收敛速度很慢的情况,所以就要根据具体的模型制定具体的求解策略,主要通过修改求解参数来完成。在求解参数中主要设置求解的控制方程、选择压力速度耦合方法、松弛因子、离散格式等。 在VOF模型中,PISO比较适合于不复杂的流体,SIMPLE和SIMPLEC适合于可压缩的流体或者处于封闭域中的流体。 ? 求解的控制方程: 在求解参数设置中,可以选择所需要求 解的控制方程。可选择的方程包括Flow(流动方 程)、Turbulence(湍流方程)、Energy(能量方 程)、Volume Fraction(体积分数方程)等。在 求解过程中,有时为了得到收敛的解,先关闭 一些方程,等一些简单的方程收敛后,再开启 复杂的方程一起计算。 ? 选择压力速度耦合方法: 在基于压力求解器中,FLUENT提供了压 力速度耦合的4种方法,即SIMPLE、 SIMPLEC(SIMPLE.Consistent)、PISO以及 Coupled。定常状态计算一般使用SIMPLE或者 SIMPLEC方法,对于过渡计算推荐使用PISO方 法。PISO方法还可以用于高度倾斜网格的定常 状态计算和过渡计算。需要注意的是压力速度 耦合只用于分离求解器,在耦合求解器中不可 以使用。 在FLUENT中,可以使用标准SIMPLE算法和SIMPLEC算法,默认是SIMPLE算法,但对于许多问题如果使用SIMPLEC可能会得到更好的结果,尤其是可以应用增加的亚松弛迭代时。 对于相对简单的问题(如没有附加模型激活的层流流动),其收敛性可以被压力速度耦合所限制,用户通常可以使用SIMPLEC算法很快得到收敛解。在SIMPLEC算法中,压力校正亚松弛因子通常设为1.0,它有助于收敛,但是,在有些问题中,将压力校正松弛因子增加到1.0可能会导致流动不稳定,对于这种情况,则需要使用更为保守的亚松弛或者使用SIMPLE算法。对于包含湍流或附加物理模型的复杂流动,只要用压力速度耦合做限制,SIMPLEC就会提高收敛性,它通常是一种限制收敛性的附加模拟参数,在这种情况下,SIMPLE和SIMPLEC 会给出相似的收敛速度。 对于所有的过渡流动计算,推荐使用PISO算法邻近校正。它允许用户使用大的时间步,而且对于动量和压力都可以使用亚松弛因子1.0。对于定常状态问题,具有邻近校正的PISO并不会比具有较好的亚松弛因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。 当使用PISO邻近校正时,对所有方程都推荐使用亚松弛因子为1.0或者接近1.0。如果只对高度扭曲的网格使用PISO倾斜校正,则要设定动量和压力的亚松弛因子之和为1.0(例如,压力亚松弛因子0.3,动量亚松弛因子0.7)。

FLUENT知识点解析(良心出品必属精品)

一、基本设置 1.Double Precision的选择 启动设置如图,这里着重说说Double Precision(双精度)复选框,对于大多数情况,单精度求解器已能很好的满足精度要求,且计算量小,这里我们选择单精度。然而对于以下一些特定的问题,使用双精度求解器可能更有利[1]。 a.几何特征包含某些极端的尺度(如非常长且窄的管道),单精度求解器可能不能足够精确地表达各尺度方向的节点信息。 b.如果几何模型包含多个通过小直径管道相互连接的体,而某一个区域的压力特别大(因为用户只能设定一个总体的参考压力位置),此时,双精度求解器可能更能体现压差带来的流动(如渐缩渐扩管的无粘与可压缩流动模拟)。 c.对于某些高导热系数比或高宽纵比的网格,使用单精度求解器可能会遇到收敛性不佳或精确度不足不足的问题,此时,使用双精度求解器可能会有所帮助。 [1] 李鹏飞,徐敏义,王飞飞.精通CFD工程仿真与案例实战:FLUENT GAMBIT ICEM CFD Tecplot[M]. 北京,人民邮电出版社,2011:114-116

2.网格光顺化 用光滑和交换的方式改善网格:通过Mesh下的Smooth/Swap来实现,可用来提高网格质量,一般用于三角形或四边形网格,不过质量提高的效果一般般,影响较小,网格质量的提高主要还是在网格生成软件里面实现,所以这里不再用光滑和交换的方式改善网格,其原理可参考《FLUENT全攻略》(已下载)。 3.Pressure-based与Density-based 求解器设置如图。下面说一说Pressure-based和Density-based 的区别: Pressure-Based Solver是Fluent的优势,它是基于压力法的求解器,使用的是压力修正算法,求解的控制方程是标量形式的,擅长求解不可压缩流动,对于可压流动也可以求解;Fluent 6.3以前的版本求解器,只有Segregated Solver和Coupled Solver,其实也是Pressure-Based Solver的两种处理方法;

详细FLUENT实例讲座翼型计算

详细FLUENT实例讲座翼型计算 部门: xxx 时间: xxx 整理范文,仅供参考,可下载自行编辑

CAE联盟论坛精品讲座系列 详细FLUENT实例讲座-翼型计算 主讲人:流沙 CAE联盟论坛总版主 1.1 问题描述 翼型升阻力计算是CFD最常规的应用之一。本例计算的翼型为 RAE2822,其几何参数可以查看翼型数据库。本例计算在来流速度0.75马赫,攻角3.19°情况下,翼型的升阻系数及流场分布,并将计算结果与实验数据进行对比。模型示意图如图1所示。 b5E2RGbCAP 1.p ng(12.13 K>2018/7/29 23:41:251.2 FLUENT前处理设置Step 1:导入计算模型 以3D,双精度方式启动FLUENT14.5。 利用菜单【File】>【Read】>【Mesh…】,在弹出的文件选择对话框中选择网格文件rae2822_coarse.msh,点击OK按钮选择文件。如图2所示。p1EanqFDPw

点击FLUENT模型树按钮General,在右侧设置面板中点击按钮Display…,在弹出的设置对话框中保持默认设置,点击Display按钮,显示网格。如图3所示。DXDiTa9E3d 2.png(11.51 K>2018/7/29 23:41:25

3.png(33.41 K>2018/7/29 23:41:253-2.png(52.04 K>2018/7/29 23:41:25Step 2:检查网格 采用如图4所示步骤进行网格的检查与显示。点击FLUENT模型树节点General节点,在右侧面板中通过按钮Scale…、Check及 Report Quality实现网格检查。 4.png(12. 10 K>RTCrpUDGiT2018/7/29 23:41:25点击按钮Check,在命令输出按钮出现如图5所示网格统计信息。从图中可以看出,网格尺寸分布: x轴:-48.97~50m

离心风机CFD模拟及改进

2005 Fluent 中国用户大会论文集 由于CFD计算可以相对准确地给出流体流动的细节,如速度场、压力场、温度场等特性,因而不仅可以准确预测流体产品的整体性能,而且很容易从对流场的分析中发现产品和工程设计中的问题,所以在国外已经逐步得到广泛的应用。另外,跨学科组合优化设计方法也已经成为复杂叶轮产品的设计平台。 如今,CFD技术运用于风机的实例在我国已不少见,但由于计算机计算能力的限制,模型过于简单。如单独一个离心叶轮的流道或单独算一个蜗壳;或运用一个流道与蜗壳迭代计算的方法研究风机内部流动,上述模型均忽略了由于蜗壳型线的非对称而导致叶轮各叶道流动呈现的非对称流动特征,而且从离心风机通道内流场分析来看,各部件间的相互影响很严 重,所以,必须充分考虑它们之间的相互影响,不能孤立地分别研究[2]。 本文应用Fluent流动分析软件,计算某型号离心通风机全流场,详细得到通风机内部流场流动情况,并根据气动流场,对叶轮前盘 形状和蜗壳出口部位等进行优化设计,同 时,运用多学科优化平台软件OPTIMUS集成流体计算软件FLUENT,优化计算通风机进口型线,比较集成优化型线与单独用Fluent 反复计算的结果,两者基本接近,说明集成优化是可信的。 将流动区域分为三部分:通风机进口部分、叶轮和蜗壳。进口部分和蜗壳是静止元件,叶轮转动,采用gambit进行参数化建模。整个通风机的网格数为80 万,网格采用四面体和六面体混合的非结构网格技术。

气体在通风机内流动时,它的气动性能在很大程度上由它本身的造型决定。由于流道形状、哥氏力和粘性力的影响,通风机内的气体流动十分复杂。一般认为气流在叶轮内的相对运动和在静止元件内的绝对运动为定常流,而且通风机内的气体压强变化不大,可忽略气体的压缩性。因此,通风机内的流动是三维、定常、不可压缩流动。求解相对稳定的、三维不可压缩雷诺平均N-S方程,湍流模型采用标准的εκ?两方程模型,采用一阶迎风格式离散方程,用SIMPLE方法求解控制方程。在OPTIMUS提供的优化算法中,采用序列二次规划算法。 3 数值计算结果与分析 3.1原通风机建模及数值模拟原有离心通风机存在风量不足、风压不均匀等问题,所以首先对原通风机模型进行数值模拟,分析其内部气流流动状况,找出问题所在。图 1 原通风机子午面的速度分布表2 原通风机回转面的速度分布叶轮出口部位的速度 分布很不均匀,在叶轮前半部分,叶轮不出风反而进风,所以此处有较多逆流存在。观察叶轮子午面上速度分布如图1所示,可以看到叶轮出口明显的逆流现象。风机出风口有较多逆流现象,如图2所示。通过上述流场仿真计算,可以确定原通风机的气动性能很不好。分析气动流场,认为性能差的原因主要基于三个方面:1叶片进口部位缺乏导流部分,气体流动的流线不能折转,所以造成叶片前半部分压强低,产生逆流。由于叶轮出口有较多逆流,导致进入蜗壳的气流速度不均匀。2原模型叶轮 采用前向叶片,叶轮的前盘采用平前盘。平前盘制造简单,但对气流的流动情况有不良影响[3]。3通风机蜗壳出口的面积过大,所以在蜗壳出口处压力过低而产生较多的逆流。3.2通风机改型优化计算优化是对通风机改型以得到较好 气动性能的过程。针对原通风机模型气动流场中存在的问题,在结构上作一些相应修改。3.2.1 改进模型A 针对原通风机模型计算中存在的叶片前半部分逆流 现象严重的问题,将叶轮前盘改为弧线型,使计算结果改善。但由于将叶轮的前盘改为弧形,而使叶轮出口宽度减小,所以为了不降低流量,将叶轮的轴向尺寸增加。叶轮出口宽度增加到252mm。如此改动后,叶轮沿子午面速度分布如图3所示。改为 弧形前盘,对气流进行导流,则气动性能改善。与原模型相比,通风机的出口风压增加24.9%,出口流量增加17%,轴功率增加9.4%,效率增加7.6%。2005 Fluent 中国用户大会论文集90 气动性能有所改善,但轴功率增加。通过观察通风机内部气流 的流动情况,叶轮进口部位的流动得到好转,但蜗壳出口部位的流动仍然不好,蜗壳

fluent求解器

Model slover slover:求解器 Pressuere based: 基于压力 Density based:基于密度 Formulation:算法 implicit: 隐式算法 explicit:显式算法 space:选择空间属性 2D:二维空间 Axisymmetric:轴对称空间 Axisymmetric swirl:轴对称旋转空间 3D:三维空间 time:时间 steady:稳态 unsteady:非稳态 (~~~~~~~~~~~~~~~~~) velocity formulation:速度属性 absolute:绝对速度 relative:相对速度 Gradient option:梯度选项 Green-Gauss Cell-Based:格林-高斯基于单元体-默认方法;解有伪扩散(求解域的拖尾现象)。伪扩散是指在平流扩散方程数值解中因平流项有限差分的截断误差引起的虚假扩散。这是解方程欧拉型模式所特有的。其大小与所用的有限差分格式有关,有时甚至完全掩盖方程中其他扩散项的作用。为克服伪扩散,须采取特殊的技术措施和各种不同的差分格式。 Green-Gauss Node-Based:格林-高斯基于节点-更精确;最小化伪扩散;推荐用在三角网格上。 least-quares cell based:基于单元体的最小二乘法插值。推荐用于多面体网格,与基于节点的格林-高斯格式具有相同的 精度和格式。 porous formulation:选择多孔算法 superficial velocity:表面速度算法 physical velocity:物理速度算法 Multiphase Model:多相流模型 Volume of Fluid:VOF模型 Mixture:混合模型 Eulerian:欧拉模型 Energy:能量方程

CFD 的Fluent后处理tecplot软件动画步骤方法

创〗tecplot 中动画制作方法。 [精华] 于 2005-11-09 09:41 个时间序列的数据读入以后利用tecplot 中的tool/Animate/选项可以创建动画。可以根据不同的需要选择contours 、zones 。在应用中一般选择zones 多一点。 主题相关图片如下: dreamoon 发帖: 13 于 2005-11-09 09:46 在zones 里有如下弹出窗口,选择起始zone 和结束zone ,然后输出即可。 此主题相关图片如下:

积分: 0 雪币: 13 dreamoon 发帖: 13 积分: 0 雪币: 13 于 2005-11-09 09:54 或者另外有一种更为方便的方法,该法可以不用一次将所有的数据文件读入,对内存和机子速度较慢的用户更实用: File/Export ,选择avi ,然后打开要输出的contour ,进行如图的操作: 此主题相关图片如下:

dreamoon 编辑于2005-11-09 10:01 dreamoon 发帖: 13 积分: 0 雪币: 13 于2005-11-09 09:56 然后: 此主题相关图片如下:

dreamoon 发帖: 13 积分: 0 雪币: 13 于2005-11-09 09:58 最后选择Finish Animation就可以了。 此主题相关图片如下:

东岸线 发帖: 361 积分: 0 雪币: 310 于2005-11-09 18:58 好 有机会试试

flyboys 发帖: 35 积分: 0 雪币: 35 于2005-11-10 22:02 楼主的数据源是来自 fluent计算获得的数据吧!我们没有用过fluent,根本不知道数据格式是什么?能否把你所作例子的数据格式呢?谢谢 dreamoon 发帖: 13 积分: 0 雪币: 13 于2005-11-11 07:32 我给的例子是一般性的数据;对于fluent 来说就是利用软件的自动编号过程将计算不同时间(或迭代步)的结果保存下来然后分别导入Tecplot就可以了,具体的方法可以参考fluent的帮助手册中关于文件的读写 的相关部分。 wilim 飞燕 发帖: 6 积分: 0 雪币: 6 于2005-11-16 20:49 直接在fluent里面做动画不就可以了,为何还要导出到tecplot中呢,不理解 dreamoon于2005-11-18 00:06

FLUENT求解器设置

FLUENT求解器设置主要包括:1、压力-速度耦合方程格式选择2、对流插值3、梯度插值4、压力插值 下面对这几种设置做详细说明。 一、压力-速度耦合方程求解算法 FLUENT中主要有四种算法:SIMPLE,SIMPLEC,PISO,FSM (1)SIMPLE(semi-implicit method for pressure-linked equations)半隐式连接压力方程方法,是FLUENT的默认格式。 (2)SIMPLEC(SIMPLE-consistent)。对于简单的问题收敛非常快速,不对压力进行修正,所以压力松弛因子可以设置为1 (3)Pressure-Implicit with Splitting of Operators (PISO)。对非定常流动问题或者包含比平均网格倾斜度更高的网格适用 (4)Fractional Step Method (FSM)对非定常流的分步方法。用于NITA格式,与PISO具有相同的特性。 二、对流插值(动量方程) FLUENT有五种方法:一阶迎风格式、幂率格式、二阶迎风格式、MUSL三阶格式、QUICK格式 (1)FLUENT默认采用一阶格式。容易收敛,但精度较差,主要用于初值计算。 (2)Power Lar.幂率格式,当雷诺数低于5时,计算精度比一阶格式要高。 (3)二阶迎风格式。二阶迎风格式相对于一阶格式来说,使用更小的截断误差,适用于三角形、四面体网格或流动与网格不在同一直线上;二阶格式收敛可能比较慢。 (4)MUSL(monotone upstream-centered schemes for conservation laws).当地3阶离散格式。主要用于非结构网格,在预测二次流,漩涡,力等时更精确。 (5)QUICK(Quadratic upwind interpolation)格式。此格式用于四边形/六面体时具有三阶精度,用于杂交网格或三角形/四面体时只具有二阶精度。 三、梯度插值梯度插值主要是针对扩散项。 FLUENT有三种梯度插值方案:green-gauss cell-based,Green-gauss node-based,least-quares cell based. (1)格林-高斯基于单元体。求解方法可能会出现伪扩散。 (2)格林-高斯基于节点。求解更精确,最小化伪扩散,推荐用于三角形网格上 (3)基于单元体的最小二乘法插值。推荐用于多面体网格,与基于节点的格林-高斯格式具有相同的精度和格式。 四、压力插值压力基分离求解器主要有五种压力插值算法。 (1)标准格式(Standard)。为FLUENT缺省格式,对大表妹边界层附近的曲线发现压力梯度流动求解精度会降低(但不能用于流动中压力急剧变化的地方——此时应该使用PRESTO!格式代替) (2)PRESTO!主要用于高旋流,压力急剧变化流(如多孔介质、风扇模型等),或剧烈弯曲的区域。 (3)Linear(线性格式)。当其他选项导致收敛困难或出现非物理解时使用此格式。 (4)second order(二阶格式)。用于可压缩流动,不能用于多孔介质、阶跃、风扇、VOF/MIXTURE多相流。 (5)Body Force Weighted体积力。当体积力很大时,如高雷诺数自然对流或高回旋流动中采用此格式。

tecplot执行fluent后处理--截面云图显示

Fluent后处理——Tecplot截面云图显示 本部分介绍一下用tecplot进行fluent后处理(等值线云图制作和矢量图制作)的一些小经验,希望能帮到和我一样在fluent数据处理途中遇到小问题的同学。如有问题,请多指正。 1、将导入的三维数据通过slices工具切出截面,并读取该截面数据。 首先,在导入fluent数据时,选择的数据类型必须是Fluent Data Loader,同时选中fluent的cas和data文件进行导入;并且保证文件目录全是英文目录,否则容易出错。 对于已导入的三维cas和data文件,选中左侧的slices,并点击slices后面的details。 弹出slice Details对话框,在slice location后有下拉菜单,可以自定义切片位置,本文选择了Z-Planes,并将Show primary planes 下方的值设为0,代表切片为Z=0截面,如图中黄色截面表示切片位置。点击close。 点击如下图所示的data—extract—extract current slices。

弹出对话框,点击extract。 切片就会变为实体,如下图所示。

为了只显示切片,点击左侧栏中的zone style,在弹出的对话框中,将slice:Z=0上面的内容全部取消选择,如下图。 选中左侧栏中的contour,云图就会显示出来(vector同理)。视角可以自由调整。至此,读取切片数据完成。

2、如何不让云图像蒙了一层雾一般,使其颜色更鲜艳?如下图,云图颜色不鲜亮。 取消下图中lighting选项,云图颜色恢复明亮。

FLUENT操作过程及全参数选择

振动流化床仿真操作过程及参数选择 1创建流化床模型。 根据靳海波论文提供的试验机参数,创建流化床模型。流化床直148mm 高1m开孔率9%孔径2mm在筛板上铺两层帆布保证气流均布。 因为实验机为一个圆形的流化床,所以可简化为仅二维模型。而实际实验中流化高度远小于1m甚至500mm所以为提高计算时间,可将模型高度缩为500mm由于筛板上铺设两层帆布以达到气流均分的目的,所以认为沿整个筛板的进口风速为均匀的。最终简化模型如下图所示: 上图为流化后的流化床模型,可以看出流化床下端的网格相对上端较密,因为流化行为主要发生的流化床下端,为了加快计算时间,所以采用这种下密上疏的划分方式。其中进口设置为velocity inlet ;出口设置为outflow ;左右两边分为设置为wall。在GAMBIT中设置完毕后,输出二维模型vfb.msh。 outflow 边界条件不需要给定任何入口的物理条件,但是应用也会有限制,大致为以下四点: 1.只能用于不可压缩流动

2.出口处流动充分发展 3.不能与任何压力边界条件搭配使用(压力入口、压力出口) 4.不能用于计算流量分配问题(比如有多个出口的问题) 2 打开FLUENT 6326,导入模型vfb.msh 点击GRID—CHECK检查网格信息及模型中设置的信息,核对是否正确,尤其查看是否出现负体积和负面积,如出现马上修改。核对完毕后,点击GRID-SCAL弹出SCALEGRID窗口,设置单位为mm 并点击change length unit 按钮。具体设置如下: 3设置求解器 保持其他设置为默认,更改TIME为unsteady,因为实际流化的过程是随时间变化的。 (1)pressure based 求解方法在求解不可压流体时,如果我们联立求解 从动量方程和连续性方程离散得到的代数方程组,可以直接得到各速

基于FLUENT的离心风机性能优化

陈一晓?龚一艳?陈小兵?等.基于FLUENT的离心风机性能优化[J].江苏农业科学?2019?47(16):250-254. doi:10.15889/j.issn.1002-1302.2019.16.055 基于FLUENT的离心风机性能优化 陈一晓?龚一艳?陈小兵?张一晓?王一果?缪友谊?刘德江 (农业部南京农业机械化研究所?江苏南京210014) 一一摘要:机械化高效施药是目前对喷药机械的基本要求?风送植保机械被普遍使用?风机是其中的关键部件?因此风 机的性能直接决定了植保机械的性能?采用正交试验方法对影响风机的主要参数进行优化?在FLUENT中对各个正交试验方案进行分析得到?叶片数对风机性能的影响最大?叶片出口安装角对风机性能的影响最小?优化后的风机方案与原风机相比?风机流量增加了3.92%?效率提高了17.07%?采用FLUENT对离心风机进行性能分析可以为优化离心风机性能提供理论依据? 一一关键词:背负式喷雾喷粉机?离心风机?叶轮?FLUENT?流场分析?正交试验?优化性能?最优方案一一中图分类号:S49一一文献标志码:A一一文章编号:1002-1302(2019)16-0250-05收稿日期:2018-03-29 基金项目:国家重点研发计划(编号:2017YFD0200303)?江苏省重点研发计划(现代农业)(编号:BE2016303)?现代农业产业技术体系建设专项-西甜瓜产业技术体系(编号:CARS-25)? 作者简介:陈一晓(1989 )?女?山东德州人?硕士?研究实习员?主要从事植保施药技术与装备研究?E-mail:chenxiao6105@163.com?通信作者:龚一艳?硕士?研究员?主要从事植保施药技术与装备研究?Tel:(025)84346241?E-mail:nnnGongyan@qq.com? 一一背负式喷雾喷粉机是一种典型的小型植保机械?因其轻 便二灵活二效率高的特点已被广泛应用于水稻二棉花二玉米二小麦二果树等大面积农作物病虫害防治?背负式喷雾喷粉机的射程二雾化效果二喷量等关键指标主要取决于离心风机的性能?因此离心风机是背负式喷雾喷粉机的关键部件?它的功 用主要是产生高速气流?将药液破碎雾化或将药粉吹散?并将之送向远方[1]? 背负式喷雾喷粉机上所使用的离心风机均为小型高速离 心风机?风机是风送植保机械的核心部件?良好的风机性能能够提高雾滴喷洒的均匀性?提高沉积量?降低飘移量[2]?虽然目前存在的风送式植保机械种类众多?但是用于风送式植保机械的风机并没有统一标准?风送式植保机械具有射程远二雾化均匀二穿透性好二靶标性好二雾滴飘移少等特点?因此被广泛应用于大田二果园等农药喷洒中[3]?配合不同的地理 位置以及作物本身?风送植保机械存在不同的种类?用于背负式喷雾喷粉机上的小型离心风机由于转速高等特点?离心风机各个参数对其性能的影响并没有明确的理论依据?因此针对不同的小型离心风机进行流场分析以得到各个参数对风机性能影响的研究很有必要?1一离心风机流场数值模拟 所有的流动都必须满足三大物理定律?即质量守恒定律二动力守恒定律以及能量守恒定律?相对应地就可以得到对应的质量守恒方程二动量守恒方程以及能量守恒方程?由于在离心风机中不需要考虑传热问题?因此能量守恒方程不需要考虑在其中? FLUENT中提供的湍流模型种类很多?但是目前还没有 适用于各种流动的湍流模型?因此要根据实际解决的问题及其对精确度的要求选择合适的湍流模型?考虑到风机实际的工作情况?本研究的湍流模型选择为K-ε模型?K-ε模型又分为标准K-ε模型二重整化(RNG)K-ε模型以及可实现K-ε模型?综合考虑风机运动的实际情况?最终选择K-ε模型中的可实现K-ε模型为本研究所用的湍流模型? FLUENT中提供了多种壁面函数处理方法?例如标准壁 面函数法二非平衡壁面函数法以及增强壁面处理?标准壁面函数法利用对数校正法提供了所必需的壁面边界条件?考虑到离心风机内部的结构比较复杂?本研究选择标准壁面函数法作为分析方法? 在流体流动中建立的基本方程为偏微分方程?在理论上可以求得其解?但是由于问题本身的复杂性?并不易得到它们的解析解或者近似解析解?因此在FLUENT中出现了离散化的概念?离散化就是将无限空间中的有效个体映射到有限的空间中?离散化的目的是将连续的偏微分方程组及其定解条件按照特定的规则在计算区域的离散网格上转换为代数方程?以得到连续系统的离散数值逼近解?在FLUENT中可以将控制方程的离散方法分为有限差分法二有限元法和有限体积法?本研究选择的是有限体积法? SIMPLE算法在1972年被提出并得到广泛的应用?是计 算不可压流场的主要方法?是后来对其算法进行改进与发展的基础?SIMPLEC算法的基本思想与SIMPLE算法一致?但是对通量的修正方法进行了改进?加快了收敛速度?本研究中求解算法采用的是SIMPLEC算法? 在FLUENT中通常认为残差小于10-3时为收敛?因此在对离心风机的内部流场进行分析时?对变量的监控指标设定为10-3? 2一离心风机内部流场模拟结果分析2.1一离心风机模型的建立 本研究的离心风机主要被应用在背负式喷雾喷粉机上? 因此采用前向叶轮和前弯式叶片?在Pro/ENGINEER软件 052 江苏农业科学一2019年第47卷第16期

相关文档
最新文档