Thin-film phase change driven by disjoining pressure difference

The thin liquid film at the contact line is gaining increasing attention due to its importance in phase-change heat transfer and wettability control. The intermolecular effects within the thin film are usually represented by of disjoining pressure. In the present work, molecular dynamics were employed to investigate the influence of disjoining pressure on thin–film evaporation and condensation in a nanoscale triple-phase system. The simulation domain was a cuboid with argon gas sandwiched between the two solid platinum walls. The two solid walls were fixed at the same temperature while the liquid films had different initial thicknesses thus corresponding to different disjoining pressures. Spontaneous evaporation and condensation were observed at the thicker and the thinner film, respectively. Disjoining pressure together with thermodynamics theories were employed to qualitatively and quantitatively explain the phenomenon. The evaporation fluxes were measured and compared to the Hertz-Knudsen-Schrage model which is based on the kinetic theory of gases. The resulting non-evaporating thickness was measured and compared to the models based on disjoining theory.     

Effect of vapor flow on the wetting behavior of unstable evaporating menisci in heated capillary tubes

The wicking height of a heated, evaporating meniscus formed by surface-wetting liquid in a vertical capillary tube with dynamic flow has been investigated. Previous experimental results and analytical models for measuring/predicting wicking heights in capillaries are also reviewed. An analytical model is presented that accounts for both major and minor vapor pressure losses along the vertical capillary tube. It is shown that during thermo-mechanical instability, vapor/meniscus interaction can become more prevalent due to increased vapor generation/pressure near the meniscus free surface. A relatively simple procedure for estimating onset of meniscus instability is presented and, when used with the vapor Reynolds number, can estimate whether vapor pressure loss is significant. By comparing the current model with the available experimental data, it is shown that the wicking height of an unstable, evaporating meniscus of n-pentane in a vertical, glass capillary tube is better estimated by considering vapor flow pressure losses – providing a 40% improvement over previous models that neglect vapor flow. In addition to vapor flow pressure loss, the dynamic contact angle and thin film profile must also be calculated to ensure accurate prediction of wicking height. Although the proposed model shows improvement, it is prone to under-predicting the actual meniscus wicking height for stable, evaporating menisci at lower relative heat loads. The proposed model can be used for predicting wicking behavior of heated, vertically-aligned liquid columns in capillary structures – which is relevant to the design of miniature heat transfer equipment/media such as wicked heat pipes, micro-channels and sintered/porous surfaces.     

Effect of disjoining pressure on the instability of a charged thin liquid film on a spherical rigid core

It is shown that the critical Rayleigh number which characterizes the stability of a thin charged viscous fluid film on the surface of a rigid spherical core develops rapidly with decrease in the film thickness to 100 nm when the effect of the disjoining pressure becomes significant. The dependence of the instability growth rate on the thickness of the fluid layer is obtained by analyzing the dispersion relation numerically. Yaroslavl’. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 102–106, January–February, 1999.     

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.     

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 Variational Model of Disjoining Pressure: Liquid Film on a Nonplanar Surface

Variational methods have been successfully used in modelling thin liquid films in numerous theoretical studies of wettability. In this article, the variational model of the disjoining pressure is extended to the general case of a two-dimensional solid surface. The Helmholtz free energy functional depends both on the disjoining pressure isotherm and on the shape of the solid surface. The augmented Young–Laplace equation (AYLE) is a nonlinear second-order partial differential equation. A number of solutions describing wetting films on spherical grains have been obtained. In the case of cylindrical films, the phase portrait technique describes the entire variety of mathematically feasible solutions. It turns out that a periodic solution, which would describe wave-like wetting films, does not satisfy Jacobi’s condition of the classical calculus of variations. Therefore, such a solution is nonphysical. The roughness of the solid surface significantly affects liquid film stability. AYLE solutions suggest that film rupture is more likely at a location where the pore-wall surface is most exposed into the pore space, and the curvature is positive.     

