瞬态空化泡演变过程的数值模拟Ξ.pdfVIP

第 14 卷　第 2 期 应 用 力 学 学 报 . 14　 . 2 V o l N o 1997 年 6 月 CH INESE JOURNAL OF APPL IED M ECHAN ICS Jun. 1997 瞬态空化泡演变过程的数值模拟 蔡悦斌　鲁传敬　何友声 (上海交通大学　上海　200030) 摘　　要 采用边界积分方程方法, 对无粘流体中三个空化泡以及自由面附近二个空化泡 相互作用的演变过程进行了数值模拟。计算中边界用二阶有限元离散, 影响系数矩阵 非对角线元素用六点高斯数值积分方法计算, 对第一类、第二类完全椭圆积分用高次 多项式近似, 对计算系数矩阵对角线元素中遇到的奇异积分进行了特殊处理。结果表 明, 在不同的给定参数下, 空化泡的溃灭形态各异, 柱状射流和环形射流都有可能发 生, 使空化泡演变成双泡或环形泡。 关键词: 边界积分方程方法; 空化泡; 射流 1　引　言 近二十年来, 国外学者对瞬态空化泡在自由液面与固壁面附近的溃灭现象进行了细致深 入的理论和实验研究。气泡在自由液面附近的成长与溃灭研究最早是由水下爆炸兴波问题引 出的, 而空化泡在固壁附近的溃灭则与空化空蚀现象有关。现认为空蚀是由靠近固体表面的空 泡崩溃时产生的高速液体射流冲击壁面造成的。已发现空化泡在柔性壁面与固壁面附近的溃 灭形态存在着很大差别[ 1 ]～ [ 7 ]。 瞬态空化泡的动力学研究, 在空蚀和空化噪声等工程问题上有重要意义, 越来越引起人们 的注意。在研究自由液面与固壁面附近单个空化泡成长溃灭过程和二个空化泡之间的相互作 [ 7 ] 用的基础上 , 本文进一步研究了自由面附近二个瞬态空化泡及多个泡之间的相互作用过程。 2　数学模型 流体中瞬态泡的演变是一个非常复杂的物理过程。但除了泡的最初成长阶段和发生崩溃 的最后阶段外, 可以忽略表面张力、相变和热力学效应等因素的影响, 假定流体是理想、不可压 的、流动无旋有势。现考虑轴对称问题, 取坐标系如图 1 所示。引入速度 , 取空泡在无界流场 中膨胀的最大半径R m 、流体密度 和 P = P ∞ - P c 为特征量, 这里P ∞ 为环境压力, P c 为饱 本项目研究得到国家自然科学基金资助 来稿日期: © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 2 应 用 力 学 学 报 第 14 卷 和蒸汽压力。无量纲化后, 满足的方程与边界条 件如下: 2 1 ( r ) + 2 = 0　　 流域 内 ( 1) r r r z D r ( ) = S b 与S f 上 2 D t D 1