饱和蒸发和热毛细力作用下低雷诺数液膜在波纹竖壁上的流动

考虑表面蒸发压力和热毛细力作用情况下,对饱和蒸发状态下低雷诺数自由降落液膜在小波幅正弦型波纹壁面上的流动进行理论分析。对控制微分方程及边界条件进行量纲一化并引入流函数,对微分方程及边界条件进行摄动展开,得到了这种情况下液膜流动的简化分析模型,求出了近似解析解。讨论了壁面波纹、表面张力、蒸发压力、热毛细力对液膜流动的影响。研究表明:液膜的波动幅度随蒸发强度和热毛细力的增大而增大;液膜波动与壁面波纹的相位差随蒸发强度增大而增大,随热毛细力增大而减小。     

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

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

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

不同壁面条件下液滴撞击铺展特性的模拟研究

液滴撞击不同润湿性壁面的传热流动问题在自然界和工业生产中广泛存在。研究采用CLSVOF方法,引入描述壁面润湿特性的动态接触角,并考虑液滴物性参数随温度的变化,建立液滴撞壁模型,模拟研究液滴撞击流动行为,通过与实验对比验证,确定模型有效性。在此基础上,对传热作用下考虑壁面润湿性的液滴撞击问题展开研究,探讨壁面传热作用对液滴撞击铺展特性的影响。研究表明,在撞击过程中,液滴先铺展后逐渐收缩,与静态接触角模型相比,采用动态接触角模型所得的液滴流动特性与实验结果更加吻合;随着接触角增大,液滴在撞壁初期不易铺展,随后则易于收缩;虽然固液传热作用会影响液滴铺展直径,但不改变液滴的运动趋势。     

上海陆域古河道分布及对工程建设影响研究

采用直接数值模拟方法, 对槽道湍流中展向振动流向传播的波动壁面的流动 控制和减阻问题进行了研究, 讨论了流向参数k_{x}对Stokes层、湍流拟序结构、湍流猝 发事件以及壁面阻力的影响, 并对此类波动壁面的湍流控制和减阻机理进行了讨论. 结果表 明, 当此类波动壁面被用来调制近壁流动时, 仅低频波对湍流流场具有显著影响, 可导致湍 流猝发事件的频率和强度的显著变化; 波数k_{x}的增大对于湍流猝发事件的频率和强度增 减的影响并不同步, 存在一个最优的波数k_{x}, 在其调制下, 固有流场对诱导流场的影响 最弱, 而诱导流场对固有流场的影响显著, 减阻效果最好.     

基于FTM的Kelvin-Helmholtz不稳定性衍化数值分析

采用FTM研究了两相流体界面中的K-H不稳定性,分析了初始振幅、波数、表面张力、密度比、速度梯度、韦伯数We对K-H不稳定性发展的影响。结果表明:增大初始扰动对K-H不稳定性具有促进作用,扰动越强,K-H波浪翻卷高度越高;波数的增加能加快K-H不稳定性的衍化过程,但对K-H波浪高度的发展影响较小;减小表面张力或密度比都能使不稳定性的发展速度增加;初始速度梯度对不稳定性的发展有较大影响,梯度越大,K-H波浪翻卷高度越高;We对K-H波浪的对称性影响较大,We越小,对称程度越高。     

用格子Boltzmann方法模拟液滴撞击固壁动力学行为

首次用格子Boltzmann方法中的伪势模型对液滴撞击固壁的动力学行为进行了数值模拟.详细研究了液滴在壁面上的流动状态以及各种因素对撞击过程的影响.通过数值模拟得到:壁面的可润湿性越小,液滴越容易发生反弹,液滴的回缩速度越快;液滴的撞击速度越大,所得到的相对直径越大,回缩速度越快;液滴的粘性越小,所得到的相对直径越大;液滴的表面张力越大,液滴越容易发生反弹现象.另外,液滴的最大相对直径与We数满足一定的线性关系,这些结果与前人的理论预测和实验结果完全吻合.     

粗糙多孔腔体自然对流的非正交MRT-LB数值模拟

