自激振荡脉冲喷嘴空化效应及其射流形态的数值分析汪朝晖１㊀胡亚男１㊀饶长健１㊀邓晓刚２１．武汉科技大学机械自动化学院,武汉,４３００８１２．重庆科技学院机械与动力工程学院,重庆,４０１３３１摘要:基于自激振荡脉冲喷嘴空化效应和多相流模型,建立了自激振荡脉冲射流空化模型.依据自激振荡腔室结构及其几何参数建立了腔室轴对称物理模型,计算得到了振荡周期１００m s 内自激振荡脉冲射流的空化泡破碎㊁腔室内两相分布㊁湍动能分布和速度分布等结果.研究表明:在１．０２~２．３７m s时,空化泡半径减小,气泡开始径向运动形成泡面加速射流;在２．６９~４．６７m s 时,空化泡面压力达到极限破碎值时气泡开始破碎;在自激振荡周期前２５m s ,主射流与空气接触边界面形成较强湍动能,自激振荡腔室中心漩涡区逐渐变大,外流场连续射流被割断成多股状射流,射流在喷射轴线附近速度达到并稳定在３０~４０m /s ;在振荡周期的４０~９０m s ,腔室内中心空化气囊形成并开始阻挡主射流运动,喷嘴出口流道出现大面积空化区域,湍动能最大区域集中在下喷嘴出口下游;在振荡后期,随着主射流与空气相互作用及射流贯穿距离增加,主射流速度逐渐趋于稳定且扩散作用减弱.关键词:自激振荡脉冲喷嘴;空化效应;射流形态;空化泡;雾化中图分类号:T P ６９D O I :１０．３９６９/j．i s s n ．１００４ １３２X．２０１７．１３．００４N u m e r i c a l A n a l ys i s o fC a v i t a t i o nE f f e c t s o f S e l f Ｇe x c i t e dO s c i l l a t i o nP u l s eN o z z l e s a n d J e t F o r m s WA N GZ h a o h u i １㊀HU Y a n a n １㊀R A O C h a n g j i a n １㊀D E N G X i a o g a n g２１．S c h o o l o fM a c h i n e r y a n dA u t o m a t i o n ,W u h a nU n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y,W u h a n ,４３００８１２．C o l l e g e o fM e c h a n i c a l a n dP o w e rE n g i n e e r i n g ,C h o n g q i n g U n i v e r s i t y of S c i e n c e a n dT e c h n o l og y ,Ch o n g qi n g,４０１３３１A b s t r a c t :B a s e do nt h ec a v i t a t i o ne f f e c t so f t h es e l f Ｇe x c i t e do s c i l l a t i o n p u l s e dn o z z l e sa n d m u l t i Ｇph a s e f l o w m o d e l ,ac a v i t a t i o n m o d e l o f t h es e l f Ｇe x c i t e do s c i l l a t i o n p u l s e d j e t sw a se s t a b l i s h e d ．T h e a x i s y mm e t r i c p h y s i c a lm o d e l o f t h e c h a m b e rw a s b u i l t a c c o r d i n gt o t h e s e l f Ｇe x c i t e d o s c i l l a t i o n c h a m b e r s t r u c t u r e s a n d i t s g e o m e t r i c p a r a m e t e r s ．T h e b r e a k u p o f c a v i t a t i o n b u b b l e s ,t w o Ｇph a s e d i s t r i b u t i o n s i n t h e c h a m b e r s ,t u r b u l e n t k i n e t i c e n e r g y d i s t r i b u t i o n s a n dv e l o c i t y di s t r i b u t i o n s i n a no s c i l l a t i o n p e r i o d o f １００m sw e r e o b t a i n e d b y nu m e r i c a l c a l c u l a t i o n s ．T h e r e s u l t s s h o wt h a t :i n １．０２Ｇ２．３７m s ,t h e r a d i i o f c a v i t a t i o nb u b b l e s d e c r e a s e r a p i d l y a n d t h e b u b b l e s b e g i n t o d o r a d i a lm o t i o n s t o f o r mt h e a c c e l e r a t i n gje t o n t h e b u b b l e s s u rf a c e ．I n ２．６９Ｇ４．６７m s ,t h e b u b b l e s s t a r t t o c r u s hw h e n t h e p r e s s u r e s o f c a v i t a t i o n b u b b l e s u r f a c e s r e a c h t h e l i m i t b r e a k i ng v a l u e ．I n th e b e gi n n i n g o f ２５m s ,t h e s t r o n g t u r b u l e n t k i n e t i c e n e r g yi s f o r m e d a t t h e c o n t a c t i n t e r f a c eb e t w e e nt h em a i n j e t a n da i r ,a n d t h e c e n t e rv o r t e xa r e a i n s e l f Ｇe x c i t e do s c i l l a t i o n c h a m b e r g r a d u a l l yg r o w s u p ．C o n t i n u o u s je t i n t h e o u tf l o wf i e l d i s c u t i n t o t h e m u l t i p l e j e t a n d t h e j e t v e l o c i t y r e a c h e s a s t a b l e v a l u e a b o u t ３０Ｇ４０m /s n e a r t h e a x i s o f i n je c t i o n ．I n ４０Ｇ９０m s ,t h e c e n t r a l c a v i t a t i o na i r b a g i nt h ec h a m b e rh a sb e e nf o r m e da n db e gi n s t ob l o c kt h e m o v e Ｇm e n t s o fm a i n j e t ,a n d t h e l a r g e c a v i t a t i o na r e aa p pe a r s i nt h en o z z l e so u t l e t c h a n n e l a n d t h em a x i Ｇm u mt u r b u l e n t k i n e t i c e n e r g yi s c o n c e n t r a t e d i n t h e d o w n s t r e a mo f t h e l o w e r n o z z l e e x i t s ．I n t h e l a t e r p e r i o do f o s c i l l a t i o n ,t h e v e l o c i t y o fm a i n je t t e n d s t ob e s t a b l e a n d t h e d if f u s i o n e f f e c t s a r ew e a k e n e d a s t h em a i n j e t i n t e r a c t sw i t h t h e a i r a n d t h e p e n e t r a t i o nd i s t a n c e s o f t h e je t i n c r e a s e ．K e y wo r d s :s e l f Ｇe x c i t e do s c i l l a t i o n p u l s e n o z z l e ;c a v i t a t i o n e f f e c t ;j e t f o r m ;c a v i t a t i o nb u b b l e ;a t Ｇo m i z a t i o n收稿日期:２０１６Ｇ０８Ｇ２５基金项目:国家自然科学基金资助项目(５１４０５３５２,５１３７６２０４)０㊀引言空化是液流系统中由于局部压力低于临界值而诱发液体内部空泡的产生㊁发展和溃灭的过程,高速射流在进入狭窄的喷嘴内部时往往伴随着复杂的湍流运动和介质密度的变化,极易形成空化[１Ｇ２].研究表明:当空化泡破碎时,泡面微射流５３５１ 自激振荡脉冲喷嘴空化效应及其射流形态的数值分析汪朝晖㊀胡亚男㊀饶长健等Copyright©博看网 . All Rights Reserved.瞬间可达到高速㊁高压和高温的状态.当喷嘴内部出现空化现象时,有效利用其空化特性可使流体雾化效果显著增强;但由于实验中无法清楚地观察空化射流形成及其喷雾形态,尤其是对空化射流的喷射流体结构和外部流场分布研究仍然存在盲区,因此利用数值计算手段研究空化效应及其射流形态具有重要意义.空化效应的形成机理极其复杂,空化泡的产生㊁溃灭以及射流不稳定性是影响空化射流形态的重要因素.S UH等[３]将喷油嘴简化成扩大的透明二维矩形空壳结构,研究喷孔在不同几何结构下空穴流动对喷雾的影响,以及空穴的发展过程,探究了空化泡对空化射流雾化的影响.A L EＧHO S S E I N等[４]和Q I N等[５]通过求解R a y l e i g hＧP l e s s e t方程模拟了空化射流中空化泡的生成㊁溃灭的变化过程,得出空化泡影响空化射流形成的结论.姚立明等[６]为研究空化喷嘴产生的空泡射流在不同深度环境下的影响因素,研究了不同环境下空泡的差异,并总结了空泡内气相体积分数和速度的变化规律.C H E N等[７]利用数值模拟手段,证实了空化射流对钢铁腐蚀的主要原因在于空化泡溃灭时产生的高强度冲击波和瞬间的高温.刘琦等[８]分析了喷孔内部气液两相流场的三维流态以及空化流动特性,发现随着喷孔直径增大,燃油高速区域扩大,空化效应增强.目前空化射流形态研究主要关注空化模型及其算法的改进.X I E等[９]建立了一种简化虚拟流体模型模拟空化射流中气泡的爆炸现象,而用另外等熵单流体空化模型描述和捕捉非定常的空化射流形态,在模拟空化射流冲击自由壁面中得到较好的计算结果.WA N G等[１０]结合拉格朗日计算方法和欧拉计算方法的优点,提出了混合计算模型,对变化剧烈的空化区域使用欧拉模型进行计算,而对于相对稳定的刚性壁面等位置通过拉格朗日网格进行计算,提高了空化模型模拟精度.曾宇杰等[１１]利用两相流的数值计算模拟血液的两相流动,得到血管壁面压力㊁壁面剪切应力㊁血液流速等血流动力学参数,弥补了单相流模型的不足.王维军等[１２]采用改进后的空化模型和湍流模型,对离心泵内部空化流动进行分析,提出了空化初生的判定准则,系统划定了离心泵的空化区域.上述计算表明:有效的空化模型及其改进算法有助于空化射流的数值模拟研究.根据H e l m h o l t z空腔模型而设计的自激振荡腔室能够产生较强的脉冲射流,其脉冲作用由封闭的振荡腔内上喷嘴出口处形成的自由剪切射流与下喷嘴碰撞壁反馈产生的压力扰动波相互作用形成[１３Ｇ１５].本文基于自激振荡脉冲喷嘴空化效应和多相流模型,建立自激振荡脉冲射流空化模型和腔室轴对称物理模型,计算得到自激振荡脉冲射流空化泡破碎㊁腔室内两相分布㊁湍动能分布和速度分布等结果.１㊀数学模型１．１㊀控制方程自激振荡脉冲射流空化效应产生机理如图１所示.高速射流中的离散涡在自激振荡腔室内的剪切层中被选择性放大,形成大尺度涡旋结构,进而形成沿腔室轴线对称分布的空化气囊.该气囊对喷嘴入口来流产生周期性的能量聚集与释放,使连续射流转变为脉冲射流,并使射流具有压力波动和一定的空化效应;同时,自激振荡腔室内存在着的大尺度漩涡以及脉冲压力振荡效应,将导致自激振荡腔室内漩涡空化和振荡空化的形成,从而加强喷嘴的空化效果.为研究空化效应下的自激振荡脉冲射流形态,对空化泡破碎㊁自激振荡脉冲腔室内场及外场进行数值分析.由于计算模型中空化泡破碎涉及空化泡㊁液体和空气两相作用,喷嘴内流场计算涉及空化现象,喷嘴外流场计算则涉及空气和射流的相互作用,故整个射流形态计算模型采用多相流模型同空化模型进行计算.图１㊀自激振荡脉冲空化效应F i g．１㊀T h e c a v i t a t i o n e f f e c t o f s e l fＧe x c i t e do s c i l l a t i o n p u l s e根据多相流计算模型,建立混合相的连续性方程和动量方程[１６]:∂ρ∂t＋Ñ(ρu)＝０(１)∂(ρu)∂t＋Ñ(ρu u)＝－Ñp＋Ñτ＋ʏδ(t)σk nδ(x－xᶄ)d S(２)式中,t为时间;u为速度;σ为表面张力系数;ρ为混合密度;k为界面曲率;p为压力;τ为黏性切应力;n为指向表面S的法向单位矢量;δ(x)为D i r a c函数.各相的体积分数满足６３５１中国机械工程第２８卷第１３期２０１７年７月上半月Copyright©博看网 . All Rights Reserved.(∂αi/∂t)＋Ñ(uαi)＝０㊀i＝l,g,v(３)αl＋αg＋αv＝１(４)ρ＝αlρl＋αgρg＋αvρv(５)μ＝αlμl＋αgμg＋αvμv(６)式中,αi为体积分数;αl为液相体积分数;αg为蒸汽相体积分数;αv为空气相体积分数;ρl为液相密度;ρg为蒸汽相密度;ρv为空气相密度;μ为混合相黏度;μl为液相黏度;μg为蒸汽相黏度;μv为空气相黏度.式(２)中的界面曲率k可根据所求解的液相体积分数αl估算,即k＝Ñ(Ñαl/|Ñαl|)(７)１．２㊀空化模型射流空化形成涉及相变传质过程,质量传输是建立空化模型所要考虑的关键问题.空化传质模型建立并不是独立的,空化模型将作为平衡方程的一部分用来描述蒸汽的产生与破碎.当流场内部压力低于饱和蒸汽压时产生蒸汽,平衡方程的其他部分解决流动和射流问题.本文采用Z w a r tＧG e r b e rＧB e l a m r i空化传质模型[１７],模型假设液体中所有的气泡具有相同的初始尺寸,采用单位体积内气泡的数目计算传质效率:R＝n(４πR２B p v d R B d t)(８)式中,n为单位体积内气泡数目;R B为气泡直径;p v为气泡内压力.忽略R a y l e i g hＧP l e s s e t方程中的二阶导数项㊁黏性项以及表面张力项可得d R B d t＝２(p B－p f)３ρl(９)式中,p B为气泡内的压力;p f为非凝结气体的部分压力.将式(９)代入式(８),可得由体积分数表示的传质效率R＝３αρvR B２(p B－p f)３ρl(１０)最终气泡蒸发与凝结的传质效率表达形式为R e＝C e３αn u c(１－αv)R Bρv ２(p v－p f)３ρl(１１)R c＝C c３αv R Bρv２(p v－p f)３ρl(１２)式中,αn u c为气核体积分数;C e为蒸发常数相;C c为凝结常数相.２㊀物理模型２．１㊀腔室结构及计算参数基于自激振荡脉冲效应的喷嘴腔室几何模型如图２所示.喷嘴结构及其主要计算参数为[１,１８]:上喷嘴入口流道长度l１＝３０mm;下喷嘴出口流道长度l２＝３０mm;上喷嘴入口流道直径d１＝８mm;下喷嘴出口流道直径d２＝１５mm;自激振荡脉冲腔室直径D＝１００mm;自激振荡脉冲腔室长度L＝６０mm;下喷嘴碰撞壁夹角α＝１２０ʎ.图２㊀自激振荡脉冲腔室结构F i g．２㊀T h e c h a m b e r s t r u c t u r e o f s e l fＧo s c i l l a t i n gp u l s e由于需要考虑外流场空气对射流的影响,故喷嘴外部构造的轴对称计算域几何模型如图３所示.计算过程中发现:当外流场设置的空气域模型是喷嘴腔室直径的３倍之后,继续增大计算区域,计算的结果基本不发生变化;同时,计算区域增大,计算量也会增大,也没有实际计算意义.由此设置外流场计算域几何模型的空气域尺寸为喷嘴腔室直径的３倍.图３㊀外流场计算域几何模型F i g．３㊀T h e g e o m e t r y m o d e l a n d c o m p u t a t i o n a ld o m a i no f t he o u tf l o wf i e l d２．２㊀有限元求解本文采用G AM B I T软件对计算区域进行网格划分,并使用F L U E N T软件对网格计算域进行求解.由于自激振荡腔室结构尺寸与外部空气域相差悬殊,故本文在划分网格时首先采用线网格划分,然后再进行面网格划分,这样可以保证自激振荡腔室内部以及靠近喷嘴出口处的网格密度较大,而远离该区域的网格密度较小.网格划分及计算区域边界条件如图４所示,其中入口１为液体压力入口,参数设置为１１０１３２５P a;入口２为空气压力入口,参数设置为１０１３２５P a;压力出口参数设置为１０１３２５P a.在计算过程中,对流体流动状态及相关计算模７３５１自激振荡脉冲喷嘴空化效应及其射流形态的数值分析 汪朝晖㊀胡亚男㊀饶长健等Copyright©博看网 . All Rights Reserved.图４㊀网格划分及计算区域边界条件F i g．４㊀T h e g i r dd i v i s i o na n d t h e b o u n d a r y c o n d i t i o no f c o m p u t a t i o n a l d o m a i n型进行如下设置:液体介质选用常温状态下的水,密度为１０００k g/m３,动力黏度为１０－３N s/m;气体介质选用常温状态下的空气,密度为１．２２５k g/m３,动力黏度为１．７９ˑ１０－５N s/m;液体饱和蒸汽压设置为３５４０P a;流体流动状态采用层流模型;计算过程采用瞬态模型;湍流模型选用可实现的kＧε模型.３㊀计算结果与分析３．１㊀空化泡破碎由于自激振荡脉冲喷嘴出口空化区域内的压力较外界大气压力小,因此空化泡由喷嘴出口进入大气环境时,空化泡外壁压力增大,在空化泡内部,相同温度条件下的饱和蒸汽压力和气体组分压力保持恒定.随着自激振荡脉冲射流空化泡半径的不断增大,对于某一确定初始半径的空化泡,存在着最大空化泡极限破碎半径,当空化泡半径达到极限值时空化泡破碎.单个空化泡在自由液面附近的破碎过程如图５所示.可以看出,在１．０２~２．３７m s时,空化泡开始溃灭,其半径开始急速减小,空化气囊开始径向运动,并且液体被带动向内挤压,形成一个泡面的加速射流.在２．６９~４．６７m s时,空化气囊加速径向运动,挤压射流的压力开始急速增大,当挤压压力达到空化泡极限破碎值时,空化泡从泡面最上端开始破碎,并且液体从空化泡面喷射.由于空化泡溃灭伴随着能量释放,因此在空化泡破碎过程中液体的湍动能增大.空化泡溃灭对自由液面附近的液体产生了破碎效果,因而当空化射流中出现大面积的空化泡溃灭时,喷射射流将会出现更加明显的雾化效果.３．２㊀两相分布根据伯努利方程计算得到喷嘴入口射流流速为４４．７６m/s,据此计算得到的自激振荡脉冲腔室频率大约为１１H z,因此本文的自激振荡脉冲腔室振荡周期为１００m s.下面分析在此一个振荡(a)０．１m s㊀㊀㊀㊀㊀㊀(b)１．０２m s(c)２．３７m s㊀㊀㊀㊀㊀㊀(d)２．６９m s(e)３．４９m s㊀㊀㊀㊀㊀㊀(f)４．６７m s图５㊀单个空化泡在自由液面附近的破碎过程F i g．５㊀T h e b r o k e n p r o c e s s o f t h e s i n g l e c a v i t a t i o nb u b b l e n e a r t h e f r e e l i q u i d s u r f ac e周期内不同时刻的两相分布.振荡周期内不同时刻的液态水相㊁蒸汽相和空气相分布如图６所示.(a)１m s㊀(b)５m s㊀(c)７．５m s㊀(d)１０m s㊀(e)２５m s(f)４０m s㊀㊀(g)５５m s㊀(h)９０m s㊀(i)１００m s㊀图６㊀振荡周期内不同时刻的两相分布F i g．６㊀T h e t w oＧp h a s e d i s t r i b u t i o na t t h e d i f f e r e n tm o m e n t s i n t h e o s c i l l a t i o n p e r i o d由图６可看出,１~１０m s为自激振荡腔室内空化气囊的形成和生长过程.空化最初出现在腔室内的分离区,伴随着射流逐渐向下喷嘴运动,分离区内的空化泡也随之运动并不断生长.当空化泡到达下喷嘴碰撞壁后,受阻的空化泡开始沿碰撞壁壁面反向运动.在２５m s时,空化泡反向运动停止,在该位置空化泡持续生长,此时的射流未被阻断,液相体积分数在出流管道内达到最大.到４０m s时,空化气囊与主射流相互作用,这种作用在腔室内表现为气液相的混合,在腔室外表现为射流形态的改变.在４０~１００m s阶段可以看到,外流场的主射流逐渐变窄,两侧的液相体积分数开始减小并逐渐趋近于０.主射流逐渐被截断,空化泡与下喷嘴碰撞壁碰撞后开始反方向运８３５１中国机械工程第２８卷第１３期２０１７年７月上半月Copyright©博看网 . All Rights Reserved.动直至腔室固定位置后停止运动,空化泡达到体积最大值.此时,主射流被中心空化气囊所阻挡,形成了射流间断,此过程持续至空化射流振荡周期结束.此阶段,出流管道内的液相体积分数变化最为复杂,它既有空气的卷吸效应,又有空化气囊在管道内的破碎.可以看出:腔室内空化气囊对射流形态产生的影响主要表现在喷雾形态的变化,主射流液柱与空气交界面出现不规则锯齿状,造成连续射流被割断成多股状射流,加剧了空化射流的扰动和不稳定性,进而促进射流雾化的产生.３．３㊀湍动能分布一个振荡周期内不同时刻的湍动能分布如图７所示.在５m s 时,腔室内部空化气囊在靠近下喷嘴位置,由于空化气囊与下喷嘴碰撞壁发生了相互作用,因而在该区域出现较强的湍动能,主射流刚进入外流场时引起了空气场扰动,所以即使在射流未穿透区域也依然可以看到比较强烈的湍(a )５m s ㊀㊀㊀㊀㊀㊀㊀㊀㊀(b )１０m s(c )２５m s ㊀㊀㊀㊀㊀㊀㊀㊀㊀(d )４０m s(e )６０m s ㊀㊀㊀㊀㊀㊀㊀㊀㊀(f )９０m s图７㊀振荡周期内不同时刻的湍动能分布F i g ．７㊀T h e t u r b u l e n t k i n e t i c e n e r g y di s t r i b u t i o na t t h e d i f f e r e n t t i m e i n t h e o s c i l l a t i o n p e r i o d动能;在１０m s 时,同样可以看到,在腔室内由于空化气囊逐渐向两侧壁面移动并在边角区域内与壁面相互作用,使得该区域内的湍动能有所增强,在外流场中,主射流与空气相互作用,在其接触的边界面形成较强湍动能,这说明该区域内空气与液柱发生了较为强烈的动量交换;在２５m s 时,主射流已经完全贯穿外流场,完整的射流形态形成,此时,在主射流与空气交界面上的湍动能继续增强,更有利于较好射流雾化效果的形成;在４０m s 时,中心空化气囊已经完全形成并开始阻挡主射流运动,由于气囊未与腔室壁面发生相互作用,所以腔室内空化区域的湍动能较弱,同时由于主射流被阻挡,喷嘴靠近出口处的湍动能较强,而外流场湍动能开始逐渐减弱,即说明在该时刻雾化效果较强的区域集中在喷嘴出口附近;在９０m s 时,腔室内空化气囊很大程度上阻挡了主射流运动,同时自激振荡喷嘴出口流道内也开始出现了大面积的空化区域,湍动能最强区域集中在下喷嘴出口下游,此时相应的雾化效果也在该区域内最强.３．４㊀速度分布一个振荡周期内不同时刻的内流场速度分布如图８所示.１m s 时,在自激振荡腔室的分离区形成涡结构,空化最先出现在漩涡中心位置,这是由于漩涡中心出现了低压区,一方面低压区内压力达到饱和蒸汽压时部分液体汽化,另一方面低压区的出现使得溶解于水体中的部分气体释放.在１０~５５m s 时,腔室中心漩涡区开始逐渐变大,直至将近布满整个半腔室;同时,由于腔室壁面进行了重构,故腔室四周边角区域不再出现次生涡.(a )１m s ㊀㊀㊀㊀㊀㊀㊀㊀(b )１０ms(c )２５m s ㊀㊀㊀㊀㊀㊀㊀㊀(d )５５m s图８㊀振荡周期内不同时刻的内流场速度分布F i g ．８㊀T h e v e l o c i t y di s t r i b u t i o no f t h e i n t e r n a l f l o wf i e l d a t t h e d i f f e r e n tm o m e n t s i n t h e o s c i l l a t i o n p e r i o d当射流由喷嘴内部进入空气时,由于高压水射流与外界大气相对速度较大,故两者的速度大小和方向不同造成了气液相的相互作用.一个振荡周期内不同时刻的外流场速度分布如图９所示.在１０m s 时,射流刚进入大气环境,该阶段内的气液相相互作用最为明显.高速运动的液柱带动处于静止状态的气体,气体流线开始偏向主射流区,这一阶段内的液柱形状也因受到气体的作用而不稳定.在射流运动过程中,外界气体运动９３５１ 自激振荡脉冲喷嘴空化效应及其射流形态的数值分析汪朝晖㊀胡亚男㊀饶长健等Copyright©博看网 . All Rights Reserved.速度开始趋于稳定.在５５m s 时,主射流区外的气体流线开始逐渐与主射流贴合.达到９０m s时,主射流与空气的相互作用基本达到稳定状态.(a )１０ms(b )５５ms(c )９０m s图９㊀振荡周期内不同时刻的外流场速度分布F i g ．９㊀T h e v e l o c i t y d i s t r i b u t i o no f t h e o u t f l o wf i e l d a t t h e d i f f e r e n tm o m e n t s i n t h e o s c i l l a t i o n p e r i o d一个振荡周期内不同时刻的外流场出口速度分布如图１０所示.横轴表示外流场出口边界与轴线的横向距离,纵轴表示边界面上某点的速度大小.在５m s 时,靠近轴线的外流场流体速度最高值达到了近９０m /s ,但是该速度并非射流的喷射速度,而是射流在没有达到稳定状态前由主射流而引起的空气湍动能增强的结果.伴随着主射流与空气相互作用逐渐稳定以及射流贯穿距离增加,主射流速度逐渐稳定.１０m s 为过渡时刻.随后射流在喷射轴线附近的速度稳定在３０~４０m /s.在４０m s 后,距轴线约１８０mm 处,射流速度开始下降并最终趋于０,这说明在一个振荡周期末射流的速度扩散作用减弱,主射流角度减小.图１０㊀振荡周期内不同时刻的外流场出口速度曲线F i g ．１０㊀T h e o u t l e t v e l o c i t y cu r v e o f t h e o u t f l o wf i e l d a t t h e d i f f e r e n t t i m e s i n t h e o s c i l l a t i o n p e r i o d４㊀结论(１)随着自激振荡脉冲射流空化泡溃灭,其半径急速减小,空化气囊开始径向运动,形成泡面加速射流.当挤压射流压力急速增大并达到空化泡极限破碎值时,空化泡从泡面最上端开始破碎.当空化射流中出现大面积的空化泡溃灭时,射流将会出现明显的雾化效果.(２)空化最初出现在自激振荡腔室的分离区,伴随着射流逐渐向下喷嘴运动,分离区内的空化泡也随之运动并不断生长.当外流场的主射流逐渐变窄时,主射流被截断,空化泡达到体积最大值,使连续射流被割断成多股状射流,加剧了空化射流的扰动和不稳定性.(３)空化气囊与下喷嘴碰撞壁发生相互作用时,出现较强湍动能.当主射流完全贯穿外流场时,主射流与空气交界面上的湍动能继续增强;在射流振荡周期末,腔室内空化气囊阻挡了主射流运动,且在喷嘴出口流道出现大面积的空化区域,湍动能最大区域集中在下喷嘴出口下游.(４)射流振荡初期气液相相互作用明显,伴随着主射流与空气相互作用逐渐稳定以及射流贯穿距离增加,主射流速度逐渐趋于稳定.在振荡周期末射流的速度扩散作用减弱,主射流角度减小.参考文献[１]㊀汪朝晖,胡亚男,廖振方,等．基于自激振荡脉冲效应的雾化喷嘴出口流道空化特性研究[J ]．机械工程学报,２０１６,５２(１４):２０４Ｇ２１２．WA N GZ h a o h u i ,HU Y a n a n ,L I A O Z h e n f a n g,e t a l ．C a v i t a t i o n C h a r a c t e r i s t i c S t u d y o n t h e O u t l e t C h a n n e l o fA u t o m i z a t i o n N o z z l eB a s e do nt h eS e l f Ｇe x c i t e dO s c i l l a t i n g P u l s eEf f e c t s [J ]．J o u r n a l o fM e Ｇc h a n i c a l E ng i n e e r i n g ,２０１６,５２(１４):２０４Ｇ２１２．[２]㊀MA R ,S L A B O C H PE ,MO R R I SSC ．F l u i d M e Ｇc h a n i c so ft h e F l o w Ｇe x c i t ed He l m h o l t z R e s o n a t o r０４５１ 中国机械工程第２８卷第１３期２０１７年７月上半月Copyright©博看网 . All Rights Reserved.[J]．J o u r n a l o fF l u i d M e c h a n i c s,２００９,６２３:１Ｇ２６．[３]㊀S UH H K,L E ECS．E f f e c t o f C a v i t a t i o n i nN o z z l e O r i f i c e o n t h eD i e s e lF u e lA t o m i z a t i o nC h a r a c t e r i sＧt i c s[J]．I n t e r n a t i o n a lJ o u r n a lo f H e a ta n d F l u i dF l o w,２００８,２９(４):１００１Ｇ１００９．[４]㊀A L E HO S S E I N H,Q I N Z．N u m e r i c a lA n a l y s i so f R a y l e i g hＧp l e s s e tE q u a t i o n f o rC a v i t a t i n g W a t e r J e t s[J]．I n t e r n a t i o n a l J o u r n a l f o rN u m e r i c a lM e t h o d s i nE n g i n e e r i n g,２００７,７２(７):７８０Ｇ８０７．[５]㊀Q I NZ,B R E MH O R S T K,A L E H O S S E I N H,e t a l．S i m u l a t i o n o f C a v i t a t i o n B u b b l e si n a C o n v e r g e n tＧd i ve r g e n tN o z z l eW a t e r J e t[J]．J o u r n a l of F l u i dM eＧc h a n i c s,２００７,５７３:１Ｇ２５．[６]㊀姚立明,赵怡,李大蔚,等．不同深度对空化射流形态影响的数值研究[J]．水动力学研究与进展A辑,２０１４,A２９(６):６７５Ｇ６８２．Y A O L i m i n g,Z HA O Y i,L I D a w e i,e t a l．N u m e r iＧc a lR e s e a r c ho f t h eE n v i r o n m e n t a l I n f l u e n c eo fD i fＧf e r e n tD e p t ho n t h e S h a p e o f C a v i t a t i o n J e t[J]．C h iＧn e s e J o u r n a l o fH y d r o d y n a m i c s,２０１４,A２９(６):６７５Ｇ６８２．[７]㊀C H E N H S,L I J,C H E N D,e t a l．D a m a g e s o nS t e e l S u r f a c ea tt h eI n c u b a t i o n S t a g e o ft h e V i b r a t i o nC a v i t a t i o nE r o s i o n i n W a t e r[J]．W e a r,２００８,２６５(５/６):６９２Ｇ６９８．[８]㊀刘琦,欧阳光耀,杨昆,等．高压喷射条件下非常态燃油喷嘴内部空化流动特性[J]．农业机械学报,２０１６,４７(６):３３３Ｇ３４０．L I U Q i,O U Y A N G G u a n g y a o,Y A N G K u n,e t a l．N o z z l e I n n e r C a v i t a t i o n F l o w C h a r a c t e r i s t i c s o fN o nＧn o r m a l F u e lB a s e do n H i g hP r e s s u r e I n j e c t i o nC o n d i t i o n[J]．T r a n s a c t i o n so ft h eC h i n e s eS o c i e t yf o rAg r i c u l t u r a lM a chi n e r y,２０１６,４７(６):３３３Ｇ３４０．[９]㊀X I E W F,L I U T G,K HO O BC．T h eS i m u l a t i o n o fC a v i t a t i n g F l o w sI n d u c e db y U n d e r w a t e rS h o c ka n dF r e eS u r f a c e I n t e r a c t i o n[J]．A p p l i e dN u m e r i c a lM a t h e m a t i c s,２００７,５７(５):７３４Ｇ７４５．[１０]㊀WA N GG,Z HA N GS,Y U M,e t a l．I n v e s t i g a t i o n o f t h eS h o c k W a v eP r o p a g a t i o nC h a r a c t e r i s t i c s a n dC a v i t a t i o n E f f e c t so f U n d e r w a t e rE x p l o s i o n N e a rB o u n d a r i e s[J]．A p p l i e dO c e a nR e s e a r c h,２０１４,４６:４０Ｇ５３．[１１]㊀曾宇杰,罗坤,樊建人,等．主动脉夹层血液两相流动数值模拟分析[J]．工程热物理学报,２０１６,３７(４):７８０Ｇ７８４．Z E N G Y u j i e,L U O K u n,F A N J i a n r e n,e ta l．A o r t i cD i s s e c t i o nB l o o dT w oＧp h a s eF l o w N u m e r iＧc a lS i m u l a t i o n S t ud y[J]．J o u r n a lo fE n g i ne e r i n gT h e r m o p h y s i c s,２０１６,３７(４):７８０Ｇ７８４．[１２]㊀王维军,王洋,刘瑞华,等．离心泵空化流动数值计算[J]．农业机械学报,２０１４,４５(３):３７Ｇ４４．WA N G W e i j u n,WA N GY a n g,L I UR u i h u a,e t a l．N u m e r i c a l C a l c u l a t i o no f C a v i t a t i o nF l o wi n aC e nＧt r i f u g a l P u m p[J]．T r a n s a c t i o n s o f t h eC h i n e s eS oＧc i e t y f o rA g r i c u l t u r a lM a c h i n e r y,２０１４,４５(３):３７Ｇ４４．[１３]㊀L I A O ZF,L I J,C H E N D S,e ta l．T h e o r y a n dE x p e r i m e n t a l S t u d y o f t h eS e l fＧe x c i t e dO s c i l l a t i o nP u l s e d J e tN o z z l e[J]．C h i n e s e J o u r n a l o fM e c h a n iＧc a l E n g i n e e r i n g,２００３,１６(４):３７９Ｇ３８３．[１４]㊀S H E R V A N IＧT A B A R M T,P A R S A S,G HO RＧB A N I M．N u m e r i c a lS t u d y o nt h e E f f e c to ft h eC a v i t a t i o nP h e n o m e n o no nt h e C h a r a c t e r i s t i c so fF u e l S p r a y[J]．M a t h e m a t i c a l a n dC o m p u t e rM o dＧe l l i n g,２０１２,５６(５/６):１０５Ｇ１１７．[１５]㊀L ID,K A N G Y,D I N G X,e t a l．E f f e c t so fA r e aD i s c o n t i n u i t y a tN o z z l e I n l e t o n t h eC h a r a c t e r i s t i c so fH i g hS p e e dS e l fＧe x c i t e dO s c i l l a t i o nP u l s e d W aＧt e r J e t s[J]．E x p e r i m e n t a lT h e r m a l&F l u i dS c iＧe n c e,２０１６,７９:２５４Ｇ２６５．[１６]㊀王福军．计算流体力学分析ＧC F D软件原理与应用[M]．北京:清华大学出版社,２００５:７Ｇ１３．WA N G F u j u n．C o m p u t a t i o n a l F l u i d D y n a m i c sA n a l y s i s:C F DS o f t w a r eP r i n c i p l ea n dA p p l i c a t i o n[M]．B e i j i n g:T s i n g h u aU n i v e r s i t y P r e s s,２００５:７Ｇ１３．[１７]㊀Z WA R TPJ,G O D I NPQ,P E N R O S EJ,e t a l．A T w oＧp h a s eF l o w M o d e lf o rP r e d i c t i n g C a v i t a t i o nD y n a m i c s[C]//F i f t hI n t e r n a t i o n a lC o n f e r e n c eo nM u l t i p h a s eF l o w．Y o k o h a m a,２００４:１Ｇ１１．[１８]㊀王循明,焦磊,王乐勤．自激脉冲发生机理数值模拟及参数影响分析[J]．浙江大学学报(工学版),２００５,３９(９):１４５０Ｇ１４５４．WA N G X u n m i n g,J I A O L e i,WA N G L e q i n．N uＧm e r i c a l S i m u l a t i o n o f S e l fＧe x c i t e d O s c i l l a t i o nP u l s e d J e t a n d A n a l y s i so fP a r a m e t e r sI n f l u e n c e[J]．J o u r n a lo fZ h e j i a n g U n i v e r s i t y(E n g i n e e r i n gS c i e n c e),２００５,３９(９):１４５０Ｇ１４５４．(编辑㊀袁兴玲)作者简介:汪朝晖,男,１９８１年生.武汉科技大学机械自动化学院教授㊁博士研究生导师.主要研究方向为计算流体动力学及其装备设计理论与方法.EＧm a i l:z h w a n g＠w u s t．e d u．c n.胡亚男,男,１９９１年生.武汉科技大学机械自动化学院硕士研究生.饶长健,男,１９９３年生.武汉科技大学机械自动化学院硕士研究生.邓晓刚,男,１９７５年生.重庆科技学院机械与动力工程学院教授.１４５１自激振荡脉冲喷嘴空化效应及其射流形态的数值分析 汪朝晖㊀胡亚男㊀饶长健等Copyright©博看网 . All Rights Reserved.。
