高雷诺数方腔流动的分子模拟.pdfVIP

第 45 卷 第 5 期 力 学 学 报 Vol. 45 , No. 5 2013 年 9 月 Chinese Journal of Theoretical and Applied Mechanics Sep. , 2013 研 究 论 文 1 高 雷 诺 数 方 腔 流 动 的 分 子 模 拟 2 费 飞 樊 菁 中 国 科 学 院 力 学 研 究 所 , 北 京 100190 2 4 摘 要 采 用 扩 散 信 息 保 存 di usive information preservation, D-IP 方 法 计 算 了 雷 诺 数 为 10 10 的 二 维 方 腔 流 动D-IP 网 和 长 步 间 时 于 对 法 方 罗 卡 特 蒙 拟 模 接 直 典 经 了 服 克 , 法 方 拟 模 子 分 的 点 观 动 运 散 扩 于 基 种 一 是 法 方 制 限 格 严 的 小 大 格, 中 算 计 在 D-IP 十 几 的 程 由 自 均 平 和 间 时 撞 碰 均 平 子 分 为 别 分 小 大 格 网 和 长 步 间 时 的 法 方 均与 所得到的方腔流线分布和旋涡的精细结构, 倍乃至几百倍, Navier--Stokes 方程数值解相符关 键 词 扩散信息保存方法, 方腔流动, DSMC 网格大小 时间步长, 方法, 中 图 分 类 号 : V211.25 文 献 标 识 码 : A DOI : 10.6052/0459-1879-13-064 引 言 于 分 子 平 均 碰 撞 时 间 、 网 格 尺 度 小 于 分 子 平 均 自 由 程 , 这 种 做 法 在 物 理 上 才 是 合 理 的以 DSMC 计 算 数 值 模 拟 方 法 可 以 大 致 分 为 两 类 : 连 续 介 质 方 [9] 的 黏 性 系 数 相 对 误 差 为 例 , Alexander 等 法 和 粒 子 方 法鉴 于 直 接 模 拟 蒙 特 卡 罗 direct simu- [10] 与 Hadjiconstantinou 的 分 析 表 明 : 它 与 时 间 步 长 lation Monte Carlo, DSMC 方 法 在 稀 薄 气 体 中 的 巨 大 t 和 网 格 大 小 x 的 平 方 成 正 比当 t 和 x 分 别 [1] 成 功 , 将 其 推 广 应 用 于 更 广 泛 的 领 域 , 是 容 易 想 与 分 子 的 平 均 碰 撞 时 间 和 平 均 自 由 程 相 当 时 , 到 , 也 很 能 引 起 兴 趣 的 工 作 , 当 前 颇 为 活 跃稀 薄 气 分 别 约 为 5% 和 10%. 在 标 准 大 气 条 件 下 , 空 气 分 体 流 动 和 连 续 介 质 流 动 通 常 以 努 森 数 Kn 0.01 分 界 8 子 的 平 均 自 由 程 和 平 均 碰 撞 时 间 分 别 约 为 6 10 m 随 着 Kn 进 一 步 增 大 或 减 小 , 都 将 有 过 渡 区 出 现 , 前 10 和 2 10 s. 由 于 Kn 反 比 于 Re 和 Ra , Re 或 Ra 很 者 通 向 自 由 分 子 流 , 后 者 通 向 湍 流在 通 向 自 由 分 子 大 , 意 味 着 Kn 很 小 , 即 流 的 动 特 征 尺 度 远 大 于 分 子 流 的 过 渡 区 , 分 子 拟 模 法 方 表 现 出 色 ; 在 通 向 湍 流 的 运 动 特 征 尺 度 , 即 使 用 当 前 最 快 的 计 算 机 , 也 很 难 转 捩 区 , 分 子 模 拟 方 法 是 否 也 能 有 所 作 为 呢 ? 这 个 满 足 DSMC 这 在 种 情 况 下 的 时 间 步 长 和 网 格 大 小 要 问 题 吸 引 了 许 多 研 究 者 , 包 括 DSMC 方 法 的 发 明 者 求[2-3] [4-5] Bird 和 分 子 动 力 学 的 发 明 者 Alder 等一 些 研 究 者 采 用 扩 散 观 点 , 将 分 子 运 动 和 碰 撞 分 子 模 拟 方 法 用 于 流 动 稳 定 性 , 从 原 理 上 讲 没 [11-12] 耦 合 在 一 起 考 虑这 种 做 法 从 源 头 上 消 除 有 问 题 , 在 无 量 纲 控 制 参 数 相 对 小 的 情 况 下 , 已 经 DSMC 对 于 时 间 步 长 和 网 格 大 小 的 苛 刻 限 制 , 其 实 取 得 一 些 新 的 认 识 , 例 如 Rayleigh--Benard 问 题 中 , 现 的 关 键 在 于 如 何 实 时 获 得 每 个 分