波纹竖板层流降膜蒸发传热的分析

对高粘度液体在等温正弦形波纹壁面上的自由降落与蒸发建立了摄动分析模型。得到了流动的分析解和蒸发传热的数值解。考察了壁面波纹的波幅和波数、液膜表面张力及贝克利数对流动与传热的影响,结果表明,加大波纹的波幅、适当选择波数、减小贝克利数可增强传热,而表面张力对蒸发传热的影响较小。     

EFFECT OF TWO-DIMENSIONAL MICROPILLAR ARRAYED TOPOGRAPHY ON SPREADING OF INSOLUBLE SURFACTANT-LADEN DROPLET

针对二维微柱阵列壁面上含不溶性活性剂液滴的铺展过程,采用润滑理论建立了液膜厚度和浓度演化模型,采用数值计算方法得到了液滴的铺展特征及相关参数的影响. 研究表明：活性剂液滴在微柱阵列壁面上铺展时,在壁面凸起处衍生出隆起结构,壁面凹槽处衍生出凹陷结构,随时间持续,隆起和凹陷均向两侧移动,且数量不断增加. 活性剂液膜流经凸起时,隆起高度呈驼峰形变化. 增大预置液膜厚度或活性剂初始浓度,铺展区域隆起和凹陷数量增多,液滴铺展速度加快. 增加凹槽深度或减小斜度会使毛细力作用增强,液膜破断可能性加大;增大凹槽宽度可加速活性剂液滴的铺展,加剧液膜表面波动幅度.     

二维微柱阵列壁面对活性剂液滴铺展的影响

针对二维微柱阵列壁面上含不溶性活性剂液滴的铺展过程,采用润滑理论建立了液膜厚度和浓度演化模型,采用数值计算方法得到了液滴的铺展特征及相关参数的影响. 研究表明:活性剂液滴在微柱阵列壁面上铺展时,在壁面凸起处衍生出隆起结构,壁面凹槽处衍生出凹陷结构,随时间持续,隆起和凹陷均向两侧移动,且数量不断增加. 活性剂液膜流经凸起时,隆起高度呈驼峰形变化. 增大预置液膜厚度或活性剂初始浓度,铺展区域隆起和凹陷数量增多,液滴铺展速度加快. 增加凹槽深度或减小斜度会使毛细力作用增强,液膜破断可能性加大;增大凹槽宽度可加速活性剂液滴的铺展,加剧液膜表面波动幅度.      

双自由面溶质?热毛细液层的不稳定性

溶质?热毛细对流是流体界面的浓度和温度分布不均导致的表面张力梯度驱动的流动, 它主要存在于空间微重力环境、小尺度流动等表面张力占主导的情况中, 例如晶体生长、微流控、合金浇筑凝固、有机薄液膜生长等. 对其流动进行稳定性分析具有重要意义. 本文采用线性稳定性理论研究了双自由面溶质?热毛细液层对流的不稳定性, 得到了两种负毛细力比(η)下的临界Marangoni数与Prandtl数(Pr)的函数关系, 并分析了临界模态的流场和能量机制. 研究发现: 溶质?热毛细对流和纯热毛细对流的临界模态有较大的差别, 前者是同向流向波、逆向流向波、展向稳态模态和逆向斜波, 后者是逆向斜波和逆向流向波. 在Pr较大时, Pr增加会降低流动稳定性; 在其他参数下, Pr增加会增强流动稳定性. 在中低Pr, 溶质毛细力使流动更加不稳定; 在大Pr时, 溶质毛细力的出现可能使流动更加稳定; 在其他参数下, 溶质毛细力会减弱流动稳定性. 流动稳定性不随η单调变化. 在多数情况下, 扰动浓度场与扰动温度场都是相似的. 能量分析表明: 扰动动能的主要能量来源是表面张力做功, 但其中溶质毛细力和热毛细力做功的正负性与参数有关.      

饱和下降液膜的稳定性研究

通过简化表面张力、蒸汽压力和液膜波动之间的关系,用线性理论分析了饱和液膜在等温竖壁上的流动稳定性,讨论了雷诺数、波数、壁面和液膜温差及流体物性变化的影响。     

纳米阵列中气体驱替液体的流动特征

