气液两相环状流诱发管道振动模型研究

根据两相流动过程中弹性恢复力、离心力、科氏力、惯性力的平衡关系,并考虑因管道弯曲产生的非线性因素,建立了气液两相环状流诱发两端简支管道振动的模型;采用Galerkin方法对模型进行离散化处理,得到了系统的非线性矩阵方程。通过对方程系数矩阵线性部分的求解,具体分析了不同折算流速以及不同弹性模量下管道前四阶模态的量纲归一化特征值的实部及虚部之间的变化规律。结果表明:当折算液速一定时,随着折算气速的增大,管道系统的前4阶固有频率会不断减小;在不同的折算液速下,随着折算气速的增大,管道系统会呈现不同的状态;管道的弹性模量对管道的振动也有着较大的影响。通过对整个方程的求解,分析了管道横向振动位移的规律。研究结果表明:弹性模量较小时,管道的振动位移情况较为复杂;而弹性模量较大时,在同样的折算流速下,管道在平衡位置做振幅较小的周期性往复运动。     

微管内气液两相流中液膜厚度的实验测量

微通道内气液两相流中气柱(plug bubble)与通道壁之间液膜厚度的实验测量,是微热管、微流动、微电子冷却以及气泡雾化等研究中普遍关注的问题.本文利用基于光学干涉和快速傅立叶变换的空间频谱分析方法,实验测量获取了含表面活性剂水中气柱在750μm 通道内运动时其与通道壁面之间的液膜厚度.实验结果表明:表面活性剂对液膜厚度的影响比较明显,即当表面活性剂浓度在一定范围内增大时,液膜厚度会减小;此外,当气柱运动速度在一定范围内增大时,液膜厚度也会减小.     

微扰动对液体螺旋槽机械密封性能的影响

针对非接触式液体螺旋槽机械密封,分析了周期性轴向微扰动和角微扰动对密封环端面间液膜厚度的影响规律,建立了单周期螺旋槽液膜模型,采用数值法求解时变雷诺方程,研究了微扰动对液膜承载力和泄漏率等密封特性参数的影响.结果表明,膜厚对液膜承载力和泄漏率的影响显著,膜厚增加,液膜承载力减小,泄漏率增大.液膜承载力、泄漏率的变化幅度和频率主要受到轴向扰动的影响.角扰动造成沿圆周方向分布的液膜承载力不均匀而使液膜稳定性变差.     

基于湍流模型的高速螺旋槽机械密封稳态性能研究

针对高速工况下的液膜润滑螺旋槽端面机械密封，建立了其湍流润滑模型，采用有限单元法结合松弛迭代技术实现了润滑方程和液膜湍流模型的数值求解，对比分析了层流模型和湍流模型下不同螺旋槽几何参数和工况参数对密封性能的影响. 结果表明:液膜湍流效应显著提升了螺旋槽机械密封端面液膜的动力润滑效应，密封的开启力、泄漏率和刚度明显大于层流模型预测值. 在不同条件下，比较而言螺旋槽内产生更加明显的湍流效应，其内液膜流动行为远不同于层流模型. 以开启力为优化目标，湍流模型获得的优化螺旋槽几何参数在螺旋角、槽深明显不同于层流模型. 在高速和低黏度介质下，机械密封的湍流效应不可忽略.      

气液两相瞬变流的流固耦合模型研究

传统的气液两相流瞬态分析和管道动力响应计算是分开的，存在一定的缺陷．针对石油工业中常见的多相混输问题，介绍了常见的气液两相流瞬态模型和流固之间存在的耦合机理，在不作薄壁管假设的前提下推导出了气液两相瞬变流的流固耦合模型．与现有相似模型的对比分析表明，这一模型比较全面地考虑了流体和管道的特性以及不同的耦合形式，可以适应实验和／或仿真研究的需要．     

搅拌槽内气液流动大涡数值模拟研究

为优化搅拌槽体型设计、提高搅拌效率,分别利用气液两相大涡模型和标准k-ε紊流模型,结合多参考系法,对搅拌槽内的气液混合过程进行了模拟计算;控制方程的离散使用了有限体积法;速度与压力耦合求解时使用了压力隐式算子分裂PISO(Pressure-Implicit with Splitting of Operators)算法;模拟自由水面采用了VOF(Volume of Fluid)法;通过计算得到了搅拌槽在相同通气率、不同转轮转速下的气液流动速度场以及不同截面含气率的分布规律。比较两种模型,模拟结果表明:叶片旋转区域紊流的各向异性随着转速的增加而明显加强;基于各向同性假设的标准k-ε紊流模型不能描述流体流速的波动变化,而大涡模拟能够捕捉叶片附近区域流场的分布规律;在通气率一定的情况下,转轮转速的大小对下叶轮附近区域的含气率影响较小,而对上叶轮附近区域含气率的影响较大。     

