垂直圆管内液氮流动沸腾的理论模型及数值模拟

分析了液氮流动沸腾过程中气液两相间动量、能量以及质量的传输规律,建立了相应的理论模型,新模型重点修正了界面面积浓度和气泡挣脱直径的计算式;采用新建立的理论模型作为封闭方程对CFX-4.3中内建的双流体模型进行了修正,并采用修正后的双流体模型模拟了液氮在垂直圆管内的流动沸腾过程.数值模拟的结果与文献中的实验数据吻合较好,证明了本文所建模型的合理性.通过数值模拟发现,两相流参数分布的不均匀性对液氮流动沸腾过程中的热质传输特性有重要影响.     

EFFECT OF NATURAL CONVECTION ON MASS TRANSFER BETWEEN SOLID PARTICLES AND FLOWING FLUIDS IN PACKED BEDS

Mass transfer between solid particles in a packed bed and a flowing fluid in laminar flow region, where both forced and natural convection are important, was analyzed theoretically. The diffusion equation based on the free surface model proposed by Happel was solved and an approximate solution relating the average Sherwood number to the Peclet number, the Rayleigh number and the void fraction of a packed bed was obtained. The published experimental mass transfer data in packed beds were found to agree well with the theoretical predictions.     

Improved Void Fraction Correlation for Two‐Phase Flow in Large‐Diameter Annuli

A flow pattern‐independent void fraction correlation for gas‐liquid two‐phase flow in vertical large‐diameter annuli is established. Two equations are proposed for the parameters of a drift flux model‐based correlation: the distribution parameter and the drift flux velocity. These equations are expressed as a function of two‐phase flow variables including void fraction, fluid properties, pipe geometry, and phase flow rates. Experiments were performed to study the void fraction of vertical air‐water two‐phase flow in large‐diameter annuli. The obtained experimental data along with the literature data of Caetano are used to verify the performance of the proposed void fraction correlation. The accuracy of this correlation is compared with nineteen frequently used correlations in literature. The proposed correlation was found to predict the void fraction consistently with a better accuracy.     

基于MUSIG模型的竖直圆管内液氮过冷沸腾数值模拟研究(英文)

Multiple size group （MUSIG） model combined with a threedimensional twofluid model were em ployed to predict subcooled boiling flow of liquid nitrogen in a vertical upward tube. Based on the mechanism of boiling heat transfer, some important bubble model parameters were amended to be applicable to the modeling of liquid nitrogen. The distribution of different discrete bubble classes was demonstrated numerically and the distribu tion patterns of void fraction in the wallheated tube were analyzed. It was found that the average void fraction in creases nonlinearly along the axial direction with wall heat flux and it decreases with inlet mass flow rate and sub cooled temperature. The local void fraction exhibited a Ushape distribution in the radial direction. The partition of the wall heat flux along the tube was obtained. The results showed that heat flux consumed on evaporation is the leading part of surface heat transfer at the rear region of subcooled boiling. The turning point in the pressure drop curve reflects the instability of bubbly flow. Good agreement was achieved on the local heat transfer coefficient aalnst experimental measurements, which demonstrated the accuracy of the numerical model.     

加热上升管内过冷流动沸腾数值模拟

采用计算流体动力学（CFD）程序CFX4.4对加热上升管内过冷流动沸腾工况下气水两相流动局部两相流参数（空泡份额和汽泡尺寸）进行了数值模拟。对数值差分方法、相关模型（界面力和气泡诱导的紊流）和汽泡尺寸进行了敏感性分析。空泡份额分布计算结果与实验结果比较表明,在低空泡份额工况下,两者符合较好,在高空泡份额工况下两者存在一定偏差,并且气相速度和汽泡尺寸的计算结果不理想。计算结果与实验结果之间的差异说明程序模型对于加热上升管内过冷流动沸腾模拟并不完善,建立更为合理的汽泡尺寸模型,考虑汽泡的合并和撕裂是必要的。     

Population balance modelling for bubbly flows with heat and mass transfer

