2023年第三章湍流模型.docx

第三章湍流模型 第三章 湍流模型 第一节 前言 湍流流动模型很多，但大致可以归纳为以下三类： 第一类是湍流输运系数模型，是Boussinesq于1877年针对二维流动提出的，将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积即： 3－1 推广到三维问题，若用笛卡儿张量（笛卡尔坐标系）表示，即有： 3－2 为DELT函数，一般i=j时为1，否则为0. 模型的任务就是给出计算湍流粘性系数的方法根据建立模型所需要的微分方程的数目，可以分为零方程模型（代数方程模型），单方程模型和双方程模型 （模拟大空间建筑空气流动）　 　μt=0.038 74ρvl （模拟通风空调室内的空气流动）　　　　　 比例系数由直接数值模拟的结果拟合而得，其中：v为当地时均速度，l为当地距壁面最近的距离 第二类是抛弃了湍流输运系数的概念，直接建立湍流应力和其它二阶关联量的输运方程 第三类是大涡模拟前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。

实际求解中，选用什么模型要根据具体问题的特点来决定选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性参见：湍流模型的选择资料 FLUENT提供的湍流模型包括：单方程（Spalart-Allmaras）模型、双方程模型（标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型）及雷诺应力模型和大涡模拟 Zero-Equation Models One-Equation Models Spalart-Allmaras Two-Equation Models Standard k-e RNG k-e Realizable k-e Reynolds-Stress Model Large-Eddy Simulation Direct Numerical Simulation 包含更多 物理机理 每次迭代 计算量增加 FLUENT提供的模型选择 RANS-based models 湍流模型种类示意图 大涡模拟启动需要用命令：(rpsetvar ＇les-2d? #t) 第二节 平均量输运方程 输运过程的粘滞系数、扩散系数和热传导率，故称为输运方程 雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。

对于速度，有： 3－3 其中，和分别是平均速度和脉动速度（i=1,2,3） 类似地，对于压力等其它标量，我们也有： 3－4 其中，表示标量，如压力、能量、组分浓度等 把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式：3－5 3－6 上面两个方程称为雷诺平均的Navier-Stokes（RANS）方程他们和瞬时Navier-Stokes方程有相同的形式，只是速度或其它求解变量变成了时间平均量额外多出来的项是雷诺应力，表示湍流的影响如果要求解该方程，必须模拟该项以封闭方程 如果密度是变化的流动过程如燃烧问题，我们可以用法夫雷（Favre）平均这样才可以求解有密度变化的流动问题法夫雷平均就是出了压力和密度本身以外，所有变量都用密度加权平均变量的密度加权平均定义为： 3－7 符号～表示密度加权平均；对应于密度加权平均值的脉动值用表示，即有：很显然，这种脉动值的简单平均值不为零，但它的密度加权平均值等于零，即： ， Boussinesq近似与雷诺应力输运模型 为了封闭方程，必须对额外项雷诺应力进行模拟。

一个通常的方法是应用Boussinesq假设，认为雷诺应力与平均速度梯度成正比，即： 3－8 Boussinesq假设被用于Spalart-Allmaras单方程模型和双方程模型Boussinesq近似的好处是与求解湍流粘性系数有关的计算时间比较少，例如在Spalart-Allmaras单方程模型中，只多求解一个表示湍流粘性的输运方程；在双方程模型中，只需多求解湍动能k和耗散率ε两个方程，湍流粘性系数用湍动能k和耗散率ε的函数Boussinesq假设的缺点是认为湍流粘性系数是各向同性标量，对一些复杂流动该条件并不是严格成立，所以具有其应用限制性 另外的方法是求解雷诺应力各分量的输运方程这也需要额外再求解一个标量方程，通常是耗散率ε方程这就意味着对于二维湍流流动问题，需要多求解4个输运方程，而三维湍流问题需要多求解7个方程，需要比较多的计算时间，对计算机内存也有更高要求 在许多问题中，Boussinesq近似方法可以得到比较好的结果，并不一定需要花费很多时间来求解雷诺应力各分量的输运方程但是，如果湍流场各向异性很明显，如强旋流动以及应力驱动的二次流等流动中，求解雷诺应力分量输运方程无疑可以得到更好的结果。

粘性面板参数说明： Cmu ：(only for the standard or RNG - model or the RSM) is the constant that is used to compute . C1-Epsilon ：(only for the standard or RNG - model or the RSM) is the constant used in the transport equation for . C2-Epsilon ：(only for the standard, RNG, or realizable - model or the RSM) is the constant used in the transport equation for . C1-PS (only for RSM) is the constant in Equation . C2-PS (only for RSM) is the constant in Equation  C1＇-PS (only for RSM) is the constant in Equation . C2＇-PS (only for RSM) is the constant in Equation  Prandtl Number (only for the Spalart-Allmaras model) is the constant in Equation  10.3-1. TKE Prandtl Number (only for the standard or realizable - model, the standard or SST - model, or the RSM) is the effective ``Prandtl＇＇ number for transport of turbulence kinetic energy . This effective Prandtl number defines the ratio of the momentum diffusivity to the diffusivity of turbulence kinetic energy via turbulent transport. TKE (Inner) Prandtl # (only for the SST - model) is the effective ``Prandtl＇＇ number for the transport of turbulence kinetic energy, , inside the near-wall region. See Section  10.5.2 for details. TKE (Outer) Prandtl # (only for the SST - model) is the effective ``Prandtl＇＇ number for the transport of turbulence kinetic energy, , outside the near-wall region. See Section  10.5.2 for details. TDR Prandtl Number is the effective ``Prandtl＇＇ number for transport of the turbulent dissipation rate, , for the standard or realizable - model or the RSM. This effective Prandtl number defines the ratio of the momentum diffusivity to the diffusivity of turbulence dissipation via turbulent transport. For the standard - model, the TDR Prandtl Number is the effective ``Prandtl＇＇ number for the transport of the specific dissipation rate, . SDR (Inner) Prandtl # (only for the SST - model) is the effective ``Prandtl＇＇ number for the transport of the specific dissipation rate, , inside the near-wall region. See Section  10.5.2 for details. SDR (Outer) Prandtl # (only for the SST - model) is the effective ``Prandtl＇＇ number for the transport of the specific dissipation rate, , outside the near-wall region. See Section  10.5.2 for details. Dispersion Prandtl Number (only for the - multiphase models) is the effective ``Prandtl＇＇ number for the dispersed phase, . See Section  22.4.8 for details. Energy Prandtl Number (for any turbulence model except the RNG - model) is the turbulent Prandtl number for energy, Pr , in Equation  10.6-20. (This item will not appear for premixed or partially premixed combustion models.) Wall Prandtl Number (for all turbulence models) is the turbulent Prandtl number at the wall, Pr in Equation  10.8-5. (This item will no。