可压缩湍流化学反应流的数值模拟研究
流体流动数值模拟
流体流动数值模拟流体流动现象普遍存在于⾃然界及多种⼯程领域中。
所有这些流动过程都遵循质量守恒、动量守恒、能量守恒和组分守恒等基本物理定律；⽽且流动若处于湍流状态，则该流动系统还要遵守附加的湍流输运⽅程。
本讲座将依据流体运动的特性阐述计算流体动⼒学的相关基础知识及任务；在流体运动所遵循的守恒定律及其数学描述的基础上，介绍数值求解这些基本⽅程的思想及其求解过程。
第⼀节计算流体动⼒学概述计算流体动⼒学（CFD）技术⽤于流体机械部流动分析及其性能预测，具有成本低，效率⾼，⽅便、快捷⽤时少等优点。
近年来随着计算流体⼒学和计算流体动⼒学及计算机技术的发展, CFD技术已成为解决各种流体运动和传热，以及场问题的强有⼒、有效的⼯具，⼴泛应⽤于⽔利、⽔电，航运，海洋，冶⾦，化⼯，建筑，环境，航空航天及流体机械与流体⼯程等科学领域。
利⽤数值计算模拟的⽅法对流体机械的部流动进⾏全三维整机流场模拟，进⽽进⾏性能预测的⽅法越来越⼴泛地被从事流体机械及产品性能取决于各种场特性的设计、科研等科技⼈员所使⽤；过去只有通过实验才能获得的某些结果或结论，现在完全可借助CFD模拟的⼿段来准确地获取。
这不仅既可以节省实验资源,还可以显⽰从实验中不能得到的许多场特性的细节信息。
⼀、什么是计算流体动⼒学计算流体动⼒学(Computational Fluid Dynamics，简称CFD)是通过计算机数值计算和图像显⽰，对包含流体流动和有热传导等相关物理现象的系统所做的分析。
CFD的基本思想可以归结为：把原来在时间域及空间域上连续的物理场（如速度场和压⼒场，以及热⼒场等），⽤⼀系列有限个离散点上变量值的集合来代替；并通过⼀定的原则和规律建⽴起关于这些离散点上的场变量之间关系，从⽽组成这些场变量之间关系的代数⽅程组；然后求解这种代数⽅程组，来获得这些场变量的近似值[1-3]；这就是流动的数值计算。
或者直观地说，通过数值计算中的各种离散⽅法，把描述连续流体运动的控制偏微分⽅程离散成代数⽅程组，由此建⽴该流动的数值模型；再根据问题的具体情况，设定边界条件和初始条件封闭⽅程组；然后通过计算机数值计算求解这种代数⽅程组，从⽽获得描述该流场场变量的某些运动参数的数值解。
(完整word版)湍流模型理论
湍流模型理论§3.1 引言自然界中的实际流动绝大部分是三维的湍流流动，如河流,血液流动等。
湍流是流体粘性运动最复杂的形式，湍流流动的核心特征是其在物理上近乎于无穷多的尺度和数学上强烈的非线性，这使得人们无论是通过理论分析、实验研究还是计算机模拟来彻底认识湍流都非常困难。
回顾计算流体力学的发展,特别是活跃的80年代，不仅提出和发展了一大批高精度、高分辨率的计算格式,从主控方程看相当成功地解决了Euler方程的数值模拟,可以说Euler方程数值模拟方法的精度已接近于它有效使用范围的极限；同时还发展了一大批有效的网格生成技术及相应的软件，具体实现了工程计算所需要的复杂外形的计算网格；且随着计算机的发展，无论从计算时间还是从计算费用考虑，Euler方程都已能适用于各种实践所需。
在此基础上，80年代还进行了求解可压缩雷诺平均方程及其三维定态粘流流动的模拟。
90年代又开始一个非定常粘流流场模拟的新局面，这里所说的粘流流场具有高雷诺数、非定常、不稳定、剧烈分离流动的特点，显然需要继续探求更高精度的计算方法和更实用可靠的网格生成技术.但更为重要的关键性的决策将是，研究湍流机理,建立相应的模式，并进行适当的模拟仍是解决湍流问题的重要途径。
要反映湍流流场的真实情况,目前数值模拟主要有三种方法：1。
平均N-S方程的求解，2。
大涡模拟（LES），3。
直接数值模拟（DNS)。
但是由于叶轮机械内部结构的复杂性以及目前计算机运算速度较慢，大涡模拟和直接数值模拟还很少用于叶轮机械内部湍流场的计算,更多的是通过求解平均N-S方程来进行数值模拟。
因为平均N-S方程的不封闭性，人们引入了湍流模型来封闭方程组,所以模拟结果的好坏很大程度上取决于湍流模型的准确度。
自70年代以来，湍流模型的研究发展迅速，建立了一系列的零方程、一方程、两方程模型和二阶矩模型，已经能够十分成功的模拟边界层和剪切层流动。
但是,对于复杂的工业流动，比如航空发动机中的压气机动静叶相互干扰问题，大曲率绕流，激波与边界层相互干扰，流动分离，高速旋转以及其他一些原因，常常会改变湍流的结构，使那些能够预测简单流动的湍流模型失效，所以完善现有湍流模型和寻找新的湍流模型在实际工作中显得尤为重要。
高精度差分格式及湍流数值模拟(三)
~m; (~10m)
最小尺度涡： 微米量级； 中小尺度涡： 0.1mm~ 1mm量级
直接数值模拟： 分辨出最小尺度涡； 网格量极为巨大： 计算量、存储量大
对高性能计算需求强烈
二、 高精度CFD软件 OpenCFD
OpenCFD: 作者开发的一套高精度、开放的CFD程序
1） OpenCFD-SC : 高精度差分 2） OpenCFD-EC: 有限体积 3） OpenCFD-Comb: 化学反应计算
computational region
200
z 300
outflow boundary
buffer region
400
500
示例2：有攻角小头钝锥边界层转捩的DNS
3） 计算结果验证
a.不同网格之间的比较 （网格收敛性） b. 不同扰动形式的结果比较 （模型正确性） d. 与Horvath静风洞实验比较 （与实验结果比较） c. 与 eN方法比较 （与理论结果比较） e. 与Stenson实验比较 （与实验结果比较）
面向工程计算开放的CFD代码(Open CFD code for Engineering Computing)
算法： 多块结构网格有限体积 + RANS 特点： 适用于复杂外形工程计算
差分-有限体积混合方法
3） 高精度化学反应模拟程序—— OpenCFD-Comb
面向化学反应的开放CFD程序： An OpenCFD code for Combustion
高精度差分格式及 湍流数值模拟 （三）
Part 3 可压缩湍流DNS
1. 背景 2. OpenCFD软件简介 3. 典型可压缩湍流的直接数值模拟示例 4. 湍流模型的评估及改进 5. 湍流燃烧的DNS 6. 小结
可压缩平板湍流边界层拟序结构的直接数值模拟
・３ ・
可 压 缩 平 板 湍 流 边 界层 拟序 结构 的直 接 数 值 模 拟
Ａ Ｄｉ ｒ ｅ ｃ ｔ Ｎｕ ｍｅ ｒ ｉ ｃ ａ ｌ Ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ｏ ｆ ａ Ｃｏ ｍｐ ｒ ｅ ｓ ｓ ｉ ｂ ｌ ｅ Ｆｌ ａ ｔ － Ｐｌ ａ ｔ ｅ Ｔｕ ｒ ｂ ｕ ｌ ｅ ｎ ｔ Ｂｏ ｕ ｎ ｄ ａ ｒ ｙ Ｌａ ｙ ｅ ｒ Ｃｏ ｈ ｅ ｒ ｅ ｎ ｔ Ｓ ｔ ｒ ｕ ｃ ｔ ｕ ｒ ｅ ｓ
着上抛和 下扫现 象， 这是壁湍能生成和耗散 的关键 。
Ａｂ ｓ ｔ ｒ ａ ｃ ｔ ：Ａ Ｄ ｉ ｒ ｅ ｃ ｔ ｎ ｕ ｍｅ ｉ ｒ ｃ ａ ｌ ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ｏ ｆ ａ ｓ ｐ ａ ｔ ｉ ａ ｌ ｌ ｙ ｅ ｖ ｏ ｌ ｖ ｉ ｎ ｇ ｃ ｏ ｍｐ ｒ ｅ ｓ ｓ ｉ ｂ ｌ ｅ ｆ ｌ ａ ｔ — ｐ ｌ ａ ｔ ｅ ｔ ｕ ｒ ｂ ｕ ｌ ｅ ｎ ｔ ｂ ｏ ｕ ｎ ｄ ａ ｒ ｙ ｌ ａ ｙ ｅ ｒ ｗｉ ｔ ｈ ｉ ｎ ｃ ｏ ｍｉ ｎ ｇ Ｍａ ｃ ｈ
ｎ ｕ ｍｂ ｅ ｒ（ Ｍ∞） ｏ ｆ ２ ． ２ ５ ｗ ａ ｓ ｐ ｅ ｆ ｒ ｏ ｒ ｍ ｅ ｄ ｉ ｎ ｔ ｈ ｉ ｓ ｐ ａ ｐ ｅ ｒ ． Ｔ ｈ ｅ ｃ ｏ ｈ ｅ ｒ ｅ ｎ ｔ ｓ ｔ ｒ ｕ ｃ ｔ ｕ ｒ ｅ ｆ ｏ ｒ ｍｉ ｎ ｇ ｐ ｒ ｏ ｃ ｅ ｓ ｓ ｉ ｓ ｓ ｔ ｕ ｄ ｉ ｅ ｄ ． Ｔ ｈ ｅ ｌ ｏ ｗ － ｓ ｐ ｅ ｅ ｄ ｓ ｔ ｉ ｒ ｐ ａ ｐ ｐ ｅ ａ ｒ ｉ ｎ ｇ ｉ ｎ ｔ ｈ ｅ
ｅ ｎ ｅ ｒ ｇ ｙ ｇ ｅ ｎ ｅ ｒ ａ ｔ ｉ ｏ ｎ ａ ｎ ｄ ｄ ｉ ｓ ｓ ｉ ｐ ａ ｔ ｉ ｏ ｎ ．
湍流的几种数值模拟方法
LES特点
抓大不放小 非常有利，有力的工具 是最近,可预见未来流体 力学研究和应用的热点 近来又出现了VLES, DES等在LES上发展而 来的工具
Will RANS survive LES? Hanjalic自问自答
会。Journal of Fluids Engineering -V127, 5, pp. 831-839 (Will RANS
Prandtl(1925)混合长度模型
也被称作零方程模型 还在被广泛应用 廉价，易收敛 基本在流场比较简单，或者对计算结果 精度要求不高或者流场形状比较复杂的 行业中，比如暖通空调，流体机械等。
Prandtl混合长度模型 缺点
最明显的缺点是:当速度梯度 为零的 时候， 消失， 这与事实不符
Launder and Li(1994), Craft and Launder (1995)
目前有很多学者在继续此方面的工作
Brian E. Launder
本科Imperial College, London 硕博 MIT 实验流体力学 1964-1976 Imperial College 讲师
涡流粘度
Eddy viscosity or turbulent viscosity
二维流场分子粘性力
为描述雷诺应力，Boussinesq 1887 定义了与之相对应的
RANS模型的核心在于给出 的数 学表达式，要求精度高，适用范围广
涡流粘度,
Prandtl 1925 Prandtl 1945 Bradshaw 1968 Kolmogorov, 1942 Hanjalic 1970 Rotta 1951 Chou 1945 Davidov 1961
湍流流动的可压缩性与不可压缩性研究
湍流流动的可压缩性与不可压缩性研究引言湍流是自然界和工程领域中普遍存在的一种现象，具有复杂、非线性和随机的特点。
在流体力学中，湍流流动的研究一直是一个重要的课题。
湍流流动可以分为可压缩性和不可压缩性两种情况，对湍流的可压缩性和不可压缩性研究对于理解和掌握湍流现象具有重要意义。
可压缩性湍流流动的研究可压缩性湍流流动是湍流中存在压力变化的一种情况。
在高速气流和超音速气流中，湍流流动通常是可压缩的。
可压缩性湍流流动的研究对于航空航天工程、气象学和天体物理学等领域具有重要应用价值。
可压缩性湍流流动的特点可压缩性湍流流动具有以下特点：1.压缩性效应：气体在可压缩性湍流流动中经历压力和密度的变化，相比于不可压缩性湍流流动，具有更强的压缩性效应。
2.高速性：可压缩性湍流流动通常发生在高速气流和超音速气流中，流动速度较大。
3.非线性性：可压缩性湍流流动是非线性的，无法用简单的线性模型描述。
4.随机性：可压缩性湍流流动具有随机性，其流动结构在空间和时间上都是不规则和随机的。
可压缩性湍流流动的数值模拟可压缩性湍流流动的数值模拟是研究该领域的重要方法之一。
通过数值模拟，可以对湍流流动的可压缩性进行分析和预测，提供工程设计和科学研究的依据。
数值模拟可压缩性湍流流动的方法主要包括直接数值模拟（DNS）和雷诺平均Navier-Stokes方程(RANS)求解。
DNS方法适用于小尺度和短时间尺度的流动问题，可以提供精确的流场信息，但计算量较大。
RANS方法通过对流场进行时间平均和空间平均，提供了湍流的统计特性信息，计算效率较高。
数值模拟中的湍流模型是可压缩性湍流流动研究的关键。
目前常用的湍流模型有标准k-ε模型、雷诺应力模型和Reynolds-Stress模型等。
这些模型基于不同的假设和参数化方法，对湍流流动进行描述。
可压缩性湍流流动的实验研究除了数值模拟，实验研究也是可压缩性湍流流动研究的重要手段之一。
实验室和风洞中的试验可以模拟可压缩性湍流流动的场景，通过对流场的测量和分析，获取流动的统计特性。
大型客机复杂可压缩流的大涡模拟主要研究方法
大型客机复杂可压缩流的大涡模拟主要研究方法一、大涡模拟基础1. 大涡模拟简介大涡模拟是一种将流场分解成小尺度湍流和大尺度湍流的方法。
在LES中，大尺度结构通过直接数值模拟来求解，而小尺度结构则通过子网格模型（sub-grid model）进行建模。
由于小尺度结构不再需要直接求解，因此可以使用更粗的网格来进行计算，从而减少计算量。
同时，LES还能够提供更加真实的湍流统计数据，如湍流强度、湍流长度等。
2. LES的优点和局限性与其他流体力学方法相比，LES有以下几个优点：（1）能够考虑湍流中的时间和空间尺度差异，提供更加真实的湍流信息；（2）计算结果对于网格的依赖性相对较小，使得计算可以在较粗的网格上进行；（3）LES能够模拟复杂流场，如湍流燃烧、多相流等。
虽然LES具有很多优点，但它也有一些局限性：（1）计算量较大，需要使用高性能计算机进行计算；（2）由于需要建立子网格模型，LES的结果可能受到模型误差的影响；（3）由于直接数值模拟只考虑了大尺度结构，因此对于小尺度结构的预测可能存在误差。
二、大涡模拟在大型客机流场研究中的应用1. 大涡模拟在飞行器气动力学研究中的应用大型客机的外形复杂，流场也非常复杂。
对于这样的流场，传统的计算流体力学方法可能无法准确地预测气动力学行为。
因此，大涡模拟成为研究大型客机流场的一种重要方法。
在大涡模拟中，通过将流场分解成大尺度结构和小尺度结构，可以更加准确地模拟大型客机流场中的湍流现象。
大涡模拟还能够提供更加真实的气动力学数据，如升阻比、气动力矩等。
这些数据对于飞机设计和优化非常重要。
2. 大涡模拟在飞行器噪声研究中的应用随着人们对噪声污染的关注度不断提高，飞机噪声研究也越来越受到关注。
大型客机飞行时产生的噪声主要来自于引擎和机翼表面的湍流。
由于湍流现象非常复杂，传统的计算流体力学方法无法准确地预测噪声的产生和传播。
因此，大涡模拟成为研究飞机噪声的一种重要方法。
通过大涡模拟，可以更加准确地模拟湍流现象，从而预测噪声的产生和传播方式。
湍流的数值模拟方法进展
3 大涡模拟（LES ）湍流大涡数值模拟（LES ）是有别于直接数值模拟和雷诺平均模式的一种数值模拟手段.利用次网格尺度模型模拟小尺度湍流运动对大尺度湍流运动的影响即直接数值模拟大尺度湍流运动, 将N —S 方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。
3。
1 基本思想很多尺度不同的旋涡一起组成了湍流运动平均流动主要取决于大漩涡的流动，大尺度运动则受到小旋涡的影响。
流动中的大涡实现了动量、能量质量、热量的交换，耗散主要是由于小涡作用的。
大旋涡中受到流场形状、阻碍物的影响，,使大漩涡的各向异性更加明显。
然而小漩涡之间各项同性，相互没有太大的区别，所以建立统一的模型比较容易一些.综上所述,大涡模拟将湍流瞬时运动量通过滤波将运动分成小尺度和大尺度.大尺度的运动受到小尺度的运动的影响可以通过应力项（类似于雷诺应力项）来表示，即为亚格子雷诺应力，以建立这种模型的方法来模拟。
而大尺度则是求解运动微分方程而计算出来的，也就是说大涡模拟,要先过滤掉小尺度的脉动,然后再推出小尺度的运动封闭方程以及大尺度的运动控制方程。
3。
2 滤波函数正如上面提到，大涡模拟要先将流动变量分解成小尺度量和大尺度量，我们把这个作用叫做滤波.滤波运算就是在一区域内按照一定的条件对函数进行加权平均，作用是将高波数滤掉，使低波数保留，滤波函数的特征尺度决定了截断波数的最大波长,下面三种滤波函数是最为常用的主要有以下三种：盒式、富氏截断以及高斯滤波函数.不可压常粘性系数的湍流运动控制方程为N-S 方程:j ij i j j i i x S x P x u u t u ∂⋅∂+∂∂-=∂∂+∂∂)2(1γρ式中:S 拉伸率张量，表达式为：2/)//(i j j i ij x u x u S ∂∂+∂∂=；γ分子粘性系数；ρ流体密度。
设将变量i u 分解为方程（11）中i u 和次网格变量(模化变量）'i u ,即'+=i i i u u u ，i u 可以采用Leonard 提出的算式表示为:(11）式中)(x x G '-称为过滤函数，显然G(x)满足x d x u x x G x u i i '''-=⎰+∞∞-)()()(⎰+∞∞-=1)(dx x G3.3 控制方程将过滤函数作用与N —S 方程的各项，得到过滤后的湍流控制方程组:由于无法同时求解出变量i u 和j i u u ,所以将j i u u 分解成i j i j ij u u u u τ=⋅+，ij τ即称为次网格剪切应力张量(亦称为亚格子应力)。
可压缩全速湍流流动的边界数值处理方法研究
对每一控制方程分别 阐述特定的数 值处理方法。应用提 出的边界处理 方法对单 圆弧 凸包通 道进行数值仿真获 得 了合理 的 结果 ， 音速和超音速情形下准确地计算 出了流场中存 在的激波。 跨
（ ．西 北 工 业 大 学 ， 西 西 安 ７ ０７ ２ １ 陕 １０ ２；．中国 燃 气 涡 轮研 究 院 ， 四川 成 都 ６０ ０ ） ｌ５ ０
摘要 ： 计算流体力学问题 的边界条件处理方法关 系到数值仿真结果 的精确度 。为解 决算法的精度 ， 出了三维可压缩 湍流 提
流动的边界条件数值处理方法 ， 对所研究 的边界类型包括进 Ｉ边界 、 口边界 和固体 壁面 ， ： １ 出 流动 的速度 范围涉及亚音 速、 跨
ａ ｄ ｗａｌ ｏ ｎ ａｙ ａ ｅ ｉｃ ｕ ｅ ｎ ｌ ｂ ｕ ｄ ｒ ｒ ｎ ｌ ｄ ｄ，ａ ｄ ｔ ｅ ｆ ｗ ｐ ｅ ｏ ｅ ｓｓ ｂ ｏ ｉ ｔａ ｓ ｎ ｃ ａ ｄ ｓ ｐ ｒ ｏ ｉ ．Ｓ ＭＰ Ｅ ａｇ ｒ ｈ ｉ ｎ ｈ ｏ ｓ ｅ ｄ ｃ ｖ ｒ ｕ ｓ ｎｃ， ｒ ｎ ｏ ｉ ｎ ｕ ｅ ｓ ｎ ｃ Ｉ Ｌ ｏ ｉ ｍ ｌ ｌ ｔ ｓ
关键词 ： 界条件 ； 边 可压缩湍流； 全速流动 ； 压力耦合算法
中图 分 类 号 ：３７． ０５ １ 文献 标 识 码 ： Ａ
Ｎｕｍ ｅ ｉａ ｅ ｈ ｄ ｆＢｏ ｎｄ ｒ ｅ ｔ ｅ ｔｆ ｒ Ｃｏ ｐ ｅ ｓｂｌ ｒｃ ｌＭ ｔ ｏ ｓ ｏ ｕ ａ ｙ Ｔｒ ａ ｍ ｎ ｏ ｍ ｒ ｓ ｉ ｅ Ｔｕｒ ｕｌｎｔＦｌｗ ｔＡ１ Ｓ ｅ ｓ ｂ ｅ ｏ ａ １ ｐｅ ｄ
不同湍流模型在管道流动数值模拟中的适用性研究
不同湍流模型在管道流动数值模拟中的适用性研究邵杰;李晓花;郭振江;刘瑞璟;田晓亮【摘要】Currently numerical simulation has been applied in thefields of scientific research and engineering in large scale. Turbulent model is often used in simulation. But different turbulent model has its applicable scope respectively. In this article, by using some common turbulent models provided in CFD software FLUENT, the numerical simulation of turbulentflow in pipe was carried out and the frictional drag resulted from simulation was compared with that obtained in experiment. It was shown from the results of analysis that Spalart-Allmaras model,k-ε (EWT) model and Reynolds stress (EWT) model are suitable for hydraulically smooth pipe with laminarflow, butk-ε model is suitable both of laminar and turbulentflows; for hydraulically smooth pipe with laminarflow, the highest precision can be reached by use of Spalart-Allmaras model; for coarse surface pipe with laminarflow, coarse degree should be adjusted in use ofk-ε model.%针对数值模拟在科学研究和工程实践领域中的大规模应用，湍流模型是数值模拟中常用的模型，不同湍流模型有自己的适用范围。
湍流直接数值模拟
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
混淆误差产生原因： 当用伪谱方法计算时，谱空分量用，它等于：
利用三角级数公式
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
在上式中，p、q、k的取值范围都是{一N/2，N/2-1}，所 以取值不等于零的波数组合情况只有两种：p+q-k＝0，或 p+q-k＝ N，于是有
Kolmogorov耗散尺度
动速度均方根值)，将以上关系代入式，可得：
1/4 3
，而 ~ u '3 / l （ u ' 是脉
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率
那么三维总网格数N：
＝10４
这是一个天文数字的估计 ，假设 Rel ＝104，就要求网格数为109， 考虑到计算的流动变量数，需要约1010字长的计算机内存。直接数值 模拟实际工程湍流运动时，对网格分辨率的要求更高。 应当指出选定最小的网格长度还和数值方法有关。谱方法的数 值精度最高，差分法的精度和差分格式有关。
边界条件的提法： (3)渐近条件 对于湍流边界层或其他薄湍流切变层，在远离薄层和物面的渐近区 域，速度场趋近于无旋的均匀场，因此对于不可压缩流体可以采用如下：
刚盖假定：数值方法只能计算有限域内的流动，渐近条件只能采用 近似形式，一种方法是在离开薄层或物体横向一定距离的平面上设 置“虚拟边界”，在虚拟边界y＝H满足以下条件，称为刚盖假定：
式中，第一项是乘积w＝uv在谱空间中的投影，第二项是 在伪谱运算中产生的误差，故称混淆误差。
应用GAO-YONG可压缩湍流模式数值模拟RAE2822翼型绕流
粘性绕流 问题模拟 的状况 看 ， 有几个 问题 没有很好 解决 ： 由于 激波 引起 的逆压 梯 度导 致 的分离 流 动 ； 由 ① ② 于对各项 异性 的湍 流粘性 考虑不 充分 ， 导致 转捩 点 、 波位 置 的捕获不 够 理想 ； 激 ③激 波湍 流边 界 层干 扰 中的
非定 常性 和热载荷 很难 准确模 拟 ．
应 用 ＧＡ Ｙ Ｇ 可 压 缩 湍 流 模 式 Ｏ— ＯＮ 数 值 模 拟 Ｒ Ｅ ８２翼 型 绕 流 Ａ ２２
闫文辉 ， 闫 巍 ， 高 歌
（ ． 京 航 空航 天大 学 能 源 与动 力 学 院 热 动 力 重点 实 验 室 ，北京 １北 １０ ８ ） ００３
［ 摘 要 ］ 应 用 Ｇ ｏ ｏｇ可压 缩 湍 流 模 式 ，数 值 模 拟 Ｒ Ｅ ８２二 维 翼 型在 两 种 不 同来 流 情况 下 的跨 音 速 粘 性 绕 ａ— ｎ Ｙ Ａ ２２
本文基 于 Ｇ Ｏ Ｙ Ｎ Ａ ． Ｏ Ｇ可压缩湍 流模式 ， 控 制方 程 时 间离散 采用 五 步 Ｒ ｎｅＫ ｔ 对 ｕ ｇ． ｕｔ ａ显示 时 间推进 、 空间 离散采用 单元 中心型有 限体积法 ， 无粘 通 量采 用 失 通量 差分 分 裂 Ｒ ｅ 式 ， ｏ格 粘性 通 量 采用 二 阶 中心差 分 格
湍流的理论与实验研究
湍流的理论与实验研究湍流的理论与实验研究湍流是流体力学界公认的难题，被认为是经典物理学中最后一个未被解决的问题。
自然界和工程领域的绝大多数流动都是湍流，因此湍流研究具有重大意义。
近年来，随着实验测量技术和数值模拟能力的不断增强，学术界对高雷诺数和高马赫数湍流有了许多新的认识。
我国科学界也结合国家重大战略需求和学科发展前沿，分析国际上湍流研究的特点、现状和发展趋势，希望对湍流产生机制和流动本质进行深入研讨，加强与航空、航天、航海等相关单位和部门间的沟通与联系，推动湍流研究的发展。
针对国内学科发展现状，尤其是实验研究相对薄弱的特点，国家自然科学基金委员会数理科学部、工程与材料科学部和政策局，于2014年3月20-21日在北京联合举办了第110期双清论坛，论坛主题为“湍流的理论与实验研究”。
来自全国15个单位的近50位流体力学与工程领域的专家学者应邀出席。
与会专家通过充分而深入的研讨，凝练了该领域的重大关键科学问题，探讨了前沿研究方向和科学基金资助战略。
本期特刊登此次论坛学术综述。
一、湍流研究的重要意义自1883年雷诺（Reynolds）发现湍流以来，湍流问题的研究一直困扰着众多学者。
著名物理学家费曼曾说，湍流是经典物理学中最后一个未被解决的难题；2005年《科学》杂志在其创刊125周年公布的125个最具挑战性的科学问题中，其中至少两个问题与湍流相关。
在我们日常生活中，湍流无处不在。
自然界和工程应用中遇到的流动，绝大部分是复杂的湍流问题。
在自然界，从宇宙星系的时空演化，到星球内部的翻滚流动，从大气环流的全球运动，到江河湖泊的区域流动，都有湍流的身影。
在工程领域，从陆地、海洋、空天等交通运载工具，到原子弹、氢弹、导弹、战斗机、舰船等国防武器的设计；从全球气象气候的预报，到地区水利工程的设计；从传统行业如叶轮机械、房桥建筑、油气管道，到新兴行业如能源化工、医疗器械、纳米器件的设计，都需要了解和利用湍流。
因此，湍流流动的研究不仅仅是一个学科发展的问题，更具有重要的工程应用价值。
蒸汽喷射器三维流场的数值模拟计算与分析
大连理工大学硕士学位论文蒸汽喷射器三维流场的数值模拟计算与分析姓名：***申请学位级别：硕士专业：热能工程指导教师：李素芬;沈胜强20000601摘要ｒ气体喷射器作为一种节能装置，可回收大量余热，起到了节能和环保的双重作用，在工业部门中得到广泛应用。
其内部经历着复杂的多维湍流流动过程，而其中喷嘴更是决定喷射器是否正常工作的关键部件。
ｊ本文在详细分析喷射器内部流动的基础上，建立了三维湍流流动的数值模拟计算模型，并主要对喷嘴的流场进行了详细的计算分析。
本文主要内容有：１、深入分析了ＫＩＶＡ系列程序与相关的ＣＦＤ理论方法，结合气体喷射器喷嘴的流动特点，建立了喷射器喷嘴复杂流场结构的三维数值模拟计算模型和计算方法，并应用于喷射器喷嘴稳态流场的数值模拟计算中。
２、根据气体喷射器结构和特点建立了喷射器整体及喷嘴通用计算网格的生成方法，并编制了相应的计算网格生成程序。
其网格生成方法及程序适用于各种结构及尺寸的喷嘴和喷射器，充分体现了其灵活性和实用性。
３、运用本文开发的通用计算网格生成程序结合三维流场数值模拟计算程序，针对不同的边界条件和结构尺寸的喷嘴流场，进行了数值模拟计算，考察了以上各特性参数对喷射器内部流动的影响，并根据计算结构的分析提出了喷射器喷嘴设计的建议。
４、比较全面地考虑了各种不可逆因素（如摩擦、散热等）对流场各参数的影响，进一步完善了喷射器的研究■一一关键词：喷嘴、数值模拟、流场ＡＢＳＴＲＡＣＴＡｓａｋｉｎｄｏｆｄｅｖｉｃｅ，ｔｈｅｓｔｅａｍｅｊｅｃｔｏｒｃａｒｌｒｅｃｙｃｌｅａｇｒｅａｔｄｅａｌｏｆｅｎｅｒｇｙ，ａｎｄａｔｔｈｅｓａｍｅｔｉｍｅ，ｉｔｐｌａｙａｇｒｅａｔｒｏｌｅｏｆｅｎｖｉｒｏｎｍｅｎｔｐｒｏｔｅｃｔｉｏｎ，ＳＯｉｔｉｓａｐｐｌｉｅｄｉｎｍａｎｙｉｎｄｕｓｔｒｙｄｅｐａｒｔｍｅｎｔｓ．Ｉｔｓｆｌｏｗｆｉｅｌｄｉｓｍｕｌｔｉ—ｄｉｍｅｎｓｉｏｎｓ，ｔｒａｎｓｉｅｎｔ，ｔｕｒｂｕｌｅｎｔ，ｓｕｂｓｏｎｉｃａｎｄｓｕｐｅｒｓｏｎｉｃｆｌｏｗｓ．Ａｎｄｔｈｅｎｏｚｚｌｅｉｓｔｈｅｋｅｙｏｆｔｈｅｅｊｅｃｔｏｒ．Ｏｎｔｈｅｂａｓｅｏｆｅｘｐａｔｉａｔｉｎｇｏｎｔｈｅｆｌｏｗｓｉｎｓｉｄｅｔｈｅｓｔｅａｍｅｊｅｃｔｏｒ，ａｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａｌ，ｔｕｒｂｕｌｅｎｔ，ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｃｏｍｐｕｔａｔｉｏｎａｌｍｏｄｅ】．ａｎｄｔｈｅｍｅｔｈｏｄｉＳｕｔｉｌｉｚｅｄｅｍｐｈａｔｉｃａｌｌｙｏｎａｎａｌｙｓｉｓａｎｄｃａｌｃｕｌａｔｉｏｎｔｈｅｆｌｏｗｆｉｅｌｄｏｆｔｈｅｎｏｚｚｌｅ．Ｔｈｅｍａｉｎｗｏｒｋｓａｒｅｓｕｍｍａｒｉｚｅｄａｓｆｏｌｌｏｗｓ：１．ＡｎａｌｙｚｅｔｈｅｐｒｉｎｃｉｐｌｅｏｆｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｏｆｔｈｅＫＩＶＡ一３ｃｏｄｅａｎｄＣＦＤｍｅａｎｓ．ｃｏｍｂｉｎｉｎｇｔｈｅｆｌｏｗｉｎｇｃｈａｒａｃｔｅｒｉｓｔｉｃｏｆｎｏｚｚｌｅ，ａｐｒｏｇｒａｍｓｕｉｔａｂｌｅｔｏｃｏｍｐｕｔｅｔｈｉｓｋｉｎｄｏｆｆｌｏｗｆｉｅｌｄｂｙｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｍｅｔｈｏｄｉｓｃｏｍｐｉｌｅｄ．２．Ａｐｐｌｙｔｈｅｍｅｔｈｏｄｏｆｂｏｄｙ·ｆｉｔｔｅｄｍｅｓｈｇｅｎｅｒａｔｉｏｎａｎｄｔｈｅｂｌｏｃｋ—ｓｔｒｕｃｔｕｒｅｄｍｅｔｈｏｄ，ａｃｏｍｍｏｎｐｒｏｇｒａｍｉｓｃｏｍｐｉｌｅｄ．Ｉｔｃａｎｂｅｎｏｔｏｎｌｙｕｔｉｌｉｚｅｄｏｎｔｈｅｅｊｅｃｔｏｒ，ｂｕｔｍａｎｙｃｏｍｐｌｉｃａｔｅｓｔｒｕｃｔｕｒｅｆｌｏｗｆｉｅｌｄｓ．３．Ｍｏｂｉｌｉｚｉｎｇｔｈｅｃｕｒｒｅｎｔｇｒｉｄｄｉｎｇｐｒｏｇｒａｍａｎｄｔｈｅｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｃｏｍｐｕｔａｔｉｏｎａｌｍｏｄｅｌ，ａｎａｌｙｚｅａｎｄｃａｌｃｕｌａｔｅｔｈｅｆｌＯＷｆｉｅｌｄＯｆｔｈｅｎｏｚｚｌｅ，ｄｉｓｃｎｓｓｔｈｅｅｆｆｅＣｔｓｏｎｔｈｅｆｌＯＷ０ｆｖａｒｊ０ＵＳｂｏｕｎｄａｒＹＣＯｎｄｉｔｉｏｎｓ，ｓｔｒｕｃｔｕｒｅＳｉｚｅ．Ｔｈｅｒｅｓｕｌｔｓｐｒｅｓｅｎｔｐａｒｔｉｃｕｌａｒｓｕｇｇｅｓｔｉｏｎｆｏｒｔｈｅｏｐｔｉｍｉｚｉｎｇｄｅｓｉｇｎｏｆｔｈｅｎｏｚｚｌｅ．４．ＧｅｎｅｒａｌｌｙｃＯｎｓｉｄｅｒｔｈｅｉｎｆｅＣｔｉＯｎＳｏｆｍａｎＹｋｉｎｄＳＯｆｕｎｒｅｖｅｒｓｉｂｌｅｆａｃｔｏｒｓ（ｆｒｉｃｔｉｏｎ，ｈｅａｔｄｉｓｐｅｒｓｉｏｎ），ａｎｄｍａｋｅｔｈｅｒｅｓｅａｒｃｈｏｆｎｏｚｚｌｅｏｒｅｊｅｃｔｏｒｍｏｒｅｐｅｒｆｅｃｔ．Ｋｅｙｗｏｒｄｓ：ｎｏｚｚｌｅ，ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ，ｆｌｏｗｆｉｅｌｄ第一章绪论第一章绪论本章在查阅大－￥－文献的基础上．ｘｅａ－喷射器及：ｇ－数值－｝－Ｉ－算等研究领域的发展和概况进行了详细的综述，并概括出本文的主要内容。
CFD软件及数值模拟湍流理论
2.2.4能量利用系数在对气流进行综合评价的时候需要用到能量利用系数，利用该系数可以解决很多在能耗方面的问题。
在空调系统中存在室内气流的影响，在这种影响下可以进行气流研究。
在对气流进行研究的时候需要对空调的工作区等部分进行调查，需要使得空调进行合理的送风而不是将空调预热带向错误的位置，通过这种方式可以提高空调排热效率，进而减少总系统的功耗，使得空调更具有节能减排的经济性。
在对其进行研究的过程中用到了能量利用系数，将其用η代表，那么其定义式如下[8]；p o n o t t t t η-=-式（2-10）式中 o t ：送风温度，℃，p t ：排风温度，℃，n t ：室内工作区的平均温度，℃， 能量利用系数在一般情况下是由默认值的，在混合通风系统中一般默认为1.0，但是在有些情况下如下送风将会大于1.0。
事实上，能量利用系数显示的是室内的热力分层特性，即室内温度的梯度变化。
对于通风系统或者全新风系统，节能潜力随着能量利用系数的增大而升高。
有时，能量利用系数越高，表明室内温度的梯度越大，所以不能够太过分地强调能量利用系数，可能会影响到舒适性。
3 CFD软件及数值模拟湍流理论气流组织的模拟运用数值计算的方法对大空间速度场、温度场的分布规律进行研究,数值方法是模拟的基础,它对实现气流组织的模拟有着重要的意义。
因此下面对数值模拟的相关理论作以详细介绍[27,28]。
3.1 CFD软件简介CFD是英文Computational Fluid Dynamics（计算流体力学）的简称[29]，其伴随数值计算及计算机技术的发展而发展。
通俗地讲，CFD是一种虚拟的实验，他通过使用计算机科技技术来虚拟实验数据，通过这些数据来模拟我们所需要的流动状况，从而进行研究。
这样的话可以利用网络技术减少很多不必要的繁琐工作，他利用了相关的数学微分公式，通过这种方程式的技术来进行近似模拟。
其基本结构包括三大模块，即前处理、求解器和后处理，每个模块都有其独特的作用。
可压缩平板湍流边界层拟序结构的直接数值模拟
可压缩平板湍流边界层拟序结构的直接数值模拟高慧【摘要】文中直接数值模拟了来流马赫数(M肄)为2.25的空间发展的可压缩平板湍流边界层.研究了可压缩平板湍流边界层内的拟序结构的形成过程.指出在粘性底层出现低速条带是触发拟序结构的第一个信息.条带升起，开始发生振动，然后破裂，并伴随着产生强烈的湍流脉动，这便形成了湍流的猝发过程.因此壁湍流拟序结构的重要意义在于它是生成湍流的重要机制.条带的破裂还伴随着上抛和下扫现象，这是壁湍能生成和耗散的关键.%A Direct numerical simulation of a spatially evolving compressible flat-plate turbulent boundary layer with incoming Mach number (M∞) of 2.25 was performed in this paper. The coherent structure forming process is studied. The low-speed strip appearing in the viscous sublayer is the first information to trig coherent structure. Strip rises, begin to vibrate, and rupture, and accompanied by strong turbulence, which will form the turbulent bursting process. Therefore the wall turbulent coherent structure is significant in that it is the important mechanism to generate turbulence. Strip rupture accompanied by up and down sweep phenomenon, it is the key of wall turbulent energy generation and dissipation.【期刊名称】《价值工程》【年(卷),期】2013(000)010【总页数】2页(P3-4)【关键词】可压缩;平板湍流边界层;直接数值模拟;拟序结构【作者】高慧【作者单位】中国石油大学石油工程学院，青岛266555【正文语种】中文【中图分类】O350 引言对湍流拟序结构的研究，使得人们在揭示湍流的物理本质方面寻觅到了一个新的有效的方向，从而得以有能力预测和控制湍流。
湍流氢氧预混燃烧流场的二维数值模拟
燃烧 室 流场 中湍 流与 燃烧 两种 物理 现象 的强 烈 耦 合， 使得 燃烧 流动 现象 异 常复 杂 ， 伴有 剧烈 的放 并
热化 学反 应 。 因此 ， 烧 室 内湍 流 燃烧 流 动是 一 燃
个非 常复 杂 的现 象 ， 目前 其 工 作 机 理 尚未 被 充 分
定 常 可 压 缩 燃 烧 流 动 的 解 算 器 Ｓｅｄ ｒａｔ ｇ ｔ ｙｅｃｉ — ａ ｎ
Ｆｏ ｍ ｏ ａ
１ ５ 计 算模 型 与计算 网格 ．
控 制 方程 中 的瞬 态项 采 用 Ｌ ．Ｇ Ｕ Ｓ Ｓ隐 式迭 代 法 求
解 ， 流项 采用 二 阶迎风 的高斯差 分格 式 ， 对 扩散 项 采用 二 阶线性 修正 的高 斯差 分格 式 。湍流模 型采 用形 式简 单 、 用 方法 的 ｋＥ模 型 ， 使 － 能得 出相 当 ห้องสมุดไป่ตู้
燃 烧 室 是 发 动机 的重 要 组成 部 分 ， 烧 室 内 燃 的燃 烧过 程 基本 上 都 是 湍 流 燃 烧 ， 由于 湍 流 自身 空 间上 的尺 度多重 性 和时 间上 的高 频脉 动性 以及
足 动力 学方 程 的湍 流瞬 时运 动分 解为平 均 运 动 和
脉 动 运动 两部 分 ， 后 把 脉 动 运 动这 部 分 对 平 均 然
湍 流流 动计 算 方 法 。Ｄ Ｓ的计 算 量 非 常大 ， Ｎ 以 消耗最 多 的计 算 资 源 来 捕 捉 流 动 中 最 细 微 的 变 化 ， 揭示 湍流 流动 现象 的 机理 方面 意义 重大 ， 在 截 至 目前 乃 至未来 相 当长 一 段 时 间 里 ， Ｎ Ｄ Ｓ方法 对 于大多 数实 际 工 程 的 应 用都 具 有 相 当 的局 限性 。
基于可压缩SIMPLE算法的叶栅通道湍流流场的数值模拟
关键词 ： 流动控 制方 程 ； 差分格式 ； 网格 分布 ； 数值模拟 分类号 ：Ｈ ５ Ｔ ４ 文献标识码 ： Ａ 文章编号 ： ０ —８ ４ ２０ ）４０ ６ －３ １ １５８ （０ ７ ０ － ６０ ０ ２
Ｎｕ ｒ ａ ｎ ｅｔ ａｉｎｏ ｕｂ ｌｎ ｌｗ ｆｒＣ ｓａ ｅＰ ｓａｅＢ ｓｄ ｍｅ ｃｌＩｖ ｓ ｇｔ ｆＴ ｒｕｅ ｔＦｏ ｏ ａｃ ｄ ａｓｇ ａｅ ｉ ｉ ｏ ｏ ｏ ｒｓｉｌ Ｉ Ｌ ｒｔｍｅｉ ｎ Ｃ ｍｐｅｓｂｅＳＭＰ Ｅ Ａ ｈ ｔ ｉ ｃ
周红梅 ， 苏莫 明， 任 雄
可压缩流体流动的建模与仿真
可压缩流体流动的建模与仿真引言流体力学是研究流体的运动规律和力学性质的学科。
在实际工程应用中，流体力学模拟与仿真对于分析和优化流体流动过程具有重要意义。
在许多情况下，需要考虑流体的可压缩性，即流体在流动过程中密度变化较大的情况。
本文将介绍可压缩流体流动的建模与仿真方法，旨在提供一个综合的理论框架和实用的工具。
1. 可压缩流体的基本概念可压缩流体与不可压缩流体在性质上存在一定差异。
可压缩流体在流动过程中，流速、压力和密度都发生变化。
同时，可压缩流体流动的数学模型也需要包括能量方程，以考虑流体的内能变化。
由于可压缩流体的复杂性，建模和仿真方法需要结合流体力学基本方程、状态方程等知识。
1.1 可压缩流体的基本方程可压缩流体的运动可以通过质量守恒、动量守恒和能量守恒方程来描述。
质量守恒方程表明流体的质量在流动过程中不会发生变化。
动量守恒方程根据牛顿第二定律描述流体的加速度与作用力之间的关系。
能量守恒方程考虑了流体的内能变化。
1.2 状态方程与压缩性状态方程描述了流体的压缩性质，它给出了流体密度与压力、温度之间的关系。
不同的流体可以采用不同的状态方程来描述其性质。
常见的状态方程有理想气体状态方程、热力学状态方程等。
2. 可压缩流体流动的建模方法可压缩流体流动的建模方法依赖于流体力学的基本方程和状态方程。
在实际工程中，常用的建模方法包括有限体积法、有限元法和高阶精确方法。
2.1 有限体积法有限体积法是一种将流体流动区域划分为有限体积单元，并在每个单元内求解基本方程的方法。
该方法通过对每个体积单元内的质量、动量和能量进行离散，建立了离散形式的基本方程组。
通过迭代求解，可以得到流场的数值解。
2.2 有限元法有限元法是一种以将流动区域离散成有限个小区域，并在每个小区域内建立适当的数学模型的方法。
该方法采用变分原理和试验函数的思想，将连续的问题转化为离散的问题。
通过求解离散形式的基本方程，得到流场的数值解。
2.3 高阶精确方法高阶精确方法是一种通过提高离散格式的阶次来提高数值解的精度和收敛性的方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
可压缩湍流化学反应流的数值模拟研究
可压缩湍流化学反应流的数值模拟研究在化学工程领域具有重要
意义，能够帮助工程师和研究人员更好地理解和优化工业过程中的化
学反应流动。
本文将介绍湍流化学反应流的数值模拟研究的基本原理、方法和应用。
在可压缩湍流化学反应流的数值模拟研究中，首先需要建立数学
模型来描述流体运动和化学反应。
通常采用的是基于Navier-Stokes
方程和化学动力学方程的耦合模型。
Navier-Stokes方程描述了流体的运动，包括惯性力、粘性力和压力力等因素；化学动力学方程则描述
了化学物质在流动中的传输和反应。
通过求解这些方程组，可以得到
流体流动和化学反应的详细信息。
数值模拟研究中，对于湍流流动的模拟通常采用雷诺平均N-S方
程（RANS）和大涡模拟（LES）等方法。
RANS方法适用于大规模湍流结构较为简单的情况，而LES方法则能够模拟湍流中的小尺度涡旋，更
适用于复杂流动的模拟。
在化学反应方面，通常采用有限元法或有限
体积法求解化学动力学方程，考虑各种化学反应过程的影响。
可压缩湍流化学反应流的数值模拟研究在石油化工、环境保护、新能源开发等领域具有广泛的应用。
例如，在石油炼制过程中，通过模拟反应流动的变化，可以提高反应器的效率和产率；在环境保护方面，可以模拟大气中的化学反应过程，预测有害物质的扩散和传输；在新能源开发中，可以模拟燃烧过程，优化能源利用效率。
总之，可压缩湍流化学反应流的数值模拟研究是一项复杂而重要的工作，需要综合考虑流体力学、化学动力学和数值方法等多方面知识。
通过不断改进模型和算法，将能够更好地应用于工程实践，推动工业化学反应流动的进步和发展。