绕二维翼型空泡流地研究.pdfVIP

第六届全国流体力学会议论文集 2001上海 绕二维翼型空泡流的研究 吴磊 鲁侍敬 薛雷平 (上趣变通大学建筑工程与力学学院，上海2fD245) 摘要研究了黏性、非定常空泡流．应用两相流模型数值求解层流、中等Reynolds数、小攻 角二维翼型空泡流的Navicr．Stokes方程。利用水汽物理性质建立水相、汽相以及术汽相变的 密度和压力间本构方程：同时饭设密度的人工过渡区以数值模拟空泡面。本研究采用有限体 积数值方法，具有良好的向三维空泡流问题的扩展|生． 关键词两相流，—定盥一 引言 空化现象普遍存在于液体中；在流场中，通常首先在固壁附近，当流动导致局部压力低 于水(本文研究的流质)的饱和蒸汽压时．就会发生空化现象，微小的气核发展成肉眼可见的 气泡(空泡)，甚至发展为成片的附若于物体表面的空泡。在水下航行体、水下推进器、高速 泵以及喷管中可以观察到大量空化现象。空化导致固体表面剥蚀．产生噪声，空泡的存在影 响流体动力特性。 目前空泡流求解方法可主要归纳为两类：边界元法和求解Naviat-Stokes方程方法。前者 基于势流理论，在边界上布置源汇或涡，同时假设空泡尾流模型；忽略或简化黏性效应。 KinnasI“2l和Lee口I等用此方法进行了大量研究。这类方法数值模拟简单。并能在台理的闭合 模型下获得与实验数据吻合的结果：所不足的是它对空泡初始点位置、空泡尾流模型等提出 较高要求，同时三维应用较困难。后者应用两相流理论建立黏性流场的Navier-Stokes方程． 同时建立密度和压力之间的辅助关系封闭方程组。 Kubomet aI州基于Rayleigh-Plesset方程建立BTF(BubbleTwo-phaseFlow)方法：此方法适 合云雾空泡研究。Chen＆Heiste一51依据空泡内压力等于水的饱和蒸汽压R的一般假定，提出 辅助关系：Dp／D,=c(P一昂)．PB；C为非常大的常数(如5000．50000)以保证泡内 et a1161假设汽化和凝结瞬态发生，建立密度和压力的函数关系 压力等于饱和蒸汽压。Delannoy P只+△P 口= P只一△尸 (I) elsewhere 式中dm为最小声速，△P=1I 求解能量方程，依靠水汽性质的大量实验数据，获得密度、黏性、温度、热传导率以及定压 比热髓压力和焓的变化关系，求解定常或非定常空泡流问题；Reynolds数取103量级：汽相 在空泡流研究中的地位得到进一步肯定。 在本文中，我们忽略温度场的细微变化，根据水汽变化的等温变化曲线建立密度和压力 的辅助关系：在全流场建立Navier-Stokes方程，求解定常和非定常黏性空泡流。由于空泡面 附近的高马赫数，整个流场是局部可压缩的：另一方面，可压缩系数k=每p／劫在空泡面附 近(过渡区)取较高值，在其它区域取较低值或零。这样我们碍到整个流场的四个区域：泡 内区，高可压缩过渡区，低可压缩区和不可压缩区。 1 控制方程和求解过程 1 1控制方程 守恒型Navier-Stokes方程具有以下偏微分形式 挲+三(肌)=0 (2) df GXi ㈣ 导c肌，+毒(puju；)=-筹+毒[F一一；％如)] 其中e“：挈+挈 “‘ ∞J 采用积分形式并化简，可得 (4) ，I，vO(aPr