Fluent_入门基础培训.ppt

Fluent基础 主要内容 发展历史基本功能软件结构及常用文件类型安装方法 基于6 3版本 基本概念 求解器的选择离散格式初始化边界条件的设定可用的参考资料 发展历史 1975年谢菲尔德大学 UK 开发了Tempest1983年美国的流体技术服务公司creature推出fluent1988年FluentInc 成立1995年收购最大对手FDI公司 FIDAP 1997年收购Polyflow公司 粘弹性和聚合物流动模拟 2006年被ANSYS收购在被ANSYS收购后为6 3版本2009年6月发布12 0版本2010年底发布13 0版本2011年底发布14 0版本 基本功能 1 可压缩与不可压缩流动问题 稳态和瞬态流动问题 无粘流 层流及湍流问题 牛顿流体及非牛顿流体 对流换热问题 包括自然对流和混合对流 导热与对流换热耦合问题 辐射换热 惯性坐标系和非惯性坐标系下的流动问题模拟 基本功能 2 多运动坐标系下的流动问题 化学组分混合与反应 可以处理热量 质量 动量和化学组分的源项 用Lagrangian轨道模型模拟稀疏相 颗粒 水滴 气泡等 多孔介质流动 一维风扇 热交换器性能计算 两相流问题 复杂表面形状下的自由面流动 软件结构及常用文件类型 1 FLUENT包应该包括以下几个部分 FLUENT解法器prePDF 用于模拟PDF燃烧过程GAMBIT 网格生成TGrid 额外的处理器 用于从现有的边界网格生成体网格 Filters Translators 转换其它程序生成的网格 用于FLUENT计算 可以接口的程序包括 ANSYS I DEAS NASTRAN PATRAN等 软件结构及常用文件类型 2 基本程序结构示意图 软件结构及常用文件类型 3 FLUENT软件文件类型 jou文件 日志文档 可以编辑运行 dbs文件 Gambit工作文件 若想修改网格 可以打开进行再编辑 可以打开进行再编辑 msh文件 Gambit输出的网格文件 cas文件 msh文件经过Fluent处理后得到的文件 dat文件 FLUENT计算数据结果的数据文件 安装方法 Fluent最初只支持UNIX操作系统 面向高端用户 随着Windows的普及 Fluent也开始支持Windows系统 但其前处理软件Gambit仍然需要在UNIX系统下运行 所以在安装Gambit前 要安装一个模拟UNIX系统的软件Exceed 基于Fluent6 3的软件安装及破解方法 1 安装之前 时间调整到2005年 2 先装exceed 安装文件在X Exceed v9 0 Exceed SETUP 再装gambit 最后装fluent 3 按照提示 复制相应的许可证license到fluent和gambit的目录里 4 设置初始环境 5 设置完初始环境之后 把时间调整过来 重启电脑 基本概念 1 理想与实际流体 根据是否考虑流体的粘性 可将流体分为理想流体和实际流体 可压缩流体与不可压缩流体 根据流体压缩性的大小 可将流体分为可压缩流体与不可压缩流体 密度随压强变化大且不可视为常数的为可压缩流体 反之 称为不可压缩流体 正常情况下 液体和低速气体 50m s 可视为不可压缩流体 牛顿流体和非牛顿流体 根据流体是否满足牛顿内摩擦定律 可将流体分为牛顿流体和非牛顿流体 符合牛顿内摩擦定律 也就是切应力与速度梯度成正比的流体为牛顿流体 水 酒精等大多数纯液体 轻质油 低分子化合物溶液以及低速流动的气体等均为牛顿流体 非牛顿流体 油漆 稀润滑脂 高分子溶液 乳化液 淀粉糊 石灰浆 牛奶 牙膏 血液 泥浆等 基本概念 2 层流和湍流 层流是流体运动规则 各部分分层流动互不掺混 质点的轨线是光滑的 而且流动稳定 湍流的特征则完全相反 流体运动极不规则 各部分激烈掺混 质点的轨线杂乱无章 而且流场极不稳定 这两种截然不同的运动形态在一定条件下可以相互转化 可用雷诺数进行判定 雷诺数是流体惯性力与黏性力比值的量度 它是一个无量纲量 雷诺数较小时 黏滞力对流场的影响大于惯性力 流场中流速的扰动会因黏滞力而衰减 流体流动稳定 为层流 反之 若雷诺数较大时 惯性力对流场的影响大于黏滞力 流体流动较不稳定 流速的微小变化容易发展 增强 形成紊乱 不规则的紊流流场 基本概念 3 定常流动和非定常流动 以时间为标准 根据流体流动的物理量 如速度 压力 温度等 是否随时间变化 将流动分为定常与非定常两大类 当流动的物理量不随时间变化 为定常流动 反之称为非定常流动 许多流体机械在起动或关机时的流体流动一般是非定常流动 而正常运转时可看作是定常流动 基本概念 4 亚音速流动与超音速流动 当气流速度很大 或者流场压力变化很大时 流体就受到了压缩性的影响 马赫数是流体速度与当地音速之比值 而音速在不同高度 温度等状态下又有不同数值 因此无法将Ma的数值换算为固定的km hr或mph等单位 Ma1 超音速流动 基本概念 5 热传导及扩散 当流体中存在温度差时 温度高的地方将向温度低的地方传送热量 这种现象称为热传导 同样地 当流体混合物中存在组元的浓度差时 浓度高的地方将向浓度低的地方输送该组元的物质 这种现象称为扩散 理想流体忽略了粘性 即忽略了分子运动的动量输运性质 因此在理想流体中也不应考虑质量和能量输运性质 扩散和热传导 因为它们具有相同的微观机制 基本概念 6 有旋流动和无旋流动 有旋流动是指流场中各处的旋度 流体微团的旋转角速度 不等于零的流动 无旋流动是指流场中各处的旋度都为零的流动 流体运动是有旋流动还是无旋流动 取决于流体微团是否有旋转运动 与流体微团的运动轨迹无关 流体流动中 如果考虑粘性 由于存在摩擦力 这时流动为有旋流动 如果粘性可以忽略 而流体本身又是无旋流 这时流动为无旋流动 例如 均匀气流流过平板 在紧靠壁面的附面层内 需要考虑粘性影响 因此 附面层内为有旋流动 附面层外的流动 就性可以忽略 为无旋流动 求解器选择 1 FLUENT求解器 1 FLUENT2d 二维单精度求解器 2 FLUENT3d 三维单精度求解器 3 FLUENT2ddp 二维双精度求解器 4 FLUENT3ddp 三维双精度求解器 求解器选择 2 FLUENT中有两种求解器 压力基和密度基 压力基求解器以动量和压力为基本变量通过连续性方程导出压力和速度的耦合算法只采用隐式方式对控制方程进行线性化压力基求解器有两种算法分离求解器 压力修正和动量方程顺序求解 耦合求解器 PBCS 压力和动量方程同时求解 用哪种压力梯度来计算控制方程中的导数项 求解器选择 3 压力基求解器通过连续性方程和动量方程导出压力方程或压力修正方程压力 速度耦合方程求解算法Semi ImplicitMethodforPressure LinkedEquations SIMPLE 默认算法 稳健性好 SIMPLE Consistent SIMPLEC 对简单问题 收敛更快 如层流Pressure ImplicitwithSplittingofOperators PISO 对非稳态流动或者高扭曲度网格有用Coupled压力基耦合求解器 求解器选择 4 密度基耦合求解器以矢量方式求解连续性方程 动量方程 能量方程和组分方程通过状态方程得到压力其他标量方程按照分离方式求解 DBCS可以显式或隐式方式求解隐式 使用高斯赛德尔方法求解所有变量显式 用多步龙格库塔显式时间积分法 求解器选择 5 求解器选择 6 求解器选择 7 压力基求解器应用范围覆盖从低压不可压缩流到高速压缩流需要的内存少求解过程灵活压力基耦合求解器 PBCS 适用于大多数单相流 比分离求解器性能更好不能用于多相流 欧拉 周期质量流和NITA比分离求解器多用1 5 2倍内存密度基耦合求解器 DBCS 适用于密度 能量 动量 组分间强耦合的现象例如 伴有燃烧的高速可压缩流动 超高音速流动 激波干扰隐式方法一般优于显式 因为其对时间步有严格的限制显式方法一般用于流动时间尺度和声学时间尺度相当的情况 如高马赫激波的传播 离散格式 1 在对指定问题进行CFD计算之前 首先要将计算区域离散化 即对空间上连续的计算区域进行划分 把它划分成许多个子区域 并确定每个区域中的节点 从而生成网格 然后将控制方程在网格上离散 即将偏微分方程格式的控制方程转化为各个节点上的代数方程组 然后求解离散方程组 得到各个节点上的解 节点之间的近似解 一般认为光滑变化 原则上可以应用插值方法确定 从而得到整个计算域上的近似解 这种插值方式被称为离散格式 discrezitationscheme 离散格式 2 离散格式 3 对流项的插值方法有 First OrderUpwind 易收敛 一阶精度PowerLaw 对低雷诺数流动 Recell 5 比一阶格式更精确Second OrderUpwind 尤其适用流动和网格方向不一致的四面体 三角形网格 二阶精度 收敛慢MonotoneUpstream CenteredSchemesforConservationLaws MUSCL 对非结构网格 局部三阶精度 对二次流 旋转涡 力等预测的更精确QuadraticUpwindInterpolation QUICK 适用于四边形 六面体以及混合网格 对旋转流动有用 在均匀网格上能达到三阶精度 离散格式 4 离散格式 5 为了得到扩散通量 速度导数 以及高阶离散格式 都需要求解变量的梯度单元中心的变量梯度由以下三种方法得到 Green GaussCell Based 可能会引起伪扩散Green GaussNode Based 更精确 更少伪扩散 建议对三角形 四面体网格采用Least SquaresCell Based 建议对多面体网格采用 精度和属性同Node based面上的梯度用多级泰勒级数展开求得 离散格式 6 使用分离算法时 计算面上压力的插值方法有 Standard 默认格式 对于近边界的沿面法向存在大压力梯度流动 精度下降 如果存在压力突变 建议改用PRESTO PRESTO 用于高度旋流 包括压力梯度突变 多孔介质 风扇模型等 或者计算域存在大曲率的面Linear 当其他格式导致收敛问题或非物理解时使用Second Order 用于压缩流 不适用多孔介质 风扇 压力突变以及VOF Mixture多相流BodyForceWeighted 用于大体积力的情况 如高瑞利数自然对流或高旋流 离散格式 7 初始化 1 FLUENT要求所有的求解变量有初始值更真实的初值能提高收敛稳定性 加速收敛过程 有些情况需要一个好的初值在特定区域对特定变量单独赋值自由射流 喷射区高速 燃烧问题 高温激活反应 单元标注 自适应 FMG初始化 2 FullMultiGrid FMG 能用来创建更好的初场FMG初始化对包括大的压力梯度和速度梯度的复杂流动有用在粗级别网格上求解一阶欧拉方程可用于压力基或密度基求解器 但限于稳态问题启动FMG初始化压力基求解器 solve init fmg initialization密度基求解器 当选择密度基求解器后在GUI里可见FMG在粗网格上用多重网格求解通过TUI命令来设置 solve init set fmg initialization 边界条件的设定 1 速度进口边界条件 velocity inlet 给出进口速度及需要计算的所有标量值 该边界条件适用于不可压缩流动问题 边界条件的设定 2 压力进口边界条件 pressure inlet 压力进口边界条件通常用于给出流体进口的压力和流动的其它标量参数 对计算可压和不可压问题都适合 压力进口边界条件通常用于不知道进口流率或流动速度时候的流动 这类流动在工程中常见 如浮力驱动的流动问题 压力进口条件还可以用于处理外部或者非受限流动的自由边界 压力边界条件的设置如图 其中第一项的表压强与绝对压强 参考压强有如下关系 Operatingpressure输入 Define operatingconditions 边界条件的设定 3 如果流动是亚声速的 则会忽略Supersonic InitialGaugePressure 它是由指定的驻点值来计算的 这里给出的表压强的大小 是入口边界上的总压 边界条件的设定 4 质量流量进口边界条件 mass flow inlet 给定入口边界上的质量流量 主要用于可压缩流动问题 对于不可压缩问题 由于密度是常数 可以使用速度入口条件 如果压力边界条件和质量边界条件都适合流动时 优先选择用压力进口条件 压力出口边界条件 pressure outlet 需要给定出口静压 表压 而且 该压力只用于亚音速计算 如果局部变成超音速 则根据前面来流条件外推出口边界条件 需要特别指出的是 这里的压力是相对于前面给定的工作压力 如果有回流出现 给的表压将视为总压 所以不必给出回流压力 回流流动方向与出口边界垂直 边界条件的设定 5 边界条件的设定 6 压力远场边界条件 pressure far field 如果知道来流的静压和马赫数 Fluent提供了的压力远场边界条件来模拟该类问题 该边界条件只适合用理想气体定律计算密度的问题 而不能用于其它问题 为了满足压力远场条件 需要把边界放到我们关心区域足够远的地方 给定边界静压和温度及马赫数 可以是亚音速 跨音速或者超音速 并且需要给定流动方向 如果有需要还必须给定湍流量等参数 边界条件的设定 7 边界条件的设定 8 自由流出边界条件 outf