离散空间直接建模的计算流体力学方法
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
520
描述流体运动 , 对应最小尺度则是分子的平均自由 程 . 在时间方向也存在类似的尺度问题 , 这里不再 赘述. 需要注意的是, 当物理模型由微分方程来描述 后, 数学意义上的最小可分辨尺度变成了无限小(零 ). 如果采用离散模型 , 则最小可分辨尺度对应于离散 尺度. 换句话说, 物理模型表现出的性质会受到描述 方式的限制. 在 CFD 计算中, 数值离散本身将引入一个新的 最小可分辨尺度, 即网格尺度∆x. 小于网格尺度的流 动变化同样只能通过(数值 )模型来表达 . 最小尺度的 差异决定了数值解、 物理解和偏微分方程的解析ຫໍສະໝຸດ Baidu之 间存在差异. 为了准确捕捉最小尺度为 l 的流场脉动, 网格尺度必须满足∆x<l. 对直接离散 PDE 而建立的 传统 CFD 方法, 还必须要求 PDE 隐含的物理模型最 小尺度 lPDEl. 对于多尺度流动问题 , 一个优秀的 CFD 方法应该能按实际需要自动而高效地刻画出相 应尺度上的流动特征. 也就是说, 如果实际应用要求 给出小尺度上的流动, 在满足要求的网格尺度上, 数 值方法能捕捉到该小尺度上的物理量 ; 如果只要求 给出较大尺度上的流动特征 , 数值方法也能在相应 较粗的网格上给出合理的流场物理量. 当然, 一个比 较直观的想法是以流动的最小尺度来布置计算网格 , 从而能将大大小小的各种不同尺度的流场脉动同时 捕捉出来 . 不过 , 如高雷诺数湍流的直接数值模拟 (DNS)一样, 这种方法可能需要极其巨大的计算资源, 甚至根本不可行. 实际上, 科学研究的目的是对自然 的最有效的描述 , 如果所有真实的描述都要归结为 最细致的描述 , 那么世界上除了粒子物理之外就没 有其他科学了. 上面提到的采用不同尺度网格计算得到的流场 物理量之间的关系中, “合理”是指采用粗网格得到的 单元内物理量应该与采用细网格得到的流场按粗网 格尺度平均后得到的平均量近似相同. 同时, 为捕捉 边界条件等大尺度因素的影响 , 网格的尺度必须足 够小. 因此, 随着计算网格的逐步加密, CFD 计算得 到的流场平均量应该趋于定值 , 这就是网格收敛性 . 需要注意的是, 这里的网格收敛性与传统 CFD 计算 中的提法不同. 后者通常是指随着网格加密, 数值解 逐渐趋近于 PDE 的解析解; 而前者是指数值解与物 理解的逼近程度随网格尺度的变化 , 包含了不同尺 度下 PDE 本身对真实流动的逼近程度. 比如, 如果 网格大小和分子直径相近, CFD 应该可以看到独立的
摘要
本文提出了一种发展适合多尺度、 多物理流动的 CFD 格式的全新途径, 即直接在离散空间利用
物理模型的跨尺度演化解来建立数值格式的新方法. 与直接离散偏微分方程的传统做法相比, 基于离散 空间直接建模构造出的数值方法考虑了网格尺度和物理模型之间的匹配, 能做到不同尺度上物理模型的 连续过渡, 从而实现对多尺度流动的高效模拟. Boltzmann 方程的模型演化解为直接建模方法提供了重要 的支撑. 本文提出了直接建模构造计算流体力学方法的理论基础, 并给出了在建立统一气体动理学格式 中的成功运用, 验证了新方法的可行性和优越性. 关键词 偏微分方程数值解, 离散空间直接建模, 多尺度流动, 统一气体动理学格式, BGK 方程
徐昆等: 离散空间直接建模的计算流体力学方法
方程基于气体分子(粒子 )混沌假设, 其模型尺度是分 子平均碰撞时间和平均自由程. 虽然 Boltzmann 方程 理论上可以适用于从连续到稀薄的跨流域流动[6], 但 其要求在任何地方都能分辨出分子的微观运动尺度 , 也就是所谓的动理学尺度 (Kinetic Scale), 所以直接 把 Boltzmann 方程用于宏观流体的计算是不现实的. 这些不同的模型方程实际上是在不同尺度上描述流 体运动(下文将对此展开讨论), 因而逼近真实流动的 近似程度不同 . 在流体计算中由于网格尺度的加入 , 使得问题更加复杂 . 虽然通常认为当网格尺度足够 小时, 直接离散 PDE 得到的数值解能以任意的精度 逼近 PDE 解析解, 但它对真实物理流动的描述还受 限于 PDE 模型本身. 因此, 数值解逼近真实流动的 精度不光与计算网格有关, 而且与所选用的 PDE 模 型有关. 对于跨流域流动这种多尺度流动问题, 不同 流域需要选取不同的 PDE, 在过渡区, 甚至找不到一 个合适的 PDE 模型, 从而导致数值计算的难度太大. 一个很有意思的问题是 , 能否直接在离散空间上建 立流动的(数值 )物理模型, 从而通过数值解来直接逼 近物理解? 这也正是本文将要重点探讨的内容. 本文以多尺度流动问题为例 , 首先分析了不同 尺度下描述流动的物理模型 , 回顾了传统的基于直 接离散 PDE 的 CFD 方法, 进而提出发展 CFD 方法的 一条全新途径, 即在离散空间直接建模的方法. 以近 年来发展起来的适用于跨流域流动的统一气体动理 学格式 (UGKS) 为例 , 对新思路进行了验证和分析 . 另外, 当前传统的基于 Riemann 解的格式在高超声速 流动计算中遇到了瓶颈 , 现在是近一步思考和发展 下一代数值计算方法的关键时代 . 本文最后对此进 行了一定的讨论.
2
流动的描述尺度与物理模型
对真实流动的描述首先应该基于一定的最小可 分辨尺度 , 不同描述尺度下能看到不同的流动细节 , 以此可以建立不同的流动物理模型 , 通常为偏微分 方程. 事实上, 工程应用中关心的也是在某个有限分 辨率下的物理量, 如速度、温度和压力等, 或者说是 具有最小可分辨尺度的微元内的平均值 . 小于最小 尺度的流动变化则反映在物理模型中 , 如宏观方程 中耗散与内能的关系. 也就是说, 不同 PDE 模型隐 含的最小尺度 lPDE 不同. NS 方程基于连续介质假设, 对应宏观耗散尺度, 而 Boltzmann 方程从介观角度
3
基于直接离散 PDE 的 CFD 方法
传统 CFD 方法是通过直接离散描述流体运动的 偏微分方程建立起来的 . 借助偏微分方程数值解理 论 , 对简单流动 , 特别是具有线性变化特征的流动 , 可以比较容易地对 CFD 方法的基本性质, 如稳定性 和收敛性等进行理论分析. 实际计算中, 可以根据需 要通过加密网格 , 从而更好地逼近 PDE 的解析解 . 此时 , 数值解对物理解的逼近问题实际上是要解的 PDE 对真实流动的近似问题. 这里的理念(原理)实际 和计算关系不大 , 更像牛顿发明微积分用小单元取 极限描述连续曲线一样 . 但实际的计算一定是在有 限网格大小上进行的, 所以传统的 CFD 方法得到的 只是相应 PDE 的一个近似解, 而且在理论上通常也 不可能知道在数值上求解的确切方程. 对于多尺度流动问题 , 上述方法存在的不足就 更加明显 . 其一 , 对于非线性问题 , 基于直接离散 PDE 的数值格式的理论分析非常困难 , 至少到目前 为止 , 仍亟需分析方法的突破 . 其二 , 对于流动的不 同尺度, 需要求解不同的 PDE, 对应的数值方法也可 能完全不同 , 其间的切换和搭接非常困难 . 其三 , 对 于某些尺度上的流动 , 难以建立起有效的宏观输运 方程, 就更谈不上数值求解了. 从连续到稀薄这种跨 流域多尺度流动就是其中最为典型的例子 . 实际上 现存的大多数跨尺度计算方法采用的是直接分辨不 切实际的最小尺度, 即所谓的蛮力方法. 跨流域流动主要出现在微机电系统 (MEMS) 和 航天气动领域 . 近年来近空间飞行器研究的迅猛发 展更是对跨流域流动模拟提出了迫切的要求 . 无论 是飞船从太空(如 100 km 以上)返回地面还是飞行器
徐昆
①④
, 李启兵 *, 黎作武
②
③
① 香港科技大学数学系, 香港; ② 清华大学航天航空学院, 北京 100084; ③ 中国空气动力研究与发展中心, 绵阳 621000; ④ 北京大学湍流和复杂系统国家重点实验室, 北京 100871 *联系人, E-mail: lqb@tsinghua.edu.cn 收稿日期: 2013-09-04; 接受日期: 2013-11-20 国家自然科学基金 (批准号: 11172154)和香港大学教育资助委员会 (编号 : 621011)资助项目
中国科学: 物理学 力学 天文学
2014 年
第 44 卷
第5期
分子; 如果网格尺度在分子平均自由程附近, 看到的 就是分子之间的运动和碰撞 ; 如果网格尺度接近宏 观耗散尺度, 看到的就是大量分子的集体效应, 像流 体的传输和压力波的传播. 总之 , 计算流体力学的根本目的不是去解特定 的 PDE 方程, 而是对真实流体流动在不同网格尺度 上的直接刻画 . 或者说如果用一个连续变化的尺度 来观测流体运动, 自然界不可能只有现在孤立的 PDE 方程, 像 NS, Boltzmann 等, 而应该是一个连续 谱 . 本文的目的就是用直接计算建模的方法来捕捉 这样一个物理上连续变化的流体描述或 PDE 方程.
PACS: 47.11.-j, 02.60.-x, 51.10.+y, 47.45.-n doi: 10.1360/SSPMA2013-00054
1
引言
计算流体力学 (CFD) 通过数值方法求解流动问 题, 获得流场在离散空间中的定量描述, 进而预测流 体运动规律. 从 1922 年 Richardson 尝试采用数值计 算方式求解偏微分方程(PDE)开展天气预报工作开始, 随着高速计算机的出现, 到 20 世纪 60 年代, 传统计 算流体力学的框架体系逐步得到建立 , 进而成为与 理论解析和实验研究地位同等重要的流体力学研究 方法 [1,2] . 近半个世纪以来 , 计算机技术得到了迅猛 发展 , CFD 方法取得了长足进步 , 在航空航天等工
引用格式: 徐昆, 李启兵, 黎作武. 离散空间直接建模的计算流体力学方法. 中国科学: 物理学 力学 天文学, 2014, 44: 519–530
Xu K, Li Q B, Li Z W. Direct modeling-based computational fluid dynamics (in Chinese). Sci Sin-Phys Mech Astron, 2014, 44: 519–530, doi: 10.1360/SSPMA2013-00054
中国科学: 物理学 力学 天文学
2014 年
第 44 卷
第 5 期: 519–530 phys.scichina.com
《中国科学》杂志社
SCIENCE CHINA PRESS
SCIENTIA SINICA Physica, Mechanica & Astronomica
论文
离散空间直接建模的计算流体力学方法
业领域得到了广泛的应用. 然而, 现有数值方法与实 际应用的要求仍然存在很长的距离. 湍流、跨流域流 动等多尺度流动问题就是其中比较典型的例子 . 在 计算资源有限的条件下 , 工业领域迫切要求对流场 进行越来越细致的时空刻画 . 因此 , 发展新型高精 度、 高效的数值方法是 CFD 领域面临的重要挑战[3–5]. 传统 CFD 方法是通过直接离散描述流体运动的 偏微分方程而建立起来的, 如 Euler, Navier-Stokes(NS), Burnett, Boltzmann 方程等. 不同 PDE 的建立基础(假 设)不同, 因而具有不同的适用范围. 比如, NS 方程基 于连续介质假设, 只能适用于宏观流动. 而 Boltzmann
从地面运行到太空 , 都经历了包括稠密和稀薄大气 环境的变化, 给数值模拟带来了极大的困难. 根据飞 行器的大小和飞行高度估算出流动的稀薄程度 ( 用 Knudsen 数表示 Kn=/L, L 为流动的特征长度), 对飞 行器面临的流态进行划分. 当 Kn<0.001 时采用 NS 方 程求解; 而当 0.001<Kn<0.1 时也采用 NS 方程来近似, 同时还需要考虑飞行器表面的速度滑移和温度跳跃 . 对于更高 Knudsen 数的区域, 目前还未找到有效的宏 观方程 , 因而必须采用基于粒子的直接模拟 Monte Carlo 方法(DSMC)[7]. DSMC 方法完全基于单一尺度 ( 分子平均自由程 )来描述流体的运动 . 实际上 , 对于 同一飞行高度 , 飞行器的不同部位同样可能面临不 同的流态, 因而同样需要进行分区, 从而采用不同的 方法求解 . 这就需要在不同求解方法之间进行切换 , 并且界面上不同方法得到的数据如何交换也是一个 非常值得研究的问题 . 近年来航空航天领域关心的 近空间飞行器(20–100 km 之间)正是要在跨流域中长 时间飞行. 从理论上讲 , 适合从连续到稀薄跨流域流动的 PDE 是存在的, 那就是 Boltzmann 方程及其简化模型 方程. Boltzmann 方程本身的模型尺度是和 DSMC 的 尺度是一样的, 即分子的平均自由程. 通常说 Boltzmann 方程适合于所有尺度, 隐含的也是用蛮力 的意思, 在所有区域分辨出最小尺度. 这在绝大多数 实际应用中基本上是不可能的 . 再有 , 从 Boltzmann 方程出发通过不同的渐近展开可以得到不同的宏观 方程. 通过一阶 Chapman-Enskog 展开得到 NS 方程 是一个成功的例子 , 但在历史上还得到了很多不太 成功的方程, 比如更高阶的 Burnett 和 Super-Burnett 方程. 所以 Boltzmann 方程只是在其模型尺度上描述 流动, 本身并没有对其他尺度的适应性. 在数值上发 展一种跨尺度的格式比写出 Boltzmann 方程更困难. 对一维单原子气体流动, Boltzmann 方程可以写为
描述流体运动 , 对应最小尺度则是分子的平均自由 程 . 在时间方向也存在类似的尺度问题 , 这里不再 赘述. 需要注意的是, 当物理模型由微分方程来描述 后, 数学意义上的最小可分辨尺度变成了无限小(零 ). 如果采用离散模型 , 则最小可分辨尺度对应于离散 尺度. 换句话说, 物理模型表现出的性质会受到描述 方式的限制. 在 CFD 计算中, 数值离散本身将引入一个新的 最小可分辨尺度, 即网格尺度∆x. 小于网格尺度的流 动变化同样只能通过(数值 )模型来表达 . 最小尺度的 差异决定了数值解、 物理解和偏微分方程的解析ຫໍສະໝຸດ Baidu之 间存在差异. 为了准确捕捉最小尺度为 l 的流场脉动, 网格尺度必须满足∆x<l. 对直接离散 PDE 而建立的 传统 CFD 方法, 还必须要求 PDE 隐含的物理模型最 小尺度 lPDEl. 对于多尺度流动问题 , 一个优秀的 CFD 方法应该能按实际需要自动而高效地刻画出相 应尺度上的流动特征. 也就是说, 如果实际应用要求 给出小尺度上的流动, 在满足要求的网格尺度上, 数 值方法能捕捉到该小尺度上的物理量 ; 如果只要求 给出较大尺度上的流动特征 , 数值方法也能在相应 较粗的网格上给出合理的流场物理量. 当然, 一个比 较直观的想法是以流动的最小尺度来布置计算网格 , 从而能将大大小小的各种不同尺度的流场脉动同时 捕捉出来 . 不过 , 如高雷诺数湍流的直接数值模拟 (DNS)一样, 这种方法可能需要极其巨大的计算资源, 甚至根本不可行. 实际上, 科学研究的目的是对自然 的最有效的描述 , 如果所有真实的描述都要归结为 最细致的描述 , 那么世界上除了粒子物理之外就没 有其他科学了. 上面提到的采用不同尺度网格计算得到的流场 物理量之间的关系中, “合理”是指采用粗网格得到的 单元内物理量应该与采用细网格得到的流场按粗网 格尺度平均后得到的平均量近似相同. 同时, 为捕捉 边界条件等大尺度因素的影响 , 网格的尺度必须足 够小. 因此, 随着计算网格的逐步加密, CFD 计算得 到的流场平均量应该趋于定值 , 这就是网格收敛性 . 需要注意的是, 这里的网格收敛性与传统 CFD 计算 中的提法不同. 后者通常是指随着网格加密, 数值解 逐渐趋近于 PDE 的解析解; 而前者是指数值解与物 理解的逼近程度随网格尺度的变化 , 包含了不同尺 度下 PDE 本身对真实流动的逼近程度. 比如, 如果 网格大小和分子直径相近, CFD 应该可以看到独立的
摘要
本文提出了一种发展适合多尺度、 多物理流动的 CFD 格式的全新途径, 即直接在离散空间利用
物理模型的跨尺度演化解来建立数值格式的新方法. 与直接离散偏微分方程的传统做法相比, 基于离散 空间直接建模构造出的数值方法考虑了网格尺度和物理模型之间的匹配, 能做到不同尺度上物理模型的 连续过渡, 从而实现对多尺度流动的高效模拟. Boltzmann 方程的模型演化解为直接建模方法提供了重要 的支撑. 本文提出了直接建模构造计算流体力学方法的理论基础, 并给出了在建立统一气体动理学格式 中的成功运用, 验证了新方法的可行性和优越性. 关键词 偏微分方程数值解, 离散空间直接建模, 多尺度流动, 统一气体动理学格式, BGK 方程
徐昆等: 离散空间直接建模的计算流体力学方法
方程基于气体分子(粒子 )混沌假设, 其模型尺度是分 子平均碰撞时间和平均自由程. 虽然 Boltzmann 方程 理论上可以适用于从连续到稀薄的跨流域流动[6], 但 其要求在任何地方都能分辨出分子的微观运动尺度 , 也就是所谓的动理学尺度 (Kinetic Scale), 所以直接 把 Boltzmann 方程用于宏观流体的计算是不现实的. 这些不同的模型方程实际上是在不同尺度上描述流 体运动(下文将对此展开讨论), 因而逼近真实流动的 近似程度不同 . 在流体计算中由于网格尺度的加入 , 使得问题更加复杂 . 虽然通常认为当网格尺度足够 小时, 直接离散 PDE 得到的数值解能以任意的精度 逼近 PDE 解析解, 但它对真实物理流动的描述还受 限于 PDE 模型本身. 因此, 数值解逼近真实流动的 精度不光与计算网格有关, 而且与所选用的 PDE 模 型有关. 对于跨流域流动这种多尺度流动问题, 不同 流域需要选取不同的 PDE, 在过渡区, 甚至找不到一 个合适的 PDE 模型, 从而导致数值计算的难度太大. 一个很有意思的问题是 , 能否直接在离散空间上建 立流动的(数值 )物理模型, 从而通过数值解来直接逼 近物理解? 这也正是本文将要重点探讨的内容. 本文以多尺度流动问题为例 , 首先分析了不同 尺度下描述流动的物理模型 , 回顾了传统的基于直 接离散 PDE 的 CFD 方法, 进而提出发展 CFD 方法的 一条全新途径, 即在离散空间直接建模的方法. 以近 年来发展起来的适用于跨流域流动的统一气体动理 学格式 (UGKS) 为例 , 对新思路进行了验证和分析 . 另外, 当前传统的基于 Riemann 解的格式在高超声速 流动计算中遇到了瓶颈 , 现在是近一步思考和发展 下一代数值计算方法的关键时代 . 本文最后对此进 行了一定的讨论.
2
流动的描述尺度与物理模型
对真实流动的描述首先应该基于一定的最小可 分辨尺度 , 不同描述尺度下能看到不同的流动细节 , 以此可以建立不同的流动物理模型 , 通常为偏微分 方程. 事实上, 工程应用中关心的也是在某个有限分 辨率下的物理量, 如速度、温度和压力等, 或者说是 具有最小可分辨尺度的微元内的平均值 . 小于最小 尺度的流动变化则反映在物理模型中 , 如宏观方程 中耗散与内能的关系. 也就是说, 不同 PDE 模型隐 含的最小尺度 lPDE 不同. NS 方程基于连续介质假设, 对应宏观耗散尺度, 而 Boltzmann 方程从介观角度
3
基于直接离散 PDE 的 CFD 方法
传统 CFD 方法是通过直接离散描述流体运动的 偏微分方程建立起来的 . 借助偏微分方程数值解理 论 , 对简单流动 , 特别是具有线性变化特征的流动 , 可以比较容易地对 CFD 方法的基本性质, 如稳定性 和收敛性等进行理论分析. 实际计算中, 可以根据需 要通过加密网格 , 从而更好地逼近 PDE 的解析解 . 此时 , 数值解对物理解的逼近问题实际上是要解的 PDE 对真实流动的近似问题. 这里的理念(原理)实际 和计算关系不大 , 更像牛顿发明微积分用小单元取 极限描述连续曲线一样 . 但实际的计算一定是在有 限网格大小上进行的, 所以传统的 CFD 方法得到的 只是相应 PDE 的一个近似解, 而且在理论上通常也 不可能知道在数值上求解的确切方程. 对于多尺度流动问题 , 上述方法存在的不足就 更加明显 . 其一 , 对于非线性问题 , 基于直接离散 PDE 的数值格式的理论分析非常困难 , 至少到目前 为止 , 仍亟需分析方法的突破 . 其二 , 对于流动的不 同尺度, 需要求解不同的 PDE, 对应的数值方法也可 能完全不同 , 其间的切换和搭接非常困难 . 其三 , 对 于某些尺度上的流动 , 难以建立起有效的宏观输运 方程, 就更谈不上数值求解了. 从连续到稀薄这种跨 流域多尺度流动就是其中最为典型的例子 . 实际上 现存的大多数跨尺度计算方法采用的是直接分辨不 切实际的最小尺度, 即所谓的蛮力方法. 跨流域流动主要出现在微机电系统 (MEMS) 和 航天气动领域 . 近年来近空间飞行器研究的迅猛发 展更是对跨流域流动模拟提出了迫切的要求 . 无论 是飞船从太空(如 100 km 以上)返回地面还是飞行器
徐昆
①④
, 李启兵 *, 黎作武
②
③
① 香港科技大学数学系, 香港; ② 清华大学航天航空学院, 北京 100084; ③ 中国空气动力研究与发展中心, 绵阳 621000; ④ 北京大学湍流和复杂系统国家重点实验室, 北京 100871 *联系人, E-mail: lqb@tsinghua.edu.cn 收稿日期: 2013-09-04; 接受日期: 2013-11-20 国家自然科学基金 (批准号: 11172154)和香港大学教育资助委员会 (编号 : 621011)资助项目
中国科学: 物理学 力学 天文学
2014 年
第 44 卷
第5期
分子; 如果网格尺度在分子平均自由程附近, 看到的 就是分子之间的运动和碰撞 ; 如果网格尺度接近宏 观耗散尺度, 看到的就是大量分子的集体效应, 像流 体的传输和压力波的传播. 总之 , 计算流体力学的根本目的不是去解特定 的 PDE 方程, 而是对真实流体流动在不同网格尺度 上的直接刻画 . 或者说如果用一个连续变化的尺度 来观测流体运动, 自然界不可能只有现在孤立的 PDE 方程, 像 NS, Boltzmann 等, 而应该是一个连续 谱 . 本文的目的就是用直接计算建模的方法来捕捉 这样一个物理上连续变化的流体描述或 PDE 方程.
PACS: 47.11.-j, 02.60.-x, 51.10.+y, 47.45.-n doi: 10.1360/SSPMA2013-00054
1
引言
计算流体力学 (CFD) 通过数值方法求解流动问 题, 获得流场在离散空间中的定量描述, 进而预测流 体运动规律. 从 1922 年 Richardson 尝试采用数值计 算方式求解偏微分方程(PDE)开展天气预报工作开始, 随着高速计算机的出现, 到 20 世纪 60 年代, 传统计 算流体力学的框架体系逐步得到建立 , 进而成为与 理论解析和实验研究地位同等重要的流体力学研究 方法 [1,2] . 近半个世纪以来 , 计算机技术得到了迅猛 发展 , CFD 方法取得了长足进步 , 在航空航天等工
引用格式: 徐昆, 李启兵, 黎作武. 离散空间直接建模的计算流体力学方法. 中国科学: 物理学 力学 天文学, 2014, 44: 519–530
Xu K, Li Q B, Li Z W. Direct modeling-based computational fluid dynamics (in Chinese). Sci Sin-Phys Mech Astron, 2014, 44: 519–530, doi: 10.1360/SSPMA2013-00054
中国科学: 物理学 力学 天文学
2014 年
第 44 卷
第 5 期: 519–530 phys.scichina.com
《中国科学》杂志社
SCIENCE CHINA PRESS
SCIENTIA SINICA Physica, Mechanica & Astronomica
论文
离散空间直接建模的计算流体力学方法
业领域得到了广泛的应用. 然而, 现有数值方法与实 际应用的要求仍然存在很长的距离. 湍流、跨流域流 动等多尺度流动问题就是其中比较典型的例子 . 在 计算资源有限的条件下 , 工业领域迫切要求对流场 进行越来越细致的时空刻画 . 因此 , 发展新型高精 度、 高效的数值方法是 CFD 领域面临的重要挑战[3–5]. 传统 CFD 方法是通过直接离散描述流体运动的 偏微分方程而建立起来的, 如 Euler, Navier-Stokes(NS), Burnett, Boltzmann 方程等. 不同 PDE 的建立基础(假 设)不同, 因而具有不同的适用范围. 比如, NS 方程基 于连续介质假设, 只能适用于宏观流动. 而 Boltzmann
从地面运行到太空 , 都经历了包括稠密和稀薄大气 环境的变化, 给数值模拟带来了极大的困难. 根据飞 行器的大小和飞行高度估算出流动的稀薄程度 ( 用 Knudsen 数表示 Kn=/L, L 为流动的特征长度), 对飞 行器面临的流态进行划分. 当 Kn<0.001 时采用 NS 方 程求解; 而当 0.001<Kn<0.1 时也采用 NS 方程来近似, 同时还需要考虑飞行器表面的速度滑移和温度跳跃 . 对于更高 Knudsen 数的区域, 目前还未找到有效的宏 观方程 , 因而必须采用基于粒子的直接模拟 Monte Carlo 方法(DSMC)[7]. DSMC 方法完全基于单一尺度 ( 分子平均自由程 )来描述流体的运动 . 实际上 , 对于 同一飞行高度 , 飞行器的不同部位同样可能面临不 同的流态, 因而同样需要进行分区, 从而采用不同的 方法求解 . 这就需要在不同求解方法之间进行切换 , 并且界面上不同方法得到的数据如何交换也是一个 非常值得研究的问题 . 近年来航空航天领域关心的 近空间飞行器(20–100 km 之间)正是要在跨流域中长 时间飞行. 从理论上讲 , 适合从连续到稀薄跨流域流动的 PDE 是存在的, 那就是 Boltzmann 方程及其简化模型 方程. Boltzmann 方程本身的模型尺度是和 DSMC 的 尺度是一样的, 即分子的平均自由程. 通常说 Boltzmann 方程适合于所有尺度, 隐含的也是用蛮力 的意思, 在所有区域分辨出最小尺度. 这在绝大多数 实际应用中基本上是不可能的 . 再有 , 从 Boltzmann 方程出发通过不同的渐近展开可以得到不同的宏观 方程. 通过一阶 Chapman-Enskog 展开得到 NS 方程 是一个成功的例子 , 但在历史上还得到了很多不太 成功的方程, 比如更高阶的 Burnett 和 Super-Burnett 方程. 所以 Boltzmann 方程只是在其模型尺度上描述 流动, 本身并没有对其他尺度的适应性. 在数值上发 展一种跨尺度的格式比写出 Boltzmann 方程更困难. 对一维单原子气体流动, Boltzmann 方程可以写为