K-Epsilon 模型参考

K-Epsilon 湍流模型可用于计算湍动能 k 和湍流耗散率 ϵ ,以提供封闭的雷诺平均纳维-斯托克斯方程。

模型概述

模型名称和缩写 EB K-Epsilon(椭圆混合 K-Epsilon) EB KE
滞后 EB K-Epsilon(滞后椭圆混合 K-Epsilon) EBL KE
可实现 K-Epsilon RKE
可实现的 K-Epsilon 两层模型* RKE 2L
标准 K-Epsilon SKE
标准 K-Epsilon 低雷诺数(标准 K-Epsilon 低雷诺数) SKE LRe
标准 K-Epsilon 两层模型 SKE 2L
理论 请参见
提供方式 [物理连续体] > 模型 > K-Epsilon 湍流模型
节点路径示例 连续体 > Physics 1(物理 1) > 模型 > 可实现的 K-Epsilon 两层模型
要求 物理模型
  • 空间:轴对称、二维或三维
  • 时间:稳态、隐式非稳态、PISO 非稳态或谐波平衡(除 EB KE 和 EBL KE 以外的所有模型)
  • 材料:气体、液体、多相、多组分气体或多组分液体
  • 粘滞态:湍流
  • 湍流:雷诺平均纳维-斯托克斯
  • 雷诺平均湍流:K-Epsilon 湍流
属性 请参见属性查找。
激活 物理模型
  • 壁面距离:壁面距离*
  • 壁面处理:全 y+ 壁面处理和低 y+ 壁面处理(SKE LRe、EB KE 和 EBL KE)、高 y+ 壁面处理(RKE 和 SKE)或两层全 y+ 壁面处理(RKE 2L 和 SKE 2L)
  • 可选模型:湍流粘度用户缩放和温度通量模型 (SKE LRe)
模型控制(子节点)
初始条件 请参见 K-Epsilon 初始条件参考。
边界输入 请参见 K-Epsilon 边界参考。
区域输入 请参见 K-Epsilon 区域参考。
求解器
  • K-Epsilon 湍流(除 EB KE 和 EBL KE 以外的所有模型)
  • EB K-Epsilon 湍流 (EB KE)
  • 滞后 K-Epsilon 湍流 (EBL KE)
  • K-Epsilon 湍流粘度
请参见 K-Epsilon 求解器参考。
监视器
  • Tdr(湍流耗散率)
  • Tke(湍动能)
  • Alpha(椭圆混合函数、EB KE 和 EBL KE)
  • Phi(折算应力函数、EB KE 和 EBL KE)
场函数
  • Effective Viscosity(有效粘度)
  • 椭圆混合函数(EB KE、EBL KE)
  • Kolmogorov 长度尺度
  • Kolmogorov 时间尺度
  • 折算应力函数(EB KE、EBL KE)
  • Strain Rate Tensor Modulus(应变率张量模量)
  • 泰勒微尺度
  • 湍流壁面距离雷诺数(RKE 2L、SKE LRe、SKE 2L)
  • 湍流耗散率
  • 湍动能
  • Turbulent Viscosity(湍流粘度)
  • 湍流粘度比

请参见 K-Epsilon 场函数参考。

属性查找

下表显示不同的 K-Epsilon 湍流模型使用了哪些属性。使用上述“模型概述”中给定的缩写。

