外流问题的外边界条件设计Ξ.PDF

　A 辑第 17 卷第 6 期 Ser. A , Vol. 17 ,No. 6 　　　　　水 动 力 学 研 究 与 进 展 　　　 　　2002 年 12 月 　　　　　　　 　　　Dec. , 2002 J OURNAL OF HYDROD YNAM ICS 　　 文章编号 (2002) 外流问题的外边界条件设计 李勇 , 　　李向群 , 　　翟大熙 (北京大学力学与工程科学系 ,北京 100087) 　　摘 　要 : 　粘性外流计算原则上要求在无限大的外空间进行 , 巨大的计算区要求巨大的计算 机资源。本文利用一般绕流问题在远区流动尺度大 ,因而雷诺数大的特点对远外处的流动近似地 用位势流代替 ,从而在较小的流体边界条件上建立了“开边界条件”, 由此大大缩小了计算工作量。 关 　键 　词 : 　NS 方程 ;开边界条件 ;势流理论 中图分类号: 　O357 　　　文献标识码 :A 1 　引言 外流场计算是数值模拟物体在空气或其他流体中运动时流体本身的运动。由它可以求得 两者之间相互作用 ,从而得知物体最后的运动。粘性流体力学方程是非线性的 ,要解它原则上 必须在无限的流体外空间布置计算 ,事实上这是做不到的。在实践上计算只能到“足够大”的 流体空间上止步 ,而将本该设在无限远处的条件挪到此空间边缘。为了计算结果不致偏离太 远 ,常常要求这个空间相当大。以机翼的外流计算为例 ,一般要求计算空间的宽度至少大于 10 倍弦长。计算区域如此大 ,而有兴趣的却仅仅只有机翼附近小范围的流体运动。于是人们 设计外边界条件 ,要求计算区域缩小而机翼附近的流动“依旧”。这样的条件常称作“开边界条 件”,最常用的开边界条件就是令参数顺边界法线的微商为零。这样的开边界条件比较简单 , 用得很多 ,但能节约的计算面积不多。本文提出在粘性流计算区外边界上利用场论建立流动 参数的边条件来压缩计算区。实践表明这可以将计算区缩到很小而又不致明显影响物体附近 的流动解。这样的开边界条件的实践意义是不可言而喻的。 对于亚声速粘性流动来说 ,随着与物体的距离增加流体流动参数的变化变缓 ,亦即空间尺 度变大 ,流动呈现“大雷诺数”性态。根据这一特点 ,在远处我们可用理想流体运动来刻画。起 源于均匀流的理想流体运动是位势流 ,场论的公式可以应用。本文对粘性区开边界条件的设 计原则便是假定粘性计算区外是位势流 ,根据位势解的公式建立粘性流体计算边界上流动参 数之间的关系 ,这就是我们设计的开边界条件。我们虽然外部流动做了简化 ,但基本的流动法 则依旧遵守 ,所以它的结果比其它的开边界条件设计要好 ,在物体周围的流动计算结果影响不 大的条件下容许采用更小的计算区。 收稿日期 : ( ) 作者简介 : 　李勇 1974 ,男 ,MIT 访问学者 ,硕士。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 710 水 　动 　力 　学 　研 　究 　与 　进 　展 　　　　　　　　　2002 年第 6 期 2 　问题的提法及解题思路[1 ] 在流体力学的计算中 ,有很多属于无界区域中的粘性绕流问题。等密度常粘性流体的绕 流问题的提法为 : 控制方程 v 1 2 ( ) t + v · v = - p + R e v 1