考虑甲烷溶解效应的液膜密封稳态性能研究

针对油气混输工况密封腔内含气率变化所引起的液膜承载力不稳定问题，考虑密封腔内油气两相介质的互溶性特征，将溶解度方程引入包含Jakobsson–Floberg–Olsson (JFO)边界条件的广义稳态Reynolds方程，建立了考虑甲烷溶解效应的液膜密封润滑模型.采用有限差分法求解该溶解润滑模型，研究了液膜压力、甲烷溶解度及油相黏度之间的相互影响机制.在不同的螺旋槽结构参数与工况条件下，对比分析了甲烷溶解效应对液膜密封机理及密封性能的影响.结果表明：甲烷溶解效应在液膜高压区对油相黏度影响大；考虑甲烷溶解时所得的液膜开启力减小、空化率增大、摩擦系数增大以及泄漏量减小，且液膜动压效应越强时甲烷溶解效应对密封性能的影响越大.在高压及输送油气介质时，气体溶解对密封性能的影响不可忽略.     

基底弹性对蒸发超薄液膜去润湿过程的影响

研究了基底的弹性变形对蒸发超薄膜的稳定性和去润湿动力学过程的影响. 基于长波近 似, 得到了关于液体薄膜厚度的演化方程. 运用线性稳定性理论和数值模拟两种方法, 研究 了基底弹性、范德华力以及液体蒸发等因素对液体薄膜的稳定性和去润湿过程的影响. 研究结果表明增大 基底的弹性系数或者减小液体的表面张力, 都能加速液膜的破碎, 并且能够影响气液界面波 的波长; 液体蒸发能促进气液界面扰动的增长, 有助于液膜的破裂.     

微管内气液两相流中液膜厚度的实验测量

微通道内气液两相流中气柱（plugbubble）与通道壁之间液膜厚度的实验测量，是微热管、微流动、微电子冷却以及气泡雾化等研究中普遍关注的问题。本文利用基于光学干涉和快速傅立叶变换的空间频谱分析方法，实验测量获取了含表面活性剂水中气柱在750μm通道内运动时其与通道壁面之间的液膜厚度。实验结果表明：表面活性剂对液膜厚度的影响比较明显，即当表面活性剂浓度在一定范围内增大时，液膜厚度会减小；此外，当气柱运动速度在一定范围内增大时，液膜厚度也会减小。     

气液两相管流水力瞬变的数值计算及实验研究

由气液两相管流的基本方程出发,通过引入矢通量分裂,对传统的特征线差分做了较大的改进,形成了基于矢通量分裂的特征线差分解法。该法首先将控制方程组的特征值分解成正、负两部分,进而将控制方程中的矢通量雅可比矩阵分裂成两个亚矢量矩阵,对其按各自的迎风格式差分,从而建立了稳定的差分求解格式。该计算法适合于计算声速变化较大且计及液流速度的气液管流的瞬变。计算求解得到的各种不同初始空隙比的压力变化曲线、声速曲线、波速变化曲线、空隙比变化曲线及气体释放影响曲线,通过与不同初始空隙比时气液管流水力瞬变的实验结果对比分析,结果表明两者吻合较好,说明本文方法对于低空隙比的气液两相管流具有较普遍的适用性。     

A two-phase, two-component model for vertical upward gas-liquid annular flow

In this paper, a new two-fluid two-component computational fluid dynamics (CFD) model is developed to simulate vertical upward two-phase annular flow. The two-phase VOF scheme is utilized to model the roll wave flow, and the gas core is described by a two-component phase consisting of liquid droplets and gas phase. The entrainment and deposition processes are taken into account by source terms of the governing equations. Unlike the previous models, the newly developed model includes the effect of liquid roll waves directly determined from the CFD code, which is able to provide more detailed and, the most important, more self-standing information for both the gas core flow and the film flow as well as their interactions. Predicted results are compared with experimental data, and a good agreement is achieved.     

Modeling transient churn-annular flows in a long vertical tube

This work investigates the transient behavior of high gas fraction gas-liquid flows in vertical pipes (annular and churn flows). Hyperbolic balance equations for mass, momentum and entropy are written for the gas and liquid, which is split between a continuous film and droplets entrained in the gas core. Closure relationships to calculate the wall and interfacial friction and the rates of droplet entrainment and deposition were obtained from the literature. A finite-difference solution algorithm based on a coefficient matrix splitting method was implemented to deal with sharp variations in the spatial and temporal domains, such as pressure and phase holdup waves. The model results were compared with steady-state experimental data from eight different sources, totaling more than 1500 data points for pressure gradient, liquid film flow rate and void/core fraction. The absolute average deviation between the model and the data was 17% for the pressure gradient and 5.8% for the void fraction. A comparison of the model results with fully transient air-water data generated in a 49-mm ID, 42-m long vertical pipe is also presented. The experimental results consist of two outlet pressure-induced and two inlet mass flow rate-induced transient tests. Two main transient parameters are compared, namely the local void fraction and the pressure difference between selected points along the test section and the outlet (taken as a reference). The comparisons between the experiments and the numerical model indicate that the model was capable of describing the transient annular to churn flow transition with absolute average deviations of 14.5% and 7.9% for the pressure difference and void fraction, respectively.     

