文档库 最新最全的文档下载
当前位置:文档库 › fluent被动运动6DOF总结

fluent被动运动6DOF总结

fluent被动运动6DOF总结
fluent被动运动6DOF总结

利用CFD软件解决动网格问题,通常可分为以下两类:

(1)主动型动网格

主动型动网格问题通常指的是边界运动规律及运动状态已知,通常可由软件使用者通过函数或程序进行描述。在程序计算过程中,求解器调用边界运动轨迹描述程序实现边界运动。这类动网格例子很多,如各类泵、风扇等。

(2)被动型动网格

还有一类动网格问题,其边界运动规律往往是未知的,常常需要通过计算边界上的力或力矩,以此来求取边界的运动。在这类动网格计算设置中,网格变化规律难以预料,导致网格参数经常需要进行多次调整才能达到目的。这类例子在现实中其实也很多,比如风力发电机的叶轮、水轮机等。

解决主动型动网格问题比较容易,利用CFD软件提供的动网格模拟能力很容易解决。需要关注的地方是边界运动后,网格节点如何重新布置和生成。如在FLUENT软件中,其动网格主要包括三种网格功能:弹簧光顺、动态层及网格重构。利用网格重构功能几乎可以解决所有主动型动网格问题。

那被动型动网格问题怎么处理呢?一般来说,这类边界的运动都是由于内部流体对其压力所造成的,那么就涉及到力和力矩计算的问题。对于这类问题,在FLUENT软件中可以采用6DOF 模型进行计算。

需要注意的是,以上所有类型动网格计算均建立在边界为刚性的情况下。即不会计算由于流动产生的力的作用导致的边界变形。若要计算边界变形,则需要采用流固耦合方法,利用固体求解器计算。被动型动网格中的力和力矩均是压力对面的积分计算而来。

1、6DOF UDF宏

在FLUENT中利用6DOF是需要定义UDF宏的。该宏的定义形式如下:

DEFINE_SDOF_PROPERTIES(name, properties, dt ,time ,dtime)

函数中:

Name:宏名称

Real *properties:存储6DOF属性的数组

Dynamic_Thread *dt:存储制定的动网格属性

Real time:当前时间

Real dtime:时间步长

该UDF宏没有返回值。用户需要定义的变变量为name、properties、dt、dtime。

在利用该UDF宏的过程中,需要注意的是properties数组,其包含多种属性:

SDOF_MASS //质量,定义方式如properties[SDOF_MASS]=7.0定义质量7kg

SDOF_IXX,SDOF_IYY,SOF_IZZ,//X、Y、Z方向惯性矩

SDOF_IXY,SDOF_IXZ,SDOF_IYZ,//惯性积

SDOF_LOAD_LOCAL,//布尔值,FALSE为全局坐标系,TRUE为体坐标系(局部坐标系),默认为FALSE

SDOF_LOAD_F_X,SDOF_LOAD_F_Y,SDOF_LOAD_F_Z,//X,Y,Z方向外力

SDOF_LOAD_M_X,SDOF_LOAD_M_Y,SDOF_LOAD_M_Z,//分别为X,Y,Z方向外力矩

SDOF_ZERO_TRANS_X,SDOF_ZERO_TRANS_Y,SDOF_ZERO_TRANS_Z,//布尔值,TRUE表示该方向位移为0,默认为FALSE

SDOF_ZERO_ROT_X,SDOF_ZERO_ROT_Y,SDOF_ZERO_ROT_Z,//布尔值,TRUE表示该方向旋转自由度为0,默认值为FALSE

SDOF_SYMMETRY_X,SDOF_SYMMETRY_Y,SDOF_SYMMETRY_Z,//指定半模型的对称轴法向向量。在半模型运动时使用。

上面的变换变量可以应付绝大多数问题,然而在FLUENT中,用户还可以自定义变换变量,这也是通过6DOF属性变量实现,这些变量包括:

SDOF_CUSTOM_TRANS,//布尔值,TRUE表示使用自定义变换变量,FALUSE表示不使用SDOF_CTRANS_11,SDOF_CTRANS_12,SDOF_CTRANS_13,SDOF_CTRANS_21,SDOF_CTRANS_22,SDOF_CTRANS_23,SDOF_CTRANS_31,SDOF_CTRANS_32,SDOF_CTRANS_33,//坐标旋转变量SDOF_DTRANS_11,SDOF_DTRANS_12,SDOF_DTRANS_13,SDOF_DTRANS_21,SDOF_DTRANS_22,SDOF_DTRANS_23,SDOF_DTRANS_31,SDOF_DTRANS_32,SDOF_DTRANS_33,//衍生旋转矩阵这些旋转矩阵在实际使用中应用较少。

下面是一个简单的SDOF实例:

#include “udf.h”

DEFINE_SDOF_PROPERTIES(stage,prop,time,dtime)

{

prop[SDOF_MASS] = 800.0;

prop[SDOF_IXX] = 200;

Prop[SDOF_IYY] = 100;

Prop[SDOF_IZZ] = 100;

}

上面的UDF定义了一个名为stage的SDOF宏,且质量800kg,X方向转动惯量为200,Y方向转动惯量100,Z方向转动惯量100

2、刚体参数的获取

应用6DOF无可避免的会碰到刚体属性值的获取问题。对于简单的几何体,应用力学公式可以很方便的计算出诸如质量、转动惯量等。但是对于复杂的几何模型,要想通过数学公式计算出这些物理量几乎是不可能的。