SKE SKE LRe SKE 2L RKE RKE 2L EB KE EBL KE
α 最小值
已传输变量 α 允许的最小值。适当的值是一个大于计算机最小浮点数的较小数值。
耗散的浮力结果
确定如何计算结果项 P ε 中的系数 C ε 3 (请参见 Eqn. (1155) 和 Eqn. (1253)(EB KE 和 EBL KE)。
  • 无:将 C ε 3 设为零。
  • 边界层方向:根据 Eqn. (1177) 计算 C ε 3 。
  • 热分层:根据 Eqn. (1178) 计算 C ε 3 。
  • 常系数:将 C ε 3 作为常系数计算。此选项需要指定对应子节点 C3e 中的 C ε 3 。
C 系数
Eqn. (1173) 中的系数 C 。
C_T
Eqn. (1250) 中的系数 C T 。
C1
Eqn. (1260) (EB KE) 和 Eqn. (1261) (EBL KE) 中的系数 C 1 。
C1e
基本传输方程中的系数 C ε 1 。
C1s
Eqn. (1261) 中的系数 C 1 * 。
C2
Eqn. (1260) (EB KE) 中的系数 C 2 。
C2e
基本传输方程中的系数 C ε 2 。
C3
Eqn. (1268) 中的系数 C 3 。
C3e
Eqn. (1259) 中的系数 C ε 3 。
C3s
Eqn. (1261) 中的系数 C 3 * 。
C4
Eqn. (1261) 中的系数 C 4 。
C5
Eqn. (1261) 中的系数 C 5 。
Cd0 系数
Eqn. (1174) 中的系数 C d 0 。
Cd1 系数
Eqn. (1174) 中的系数 C d 1 。
Cd2 系数
Eqn. (1174) 中的系数 C d 2 。
Ceta
计算湍流长度尺度 L 时使用的系数 C η 。
Ck
Eqn. (1265) 中的系数 C k 。
Cl
计算湍流长度尺度 L 时使用的系数 C L 。
Cmu
湍流粘度 μ t 计算和基本传输方程中的系数 C μ 。
本构关系
控制使用的本构关系类型。
  • 线性:选择 Boussinesq 近似隐含的线性本构关系(请参见 Eqn. (1129))。对大多数模拟均使用此选择。
  • 二次:选择二次本构关系(请参见 Eqn. (1189))。此关系需要系数 C μ 、 C 1 、 C 2 和 C 3 。
    非线性 Cmu 系数
    用于提供参数 Ca0、Ca1、Ca2 和 Ca3 以计算 C μ 的子节点。
    非线性二次系数
    用于提供参数 Cnl1、Cnl2、Cnl3、Cnl6 和 Cnl7 以计算 C 1 、 C 2 和 C 3 的子节点。
  • 三次:选择三次本构关系(请参见 Eqn. (1190))。此关系需要系数 C μ 、 C 1 、 C 2 、 C 3 、 C 4 和 C 5 。
    非线性 Cmu 系数
    与二次相同。
    非线性二次系数
    与二次相同。
    非线性三次系数
    用于提供参数 Cnl4、Cnl5 以计算 C 4 和 C 5 的子节点。

有关何时使用非线性本构关系的准则,请参见考虑湍流各向异性。

对流
控制对流格式。
  • 一阶:选择一阶迎风对流格式。
  • 二阶:选择二阶迎风对流格式。
Ct
计算湍流时间尺度 T 时使用的系数 C t 。
曲率校正选项
控制是否要应用曲率校正以提供强(流线)曲率和坐标系旋转效果。
打开时,将应用曲率校正因子 f c ,以根据局部旋转速率和涡旋速率更改湍动能结果项(请参见 Eqn. (1163) 至 Eqn. (1166))。仅当具有稳定求解时,才可使用此选项以改善结果。
曲率校正参数
用于提供参数 Cr1、Cr2、Cct 和上限 ( C max ⁡ ) 以使用 Eqn. (1272) 来计算 f c 的子节点。

有关何时应用曲率校正的准则,请参见考虑强流线曲率和坐标系旋转。

D 系数
Eqn. (1170) 中的系数 D 。
E 系数
Eqn. (1170) 中的系数 E 。
自由流选项
控制 ε 传输方程中 C ε 2 * 的定义(请参见 Eqn. (1253)):
  • Variable C2e Option(变量 C2e 选项):使用变量 C ε 2 * 值(请参见 Eqn. (1266))。默认情况下设置为使用该选项。
  • 关闭:使用常数 C ε 2 * 值(请参见 Eqn. (1267))。
在某些外部空气动力情况下,变量 C2e 选项会导致自由流中出现不稳定。在这种情况下,需将自由流选项设为关闭。
法向应力项
此属性是一个显式项,根据完整的 Boussinesq 近似直接合并发散和湍动能 - 2 3 ρk I 。

打开时,应力张量建模如下:

T RANS = 2 μ t S - 2 3 ( μ t ∇ ⋅ v + ρ k ) I

湍流结果使用下式建模:

G k = μ t S 2 − 2 3 ρ k ∇⋅ v ¯ − 2 3 μ t ( ∇⋅ v ¯ ) 2

关闭时,应力张量建模如下:

T RANS = 2 μ t S - 2 3 ( μ t ∇ ⋅ v ) I

湍流结果假设为:

G k = μ t S 2 − 2 3 μ t ( ∇⋅ v ¯ ) 2

默认情况下,此属性处于关闭状态,在这种情况下,物理量 - 2 3 ρ k I 将并入压力,导致压力略有不同。但即使未显式包括此项,结果仍会相同。

在不可压缩流中,仅压力梯度起作用,因此此设置对结果没有影响。但在可压缩流中,理想气体定律 (Eqn. (630)) 中将使用压力绝对值。如果 k 初始化效果不佳,激活此选项会影响求解收敛。

Φ 最小值
已传输变量 φ 允许的最小值。适当的值是一个大于计算机最小浮点数的较小数值。
可实现性选项
选择是否将限制器应用于湍流时间尺度。可用的选项包括:
  • 无 - 不应用限制器。
  • 德宾尺度限制器 - 在 Eqn. (1151) 中激活对湍流时间尺度 T 的德宾可实现性最小约束。
  • 涡旋限制器 - 在 Eqn. (1152) 和 Eqn. (1153) 中激活湍流时间尺度的涡旋限制器。激活后,提供涡旋系数节点。仅当选择 VOF 波模型时才可用。
可实现系数
指定 Eqn. (1151) 中的可实现系数 C T 。

有关何时激活德宾可实现性约束的准则,请参见克服 K 意外大幅增长。

涡旋系数
指定用于计算湍流时间尺度 T 的涡旋系数。这是 Eqn. (1152) (SK) 和 Eqn. (1153) (RKE) 中的项 λ 2 C ϵ 2 C ϵ 1 。

有关何时激活涡旋限制器的准则,请参见自由表面波的湍流阻尼。

Sarkar
可压缩性修正 ϒ M 中的系数 C M (请参见 Eqn. (1171))。对于 EB KE 和 EBL KE,请参见 Eqn. (1264)。
二阶梯度
忽略或包括用于扩散的边界二阶梯度和/或网格面上的内部二阶梯度。
  • 打开:包括两个二阶梯度。
  • 关闭:排除两个二阶梯度。
  • 仅限内部:仅包括内部二阶梯度。
  • 仅限边界:仅包括边界二阶梯度。
Sigma_e
基本传输方程中的系数 σ e 。
Sigma_k
基本传输方程中的系数 σ k 。
Sigma_phi
Eqn. (1254) 中的系数 σ φ 。
Tdr 最小值
已传输变量 ε 允许的最小值。适当的值是一个大于计算机最小浮点数的较小数值。
Tke 最小值
已传输变量 k 允许的最小值。适当的值是一个大于计算机最小浮点数的较小数值。
两层 δ ReY
Eqn. (1181) 中的 Δ Re y 值。
两层 ReY*
Eqn. (1180) 中的 Re y * 值。
两层型
控制两层公式的类型
  • 剪切驱动 (Wolfstein):选择 Wolfstein [310] 的两层公式。对于浮力不占主导的流体,使用此选择。请参见 Eqn. (1183) 和 Eqn. (1184)。
  • 浮力驱动(徐氏):选择徐氏 [311] 的两层公式。对于浮力占主导的流体,使用此选择。请参见 Eqn. (1187) 和 Eqn. (1188)。
  • 剪切驱动 (Norris-Reynolds):选择 Norris 和 Reynolds [305] 的两层公式。对于浮力不占主导的流体,使用此选择。请参见 Eqn. (1185) 和 Eqn. (1186)。
Yap 校正
打开时,Yap 校正 ϒ y 包括在 ε 的传输方程中,请参见 Eqn. (1160) (SKE 2L) 和 Eqn. (1162) (SKE LRe)。默认设置为关闭。
Yap Cl
Yap 校正中的系数 C l (请参见 Eqn. (1172))。
Yap Cw
Yap 校正中的系数 C w (请参见 Eqn. (1172))。