高等流体力学第四章..ppt

4.5 大涡模拟（LES） 实现大涡模拟所需的工作 建立一种数学滤波函数，从湍流瞬时运动方程中将尺度比滤波函数的尺度小的涡滤掉，从而分解出描写大涡流场的运动方程； 建立附加应力项的数学模型，以体现被滤掉的小涡对大涡运动的影响。该模型称为亚格子尺度模型（Sub Grid-Scale model)，简称SGS模型。 4.5.1 大涡的运动方程 在LES方法中，通过使用滤波函数，每个变量都被分成两部分。 举例 对于瞬时变量φ，有： a．大尺度的平均分量 。该部分叫做滤波后的变量，是在LES模拟时直接计算的部分。 b．小尺度分量 。该部分是需要通过模型来表示的。 4.5.1 大涡的运动方程 注意 这里的平均分量 是滤波后得到的变量，它不是在时间域上的平均，而是在空间域上的平均。 式中，D是流动区域； 是实际流动区域中的空间坐标； X是滤波后的大尺度空间上的空间坐标； 是滤波函数。 (4-44) 4.5.1 大涡的运动方程 决定了所求解的涡的尺度，即将大涡与小涡划分开来； 换句话说， 只保留了φ在大于滤波函数 宽度的尺度上的可变性。 的表达式有多种选择，但有限体积法的离散过程本身就隐含地提供了滤波功能，即在一个控制体积上对物理量取平均值，因此，这里采用如下的表达式: 其中V是表示控制体积 所占几何空间的大小 (4-45) 4.5.1 大涡的运动方程 式（4-44）可以写成： 用式(4-46)表示的滤波函数处理瞬时状态下的Navier-Stokes方程及连续方程，有： LES方法中瞬时状态下使用的控制方程组 ，式中带有上划线的量为滤波后的场变量 (4-46) (4-47) (4-48) 4.5.1 大涡的运动方程 被定义为亚格子尺度应力(subgrid-scale streese，简称SGS应力)，它体现了小尺度涡的运动对所求解的运动方程的影响。 (4-49) 4.5.1 大涡的运动方程 滤波前后比较 滤波后的Navier-Stokes方程与RANS方程在形式上非常类似，区别在于这里的变量是滤波后的值，仍为瞬时值，而非时均值，同时湍流应力的表达式不同； 滤波后的连续方程与时均化的连续方程相比，则没有变化，这是由于连续方程具有线性特征。 4.5.2 亚格子尺度模型 目的 亚格子尺度模型简称SGS模型，是关于SGS应 力 的表达式，建立该模型是为了使方程 (2-47)、(2-48)封闭。 根据Smagorinsky的基本SGS模型，假定SGS应力具有下面的形式: (4-50) 4.5.2 亚格子尺度模型 式中， 是亚格子尺度的湍动粘度，推荐用下式计算： 其中， ， ， 式中， 代表沿 轴方向的网格尺寸， 是Smagorinsky常数。 (4-51) (4-52) 4.5.2 亚格子尺度模型 理论上， 通过Kolmogorov常数 来计算， 即 。 当 =1.5时， =0.17。但实际应用表明， 应取一个更小的值，以减小SGS应力的扩散 影响。尤其是在近壁面处，该影响尤其明显。因此，Van Driest模型建议按下式调整 ： 式中， 是到壁面的最近距离， 是半经验常数，取25.0。 是Van Driest常数，取0.1。 (4-53) 4.5.3 LES控制方程的求解 通过式(4-50)将 用相关的滤波后的场变量表示后，方程(4-47)(实际对应3个动量方程)与(4-48)构成了封闭的方程组。 在该方程组中，共包含 、 、 和 4个未知量，而方程数目正好是4个，可利用CFD的各种方法进行求解。如使用有限元法求解，也可使用有限体积法得到各变量的值。 目前多数采用有限体积法求解，如SIMPLEC算法等。 4.5.3 LES控制方程的求解 关于LES的求解过程的补充说明 如果需要对能量或组分等进行计算，需要建立其他针对滤波后的标量型变量 的控制方程。 LES方法在某种程度上属于直接数值模拟(DNS)，在时间积分方案上，应该选择具有至少二阶精度的Crank-Nicolson半隐式方案或Adams-Bashforth方案，甚至是混合方案。在基于有限体积法的空间离散格式上，为克服假扩散，应选择具有至少二阶精度的离散格