其实,几乎所有的CAD建模软件中均具有这些物理量的计算能力,不过要注意区分计算的是全局坐标系还是局部坐标系。

6DOF实例1:设置过程

对于被动型动网格计算问题,可以采用6DOF动网格模型解决。利用该模型时,需要确定计算模型中运动部件的质量、三方向转动惯量及惯性矩、重心坐标。因此在几何模型创建过程中,尽量使这些值容易获取(比如将重心位置放置在坐标原点),当然这不是必须的,在CAD 建模软件中都包含了这些物理量的查询。

如图1所示的几何体,在受到图中流动方向的流体作用下,若轴向旋转自由度未被约束,则会产生旋转位移。本次实例几何如图1所示。

在建立流体计算域之前,需要利用CAD软件获取几何的特征物理量,我们这里使用Solidworks实现此功能。

1、获取几何属性

利用solidworks载入几何体,进入【评估】>【质量属性】,选择几何体进行计算,如图2所示。

图 2 质量属性

在计算几何质量属性对话框中,可以点击【选项】按钮,设置材料参数,如图3所示。

图 3 几何参数

从图2可以看出,该几何体的重心(0,0,25.19),Z方向惯性矩IZZ=324047.793g*mm^2,由于我们的几何体约束了X,Y,Z三方向移动自由度,X,Y方向的旋转自由度,因此只需要Z 方向的惯性矩。需要注意的是,这里的单位是mm,而UDF中必须全部为国际单位制,需要做单位转换。

2、计算域模型

在workbench中的DM中创建计算域模型。主要是利用DM中方便的布尔运算功能。建立的计算域如图4所示。

图 4 计算域模型

为计算域边界命名,分别命名内部区域面、入口面、出口面、圆柱面边界为:WallRotation、VelocityInlet、Outlet、wallCylinder。

计算域建立完毕后即进行网格的划分,这里不详细描述。

3、UDF宏的编写

利用文本编辑器编写UDF宏文件。如下:

#include "udf.h"

DEFINE_SDOF_PROPERTIES(rotate,prop,dt,time,dtime)

{

prop[SDOF_MASS] = 845.154E-3; /*质量*/

prop[SDOF_IXX]= 336179E-9; /*x方向惯性矩*/

prop[SDOF_IYY] = 336179E-9; /*y方向惯性矩*/

prop[SDOF_IZZ]= 324047.793E-9; /*z方向惯性矩*/

prop[SDOF_ZERO_TRANS_X]=TRUE; /*限制x方向位移*/

prop[SDOF_ZERO_TRANS_Y]=TRUE; /*限制y方向位移*/

prop[SDOF_ZERO_TRANS_Z]=TRUE; /*限制z方向位移*/

prop[SDOF_ZERO_ROT_X]=TRUE; /*限制x方向旋转*/

prop[SDOF_ZERO_ROT_Y]=TRUE;/*限制y方向旋转*/

printf("\nstage: updated 6dof properties");

}

4、进入FLUENT中设置

在FLUENT中主要包括动网格参数的设置。

其他设置:瞬态计算、标准K-E湍流模型、标准壁面函数、工作介质为water、入口5m/s,出口静压为0。

动网格设置:

(1)、激活sixDOF

如图5所示,进入Dynamic Mesh面板,激活Dynamic Mesh,勾选Remeshing选项及SixDOF 项。

图 5

(2)、动网格参数设置

进入mesh methods下的settings按钮,进入动网格参数设置。图6、图7分别为smoothing 与Remeshing的设置。

图6 smoothing面板

图7 remeshing面板

6DOF中的settings设置对话框如图8所示,这里我们不考虑重力影响,设置三方向重力加速度为0。同时勾选write motion history项,记录保存运动中几何姿态。

图 8 6DOF设置

(3)、编译加载UDF宏

利用菜单【Define】>【user defined】>【function】>【Complied】编译宏。

图 9 编译及加载宏文件

(4)、设置动网格区域

动网格设置面板中,点击【Create/Edit…】按钮定义动网格区域。需要设置的区域如图10所示的红色框选部分。

图10运动区域设置

同时注意设置meshing Options标签页下的网格高度。

5、需要注意的问题

3D动网格问题,首先要注意的是负体积问题。前面说过,由于运动状态由受力状态控制,在计算之前是未知的,所有调整动网格参数比较困难。不过仔细耐心的修正还是可以解决的。其次,关于这类被动运动问题,一些模型较为简单的几何体,实际上是可以化被动为主动,只是计算壁面受力,运动状态通过受力状态用UDF确定。比如此例,我们可以在计算过程中通过宏计算出每次迭代加载在重心上的力和力矩,然后通过力学公式确定几何体的运动状态。这样被动运动问题就变为人为定义运动状态的主动动网格问题。

这次几何有点复杂,计算时间比较长,不过主要用于演示6DOF模型的定义。下次使用简单的模型进行演示。

第三章,湍流模型

第三章,湍流模型 第一节, 前言 湍流流动模型很多,但大致可以归纳为以下三类: 第一类是湍流输运系数模型,是Boussinesq 于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即: 2 1 21 x u u u t ??=-μρ 3-1 推广到三维问题,若用笛卡儿张量表示,即有: ij i j j i t j i k x u x u u u δρμρ32 -??? ? ????+ ??=''- 3-2 模型的任务就是给出计算湍流粘性系数t μ的方法。根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。 第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。 第三类是大涡模拟。前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes 方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。 实际求解中,选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。 FLUENT 提供的湍流模型包括:单方程(Spalart-Allmaras )模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。 湍流模型种类示意图 第二节,平均量输运方程 包含更多 物理机理 每次迭代 计算量增加 提的模型选 RANS-based models

