第15讲-湍流及转捩3推荐.ppt

针对如下Sod 激波管问题 用5阶WENO格式计算其数值解，画出t=0.14时刻密度、速度及压力的分布；并与精确解进行比较（要求数值解与精确解画在同一张图上，便于比较）。 要求： 1） 空间网格数100， 时间推进格式选用3阶Runge-Kutta,时间步长自选。 2） 结合使用Steger-Warming 流通矢量分裂，画出结果 3） 结合使用特征投影分裂（Roe平均计算Uj+1/2) ，画出结果 * Copyright by Li Xinliang 习题8.1 使用WENO格式 求解激波管问题 针对如下Shu-Osher 激波-密度扰动波干扰问题： 用5阶WENO格式计算其数值解，画出t=0.1时刻密度、速度及压力的分布 要求： 1） 空间网格数200， 时间推进格式选用3阶Runge-Kutta,时间步长自选。 2） 结合使用Steger-Warming 流通矢量分裂，画出结果 3） 结合使用特征投影分裂（Roe平均计算Uj+1/2) ，画出结果 4） 使用2000个网格点计算，其结果作为“精确解”，与其它结果画在一起，便于比较。 * Copyright by Li Xinliang 初值为: 习题8.2 使用WENO格式 求解Shu-Osher问题 * Copyright by Li Xinliang 推导7阶精度的WENO格式，给出详细的推导过程及格式的具体表达式 提示： 与5阶WENO推导思路相同，但网格基架点扩大了，模板数目也增加了 j-3 j-2 j-1 j j+1 j+2 j+3 正通量(a0)的WENO通量 使用基架点{j-3,j-2,j-1,j,j+1,j,j+2,j+3} 如上图示，将其分割为4个组（模板），每个模板上4个基架点。 1） 先构建整个网格基（7个点）上7阶迎风格式的通量表达式 2） 对于每个模板，构造逼近j点导数的4阶差分格式的通量表达式 3） 按照理想权重进行组合可得到理想格式（7阶迎风格式） 对照具体表达式，求出理想权重Ck 4） 构建WENO通量的加权表达式 5）仿照 本PPT第35页的方法，给出光滑度量因子的具体表达式 可利用小程序求系数 coeff-schemes.f 减小工作量 习题8.3 推导7阶精度的WENO格式 Copyright by Li Xinliang * 习题 10.1 运用有限体积法求解翼型绕流 计算RAE2822超临界翼型绕流的流场，给出压力分布云图及翼型表面的压力分布。 流动参数： Mach 0.729, 攻角（AoA）= 2.31°， Re= 6.5×106 （基于弦长及来流值） 2.31° 建议方法： 采用本PPT介绍的有限体积法，具体步骤见本PPT 10.3节 Copyright by Li Xinliang * 习题 11.1 网格生成 通过解椭圆型方程生成NACA0012翼型的网格 要求： 推荐采用图示的C型网格， 网格点不限； 外边界的位置不限； 可求解无源项的Laplace方程或有源项的Poisson方程； 绘制出网格，并给出具体计算说明及公式。 “翼型数据库大全” /ph/www/m-selig/ads.html /bbs/viewthread.php?tid=46851extra=page%3D1 NACA 0012 翼型（对称翼型）的拟合曲线为 （宋宇宁等 “微型飞行器的翼型拟合与模具加工”，电加工与模具，2002第5期，33-36） 习题 12.1 求解方腔问题 问题描述： 如图示边长为L的方腔，上表面流体以常速度U运动，求解里面的流场（假设流动定常）。 考虑 三种情况 要求： 数值方法不限 （人工压缩性方法、投影法、涡量-流函数方法及SIMPLE方法均可）； 空间离散采用差分法，建议采用较高阶精度的方法。 绘制出定常解的流线图。 请详细写明方程及公式的推导过程及计算流程，切勿只上交计算结果。 * Copyright by Li Xinliang Copyr