空化紊流流动的数值计算模型及其验证

本文介绍了一种空化紊流流动的数值计算模型,并计算了射流放水阀内部的空化流动现象。该模型基于均相平衡模型和液相与汽相传输方程,基本方程采用N-S方程,空化模型采用Kunz等提出的汽液质量转换模型。紊流封闭采用标准的κ-ε素流模型;稳定流动计算采用扩展的SIMPLE压力修正方法,非稳定计算采用PISO算法。为评价该数值模型,计算了绕射流放水阀的空化流动,并与实验结果进行了对比,计算结果与实验结果取得较好的一致,说明该流动模型和计算方法是可靠的。     

基于几何重建的气液两相流孔隙尺度模拟

本文基于二维micro CT扫描图像进行数值重构得到三维多孔介质几何模型,采用Shan-Chen多相格子Boltzmann方法(LBM)建立适用不规则几何边界区域内气液两相流动数值模型,在此基础上对复杂多孔介质内的两相流动过程进行孔隙尺度模拟。结果表明,气-液两相流动受限于多孔介质的复杂孔隙结构,气液表面张力、壁面润湿特性以及进出口压差对两相流动特性将产生明显影响。     

三维小空化数空泡流数值方法

采用均质平衡流空泡模型,引入基于求解液体质量份数输运方程的空泡流模型,提出新的压力密度耦合的压缩性方法和控制气液两相转换的源项,求解汽水混合介质的RANS方程和带低雷诺数修正的k-ε模型,实现了小空化数(σ=0.2~0.01)下水下航行体的空泡流数值模拟.得到清晰的空泡形态特征与内部结构,以及空泡长度和最大直径随空化数的变化规律,给出航行体运动阻力系数与空化数之间的变化关系.将计算结果与解析结果及实验数据进行比较表明,该方法保证了极小空化数下自然空泡流计算的稳定性与收敛性,提高了空泡形态特性的预报精度.     

部分流低温流体循环泵空化性能预测

为了对部分流低温流体(液氮) 循环泵空化特性进行预测. 基于软件 ANSYS-FLUENT, 计算选用Stardard k-ε 湍流模型,Simplec 压力耦合方式, 进行空化前 MRF 模型定常计算, 添加两相流参数, 选取 Singhal-et-al 空化模型, 添加两相参数, 考虑液-气密度比对低温流体液氮泵内能量的传递和交换的影响, 得到气液动量、 质量和能量守恒方程, 利用 RNGk-ε 湍流模型,Simplec 压力耦合方式, 对不同进液压力条件下, 部分流低温流体(液氮) 循环泵空化特性进行全流域空化数值计算. 进行泵空化特性试验, 在额定转速下, 随着泵前流体压力的降低呈现的空化性能, 数值计算与试验测得的泵头数值最大偏差在10% 以内, 曲线吻合性较好, 泵内流场空化发生伴有显著的压头下降, 空化过程增强, 空化区变大, 液相和汽相相互拖拽能力增强, 空化核心区由叶顶背压部分扩散到整个流道, 汽液界面不清晰, 直至断流. 本文采用的计算方法和研究结果为低温流体循环泵内部流体空化的诊断和性能优化提供了一定依据.     

绕水翼片状空化流动结构的数值与实验研究

采用数值与实验相结合的方法研究了水翼片状空化流动结构.实验采用高速录像技术观察了片状卒泡形态,应用LDV测量了翼型周围的湍动能和速度分布;采用N-S方程和基于空泡动力学方程的空化模型计算了绕水翼片状空化流场.结果表明:在片状空化阶段,翼型吸力面上附着很薄的一层透明空泡,空泡彤态呈现于指状;随着空化数的减少,空泡尾部水汽交界面相互作用增强,并且空泡尾部出现大的旋涡,影响了空泡尾部区域压力和速度分布,片状空泡尾部的水汽混合区出现不稳定现象,同时存在小的空泡团脱落.数值模拟得到的水翼片状空化流动现象和实验观察到的结果基本一致,验证了计算模型和数值方法的可靠性.     

零净液流量两相流持液率与阻力特性研究

分别以牛顿流体和非牛顿流体为液相，研究了垂直管中零净液流量气液两相流的流动特性。提出了零净液流量气液两相流动模型，应用这一模型计算了零净液流量气液两相流的持液率和压力降，模型计算结果与试验结果相符。研究结果表明，零净液流量气液两相流与常规气液两相流相比具有特殊性，表现为其持液率仅由质量平衡方程控制，其摩擦阻力压力降为负值。     

基于双重介质模型的正仲氢转化器流动规律研究

为了研究正仲氢转化器中流体的流动规律，建立了催化剂堆积床的颗粒孔隙-堆床空隙的双重介质模型。基于双重介质模型分别建立颗粒孔隙和堆床空隙中正氢和仲氢的质量守恒方程，两种介质之间通过流量交换源项相互耦合。采用控制容法离散控制方程，并数值求解离散方程组。数值算例表明，由于颗粒孔隙与堆床空隙的差异，两种介质中的流动规律存在明显差异。     