Flow boiling visualization in a vertical circular minichannel at high vapor quality

This paper reports on an experimental study of saturated flow boiling of R134a inside a circular vertical quartz tube coated with a transparent heater. The inner diameter of the tube was 1.33 mm and the heated length 235.5 mm. The flow pattern at high vapor qualities and the dryout of the liquid film were studied using a high speed CCD camera at the mass fluxes 47.4 and 124.4 kg/m² s in up flow at 6.425 bar. The heat fluxes ranged from 5 to 13.6 kW/m² for the lower mass flux and from 20 to 32.4 kW/m² for the higher mass flux.

The behavior of the flow close to dryout was found to be different at low and high mass flux. At low mass flux the location of the liquid front fluctuated with waves passing high up in the tube. In between the waves, a thin film was formed, slowly evaporating without breaking up.

At high mass flux the location of the liquid front was more stable. In this case the liquid film was seen to break up into liquid streams and dry zones on the tube wall.     

磁场力及膜曲率对磁敏感薄膜-基底界面黏附性能的影响与调控

界面黏附和脱黏的可调控在攀爬装置、黏附开关、机械抓手等方面具有重要的应用需求. 针对磁敏感薄膜-基底界面, 开展了薄膜初始曲率及外加磁场对界面黏附性能影响机制的研究. 首先实验制备了具有初始曲率的磁敏感薄膜, 分别开展了具有初始曲率的磁敏感薄膜-基底界面撕脱实验及理论研究, 研究了薄膜初始曲率、弯曲刚度和外加磁场强度对界面黏附性能的影响规律. 实验和理论结果一致表明: 具有初始曲率的磁敏感薄膜-基底界面黏附力随薄膜初始曲率的增大而减小, 而外加磁场能够有效提高界面黏附力;相比于初始零曲率薄膜-基底界面稳态撕脱力与薄膜弯曲刚度无关, 薄膜弯曲刚度减弱了具有初始曲率薄膜-基底界面的稳态撕脱力. 进一步从能量角度分析了界面等效黏附性能, 揭示了薄膜弯曲能、磁场势能、界面黏附能的相互竞争机制. 最后, 基于本文的实验及理论结果, 提出了一种磁场和薄膜初始曲率协同调控的简易机械抓手, 可连续实现物体的拾取、搬运和释放功能. 本文结果不仅有助于理解多场调控的界面可逆黏附机制, 对界面黏附可控的功能器件设计亦提供了一种新方法.      

An updated three-zone heat transfer model for slug flow boiling in microchannels

This work proposes a novel physics-based model for the fluid mechanics and heat transfer associated with slug flow boiling in horizontal circular microchannels to update the widely used three-zone model of Thome et al. (2004). The heat transfer model has a convective boiling nature and predicts the time-dependent variation of the local heat transfer coefficient during the cyclic passage of a liquid slug, an evaporating elongated bubble and a vapor plug. The capillary flow theory, extended to incorporate evaporation effects, is applied to estimate the bubble velocity along the channel. A liquid film thickness prediction method also considering bubble proximity effects, which may limit the radial extension of the film, is included. The minimum liquid film thickness at dryout is set to the channel wall roughness. Theoretical heat transfer models accounting for the thermal inertia of the liquid film and for the recirculating flow within the liquid slug are utilized. The heat transfer model is compared to experimental data taken from three independent studies. The 833 slug flow boiling data points cover the fluids R134a, R245fa and R236fa, and channel diameters below 1 mm. The proposed evaporation model predicts more than 80% of the database to within ±30%. It demonstrates a stronger contribution to heat transfer by the liquid slugs and correspondingly less by the thin film evaporation process compared to the original three-zone model. This model represents a new step towards a complete physics-based modelling of the bubble dynamics and heat transfer within microchannels under evaporating flow conditions.     

接触角滞后与流体动压润滑的相关性研究

