混合网格重叠方法在多体相对运动中的应用.pdfVIP

第 7卷 第 l期 Vo1_7 No．1 2015年 2月 Feb．2O15 混合网格重叠方法在多体相对运动中的应用 康忠 良 方媛媛 (I．中国建筑科学研 究院建研科技股份有 限公司，北京 100013； 2．北京市燃气集 团研究院，北京 100011) 【摘 要】本文发展 了混合网格的动态重叠方法。基于线性重构方法，给出了一种适用于混合 网格间的插值策略。 所建立的网格间边界定义方法的性能优秀可靠，网格装配后的插值边界光滑且 网格匹配性较好。耦合求解非定常 Ns方程和刚体 6DOF运动方程来模拟 多体相对运动 问题。外挂物分离问题的数值计算表明，所发展的动网格方法 对于处理 多体相对运动模拟 问题是快速准确可靠的。 【关键词】动网格；重叠网格；混合网格；多体运动 【中图分类号】V211．3 【文献标识码】A 【文章编号】1674—7461(2015)O1—0069—06 散。另一方面，在实际工程应用 中还应该强调动 网 1 引言 格处理的 自动化程度、最大限度地减少人工干预 ， 否则也不能体现混合网格重叠方法的优势。 本文在优化网格问边界定义的基础上 ，给出一 在计算流体力学 (CFD)应用 中，采用混合 网格 种适用于任意单元类型的插值策略，并给出一种完 方法不但对复杂构型具有强大 的几何适应能力 ，网 全 自动化的高效的网格动态装配策略。 格生成 的人工工作量少 ，而且 易于通过采用不 同的 单元类型及调整 网格 的疏密来适应不同的流场特 2 数值方法 征 ，并容易生成整体 网格 、整体求解 。但是 ，在实际 应用中，混合 网格方法也遇 到了一些 困难 ，比如对 本文算法构造基于任意多面体混合 网格单元。 于多体相对运动 问题 ，就需要对 网格进行再生 ，处 算法不考虑具体的网格拓扑，对不同的单元类型统一 理非常复杂。 处理。空间离散采用格心有限体积法 ，控制体取网格 为了解决这一 问题 ，Nakahashi¨首先提 出了非 单元 ，控制面取单元表面，流场变量存储在单元中心。 结构网格重叠方法。该方法吸收了结构 网格重叠方 2．1 控制方程 法 的优点 ，后来又被推广到混合网格重叠及动态重 舍去源项的三维非定常可压缩 Ns方程组在直 叠应用中 ，得到了较大的发展。混合 网格动态重 角坐标系下的守恒积分形式可表示为 叠方法的优点就是在模拟多体间具有大幅相对运动 和 + F)dS=0 (1) 问题时，在物体运动过程 中不需要网格再生，算法实 式中： 为控制体 ；d 为控制面 ；Q为守恒变矢量； 现简单 ，并且可以获得较高的动网格处理效率。 F 为对流通量；F 为粘性通量。各项的具体描述详 尽管如此 ，混合网格重叠方法 目前仍然面临诸 多难题，其 中，网格间边界的定义及其插值方法就 见文献 。 是其中一个重要方面 ，已有研究均基于多层 网格节 2．2 数值离散方法 点的流场变