Pure steam condensation model with laminar film in a vertical tube

A new physical model for calculating the liquid film thickness and condensation heat transfer coefficient in a vertical condenser tube is proposed by considering the effects of gravity, liquid viscosity, and vapor flow in the core region of the flow. To estimate the velocity profile in the liquid film, the liquid film was assumed to be in Couette flow forced by the interfacial velocity at the liquid–vapor interface. For simplifying the calculation procedures, the interfacial velocity was estimated by introducing an empirical power-law velocity profile. The resulting film thickness and heat transfer coefficient from the model were compared with the experimental data and the results obtained from the other condensation models. The results demonstrated that the proposed model described the liquid film thinning effect by the vapor shear flow and predicted the condensation heat transfer coefficient from experiments reasonably well.     

A separated flow model for predicting two-phase pressure drop and evaporative heat transfer for vertical annular flow

A separated flow model has been developed that is applicable to vertical annular two-phase flow in the purely convective heat transfer regime. Conservation of mass, momentum, and energy are used to solve for the liquid film thickness, pressure drop, and heat transfer coefficient. Closure relationships are specified for the interfacial friction factor, liquid film eddy-viscosity, turbulent Prandtl number, and entrainment rate. Although separated flow models have been reported previously, their use has been limited, because they were tested over a limited range of flow and thermal conditions. The unique feature of this model is that it has been tested and calibrated against a vast array of two-phase pressure drop and heat transfer data, which include upflow, downflow, and microgravity flow conditions. The agreements between the measured and predicted pressure drops and heat transfer coefficients are, on average, better or comparable to the most reliable empirical correlations. This separated flow model is demonstrated to be a reliable and practical predictive tool for computing two-phase pressure drop and heat transfer rates. All of the datasets have been obtained from the open literature.     

Comparison of film condensation models in presence of non-condensable gases implemented in a CFD Code

Several film condensation models in presence of non-condensable gases are presented. They have been implemented in a CFD code and compared with experimental data. The aim was to improve the code for simulating the gas mixing process in large containment buildings involving steam. The models based on correlation are more robust and simpler, but they work badly out of their experimental conditions. The mechanistic models, based on the diffusion layer theory, work well in numerous conditions but the algorithm are more complicated. Moreover, they run badly when the convective heat transfer is not well predicted by the code.The main material of this paper is part of the PhD Thesis Numerical Models Implemented in a Computational Fluid-dynamics Code for Hydrogen Combustion Risk Assessment within Containments, J.M. Martín-Valdepeñas, Universidad Politécnica de Madrid, 2004.     

Free- and forced-convection film condensation from a horizontal elliptic tube with a vertical plate and horizontal tube as special cases

A simple mathematical model is developed for the study of the mixed-convection film condensation with downward flowing vapors onto a horizontal elliptic tube. Analytical analysis for both the local condensate film thickness and heat transfer characteristics under simultaneous effects of interfacial vapor shear and pressure gradient has been performed by adopting a unified geometry parameter, eccentricity e. The present results for two limit cases, e = 0 (circular tube) and e = 1.0 (vertical plate) are in an excellent agreement with the earlier works. For very slow vapor flow, the present result for dimensionless mean heat transfer coefficient reduces to the same form as in the earlier works, , whose value is 0.728 for e = 0 and 0.943 for e = 1.0. As for very fast vapor flow, the dimensionless mean heat transfer coefficient, increase with increasing eccentricity under the effects of pressure gradient caused by potential flow and surface tension.     

Experimental study on condensation heat transfer of steam on vertical titanium plates with different surface energies

