本构关系

本构关系描述应力张量和用于 Boussinesq 近似的平均应变率之间的关系。

默认情况下,Boussinesq 近似表示线性本构关系。通过添加应变和旋转张量的非线性函数,非线性本构关系 [311] 可以考虑湍流各向异性。两个非线性本构关系适用于 Simcenter STAR-CCM+ 中的标准 K-Epsilon 模型变体(SKE、SKE 2L 和 SKE LRe):

本构关系 公式
二次
1. EQUATION_DISPLAY
TRANS,NL=-4μtkε{C1[SS-13I(S:S)]+C2(WS+SWT)+C3[WWT-13I(W:WT)]}
(1203)

其中:

C1=CNL1(CNL6+CNL7S¯3)Cμ

C2=CNL2(CNL6+CNL7S¯3)Cμ

C3=CNL3(CNL6+CNL7S¯3)Cμ

三次
2. EQUATION_DISPLAY
TRANS,NL=TRANS,quad-8μtk2ε2{C4[(SS)W+WT(SS)]+C5(S:S-W:WT)[S13tr(S)I]}
(1204)

其中,TRANS,quadEqn. (1203) 给出,且:

C4=CNL4Cμ2

C5=CNL5Cμ2

其中:

系数 Cμ 由以下公式给出:

3. EQUATION_DISPLAY
Cμ=Ca0Ca1+Ca2S¯+Ca3W¯
(1205)

其中:

  • Ca0Ca1Ca2Ca3模型系数
  • S¯=kε2S:S
  • W¯=kε2W:W

使用非线性本构关系时,变量系数 Cμ (Eqn. (1205)) 将替换湍流粘度 μt 关系中的常数值 Cμ。请参见 Eqn. (1163)

对于动量方程 Eqn. (665),应力张量 TRANS 定义如下:

4. EQUATION_DISPLAY
TRANS=TRANS,L+TRANS,NL
(1206)

其中,TRANS,LEqn. (1147) 给出。

模型系数

CNL1 CNL2 CNL3 CNL4 CNL5 CNL6 CNL7
0.75 3.75 4.75 -10 -2 1000 1
Ca0 Ca1 Ca2 Ca3
0.667 1.25 1 0.9