混流泵叶轮内空化流动的数值计算

本文应用三维Navier-Stokes方程和k-ε两方程湍流模型以及全空化模型,对一台喷水推进混流泵叶轮内的空化流动进行数值模拟.介绍了空化流动的模拟方法,并计算了混流泵的空化性能指标,有效地预测了叶片上空化发生的区域和空化流动的发展情况.文中还对两种进口相对流动角作了数值计算,结果表明小的相对流动角可以提高混流泵的空化性能,但却使扬程有所下降.     

涡流二极管空化流动的数值计算

为了研究涡流二极管内空化现象的机理特性及对其性能的影响,我们假设流体为气液混相均质,并考虑不可凝结气相,采用基于组分输运方程,求解了涡流二极管全流道内气液混相均质流的雷诺平均N-S方程以及气相组分输运方程。数值计算结果显示了空化形成时涡流二极管入口、出口及旋流腔内的流场形态,研究表明:涡流二极管空化现象主要发生在流体切向进入时旋流腔和中心管的中心部位;空化流是由于液体在中心旋流场低压条件下汽化,同时不可凝结气体由于亨利定律在旋转流场形成的压力梯度下而发生的输运效应综合形成的;空化流由于强旋涡的原因对涡流二极管的性能产生明显的影响。上述结论对涡流二极管的设计及其指导工程应用具有重要的价值。     

移动网格与Level Set耦合方法在气-液两相流数值模拟中的应用

本文采用自适应移动网格与Level Set函数相耦合的方法来实现气-液两相流的数值模拟与计算.作为自适应网格方法的一种,移动网格方法主要是为了解决发展方程的计算问题而设计的方法.文中给出了移动网格的生成方程,并针对方程的非线性,给出了一种半隐式的离散方法用于进行求解.本文将移动网格方法与Level Set方法相耦合,将控制流体运动的Navier-Stokes方程以及追踪相界面的Level Set方程转换到曲线坐标下,应用一套曲线坐标方程组来同时描述气、液两相流的运动规律,成功实现了对气-液两相流问题的数值模拟.通过对顶盖驱动流的计算以及对液滴沉降现象的模拟计算,验证了本文方法的可靠性.本文对常重力与微重力下两气泡融合的发展规律进行了数值模拟,通过分析对比,得到了重力对两气泡融合变形的影响规律.     

Analysis of Two-Phase Cavitating Flow with Two-Fluid Model Using Integrated Boltzmann Equations

In the present work, both computational and experimental methods are employed to study the two-phase flow occurring in a model pump sump. The two-fluid model of the two-phase flow has been applied to the simulation of the three-dimensional cavitating flow. The governing equations of the two-phase cavitating flow are derived from the kinetic theory based on the Boltzmann equation. The isotropic RNG$k-\epsilon-k_{ca}$ turbulence model of two-phase flows in the form of cavity number instead of the form of cavity phase volume fraction is developed. The RNG $k-\epsilon-k_{ca}$ turbulence model, that is the RNG$k-\epsilon$ turbulence model for the liquid phase combined with the $k_{ca}$model for the cavity phase, is employed to close the governing turbulent equations of the two-phase flow. The computation of the cavitating flow through a model pump sump has been carried out with this model in three-dimensional spaces. The calculated results have been compared with the data of the PIV experiment. Good qualitative agreement has been achieved which exhibits the reliability of the numerical simulation model.     

基于小关联时间两相湍流色噪声的近似模型及数值模拟

从一般高斯型色噪声模型出发,通过泛函导数,应用小关联时间,近似计算多维色噪声,得到有效Fokker-Planck方程.将其应用到两相湍流中得到颗粒相的概率密度函数输运方程,从而得到颗粒相的二阶矩模型.将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,并将求解结果与实验结果进行对比分析.     

轴流式水轮机内部空化流动的数值预测

将一种完整空化模型和一种混合流体两相流模型相结合,应用于轴流式水轮机全流道内定常空化流动的数值模拟。根据模拟结果,分别预测了水轮机在正常和较低尾水位情况下运行时,流道内空化发生的部位和程度,重点考察了尾水位的降低对水轮机空化性能的影响,为水轮机在特定运行工况下的空化性能预测提供参考。     

用视密度加权平均二阶矩模型模拟旋流两相流动

本文用视密度加权平均代替时平均,建立了视密度加权平均的统一二阶矩两相湍流模型方程组（ＭＵＳＭ）,其中用体积分数代替了数密度,用颗粒驰豫时间作为封闭两相脉动速度关联方程耗散项的时间尺度,并引入了颗粒视在的气体速度脉动的输运方程。用ＭＵＳＭ模型模拟了旋流数为０．４７的气粒两相流动。并和实验结果及时间平均的ＵＳＭ模型的模拟结果进行了对照,两种模型均能较好地预报的两相的轴向和切向速度,轴向和切向脉动速度。此外,ＭＵＳＭ模型可以减少所用方程数,节省计算量。因此视密度加权平均的统一二阶矩两相湍流模型是一种对时间平均的统一二阶矩模型的改进,今后可以进一步扩大应用。     