理论研究表明不同润湿性界面对流体动压润滑油膜厚度有着显著地影响，一般采用接触角(CA)来表征固液界面润湿性. 而由热力学原理推导出的界面势能垒理论模型不仅与接触角相关，也是接触角滞后(CAH)的函数. 本文作者通过对不同基体材料的滑块进行表面张力修饰，获得了不同亲和性的界面. 利用干涉法及荧光法分别测量了不同润湿性界面的流体动压润滑膜厚及油膜受剪切的流动特性，研究了接触角及接触角滞后两个界面参数对流体动压润滑油膜厚度的影响，并对势能垒与接触角滞后的关系进行了讨论. 结果表明:接触角与流体动压润滑油膜厚度的相关性较差，接触角滞后可以更好地表征界面效应对流体动压润滑油膜厚度的影响.      

STABILITY OF ULTRATHIN LIQUID FILM EVOLUTION WITH SURFACTANT

针对固体基底上厚度小于100 nm的含活性剂超薄液膜演化过程, 基于润滑理论推导出包含分离压影响的液膜厚度和活性剂浓度的演化方程, 采用正则模态法导出了描述液膜线性稳定性的特征方程, 分析了多个特征参数对线性稳定性的影响, 数值模拟了液膜厚度和活性剂浓度演化历程, 对比了模拟所得非线性结果与线性分析预测结果的一致性.结果表明:范德华力具有促进扰动增长的作用, 较强的玻恩斥力促使扰动衰减, 使液膜趋于稳定;较小的毛细力数易使液膜凹陷处发生二次失稳, 并最终导致去润湿现象发生;液膜厚度和溶于液膜内部的活性剂浓度初值越大, 液膜稳定性越强, 液膜表面活性剂浓度影响则相反;增大吸附系数不利于液膜稳定性.     

Extension of Stoney's formula to non-uniform temperature distributions in thin film/substrate systems. The case of radial symmetry

Current methodologies used for the inference of thin film stress through curvature measurements are strictly restricted to stress and curvature states which are assumed to remain uniform over the entire film/substrate system. By considering a circular thin film/substrate system subject to non-uniform, but axisymmetric temperature distributions, we derive relations between the film stresses and temperature, and between the plate system's curvatures and the temperature. These relations featured a “local” part which involves a direct dependence of the stress or curvature components on the temperature at the same point, and a “non-local” part which reflects the effect of temperature of other points on the location of scrutiny. Most notably, we also derive relations between the polar components of the film stress and those of system curvatures which allow for the experimental inference of such stresses from full-field curvature measurements in the presence of arbitrary radial non-uniformities. These relations also feature a “non-local” dependence on curvatures making full-field measurements of curvature a necessity for the correct inference of stress. Finally, it is shown that the interfacial shear tractions between the film and the substrate are proportional to the radial gradients of the first curvature invariant and can also be inferred experimentally.     

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.     

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.     

Turbulent boundary layer on a mildly curved convex surface

Extensive single point turbulence measurements made in the boundary layer on a mildly curved heated convex wall show that the turbulence heat fluxes and Stanton number are more sensitive to a change in wall curvature than the Reynolds stresses and skinfriction coefficient, and that downstream, as the flow adjusts to new curved conditions, the St/c _f ratio of Reynolds analogy is appreciably lower than in plane wall flow for the same conditions. Details of the turbulence structure in unheated flow have been documented in an earlier paper; temperature field measurements now described comprise mean temperature distributions, the streamwise variation of wall heat flux, profiles of the temperature variance, transverse and streamwise heat fluxes, and triple correlations. Turbulent diffusion of heat flux is drastically reduced even by mild curvature; changes in the heat fluxes are of the same order as changes in the shear stress, that is, an order of magnitude greater than the ratio of boundary layer thickness to wall radius of curvature. The data include plane flow measurements taken in a developed boundary layer upstream of a change in wall curvature.     