雷诺平均就是把Navier-Stokes 方程中的瞬时变量分解成平均量和脉动量两部分。对于速度,有: i i i u u u '+= 3-3 其中,i u 和i u '分别是平均速度和脉动速度(i=1,2,3) 类似地,对于压力等其它标量,我们也有: φφφ'+= 3-4 其中,φ表示标量,如压力、能量、组分浓度等。 把上面的表达式代入瞬时的连续与动量方程,并取平均(去掉平均速度i u 上的横线),我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式: 0)(=?? +??i i u x t ρρ 3-5 () j i j l l ij i j j i j i i u u x x u x u x u x x p Dt Du -?? +???????????? ????-??+????+??-=ρδμρ32 3-6 上面两个方程称为雷诺平均的Navier-Stokes (RANS )方程。他们和瞬时Navier-Stokes 方程有相同的形式,只是速度或其它求解变量变成了时间平均量。额外多出来的项j i u u ''-ρ是雷诺应力,表示湍流的影响。如果要求解该方程,必须模拟该项以封闭方程。 如果密度是变化的流动过程如燃烧问题,我们可以用法夫雷(Favre )平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外,所有变量都用密度加权平均。变量的密度加权平均定义为: ρρ/~ Φ=Φ 3-7 符号~表示密度加权平均;对应于密度加权平均值的脉动值用Φ''表示,即有: Φ''+Φ=Φ~ 。很显然,这种脉动值的简单平均值不为零,但它的密度加权平均值等于零,即: 0≠Φ'', 0=Φ''ρ Boussinesq 近似与雷诺应力输运模型 为了封闭方程,必须对额外项雷诺应力j i u u -ρ进行模拟。一个通常的方法是应用Boussinesq 假设,认为雷诺应力与平均速度梯度成正比,即: ij i i t i j j i t j i x u k x u x u u u δμρμρ)(32 ??+-??? ? ????+??=''- 3-8 Boussinesq 假设被用于Spalart-Allmaras 单方程模型和ε-k 双方程模型。Boussinesq 近似 的好处是与求解湍流粘性系数有关的计算时间比较少,例如在Spalart-Allmaras 单方程模型中,只多求解一个表示湍流粘性的输运方程;在ε-k 双方程模型中,只需多求解湍动能k 和耗散率ε两个方程,湍流粘性系数用湍动能k 和耗散率ε的函数。Boussinesq 假设的缺点是认为湍流粘性系数t μ是各向同性标量,对一些复杂流动该条件并不是严格成立,所以具有其应用限制性。

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中文全教程1-250

FLUENT 教程 赵玉新 I、目录 第一章、开始 第二章、操作界面 第三章、文件的读写 第四章、单位系统 第五章、读入和操作网格 第六章、边界条件 第七章、物理特性 第八章、基本物理模型 第九章、湍流模型 第十章、辐射模型 第十一章、化学输运与反应流 第十二章、污染形成模型 第十三章、相变模拟 第十四章、多相流模型 第十五章、动坐标系下的流动 第十六章、解算器的使用 第十七章、网格适应 第十八章、数据显示与报告界面的产生 第十九章、图形与可视化 第二十章、Alphanumeric Reporting 第二十一章、流场函数定义 第二十二章、并行处理 第二十三章、自定义函数 第二十四章、参考向导 第二十五章、索引(Bibliography) 第二十六章、命令索引 II、如何使用该教程 概述 本教程主要介绍了FLUENT 的使用,其中附带了相关的算例,从而能够使每一位使用 者在学习的同时积累相关的经验。本教程大致分以下四个部分:第一部分包括介绍信息、用户界面信息、文件输入输出、单位系统、网格、边界条件以及物理特性。第二和第三部分包含物理模型,解以及网格适应的信息。第四部分包括界面的生成、后处理、图形报告、并行处理、自定义函数以及FLUENT 所使用的流场函数与变量的定义。 下面是各章的简略概括 第一部分: z开始使用:本章描述了FLUENT 的计算能力以及它与其它程序的接口。介绍了如何对具体的应用选择适当的解形式,并且概述了问题解决的大致步骤。在本章中,我们给出

