气相湍流模拟目前大致分为三大类:微观模拟、概率密度函数模拟、统观模拟

气相湍流模拟目前大致分为三大类:微观模拟、概率密度函数模拟、统观模拟:

 

微观模拟

     目前,基于连续介质理论体系下的湍流微观模型包括直接求解N—S方程(DNS)、大涡模拟(LES)、离散涡模拟(DVS)三类。

      DNS方法是在湍流尺度下的网格尺寸内不引入任何封闭模型的前提下对N—S方程直接求解。这种方法能对湍流流动中最小尺度涡进行求解,其所得结果的误差仅是一般数值计算所引起的那些误差,并且可以据需要而加以控制,但必须采用很小的时间与空间步长以达到统计的稳定状态,因此,DNS方法仅限于很小空间内相对低的雷诺数和简单边界条件的湍流流动模拟,如槽道流动与平板边界层流动;另外由于DNS方法计算量很大,目前尚无法用于工程问题。

      LES方法的基本思想是:湍流流动是有许多大小不同尺度的涡旋组成,大尺度的涡旋对平均流动影响比较大,各种变量的湍流扩散、热量、质量、动量和能量的交换以及雷诺应力的产生都是通过大尺度涡旋来实现的,而小尺度涡旋主要对耗散起作用,通过耗散脉动来影响各种变量。由于网格尺度比湍流尺度大,可以模拟湍流发展过程中的一些细节,用于比较简单的剪切流运动和管流。由于LES方法计算精度可与DNS方法相比,LES方法可得到真实瞬态流场精度相当高,但其计算费用却便宜得多,所以近几十年来得到飞快发展。目前该模型主要用于检验统观模型。

      离散涡或随机涡模拟是把湍流流场分成一系列大尺度涡元,用涡元的随机运动来模拟湍流。它能获得混合层和尾涡中大尺度涡旋机理的本质,目前也得到了一定程度的应用[54]。总的来说,离散涡模拟(DVS)仍是比较新的求解N—S方程的一种数值方法,尚处于形成阶段。最主要的原因是定性上描述尚可而定量上精度难以确定。就已得到的结果而言,只能给出合理的现象。其主要问题是:(1)涡旋尺寸分布的确定有任意性;(2)未考虑含能涡和耗散涡的差别;(3)尚难以用于三维问题。

概率密度函数模拟

      近年来,由分子运动理论推导的概率密度函数(Probability Density Funtion)输运方程,简称PDF模型。与通常的湍流模型相比,PDF模型能够对湍流流场进行更为精确和完整的描述。应用PDF模型计算的平面混合层流流动、平面尾流、平面射流和轴对称射流四种剪切流动的扩展率、平均速度的分布和Reynolds应力的计算值与试验值吻合较好。用这种方法对湍流自由射流进行模拟,得到的射流扩张率、平均速度及Reynolds应力与试验值十分吻合。采用PDF模型模拟湍流燃烧和污染物的生成的较多。

统观模拟

      在单相湍流领域内,上述两种模拟方法在工程应用上的困难,按国内外一些学者的看法,可用于工程应用上的现实模拟方法是雷诺方程组与关联矩封闭方法。它将雷诺时均方程及湍流特征量输运方程中的高阶未知关联项用低阶未知关联项或时均量来表达,从而使雷诺时均方程封闭。对于这些关联量的处理,有各种模型:湍流粘性系数模型、代数雷诺应力模型(ASM)、雷诺应力模型(DSM)。而其中较为常见、较为有效的是k-ε双方程模型和雷诺应力模型。k-ε双方程模型相对来说较为简单,使用方便,适合工程问题的研究,目前应用得最多,对其解决各种流动问题的有效性也认识得最清楚。

      k-ε双方程模型有以下优点:(1)通过求偏微分方程考虑湍流物理量的输运过程,即通过求解偏微分方程确定脉动特征速度与平均场速度梯度的关系,而不是直接将两者联系起来;(2)特征长度不是由经验确定,而是以耗散尺度作为特征长度,并由求解相应的偏微分方程得到。

      k-ε双方程模型在应用时也有其局限性:它不能用于低雷诺数的层流条件,适合用于高雷诺数、充分发展的湍流流动区域。基于各向同性的标准k-ε双方程模型可以较成功的模拟无浮力平面射流、平壁边界层、管流及有弱旋的回流流动;而对强旋流、浮力流、重力分层流以及曲壁边界层和低雷诺数流动问题,有很大的缺陷。这主要有两方面的原因:一方面是因为k-ε模型中ε方程模拟的不准确性,另一方面是因为采用了各向同性涡粘系数的Boussinesq假设,无法预报各向异性很强的湍流流动。一般认为Sn(旋流数)﹤0.5以下,采用各向同性的湍流粘性系数模型,这时标准的k-ε双方程模型能较好地预报流场内的分布特性,当Sn﹥1时,标准的k-ε双方程模型的预报是不成功的,必须对标准的k-ε双方程模型进行修正。