纳米尺度下气体驱动液体流动特征在纳流控芯片及页岩气开发中具有广泛的应用前景. 利用管径规格为292.8 nm,206.2 nm,89.2 nm,67.0 nm,26.1 nm的氧化铝膜为纳米阵列,进行气驱水实验和单相气体流动实验,分析纳米尺度下气驱水流动特征. 实验表明,纳米阵列中气驱水时气体流量随驱动压力变化经历三个阶段:第一阶段流量缓慢增大,且比单相气体流量降低约一个数量级;第二阶段纳米阵列中的水被大量驱替出,流量迅速增大;第三阶段纳米阵列中的水全部被驱替出,流动特征与单相气体流动保持一致. 分析表明,气驱水第一阶段存在气液界面毛细管力的“钉扎”作用及固液界面相互作用力的影响,是产生非线性流动的主要原因;而一旦“钉扎”作用破坏,气体进入管道推动界面运动,气柱与液柱之间的毛细曲面曲率变化,毛细管力减小,气体流量急剧增大,其中毛细管力随驱替压力增大急剧变化,是造成第二阶段气体流量突变的主要原因.      

纳米阵列中气体驱替液体的流动特征

纳米尺度下气体驱动液体流动特征在纳流控芯片及页岩气开发中具有广泛的应用前景.利用管径规格为292.8 nm,206.2 nm,89.2 nm,67.0 nm,26.1 nm的氧化铝膜为纳米阵列,进行气驱水实验和单相气体流动实验,分析纳米尺度下气驱水流动特征.实验表明,纳米阵列中气驱水时气体流量随驱动压力变化经历三个阶段:第一阶段流量缓慢增大,且比单相气体流量降低约一个数量级;第二阶段纳米阵列中的水被大量驱替出,流量迅速增大;第三阶段纳米阵列中的水全部被驱替出,流动特征与单相气体流动保持一致.分析表明,气驱水第一阶段存在气液界面毛细管力的"钉扎"作用及固液界面相互作用力的影响,是产生非线性流动的主要原因;而一旦"钉扎"作用破坏,气体进入管道推动界面运动,气柱与液柱之间的毛细曲面曲率变化,毛细管力减小,气体流量急剧增大,其中毛细管力随驱替压力增大急剧变化,是造成第二阶段气体流量突变的主要原因.     

考虑毛细作用微梁剥离的实验和力学模型

毛细力诱发的粘附现象在自然界和工农业生产中广泛存在,例如微机电系统、微纳米自组装、油气驱替等．本文系统研究了两根微梁的毛细粘附行为,包括梁剥离过程中液桥的形貌以及剥离力-位移变化规律．试验发现,微梁在毛细力作用下的剥离部分经历了液膜粘附和液滴粘附两个阶段．考虑两个阶段的液桥形状特征,分别建立系统的能量泛函,采用变分原理推导了考虑毛细力的微梁剥离的非线性微分方程和边界条件．基于Matlab 编程求解方程,得到了剥离力-位移曲线,理论计算与试验结果具有很好的一致性．另外,参数研究表明,接触角和表面张力系数对液膜粘附的微梁剥离影响显著,而对液滴粘附的剥离影响较小;微梁刚度对两个阶段的剥离都有明显影响．本文的试验结果和理论分析对于实际工程中微结构的定量设计具有一定参考价值．     

气泡夹带对壁面油膜流动特性的影响

油-气润滑系统工作过程中,润滑油膜受微油滴冲击和压缩空气扰动影响易形成气泡夹带现象,气泡夹带行为将对壁面润滑油膜层的形成及流动过程产生重要影响。基于VOF数值模拟方法,对含气泡油膜沿倾斜壁面的流动行为进行研究,考察了气泡的存在对油膜形态和流动速度的影响规律,以及气泡破裂阶段空腔邻域内流体压力变化特性。研究表明,油膜夹带气泡的形变和迁移诱发气泡周围微流场的速度扰动现象,导致气液界面处产生非均匀速度梯度分布,进而引发油膜表面的形态波动。气泡发生破裂时,油膜空穴部位发生明显的正负压力波动现象,气泡附近壁面将承受一定的交变载荷作用。     

大尺寸液桥热毛细对流失稳性地面实验研究

开展大尺寸液桥浮力-热毛细对流地面实验, 探究流场转捩的临界条件及临界状态附近的流动情况. 通过粒子图像测速方法(PIV) 获得流体速度场, 研究液桥内部定常和转捩后的流场结构以及流体运动规律;并用红外热像仪测量液桥自由面温度分布, 研究流体流动的时空演化和温度振荡. 实验发现大尺寸半浮区液桥浮力-热毛细对流临界值与几何参数有关, 在大普朗特(Prandtl) 数情况下, 流场存在由稳定态向不稳定态再到混沌的转捩过程, 在临界马兰哥尼(Marangoni) 数附近, 流场内会出现行波现象, 流动模式也会随高径比的变化而发生变化;当继续增大马兰哥尼数, 流动会进入混沌状态.      