了一个可以在你自己计算机上运行的简单的算例。 z使用界面:本章描述了用户界面、文本界面以及在线帮助的使用方法。同时也提供了远程处理与批处理的一些方法。(请参考关于特定的文本界面命令的在线帮助) z读写文件:本章描述了FLUENT 可以读写的文件以及硬拷贝文件。 z单位系统:本章描述了如何使用FLUENT 所提供的标准与自定义单位系统。 z读和操纵网格:本章描述了各种各样的计算网格来源,并解释了如何获取关于网格的诊断信息,以及通过尺度化(scale)、分区(partition)等方法对网格的修改。本章还描述了非一致(nonconformal)网格的使用. z边界条件:本章描述了FLUENT 所提供的各种类型边界条件,如何使用它们,如何定义它们and how to define boundary profiles and volumetric sources. z物理特性:本章描述了如何定义流体的物理特性与方程。FLUENT 采用这些信息来处理你的输入信息。 第二部分: z基本物理模型:本章描述了FLUENT 计算流体流动和热传导所使用的物理模型(包括自然对流、周期流、热传导、swirling、旋转流、可压流、无粘流以及时间相关流)。以及在使用这些模型时你需要输入的数据,本章也包含了自定义标量的信息。 z湍流模型:本章描述了FLUENT 的湍流模型以及使用条件。 z辐射模型:本章描述了FLUENT 的热辐射模型以及使用条件。 z化学组分输运和反应流:本章描述了化学组分输运和反应流的模型及其使用方法。本章详细的叙述了prePDF 的使用方法。 z污染形成模型:本章描述了NOx 和烟尘的形成的模型,以及这些模型的使用方法。 第三部分: z相变模拟:本章描述了FLUENT 的相变模型及其使用方法。 z离散相变模型:本章描述了FLUENT 的离散相变模型及其使用方法。 z多相流模型:本章描述了FLUENT 的多相流模型及其使用方法。 z Flows in Moving Zones(移动坐标系下的流动):本章描述了FLUENT 中单一旋转坐标系,多重移动坐标系,以及滑动网格的使用方法。 z Solver 的使用:本章描述了如何使用FLUENT 的解法器(solver)。 z网格适应:本章描述了explains the solution-adaptive mesh refinement feature in FLUENT and how to use it 第四部分: z显示和报告数据界面的创建:本章描述了explains how to create surfaces in the domain on which you can examine FLUENT solution data z图形和可视化:本章描述了检验FLUENT 解的图形工具 z Alphanumeric Reporting:本章描述了如何获取流动、力、表面积分以及其它解的数据。 z流场函数的定义:本章描述了如何定义FLUENT 面板内出现的变量选择下拉菜单中的流动变量,并且告诉我们如何创建自己的自定义流场函数。 z并行处理:本章描述了FLUENT 的并行处理特点以及使用方法 z自定义函数:本章描述了如何通过用户定义边界条件,物理性质函数来形成自己的FLUENT 软件。 如何使用该手册 z根据你对CFD 以及FLUENT 公司的熟悉,你可以通过各种途径使用该手册 对于初学者,建议如下:

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经验总结

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可能会得到更好的结果,尤其是可以应用增加的亚松驰迭代时,具体介绍如下: 对于相对简单的问题(如:没有附加模型激活的层流流动),其收敛性已经被压力速

fluent使用基本步骤

fluent使用基本步骤 步骤一:网格 1.读入网格(*.msh) File →Read →Case 读入网格后,在窗口显示进程 2.检查网格 Grid →Check Fluent对网格进行多种检查,并显示结果。注意最小容积,确保最小容积值为正。 3.显示网格 Display →Grid ①以默认格式显示网格 能够用鼠标右键检查边界区域、数量、名称、类型将在窗口显示,本操作关 于同样类型的多个区域情形专门有用,以便快速区别它们。 4.网格显示操作 Display →Views (a)在Mirror Planes面板下,axis (b)点击Apply,将显示整个网格 (c)点击Auto scale, 自动调整比例,并放在视窗中间 (d)点击Camera,调整目标物体位置 (e)用鼠标左键拖动指标钟,使目标位置为正 (f)点击Apply,并关闭Camera Parameters 和Views窗口 步骤二:模型 1. 定义瞬时、轴对称模型

Define →models→Solver (a)保留默认的,Segregated解法设置,该项设置,在多相运算时使用。 (b)在Space面板下,选择Axisymmetric (c)在Time面板下,选择Unsteady 2. 采纳欧拉多相模型 Define→Models→Multiphase (a) 选择Eulerian作为模型 (b)假如两相速度差较大,则需解滑移速度方程 (c)假如Body force比粘性力和对流力大得多,则需选择implicit body force 通过考虑压力梯度和体力,加快收敛 (d)保留设置不变

FLUENT学习经验总结(狠珍贵,学长传授)