基于非正交多松弛系数格子Boltzmann(MRT-LB)方法建立了适用于多孔方腔自然对流计算模型,选取典型热流动问题分析了非正交转换矩阵的MRT-LB模型数值准确性和运算效率,布置两种粗糙多孔方腔模型进行研究。讨论了Rayleigh数、粗糙单元个数n、粗糙单元横纵比A等参数对方腔内流动传热的影响,同时给出了各个参数条件下腔内流场与温度场分布。结果表明,非正交MRT-LB模型具有很好的数值准确性和收敛速度,布置粗糙单元使得壁面处流线等温线发生变形,增大粗糙元单元个数n或横纵比A均会恶化方腔传热。     

浮力对混合对流流动及换热特性的影响

用热线和冷线相结合的技术测量垂直圆管内逆混合对流流体的平均速度、 温度以及它们的脉动. 较详细地研究了浮力对逆混合对流的流动特性和传热特性的影响. 评 估了实验中采用的冷线测量温度补偿速度探头温度敏感的影响. 逆混合对流的传热结果用无 量纲参数Ω (Ω= Grd / Red2 )来表示，其中，基于管道直 径的雷诺数Red变化范围为900~18000, 浮力参数Ω变化范围为 0.004899~0.5047. 研究结果表明，浮力对逆混合对流的换热有强化作用. 随着葛拉晓夫数Grd的增加，温度脉动，流向雷诺正应力和流向温度通量增 大，并且在靠近壁面的流体区域尤其明显. 热线与冷线相结合的技术适合于研究非绝热的流 动测量，可以用于研究浮力对流动和换热特性的影响.     

正弦波温度边界下多孔介质方腔内非热平衡对流传热数值模拟

采用局部非热平衡模型,在方腔左侧壁面温度正弦波变化、右侧壁面温度均一的边界条件下,通过SIM-PLER算法数值研究了固体骨架发热多孔介质方腔内的稳态非达西自然对流,主要探讨了不同正弦波波动参数N及方腔的高宽比M/L对方腔内自然对流与传热的影响规律。计算结果表明:正弦波温度边界使得方腔内的流场出现了复杂的变化,流体及固体区域左侧壁面附近出现了周期性的正负变化的温度场分布,左侧壁面局部Nusselt数出现了周期性的震荡现象;存在一个最佳温度波动参数N=1,此时多孔介质方腔内的整体散热量达到最大值;增加方腔高宽比会显著地削弱方腔内的自然对流传热过程,小高宽比也会在一定的程度上削弱多孔介质方腔内的对流传热。     

Numerical analysis on flow dynamics and heat transfer of falling liquid films with interfacial waves

Flow dynamics and heat transfer of falling liquid films with interfacial waves flowing on a vertical plate have been studied with originally proposed numerical simulation method. To discretize basic equations a staggered grid fixed on a physical space is employed. A small amplitude disturbance generated at inflow boundary develops to a solitary wave which consists of a large amplitude roll wave and small amplitude capillary waves. Instantaneous streamwise velocity profiles at the wave crest and trough are very different from a laminar flow. A circulation flow occurs in the roll wave and it affects temperature distributions, especially the strong effect is observed for high Prandtl number liquids. The interfacial wave enhances the heat transfer by two kinds of effects which are a film thinning effect and a convection effect. The dominating effect depends on the Prandtl number. Received on 23 December 1998     

Temperature pulsations and thermocapillary effects on the surface of a wavy heated liquid film

The temperature pulsations and wave characteristics in water film flow along a vertical plate with a heater are investigated. Using an infrared scanner, the temperature field on the film surface is measured for various heat flux densities on the heater. Experimental data on the variation of the temperature with time on a local segment of the liquid film surface during wave transmission are obtained. In the absence of a heat flux the data obtained are in good agreement with the results of other researchers for an isothermal liquid film. When the down-flowing liquid is heated, the thermocapillary forces lead to the formation of rivulets and a thin film between them. It is shown that in the inter-rivulet zone the relative wave amplitude increases due to the action of the thermocapillary forces.     

Heat and mass transfer of a water film falling down a tilted plate with radiant heating and evaporation

This paper has concerned the heat and mass transfer of a water film falling down a tilted plate with radiant heating and water evaporation. A cluster of physical models was developed for evaluating the properties of heat and mass transfer. A fully implicit control-volume finite-difference procedure was used to solve the coupling equations. The effects of various parameters on heat and mass transfer were investigated. The results showed that the mass fraction of water vapor in ambient atmosphere and the flow turbulence played key roles in the heat and mass transfer. The ambient atmospheric temperature dramatically affected the sensible heat flux. However its effect on the latent heat flux is negligibly small. The magnitude of solar incident flux had an intense influence on the water film temperature. Received on 29 January 1998     

