湍流大涡数值模拟.ppt

湍流大涡数值模拟 四、复杂流动的大涡数值模拟算例 四、复杂流动的大涡数值模拟算例 一般认为有代表性的湍流流动： 湍流复杂程度 无分离的二维湍流边界层、直槽和圆管湍流、远场平面射流和尾流等。 简单湍流: 准平行的平衡剪切流动 复杂湍流: 简单湍流场中附加剪切或旋转 三维湍流边界层、大曲率的湍流边界层和射流、旋转系统的剪切湍流以及钝体分离湍流等。 横向喷流湍流， Kleisery，2006 四、复杂流动的大涡数值模拟算例 算例1 平面扩压器 平面扩压器的几何外型和网格 条件 雷诺数 入口截面的平均速度 ，入口半高 网格（352×64×128）,展向采用均匀网格；流向采用局部加密网格 展向采用傅里叶展开，流向和垂向平面上应用有限体积法离散。 标准Smagorinsky模型为基准的动力模式。 亚格子模型 四、复杂流动的大涡数值模拟算例 算例1 平面扩压器 计算结果 四、复杂流动的大涡数值模拟算例 算例1 平面扩压器 计算结果 四、复杂流动的大涡数值模拟算例 算例1 平面扩压器 计算结果 平均速度分布和实验结果符合良好。数值计算的分离点在 ，实验结果是 ，两者略有差别。计算的分离泡长度等于52δ ，实验结果是47δ 。 LES计算的峰值大于实验结果l0％～20％。雷诺应力分市的总体情况符合良好，误差在10％～20％。所有统计量在扩压器出口有较大误差。 结果还表明，LES能够捕捉流动特征，但是在计算的准确度上还有待改善。特别是扩压器二次流的影响。展向的周期性条件是否能准确计算实际扩压器的二次流值得研究。 四、复杂流动的大涡数值模拟算例 算例2 绕圆柱流动 几何外型和网格 条件 流动雷诺数 O型网格,圆柱体径向的计算域长度等于60倍的圆柱直径，(Kravchenko和Moin，2000) 分别采用Smagorinsky亚格子模型和动力亚格子模型。 亚格子模型 方法 一般有限体积法，对流项常用迎风格式，2阶守恒划中心差分格式或H—Spline迦辽金方法 四、复杂流动的大涡数值模拟算例 算例2 绕圆柱流动 网格影响结果 当雷诺数等于3900时，通过三组稀密程度不同网格的计算比较，考察计算结果对网格的依赖性。稀网格算例的网格点总数约50万，中等密网格的网格点总数约130万，密网格的网格点总数约240万。三种网格的计算结果几乎一样，说明中等密网格的结果可以采用。 (Kravchenko和Moin，2000) 数值格式影响结果 低阶迎风格式的数值耗散过大，不能获得正确的非定常结果，应当采用高精度格式。 (Mi.al和Moin，1997) (Kravchenko等，1999) 国防科技大学 * 李 桦 《湍流》8 大涡数值模拟的基本思想是直接计算大尺度脉动，而只对小尺度脉动做模式。所以，实现大涡数值模拟的第一步是把小尺度脉动过滤掉。本节首先介绍过滤方法，然后再导出大尺度运动的控制方程，并介绍小尺度脉动项的封闭模型。 大涡模拟的基本思想 二、大涡模拟的控制方程和亚格子应力 一、脉动的过滤 三、亚格子模型的检验 四、复杂流动的大涡数值模拟算例 五、关于大涡模拟的几个问题 湍流大涡数值模拟 内容： 一、脉动的过滤 一、脉动的过滤 湍流速度场 一、脉动的过滤 常用的均匀过滤器 ? 谱空间低通滤波 ? 物理空间的盒式滤波器 ? 高斯过滤器 若 k 为湍流脉动波数，物理空间的湍流脉动在谱空间的投影为 ，则谱空间过滤后的脉动可表示为： 为各向同性低通滤波 其中 是截断函数， 　是截断波数。 一维低通滤波器 谱空间的过滤就是令高波数的脉动等于零，相当于对脉动信号做低通滤波。 一、脉动的过滤 常用的均匀过滤器 ? 谱空间低通滤波 ? 物理空间的盒式滤波器 ? 高斯过滤器 在尺度 l 上进行的滤波函数记作 ，则任意湍流脉动 的过滤为 盒式过滤器 对于复杂流动不可能采用谱空间低通滤波，物理空间的低通过滤可以用积分方法实现。 物理空间的滤波器必须满足正则条件 是过滤的空间体积。 一、脉动的过滤 常用的均匀过滤器 ? 谱空间低通滤波 ? 物理空间的盒式滤波器 ? 高斯过滤器 高斯过滤器 将过滤函数取作高斯函数。 ? 只有高斯滤波器在Fourier积分变换时保持高斯函数形式。 应当注意： ? 物理空间的盒式过滤，在 的高波数区有微小的泄漏(Gibbs效应)。 ? 谱空间的盒式过滤器变换到物理空间时并不是盒式过滤器，而在 的盒子以外，过滤函数仍有微小的震荡。 一、脉动的过滤 脉动速度的分解 经过