Do huge waves exist in horizontal gas-liquid pipe flow?

Huge waves are periodic interfacial structures which are observed in vertical co-current gas-liquid two-phase flow under churn and the transition between churn and annular flows. Published data examining vertical gas-liquid flow indicate that a huge wave has either a continuous gas core surrounded by a large-scale interfacial wave or a core with a highly-agitated mixture of gas and liquid.Employing a Wire-Mesh Sensor (WMS), the spatio/temporal investigation of high flow rate horizontal air-water flow divulged some recurrent liquid structures (one may call pseudo-slugs) analogous to huge waves of (vertical) churn flow. In both cases, the blow-through (penetration of gas into the liquid structure) was the most manifest feature.Different qualitative and quantitative methods were employed to compare the behavior of pseudo-slug to churn flow. The quantitative measures included Probability Density Function analysis (PDF), distribution coefficient in drift flux model, structural velocity, core average velocity, interfacial friction factor, and slippage number. Both flow regimes demonstrated similar behavior.     

A multiphase,micro-scale PIV measurement technique for liquid film velocity measurements in annular two-phase flow

Prediction methods for two-phase annular flow require accurate knowledge of the velocity profile within the liquid film flowing at its perimeter as the gradients within this film influence to a large extent the overall transport processes within the entire channel. This film, however, is quite thin and variable and traditional velocimetry methods have met with only very limited success in providing velocity data. The present work describes the application of Particle Image Velocimetry (PIV) to the measurement of velocity fields in the annular liquid flow. Because the liquid is constrained to distances on the order of a millimeter or less, the technique employed here borrows strategies from micro-PIV, but micro-PIV studies do not typically encounter the challenges presented by annular flow, including very large velocity gradients, a free surface that varies in position from moment to moment, the presence of droplet impacts and the passage of waves that can be 10 times the average thickness of the base film. This technique combines the seeding and imaging typical to micro-PIV with a unique lighting and image processing approach to deal with the challenges of a continuously varying liquid film thickness and interface. Mean velocity data are presented for air–water in two-phase co-current upward flow in a rectangular duct, which are the first detailed velocity profiles obtained within the liquid film of upward vertical annular flow to the authors’ knowledge. The velocity data presented here do not distinguish between data from waves and data from the base film. The resulting velocity profiles are compared with the classical Law of the Wall turbulent boundary layer model and found to require a decreased turbulent diffusivity for the model to predict well. These results agree with hypotheses previously presented in the literature.     

Prediction of amount of entrained droplets in vertical annular two-phase flow

Prediction of amount of entrained droplets or entrainment fraction in annular two-phase flow is essential for the estimation of dryout condition and analysis of post dryout heat transfer in light water nuclear reactors and steam boilers. In this study, air–water and organic fluid (Freon-113) annular flow entrainment experiments have been carried out in 9.4 and 10.2 mm diameter test sections, respectively. Both the experiments covered three distinct pressure conditions and wide range of liquid and gas flow conditions. The organic fluid experiments simulated high pressure steam–water annular flow conditions. In each experiment, measurements of entrainment fraction, droplet entrainment rate and droplet deposition rate have been performed by using the liquid film extraction method. A simple, explicit and non-dimensional correlation developed by Sawant [Sawant, P.H., Ishii, M., Mori, M., 2008. Droplet entrainment correlation in vertical upward co-current annular two-phase flow. Nucl. Eng. Des. 238 (6), 1342–1352] for the prediction of entrainment fraction is further improved in this study in order to account for the existence of critical gas and liquid flow rates below which no entrainment is possible.Additionally, a new correlation is proposed for the estimation of minimum liquid film flow rate at the maximum entrainment fraction condition. The improved correlation successfully predicted the newly collected air–water and Freon-113 entrainment fraction data. Furthermore, the correlations satisfactorily compared with the air–water, helium–water and air–genklene experimental data measured by Willetts [Willetts, I.P., 1987. Non-aqueous annular two-phase flow. D.Phil. Thesis, University of Oxford]. However, comparison of the correlations with the steam–water data available in literature showed significant discrepancies. It is proposed that these discrepancies might have been caused due to the inadequacy of the liquid film extraction method used to measure the entrainment fraction or due to the change in mechanism of entrainment under high liquid flow conditions.     

