第三章, 湍流模型 第一节, 前言 湍流流动模型很多, 但大致可以归纳为以下三类 : 第一类是湍流输运系数模型, 是 Boussnesq 于 877 年针对二维流动提出的, 将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积 即 : u ρ u u = - 推广到三维问题, 若用笛卡儿张量表示, 即有 : u u ρuu = + ρδ - 模型的任务就是给出计算湍流粘性系数 的方法 根据建立模型所需要的微分方程的数 目, 可以分为零方程模型 ( 代数方程模型 ), 单方程模型和双方程模型 第二类是抛弃了湍流输运系数的概念, 直接建立湍流应力和其它二阶关联量的输运方程 第三类是大涡模拟 前两类是以湍流的统计结构为基础, 对所有涡旋进行统计平均 大涡模拟把湍流分成大尺度湍流和小尺度湍流, 通过求解三维经过修正的 Naver-Soes 方程, 得到大涡旋的运动特性, 而对小涡旋运动还采用上述的模型 实际求解中, 选用什么模型要根据具体问题的特点来决定 选择的一般原则是精度要高, 应用简单, 节省计算时间, 同时也具有通用性 FLUENT 提供的湍流模型包括 : 单方程 (Spalar-Allmaras) 模型 双方程模型 ( 标准 κ- 模型 重整化群 κ- 模型 可实现 (Realzable)κ- 模型 ) 及雷诺应力模型和大涡模拟 Zero-Equaon Models RANS-based models 包含更多物理机理 One-Equaon Models Spalar-Allmaras Two-Equaon Models Sandard - RNG - Realzable - Reynolds-Sress Model Large-Eddy Smulaon FLUENT 提供的模型选择 每次迭代计算量增加 Drec Numercal Smulaon 湍流模型种类示意图 4
第二节, 平均量输运方程 雷诺平均就是把 Naver-Soes 方程中的瞬时变量分解成平均量和脉动量两部分 对于速度, 有 : u = u + u - 其中, u 和 u 分别是平均速度和脉动速度 (=,,) 类似地, 对于压力等其它标量, 我们也有 : φ φ + φ = -4 其中,φ 表示标量, 如压力 能量 组分浓度等 把上面的表达式代入瞬时的连续与动量方程, 并取平均 ( 去掉平均速度 u 上的横线 ), 我 们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式 : ρ + ( ρu ) = 0 p u u u l = + + δ + ( ρu u l Du ρ ) D -6 上面两个方程称为雷诺平均的 Naver-Soes(RANS) 方程 他们和瞬时 Naver-Soes 方程有相同的形式, 只是速度或其它求解变量变成了时间平均量 额外多出来的项 ρ u u 是雷诺应力, 表示湍流的影响 如果要求解该方程, 必须模拟该项以封闭方程 如果密度是变化的流动过程如燃烧问题, 我们可以用法夫雷 (Favre) 平均 这样才可以求解有密度变化的流动问题 法夫雷平均就是出了压力和密度本身以外, 所有变量都用密度加权平均 变量的密度加权平均定义为 : ~ Φ = ρ Φ / ρ -7 符号 ~ 表示密度加权平均 ; 对应于密度加权平均值的脉动值用 Φ 表示, 即有 : Φ = Φ ~ + Φ 很显然, 这种脉动值的简单平均值不为零, 但它的密度加权平均值等于零, 即 : Φ 0, ρ Φ = 0 Boussnesq 近似与雷诺应力输运模型 为了封闭方程, 必须对额外项雷诺应力 ρ u u 进行模拟 一个通常的方法是应用 Boussnesq 假设, 认为雷诺应力与平均速度梯度成正比, 即 : u u u ρ uu = + ( ρ + ) δ -8 Boussnesq 假设被用于 Spalar-Allmaras 单方程模型和 双方程模型 Boussnesq 近似的好处是与求解湍流粘性系数有关的计算时间比较少, 例如在 Spalar-Allmaras 单方程模型中, 只多求解一个表示湍流粘性的输运方程 ; 在 双方程模型中, 只需多求解湍动能 和耗散率 两个方程, 湍流粘性系数用湍动能 和耗散率 的函数 Boussnesq 假设的缺点是认为湍 -5 5
流粘性系数 是各向同性标量, 对一些复杂流动该条件并不是严格成立, 所以具有其应用限 制性 另外的方法是求解雷诺应力各分量的输运方程 这也需要额外再求解一个标量方程, 通常是耗散率 方程 这就意味着对于二维湍流流动问题, 需要多求解 4 个输运方程, 而三维湍流问题需要多求解 7 个方程, 需要比较多的计算时间, 对计算机内存也有更高要求 在许多问题中,Boussnesq 近似方法可以得到比较好的结果, 并不一定需要花费很多时间来求解雷诺应力各分量的输运方程 但是, 如果湍流场各向异性很明显, 如强旋流动以及应力驱动的二次流等流动中, 求解雷诺应力分量输运方程无疑可以得到更好的结果 第三节, 湍流模型.. 单方程 (Spalar-Allmaras) 模型 Spalar-Allmaras 模型的求解变量是 ~ ν, 表征出了近壁 ( 粘性影响 ) 区域以外的湍流运动粘性系数 ~ ν 的输运方程为 : D ~ ν ~ ν ν ρ Gν ρν Cb ρ ~ = + ( + ~ ) + Yν -9 D σ ~ ν 其中, G 是湍流粘性产生项 ; Y 是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少 ; σ ~ 和 C b ν 是常数 ;ν 是分子运动粘性系数 湍流粘性系数用如下公式计算 : ρν~ f = ν f ν χ + C ν G ν 其中, 是粘性阻尼函数, 定义为 : ν ν f χ =, 并且 χ ~ ν ν 湍流粘性产生项, 用如下公式模拟 : ~ Gν = C ρ ~ b Sν -0 ~ ~ ν χ 其中, S S + f ν, 而 fν = 其中, C b 和 是常数,d 是计算点到 d + χf Ω 壁面的距离 ;S Ω Ω 定义为 : u u Ω = - 由于平均应变率对湍流产生也起到很大作用,FLUENT 处理过程中, 定义 S 为 : S Ω + C mn( 0, S Ω ) - prod 其中, C prod =.0, Ω ΩΩ, S S S, 平均应变率 S 定义为 : S u = u + ν - ν 6
在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小 这适合涡流靠近涡旋中心的区域, 那里只有 单纯 的旋转, 湍流受到抑止 包含应变张量的影响更能体现旋转对湍流的影响 忽略了平均应变, 估计的涡旋粘性系数产生项偏高 湍流粘性系数减少项 Y ν 为 : ~ ν Yν = Cw ρf w -4 d / 6 6 + C w 其中, f w = g 6 6-5 g + Cw g = r + C ( 6 w r r) -6 ~ ν r ~ -7 S d ~ ~ ν 其中, Cw, Cw, Cw 是常数, S S + f ν 在上式中, 包括了平均应变率对 S d 的影响, 因而也影响用 S ~ 计算出来的 r 上面的模型常数在 FLUENT 中默认值为 : C b = 0. 5, C b = 0. 6, σ ~ ν = /, Cν= 7., C w = Cb / + ( + Cb ) / σ ~ ν, C w = 0., C w =. 0, = 0.4 壁面条件 在壁面, 湍流运动粘性 ~ ν 设置为零 当计算网格足够细, 可以计算层流底层时, 壁面切应力用层流应力 - 应变关系求解, 即 : u ρu τ y = -8 u τ 如果网格粗错不能用来求解层流底层, 则假设与壁面近邻的网格质心落在边界层的对数区, 则根据壁面法则 : u ρu = τ y ln E -9 uτ 其中,=0.49,E=9.79 对流传热传质模型 在 FLUENT 中, 用雷诺相似湍流输运的概念来模拟热输运过程 给出的能量方程为 : c p T ( ρ E) + [ u ( ρe + p)] = + u + S h + ( τ ) -0 Pr 式中,E 是总能量, (τ 是偏应力张量, 定义为 : ) 7
u u u ( τ ) = ( + ) δ - 其中, ) ) (τ 表示粘性加热, 耦合求解 如果默认为分开求解,FLUENT 不求解处 (τ 但是可以通过变化 粘性模型 面板上的湍流普朗特数 (Pr), 其默认值为 0.85 湍流质量输运与热输运类似, 默认的 Schmd 数是 0.7, 该值同样也可以在 粘性模型 面板上调节 标量的壁面处理与动量壁面处理类似, 分别选用合适的壁面法则 综上所述,Spalar-Allmaras 模型是相对简单的单方程模型, 只需求解湍流粘性的输运方程, 并不需要求解当地剪切层厚度的长度尺度 该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果, 在透平机械湍流模拟方面也有较好结果 Spalar-Allmaras 模型的初始形式属于对低雷诺数湍流模型, 这必须很好解决边界层的粘性影响区求解问题 在 FLUENT 中, 当网格不是很细时, 采用壁面函数来解决这一问题 当网格比较粗糙时, 网格不满足精确的湍流计算要求, 用壁面函数也许是最好的解决方案 另外, 该模型中的输运变量在近壁处的梯度要比 中的小, 这使得该模型对网格粗糙带来数值误差不太敏感 但是,Spalar-Allmaras 模型不能预测均匀各向同性湍流的耗散 并且, 单方程模型没有考虑长度尺度的变化, 这对一些流动尺度变换比较大的流动问题不太适合 比如, 平板射流问题, 从有壁面影响流动突然变化到自由剪切流, 流场尺度变化明显.. 标准 模型 标准 模型需要求解湍动能及其耗散率方程 湍动能输运方程是通过精确的方程推导 得到, 但耗散率方程是通过物理推理, 数学上模拟相似原形方程得到的 该模型假设流动为完 全湍流, 分子粘性的影响可以忽略 因此, 标准 模型只适合完全湍流的流动过程模拟 标准 模型的湍动能 和耗散率 方程为如下形式 : D ρ = + G + Gb YM D + ρ σ - D ρ = + C ( G + C Gb ) C ρ D + σ - 在上述方程中, G 表示由于平均速度梯度引起的湍动能产生, G 是用于浮力影响引起 的湍动能产生 ; Y M 可压速湍流脉动膨胀对总的耗散率的影响 湍流粘性系数 b ρc = 在 FLUENT 中, 作为默认值常数, C =.44, C =.9, C = 0. 09, 湍动能 与耗散 率 的湍流普朗特数分别为 σ =.0, σ =. 可以通过调节 粘性模型 面板来调节这些 常数值 8
.. 重整化群 κ- 模型 重整化群 κ- 模型是对瞬时的 Naver-Soes 方程用重整化群的数学方法推导出来的模 型 模型中的常数与标准 κ- 模型不同, 而且方程中也出现了新的函数或者项 其湍动能与 耗散率方程与标准 κ- 模型有相似的形式 : D ρ = ( α ) + G + Gb ρ YM D -4 D ρ = ( α ) + C ( G + C Gb ) C ρ R D -5 G 表示由于平均速度梯度引起的湍动能产生, G 是用于浮力影响引起的湍动能产生 ; Y 可 压速湍流脉动膨胀对总的耗散率的影响, 这些参数与标准 κ- 模型中相同 α 和 α 分别是 湍动能 和耗散率 的有效湍流普朗特数的倒数 湍流粘性系数计算公式为 : ρ ~ ν d =.7 d ~ ν ~ ν Cν 其中, ~ ν = /, 00 C ν b M -6 对上面方程积分, 可以精确得到有效雷诺数 ( 涡旋尺度 ) 对湍流输运的影响, 这有助于处理低雷诺数和近壁流动问题的模拟 对于高雷诺数, 上面方程可以给出 : = ρc, C = 0. 0845 这个结果非常有意 = 0. 思, 和标准 κ- 模型的半经验推导给出的常数 C 09 非常近似 在 FLUENT 中, 如果是默认设置, 用重整化群 κ- 模型时候是针对的高雷诺数流动问题 如果对低雷诺数问题进行数值模拟, 必须进行相应的设置 重整化群 κ- 模型有旋修正 通常, 平均运动有旋时候对湍流有重要影响 FLUENT 中重整化群 κ- 模型通过修正湍流粘性系数来考虑了这类影响 湍流粘性的修正形式为 : = f ( α, Ω, ) -7 0 s 其中, 0 是不考虑有旋计算出来的湍流粘性系数 ;Ω 是 FLUENT 计算出来的特征旋流数 ; α s 是旋流常数, 不同值表示有旋流动的强度不同 流动可以是强旋或者中等旋度的 FLUENT 默认设置 α s =0.05, 针对中等旋度的流动问题, 对于强旋流动, 可以选择较大的值 湍动能及其耗散率的有效湍流普朗特数倒数的计算公式为 : 9
式中, 0 0.6 0.679 α.99 α +.99 mol = α 0.99 α 0 +. 99 α =, 在高雷诺数流动问题中, /, α =. 9 mol α = -8 湍流耗散率方程右边的 R 为 : C ρη ( η / η0 ) R = -9 + 其中, βη η S /, η 0 = 4. 8, = 0. 0 β 为了更清楚体现 R 对耗散率的影响, 我们把耗散率输运方程重写为 : D * ρ = ( α ) + C ( G + C Gb ) C ρ C ρ D -0 C * ρη ( η / η0 ) 则 : C = C + + βη - * η < 的区域,R 的贡献为正 ; C 大于 C 以对数区为例, 在 η0 和标准 κ- 模型中给出的 C 群 κ- 模型给出的结果比标准 κ- 模型的结果要大 重整化群模型中, C =. 4, C =.68 * η, C.0, 这 =.9 接近 因此, 对于弱旋和中等旋度的流动问题, 重整化..4 可实现 κ- 模型 可实现 κ- 模型的湍动能及其耗散率输运方程为 : D ρ = + G + Gb YM D + ρ σ - D ρ = ρc S ρc + C C Gb D + + σ + ν - η 其中, C = ma 0.4,, η = S / η + 5 在上述方程中, G 表示由于平均速度梯度引起的湍动能产生, G 是用于浮力影响引起 的湍动能产生 ; 可压速湍流脉动膨胀对总的耗散率的影响 C 和 C 是常数 ; σ, σ 分 YM 别是湍动能及其耗散率的湍流普朗特数 在 FLUENT 中, 作为默认值常数, =.44, C =.9, σ =.0, σ =. b C 可实现 κ- 模型的湍动能的输运方程与标准 κ- 模型和重整化群 κ- 模型有相同的形式, 只是模型参数不同 但耗散率方程有较大不同 首先耗散率产生项 ( 方 0
G 程右边第二项 ) 不包含湍动能产生项, 现在的形式更能体现能量在谱空间的传输 另外的特色在于耗散率减少项中, 不具有奇异性 并不象标准 κ- 模型模型那样把 K 放在分母上 该模型适合的流动类型比较广泛, 包括有旋均匀剪切流, 自由流 ( 射流和混合层 ), 腔道流动和边界层流动 对以上流动过程模拟结果都比标准 κ- 模型的结果好, 特别是可实现 κ- 模型对圆口射流和平板射流模拟中, 能给出较好的射流扩张角 湍流粘性系数公式为 = ρc, 这和标准 κ- 模型相同 不同的是, 在可实现 κ- 模型中, C 不再是个常数, 而是通过如下公式计算 : C = -4 * U K A0 + As * ~ ~ ~ 其中, U = S S + ΩΩ, Ω = Ω ω, Ω = Ω ω, Ω 是 s he mean rae-of roaon ensor vewed n a roang reference frame wh he angular velocy ω 模型常数 A 0 = 4. 04, A s = 6 cosφ, 而 : S S S ~ u u φ = arccos( 6W ), 式中 W= ~, S = S S, S = ( + ) S 我们可以发现, C 是平均应变率与旋度的函数 在平衡边界层惯性底层, 可以得到 C = 0.09, 与标准 κ- 模型中采用底常数一样 双方程模型中, 无论是标准 κ- 模型 重整化群 κ- 模型还是可实现 κ- 模型, 三个模型有类似的形式, 即都有 κ 和 的输运方程, 它们的区别在于 :, 计算湍流粘性的方法不同 ;, 控制湍流扩散的湍流 Prandl 数不同 ;, 方程中的产生项和 G 关系不同 但都包含了相同的表示由于平均速度梯度引起的湍动能产生 G, 用于浮力影响引起的湍动能产生 Gb ; 可压速湍流脉动膨胀对总的耗散率的影响 Y M 湍动能产生项 u G = ρ u u -5 T Gb = β g -6 Pr 式中,Pr 是能量的湍流普特朗数, 对于可实现 κ- 模型, 默认设置值为 0.85 对于重整 ρ 化群 κ- 模型, Pr = / α, α = / Pr = / C p 热膨胀系数 T β =, 对于理想 ρ p 气体, 浮力引起的湍动能产生项变为 :
G b = g ρ ρ Pr -7 在 FLUENT 程序中, 如果有重力作用, 并且流场里有密度或者温度的梯度, 浮力对湍动能的影响都是存在的 浮力对耗散率的影响不是很清楚, 因此, 默认设置中, 耗散率方程中的浮力影响不被考虑 如果要考虑浮力对耗散率的影响, 用 粘性模型 面板来控制 浮力对耗散率影响是用 C 来体现 但 C 并不是常数, 而是如下的函数形式 : v C = anh -8 u v 是平行于重力方向的速度分量 ;u 是垂直于重力方向的速度分量 如果流动速度与重力方向相同的剪切流动, C =, 对于流动方向与重力方向垂直的剪切流, =0 C 对于高马赫数的流动问题, 可压速性对湍流影响在 Y = ρ M M Y M 中体现 其中, M 是马赫数, 定义为 : M = ( a γrt 是声速 ) a 默认设置中, 只要选择可压速理想气体, 可压速效应都是考虑的 在上述的双方程模型中, 对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的 能量方程形式为 : T ( ρ E) + [ u ( ρe + p) ] = u + S h + ( τ ) -9 式中,E 是总的能量, 是有效导热系数 ; (τ 是偏应力张量, 定义为 : ) u u u τ ( ) = + δ -40 (τ 表示的是粘性加热, 耦合求解时总是计算 如果不是耦合求解时候, 作为默认设置, ) 并不求解该量 如果有需要, 需在 粘性模型 面板中设置 对于重整化群 κ- 模型, 有效导热系数为 : = α -4 c p α 用 (-8) 计算, 式中, 0 = / Pr = C p 事实上,α 随着 mol α / / 的变化而 变化, 这是重整化群 κ- 模型的一个优点, 因为实验中证明, 湍流普朗特数随分子普朗特数及湍流而变化 湍流质量输运处理过程与能量输运过程类似 对于标准 κ- 模型和可实现的 κ- 模型, 默认的 Schmd 数是 0.7, 重整化群模型中, 是通过方程 -8 来计算的, 其中, α 0 = / Sc, Sc 是分子 Schmd 数
..5 雷诺应力模型 (RSM) 雷诺应力模型是求解雷诺应力张量的各个分量的输运方程 具体形式为 : ( ρ uu ) + ( ρu uu ) = 对流项 C [ ρ u u u + p( δ u + δ u )] + u u D T 湍流扩撒项分子扩散 L D ( g u θ g u θ ) U U ρ u u + u u ρβ + P 应力产生项浮力产生项目 u u p 压力应变项 Φ + + u u 耗散项 ( u u u u ) ρ Ω + -4 m m 系统旋转产生项 F m m G L T 上面方程中, C, D, P, F 不需要模拟, 而 D, G, Φ, 需要模拟以封闭 方程 下面简单对几个需要模拟项的模拟 D 可以用 Delay and Harlow [L8] 的梯度扩散模型来模拟, 即 : T T uu uu l D = Cs ρ -4 l 但这个模型会导致数值不稳定, 因此 FLUENT 程序中采用标量湍流扩散模型 : T uu D = -44 σ 式中, 湍流粘性系数用 = ρc 来计算, 根据 Len and Leschzner [L98], σ = 0. 8, 这和标准 κ- 模型中选取.0 有所不同 根据 Gbson and Launder [L58], Fu [L55], Launder [L88,L89], 压力应变项可以分解为三 项, 即 : = Φ Φ + Φ,, + Φ -45 w Φ
Φ, Φ 和 Φ w 分别是慢速项, 快速项和壁面反射项,, Φ, = Cρ uu δ, 常数 8 C =. Φ, = C ( P + F + G C ) δ ( P + G C), C = 0. 60, P = P, G = G, C = C 壁面反射项用于重新分布近壁的雷诺正应力分布, 主要是减少垂直于壁面的雷诺正应力, 增加平行于壁面的雷诺正应力 该项模拟为 : / w Φ = C uumn nmδ uu n n u u nn C d / + C Φ m,n nmδ Φ,n n Φ,nn -46 Cld 式中, C =0.5, C = 0., n 是 在垂直于壁面方向上的单位分量,d 是到壁面的距 4 离 ; C l = C / /, C = 0. 09,=0.4 默认设置时候,FLUENT 不计算 Φ 线性压力应变模型 w l 如果需要计算时候, 在 粘性模型 面板中设置 对于小雷诺数流动, 特别是用双层模型求解近壁流动问题时,FLUENT 中通过改进模型常数 C, C, C 和 C 来改进压力应变项 Launder [L9] 这一过程只有在选择双层流模型时候, 在 粘性模型 面板上调节 C = +.58A A e (0.0067 Re ) C = 0. 75 A C = C +. 67 C C = 6, 0 ma C 其中, Re = ρ /( ), 参数 A 和张量不变量 A, A 定义为 : 9 A ( A A ) 8 A α α A α α α 4
式中, α 是雷诺应力张量各向异性部分, 定义为 : ρuu + ρδ α = ρ -47 二阶压力应变模型 二阶压力应变模型由 Spezble {L57} 等人提出 * * Φ = ( Cρ + C P) b + C ρ( bb bmnbmnδ ) + ( C C bb ) ρs + C4 ρ ( b S + b S bmns mnδ ) + C5ρ( bω + b Ω ) -48 式中, b 是雷诺应力各向异性张量, 定义为 : ρuu + ρδ b = -49 ρ u 平均应变率 S 定义为 : u + u S = ; u Ω = ; 模型常数 C =.4, C * =. 8, C = 4., C 8 = 0., C * =., C =. 5 4, C 5 = 0.4 二阶压力应变模型不需要考虑壁面反射影响去模拟对数区湍流边界层过程 浮力对湍流的影响 浮力引起的产生项模拟为 : T T G + = β g g -50 Pr 其中,Pr 是能量的湍流普朗特数, 默认设置值为 0.85 对于理想气体, 把热膨胀系数的定义代入上式, 得 : G = g ρ Pr ρ + g ρ -5 耗散项 的模拟 5
耗散张量 模拟为 : = δ ( ρ + YM ) -5 式中, Y ρ M, M 是马赫数 ; 标量耗散率 用标准 - 模型中的采用的耗散率输运方 程求解 M = 雷诺应力模型的边界条件 在流场进口, 雷诺应力模型需要各个雷诺应力分量和湍动能耗散率的值 这些值可以直接输入, 也可以湍流强度和特征长度来计算 在壁面, 雷诺应力模型通过壁面函数, 给出各个雷诺应力分量和耗散率的值 雷诺应力模型的能量与质量输运方程 在雷诺应力模型中, 对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的 能量方程形式为 : c p T ( ρ E) + [ u ( ρe + p) ] = u + S h ( + ) + ( τ ) -5 Pr 式中,E 是总的能量 ; (τ 是偏应力张量, 定义为 : ) u u u τ ( ) = + δ -54 (τ 表示的是粘性加热, 耦合求解时总是计算 如果不是耦合求解时候, 作为默认设置, ) 并不求解该量, 并且 Pr=0.85 如果有需要, 需在 粘性模型 面板中设置..6 大涡模拟 (LES) 湍流中包含了不同时间与长度尺度的涡旋 最大长度尺度通常为平均流动的特征长度尺度 最小尺度为 Komogrov 尺度 LES 的基本假设是, 动量 能量 质量及其它标量主要由大涡输运 ;, 流动的几何和边界条件决定了大涡的特性, 而流动特性主要在大涡中体现 ;, 小尺度涡旋受几何和边界条件影响较小, 并且各向同性 ; 大涡模拟过程中, 直接求解大涡, 小尺度涡旋模拟, 从而使得网格要求比 DNS 低..6. 大涡模拟的控制方程 6
LES 的控制方程是对 Naver-Soes 方程在波数空间或者物理空间进行过滤得到的 过滤的过程是去掉比过滤宽度或者给定物理宽度小的涡旋, 从而得到大涡旋的控制方程 过滤变量 ( 上横线 ) 定义为 : φ ( ) = D φ( ) G(, ) d -55 其中,D 表示流体区域 ;G 是决定涡旋大小的过滤函数 在 FLUENT 中, 有限控制体离散本身暗中包括了过滤运算, φ ( ) = V φ( ) d, V -56 V 其中 V 是计算控制体体积, 过滤函数为 : / V V G (, ) = -57 0 V 目前, 大涡模拟对不可压流动问题得到较多应用, 但在可压缩问题中的应用还很少, 因此这里涉及的理论都是针对不可压流动的大涡模拟方法 在 FLUENT 中, 大涡模拟只能针对不可压流体 ( 当然并非说是密度是常数 ) 的流动 过滤不可压的 Naver-Soes 方程后, 可以得到 LES 控制方程 : ρ ρu + = 0 ρ u ) + ( ρu u ) = ( ) -58 u p τ ( -59 其中, τ 为亚网格应力, 定义为 : τ = ρu u ρu u -60 很明显, 上述方程与雷诺平均方程很相似, 只不过大涡模拟中的变量是过滤过的量, 而非时间平均量, 并且湍流应力也不同..6. 亚网格模型 由于 LES 中亚网格应力项是未知的, 并且需要模拟以封闭方程 目前, 采用比较多的亚网格模型为涡旋粘性模型, 形式为 : τ τ δ = S -6 式中, 是亚网格湍流粘性系数 ; S 是求解尺度下的应变率张量, 定义为 : u u S = + -6 求解亚网格湍流粘性系数 时,FLUENT 提供了两种方法 第一,Smagornsy-Llly 模型 ; 第二, 基于重整化群的亚网格模型 最基本的亚网格模型是 Smagornsy [L45] 最早提出的,Llly [L99] 把它进行了改善, 这就 7
是今天的 Smagornsy-Llly 模型 该模型的涡粘性计算方程为 : = ρ S -6 L s 式中, Ls 是亚网格的混合长度 ; L s 可以用下式计算 S S S C s 是 Smagornsy 常数, 则亚网格混合长度 / Ls = mn( d, CsV ) -64 其中,=0.4,d 是到最近壁面的距离,V 是计算控制体体积 Llly 通过对均匀各向同性湍流惯性子区湍流分析, 得到了 =0. 但是研究中发现, 对于有平均剪切或者过渡流动中, 该系数过高估计了大尺度涡旋的阻尼作用 因此, 对于比较多的流动问题, C s =0. 有比较好的模拟结果, 该值是 FLUENT 的默认设置值 我们再来看看基于重整化群思想的亚网格模型 人们用重整化群理论推导出了亚网格涡旋粘性系数 [L8], 该方法得到的是亚网格有效粘性系数, =, 而 / C s + s = + H C -65 式中, = ( C V / ) S S,H() 是 Heavsde 函数, s rng > 0 H ( ) = -66 0 0 V 是计算控制体体积 ; 重整化群常数 = 0.57, 而常数 C=00 对于高雷诺数流动 ( >> C rng ),, 基于重整化群理论的亚网格模型就与 Smagornsy-Llly 模型相同, 只是模型常数有区别 在流动场的低雷诺数区域, 上面的函数就小于零, 从而只有分子粘性起作用 所以, 基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果..6. 大涡模拟的边界条件 对于给定进口速度边界条件, 速度等于各个方向分量与随机脉动量的和, 即 : u =< u > + Iψ u 其中,I 是脉动强度,ψ 是高斯随机数, 满足 ψ = 0, ψ = 如果网格足够密并可以求解层流底层的流动的话, 壁面切应力采用线性应力应变关系, 即 : u ρu τ y = -67 uτ 如果网格不够细, 则假定与壁面邻近网格质心落在边界层对数区内, 则 : u ρu ln ( τ y = E ) -68 uτ 其中,=0.48,E=9.79 8
表 -, 雷诺平均模型的比较 模型名字 优点 缺点 Spalar-Allmaras 计算量小, 对一定复杂程度的边界层问题有较好效果 计算结果没有被广泛测试, 缺少子模型, 如考虑燃烧或 标准 应用多, 计算量合适, 有较多数据 积累和相当精度 RNG 能模拟射流撞击, 分离流, 二次流, 旋流等中等复杂流动 Realzable 和 RNG 模型差不多, 还可以模拟圆 口射流问题 雷诺应力模型 考虑的物理机理更仔细, 包括了湍流 各向异性影响 浮力问题对于流向有曲率变化, 较强压力梯度有旋问题等复杂流动模拟效果欠缺受到涡旋粘性各向同性假设限制 受到涡旋粘性各向同性假设限制 CPU 时间长 (~ 倍 ), 动量和湍流量高度耦合 第四节, 湍流模型算例及其设置 湍流模型设置命令 :Defne-model-vscous 9
无粘, 层流和湍流湍流模型选项近壁处理方法选择附加湍流选项 40
算例分析 : 有换热的腔道流动问题 adabac wall cold ar V = 50 fpm T = 0 F nsulaon P f f consan emperaure wall T = 00 F 0 f 步骤 :, 检查是否湍流 : Re Dh = 5980, 选择低雷诺数湍流模型 -RNG 模型 ; 壁面处理用非平衡壁面函数, 考虑压力梯度影响, 网格划分 :() 四边形网格 ;() 由于在靠近水平板处, 垂直方向梯度较大, 则近壁网格加密, 并保证第一个控制体在对数区内 ;() 变化流向网格间距, 用于捕捉边界层发展 ;(4), 根据计算结果, 自动调节网格, 用于进一步计算温度梯度 计算结果 : Velocy conours BLs on upper & lower surfaces accelerae he core flow Temperaure conours P Imporan ha hermal BL was accuraely resolved as well 算例二, 圆柱绕流 4
wall f ar V = 4 fps f f 5 f 4.5 f wall Compue drag cocen of he cylnder 步骤 :, 确定雷诺数, Re D = 4600, 钝体绕流, 后面有不稳定的涡旋脱落 采用 RNG 模型, 壁面处理是双层区模型 ; +, 网格处理 : 近壁网格加密, 由于是双层区模型, 需要网格划分到 y = 4
计算圆柱绕流的涡旋脱落过程 Conours of ecve vscosy = + C D = 0.5 Srouhal Number = 0.97 where D S τ U 4