2012
年
7
月
第
33
卷第
4
期
郑州大学学报(工学版)
Jul. 2012
Vol.33
No.4
Journal of Zhengzhou University (Engineering Science)
文章编号:
1671
-6833(2012)04
-0060
-05
空化水射流空泡溃灭过程的数值分析
杨博凯,卢义玉,杨晓峰,冯明涛,张
赛
(重庆大学资源及环境科学学院,重庆
400030)
摘
要:空化空泡发育与溃灭的计算是研究空化水射流的一个重点和难点对比分析了空泡运动方程
即
Ray
leigh-Plesset
方程的不同数值解法,并针对固定时间步长的缺点,提出了变步长法解
Rayleigh
Plesset
方程的思路.通过优化系数
λ
=
R,IR,_1
,得出如下分析结果:在空泡内外压力变化的不同情况下,
变步长法均优于固定步长法且没有奇异点.变步长法在节约计算量,保证结采精度等方面都有较好表
现,该方法为研究空化水射流提供新的思路和算法.
关键词:空泡;
Rayleigh-Plesset
方程;固定时间步长;交步长法
中图分类号:
TB126
文献标志码:
A
doi:lO.
3969/j.
issn.1671-6833.2012.
04. 014
0
51
言
目前,高压空化水射流已经应用于很多领域,
如岩石破碎、海上石油钻井、船舶清洗、空化水射
流处理有机物废水以及医疗技术等,但目前对于
空泡生成和愤灭的机理还不清楚.
对于空化水射流而言,其空泡半径与压力、时
间存在确定的函数关系,而三者的关系直接影响
空化射流的效率以及能量释放程度.因此,旨在研
究它们之间关系的空泡动力学方程在空化理论中
占有重要的地位.
对于一个包含空泡的液体,如果已知射流中
单个空泡的内外压力,就可以计算它的半径
R
[I
J '
同时能够建立一个
Rayleigh-Plesset
方程来近似求
得气泡的运动状态[
2 J
.理论上这个方程可以使用
常见的任意一种有限元差分法解决,但由于空泡
半径随时间变化的过程是非线性的,常用的迭代
方法在求解
Rayleigh-Plesset
方程时存在困难,计
算量大、精确度低,特别是在空泡渍灭阶段会产生
振荡.笔者根据变步长法的基本思想,结合自主研
发的空化水射流对各种迭代法的计算效率进行分
析,提出了变步长的
Rayleigh-Plesset
方程计算方
法.变步长法能够更有效地模拟空泡最终愤灭阶
收稿日期:
2012
-01
-08
;修订日期:
2012-03-10
段的能量释放过程,而这个能量释放过程正是研
究其对附近边界产生破坏作用机理的关键[
3
叫.
1
Rayleigh-Plesset
空泡运动方程
Rayleigh-Plesset
方程是解决空泡运动方程的
经典方程,在空化水射流的研究中,经常需要确定
空泡随时间的变化以及最终愤灭放能情况,而通
过
Rayleigh-Plesset
方程可以较为准确地得到空泡
半径与时间的关系.它是一个二阶微分方程[
5 J •
『
d2
R 3 dR
21
吁,,
4udR
Pl
R
ci7
+言(忑
)
j=pn-p
.,-
立一
~dt
(1)
式中:
p
为液体密度;
μ
为液体运动教度;
t
为时
间;
p
.,
为液体静压;
Pn
为空泡内的压力.
式(
1
)可以转换为以下的一阶方程
去
= y
,川=
Xo (
2)
去=
f(
t,
x,
y)
, y ( t0)
=儿(
3)
当
z
为因变量时,令
R
=x
,可同时求解式(
2)
和(
3
)以得到方程的解.令:
h = Llt
=ti
- ti I
(4)
Llx,
x i
Llt
i
(5)
基金项目:国家自然科学基金委专项创新研究群体基金资助项目(
50621403
);重庆市杰出青年基金资助项目
( 2009 BA604
7
);重庆大学研究生科技创新基金资助项目(
200802Bl
80010289).
作者简介:杨博凯(
1985 -
),男,天津人,重庆大学研究生,主要从事高压水射流方面的研究.
E-mail: 20046013@ cqu.
edu. cn.
2012
年
7
月
第
33
卷第
4
期
郑州大学学报(工学版)
Jul. 2012
Vol.33
No.4
Journal of Zhengzhou University (Engineering Science)
文章编号:
1671
-6833(2012)04
-0060
-05
空化水射流空泡溃灭过程的数值分析
杨博凯,卢义玉,杨晓峰,冯明涛,张
赛
(重庆大学资源及环境科学学院,重庆
400030)
摘
要:空化空泡发育与溃灭的计算是研究空化水射流的一个重点和难点对比分析了空泡运动方程
即
Ray
leigh-Plesset
方程的不同数值解法,并针对固定时间步长的缺点,提出了变步长法解
Rayleigh
Plesset
方程的思路.通过优化系数
λ
=
R,IR,_1
,得出如下分析结果:在空泡内外压力变化的不同情况下,
变步长法均优于固定步长法且没有奇异点.变步长法在节约计算量,保证结采精度等方面都有较好表
现,该方法为研究空化水射流提供新的思路和算法.
关键词:空泡;
Rayleigh-Plesset
方程;固定时间步长;交步长法
中图分类号:
TB126
文献标志码:
A
doi:lO.
3969/j.
issn.1671-6833.2012.
04. 014
0
51
言
目前,高压空化水射流已经应用于很多领域,
如岩石破碎、海上石油钻井、船舶清洗、空化水射
流处理有机物废水以及医疗技术等,但目前对于
空泡生成和愤灭的机理还不清楚.
对于空化水射流而言,其空泡半径与压力、时
间存在确定的函数关系,而三者的关系直接影响
空化射流的效率以及能量释放程度.因此,旨在研
究它们之间关系的空泡动力学方程在空化理论中
占有重要的地位.
对于一个包含空泡的液体,如果已知射流中
单个空泡的内外压力,就可以计算它的半径
R
[I
J '
同时能够建立一个
Rayleigh-Plesset
方程来近似求
得气泡的运动状态[
2 J
.理论上这个方程可以使用
常见的任意一种有限元差分法解决,但由于空泡
半径随时间变化的过程是非线性的,常用的迭代
方法在求解
Rayleigh-Plesset
方程时存在困难,计
算量大、精确度低,特别是在空泡渍灭阶段会产生
振荡.笔者根据变步长法的基本思想,结合自主研
发的空化水射流对各种迭代法的计算效率进行分
析,提出了变步长的
Rayleigh-Plesset
方程计算方
法.变步长法能够更有效地模拟空泡最终愤灭阶
收稿日期:
2012
-01
-08
;修订日期:
2012-03-10
段的能量释放过程,而这个能量释放过程正是研
究其对附近边界产生破坏作用机理的关键[
3
叫.
1
Rayleigh-Plesset
空泡运动方程
Rayleigh-Plesset
方程是解决空泡运动方程的
经典方程,在空化水射流的研究中,经常需要确定
空泡随时间的变化以及最终愤灭放能情况,而通
过
Rayleigh-Plesset
方程可以较为准确地得到空泡
半径与时间的关系.它是一个二阶微分方程[
5 J •
『
d2
R 3 dR
21
吁,,
4udR
Pl
R
ci7
+言(忑
)
j=pn-p
.,-
立一
~dt
(1)
式中:
p
为液体密度;
μ
为液体运动教度;
t
为时
间;
p
.,
为液体静压;
Pn
为空泡内的压力.
式(
1
)可以转换为以下的一阶方程
去
= y
,川=
Xo (
2)
去=
f(
t,
x,
y)
, y ( t0)
=儿(
3)
当
z
为因变量时,令
R
=x
,可同时求解式(
2)
和(
3
)以得到方程的解.令:
h = Llt
=ti
- ti I
(4)
Llx,
x i
Llt
i
(5)
基金项目:国家自然科学基金委专项创新研究群体基金资助项目(
50621403
);重庆市杰出青年基金资助项目
( 2009 BA604
7
);重庆大学研究生科技创新基金资助项目(
200802Bl
80010289).
作者简介:杨博凯(
1985 -
),男,天津人,重庆大学研究生,主要从事高压水射流方面的研究.
E-mail: 20046013@ cqu.
edu. cn.