1对于刚接触到FLUENT新手来说,面对铺天盖地的学习资料和令人难读的FLUENT help,如何学习才能在最短的时间内入门并掌握基本学习方法呢? 答:学习任何一个软件,对于每一个人来说,都存在入门的时期。认真勤学是必须的,什么是最好的学习方法,我也不能妄加定论,在此,我愿意将我三年前入门FLUENT心得介绍一下,希望能给学习FLUENT的新手一点帮助。 由于当时我需要学习FLUENT来做毕业设计,老师给了我一本书,韩占忠的《FLUENT流体工程仿真计算实例与应用》,当然,学这本书之前必须要有两个条件,第一,具有流体力学的基础,第二,有FLUENT 安装软件可以应用。然后就照着书上二维的计算例子,一个例子,一个步骤地去学习,然后学习三维,再针对具体你所遇到的项目进行针对性的计算。不能急于求成,从前处理器GAMBIT,到通过FLUENT进行仿真,再到后处理,如TECPLOT,进行循序渐进的学习,坚持,效果是非常显著的。如果身边有懂得FLUENT的老师,那么遇到问题向老师请教是最有效的方法,碰到不懂的问题也可以上网或者查找相关书籍来得到答案。另外我还有本《计算流体动力学分析》王福军的,两者结合起来学习效果更好。 2 CFD计算中涉及到的流体及流动的基本概念和术语:理想流体和粘性流体;牛顿流体和非牛顿流体;可压缩流体和不可压缩流体;层流和湍流;定常流动和非定常流动;亚音速与超音速流动;热传导和扩散等。A.理想流体(Ideal Fluid)和粘性流体(Viscous Fluid): 流体在静止时虽不能承受切应力,但在运动时,对相邻的两层流体间的相对运动,即相对滑动速度却是有抵抗的,这种抵抗力称为粘性应力。流体所具备的这种抵抗两层流体相对滑动速度,或普遍说来抵抗变形的性质称为粘性。粘性的大小依赖于流体的性质,并显著地随温度变化。实验表明,粘性应力的大小与粘性及相对速度成正比。当流体的粘性较小(实际上最重要的流体如空气、水等的粘性都是很小的),运动的相对速度也不大时,所产生的粘性应力比起其他类型的力如惯性力可忽略不计。此时我们可以近似地把流体看成无粘性的,这样的流体称为理想流体。十分明显,理想流体对于切向变形没有任何抗拒能力。这样对于粘性而言,我们可以将流体分为理想流体和粘性流体两大类。应该强调指出,真正的理想流体在客观实际中是不存在的,它只是实际流体在某些条件下的一种近似模型。 B.牛顿流体(Newtonian Fluid)和非牛顿流体(non-Newtonian Fluid): 日常生活和工程实践中最常遇到的流体其切应力与剪切变形速率符合下式的线性关系,称为牛顿流体。而切应力与变形速率不成线性关系者称为非牛顿流体。图2-1(a)中绘出了切应力与变形速率的关系曲线。其中符合上式的线性关系者为牛顿流体。其他为非牛顿流体,非牛顿流体中又因其切应力与变形速率关系特点分为膨胀性流体(Dilalant),拟塑性流体(Pseudoplastic),具有屈服应力的理想宾厄流体(Ideal Bingham Fluid)和塑性流体(Plastic Fluid)等。通常油脂、油漆、牛奶、牙膏、血液、泥浆等均为非牛顿流体。非牛顿流体的研究在化纤、塑料、石油、化工、食品及很多轻工业中有着广泛的应用。图2-1(b)还显示出对于有些非牛顿流体,其粘滞特性具有时间效应,即剪切应力不仅与变形速率有关而且与作用时间有关。当变形速率保持常量,切应力随时间增大,这种非牛顿流体称为震凝性流体(Rheopectic Fluid)。当变形速率保持常量而切应力随时间减小的非牛顿流体则称为触变性流体(Thixotropic Fluid)。 C.可压缩流体(Compressible Fluid)和不可压缩流体(Incompressible Fluid): 在流体的运动过程中,由于压力、温度等因素的改变,流体质点的体积(或密度,因质点的质量一定),或多或少有所改变。流体质点的体积或密度在受到一定压力差或温度差的条件下可以改变的这个性质称为压缩性。真实流体都是可以压缩的。它的压缩程度依赖于流体的性质及外界的条件。例如水在100个大气压下,容积缩小0.5%,温度从20°变化到100°,容积降低4%。因此在一般情况下液体可以近似地看成不可压的。但是在某些特殊问题中,例如水中爆炸或水击等问题,则必须把液体看作是可压缩的。气体的压缩性比液体大得多,所以在一般情形下应该当作可压缩流体处理。但是如果压力差较小,运动速度较小,并且没有很大的温度差,则实际上气体所产生的体积变化也不大。此时,也可以近似地将气体视为不可压缩的。 在可压缩流体的连续方程中含密度,因而可把密度视为连续方程中的独立变量进行求解,再根据气体的状态方程求出压力。不可压流体的压力场是通过连续方程间接规定的。由于没有直接求解压力的方程,不可压流体的流动方程的求解具有其特殊的困难。 D. 层流(Laminar Flow)和湍流(Turbulent Flow):

FLUENT中常用的湍流模型

The Spalart-Allmaras模型 对于解决动力漩涡粘性,Spalart-Allmaras 模型是相对简单的方程。它包含了一组新的方程,在这些方程里不必要去计算和剪应力层厚度相关的长度尺度。Spalart-Allmaras 模型是设计用于航空领域的,主要是墙壁束缚流动,而且已经显示出很好的效果。在透平机械中的应用也愈加广泛。 在原始形式中Spalart-Allmaras 模型对于低雷诺数模型是十分有效的,要求边界层中粘性影响的区域被适当的解决。在FLUENT中,Spalart-Allmaras 模型用在网格划分的不是很好时。这将是最好的选择,当精确的计算在湍流中并不是十分需要时。再有,在模型中近壁的变量梯度比在k-e模型和k-ω模型中的要小的多。这也许可以使模型对于数值的误差变得不敏感。想知道数值误差的具体情况请看5.1.2。 需要注意的是Spalart-Allmaras 模型是一种新出现的模型,现在不能断定它适用于所有的复杂的工程流体。例如,不能依靠它去预测均匀衰退,各向同性湍流。还有要注意的是,单方程的模型经常因为对长度的不敏感而受到批评,例如当流动墙壁束缚变为自由剪切流。 标准k-e模型 最简单的完整湍流模型是两个方程的模型,要解两个变量,速度和长度尺度。在FLUENT中,标准k-e模型自从被Launder and Spalding提出之后,就变成工程流场计算中主要的工具了。适用范围广、经济,有合理的精度,这就是为什么它在工业流场和热交换模拟中有如此广泛的应用了。它是个半经验的公式,是从实验现象中总结出来的。 由于人们已经知道了k-e模型适用的范围,因此人们对它加以改造,出现了RNG k-e模型和带旋流修正k-e 模型。k-ε模型中的K和ε物理意义:k是紊流脉动动能(J),ε是紊流脉动动能的耗散率(%);k越大表明湍流脉动长度和时间尺度越大,ε越大意味着湍流脉动长度和时间尺度越小,它们是两个量制约着湍流脉动。 RNG k-e模型 RNG k-e模型来源于严格的统计技术。它和标准k-e模型很相似,但是有以下改进: ?RNG模型在e方程中加了一个条件,有效的改善了精度。 ?考虑到了湍流漩涡,提高了在这方面的精度。 ?RNG理论为湍流Prandtl数提供了一个解析公式,然而标准k-e模型使用的是用户提供的常数。 ?然而标准k-e模型是一种高雷诺数的模型,RNG理论提供了一个考虑低雷诺数流动粘性的解析公式。这些公式的效用依靠正确的对待近壁区域 这些特点使得RNG k-e模型比标准k-e模型在更广泛的流动中有更高的可信度和精度。 带旋流修正的k-e模型 带旋流修正的k-e模型是近期才出现的,比起标准k-e模型来有两个主要的不同点。 ?带旋流修正的k-e模型为湍流粘性增加了一个公式。 ?为耗散率增加了新的传输方程,这个方程来源于一个为层流速度波动而作的精确方程。 术语“realizable”,意味着模型要确保在雷诺压力中要有数学约束,湍流的连续性。带旋流修正的k-e模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。带旋流修正的k-e模型和RNG k-e模型都显现出比标准k-e模型在强流线弯曲、漩涡和旋转有更好的表现。由于带旋流修正的k-e模型是新出现的模型,所以现在还没有确凿的证据表明它比RNG k-e模型有更好的表现。但是最初的研究表明带旋流修正的k-e模型在所有k-e模型中流动分离和复杂二次流有很好的作用。带旋流修正的k-e模型的一个不足是在主要计算旋转和静态流动区域时不能提供自然的湍流粘度。这是因为带旋流修正的k-e模型在定义湍流粘度时考虑了平均旋度的影响。这种额外的旋转影响已经在单一旋转参考系中得到证实,而且表现要好于标准k-e模型。由于这些修改,把它应用于多重参考系统中需要注意。 标准k-ω模型 标准k-ω模型是基于Wilcox k-ω模型,它是为考虑低雷诺数、可压缩性和剪切流传播而修改的。Wilcox k-ω模型预测了自由剪切流传播速率,像尾流、混合流动、平板绕流、圆柱绕流和放射状喷射,因而可以应用于墙壁束缚流动和自由剪切流动。标准k-e模型的一个变形是SST k-ω模型,它在FLUENT中也是可用的,将在10.2.9中介绍它。 剪切压力传输(SST)k-ω模型