Visual experiments were employed to investigate heat transfer characteristics of steam on vertical titanium plates with/without surface modifications for different surface energies. Stable dropwise condensation and filmwise condensation were achieved on two surface modification titanium plates, respectively. Dropwise and rivulet filmwise co-existing condensation form of steam was observed on unmodified titanium surfaces. With increase in the surface subcooling, the ratio of area (η) covered by drops decreased and departure diameter of droplets increased, resulting in a decrease in condensation heat transfer coefficient. Condensation heat transfer coefficient decreased sharply with the values of η decreasing when the fraction of the surface area covered by drops was greater than that covered by rivulets. Otherwise, the value of η had little effect on the heat transfer performance. Based on the experimental phenomena observed, the heat flux through the surface was proposed to express as the sum of the heat flux through the dropwise region and rivulet filmwise region. The heat flux through the whole surface was the weighted mean value of the two regions mentioned above. The model presented explains the gradual change of heat transfer coefficient for transition condensation with the ratio of area covered by drops. The simulation results agreed well with the present experimental data when the subcooling temperature is lower than 10 °C.     

Forced convective flow and heat transfer of upward cocurrent air–water slug flow in vertical plain and swirl tubes

This experimental study comparatively examined the two-phase flow structures, pressured drops and heat transfer performances for the cocurrent air–water slug flows in the vertical tubes with and without the spiky twisted tape insert. The two-phase flow structures in the plain and swirl tubes were imaged using the computerized high frame-rate videography with the Taylor bubble velocity measured. Superficial liquid Reynolds number (Re_L) and air-to-water mass flow ratio (AW), which were respectively in the ranges of 4000–10000 and 0.003–0.02 were selected as the controlling parameters to specify the flow condition and derive the heat transfer correlations. Tube-wise averaged void fraction and Taylor bubble velocity were well correlated by the modified drift flux models for both plain and swirl tubes at the slug flow condition. A set of selected data obtained from the plain and swirl tubes was comparatively examined to highlight the impacts of the spiky twisted tape on the air–water interfacial structure and the pressure drop and heat transfer performances. Empirical heat transfer correlations that permitted the evaluation of individual and interdependent Re_L and AW impacts on heat transfer in the developed flow regions of the plain and swirl tubes at the slug flow condition were derived.     

Experimental study of inverted annular film boiling in a vertical tube cooled by R-134a

An experimental investigation of inverted annular film boiling heat transfer has been performed for vertical up-flow in a round tube. The experiments used R-134a coolant and covered a pressure range of 640–2390 kPa (water equivalent range: 4000–14,000 kPa) and a mass flux range of 500–4000 kg m⁻² s⁻¹ (water equivalent range: 700–5700 kg m⁻² s⁻¹). The inlet qualities of the tests ranged from −0.75 to −0.03. The hot-patch technique was used to obtain the subcooled film boiling measurements. It was found that the heat transfer vs. quality curve can be divided into four different regions, each characterized by a different mechanisms and trends. These regions are dependent on pressure, mass flux and local quality. A detailed examination of the parametric trends of the heat transfer coefficient with respect to mass flux, inlet quality, heat flux and pressure was performed; reasonably good agreement between observed trends and those reported in the literature were noted.     

An experimental technique for measuringtransient natural/forced convective heatfluxes in a vertical channel

A series of systematic experiments for measuring transient natural andforced convected heat fluxes in a one-sided heated vertical channel have beenconducted. The total heat input in transient natural and forced experiments iscomposed of four kinds modes: radiative heat loss, conductiveheat lossm, thermal capacity of the plate, and convective heat transfer in thetest channel. In the transient periods, the generalized correlations for estimatingtransient convective heat flux are proposed in natural and forced convectionexperiments. By using the proposed q_c^″ distributions, satisfactory agreementsare achieved between the transient Nu_x predictions and experimentaldata.     

Turbulent forced convective heat transfer in two-phase evaporating droplet flow through a vertical pipe

A physical model was developed to study heat transfer in turbulent dispersed flow at very high vapor quality in a vertical pipe by numerically solving the coupling governing differential equations for both phases. Major heat transfer mechanisms included in the model were the thermal nonequilibrium effects, droplet vaporization, droplet deposition on the duct wall and thermal radiative transfer. The predicted results indicated that vapor superheating is dominant for the cases with high wall superheat, otherwise droplet vaporization dominates the energy transport processes. Heat transfer during the droplet-wall interaction only exists at low wall superheat but in small amounts.     

MHD effects on free convective flow over moving semi-infinite vertical cylinder with temperature oscillation

