空泡群动态响应的直接数值模拟和分析.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 空泡群动态响应的直接数值模拟与分析· 陈晖1，李胜才2，左志刚2 《1两安航天动力研究所，中国两安710loo，Emajl：曼h塑血L主；@x垒hQ垒：曼Qm：鲤 2华威人学T程学院，英困考义垂CV47AL) 摘要：为了在微观层面上探索空化的动力学行为，对空泡群的动态响应进行了直接数值 模拟．与目前多数计算研究不同，该直接数值模拟技术由于采用界面追踪的方法，直接解算 果表明：外界扰动会诱发气液介质中形成压力波，而且由于空泡问的相互作用，形成了一个 具有某种低频的总体积脉动，它反映了空泡群具有整体性的动力学特性．该研究对于进一步 理解空化流路系统中低频脉动现象的机理提供了依据。 关键词：空化；数值模拟；动力学；脉动 l引言 两相流动系统中低频压力脉动的研究在很多领域十分重要【l】。20多年前，李胜才等在密 歇根人学文凸利实验台上首次发现了空化共振现象12弓J，此后各国研究人员在空泡群动力学特 得计算对象仅限于低气体体积分数(1％～2％)【4l，成为了阻碍该领域研究深入的瓶颈。本 研究采用的是一种基于界面追踪方法的直接数值模拟(DNS)技术，其理论基础是由 mggv弱onl5l教授等发展的。它的土要特点是同时采用了同定和动态两套网格，并且引入了指 示函数的概念，使得在整个运算过科中，既能保证气液界面具有一定的锐度，义能保证数值 模型对最人气体体积分数的限制，使得计算对象更接近真实。 2物理模型 为了模拟肇面附着空化层对压力扰动的动态响应，计算域采用的是一个矩形梓体，里面 充满了含有一定数鼙空泡的液体(图1)。根据前人经验，一般非定常空化云的主要部分气体体 ’英国EPSRCsWan^rick-IMRC(R．ESCM．9001，R．ESCM．9219)项目资助． 一233． 第九届全国水动力学学术会议鬻第一十二届全国水动力学研讨会文集 积分数常常任10％以上，所以在初始状态F．这些空泡的形状假最为人小辑异的球体，以 a=11％至口=27％的气体体积分数井以不同分布方式悬浮在液体中，假设住辂个过程中空 泡内的气体肝力遵守等温变化规律。 幽1计算域和边界条件 在此l×l×4柱体计算域中，四个横向边界面为周划性边界条件，使得实际研究区域在 横向是无限拓艘的；{}体底部边界条件殴置为刚体壁面；压力激励△p则是沿z向施加在该辑! 体的顶部表面(削1)。 为了分析以及推J_方便，程序中的主要参数，如计算域尺寸、空泡火小、时间、密度和 黏度等均作无量纲化处理。 3计算方法 控制^胖足1|定常N“1erStokes办拌，其守恒形式为 F(7Ⅳ+9‘Ⅳ)+腭+f 蔷(删)+V(删“)=一即十9 』℃thH 琏鹰父繁．P 密度，P～】_力：p黏度：g一重力加速度 F^一表丽jK力 引人，、指小函数，h)米表币-r泄的界面位置，这}r通过，h)就可以方便地计算山各个 网格_J．的f。J削霄瞍_}|：孙瘦值 第九届全国水动力学学术会议暨第二十一届全国水动力学研讨会文集 P(x)=p。+∞。一p。)Jb) 2) ∥(x)=∥。+∞。一比)，(z) ：3) 式中，下标b平¨o分别H{来指示空泡内的气体和空泡外的液体。 采用快速泊松求解器来求解指示函数： V2』(x)=VG0) (4 式中，Gh)一每个静ll：网格点上指示函数的梯度。 因此，每个网格点上的表面张力可以表示为： c(x)=∑Db—x“’卜“’ (5) 式中，D一分布函数，界面网格点上的参昔转移到临近静Ir网