尺度自适应模拟(SAS)模型

尺度自适应模拟(SAS)模型

本节介绍尺度自适应模拟(SAS)模型的理论基础。关于在Ansys Fluent中使用该模型的详细信息,请参阅用户指南中的“设置尺度自适应模拟(SAS)建模”部分。

更多信息,请参见以下章节:

4.12.1 概述

4.12.2 SST-SAS模型的传输方程

4.12.3 与其他基于ω​的湍流模型结合的SAS

4.12.1 概述

尺度自适应模拟(SAS)是一种改进的URANS公式,它允许在不稳定流动条件下解析湍流谱。图4.6:横流中圆柱的解析结构(顶部:URANS;底部:SST-SAS)(第107页)显示了使用SST模型(URANS)和SST-SAS模型计算的横流中圆柱的等值面,其中S−Ω=1051/s2(Re=3.6∗106​)。URANS模拟仅产生大尺度非定常性,而SST-SAS模型以动态方式适应已解析的尺度,并允许在分离区域发展湍流谱。图4.6:横流中圆柱的解析结构(顶部:URANS;底部:SST-SAS)

SAS概念基于将von Kármán长度尺度引入湍流尺度方程。von Kármán长度尺度提供的信息使SAS模型能够动态适应URANS模拟中的已解析结构,从而在流动场的不稳定区域产生类似LES的行为。同时,该模型在稳定流动区域提供标准的RANS功能。

关于入口合成湍流生成的信息,请参阅“尺度解析模拟的入口边界条件”(第126页)。

4.12.2 SST-SAS 模型的输运方程

从基本原理的角度来看,目前使用的所有两方程模型都缺乏一个基础的精确输运方程,该方程可以作为逐项模型开发的指导。这种缺陷的原因在于,ε 的精确方程描述的不是大尺度,而是耗散尺度。然而,两方程模型的目标是对大尺度运动对平均流动的影响进行建模。由于缺乏精确方程,ε 和 ω 方程是根据湍流动能 k 的方程,使用纯启发式论据进行建模的。Rotta(1968, 1972)开发了一种更一致的方法来制定尺度方程。Rotta 没有使用纯启发式和量纲论据,而是为湍流动能乘以长度尺度 kL 制定了一个精确的输运方程。Rotta 的方程代表了湍流的大尺度,因此可以作为逐项建模的基础。Ansys Fluent 中实现的 SST-SAS 模型的输运方程基于将 Rotta 的方法转换为 k−ω(SST)并定义为:

∂t∂ρk​+∂xi​∂​(ρui​k)=Gk​−ρcμ​kω+∂xj​∂​[(μ+σk​μt​​)∂xj​∂k​](4.263)

∂t∂ρω​+∂xi​∂​(ρui​ω)​=αkω​Gk​−ρβω2+QSAS​+∂xi​∂​[(μ+σω​μt​​)∂xj​∂ω​]+(1−F1​)σω,2​2ρ​ω1​∂xj​∂k​∂xj​∂ω​​(4.264)

详细推导请参阅 Egorov 和 Menter [431](第 1082 页)。SST-SAS 模型的传输方程(方程 4.263,第 107 页)和方程 4.264(第 107 页)与 SST RANS 模型 [428](第 1081 页)的不同之处在于,在湍流涡频率 ω 的传输方程(方程 4.264,第 107 页)中增加了 SAS 源项 QSAS​。在方程 4.264(第 107 页)中,σω,2​ 是 SST 模型中 k−ε 区域的 σω​ 值。

额外的源项 QSAS​ 如下所示(详细内容请参阅 Egorov 和 Menter [431],第 1082 页):

QSAS​=max[ρη2​κS2(Lνκ​L​)2−C⋅σϕ​2ρk​max(ω21​∂xj​∂ω​∂xj​∂ω​,k21​∂xj​∂k​∂xj​∂k​),0](4.265)

这个SAS源项源自Rotta输运方程[430](第1082页)中的二阶导数项。在SAS源项方程4.265(第108页)中的模型参数是...

η2​=3.51

σϕ​=2/3

C=2

这里,L 代表所模拟湍流的长度尺度。

L=k​/(cμ1/4​⋅ω)(4.266)

冯·卡门长度尺度 ( L_{vK} ) 是经典边界层定义的三维推广,即 ( L_{BL}^{vK} = \kappa \frac{\partial U}{\partial y}/\frac{{\partial }^{2}U}{\partial {y}^{2}} )。

LvK​=κ​U′′U′​​(4.267)

速度的第一个导数 ∂y∂U​ 在 LvK​ 中由 U′ 表示,其等于 S,即应变率张量 Sij​ 的标量不变量。

U′=S=2⋅Sij​Sij​​;Sij​=21​(∂xj​∂Ui​​+∂xi​∂Uj​​)(4.268)

请注意,同样的 S 也直接参与了 QSAS​(公式 4.265(第 108 页))以及湍流生成项 Pk​=μt​S2 中。第二个速度导数 U′′ 通过速度拉普拉斯算子的模来推广到三维:

U′′=∂xk2​∂2Ui​​∂xj2​∂2Ui​​​(4.269)

因此定义的 L 和 LvK​ 在边界层的对数部分都等于 (κy),其中 κ=0.41 是冯·卡门常数。

该模型还提供了对高波数阻尼的直接控制。这是通过以下方式对 LvK​ 值设置较低约束来实现的:

LvK​=max(κ​U′′U′​​,CS​(β/cμ​)−ακη2​​​⋅Δ),Δ=ΩCV1/3​(4.270)

此限制器与网格单元尺寸 Δ 成正比,该尺寸计算为控制体积大小 ΩCV​ 的立方根。此限制器的目的在于控制对最精细解析湍流脉动的阻尼。限制器的结构来源于对 SST-SAS 模型 [431](第 1082 页)平衡涡粘性的分析。

4.12.3 与其他基于 ω 的湍流模型的 SAS

SAS 方法不仅可以与 SST 模型结合,还可以与几种其他基于 ω 的湍流模型结合:

标准 k−ω、BSL k−ω 和广义 k−ω(GEKO)

过渡 SST

基于 ω 的雷诺应力模型(RSM)和显式代数雷诺应力模型(EARSM)

将 SAS 功能添加到这些模型之一中,将在基于 ω 的方程中启用额外的源项 QSAS​(如公式 4.265(第 108 页)所定义)。如果 SAS 功能与 BSL(BSL k−ω、Stress-BSL、WJ-BSL-EARSM)结合,则公式 4.270(第 108 页)中的系数 α 和 β 是 BSL 模型的系数。如果与 GEKO 结合,则 α 和 β 由 GEKO 模型系数确定。对于所有其他组合,使用 SST 模型的系数。此外,在 SAS 项中,这些系数不包括压缩性效应。

相关

[衣]头版头条:飞利浦大金刚和小金刚哪个好?分享区别如何?深度体验感受!![复制链接]
揭秘网站品质认证:如何选择可靠保证机构,守护你的网络安全与信任
英雄联盟补丁包怎么用?LOL手动补丁包怎么用