Numerical solutions of magnetodynamics(MHD) effects on the free convective flow of an incompressible viscous fluid past a moving semi-infinite vertical cylinder with temperature oscillation are presented.The dimensionless,unsteady,non-linear,and coupled governing partial differential equations are solved by using an implicit finite difference method of the Crank-Nicolson type.The velocity,temperature,and concentration profiles are studied for various parameters.The local skin-friction,the average skin-fr...     

Flow condensation heat transfer characteristics of hydrocarbon refrigerants and dimethyl ether inside a horizontal plain tube

Flow condensation heat transfer coefficients (HTCs) and pressure drop of R22, propylene, propane, DME and isobutane are measured on a horizontal plain tube. The main test section in the experimental flow loop is made of a plain copper tube of 8.8 mm inner diameter and 530 mm length. The refrigerant is cooled by passing cold water through the annulus surrounding the test section. Tests are performed at a fixed refrigerant saturation temperature of 40 ± 0.2 °C with mass fluxes of 100, 200, and 300 kg/m² s and heat flux of 7.3–7.7 kW/m². The heat transfer and pressure drop data are obtained in the vapor quality range of 10–90%. Test results show that for a given mass flux the flow condensation HTCs of propylene, propane, DME and isobutane are higher than those of R22 by up to 46.8%, 53.3%, 93.5% and 61.6%, respectively. Also well-known correlations developed based upon conventional fluorocarbon refrigerants predict the present data within a mean deviation of 33%. Finally, the pressure drop increases as the mass flux and quality increase and isobutane shows the highest pressure drop due to its lowest vapor pressure among the fluids tested.     

Characteristics of film condensation of supersaturated steam–air mixture on a flat plate

Experimental and numerical work was performed for the laminar film condensation of steam–air mixture flow over a flat plate. For small temperature difference between the gas mixture and the cold wall, the gas mixture in the boundary layer can be treated as superheated gas. When the temperature difference is large, the gas mixture becomes supersaturated near the interface. In that case, mist formed near the interface, the temperature profile of the gas mixture was greatly concaved toward the interface and the heat transfer was enhanced. However the velocity profile measured by the laser Doppler anemometer (LDA) showed the same trend without mist formation. A calculation model is proposed and compared with the experimental data and previous models for the superheated or the saturated conditions.     

滑滚表面间的流体膜壁面滑移和流体动力学

经典雷诺润滑理论建立在无壁面滑移的假设基础之上。近年来许多试验报告了发生在流体膜流动的壁面滑移证据。本文研究了两固体表面间的流体膜流动特性和流体动力学,发现壁面滑移显著影响膜的流体动力学问题,流体动压力不仅受黏度和几何间隙的影响,而且还由壁面滑移和表面运动强力控制,通过控制表面的吸附性质,甚至可以得到零摩擦表面。另一方面,如果两个表面具有相同的滑移特性,存在一个临界滑动速度使得流体动压效应完全消失;但是在纯滚动条件下,即使界面极限剪应力很小,仍然有相当可观的流体动压效应。     

Experimental investigation of heat transfer and pressure drop characteristics of flow through spirally fluted tubes

Spirally fluted tubes are used extensively in the design of tubular heat exchangers. In previous investigations, results for tubes with flute depths e/D_vi < 0.2 were reported, with most correlations applicable for Re ≥ 5000. This paper presents the results of an experimental investigation of the heat transfer and pressure drop characteristics of spirally fluted tubes with the following tube and flow parameter ranges: flute depth e/D_vi = 0.1−0.4, flute pitch p/D_vi = 0.4−7.3, helix angle θ/90° = 0.3−0.65, Re = 500−80,000, and Pr = 2−7. The heat transfer coefficients inside the fluted tube were obtained from measured values of the overall heat transfer coefficient using a nonlinear regression scheme. The friction factor data obtained consisted of 507 data points. The proposed correlation for the friction factor predicts 96% of the database within ±20%. The heat transfer correlation for the range 500 ≤ Re ≤ 5000 predicts 76% of the database (178 data points) within ±20%, and the correlation for the higher Re range predicts 97% of the 342 data points within ±20%. Comparison of heat transfer and friction data show that these tubes are most effective in the laminar and transition flow regimes. The present results show that the increase of flute depth in the range considered does not improve heat transfer.     

