湍流与燃烧课程相关报告.docVIP

《湍流与燃烧》课程 相关报告 王春海 硕0902班 2010年4月13日 目录 《湍流与燃烧》课程 1 相关报告 1 目录 2 《湍流与燃烧》课程相关报告 4 第一部分 流动数值模拟过程概述 4 第二部分 流动基本控制方程 5 1 连续性方程 5 1.1 方程物理意义 5 1.2 方程的形式 5 1.3 方程的推导思路及过程 6 1.4 微分方程的主要形式及其简化 8 2 动量方程 9 2.1 方程物理意义 9 2.2 方程的推导思路及过程 9 2.3 微分方程的主要形式及其简化 11 3 能量方程 13 3.1 方程物理意义 13 3.2 方程的推导思路及过程 13 3.3 微分方程的主要形式及其简化 14 4 组分方程 16 4.1 方程物理意义 16 4.2 方程的推导思路及过程 16 第三部分 湍流模拟——大涡模拟 18 1大涡模拟的产生背景 18 2 大涡模拟的基本思想 18 3 大涡模拟的关键过程 19 3.1 过滤过程[5] 19 3.2 亚网格尺度模型的建立 21 4 大涡模拟的适用范围及优缺点分析 21 5 大涡模拟的优化和改进 22 5.1 Smargorinsky模型 22 5.2 Germano模型（动态Smargorinsky模型） 23 5.3 结构相似性模型 23 第四部分 湍流燃烧模型——涡团破碎模型 24 1涡团破碎模型概述 24 1.1 涡团破碎模型产生背景 24 1.2 涡团破碎模型的基本思想 24 1.3 涡团破碎模型的基本假设 24 2涡团破碎模型的推导思路及过程 24 3涡团破碎模型的优缺点及适用范围 25 4模型的修正 26 4.1 EDM模型（Eddy Dissipation Model） 26 4.2 EDC模型（Eddy Dissipation Concept Model） 26 4.3 特征时间模型（混合模型） 27 4.4 拉切滑模型 27 4.5 EBU系列模型的比较 27 参考文献： 29 《湍流与燃烧》课程相关报告 硕0902 王春海 摘要：本文对流动数值模拟的几个关键部分进行了介绍，主要内容包括基本控制方程的推导、湍流大涡数值模拟基本思想的介绍及涡团破碎系列模型的介绍三个部分。 关键词：基本控制方程 湍流大涡数值模拟 涡团破碎模型 第一部分 流动数值模拟过程概述  第二部分 流动基本控制方程 流动的基本控制方程的建立均是基于连续介质假设，即物质连续地无间隙地分布于物质所占有的整个空间，流体宏观物理量是空间点及时间的连续函数。 1 连续性方程 1.1 方程物理意义 流体运动连续性方程是质量守恒定律在流体力学中的应用。 其具体含义是包含在一流体系统中的流体质量在运动过程中保持不变。另一种说法是，在一固定空间中的流体质量的减少率等于在此期间通过其表面的质量通量[1]。 1.2 方程的形式 流动基本控制方程的分类如下所示： 具体到连续方程，其分类如下： 其中微分形式的欧拉型方程在流动计算过程中应用较多，本文中的基本控制方程均指这一类。 1.3 方程的推导思路及过程 1.3.1 基于雷诺输运定理的推导方式 根据雷诺输运定理[1]有： (2.1.1) 式2.1.1的具体推导过程见参考文献[1]。 应用质量守恒定律有： (2.1.2) 其中 (2.1.3) 又由Gauss公式[2] (2.1.4) 由上述三式可得： (2.1.5) 故 (2.1.6) 经数学变换可得： (2.1.7) 1.3.2 基于六面体外节点（图1-1中A点）的推导方式 图 1-1 控制体示意图 如图1-1所示，在流场中取一固定不动的微平行六面体，在直角坐标系Oxyz中，微六面体边长分别为dx，dy，dz。依据质量守恒定律，建立连续性方程： t时刻 A点流体密度ρ(x,y,z,t)，速度ν(x,y,z,t)，沿x, y, z轴分量是u, v, w。 dt时间内从左边微元面积dydz流入的流体质量为： （2.1.8） 从右边微元面积dydz流出的流体质量为： （2.1.9） dt时间内沿x轴方向流体质量的变化（净流出）为： （2.1.10） 同理，dt时间内沿y、z轴方向流体质量的变化（净流出）分别为： （2.1.11） （2.1.12） 所以，dt时间内经微六面体表面的流体质量的尽通量为： （2.1.13） dt时间内微六面体内质量的变化为： （2.1.14） 根据质量守恒有： （2.1.15） 即 （2.1.16） 即： （2.1.17） 1.3.3 基于六面体中心点（图1-1中六面体中心点P）的推导方式 除了1.