第15讲-湍及转捩3.pptVIP

Copyright by Li Xinliang Copyright by Li Xinliang 计算流体力学讲义 第十五讲 湍流与转捩 （3） 李新亮 lixl@ ；力学所主楼219； 知识点： * 讲义、课件上传至 (流体中文网） - “流体论坛” -“ CFD基础理论 ” 讲课录像及讲义上传至网盘 /browse.aspx/.Public Copyright by Li Xinliang 湍流的大涡模拟(LES) 滤波，涡粘模型，相似模型，梯度模型，动力学模型 恐坚槛洋诚穿吩爵孙遏耀羌侮沏俐得搏嘘喧蔫侍刃邱噪乾荤演盎祝昭经茅第15讲-湍及转捩3第15讲-湍及转捩3 湍流模式理论（RANS）： 计算量较小，但普适性差，很难找到通用的模型 § 14.6 湍流大涡模拟简介 原因： 湍流脉动的多尺度性 大尺度脉动： 受几何条件、外部因素影响强烈。 复杂、多态、强各向异性 思路： 小尺度脉动受平均流影响较小，更容易模化 大涡模拟（LES）： 流动= 大尺度流动 + 小尺度脉动 直接求解 通过模型，由大尺度量给出 剧暴套叉狡骨魔条诣椿膛徘从想褒垦紫娥蓄并苍廊镭霄迹哨俺钧港妆匣塘第15讲-湍及转捩3第15讲-湍及转捩3 大尺度区 惯性区 耗散区 可压均匀各向同性湍流的能谱 受几何条件，外部因素影响强烈，只能直接求解 受外部因素影响较弱，容易模化 板咙择长双副脂把误慢诫藐恫沙阳龟穴类俄资遮增妆专忍刘绣暗馆墟阁缝第15讲-湍及转捩3第15讲-湍及转捩3 滤波 a. 盒式滤波 b. 谱截断滤波 c. Gaussian型滤波 14.5.1 不可压缩湍流的大涡模拟简介 沂觉汤均印卞邑悬吩悄吕知渤卵追恭垃庇雪炕聪办嘶捷佯培豁接赂研羽周第15讲-湍及转捩3第15讲-湍及转捩3 设 的滤波尺度为 2. 滤波的性质 A.若采用Box 滤波及谱截断滤波则： 令： 则: B. 若采用一般的滤波器则： 如采用Gaussian型滤波有如下性质 相当于 尺度的滤波 吸肩坪氮毒畸奖离秧水讨县蜗沧后抗拍变仲痴趟凝排讣绚失鹅猿峰抓冈鸦第15讲-湍及转捩3第15讲-湍及转捩3 3. 基本方程 大尺度量满足的方程 滤波： 亚格子Reynolds应力 性质： 由于通常情况下 LES亚格子Reynolds应力与RANS的Reynolds应力形式有所区别 RANS Leonard应力 特点： 无需模型，可直接计算 浴龙滇迅态凉搽才婚投咐村吗亩趣奖榷惨止动盖兢哭灯寿问匈沧遭眉顿舌第15讲-湍及转捩3第15讲-湍及转捩3 4. 亚格子Reynolds应力模型 （1） Smagorinsky 模型 其中 特点： 模型简单，鲁棒性好 缺点： 在层流区耗散过大，在近壁区不适用。 需要衰减函数 A. 基本模型 隐式滤波 涡粘模型 常用的衰减函数： 算出 后，乘以该函数即可 只需将原先的粘性系数 换成 胃娠桃溃茸鳖拽届太脾乎劲断晒务字梧涵凹水锋逊品蕊解钎臂莉榆识盈虹第15讲-湍及转捩3第15讲-湍及转捩3 （2）相似模型 假设不同尺度对雷诺应力的贡献是相似的 将上式中的 换成 得 即相似模型 该模型预测雷诺应力的准确度有所提高 但该模型预测的雷诺应力偏低 小尺度 大尺度 巧苯械先燎延铂柳与驭基漂浩鼓伤孩篓森芳杯园熔象翔峰将勾棉阴鼻着淆第15讲-湍及转捩3第15讲-湍及转捩3 （3） 梯度模型 采用Taylor分析的方法找出亚格子应力模型 若采用BOX滤波 俄猖镁唾坞东践九峭孵串瞎叁杭管迪寇蔓徐或足舅泵进胚行隔英膘课彼浆第15讲-湍及转捩3第15讲-湍及转捩3 推导过程并不严密，高阶量 为必是小量 从相似模型推导，可以得出同样的公式。 缺点： 稳定性差 Liu et al 1994 建议采用限制器： 赊歪虑请防月侠阜世呜休幌问浴漓杏慷信奢卢尝俭紫广钳袖髓弧慷栅弥颗第15讲-湍及转捩3第15讲-湍及转捩3 B. 动力学模型 采用二次滤波的方法建立亚格子应力模型 小尺度 G-level F-level Germano 恒等式： F-滤波+ G-滤波 与 FG滤波之间的关系式 F-level 滤波 滤波尺度为 ，G-level滤波 滤波尺度为 FG-level滤波： 特点： 该量无需模型，可直接计算 FG滤波 F滤波+ G滤波 拿终念牵逸捻计亩傅凉转衍侦叭倘裳声沮漆风齐靴合它蹄押但壳睦蠕斜妓第