一本目前为止最好的fluent学习书本 第一章 流体力学基础与FLUENT简介 第一节 概论 一、流体的密度、重度和比重 二、流体的黏性——牛顿流体与非牛顿流体 三、流体的压缩性——可压缩与不可压缩流体 四、液体的表面张力 第二节 流体力学中的力与压强 一、质量力与表面力 二、绝对压强、相对压强与真空度 三、液体的汽化压强 四、静压、动压和总压 第三节 能量损失与总流的能量方程 一、沿程损失与局部损失 二、总流的伯努里方程 三、人口段与充分发展段 第四节 流体运动的描述 一、定常流动与非定常流动 二、流线与迹线 三、流量与净通量 四、有旋流动与有势流动 五、层流与湍流 第五节 亚音速与超音速流动 一、音速与流速 二、马赫数与马赫锥 三、速度系数与临界参数 四、可压缩流动的伯努里方程 五、等熵滞止关系式 第六节 正激波与斜激波 一、正激波 二、斜激波 第七节 流体多维流动基本控制方程 一、物质导数 二、连续性方程 三、N—S方程 第八节 边界层与物体阻力 一、边界层及基本特征 二、层流边界层微分方程 三、边界层动量积分关系式 四、物体阻力 第九节 湍流模型 第十节 FLUENT简介 一、程序的结构 二、FLUENT程序可以求解的问题 三、用FLUENT程序求解问题的步骤 四、关于FLUENT求解器的说明 五、FLUENT求解方法的选择 六、边界条件的确定 第二章 二维流动与传热的数值计算 第一节 冷、热水混合器内部二维流动 一、前处理——利用GAMBIT建立计算模型 第1步 确定求解器 第2步 创建坐标网格图 第3步 由节点创建直线 第4步 创建圆弧边 第5步 创建小管嘴 第6步 由线组成面 第7步 确定边界线的内部节点分布并创建结构化网格 第8步 设置边界类型 第9步 输出网格并保存会话 二、利用FLUENT进行混合器内流动与热交换的仿真计算 第1步 与网格相关的操作 第2步 建立求解模型 第3步 设置流体的物理属性 第4步 设置边界条件 第5步 求解 第6步 显示计算结果 第7步 使用二阶离散化方法重新计算 第8步 自适应性网格修改功能 小结 课后练习 第二节 喷管内二维非定常流动 一、利用GAMBIT建立计算模型 第1步 确定求解器 第2步 创建坐标网格图和边界线的节点 第3步 由节点创建直线 第4步 利用圆角功能对I点处的角倒成圆弧 第5步 由边线创建面 第6步 定义边线上的节点分布 第7步 创建结构化网格 第8步 设置边界类型 第9步 输出网格并保存会话 二、利用FLUENT进行喷管内流动的仿真计算 第1步 与网格相关的操作 第2步 确定长度单位 第3步 建立求解模型 第4步 设置流体属性 第5步 设置工作压强为0 atm 第6步 设置边界条件 第7步 求解定常流动 第8步 非定常边界条件设置以及非定常流动的计算 第9步 求解非定常流 第10步 对非定常流动计算数据的保存与后处理 小结 课后练习 第三节 三角翼的可压缩外部绕流 一、利用GAMBIT建立计算模型 第1步 启动Gambit,并选择求解器为FLUENT5/6 第2步 创建节点 第3步 由节点连成线 第4步 由边线创建面 第5步 创建网格 第6步 设置边界类型 第7步 输出网格文件 二、利用FLUENT进行仿真计算 第1步 启动FLUENT 2D求解器并读入网格文件 第2步 网格检查与确定长度单位 第3步 建立计算模型 第4步 设置流体材料属性 第5步 设置工作压强 第6步 设置边界条件 第7步 利用求解器进行求解 第8步 计算结果的后处理 小结 课后练习 第四节 三角翼不可压缩的外部绕流(空化模型应用) 第1步 启动FLUENT 2D求解器并读入网格文件 第2步 网格检查与确定长度单位 第3步 设置求解器 第4步 设置流体材料及其物理性质 第5步 设置流体的流相 第6步 设置边界条件 第7步 求解 第8步 对计算结果的后处理 小结 课后练习 第五节 VOF模型的应用 一、利用GAMBIT建立计算模型 第1步 启动GAMBIT并选择FLUENT5/6求解器 第2步 建立坐标网格并创建节点 第3步 由节点连成直线段 第4步 创建圆弧 第5步 创建线段的交点G 第6步 将两条线在G点处分别断开 第7步 删除DG直线和FG弧线 第8步 由边创建面 第9步 定义边线上的节点分布 第10步 在面上创建结构化网格 第11步 设置边界类型 第12步 输出网格文件并保存会话 二、利用FLUENT 2D求解器进行求解 第1步 读入、显示网格并设置长度单位 第2步 设置求解器 第3步 设置流体材料及属性 第4步 设置基本相和第二相 第5步 运算环境设置 第6步 设置边界条件 第7步 求解 第8步 计算结果的后处理 小结 第六节 组分传输与气体燃烧 一、利用GAMBIT建立计算模型 第1步 打开GAMBIT 第2步 对空气进口边界进行分网 第3步 设置边界条件 第4步 输出2D网格 二、利用FLUENT-2D求解器进行 模拟计算 第1步 与网格相关的操作 第2步 设置求解模型 第3步 流体材料设置 第4步 边界条件设置 第5步 使用常比热容的初始化并求解 第6步 采用变比热容的解法 第7步 后处理 第8步 NOx预测 小结 第三章 三维流动与传热的数值计算 第一节 冷、热水混合器内的三维流动与换热 一、利用GAMBIT建立混合器计算模型 第1步 启动GAMBIT并选定求解器(FLUENT5/6) 第2步 创建混合器主体 第3步 设置混合器的切向入流管 第4步 去掉小圆柱体与大圆柱体相交的多余部分,并将三个圆柱体联结成一个整体 第5步 创建主体下部的圆锥 第6步 创建出流小管 第7步 将混合器上部、渐缩部分和下部出流小管组合为一个整体, 第8步 对混合器内区域划分网格 第9步 检查网格划分情况 第10步 设置边界类型 第11步 输出网格文件(.msh) 二、利用FLUENT 3D求解器进行求解 第1步 检查网格并定义长度单位 第2步 创建计算模型 第3步 设置流体的材料属性 第4步 设置边界条件 第5步 求解初始化 第6步 设置监视器 第7步 保存Case文件 第8步 求解计算 第9步 保存计算结果 三、计算结果的后处理 第1步 读入Case和Data文件 第2步 显示网格 第3步 创建等(坐标)值面 第4步 绘制温度与压强分布图 第5步 绘制速度矢量图 第6步 绘制流体质点的迹线 第7步 绘制XY曲线 小结 课后练习 第二节 粘性流体通过圆管弯头段的三维流动 一、前处理——利用GAMBIT建立计算模型 第1步 确定求解器 第2步 创建圆环 第3步 创建立方体 第4步 移动立方体 第5步 分割圆环 第6步 删除3/4圆环 第7步 建立弯管直段 第8步 移动弯管直段 第9步 整合弯管和直段 第10步 边界层的设定 第11步 划分面网格 第12步 划分体网格 第13步 定义边界类型 第14步 输出网格文件 二、利用FLUENT 3D求解器进行模拟计算 第1步 启动FLUENT,进入3D模式 第2步 读入网格数据 第3步 网格检查 第4步 显示网格 第5步 建立求解模型 第6步 设置标准湍流模型 第7步 设置流体的物理属性 第H步 设置边界条件 第9步 求解控制 第10步 求解 第11步 显示初步计算结果 第12步 流线显示 小结 第三节 三维稳态热传导问题 一、利用GAMBIT进行网格划分 第1步 导入几何模型 第2步 选取求解器 第3步 网格划分 第4步 边界条件设置 第5步 网格检查 第6步 输出网格 二、利用FLUENT-3D求解器进行数值模拟计算 第1步 在FLUENT中读入网格文件 第2步 选取求解器 第3步 材料设置 第4步 边界条件 第5步 求解控制 第6步 后处理 小结 第四节 动网格问题 一、利用FLUENT-3D进行计算 第1步 与网格有关的操作 第2步 模型没置 第3步 材料设置 第4步 边界条件设置 第5步 网格运动设置 第6步 求解 二、利用FLUENT—3D进行后处理 第1步 检查最后一个时间步(BDC)的解 第2步 检查上死点的解 第3步 回放温度等高线动画 第4步 显示上死点时缸内的流动矢量切面 小结 第五节 叶轮机械的Mixing Plane模型 —、利用FLUENT-3D求解器进行计算 第1步 网格 第2步 单位设置 第3步 计算模型设置 第4步 混合面(Mixing Plane)设置 第5步 流体材料设置 第6步 边界条件设置 第7步 求解 二、利用FLUENT-3D进行后处理 第1步 生成后处理的—个等值画 第2步 显示速度矢量 第3步 平面x=0上绘全压的周向平均量 第4步 显示全压的等高线图 小结 附录 参考文献
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值