fluent 计算错误汇总

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: 改成:x: 2d -driver msw 如果还有三维的,可以再建立一个快捷方式改成: x: 3d -driver msw 这就可以直接调用了。如果不是以上原因引起的话,也有可能是和别的软件冲突,如MATLAB等,这也会使fluent无法显示图像。 Q1:GAMBIT安装后无法运行,出错信息是“unable find Exceed X Server” A. GAMBIT需要装EXCEED才能用。 gambit的运行:先运行命令提示符,输入gambit,回车 fluent的运行:直接在开始-程序-Fluent Inc里面 Q2:Fluent安装后无法运行,出错信息是甥?挱湵扡敬映湩層漯数? A. FLUENT和GAMBIT需要把相应文件拷贝到license目录下 文件?gambit时提示找不到gambit出错信息:运行Q3: 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,重起检查 Q5:Gambit运行失败,出错信息“IDENTIFIER default_ Server ” 等文件default_id.*的缺省文件已经打开,到用户默认目录删除gambitA.

Fluent 湍流模型小结

Fluent 湍流模型小结湍流模型目前计算流体力学常用的湍流的数值模拟方法主要有以下三种: 直接模拟(direct numerical&Oσλαση; simulation, DNS) 直接数值模拟(DNS)特点在湍流尺度下的网格尺寸内不引入任何封闭模型的前提下对Navier-Stokes方程直接求解。这种方法能对湍流流动中最小尺度涡进行求解,要对高度复杂的湍流运动进行直接的数值计算,必须采用很小的时间与空间步长,才能分辨出湍流中详细的空间结构及变化剧烈的时间特性。基于这个原因,DNS目前仅限于相对低的雷诺数中湍流流动模型。另外,利用DNS模型对湍流运动进行直接的数值模拟对计算工具有很高的要求,计算机的内存及计算速度要非常的高,目前DNS模型还无法应用于工程数值计算,还不能解决工程实际问题。 大涡模拟(large&Oσλαση; eddy simulation, LES) 大涡模拟(LES)是基于网格尺度封闭模型及对大尺度涡进行直接求解N-S方程,其网格尺度比湍流尺度大,可以模拟湍流发展过程的一些细节,但其计算量仍很大,也仅用于比较简单的剪切流运动及管流。大涡模拟的基础是:湍流的脉动与混合主要是由大尺度的涡造成的,大尺度涡是高度的非各向同性,而且随流动的情形而异。大尺度的涡通过相互作用把能量传递给小尺度的涡,而小尺度的涡旋主要起到耗散能量的作用,几乎是各向同性的。这些对涡旋的认识基础就导致了大涡模拟方法的产生。Les大涡模拟采用非稳态的N-S方程直接模拟大尺度涡,但不计算小尺度涡,小涡对大涡的影响通过近似的模拟来考虑,这种影响称为亚格子Reynolds应力模型。大多数亚格子Reynolds模型都是将湍流脉动所造成的影响用一个湍流粘性系数,既粘涡性来描述。LES对计算机的容量和CPU的要求虽然仍然很高,但是远远低于DNS方法对计算机的要求,因而近年来的研究与应用日趋广泛。 应用Reynolds时均方程(Reynolds-averaging&Oσλαση; equations)的模拟方法 许多流体力学的研究和数值模拟的结果表明,可用于工程上现实可行的湍流模拟方法仍然是基于求解Reynolds时均方程及关联量输运方程的湍流模拟方法,即湍流的统观模拟方法。统观模拟方法的基本思想是用低阶关联量和平均流性质来模拟未知的高阶关联项,从而封闭平均方程组或关联项方程组。虽然这种方法在湍流理论中是最简单的,但是对工程应用而言仍然是相当复杂的。即便如此,在处理工程上的问题时,统观模拟方法仍然是最有效、最经济而且合理的方法。在统观模型中,使用时间最长,积累经验最丰富的是混合长度模型和K-E模型。其中混合长度模型是最早期和最简单的湍流模型。该模型是建立在层流粘性和湍流粘性的类比、平均运动与湍流的脉动的概念上的。该模型的优点是简单直观、无须增加微分方程。缺点是在模型中忽略了湍流的对流与扩散,对于复杂湍流流动混合长度难以确定。 到目前为止,工程中应用最广泛的是k-ε模型。另外针对k-ε模型的不足之处,许多学者通过对K-E模型的修正和发展,开始采用雷诺应力模型(DSM)和代数应力模型(ASM)。近年来,DSM模型已用来预报燃烧室及炉内的强旋及浮力流动。很多情况下能够给出优于k-ε模型的结果。但是该模型也有不足之处,首先它对工程预报来说太复杂,其次经验系数太多难以确定,此外,对压力应变项的模拟还有争议。更主要的是,尽管这一模型考虑了各种应变效应,但是其总精度并不总是高于其它模型,这些缺点导致了DSM模型没有得到广泛的应用。总之,虽然从本质上讲DSM模型和ASM模型比k-ε模型对湍流流场的模拟更加合理,但DSM和ASM中仍然采用精度不高的E方程,模型中常数的通用性还没有得到广泛的验证,边界条件不好给定,计算也比较复杂。正因为如此,目前用计算解决湍流问题时仍然采用比较成熟的K-E模型。 需要注意的是: 1、大涡模拟有自己的亚格子封闭模型,这和k-ε模型完全是两回事。LES的亚格子模型表