Population balance equations combined with a three-dimensional two-fluid model are employed to predict bubbly flows with the presence of heat and mass transfer processes. Subcooled boiling flow belongs to this specific category of bubbly flows is considered. The MUSIG (MUltiple-SIze-Group) model implemented in CFX4.4 is further developed to account for the wall nucleation and condensation in the subcooled boiling regime. Comparison of model predictions against local measurements near the test channel exit is made for the radial distribution of the bubble Sauter diameter, void fraction, interfacial area concentration and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Additional comparison was also performed against existing boiling model in CFX4.4 and the modified model developed in our previous work (Int. J. Heat Mass Transfer 45 (2002) 1197). Good agreement is better achieved with the local radial bubble Sauter diameter, void fraction, interfacial area concentration and liquid velocity profiles against measurements using the newly formulated MUSIG boiling model over the simpler boiling models. However, significant weakness of the model is still evidenced in the prediction of the vapour velocity. Work is in progress to circumvent the deficiency of the model by the consideration of additional momentum equations or an algebraic slip model to account for bubble separation.     

Void fraction in vertical gas-liquid slug flow: Influence of liquid slug content

In this study we develop a model for computing the mean void fraction and the liquid slug void fraction in vertical upward gas-liquid intermittent flow. A new model for the rate of gas entrained from the Taylor bubble to the liquid slug is formulated. It uses the work done by the pressure force at the rear of the Taylor bubble. Then an iterative approach is employed for equating the gas entrainment flux and the gas flux obtained via conservation equations. Model predictions are compared with experimental data. The developed iterative method is found to provide reasonable quantitative predictions of the entrainment flux and of the void fraction at low and moderate liquid slug void fraction conditions. However, with an increased liquid slug void fraction experimental data indicate that the flow in the liquid slug transits to churn-heterogeneous bubbly flow thus gas entrainment flux tends to zero. Considering this effect in the iterative model significantly improved the predictions for large liquid slug void fraction conditions.     

Effect of the liquid entrainment on the predicted pressure drop in vertical flow-boiling

Flow boiling in a vertical uniformly heated tube has been modelled employing the drift flux and two-fluid formulations. The governing equations are solved numerically and the results presented in the form of pressure drop versus mass flux curves are compared to previously obtained experimental results. The effect of the entrained liquid fraction at the transition from churn-turbulent to annular flow when a switch from the drift flux model to the two-fluid formulation is made is shown to affect the results and could improve the pressure drop predictions.     

Hydrodynamics of an external-loop gas-lift system with restrictions located in the downcomer

Experiments were conducted to study the fluid dynamics in the case that slug flow occurs in the riser of an external-loop gas-lift system with a restriction section located in the downcomer. Complex fluctuation behaviors of the liquid circulation velocity and the wall shear stress in the riser were observed and discussed. Based on the slug flow hydrodynamic behaviors and the balance of momentum and pressure drop over the circulating loop, a model was developed to predict the main parameters of the system: the liquid circulation velocity, the void fraction, the length and velocities of Taylor bubbles and liquid slugs. The predicted results of these parameters were compared with the experimental data and a good agreement was obtained.     

Modeling the pressure drop of wet gas in horizontal pipe

A model for gas-liquid annular and stratified flow through a horizontal pipe is investigated,using the two-phase hydrokinetics theory.Taking into consideration the flow factors including the void fraction,the friction between the two phases and the entrainment in the gas core,the one-dimensional momentum equation for gas has been solved.The differential pressure of the wet gas between the two tapings in the straight pipe has been modeled in the pressure range of 0.1-0.8 MPa.In addition a more objective iteration approach to determine the local void fraction is proposed.Compared with the experimental data,more than 83％ deviation of the test data distributed evenly within the band of ± 10％.Since the model is less dependent on the specific empirical apparatus and data,it forms the foundation for further establishing a flow measurement model of wet gas which will produce fewer biases in results when it is extrapolated.     

摇摆条件下圆管内过冷沸腾局部空泡时空分布特性

过冷流动沸腾在核能、动力、化工等工业领域广泛存在,局部空泡分布特性的准确预测对构建两相流数理模型及两相流动压降和传热特性计算均具有重要意义。实验借助光学探针测量手段研究摇摆条件下圆管内过冷沸腾局部空泡分布特性,并提出多周期同相位叠加方法获得摇摆运动条件下局部瞬时空泡份额。实验结果表明:摇摆条件下圆管内局部空泡分布呈现周期性波动规律,摇摆条件下圆管轴心区及近壁区均会出现局部空泡份额峰值,且流道近壁区局部空泡份额波动幅度最大;与静止条件下不同,摇摆运动周期内过冷沸腾局部空泡份额空间分布形态随时间发生变化;摇摆周期和摇摆角度对过冷沸腾局部空泡分布特性也存在一定影响。     

