【CFD理论】k-ω模型

目录

           k-ω模型

标准k-ω模型

低雷诺数修正

剪切流校正

输运方程

BSL k-ω模型

输运方程

混合函数

SST k-ω模型

湍流选项

Buoyancy Effects

Viscous Heating

Curvature Correction（曲率修正）

Corner Flow Correction（拐角流修正）

Production Limiters

---

本文基于fluent theory guide和user guide，结合自己的经验，对k-ω模型进行以下总结。

k-ω模型

标准k-ω模型是基于湍流动能（k）和比耗散率（ω）的模型输运方程的经验模型，也可以认为是 ε与 k 的比值。

该模型包含对低雷诺数效应、可压缩性和剪切流扩展的修改。多年来， 随着模型的修改， 湍流动能在和比耗散率在输运方程中都添加了生产项，这提高了模型预测自由剪切流的准确性。

标准k-ω模型

粘性项建模如下：

低雷诺数修正

与k-ε方程不同，此处使用α进行低雷诺数修正，可以在使用k-ω时勾选上Low-Re Corrections来启动对于低雷诺数的湍流粘度修正。

在低雷诺数时：

在高雷诺数形式中：

此修正可以用于标准、BSL、SST模型之中

剪切流校正

在标准k-ω模型中，可以勾选上Shear Flow Corrections进行剪切流校正，默认是启动的。

输运方程

湍流动能k和比耗散率 ω由以下输运方程计算得出：

有效扩散率由下式给出：

一般来说不建议使用标准 k-ω模型，下面的两个k-ω都比它更好。

BSL k-ω模型

相比于标准K-ω模型，BSL k-ω模型能够有效地将近壁区域 k- ω模型的鲁棒性和精确与k-ε模型在远场中的自由流独立性相结合。主要变化如下：

• 使用混合函数将模型结合，混合函数在近壁区域激活k- ω模型，当远离壁面时混合函数为0，激活k-ε模型

• BSL模型在 ω方程中增加了一个Dω项

• 建模常数不同

输运方程

输运方程基本和标准模型相似，只在ω方程中增加了Cross-Diffusion Modification项。

BSL 模型 基于标准 K-ω 模型和标准 K-ε 模型。为了将这两个模型混合在一起，标准 K-ε 模型被转换为基于 K和ω的方程，这导致了交叉扩散项的引入。该项具体构成如下：

混合函数

有效扩散率依然由下式给出

湍流动能和比耗散率的湍流普朗特数如下。其中加入了F1以控制它们的值，简单的理解就是：当靠近壁面时F1为1，远离时为0。

混合函数F1具体计算式如下所示：

SST k-ω模型

SST k-ω模型包含了BSL k-ω模型的所有改进，此外还在湍流粘度的定义中考虑了湍流剪切应力的传递。可以说SST模型比标准模型和BSL 模型更加准确和可靠，适用于更广泛的流态。

湍流粘度建模如下：

湍流选项

当你不知道你为什么要使用这些选项时，一般保持默认即可。当有需要时，可以针对性的使用。

Buoyancy Effects

可以使用以下选项

• Buoyancy Effects: Off，此选项忽略了浮力对湍流的影响

• Buoyancy Effects: Only Turbulence Production

• Buoyancy Effects: Full

当重力场和温度梯度同时存在时，由于浮力会产生湍流动能k。

在湍流动能k的输运方程中，由于浮力引起的湍流由下式给出：

在k的输运方程中，可以看出，当在不稳定分层中，湍流动能趋于增强（ Gb>0），对于稳定的分层，浮力倾向于抑制湍流（Gb<0 ）。

然而浮力对于ε项的影响不是很明晰，一般也是用Gb表示，见下式输运方程：

至于ω方程中的浮力项，由以下关系推导：

最终得到浮力项如下：

默认情况下，即选择Only Turbulence Production时，仅包含上式的第二部分，而仅当在选择Buoyancy Effects: Full，才包含第一部分。

Viscous Heating

当使用压力基的求解器时，Fluent能量方程中不会计算粘性耗散项（粘性剪切产生的热能）。

这部分热量通常很小，但是当布林克曼数Br接近1时，就值得关注了。此时可以勾选Viscous Heating进行计算。

对于密度基求解器时，求解能量方程时始终包括这一项。

Curvature Correction（曲率修正）

流线曲率和系统旋转在许多具有实际意义的湍流中起着重要作用，然而涡流粘度模型对它们并不敏感。所有Fluent提供了曲率修正，勾选曲率修正后，可以指定 CCURV 的值，以影响特定流所需的曲率校正强度。

值得注意的是，RNG和Realizable湍流模型本身就包括旋转或涡流效应修正。最好谨慎选择

Corner Flow Correction（拐角流修正）

对于拐角流，可以勾选Corner Flow Correction，然后设置系数 CCORNER。

Production Limiters

标准双方程湍流模型的一个缺点是在停滞点附近产生过多的湍流能量Gk，为了避免湍流动能在停滞区积聚，有两种方式可用于限制湍流动能方程中的生产项。

• Production Kato-Launder

注意到湍流动能产生的过量水平是由停滞区非常高的剪切应变率S引起的。由于停滞点附近的流场几乎是无旋转的，涡度率Ω非常小，故提出Gk替代项如下：

在简单的的剪切流中， S和 Ω相等，方程等价于原来的形式。该方法可用于k-ω模型和K-ε模型（除了realizable）

• Production Limiter

限制公式如下，其中Clim默认值为10：

可用于K-ε和k-ω模型，默认情况下，此限制器为基于 k-ω方程的所有湍流模型启用。

值得注意的是，这两种方法可以一同使用