2气泡群溃灭过程的数值模拟.doc

第八届全国流体力学学术会议 2014年9月18～21日 甘肃兰州 文章编号： CSTAM2014-B01-0130 标题：气泡群溃灭过程中的主控参数及影响规律研究 作者：杜特专，黄晨光，王一伟，于娴娴，吴小翠 单位：中国科学院力学研究所流固耦合系统力学重点实验室 Copyright ? 2014 版权所有 中国力学学会地址: 北京市北四环西路15号　邮政编码:100190 　Address: No.15 Beisihuanxi Road, Beijing 10 杜特专*，黄晨光*，王一伟*，于娴娴*，吴小翠* *（中国科学院力学研究所流固耦合系统力学重点实验室，北京海淀区 100190）获空化[1]是水下高速航行体和高速推进器面临的主要问题之一，其可能诱发系统的异常运行状态[2]，结构的变形[3]、破坏和噪声[4]等。本质上是由大量微小空泡组成空泡群Reisman等人[5][6]通过实验测量了水翼表面云状空化溃灭的压力及其传播，其典型溃灭压力幅值和作用时间分别[7]首先采用混合物连续性方程和动量方程与Rayleigh-Plesset方程[8][9]耦合的方法计算了一层气泡群在固壁附近的溃灭过程。随后Chahine[10]，Noordzij和van wijngaarden[11]，Kameda和Matsumoto[12]在此基础上分析了压力波的传播，d’Agostino和Brennen[13]分析了球状气泡群的自然频率。目前大部分的工作主要集中在线性分析和弱非线性分析方面，对于强非线性的研究目前还很少。Wang和Brennen[14][15]采用混合物模型结合Rayleigh-Plesset方程研究了气泡大小均匀的球形空泡群的非线性生长和溃灭，引入?=?0A02/R02来表征气泡间互相影响的大小。?越小，气泡间的相互影响越小，气泡群的自然频率越接近于单气泡。当?大于1时，气泡间的相互影响就不能忽略，此时气泡群的自然频率低于单个气泡自然频率。Wang[16]对包含不同尺寸气泡的气泡群进行了研究，发现气泡大小分布会对溃灭压力波传播规律和气泡群宏观动力学特性产生影响。Kumar和Brennen[17]通过研究发现小气泡组成的气泡群的自然频率远高于大气泡组成的气泡群。国内对气泡群的研究起步于80年代末，黄建波，倪汉根[18]基于分层溃灭的假定分析了固壁附近空泡群的溃灭过程。张阿漫，姚熊亮等[19]采用边界元方法研究了三维气泡群中气泡之间的相互作用。张凌新[20]，邵雪明[21]等基于气泡群的直接模拟研究了气泡间的相互作用对空泡群溃灭的影响。 本文以规则排列的气泡群为研究对象，通过量纲分析获得了气泡群溃灭过程中的主控无量纲参数，并通过数值模拟分析了主控参数对气泡群溃灭的影响规律，进而得到了气泡群溃灭速度的经验公式。 1 气泡群溃灭过程的量纲分析 考虑特征尺度为a的气泡群在远场压力作用下的溃灭过程，假设初始气泡为球形，泡内压力为，如图1所示： 气泡群示意图 针对这样的气泡群溃灭过程，控制参数包括（中括号内为物理量的量纲）： 时间：； 几何参数：初始泡径，控制体特征长度； 控制体内气相的气泡数密度； 物性参数：液体密度，液体粘性，表面张力系数，气体参考密度，气体粘性，气体比热比； 外界压力为，泡内压力为。 则气泡群的体积变化率为： ……(1) 取为基本量，进行无量纲化，同时考虑到气泡溃灭过程中驱动力是以压差的形式出现，可将改写成，则有： ……(2) 其中可将改写为，表示气泡在控制体内所占的体积分数；可改写为，表示控制体内的气泡数；表征单气泡溃灭特征时间（），用替代，整理可得： ……(3) 其中与为常数，对于气泡溃灭过程，主要是压差和惯性力起主要作用，可忽略粘性和表面张力的影响，因此（3）式可进一步简化为： ……(4) 由于可能为零，为了方便可用替代；若体积变化率为溃灭过程的平均值，则不显含t，即： ……(5) （5）式中代表单个球泡的平均体积变化率，所以表征气泡群的无量纲平均体积变化率。为气泡群中气相的体积分数，N表示气泡群中的气泡数。相同的条件下，越大则气泡间距越小，气泡间相互作用越强；越小，气泡间距越大，气泡间的相互作用越弱。为无量纲压力，该数越接近于零表示驱动力越小，气泡群处于缓慢压缩的过程；相反的，无量纲压力越接近于一说明驱动力越大，气泡群溃灭越剧烈。 为了确定（5）式的具体形式，可设…..(6) 其中为经验参数，下面我们将通过对气泡群溃灭过程的数值模拟来定出这些参数。 2 气泡群溃灭过程的数值模拟 2.1 计算模型 如图2所示，