刚性圆管中血液周期振荡流的切应力分布

本文通过求解圆管内血液振荡流的基本方程，求得圆管内血液流的压力梯度与切应力之间的关系式。在此基础上，详细讲座了圆管中轴向流速和切变率谐波的变化规律，指出流速谐波和切变率谐波的幅值都将随着谐波次数的增大而逐渐减小。为了使所得结果便于应用。文章通过管轴向中心线流速与压力梯度之间的关系式，进一步给出一种利用管轴向中心线流速计算管内切应力分布的简便方法。该方法用于检测活体血管内血液振荡流的切应力分布，具有操作简单，精度较高的优点。最后，以人体颈动脉为例，讨论血液周期振荡流的切应力的分布特性。发现在任意时刻，除了邻近管壁处切应力急剧增大到一定数值之外，沿管截面切应力分布相当均匀且接近于零，呈现出与定常流不同的切应力分布特征。     

Film thickness measurement with an ultrasonic transducer

A technique for measuring condensate film thickness using an ultrasonic transducer is described. In the experiment, the condensate film thickness with R-113 and FC-72 (a fluorinert compound developed by the 3M Company) condensing on the horizontal lower surface of a rectangular duct was measured at several locations. From the measured values a power law relation between the condensate film thickness and the axial distance from the leading edge of the condensing surface was derived by regression analysis. Assuming a linear temperature profile in the condensate film, local and average heat transfer coefficients were computed from the condensate film thickness. The average heat transfer coefficients were compared with the values obtained by measuring the heat transfer rate to the coolant. The two values were within ±12% of each other. As yet there is no satisfactory analytical model to predict the local heat transfer coefficient even in the annular condensation regime. One of the main difficulties in modeling the condensation is the lack of a suitable model to predict the interfacial shear stress. With the measurement of the film thickness it is possible to determine the interfacial shear stress. It is hoped that the shear stresses so determined will lead to the development of a satisfactory model for interfacial shear stress with condensation.     

全长粘结化学植筋的载荷传递机理

运用弹性力学空间轴对称问题与剪切位移法的基本假定和解答方法, 建立分析模型和界面应力基本方程, 得到化学植筋的界面应力和植筋轴力的弹性解析解答, 结合算例说明了界面应力的衰减曲线分布规律, 提出载荷传递系数$k$并简要分析了它对载荷传递规律影响.     

Forced convection from an isolated heat source in a channel with porous medium

A numerical study is made of flow and heat transfer characteristics of forced convection in a channel that is partially filled with a porous medium. The flow geometry models convective cooling process in a printed circuit board system with a porous insert.The channel walls are assumed to be adiabatic. Comprehensive numerical solutions are acquired to the governing Navier-Stokes equations, using the Brinkman-Forchheimer-extended Darcy model for the regions of porous media. Details of flow and thermal fields are examined over ranges of the principal parameters; i.e., the Reynolds number Re, the Darcy number Da (≡K/H²), the thickness of the porous substrate S, and the ratio of thermal conductivities R_k (≡k_eff/k). Two types of the location of the porous block are considered. The maximum temperature at the heat source and the associated pressure drop are presented for varying Re, Da, S, and R_k. The results illustrate that as S increases or Da decreases, the fluid flow rate increases. Also, as R_k increases for fixed Da, heat transfer rates are augmented. Explicit influences of Re on the flow and heat transport characteristics are also scrutinized. Assessment is made of the utility of using a porous insert by cross comparing the gain in heat transport against the increase in pressure drop.     

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明,模型方程在小到中等Reynolds数下都适用,并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加,展向不稳定性增长率则随着Marangoni数增加而增加,随着Reynolds数增加而减小,流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明,在等热流密度加热条件下,液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构,液膜与气体的换热也因溪流状结构的出现而加强;在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强,液膜会提前破裂;在随机扰动的基础上引入占优势地位的流向最不稳定扰动时,液膜的换热会增强,但不会提前破裂;在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时,换热会加强且液膜会提前破裂。