CFD经验漫谈+湍流初步.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 北京航空航天大学 能源与动力工程学院 School of Propulsion, BUAA 北京航空航天大学 能源与动力工程学院 School of Propulsion, BUAA CFD经验漫谈及湍流模化初步 CFD的基本流程 前处理 建模 Autocad, CATIA, UG, Solidworks 网格生成 Gridgen, ICEM-CFD, Gridpro 求解器 设置计算参数 后处理 Tecplot Ensight 网格类型 结构网格 多块对接网格 多块重叠网格 非结构网格 笛卡尔网格 网格划分原则 1.第一层网格距壁面的距离 层流计算:y(1)~1/sqrt(Re) 湍流计算:y+1 壁函数 30-200 2.网格法向间距拉伸率 工程应用:r1.25~1.3 科学研究:r1.2，ideal: r1.1 拉伸不是无休止的，deltay到边界层厚度/10停止(大约位于边界层中心) 3.流向网格间距 RANS计算: 找感觉。个人经验，Re=10e.6, 单位长的壁面附近布置80+点足够 Grid independent Study:粗中细三套网格，比较结果(特别是摩阻) LES计算： 理想LES三方向的网格间距应相当 法向与RANS相当，流向deltax+, deltaz+~8000(Spalart)，但不得超过边界层厚度/20 Grid Refinement Study 4. 网格品质要求 正交性 物面网格线垂直于壁面 空间网格线正交 软件带有网格正交化功能 光滑性 疏密过渡要平缓 考虑流场特性 流动变化剧烈的区域，如激波、分离区域局部加密 Examples Grid independent Study y+(1)大小的影响 求解器 求解参数 时间步长 隐式求解LUSGS, LU-ADI: CFL=5~10 显式求解Runge-Kutta: Dt很小 情况变化(Mach提高，高级湍流模型)，收敛困难 手工减小时间步长，计算稳定后逐渐提高 构造函数，使时间步长自动缓慢增长 设置合理初场 使用当地时间步技术 限制器 无激波，无需限制器，否则会引入不必要的粘性 出现激波，限制器顺序min mod(高Mach), Smooth, Van Albada, double minmod (低Mach) 推荐应用多重网格，极大加速收敛，并促进计算稳定 无粘通量空间计算方法 Roe AUSM 非定常计算 Dtdelta/Umax NaN错误 非物理解修正 Stack Overflow错误 数组别都放栈里，建议用动态数组 建议 学习一个开源代码，希望编程模块化，输入输出接口标准化，参数设置简单化 计算环境使用Linux 编译器使用intel fortran/C++ 如有可能，尽量抛弃Windows 使用Linux，办公效率提高几十倍！ 对新生寄予厚望 什么是湍流？ Leonardo da Vinci (1452-1519) 在水中放置障碍物后观察流动现象。 湍流：极不规则的流动现象。 葛饰北斋 (1760—1849)：神奈川巨浪 歌川広重 (1797—1858)：鸣门的风浪 艺术家笔下的湍流 云层 自然界的湍流 河流 香烟的烟雾 海浪 减阻 控制转捩：减少阻力 燃烧 强化湍流：加强扩散 噪声 抑制湍流：降噪 气象 预测大气湍流 湍流的研究意义 目的：预测和控制湍流 雷诺实验 湍流的形成 雷诺数： 层流 湍流 O. Reynolds (1842-1912) 层流向湍流的发展过程 阶段一：当Re很小时，流体规则地流过圆柱 阶段二： 随着Re增加，尾涡不断增大 阶段三：经过Re阀值后，变为极不规则的运动 雷诺数： 湍流的特征 特征一：不规则性 射流 特征二：湍流脉动的尺度具有宽广的范围 特征三：拟序结构。 剪切湍流存在拟序结构 湍流脉动并非完全不规则的随机过程 有序大尺度结构包含于不规则脉动中 F. N. M. Brown, A. Roshko, 1974 湍流的分类 均匀各向同性湍流 剪切湍流： 壁湍流：槽道流，圆管流，边界层湍流 自由剪切湍流：射流，混合层，远场尾流 平均剪切 固壁 均匀和各向同性 平移几何构形 平移加转动几何构形 均匀：任意n点的空间几何构形在空间平移时，脉动量的n阶相关函数不变。 均匀各向同性：任意n点的空间几何构形在空间平移和刚体转动时，脉动量的n阶相关函数不变。 均匀：脉动量的任意阶统计相关函数和坐标系的平移无关。 均匀各向同性：脉动量的任意阶统计