薄膜沿加热平板下落的稳定性分布

本文分析了薄膜沿加热平板下落的稳定性。在时间模式下，发现流动的不稳定性是由表面波不稳定和加毛细不稳定构成的，同时当流体的热扩散越大以及界面热量损失越小时，热毛细不稳定越剧烈，在时空模式下，流动随着Marangoni数的增大。流动有可能从对流不稳定过渡到绝对不稳定，这一结论尚待实验验证。     

An experimental investigation on the effects of surface gravity waves on the water evaporation rate in different air flow regimes

Estimating rate of evaporation from undisturbed water surfaces to moving and quiet air has been the topic a vast number of research activities. The obvious presence of various shapes of gravity waves on the water body surfaces was the motivation of this experimental investigation. In this investigation experimental measurements have been done to quantify evaporation rate from wavy water surfaces in free, mixed and forced convection regimes. The effects of a wide range of surface gravity waves from low steepness, round shaped crest with slow celerity, to steep and very slight spilling crest waves, on the water evaporation rate have been investigated. A wide range of ${\text{Gr}}/{\text{Re}}^{2} (0.01 \le {\text{Gr}}/{\text{Re}}^{2} \le 100)$ was achieved by applying different air flow velocities on a large heated wave flume equipped with a wind tunnel. Results reveal that wave motion on the water surface increase the rate of evaporation for all air flow regimes. For free convection, due to the effect of wave motion for pumping rotational airflows at the wave troughs and the dominant effect of natural convection for the air flow advection, the maximum evaporation increment percentage from wavy water surface is about 70 %. For mixed and forced convection, water evaporation rate increment is more sensitive to the air flow velocity for the appearance of very slight spilling on the steep wave crests and the leeward air flow structures.     

Evaporation of Falling and Shear-Driven Thin Films on Smooth and Grooved Surfaces

One of the most important tasks in development of modern gas turbine combustors is the reduction of NO_x emissions. An effective way to reduce the NO_x emission is using the lean premixed prevaporization (LPP) concept. An important phenomenon taking place in LPP chambers is the evaporation of thin fuel films. To increase the fuel evaporation rate, the use of microstructured walls has been suggested. The wall microstructures make use of the capillary forces to evenly distribute the liquid fuel over the wall, so that the appearance of uncontrolled dry patches can be avoided. Moreover, the wall structures promote the thin film evaporation characterized by ultra-high evaporation rates. An experimental setup was built for the investigation of thin liquid films falling down on the outer surface of vertical tubes with either a smooth or structured surface. In the first testing phase water is used, fuel like liquids will be used later on. The thin film can be heated from both sides, by hot oil flowing inside the tube, and by hot compressed air flowing in co-current direction to the thin film. The film is partly evaporated along the flow. Results for the wavy film structure at different Reynolds numbers are reported. For theoretical investigations a model describing the hydrodynamics and heat transfer due to evaporation of the gravity- and shear-driven undisturbed liquid film on structured surfaces was developed. For low Reynolds numbers or low liquid mass fluxes the wall surface is only partly covered with liquid and the heat transfer is shown to be governed by the evaporation of the ultra-thin film in the vicinity of the three-phase contact line. A numerical model for the solution of a two-dimensional free-surface flow of a liquid film over a structured wall was also developed. The Navier–Stokes equations are solved using the Volume of Fluid (VOF) technique. The energy equation is included in the model. The model is verified by comparison with data from the literature showing favorable agreement. In particular, the proposed model predicts the formation of capillary waves observed in the experiments. The model is used to investigate the flow of liquid on a structured wall. This calculation is the first step towards the modeling of a three-dimensional wavy flow of a gravity- and shear-driven film along a wall with longitudinal grooves. It is found that due to the Marangoni effect, a circulating flow arises within the cavity, thereby leading to an enhancement in the evaporation rate.     

Numerical investigations on vortical structures of viscous film flows along periodic rectangular corrugations

