若干预处理方法在低速粘性流动中的应用比较.pdfVIP

江　苏　航　空 2010增刊　 若 干 预 处 理 方 法 在 低 速 粘 性 流 动 中 的 应 用 比 较 胡丽燕　陈红全 (南京航空航天大学航空宇航学院, 南京, 2 10016) 摘 要: 选 用W e iss-Sm ith, Ch o i and M erk le 和 P le tch er and Ch en 三种 预处 理方 法, 成功 地使 求 解可 压 N av ier- Stok e s方程 的有限体积方法能用于低速粘性流动 问题 。文 中给出了圆柱等低雷诺数粘性流动 的数值结果, 整体 上能与文 献结果相吻合, 文 中同时对三种 预处理结果进行了分析与 比较, 展示出不同预处理方法 所对应的收 敛 历程及预处理效果。 关键词: 预处理; 低速; 粘性; W eiss-Sm ith; Ch o i an d M erk le; P le t ch e r and C h en (Qv )具体表达式参见文献 [ 1 ], 其 中 Γ是预处理矩 引　　言 阵, 可 以有多种选择, 本文采用 - 预处 W eiss Sm ith 在流体力学中, 对于包含低速、亚声速 、跨声速 理矩阵[2 ], Ch oi an d M erk le预处理矩阵 [1 ], P letch - 和超声速 的全速度流场, 以马赫数为衡量尺度, 把 er an d Ch en 预处理矩阵 [3 ]具体表达式如下: 流体分为不可压流和可压流, 一般来说, 在 M a 1 0 0 0 0. 3 时, 认为是不可压, M a 在 0. 3 以上则属于可压 Xad c 缩流。描述可压和不可压流动的控制方程, 从数学 Γw = 0 1 0 0 ( 2) 性质上来讲有很大 的不 同, 可压流控制方程是抛物 0 0 1 0 2 -双 曲型的, 不可压流动控制方程是椭圆型的, 由此 - c 0 0 1 两者的数值求解算法有很大差异 。 1 0 0 0 采用解初值 问题的时间推进方法—— 时间相 0 1 0 0 关法求解问题时, 马赫数大于 0. 3时压缩性较为 明 Γc = 0 0 1 0 ( 3) 显, 可 以得到较为满意的结果, 随着马赫数趋近于 0 0 0 1