fluent udf学习总结

fluent udf 阶段性小结 ——Flying_U 因工作需要,最近开始学习fluent二次开发功能。现在,根据工作日志将这一段时间主要的学习过程和总结的经验整理如下。 学习计划:从4月5号开始,计划花上一个月的时间了解和学习fluent udf的基本知识。目标是能够运用udf初步实现物理模型简化、掌握udf的基本用法并能根据工作需要实现相关udf功能。 4.5-4.6:浏览网站尽可能更多了解udf的知识,结合自己的实际情况分析那些知识是自己需要进一步深入学习的。 此阶段总结:1.udf是用户自定义函数的简称,其通过与fluent接口连接实现扩展fluent功能的作用。udf的主要功能有: ●定制边界条件、材料属性、表面和体积反应率、fluent输运方程的源项、用户自定义的 标量方程的源项、扩散函数等 ●调整每次迭代后的计算结果 ●初始化流场的解 ●在需要时进行udf的异步执行 ●强化后处理功能 ●强化现有的udf模型 ●传送返回值、修改fluent变量、操作外部文件案例和data文件 2. 自己现在想要实现的是udf功能是定制边界条件、定制fluent输运方程的源项、初始化流场的解和强化后处理功能;(刚开始自己也不太明确自己到底想用udf来做什么,对应上udf的主要功能是哪一部分,然后对自己不懂没理解的功能一一查询。) 3. 有相关资料的渠道有:百度知道,百度文库和doc88。其中,百度文库各种教程最多,百度知道能够快速定位回答具体的问题,doc88资料觉得更深入一些。(对搜集的资料进行及时的整理和归纳对自己学习有很大助力,很多资料都是不完全的或者自己当时没有完全理解的需要不同版本或者前后不同时间段对照着学习。) 4.6-4.9 根据自己的需求在udf帮组手册中查找实例并尝试按实例进行对照练习,初步了解udf相关知识,打通udf实现的过程(udf编写、编译和连接)。主要目的是了解udf的基本用法,初步了解udf宏命令。 此阶段总结:1.udf帮助手册里的实例对初学者特别有用,例子难度小,侧重流程和用法。 2. 解释型udf使用基本过程:将物理模型简化,确定udf实现的功能;编制c语言源代码;启动fluent并完成相关设置;加载并解释c语言代码;将udf与fluent相关联。初次尝试觉得难度很小,只是对c语言源程序的边界有点模糊,特别是在语言结构方面。 3. 进行编译型udf使用尝试,根据需要需先安装一个c编辑器,因手头有vs2012所以直接安装vs2012。安装vs2012的过程有点艰辛,因为自己电脑上装过其他的VS版本(未完全卸载),以及其他的一些开源c编辑器(一直在用),安装vs2012过程老是不成功(这里推荐使用vs2008,网上教程较多,出现相关的问题都能较容易地找到方法解决)。最后的解决办法是重装系统(各种尝试都未能解决安装问题后,可能与以前安装过的vs卸载不完全有关),并逐步尝试环境变量的设置。vs2012与fluent环境变量成功设置如下(win7 64位系统):

fluent湍流模型

