一种全耦合多相流分析的并行计算方法力学学报.pdfVIP

第 31 卷 第 3 期 力 　学 　学 　报 Vol. 31 , No. 3 1999 年 5 月 ACTA MECHANICA SINICA May ,1999 一种全耦合多相流分析的并行计算方法1) 王希诚 (大连理工大学力学系工业装备结构分析国家重点实验室 ,大连　116023) 摘要 　研究了孔隙介质中热、水和汽流全耦合分析的并行计算方法. 模型中采用了考虑毛细压 力关系的修正有效应力概念 ,并考虑了相变和潜热传递. 基本变量为位移、毛细压力、汽压和温 ( ) ( 度. 并行程序是在国家高性能计算中心 北京 的曙光 1000A 上借助 PVM Parallel Virtual Machine) 软件系统实现的 ,考题显示出较高的并行加速比和效率. 关键词 　并行计算 ,多相流 ,可变形孔隙介质 ,相变 引　言 多相流问题的研究对于石油和天然气的开采 ,地热的开发 ,环境保护 ,海底隧道及大型地 下建筑等工程领域具有十分重大的意义 ,成为近年来一个热门研究领域[1 ] . 由于问题的非线性 以及求解的时间跨度较长,多相流分析具有高度计算复杂性 ,高速并行计算是解决这一困难的 关键. 然而 ,在多相流分析领域中有效并行算法的研究尚不多见. 本文在分析该问题的内在并 行性的基础上 ,提出一种适于非线性分析的并行波前算法 ,用于求解孔隙介质中包括热和质量 传递的全耦合多相流问题. 本算法是将问题的求解区域划分成一些子区域 ,在子区域内利用波 前法[2 ,3 ] 并行进行组装和内部自由度的消元过程. 当每一个波前前进到各自的边界时 ,从波前 活动变量核心数组提取相应系数组装界面方程. 求解界面方程得出界面节点自由度值 ,再返回 各子区域利用波前回代过程求出域内节点的自由度值. 与一般基于带宽格式的区域分裂法相 比,本方法不限制有限元节点的编号顺序 , 内存需要量少,编程容易. 本算法已在国家高性能计算 ( ) 中心 北京 的曙光 1000A 上借助 PVM 系统予以实现 ,显示出较高的并行加速比和计算效率. 1 问题的数学模型 1 . 1 本构关系 将饱和度表示为毛细压力和温度的函数 ( ) ( ) S w = S w p c , T 1 毛细压力 p c 为 ( ) p c = p g - p w 2 ( ) p g , p w 分别为潮湿空气和流体水压力 ,孔隙系统中的潮湿空气假定为干空气 ga 与水蒸汽 (gw) 两种理想气体的完全混合 ,服从理想气体状态方程. 　1997 - 04 - 29 收到第一稿 , 1998 - 09 - 21 收到修改稿. 1) 国家自然科学基金和国家高性能计算基金资助项目. 第 3 期 王希诚 :一种全耦合多相流分析的并行计算方法 　　　 277 ρ ρ ( ) p g