Thermal characteristics of a high heat flux micro-evaporator

Our purpose is to design a high heat flux micro-evaporator that can remove more than 100 W/cm². For this purpose a thin liquid film is evaporized. The liquid film is stabilized in micro-channels by capillary forces. The micro-channels are fabricated by chemical etching on silicon to reduce thermal resistance. For the experiments, the channel plate is heated by an ITO thin film heater deposited on the opposite side of the channel plate. Influence of heat flux, coolant flow rate, and inlet temperature on the temperature of the heater element are investigated. Water is used as working fluid. A maximal heat flux of 125 W/cm² could be achieved for water inlet temperature of 90 °C and flow rate of 1.0 mL/min. The temperature of the heater element is kept constant at about 120 °C with fluctuations within 8 °C. The measured pressure drop is less than 1000 Pa.     

Non-uniform, axisymmetric misfit strain： in thin films bonded on plate substrates/substrate systems： the relation between non-uniform film stresses and system curvatures

Current methodologies used for the inference of thin film stress through curvature measurements are strictly restricted to stress and curvature states which are assumed to remain uniform over the entire film/substrate system. By considering a circular thin film/substrate system subject to non-uniform, but axisymmetric misfit strain distributions in the thin film, we derived relations between the film stresses and the misfit strain, and between the plate system’s curvatures and the misfit strain. These relations feature a ‘‘local’’ part which involves a direct dependence of the stress or curvature components on the misfit strain at the same point, and a ‘‘non-local’’ part which reflects the effect of misfit strain of other points on the location of scrutiny. Most notably, we also derived relations between the polar components of the film stress and those of system curvatures which allow for the experimental inference of such stresses from full-field curvature measurements in the presence of arbitrary radial non-uniformities. These relations also feature a ‘‘non-local’’ dependence on curvatures making a full-field measurement a necessity. Finally, it is shown that the interfacial shear tractions between the film and the substrate are proportional to the radial gradients of the first curvature invariant and can also be inferred experimentally.     

纳米压痕法测磁控溅射铝薄膜屈服应力

为了在考虑残余应力下测量出磁控溅射铝薄膜的屈服应力,提出了一种实验测量方法,通过曲率测试法和球形压头纳米压痕法测出磁控溅射铝薄膜的屈服应力．建立球形压痕力学模型,并用ANSYS对球形压痕进行力学有限元仿真,利用直流磁控溅射技术在硅基上淀积一层1 μm厚的铝薄膜,首先通过曲率测试法测量膜内等双轴残余应力,再利用最小二乘曲线拟合法从薄膜/基底系统的球形压头纳米压痕实验数据中提取出铝薄膜的屈服应力,测得磁控溅射铝薄膜的屈服应力为371 MPa．该方法也可以用来研究其他材料的薄膜和小体积材料的力学特性．     

Algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow

The study considers algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow, focusing in particular on momentum and heat transfer (so-called ‘convective boiling’) through the annular liquid film. In contrast with single-phase wall-bounded flow theory, shear-driven annular liquid films are assumed here to behave as fluid-bounded flows, mostly interacting with the shearing gas-entrained droplets core flow. Besides providing velocity and temperature profiles through the liquid film, the turbulence model proposed here predicts key parameters such as the average liquid film thickness, the void fraction and the convective boiling heat transfer coefficient with accuracies comparable or better than those of leading design correlations. This turbulence model is part of a unified annular flow modeling suite that includes methods to predict the entrained liquid fraction and the axial frictional pressure gradient. The underlying heat transfer database covers nine fluids (water, two hydrocarbons and six refrigerants) for vertical and horizontal tubes of 1.03-14.4 mm diameter and pressures of 0.1-7.2 MPa. Importantly, this study shows that there appears to be no macro-to-microscale transition when it comes to annular flow. Simply better physical modeling is required to span this range.