用确定性涡方法模拟二维钝体绕流.pdfVIP

第 26卷 第 4期 空 气 动 力 学 学 报 Vo1．26。No．4 2OO8年 l2月 ACTA AERoDYNAM ICA SINICA Dec．，2OO8 文章编号：O258一l825(2OO8)O4·O526一O6 用确定性涡方法模拟二维钝体绕流 杨素哲，周志勇，陈艾荣 (同济大学桥梁工程系，上海 2O0092) 摘 要：以离散涡方法为基础，采用粒子强度交换法模拟涡量扩散方程，对粒子场分布产生扭 曲变形时采用粒子强 度重分配，数值模拟了突然起动圆柱、风攻角为的方形截面和大海带东桥的流场，计算结果与解析解及其他数值方 法计算的结果均吻合的很好，数值结果显示确定性涡方法的高精度的特性使其在模拟钝体扰流问题方面具有很大 的优势 。 关键词 ：离散涡方法 ；粒子强度交换法；粒子强度重分配 ；钝体 ；三分力系数 ；斯托洛哈数 中图分类号；V2l1．3 文献标识码 ：A 强度交换法 尚未应用到工程中，且 国内对这种方法的 0 引 言 研究比较少，对其在工程上的应用更是空白。本文对 粒子强度交换法进行了研究并程序化，将其初步应用 涡方法是模拟非定常流体的拉格朗 日方法，可以 到桥梁截面的静风分析中。 忽略离散误差 ；涡方法具有很好的守恒特性。这两个 特点使得涡方法成为模拟流场 的有效的工具。二维 1 确定性涡方法 涡方法发展最成熟 的是 由Chorin[1]提出的随机走步 涡团法，它是用随机走步来模拟涡量的粘性扩散，具 1．1 涡方法 有统计平均意义，该方法已经成功地运用到了桥梁断 涡量场的二维 Navier-Stokes方程可以汇总为如 面的流场模拟和抗风分析 中，如 AllanLarsen[2和 下数学模型 Zhou[3]。但是它 的缺点不能忽视 ，对长时间历程 问 +l1． ∞一 v 。ccJ inD (1a) 题 ，误差急剧增长 ，随机走步模拟扩散精度较低 ，且用 新生涡团处理无滑移边界条件有较高的数值噪声，确 V =一 inD (1b) 定性涡方法克服了随机涡方法的缺点，而不改变涡方 lf— × inD (1c) 法作为Lagrange方法的优点。过去的二十年里，在 (z，O)= ∞0() inD (1d) 涡方法领域内，已经取得 了较大的发展 ：快速多极子 lI— U (￡)+ n(￡)P：× (，一z6) onaD (1e) 展开 的应用 (Greedgard和 Rohklin4【]，Barnes和 at。。 (1f) Hut[5])用粒子强度交换法精确地模拟粘性扩散(De— 其中∞为涡量 ID分别为流体的运动粘性系数 gond和Mas—Ga1lic[6])，以及精确施加粘性边界条件 和流体 的密度 ， 为流函数，固体 的平动速度为 (K0um0utsakos，etal[ ，Mas—Gallic[8。，Leonardet (￡)，绕其形心 的角速度为n(￡)，U 是无穷远处 al[，Ploumhans，etal[0。，Benhaddouch[1)。所有这 的速度 ，式(1b)为泊松方程 ，(1c)为不可压缩条件 ， 些使得涡方法能够更高精度地模拟具