Effects of Interphase Heat Transfer and Surface Tension on Nonlinear Wave Propagation in a Gas-Liquid Mixture

An evolutionary equation with two nonlinearities (basic hydrodynamic and an additional nonlinearity attributable to interphase heat transfer) describing the wave dynamics in a mixture in the presence of viscosity is obtained. For a quasi-adiabatic regime, depending on the heat transfer laws determined by the relations between the thermophysical parameters of the mixture two exact particular solutions predicting different pressure behaviors in the wave are given. Nonmonotonous pressure behavior associated with an increase in the wave amplitude is interpreted as wave amplification as compared with the case of monotonous behavior. Depending on the thermophysical parameters, an interval of permissible Weber numbers is found for which the amplification effect is observed. The surface tension reduces the wave amplitude and width and results in a nonmonotonous dependence of the wave velocity on the Weber number.     

Evaporation of pure liquids with increased viscosity in a falling film evaporator

The present study investigated fluid dynamics and heat transfer of viscous pure liquids in a falling film evaporator. This is of special benefit as it avoids mass transfer effects on the evaporation behaviour. Experiments at a single-tube glass falling film evaporator were conducted. It allowed a full-length optical film observation with a high-speed camera. Additionally the evaporator was equipped with a slotted weir distribution device. Test fluids provided viscosities ranging from μ = 0.3 to 41 mPa s. The Reynolds number was between 0.7 and 1,930. Surface evaporation and the transition to nucleate boiling were studied to gain information about the film stability at maximum wall superheat. A reliable database for laminar and laminar-wavy viscous single component films was created. The experimental results show a significant enhancement in the wave development due to the film distribution. A wavy flow with different wave velocities was superposed to the film in each liquid load configuration without causing a film breakdown or dry spots on the evaporator tube. It was found that nucleate boiling can be allowed without causing film instabilities over a significant range of wall superheat.     

Numerical simulation of convective heat transfer to a radial free surface jet impinging on a hot solid

The convective heat transfer between a circular free surface impinging jet and a solid surface has been studied numerically. The thin liquid film formed on the surface has been assumed to be in non-turbulent free surface flow. The effects of surface tension, viscosity, gravity and heat transfer between the film flow and the solid surface have been taken into account. The flow structure on a non-heated surface has been investigated first. Next, the steady-state flow structure in the liquid film as well as the heat transfer has been examined. The predicted results have been compared with experimental data for the purpose of validating the analysis. The hydrodynamics of the liquid film and the heat transfer processes have been investigated numerically to understand the physics of the phenomena. Received on 5 October 1998     

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

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

Study of the thermocapillary effect on a wavy falling film using a fiber optical thickness probe

The wave flow of a water film down a vertical plate with a 150×150 mm heater has been experimentally studied. The effect of the heat flux on the film flow leads to the formation of periodically flowing rivulets and thin film between them due to the action of thermocapillary forces in spanwise direction. The local film thickness between rivulets is measured by means of a noncontact fiber optical probe. As the heat flux grows, the average film thickness continuously decreases but upon reaching about 50% of the initial thickness, the film spontaneously breaks down. It is found that the decrease of the wave amplitude between rivulets is caused by the reduction of the local Reynolds number and is in a qualitative agreement with the laws of the hydrodynamics for the isothermal case. That is, no appreciable effect of streamwise thermocapillary forces on the wave amplitudes is detected. The experimental results are in good agreement with recently published data obtained by the capacitance method.     

Wärmeübergang bei der Fallfilmverdampfung

At film evaporation in a vertical tube the heat transfer depends on the flow pattern of the falling film and the shear stress at film surface due to a significant flow in the gas phase. For heat transfer without superimposed gas flow equations exist which are checked experimentally. These equations have general validity for film evaporation and film condensation. By describing the influence of gas flow great differences occur between condensation and evaporation. The present work deals with a calculation method which was verified by condensation data. This method is compared and discussed with experimental results of film evaporation. It is shown that the presented method describing the shear stress influence is also applicable.