Fluctuation characteristics of two-phase flow splitting at a vertical impacting T-junction

In order to investigate the fluctuation characteristics of two-phase flow splitting at a T-junction, particular attention was paid on Churn flow which had the strongest fluctuation comparing with bubble flow and annular flow. The main tube of the T-junction was vertical and the two branches were horizontal. All three pipes connecting to the junction were of 15 mm inner diameter. A statistical analysis based on Root Mean Square (RMS) was applied to temporal differential pressure signals and gas flow rate signals. The Power Spectral Density (PSD) was also employed to reveal their peculiar features in frequency domain as well. The effects of the extraction flow ratio and the gas and liquid superficial velocity upstream on fluctuation characteristics of gas-liquid two-phase flow splitting at the T-junction were investigated in detail. It is found that there is a wide fluctuation in both differential pressure and gas flow rate downstream at every extraction ratio (W₃/W₁) and the fluctuation intensity increases as W₃/W₁ increasing. It is also made clear that increasing either water superficial velocity or gas superficial velocity in inlet causes fluctuation to become more intensive.     

Characteristics of flow and liquid distribution in a gas–liquid vortex separator with multi spiral arms

Fischer–Tropsch (F–T) synthesis is an important route to achieve the clean fuel production. The performance of gas–liquid separation equipment involving in the progressive condensation and separation of light and heavy hydrocarbons in the oil-gas products has become a bottleneck restricting the smooth operation of the F–T process. In order to remove the bottleneck, a gas–liquid vortex separator with simple structure, low pressure drop and big separation capacity was designed to achieve the efficient separation between gas and droplets for a long period. The RSM (Reynolds Stress Model) and DPM (Discrete Phase Method) are employed to simulate the flow characteristics and liquid distribution in the separator. The results show that the separation efficiency is influenced by the flow field and liquid phase concentration in the annular zone. The transverse vortex at the top of spiral arm entrains the droplets with small diameter into the upper annular zone. The entrained droplets rotate upward at an angle of about 37.4°. The screw pitch between neighbor liquid threads is about 0.3 m. There is a top liquid ring in the top of annular zone, where the higher is the liquid phase concentration, the lower is the separation efficiency. It is found that by changing the operating condition and the annular zone height the vortex can be strengthened but not enlarged by the inlet velocity. The screw pitch is not affected by both inlet velocity and annular zone height. The liquid phase concentration in the top liquid ring decreases with both the increases of inlet velocity and annular zone height. The total pressure drop is almost not affected by the annular zone height but is obviously affected by the inlet velocity. When the height of annular zone is more than 940 mm, the separation efficiency is not changed. Therefore, the annular zone height of 940 mm is thought to be the most economical design.     

A model for droplet entrainment in heated annular flow

The ability to accurately predict droplet entrainment in annular two-phase flow is required to effectively calculate the interfacial mass, momentum, and energy transfer, which characterizes nuclear reactor safety, system design, analysis, and performance. Most annular flow entrainment models in the open literature are formulated in terms of dimensionless groups, which do not directly account for interfacial instabilities. However, many researchers agree that there is a clear presence of interfacial instability phenomena having a direct impact on droplet entrainment. The present study proposes a model for droplet entrainment, based on the underlying physics of droplet entrainment from upward co-current annular film flow that is characteristic to light water reactor safety analysis. The model is developed based on a force balance and stability analysis that can be implemented into a transient three-field (continuous liquid, droplet, and vapor) two-phase heat transfer and fluid flow systems analysis computer code.     

Algebraic turbulence modeling in adiabatic gas–liquid annular two-phase flow

The study considers algebraic turbulence modeling in adiabatic gas–liquid annular two-phase flow. After reviewing the existing literature, two new algebraic turbulence models are proposed for both the liquid film and the droplet laden gas core of annular two-phase flow. Both turbulence models are calibrated with experimental data taken from the open literature and their performance critically assessed. Although the proposed turbulence models reproduce the key parameters of annular flow well (average liquid film thickness and pressure gradient) and the predicted velocity profiles for the core flow compare favorably with available core flow velocity measurements, a more accurate experimental database is required to further improve the models accuracy and range of applicability.     

Entrainment in condensing annular flow

Experiments were carried out on low pressure, steam-water, condensing annular flow in a 38.1 mm i.d. horizontal tube. The velocity of the steam at inlet was in the range 97–186 m/s.Measurements of the liquid film flow rate at the end of the test section, which arose as a result of condensation, entrainment, and deposition, were made by extracting the film through a porous sinter bush. The liquid flow rate in the vapour core at exit was deduced from these measurements together with a heat balance on the condenser section.These results were compared with three correlations for entrainment developed from air-water studies. On the basis of the experimental data available, there was sufficient agreement in one case to warrant further investigation.