单液滴撞击球面液膜水花形成特征分析

应用VOF(volume of fluid)方法对液滴撞击球面液膜的飞溅过程进行数值模拟,由计算结果分析了不同液滴撞击速度和液膜厚度对撞击产生的水花形态、高度和直径的影响。结果表明:随着液滴撞击速度和液膜厚度的增加,水花高度增大;随着液滴撞击速度增加,水花直径增大,但当液膜厚度增加时,水花直径减小。     

单液滴撞击平面液膜飞溅过程的CLSVOF模拟

采用CLSVOF(coupled level set and volume of fluid)方法对单液滴撞击平面液膜初期的飞溅过程进行了数值模拟,探讨了撞击速度和液膜厚度对撞击后形成冠状水花形态及扩展直径的影响.结果显示,扩展直径随撞击速度的增大而增大,随液膜厚度的增大而减小,CLSVOF方法可较好地用于液滴撞击过程的数值计算.     

Evolution and heat transfer after droplet impact on heated liquid film with vapor bubbles inside

Abstract

To better understand the droplet impact on the liquid film with vapor bubbles in spray cooling, a two-dimensional symmetric numerical model is set up using the Coupled Level Set and Volume of Fluid method (CLSVOF). Three simulative cases are taken, considering the effects of film thickness and the presence of vapor bubbles or not. The main purposes of this paper are to investigate the evolution of vapor bubbles during droplet impact and to identify the effect of vapor bubbles on convection heat transfer. The results indicate that vapor bubbles will detach from the wall and break up at the surface of the liquid film during droplet impact, for a thinner film, later a “sawtooth” liquid film appears at the non-impact region. However, for a thicker film, no bubbles rupture and the detached bubbles will flow inside the liquid film and then some will merge into larger bubbles. In the presence of vapor bubbles, the crater radius is larger for a thicker liquid film. The presence of vapor bubbles will facilitate the subcooled droplet to spread to the heated wall, leading to a substantial increase in surface heat flux.     

Flow and heat transfer around a single row of circular cylinders

A numerical investigation of flow and heat transfer around a single row of circular cylinders was conducted using a boundary-fitted coordinate system. Numerical calculations for center-to-center distance between cylinders L/d=2.0, 2.5, 3.3, and ∞ were made of water flows in the Reynolds number range from 75 to 500. Numerical values of average Nusselt number for uniformly heated cylinders are in relatively good agreement with those obtained from experiments in water (Prandtl number Pr ≒ 8). The interaction of wake flows behind cylinders, observed in the experiments, was also found to occur with decreasing cylinder spacing L/d. © 1997 Scripta Technica, Inc. 25(3): 192–200, 1996     

A coupled level set and volume-of-fluid simulation for heat transfer of the double droplet impact on a spherical liquid film

A coupled level set and volume-of-fluid method is applied to investigate the double droplet impact on a spherical liquid film. The method focuses on the analysis of surface curvature, droplet diameter, impact velocity, double droplets vertical spacing, the thickness of the liquid film of two liquid droplets after the impact on a spherical liquid film, and the influence of flow and heat transfer characteristics. The results indicate that the average wall heat flux density of the double liquid droplet impact on a spherical liquid film is greater than that of a flat liquid film. Average wall heat transfer coefficient increases with the increase in the liquid film’s spherical curvature. When the liquid film thickness is smaller, the average wall heat flux density of the liquid film is significantly reduced by the secondary droplets generated from the liquid film. When the liquid film thickness is larger, the influence of liquid film thickness on the average wall heat flux density gradually decreases. The average wall heat flux density increases with the increase in impact velocity and the droplet diameter; it also decreases with the increase in double droplets vertical spacing.     

强化凝结与液膜传热过程的板壳式换热器波纹方案

板壳式换热器敝开的板外侧通道可以实现多种板式换热器所不能实现的传热过程，特别是液膜传热伟质过程。针对不同过程的特点设计波纹的波型是板壳式换热器推广应用的关键，文章介绍了适合于凝结过程的同向双尺度波纹板片和适合于吸收，蒸发等液膜传热传质过程的交叉双尺度波纹板片方案。     

Characteristics of transient heat transfer during quenching of a vertical hot surface with a falling liquid film

An experimental study has been conducted to elucidate characteristics of transient heat transfer during quenching of a vertical hot surface with a falling liquid film. The experiment was done at atmospheric pressure for the following conditions: an initial surface temperature from 200 to 400^°C, a subcooling of 20– 80 K, average velocity of 0.52– 1.24 m/s, and the block material is copper and carbon steel. The surface temperature and heat flux are estimated from the measured temperatures in the block during the quench by a two‐dimensional inverse solution. It follows that as the position of wetting advances downward, the position at which the heat flux becomes a maximum also advances downward. The time at which the position of maximum heat flux begins to move is one of the most important parameters and can be predicted by a proposed correlation. In addition, it is revealed that the maximum heat flux for copper depends on the length to which it occurs from the leading edge. © 2007 Wiley Periodicals, Inc. Heat Trans Asian Res, 36(6): 345– 360, 2007; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.20167     