Numerical simulation of pulsation processes in hydraulic turbine based on 3D model of cavitating flow

A new approach was proposed for simulation of unsteady cavitating flow in the flow passage of a hydraulic power plant. 1D hydro-acoustics equations are solved in the penstock domain. 3D equations of turbulent flow of isothermal compressible liquid-vapor mixture are solved in the turbine domain. Cavitation is described by a transfer equation for liquid phase with a source term which is responsible for evaporation and condensation. The developed method was applied for simulation of pulsations in pressure, discharge, and total energy propagating along the flow conduit of the hydraulic power plant. Simulation results are in qualitative and quantitative agreement with experiment. The influence of key physical and numerical parameters like discharge, cavitation number, penstock length, time step, and vapor density on simulation results was studied.     

On the Computations of Gas-Solid Mixture Two-Phase Flow

This paper presents high-resolution computations of a two-phase gas-solid mixture using a well-defined mathematical model. The HLL Riemann solver is applied to solve the Riemann problem for the model equations. This solution is then employed in the construction of upwind Godunov methods to solve the general initial-boundary value problem for the two-phase gas-solid mixture. Several representative test cases have been carried out and numerical solutions are provided in comparison with existing numerical results. To demonstrate the robustness, effectiveness and capability of these methods, the model results are compared with reference solutions. In addition to that, these results are compared with the results of other simulations carried out for the same set of test cases using other numerical methods available in the literature. The diverse comparisons demonstrate that both the model equations and the numerical methods are clear in mathematical and physical concepts for two-phase fluid flow problems.     

A modified PANS model for computations of unsteady turbulence cavitating flows

A modification to the PANS (partially averaged Navier-Stokes) model is proposed to simulate unsteady cavitating flows. In the model, the parameter f _k is modified to vary as a function of the ratios between the water density and the mixture density in the local flows. The objective of this study is to validate the modified model and further understand the interaction between turbulence and cavitation around a Clark-Y hydrofoil. The comparisons between the numerical and experiment results show that the modified model can be improved to predict the cavity evolution, vortex shedding frequency and the lift force fluctuating in time fairly well, as it can effectively modulate the eddy viscosity in the cavitating region and various levels of physical turbulent fluctuations are resolved. In addition, from the computational results, it is proved that cavitation phenomenon physically influences the turbulent level, especially by the vortex shedding behaviors. Also, the mean u-velocity profiles demonstrate that the attached cavity thickness can alter the local turbulent shear layer.     

卧式螺旋管内汽液两相流不稳定性试验研究

本文在中、低压汽水试验台上对卧式螺旋管内汽液两相流动不稳定性进行了详细的试验研究，获得了各类脉动发生的界限及各主要参数对脉动的影响规律，并在无固次分析的基础上，给出了密度波脉动起始边界的预报关系式。     

DSM-LPDF两相湍流模型及旋流两相流动的模拟

本文由流体－颗粒速度的拉氏联合概率密度函数（ＰＤＦ）输运方程出发,用Ｓｉｍｏｎｉｎ建议的Ｌａｎｇｅｖｉｎ模型封闭颗粒所遇到流体瞬时速度的条件期望项,并用Ｍｏｎｔｅ Ｃａｒｌｏ方法直接求解 ＰＤＦ输运方程,将其和求解流体雷诺应力方程模型的有限差分方法结合,建立了雷诺应力－拉氏ＰＤＦ（ＤＳＭ－ＬＰＤＦ,简称ＤＬ）两相湍流模型．用此模型模拟了旋流数为０．４７的突扩旋流气粒两相流动,并与文献中ＰＤＰＡ实验和用类似于单相流动湍流模型封闭方法的时平均统一二阶矩（ＵＳＭ）模型的预报进行了对比．     

Turbulent gas-dispersed flow in a pipe with sudden expansion: numerical simulation

A mathematical model was developed to simulate two-phase gas-dispersed flow moving through a pipe with axisymmetric sudden expansion. In the model, the two-fluid Euler approach was used. The model is based on solving Reynolds-averaged Navier — Stokes equations for a two-phase stream. In calculating the fluctuating characteristics of the dispersed phase, equations borrowed from the models by Simonin (1991), Zaichik et al. (1994), and Derevich (2002) were used. Results of a comparative analysis with previously reported experimental and numerical data on two-phase flows with separation past sudden expansion in a plane channel and in a pipe are given. This work was supported by the President of the Russian Federation through the Foundation for Young Candidates of Sciences under Grant MK-186.2007.8 and by the Russian Foundation for Basic Research (Grants Nos. 05-08-33586 and 06-08-00967).