1、 基础第1页，共43页。主要内容n 发展历史n 基本功能n 软件结构及常用文件类型n 安装方法（基于6.3版本）n 基本概念n 求解器的选择n 离散格式n 初始化n 边界条件的设定n 可用的参考资料第2页，共43页。发展历史n 1975年 谢菲尔德大学（）开发了n 1983年 美国的流体技术服务公司推出 n 1988年 .成立 n 1995年 收购最大对手公司（）n 1997年 收购公司（粘弹性和聚合物流动模拟）n 2006年 被收购 n 在被收购后为6.3版本n 2009年6月发布12.0版本n 2010年底发布13.0版本n 2011年底发布14.0版本第3页，共43页。基本功能（1）n
2、可压缩与不可压缩流动问题；n 稳态和瞬态流动问题；n 无粘流、层流及湍流问题；n 牛顿流体及非牛顿流体；n 对流换热问题（包括自然对流和混合对流）；n 导热与对流换热耦合问题；n 辐射换热；n 惯性坐标系和非惯性坐标系下的流动问题模拟；第4页，共43页。基本功能（2）n 多运动坐标系下的流动问题；n 化学组分混合与反应；n 可以处理热量、质量、动量和化学组分的源项；n 用 轨道模型模拟稀疏相（颗粒，水滴，气泡等）；n 多孔介质流动；n 一维风扇、热交换器性能计算；n 两相流问题；n 复杂表面形状下的自由面流动。第5页，共43页。软件结构及常用文件类型（1）n 包应该包括以下几个部分：n 解法器
3、n，用于模拟燃烧过程n，网格生成n，额外的处理器，用于从现有的边界网格生成体网格。n()，转换其它程序生成的网格，用于计算。可以接口的程序包括：,等。第6页，共43页。软件结构及常用文件类型（2）GAMBIT设置几何形状生成2D或3D网格其它软件包，如CAD，CAE等FLUENT网格输入及调整物理模型边界条件流体物性确定计算结果后处理TGrid2D三角网格3D四面体网格2D和3D混合网格prePDFPDF查表2D或3D网格几何形状或网格PDF程序网格边界和（或）体网格边界网格基本程序结构示意图第7页，共43页。软件结构及常用文件类型（3）n 软件文件类型：n 文件：日志文档，可以编辑运行。n
4、文件：工作文件，若想修改网格，可以打开进行再编辑,可以打开进行再编辑。n 文件：输出的网格文件。n 文件：文件经过处理后得到的文件。n 文件：计算数据结果的数据文件。第8页，共43页。安装方法n最初只支持操作系统，面向高端用户。随着的普及，也开始支持系统。但其前处理软件仍然需要在系统下运行。所以在安装前，要安装一个模拟系统的软件。n 基于6.3的软件安装及破解方法：n 1.安装之前，时间调整到2005年；n 2.先装（安装文件在X:9.0），再装，最后装；n 3.按照提示，复制相应的许可证到和的目录里；n 4.设置初始环境；n5.设置完初始环境之后，把时间调整过来，重启电脑！第9页，共43页。
5、基本概念（1）n 理想与实际流体：根据是否考虑流体的粘性，可将流体分为理想流体和实际流体。n 可压缩流体与不可压缩流体：根据流体压缩性的大小，可将流体分为可压缩流体与不可压缩流体。密度随压强变化大且不可视为常数的为可压缩流体，反之，称为不可压缩流体。n 正常情况下，液体和低速气体（50）可视为不可压缩流体。n 牛顿流体和非牛顿流体：根据流体是否满足牛顿内摩擦定律，可将流体分为牛顿流体和非牛顿流体。符合牛顿内摩擦定律，也就是切应力与速度梯度成正比的流体为牛顿流体。n 水、酒精等大多数纯液体、轻质油、低分子化合物溶液以及低速流动的气体等均为牛顿流体；n 非牛顿流体：油漆、稀润滑脂、高分子溶液、乳化
6、液、淀粉糊、石灰浆、牛奶、牙膏、血液、泥浆等。第10页，共43页。基本概念（2）n 层流和湍流：层流是流体运动规则，各部分分层流动互不掺混，质点的轨线是光滑的，而且流动稳定。湍流的特征则完全相反，流体运动极不规则，各部分激烈掺混，质点的轨线杂乱无章，而且流场极不稳定。这两种截然不同的运动形态在一定条件下可以相互转化。n 可用雷诺数进行判定。雷诺数是流体惯性力与黏性力比值的量度，它是一个无量纲量。n 雷诺数较小时，黏滞力对流场的影响大于惯性力，流场中流速的扰动会因黏滞力而衰减，流体流动稳定，为层流；反之，若雷诺数较大时，惯性力对流场的影响大于黏滞力，流体流动较不稳定，流速的微小变化容易发展、增强
7、，形成紊乱、不规则的紊流流场。第11页，共43页。基本概念（3）n 定常流动和非定常流动：以时间为标准，根据流体流动的物理量（如速度、压力、温度等）是否随时间变化，将流动分为定常与非定常两大类。当流动的物理量不随时间变化，为定常流动；反之称为非定常流动。n 许多流体机械在起动或关机时的流体流动一般是非定常流动，而正常运转时可看作是定常流动。第12页，共43页。基本概念（4）n 亚音速流动与超音速流动：当气流速度很大，或者流场压力变化很大时，流体就受到了压缩性的影响。n 马赫数是流体速度与当地音速之比值，而音速在不同高度、温度等状态下又有不同数值，因此无法将 的数值换算为固定的 或 等单位。n
8、1时，为亚音速流动；n 1，超音速流动。第13页，共43页。基本概念（5）n 热传导及扩散：当流体中存在温度差时，温度高的地方将向温度低的地方传送热量，这种现象称为热传导。同样地，当流体混合物中存在组元的浓度差时，浓度高的地方将向浓度低的地方输送该组元的物质，这种现象称为扩散。n 理想流体忽略了粘性，即忽略了分子运动的动量输运性质，因此在理想流体中也不应考虑质量和能量输运性质扩散和热传导，因为它们具有相同的微观机制。第14页，共43页。基本概念（6）n 有旋流动和无旋流动：有旋流动是指流场中各处的旋度(流体微团的旋转角速度)不等于零的流动，无旋流动是指流场中各处的旋度都为零的流动。n 流体运动
9、是有旋流动还是无旋流动，取决于流体微团是否有旋转运动，与流体微团的运动轨迹无关。流体流动中，如果考虑粘性，由于存在摩擦力，这时流动为有旋流动：如果粘性可以忽略，而流体本身又是无旋流，这时流动为无旋流动。n 例如，均匀气流流过平板，在紧靠壁面的附面层内，需要考虑粘性影响，因此，附面层内为有旋流动，附面层外的流动，就性可以忽略，为无旋流动。第15页，共43页。求解器选择（1）n 求解器：n(1)2d二维单精度求解器；n(2)3d三维单精度求解器：n(3)2二维双精度求解器；n(4)3三维双精度求解器。第16页，共43页。求解器选择（2）n中有两种求解器 压力基和密度基。n压力基求解器以动量和压力为
10、基本变量n通过连续性方程导出压力和速度的耦合算法n只采用隐式方式对控制方程进行线性化n压力基求解器有两种算法n分离求解器 压力修正和动量方程顺序求解。n耦合求解器()压力和动量方程同时求解 用哪种压力梯度来计算控制方程中的导数项。第17页，共43页。求解器选择（3）n压力基求解器通过连续性方程和动量方程导出压力方程或压力修正方程n压力-速度耦合方程求解算法n ()n默认算法，稳健性好 ()对简单问题，收敛更快，如层流 ()对非稳态流动或者高扭曲度网格有用 压力基耦合求解器第18页，共43页。求解器选择（4）n密度基耦合求解器n以矢量方式求解连续性方程、动量方程、能量方程和组分方程n通过状态方程
11、得到压力n其他标量方程按照分离方式求解n 可以显式或隐式方式求解n隐式:使用高斯赛德尔方法求解所有变量n显式:用多步龙格库塔显式时间积分法。第19页，共43页。求解器选择（5）第20页，共43页。求解器选择（6）第21页，共43页。求解器选择（7）n压力基求解器应用范围覆盖从低压不可压缩流到高速压缩流n需要的内存少n求解过程灵活n压力基耦合求解器()适用于大多数单相流，比分离求解器性能更好n不能用于多相流（欧拉）、周期质量流和 n比分离求解器多用1.52倍内存n密度基耦合求解器()适用于密度、能量、动量、组分间强耦合的现象n例如:伴有燃烧的高速可压缩流动，超高音速流动、激波干扰n隐式方法一般优
12、于显式，因为其对时间步有严格的限制n显式方法一般用于流动时间尺度和声学时间尺度相当的情况（如高马赫激波的传播）第22页，共43页。离散格式（1）n 在对指定问题进行计算之前，首先要将计算区域离散化，即对空间上连续的计算区域进行划分，把它划分成许多个子区域，并确定每个区域中的节点，从而生成网格。然后将控制方程在网格上离散，即将偏微分方程格式的控制方程转化为各个节点上的代数方程组，然后求解离散方程组，得到各个节点上的解。n 节点之间的近似解，一般认为光滑变化，原则上可以应用插值方法确定，从而得到整个计算域上的近似解。n 这种插值方式被称为离散格式（）第23页，共43页。离散格式（2）第24页，共4
13、3页。离散格式（3）n对流项的插值方法有:n 易收敛，一阶精度n 对低雷诺数流动(第33页，共43页。边界条件的设定（3）表压表压operating pressure压力压力水平水平operating pressure绝对压力绝对压力真空真空如果流动是亚声速的，则会忽略如果流动是亚声速的，则会忽略 ，它是由指定的驻点值来计算的。，它是由指定的驻点值来计算的。这里给出的表压强的大小，是入口边界上的总压。这里给出的表压强的大小，是入口边界上的总压。第34页，共43页。边界条件的设定（4）n 质量流量进口边界条件()：给定入口边界上的质量流量。主要用于可压缩流动问题，对于不可压缩问题，由于密度是常数
14、，可以使用速度入口条件。如果压力边界条件和质量边界条件都适合流动时，优先选择用压力进口条件。n 压力出口边界条件()：需要给定出口静压（表压）。而且，该压力只用于亚音速计算。如果局部变成超音速，则根据前面来流条件外推出口边界条件。需要特别指出的是，这里的压力是相对于前面给定的工作压力。n 如果有回流出现，给的表压将视为总压，所以不必给出回流压力。回流流动方向与出口边界垂直。第35页，共43页。边界条件的设定（5）第36页，共43页。边界条件的设定（6）n 压力远场边界条件():如果知道来流的静压和马赫数，提供了的压力远场边界条件来模拟该类问题。n 该边界条件只适合用理想气体定律计算密度的问题，
15、而不能用于其它问题。为了满足压力远场条件，需要把边界放到我们关心区域足够远的地方。n 给定边界静压和温度及马赫数。可以是亚音速，跨音速或者超音速，并且需要给定流动方向，如果有需要还必须给定湍流量等参数。第37页，共43页。边界条件的设定（7）第38页，共43页。边界条件的设定（8）n 自由流出边界条件():不知道流出口的压力或者速度，这时候可以选择流出边界条件。n 如下列情况，不能用流出边界条件：n 包含压力进口条件n 可压速流动问题n 有密度变化的非稳定流动问题（即使是不可压速流动）第39页，共43页。边界条件的设定（9）n 固壁边界条件（）:对于粘性流动问题，默认设置是壁面无滑移条件。n 壁面热边界条件包括固定热通量、固定温度、对流换热系数、外部辐射换热、外部辐射换热与对流换热等。第40页，共43页。可用的参考资料n 参考书：n 于勇：入门与进阶教程n 韩占忠：流体工程仿真计算实例与分析n 李进良：精通 6.3流场分析n 王福军：计算流体动力学分析软件原理与应用n n 论坛：第41页，共43页。第42页，共43页。演讲完毕，谢谢观看！第43页，共43页。
