湍流理论和湍流模型(博士课程)资料.ppt

S-A模型的数值求解 湍流模型求解与N-S方程组的求解可采用“松耦合”的方式，即在同一次时间推进中它们的求解是相对独立的。S-A湍流模型方程的对流项和耗散项采用中心格式进行有限体积离散。 这样，在第 i 个网格单元上有 （左端第一项） （左端第二项） （右端第一项） （右端第二项） （右端第三项） （右端第四项） 经过上述空间离散后，湍流模型方程可以写为下面的半离散形式 其中， 是湍流模型方程的残值项。 采用隐式欧拉方法进行时间离散求解，表达式如下， 将方程右端的残值项进行一阶泰勒展开，略去高阶项后得到 式中的偏导数 可以通过数值求解雅可比矩阵的方法求得，即 其中 为一个很小的正数，一般取 经过整理后，方程变为 这个方程采用高斯-赛德尔迭代进行求解。 边界条件方面，固壁面边界上 ；自由流处 。 初始条件可将 取定为自由流的值。 （三）二方程模式(标准 k-ε 模式) k- ε 模式是在涡粘模式的基础上发展起来的，它和代数模式的主要区别是在于k- ε模式的涡粘系数包含部分历史效应。具体来说，它把涡粘系数和湍动能及湍动能耗散联系在一起。 湍动能耗散率的输运方程为 用量纲分析，涡粘系数可以写成 在湍流模式中，湍动能耗散是最难构造准确模型的。 通常采用的ε模式是依据类比方法，基本思想是：湍动能耗散的生成、扩散以及消耗等项与湍动能方程中的对