Two-fluid modelling of the rewetting and refilling of hot horizontal tubes

This paper deals with the use of a two-fluid model for predicting the thermal-hydraulic characteristics of the rewetting and refilling of hot horizontal tubes. The two-fluid model equations were simplified and solved using an explicit finite difference scheme. A set of constitutive equations is used for closure. A new film boiling model and a model for evaluating the rate of latent to total heat in film boiling were incorporated in the two-fluid model. A newly proposed rewetting criterion, based on vapour film collapse, was used to mark the end of the film boiling regime. The model is capable of predicting the rewetting and refilling process in horizontal tubes under conditions of stratified or inverted annular refilling flow. The model predictions were found to be in good agreement with the experimental results.     

γ射线法测量高压管束间气液两相流的截面含气率分布

在立式蒸汽发生器垂直管束间的气液两相流中，截面含气率是其中一个重要参数。使用γ射线法对高温高压下垂直管束间气液两相流截面含气率的分布规律进行了实验研究。实验压力分别为5、7、9 MPa，质量流速为300 kg/(m²?s)，热力学干度的范围为0.003 ~ 0.4。实验得到了垂直管束间截面含气率随热力学干度、体积含气率和压力的变化关系；并与经典公式的计算结果对比发现，在低干度区域，实验结果与Miropolskii模型、Smith模型和Armand模型偏差较大，均大于30%，在高干度区域偏差较小；基于Armand理论，通过多元线性回归法拟合出本文工况下平均截面含气率的计算关联式，与日本核动力工程公司(NUPEC)的实验数据偏差小于15%。本研究对蒸汽发生器的结构设计和流动特性研究具有重要意义。     

竖直环形通道内液氮流动沸腾的数值模拟

在多尺寸组模型的基础上,从加热壁面上脱离汽泡的受力分析人手,对液氮过冷流动沸腾模型进行了修正.将新模型应用于环形通道内液氮过冷流动沸腾的数值模拟,同时为了比较,采用基于Kirichenko,Fritz汽泡脱离直径公式的多尺寸组模型对同一管道内液氮过冷流动进行了数值模拟.结果表明:结合脱离汽泡受力分析模型的多尺寸组模型可...     

Continuous production of solid polystyrene in back‐mixed and linear‐flow reactors

A mathematical model has been developed to predict the steady state performance of a continuous bulk styrene polymerization process with catalytic initiation for solid polystyrene. The polymerization section contains one boiling CSTR, followed by multiple linear‐flow reactors. The devolatilization section consists of two polymer pre‐heaters and two high‐solids flashes. The polymer moment equations were solved simultaneously with the reactor modeling equations. The non‐linear algebraic equations were solved by a Newton‐Raphson iteration technique to give the steady‐state styrene monomer weight fraction in a CSTR. The coupled, non‐linear ordinary differential equations were numerically integrated using a single‐step, 4th‐order Runge‐Kutta technique, followed by a multi‐step Adams‐Moulton technique. The resulting computer simulation model is capable of evaluating how the production rate and product quality are affected by feed composition, temperature, initiator type, initiator concentration, and residence time. Several case studies were given for commercially important crystal‐clear and impact‐resistant resins. A binary initiation system gives a good balance of monomer conversion, polymer molecular weights, and rubber grafting compared to a single initiation system. The styrene dimer/trimer occur in low concentrations but can be substantially reduced with a low temperature initiator. The ideal mean residence time is approximately one minute or less in a shell‐and‐tube devolatilization pre‐heater. Low flash chamber vacuum is more effective than high polystyrene melt temperature to reduce the volatile content of the final product. The water injected to the low volatile melt shows promising improvement in the second‐stage polystyrene devolatilization.     

矩形通道干涸后最小膜沸腾传热试验研究

