06-湍流模式和大涡模拟理论汇编.ppt

流体力学讲义 第六讲 湍流模式和大涡模拟理论 张凌新 浙江大学航空航天学院 流体工程研究所 Institute of Fluid Engineering 内容 雷诺平均统计模式 建立统计模式的一般原理 湍流涡粘模式 代数模式 标准双方程模式 非线性双方程模式 壁函数 低雷诺数修正 二阶矩模式 湍流大涡模拟 脉动的过滤 控制方程和亚格子应力 控制方程 亚格子模型 亚格子模型的检验 先验 后验 建立统计模式的一般原理（1） 建立统计模式的一般原理（2） 建立统计模式的一般原理（3） 湍流的涡粘模式（1） 湍流的涡粘模式（2） 湍流的涡粘模式（3） 湍流的涡粘模式（4） 湍流的涡粘模式（5） 湍流的涡粘模式（6） 大涡模拟-脉动的过滤（1） 大涡模拟-脉动的过滤（2） 控制方程和亚格子应力（1） 控制方程和亚格子应力（2） 控制方程和亚格子应力（3） 亚格子模型的检验（1） 亚格子模型的检验（2） 亚格子模型的检验（3） IFE , Zhejiang University 浙江大学航空航天学院流体工程研究所 * 高等流体力学 尺度问题 c l c l l 对应的是含能尺度 统计模式是对尺度l 以下尺度的脉动进行模化《 模式思想 统计模式的思想是建立高阶矩与低阶矩之间的关系式 有两种方法来对雷诺应力进行模化 一种：是建立二阶矩（雷诺应力）与一阶矩（平均速度）之间的关系 另一种：是建立三阶矩与二阶矩之间的关系《 近似准则1-涡粘模式 脉动场与平均速度场、脉动速度的初始条件和边界条件有关《 脉动速度的一般形式 雷诺应力的一般形式 从统计的角度，初始扰动的影响可以排除 而边界的影响则很大《 这是二阶矩由一阶矩近似的准则《 近似准则2-二阶矩模式 这是三阶矩、二阶矩由速度二阶矩近似的准则《 涡粘模式 S为流体平均速度的应变率张量 湍流粘度就是模式理论研究的关键 实际上更为一般的式子为 代数涡粘模式 普朗特混合长度理论 Smagorinsky模型 三维湍流粘度公式《 Baldwin-Lomax模型 采用涡量来计算湍流粘度《 分内外层来分别计算湍流粘度《 普朗特湍流粘度公式《 湍动能方程 标准双方程涡粘模式 c l 标准双方程涡粘模式 耗散率方程最难模化 耗散率模化方程是比照湍动能模化方程而来《 耗散率方程 标准双方程涡粘模式-系数的确定 雷诺应力的生成项与耗散项相平衡《 雷诺应力的表达式《 测量确定《 在湍流边界层的等应力区1 标准双方程涡粘模式-系数的确定 各向同性湍流只有耗散《 对于早期衰变，n=1.3《 或 对各向同性湍流 重整化群双方程涡粘模式 RNG模式 非线性双方程涡粘模式 应变率、涡量的二次项 历史效应《 后台阶流动的比较 标准双方程模式 非线性双方程模式 壁函数 粘性底层 对数律层 等应力层 粘性底层 对数律层 低雷诺数修正 标准模式在粘性底层过高估计了湍流粘度《 c l c l 过滤的尺度 大涡模拟是对l c以下的尺度脉动进行模化，从而计算大尺度脉动《 脉动的过滤-谱空间的过滤 低通滤波器《 脉动的过滤-物理空间的过滤 滤波函数《 低通脉动和剩余脉动《 滤波函数与低通滤波器《 平均的再平均不等于前面的平均《 亚格子应力 控制方程 改写 控制方程与雷诺平均方程相似《 为什么不写成雷诺应力的形式《 亚格子模型-Smagorinsky涡粘模型 系数Cs的确定 或 与RANS模型相似，可以很好与雷诺统计模式兼容 与RANS模型相似，也存在涡粘模式的固有缺点《 并未用到脉动的过滤 亚格子模型-尺度相似模型 亚格子应力 二次过滤的亚格子应力 尺度相似模型 耗散太小《 混合模型《 亚格子模型-动力模型 细网格上过滤的亚格子应力《 粗网格上过滤的亚格子应力《 两次过滤的亚格子应力《 其中 M的表达式利用了Smagorinsky模型《 L是已知的《 动力模型就是基于其他模型来确定Cs《 先验比较 DNS的结果进行过滤，计算的亚格子应力与LES的模型进行比较《 Smagorinsky模型的亚格子应力相关系数0.25， 尺度相似模型和混合模型相关系数0.6 先验比较并不需要LES的计算结果《 逐点比较没有意义，进行相关性的比较 但与此同时Smagorinsky模型的亚格子耗散的相关系数0.5-0.7， 而尺度相似模型的相关性很差 后验比较 用LES的计算结果与实验或者DNS的进行比较《 M1-Smagorinsky模型 M2-尺度相似模型 M3-基于Smagorinsky的动力模型 M4-基于混合模型的动力模型 混合层的计算结果比较《 IFE , Zhejiang University 浙江大学航空航天学院流体工程研究所 * * * *