不可压缩流场多体运动问题的两种数值解法.pdfVIP

第23卷第5期 　 计 算 力 学 学 报 　 Vol . 2 3 ,No . 5 2006 年 10 月 Chinese Journal of Computational Mechanics October 2006 文章编号 (2006) 不可压缩流场多体运动问题的两种数值解法 1 2 于 　军 , 　聂义勇 ( 1. 中国科学院 机器人学重点实验室 ,沈阳 1100 16 ;2 . 中国科学院 研究生院 ,北京 100039) 摘 　要 :不可压缩流场中多体运动的数值模拟是计算流体动力学的前沿课题 , 目前可以求解此问题的两种高精度 方法是重叠网格法和移动网格法 。本文详细讨论了两种方法的理论基础 ,典型算法的实现过程 ,并采用两种方法 对单圆柱和多圆柱运动问题进行了数值试验 。通过分析比较 ,本文认为重叠网格法计算效率高于移动网格法 ,重 叠网格法比移动网格法更适合于求解物体大位移运动问题 ,尤其当计算区域几何形状复杂时重叠网格法的优势 更加明显 。 关键词 :不可压缩流场 ;多体运动 ;重叠网格方法 ;移动网格方法 中图分类号 :O242 . 2 1 　　　文献标识码 :A 1 　引 言 解决物体在流体中运动问题时 ,在一个时间步内分 ( 为以下步骤 :根据流体内边界节点的位置 由物体 应用差分法或者有限元法求解微分方程初边 ) 在上一时刻的位移确定 采用网格移动方法计算流 值问题时一般对其求解区域剖分形成一套网格 ,在 体内网格节点位置 ,重新生成网格 ;采用物理守恒 网格上离散微分方程 ,计算出结果并进行分析处 量重映方法将旧网格上的物理量映射到新网格上 , 理 。当流体中存在运动物体时 ,通常采用坐标系变 并将其作为该时间步上物理量初值 ;在新网格上对 换消去物体的运动变量 ,使物体好像是固定在区域 流体运动控制方程组离散 , 以物体速度作为流体运 内,类似风洞试验 。这时仅有流体控制方程组 ,求 动内边界条件 ,结合外边界条件 ,求解出流场速度 解区域不会随时间变化 ,单求流体运动就可以解决 和压力 ;对物体表面流体的张力和压力积分 ,计算 该问题[ 1 ,2 ] 。对于本文讨论的多体运动问题 ,很难 其受到的流体作用力 ,叠加其他作用力得到其所受 进行坐标系变换以消去物体的运动变量 ,因为物体 合力 ;根据刚体运动方程计算物体位移 ,得到物体 之间存在相对运动 ,并且物体运动不能预先确定 , 下一时刻位置和速度 ,移动物体 。依次重复这些步 ( ) 所以物体 主动者 运动和流体运动必须耦合求解 。