第十章湍流模型 本章主要介绍Fluent所使用的各种湍流模型及使用方法。 各小节的具体内容是: 10.1 简介 10.2 选择湍流模型 10.3 Spalart-Allmaras 模型 10.4 标准、RNG和k-e相关模型 10.5 标准和SST k-ω模型 10.6 雷诺兹压力模型 10.7 大型艾迪仿真模型 10.8 边界层湍流的近壁处理 10.9 湍流仿真模型的网格划分 10.10 湍流模型的问题提出 10.11 湍流模型问题的解决方法 10.12 湍流模型的后处理 10.1 简介 湍流出现在速度变动的地方。这种波动使得流体介质之间相互交换动量、能量和浓度变化,而且引起了数量的波动。由于这种波动是小尺度且是高频率的,所以在实际工程计算中直接模拟的话对计算机的要求会很高。实际上瞬时控制方程可能在时间上、空间上是均匀的,或者可以人为的改变尺度,这样修改后的方程耗费较少的计算机。但是,修改后的方程可能包含有我们所不知的变量,湍流模型需要用已知变量来确定这些变量。 FLUENT 提供了以下湍流模型: ·Spalart-Allmaras 模型 ·k-e 模型 -标准k-e 模型 -Renormalization-group (RNG) k-e模型 -带旋流修正k-e模型 ·k-ω模型 -标准k-ω模型 -压力修正k-ω模型 -雷诺兹压力模型 -大漩涡模拟模型 10.2 选择一个湍流模型 不幸的是没有一个湍流模型对于所有的问题是通用的。选择模型时主要依靠以下几点:流体是否可压、建立特殊的可行的问题、精度的要求、计算机的能力、时间的限制。为了选择最好的模型,你需要了解不同条件的适用范围和限制 这一章的目的是给出在FLUENT中湍流模型的总的情况。我们将讨论单个模型对cpu 和内存的要求。同时陈述一下一种模型对那些特定问题最适用,给出一般的指导方针以便对于你需要的给出湍流模型。 10.2.1 雷诺平均逼近vs LES 在复杂形体的高雷诺数湍流中要求得精确的N-S方程的有关时间的解在近期内不太可能实现。两种可选择的方法用于把N-S方程不直接用于小尺度的模拟:雷诺平均和过滤。

fluent学习心得

1. 分离式求解器和耦合式求解器:都适用于从不可压到高速可压的很大范围的流动,总得来说,计算高速可压时,耦合式求解器更有优势;分离式求解器中有几个模型耦合式求解器中没有,如VOF,多项混合模型等。 2. 对于绝大多数问题,选择1st-Order Implicit就已经足够了。精度要求高时,选择2st-Order Implicit.而Explicit选项只对耦合显式求解器有效。 3. 压力都是相对压力值,相对于参考压力而言。对于不可压流动,若边界条件中不包含有压力边界条件时,用户应设置一个参考压力位置。计算时,fluent强制这一点的相对压力值为0. 4. 选择什么样的求解器后,再选择什么样的计算模型,即通知fluent是否考虑传热,流动是无粘、层流还是湍流,是否多相流,是否包含相变等。默认情况,fluent只进行流场求解,不求解能量方程。 5. 多相流模型:其中vof模型通过单独的动量方程和处理穿过区域的每一流体的容积比来模拟两种或三种不能混合的流体。 6. 能量方程:选中表示计算过程中要考虑热交换。对于一般流动,如水利工程及水力机械流场分析,可不考虑传热;气流模拟时,往往要考虑。默认状态下,fluent在能量方程中忽略粘性生成热,而耦合式求解器包含有粘性生成热。 7. 粘性模型:inviscid无粘计算;Laminar模型,层流模型;k-epsilon(2 eqn)模型,目前常用模型。 8. 材料定义:比较简单 9. 边界条件:见P210-211 10. 给定湍流参数:在计算区域的进口、出口及远场边界,需给定输运的湍流参数。Turbulence specification Method项目,意为让用户指定使用哪种模型来输入湍流参数。用户可任选其一,然后按公式计算选定的湍流参数,并作为输入。 湍流强度,湍动能k,湍动耗散率e。 11. 常用的边界条件: 压力进口:适用于可压和不可压流动,用于进口的压力一直但流量或速度未知的情况。Fluent 中各种压力都是相对压力值。 速度入口:用于不可压流,如果用于可压流可能导致非物理结果。 质量进口:规定进口的质量。 压力出口:需要在出口边界处设置静压。静压只用于亚音速流动。在fluent求解时,当压力出口边界上流动反向时,就是用这组回流条件。出口回流有三种方式:垂直与边界,给定方向矢量,来自相邻单元。 出流:用于模拟求解前流速和压力未知的出口边界。适用于出流面上的流动情况由区域内外推得到,且对上游没影响。不用于可压流动,也不能与压力进口边界条件一起是用。 压力远场:只适用于可压气体流动,气体的密度通过理想气体定律来计算。 12. 设置求解控制参数:为了更好的控制求解过程,需要在求解器中进行某些设置,内容包括选择离散格式、设置欠松弛因子、初始化场变量及激活监视变量等。 Fluent允许用户对流项选择不同的离散格式。默认情况下,当是用分离式求解器时,所有方程中的对流相一阶迎风格式离散;耦合式求解时,二阶精度格式,其他仍一阶。对于2D三角形和3D四面体网格,注意要是用二阶精度格式。一般,一阶容易收敛,精度差。 欠松弛因子:为了加速收敛,在迭代10次左右后,检查残差是增加还是减小,若增大,则减小欠松弛因子的值;反之,增大它。 Pressure-velocity coupling:包含压力速度耦合方式的列表。该项只在分离式求解器中出现。可选SIMPLE、SIMPLEC、PISO。多数选择simplec,piso算法主要用于瞬态问题的模拟,

相关文档