[发明专利]一种Reynolds时均法建立湍流模型的构建方法

权利要求书

1.一种Reynolds时均法建立湍流模型的构建方法，其特征在于，它包括如下步骤：

(1)采用Reynolds时均法把湍流的瞬时运动看作是由平均值和脉动值所对应的两个流动叠加而成，在时均系下，在直角坐标形式的绝对参照系下，建立瞬时流动状态的控制方程为：

∂∂t(ρui)+∂∂xj(ρuiuj)=-∂p∂xi+∂∂xj(μ∂ui∂xj-ρui′uj′‾)+Si]]>

(2)引入涡粘系数μ_t，通过涡粘系数μ_t建立了Reynolds应力与平均速度梯度的关系：

-ρui′uj′‾=μt(∂ui∂xj+∂uj∂xj)-23(ρk+μt∂ui∂xj)δij]]>

式中，μ_t为涡粘系数，u_i为时均速度，δ_ij是“Kronecker delta”符号，k为湍动能；

(3)依据确定涡粘系数μ_t的微分方程数目，提出了标准k-ε模型；在标准k-ε模型中，湍动能k和耗散率ε是两个基本未知量，对应的输运方程为：

∂(ρk)∂t+∂(ρkui)∂xi=∂∂xj[(μ+μtσk)∂k∂xj]+Gk+Gb-ρϵ-YM+Sk]]>

∂(ρϵ)∂t+∂(ρϵui)∂xi=∂∂xj[(μ+μtσk)∂ϵ∂xj]+C1ϵϵk(Gk+C3ϵGb)-C2ϵρϵ2k+Sϵ]]>

其中，G_k是由平均速度梯度引起的湍动能k的产生项，G_b是由浮力引起的湍动能k的产生项，Y_M代表可压湍流中脉动扩张的贡献，C_1ε、C_2ε和C_3ε为经验常数，σ_k和σ_ε分别是与湍动能k和耗散率ε对应的Prandtl数，S_k和S_ε是用户定义的源项；

(4)将步骤(1)(2)(3)方程离散为代数方程，建立离散方程，有

(aPa)φP=Σnbanbφnb+b+(1-a)aPaφPo]]>

式中，φ_nb是φ_p的邻点速度，a_P和a_nb为系数，φ_P^o为上一层次之解，a为松弛因子，是代数方程源项，b为源项中的常数部分；

(5)采用分离式解法，顺序地、逐个地求解离散方程中各变量代数方程组；

(6)检查结果是否收敛。若不收敛，以得到的结果作为新的猜测，重复该过程。