Numerical analysis on heat transfer enhancement by waves on falling liquid film

IntroductionLiquid films flowing on a vertical or inclined wall bythe gravitational force are encountered in the wideindustrial and engineering fields['1, such as condensatefilm, evaporating falling film, gas absorb by liquid film,etc. In the case of Indnar film flow with smooth surface,heat transfer through the liquid film is mainly carried outby conduction, and it is sufficiently explained by theNusselt's theory. On the other hand, the heat transfer isfairly enhanced for films generating su…     

Flow and heat transfer of liquid plug and neighboring vapor slugs in a pulsating heat pipe

A model of fluid flow and heat transfer on liquid slug and neighboring vapor plugs in a pulsating heat pipe (PHP) is proposed. A new energy equation for the liquid slug is built by aid of Lagrange method. The shear stress term related with the fluid flow state is included in the motion equation of the liquid slug. A sensitive heat term is replaced by a phase change term in the energy equation of the vapor plug. Based on our analysis on the displacement variation of the liquid slug with time, it is known that the harmonic force acting on the liquid slug in PHPs is the pressure difference between the vapor plugs. The flow oscillation can be considered as a forced damping vibration of one degree of freedom system. The phase difference of the oscillating flow between with and without the gravity effect can reach 45°. The amplitude and angular frequency of flow oscillation is irrespective with the initial displacement of liquid slug. If the flow pattern remains strictly slug flow in the entire system, the contribution of the sensible heat exchange to the total heat transfer of the PHP is about 80%.     

Flow regime-based modeling of heat transfer and pressure drop in microchannel flow boiling

Local heat transfer coefficients and pressure drops during boiling of the dielectric liquid fluorinert FC-77 in parallel microchannels were experimentally investigated in recent work by the authors. Detailed visualizations of the corresponding two-phase flow regimes were performed as a function of a wide range of operational and geometric parameters. A new transition criterion was developed for the delineation of a regime where microscale effects become important to the boiling process and a conventional, macroscale treatment becomes inadequate. A comprehensive flow regime map was developed for a wide range of channel dimensions and experimental conditions, and consisted of four distinct regions – bubbly, slug, confined annular, and alternating churn/annular/wispy-annular flow regimes. In the present work, physics-based analyses of local heat transfer in each of the four regimes of the comprehensive map are formulated. Flow regime-based models for prediction of heat transfer coefficient in slug flow and annular/wispy-annular flow are developed and compared to the experimental data. Also, a regime-based prediction of pressure drop in microchannels is presented by computing the pressure drop during each flow regime that occurs along the microchannel length. The results of this study reveal the promise of flow regime-based modeling efforts for predicting heat transfer and pressure drop in microchannel boiling.     

Mass transfer enhancement through Marangoni instabilities during single drop formation

The impact of Marangoni convection on the extraction efficiency during the drop formation stage is investigated in the system toluene/acetone/water for different initial solute concentrations and different drop diameters. Both mass transfer directions of the solute have been considered. Marangoni instabilities are supposed to increase the internal mixing and thus enhance mass transfer coefficients. Experimental results show a strong dependency on the mass transfer direction. The amount of solute extracted is between 19% and 55%. The total transferred mass M_A increases with drop diameter and initial concentration. Present models from the literature which predict extraction efficiencies do not take into account interfacial effects like Marangoni convection. A correlation is proposed introducing an effective diffusivity which depends on the initial solute concentration. The diffusivity factor increases linearly with initial solute concentration and is more sensitive in the mass transfer direction c → d.     

Mass transfer between a single drop and a continuous phase

The time dependent convective-diffusion equations for mass transfer between a drop and a continuous phase are solved in two cases: (1) the case of small Reynolds numbers and (2) the case of potential flow. The equations are solved by means of a similarity variable η_i = y/δ_i(θ, t) which enable their transformation into an ordinary differential equation for the concentration c_i = c_i(η_i) and a first order equation with partial derivatives for δ_i = δ_i(θ, t). Equations for the mass-transfer coefficient for the unsteady and steady states are obtained. The time in which the steady state is reached is evaluated.     

Numerical Study of Heat Transfer Enhancement by Liquid Film on the Walls

ＮｕｍｅｒｉｃａｌＳｔｕｄｙｏｆＨｅａｔＴｒａｎｓｆｅｒＥｎｈａｎｃｅｍｅｎｔｂｙＬｉｑｕｉｄＦｉｌｍｏｎｔｈｅＷａｌｌｓ￥ＣｈｕｎｌｉｎＸｉａ（ＤｅｐａｒｔｍｅｎｔｏｆＰｏｗｅｒＥｎｇｉｎｅｅｒｉｎｇ，ＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏ...     

Laminar flow and heat transfer of non-newtonian falling liquid film on a horizontal tube with variable surface heat flux

An analysis is performed to study the laminar flow and heat transfer of non-Newtonian falling liquid film on a horizontal tube for the case of variable surface heat flux. The inertia and convection terms are taken into account. The governing boundary layer equations are solved numerically using an implicit finite difference method. Of particular interest are the effects of the mass flow rate Γ, the concentration C of carboxymethylcellulose (CMC) solutions, the exponent m for the power-law surface heat flux, and the tube diameter D on the film thickness profiles, as well as on the local and average Nusselt numbers. It was found that an increase in the mass flow rate Γ and exponent value m increases the local and average heat transfer rates. Finally, the present simulation is found to be in good agreement with previous experimental and numerical results for Newtonian films.     

液滴撞击液膜过程气泡卷吸的FTM数值模拟

采用界面追踪法(front tracking method,FTM)对液滴撞击液膜形成的气泡卷吸现象进行直接数值模拟,分析不同无量纲液膜厚度(H*)和不同Bo对卷吸气泡出现时间和消失时间的影响。结果表明:在Bo=4.00时,H*的增加会使得卷吸气泡出现和消失的时间提前,并且在H*0.415时,卷吸气泡不再出现;在H*分别为0.250和0.150时,卷吸气泡不再出现的临界Bo分别为3.58和3.19,但不同的H*对卷吸气泡出现和消失的时间的影响与固定Bo时趋势一致。     

Research on a new heat transfer model for liquid film convective evaporation in annular flow

Based on the phenomenon of turbulence restraint in liquid‐vapor interface, an analytical model is proposed for annular flow with a velocity distribution. The liquid‐vapor interface affecting district mixing length model was amended, and a new liquid film convective evaporation heat transfer model at the annular flow was developed. Compared with the experimental data, the results show that the new model is better than the model based on full tube flow velocity distribution. © 2002 Wiley Periodicals, Inc. Heat Trans Asian Res, 31(7): 524–530, 2002; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10051     

The effect of liquid film evaporation on flow boiling heat transfer in a micro tube

Flow boiling in micro channels is attracting large attention since it leads to large heat transfer area per unit volume. Generated vapor bubbles in micro channels are elongated due to the restriction of channel wall, and thus slug flow becomes one of the main flow regimes. In slug flow, sequential bubbles are confined by the liquid slugs, and thin liquid film is formed between tube wall and bubble. Liquid film evaporation is one of the main heat transfer mechanisms in micro channels and liquid film thickness is a very important parameter which determines heat transfer coefficient. In the present study, liquid film thickness is measured by laser focus displacement meter under flow boiling condition and compared with the correlation proposed for an adiabatic flow. The relationship between liquid film thickness and heat transfer coefficient is also investigated. Initial liquid film thickness under flow boiling condition can be predicted well by the correlation proposed under adiabatic condition. Under flow boiling condition, liquid film surface fluctuates due to high vapor velocity and shows periodic pattern against time. Frequency of periodic pattern increases with heat flux. At low quality, heat transfer coefficients calculated from measured liquid film thickness show good accordance with heat transfer coefficients obtained directly from wall temperature measurements.     

Single-phase pressure drop and heat transfer characteristics of turbulent liquid nitrogen flow in micro-tubes

Experiments have been conducted to investigate the single-phase pressure drop and heat transfer characteristics of liquid nitrogen in four micro-tubes with the diameters of 1.931, 1.042, 0.834 and 0.531 mm. The friction factors are compared with the conventional correlations over a Reynolds number range of 10,000–90,000. The effect of the variable thermal properties of liquid nitrogen, i.e., viscosity and thermal conductivity, on the flow and local heat transfer in the micro-tubes is clarified. The average Nusselt numbers are determined and compared with the correlations for the conventional channels and micro-channels, respectively. It is found that large roughness of the micro-tube causes high friction factor, and the modified Colebrook correlation can well predict the experimental friction factors by using the measured surface roughness. With the increase of liquid nitrogen temperature, the pressure drop decreases as a result of the lower viscosity. Opposite to water, the local heat transfer coefficient of liquid nitrogen flow in the micro-tube drops by 12.5% along the tube. The experimental data show that the average Nusselt numbers for the micro-tubes are higher than those predicted by the correlations for the conventional channels. Taking into account the effect of surface roughness of the micro-tubes on the heat transfer, the modified Gnielinski correlation enables to predict the experimental Nusselt numbers with a mean absolute error (MAE) of 6.4%.