Two-dimensional gravity-driven film flows along a substrate with rectangular corrugations are studied numerically by using Finite Volume Method. The volume of fluid (VOF) method is utilized to capture the evolution of free surfaces. The film flows down an inclined plate are simulated to validate the numerical implementation of the present study. Results obtained indicate that the phase shift between the surface wave and the wall corrugation increases as the Reynolds number. The parametric studies on the interesting resonant phenomenon indicate that the peak Reynolds numbers increase as the raise of the wall depth or the decline of the inclination angle. The dependence of the flow fields is analyzed on the Reynolds numbers and wall depth in details. It is found that the vortical structures in the steady flows, either produced by the interaction between capillary wrinkling and inertia, or by the rectangular geometry, are closely related to the remarkable deformation of the free surfaces. This conclusion is also confirmed by the transient flow development of two typical simulations, i.e., flows in capillary–inertial regime and in inertial regime.     

Investigations of the Marangoni effect on the regular structures in heated wavy liquid films

Formation and development of quasi-regular metastable structures within laminar-wavy falling films were studied. These structures emerge within the residual layer between large waves and could be one reason for the break up of the falling film. The temperature field of the film surface was visualised using IR-thermography. The film thickness was obtained from point measurements with the chromatic confocal imaging method and converted into a film thickness field, based on a quasi-steady assumption and IR thermography images. The thermo-capillary nature (Marangoni effect) of the regular structures was proven experimentally.     

PIV measurements of flow in drying polymer solutions during solvent casting

An experimental method based on confocal microscopy and particle image velocimetry (PIV) is used to characterize the flow in a polymer solution during solvent casting. The flow inside a 200-μm-thick film of a poly(vinyl alcohol) (PVA) solution is visualized near a vertical wall of a mold using confocal microscopy of seed particles during solvent evaporation at 25, 35, and 45°C, and the corresponding velocity vector fields are determined from projections of the confocal images. Flow toward the vertical wall is observed inside the film as well as a slower Marangoni-type counter flow at the film surface during the initial phase of solvent evaporation, resulting from a polymer concentration gradient along the film due to a local variation in evaporation rate. Total volume of the polymer solution in the observation volume as well as solvent evaporation rate are determined as a function of time, both revealing close correlation to average horizontal velocity data from PIV. The PIV measurements show significant differences in the flow velocity fields at different temperatures. The PIV measurements correlate with the solvent evaporation rates as well as the final polymer thicknesses on the vertical wall of the mold. Surface tension and viscosity measurements are taken for different concentrations of PVA solution.     

Nonlinear vibrations of a vapor film in the presence of intense heat fluxes

Waves propagating along the interface between a thin vapor film and a liquid layer in the presence of a heat flux are investigated. The boundary conditions on the vapor-liquid phase surface take into account the temperature dependence of the pressure and the possibilities of formation of the metastable state of the superheated liquid and mass flow. Variations in the saturation pressure as functions of the temperature and mass flux lead to generation of weakly damped periodic waves of low amplitude whose velocity can be much higher than the velocity of the gravity waves. The waves ensure stability of the vapor film beneath the liquid layer in the gravity field. The finite-amplitude waves on the surface of the vapor film differ from the Stokes surface waves on the free surface of isothermal fluid. Instability regimes related with superheating of the liquid ant its explosive boiling when the amplitude of an initially small wave increases to infinity in a finite time can develop in a certain working-parameter regime.     

生物芯片中周期性电渗驱动液体薄膜的流动特性

研究了二维周期性电渗驱动液体薄膜的流动特性. 以Debye-Hückel 假设近似下线性化的Poisson-Boltzmann方程描述双电层电动势分布和电荷密度的分布关系, 与黏性不可压缩流体Navier-Stokes方程相耦合, 得到流体在自由面与固壁之间的周期电渗流流场的精确解. 结果显示, 薄膜内速度振幅与流体黏性密切相关, 雷诺数越大, 速度振幅就越小. 该文还细致分析了雷诺数和自由面ζ电势对自由面的流速振幅和薄膜内速度相位差的影响.     

Viscous liquid film flow down an inclined corrugated surface. Calculation of the flow stability to arbitrary perturbations using an integral method

Viscous liquid film flow along an inclined corrugated (sinusoidal) surface has been studied. Calculations were performed using an integral model. The stability of nonlinear steady-state flows to arbitrary perturbations was examined using the Floquet theory. It has been shown that for each type of corrugation there is a critical Reynolds number for which unstable perturbations occur. It has been found that this value greatly depends on the physical properties of the liquid and geometric parameters of the flow. In particular, in the case of film flow down a smooth wall, the critical waveformation parameter depends only on the angle of inclination of the flow surface. The values of the corrugation parameters (amplitude and period) were obtained for which the film flow down a wavy wall is stable to arbitrary perturbations up to moderate Reynolds numbers. Such parameter values exist for all investigated angles of inclination of the flow surface.