在中国核动力研究设计院流动传热基础试验平台上进行了矩形通道干涸后最小膜沸腾传热试验。通过对各种热工水力参数的试验研究,得出结论:随着进口含汽质量分数的增加,最小膜沸腾热流密度减小,含汽质量分数增加,壁面温度降低,传热系数减小;随着质量流速的增大,最小膜沸腾热流密度增大,含汽质量分数减小,壁面温度升高,传热系数增大;随着系统压力的升高,最小膜沸腾热流密度增大,含汽质量分数增加,壁面温度升高,传热系数增大;相对于偏离泡核沸腾(DNB)后最小膜沸腾,干涸后最小膜沸腾的热流密度更大,壁面过热度较低,传热系数更大。     

蒸发器沸腾两相段熵产分析

以管内侧制冷剂,管外侧水的蒸发换热器为研究对象,基于热力学第一、第二定律建立了蒸发器沸腾两相流段的熵产计算模型,用量纲1熵产数N_s表示管内外流体因温差传热和压降损失引起的不可逆损失。分析水侧、制冷剂侧质量流量,水侧入口温度和管长尺寸对系统熵产数的影响。结果表明：系统熵产主要由制冷剂侧两相压降引起;系统熵产数随水侧质量流率的增大迅速降低并趋于稳定,随制冷剂侧质量流率的增大呈指数增长,随水侧进口温度和管长尺寸的增大而增加。     

基于ECT和蚂蚁算法的油气两相流空隙率在线测量

基于电容层析成像系统（ECT）和蚂蚁算法,提出了一种油气两相流空隙率在线测量的新方法。该方法利用电容层析成像系统12电极电容传感器所获取的66个测量电容信息,首先根据电容层析成像系统所获取的流型辨识结果确定对应流型下的实际空隙率测量模型参数f和b,然后利用蚂蚁算法的信息素信息,找到当前测量状态下对空隙率起主要作用的组合电容集合和相应的权重系数,从而实现空隙率测量。与流型相对应的空隙率测量线性模型参数f和b基于先验数据通过最小二乘方法确定。油气两相流的实验结果表明,该方法对空隙率的在线测量是有效的,避免了复杂的图像重建计算,实时性能佳,测量时间小于0.08 s。在几种典型流型下,提出的空隙率测量方法与常用的快关阀方法相比最大测量偏差小于5.5%。     

Estimation of void boundaries in flow field using expectation–maximization algorithm

Phase distribution in the flow field provides an insight into the hydrodynamics and heat transfer between the fluids. Void fraction, which is one of the key flow parameters, can be determined by estimating the phase boundaries. Electrical impedance tomography (EIT), which has high temporal characteristics, has been used as an imaging modality to estimate the void boundaries, using the prior knowledge of conductivities. The voids formed within the process vessel are not stable and their movement is random in nature, thus dynamic estimation schemes are necessary to track the fast changes. Kalman-type estimators like extended Kalman filter (EKF) assume the knowledge of model parameters, such as the initial states, state transition matrix and the covariance of process and measurement noise. In real situations, we do not have the prior information of the model parameters; therefore, in such circumstances the estimation performance of the Kalman-type filters is affected. In this paper, the expectation–maximization (EM) algorithm is used as an inverse algorithm to estimate the model parameters as well as non-stationary void boundary. The uncertainties caused in Kalman-type filters, due to the inaccurate selection of model parameters are overcome using an EM algorithm. The performance of the method is tested with numerical and experimental data. The results show that an EM has better estimation of the void boundary as compared to the conventional EKF.     

基于截面气含率的文丘里湿气压降模型

运用两相流理论对湿气中的气相与液相流动进行分析,在分层流与环雾状流的条件下推导了两相流通过水平标准文丘里流量计的理论模型.通过分别考虑截面气含率、相间摩擦以及液滴夹带等因素,在文丘里轴向对湿气流动中的气相动量方程进行求解.通过对水平直管中的截面气含率公式进行修正,建立了适用于收缩管道的截面气含率模型,并在此基础上模拟了湿气流经标准文丘里时其两个取压孔之间的轴向静压分布.实验证实,使用修正后的截面气含率公式将使模型对文丘里压降的预测准确度明显提升,其相对误差在15%以内.该模型以湿气两相流在水平文丘里中的流动形态为依据,具有较充分的物理背景,而且在推导过程中较少依赖特定实验装置与数据,为建立具有一定普适性的文丘里湿气计量模型奠定了基础.