Crisis phenomena in falling liquid films at periodic heat loads

Experimental results on development of crisis phenomena in a falling film of cryogenic liquid at alternating pulse heat release are presented. Experimental data on local temperature evolution along a heat releasing surface are obtained. It is shown that amplitude of heater temperature pulsations depends significantly on the heat flux density and coordinates along the liquid film flow. New experimental data on the critical heat fluxes corresponding to formation of stable dry spots and drying crisis are presented depending on duration of heat release pulses. It is shown that parameters of the forming metastable regular structures and critical drying parameters of the heat-releasing surface are determined by dynamics of the movable boundaries of wetting during self-organization of a system.     

Measuring deformations of the heated liquid film by the fluorescence method

The fluorescence method was used to measure the instantaneous thickness field of the falling nonisothermal water film. The process of rivulet formation in a heated film was registered. Measurement averaging allowed determination of the degree of transverse deformation of a film. In the lower half of the heater within the interrivulet zone of the non-isothermal film, the wave amplitude decreases with a rise of the heat flux and reduction of the average thickness. Two zones of the heat flux effect on liquid film deformation were distinguished. At low heat fluxes, the film flow is weakly deformed. At high heat fluxes the thermal-capillary forces provide formation of rivulets and a thin film between them. The work was financially supported by the President of RF (NSh-6749.2006.8), Russian Foundation for Basic Research (Grants Nos. 05-08-33325-a, 06-01-00360-a), National Center on Science and Innovations (State contract No. 02.438.11.7002), and SB RAS (Interdisciplinary Integration Project No. 111).     

The features of rewetting dynamics of the overheated surface by a falling film of cryogenic liquid

The dynamic process of rewetting of the overheated surface by gravitationally falling film of cryogenic liquid was firstly modeled numerically with consideration of local distribution of heat transfer coefficient in the wetting zone along the 2D front. The front shape corresponding to self-organizing regular structures observed in experiments was obtained in the numerical experiment. Evolution of the front shape was studied. It was shown that local motion velocities of different areas of the 2D wetting front differed significantly. Total time of transitional process was determined by the minimal velocity of evaporating liquid boundaries in the front zones between boiling jets. This model allows quantitative determination for the wetting front velocity, variable in time and space, and temperature fields in the heater. Reliability of calculation results was proved by direct comparison with experimental data.     

Measurement of the wave characteristics of a film under the effect of external perturbations

The effect of artificial perturbations on structure formation in the water film flow over a vertical plate with a heater was studied experimentally. To measure the film thickness an eight-channel capacitance probe was used. It is shown that artificial perturbations on the liquid film surface can change the distance between rivulets from the values corresponding to the thermocapillary-wave regime of rivulet formation to the values related to the thermocapillary regime. The distance between the rivulets can be changed at Reynolds numbers higher than those corresponding to the thermocapillary regime. Artificial perturbations do not change significantly the relative amplitude of large waves, but they affect the character of amplitude dependence on the heat flux density. The work was financially supported by the Russian Foundation for Basic Research (Grants Nos. 05-08-33325-a, 06-01-00360-a) and SB RAS (Interdisciplinary Integration Project No. 111).     

Heat transfer and crisis phenomena with intense boiling in the falling wave liquid films

Experimental data on heat transfer with intense evaporation in the falling films of liquid nitrogen were analysed. According to data generalization, heat transfer at evaporation becomes more intense under the precrisis modes at high heat fluxes for two studied boundary conditions on the heat-releasing surface: T _w ≈ const and q _w ≈ const. The relative contributions of conductive and convective components of heat transfer for different heat fluxes were estimated due to statistical treatment of the wave characteristics carried out by the capacitance probes for measurement of the local liquid film thickness. It was found out that heat transfer intensification is mainly caused by a drastic decrease in thermal resistance of the local zones with intensely evaporating residual layer between large waves. At that, the convective component of heat transfer related to wave perturbations on a free surface of a liquid film decreases significantly with a rise of heat fluxes. New data on pulsations of the local temperature of the heat-releasing surface were obtained at different points along the flow with the modes of “dry spot” formation. The work was financially supported by the Russian Foundation for Basic Research (grant No. 03-02-04027-NNIO-a) and DFG (Deutsche Forschungsgemeinschaft, project Nos. Re-463-37-1 and SFB-540).     

Influence of thermocapillary effects on wave characteristics of the heated liquid film

Wave characteristics of the water film flow over a vertical plate with a heater were studied. Without a heat flux, data obtained perfectly coincide with results of other authors obtained for an isothermal liquid film. When heating the falling liquid, thermocapillary forces form the rivulets and thin film between them. It was found out that an increase in the heat flux causes a rise of the phase velocity and decreases the frequency of 3D waves. It is shown that in contrast to the known data for an isothermal liquid film, the average relative wave amplitude increases with a rise of dimensionless complex Re_loc/Ka _loc ^1/11 at relatively high heat flux densities in the interrivulet. The growth of wave amplitude under the action of thermocapillary forces was revealed in the interrivulet, what corresponds to calculations. The work was financially supported by the Russian Foundation for Basic Research (Grants Nos. 05-08-33325-a and 06-01-00360-a), Federal Target R&D Program (State contracts Nos. 02.438.11.7002 and NSh-6749.2006.8), and SB RAS (International interdisciplinary project No. 111).     

Influence of the plate inclination angle on rivulet formation and breakdown of non-isothermal liquid film

The influence of the plate inclination angle relative to the horizon on a distance between formed rivulets and breakdown of a falling water film was studied experimentally on a heater of 150×150 mm. Dependences of the distance between rivulets on the heat flux density and plate inclination angle were examined. In experiments two zones of influence of the heat flux on a distance between rivulets were revealed: at θ ≥ 20° distances between the rivulets do not depend on the heat flux, and at θ ≤ 15° they decrease with a rise of the heat flux. Data on film breakdown were generalized with consideration of evaporation effect. It is shown that there is almost no effect of inclination angle on film breakdown in the whole range of studied Reynolds numbers. The work was financially supported by the President of RF (No. NSh-6749.2006.8), Russian Foundation for Basic Research (Grants Nos. 05-08-33325-a, 06-01-00360-a) and SB RAS (Interdisciplinary integration project No. 111).     

Interaction of a premixed flame with a liquid fuel film on a wall

In piston engines and in gas turbines, the injection of liquid fuel often leads to the formation of a liquid film on the combustor wall. If a flame reaches this zone, undesired phenomena such as coking may occur and diminish the lifetime of the engine. Moreover, the effect of such an interaction on maximum wall heat fluxes, flame quenching, and pollutant formation is largely unknown. This paper presents a numerical study of the interaction of a premixed flame with a cold wall covered with a film of liquid fuel. Simulations show that the presence of the film leads to a very rich zone at the wall in which the flame cannot propagate. As a result, the flame wall distance remains larger with liquid fuel than it is for a dry wall, and maximum heat fluxes are smaller. The nature of the interaction of flame wall interaction with a liquid fuel is also different from the classical flame/dry wall interaction: it is controlled mainly by chemical mechanisms and not by the thermal quenching effect observed for flames interacting with dry walls: the existence of a very rich zone created above the liquid film is the main mechanism controlling quenching.     

Role of hydroxyl production and heat release in the two-zone fuel-rich adiabatic dimethyl ether/air flames at atmospheric pressure

Fuel-rich laminar adiabatic flames of premixed dimethyl ether/air mixtures at a high initial temperature and atmospheric pressure have been studied by numerical simulation and sensitivity analysis. These flames, having two heat release zones, are of great interest as an unusual and little-studied subject. We have investigated the chemical processes occurring in the two zones and analysed the mechanism of heat release in the flame. It has been found that the key reactions that have a significant influence on the flame speed are those involving dimethyl ether and the products of its incomplete oxidation. Calculation of the heat release rate confirms the presence of two heat release zones in the flame. A comparison of the reactions making a major contribution to the heat release with those significantly affecting the flame speed indicates that the main factor determining the flame speed is the formation of hydroxyls, rather than heat release. Analysis of the flame speed sensitivity shows that in the case of a two-zone structure of the flame, its speed is mainly determined by the reactions taking place in the low-temperature zone. That is, the cool zone with a higher temperature gradient is the leading one.     

液滴撞击液膜的流热耦合界面追踪方法数值模拟

采用界面追踪方法研究蒸馏过程中液滴撞击高温薄液膜的流动和传热特性,将数值结果与解析解和实验进行比较验证模型的正确性,研究气液界面和热流分布的演变过程.同时,分析液滴We数和无量纲液膜厚度对传热的影响.液滴撞击后的热流密度分布显示：液膜可分为撞击区、过渡区和静态区.由于液滴的撞击作用,强制对流是撞击区内主要的传热机制.增大液滴的韦伯数或减小无量纲液膜厚度会加强热量传递.随着液滴韦伯数的增加,冲击引起的扰动增强,在动量和能量共同作用下,平均热流密度明显增大,撞击区冠状水花越明显.无量纲液膜厚度越小,平均热流密度越大,且有更长的时间保持高热流密度换热.     

Investigation of transient processes at liquid boiling under nonstationary heat generation conditions

The paper presents results on experimentally investigated dynamics of boiling development and formation of film boiling zones under stepwise heat generation on a horizontally and vertically oriented cylindrical surface in a large volume of Freon-21. Experimental data on the expectation time and boiling temperature, the propagation velocity and structure of evaporation and boiling fronts for different heat flux density both in saturated liquid and in subcooling conditions are obtained. Results of experiments on investigating the nucleation forms under development of nonstationary heat release crisis caused by heat loading on the vertical heater immersed into the volume of liquid (water, ethanol) subcooled to saturation temperature are presented. A calculation ratio for determining the expectation time to the beginning of intense vaporization in water is proposed and compared with experimental data obtained on surfaces with different-size roughness. Peculiarities of evolution of evaporation fronts from incipient bubbles are investigated in the experiments with ethanol. Data on the evaporation front velocity as a function of wall overheating are obtained. The obtained experimental data on the propagation velocity of self-sustained evaporation fronts are compared with the known calculated data.     

Influence of the intensity of chemical heat-release source on the heat exchange at a flow of pseudo-plastic liquid in the initial interval of coaxial duct

The critical regimes of heat exchange at a laminar steady flow of a pseudo-plastic liquid in the initial interval of a coaxial duct are investigated with regard for both dissipative and chemical sources in the Arrhenius representation [5] under the conditions of an insignificant variation of the concentrations of reacting substances. The work was financially supported by the Russian Foundation for Basic Research (Grant No. 05-08-50043), FANI (State Contract No. 02.434.11.5009), and the Presidium of RAS (Program P-09).     

Growth criterion for small dry spots in the falling liquid films

In thin films of liquid falling over an inclined surface random (or non-random) small dry spots shut and disappear or expand into the large-scale dry spots depending on the initial spot sizes and parameters of the liquid film. On consideration of the gravity forces, inertial forces of the film and surface tension, affecting the elements of a roller surrounding the dry spot, we have formulated a criterion, which helped us to answer the following question: if this spot shuts or not? This criterion includes the numbers of Bond and Reynolds and the contact wetting angle as the determining parameters. It is generalization of the known Hartley — Murgatroid criterion, which allows determination of the zone of falling film metastability, i. e., determination of the critical film thickness, below whose value the formation of stable dry spots is possible. It is shown that the critical values of parameters of the initial dry spot depend on the shape and sizes of a roller surrounding this spot.     

Numerical Simulation for Falling Film Flow Characteristics of Refrigerant on the Smooth and Structured Surfaces

Flow characteristics of a liquid film flowing over a smooth surface and structured surface with the Reynolds number range from 10 to 1121 are studied. The mixture of R21 and R114 refrigerants is used as the test liquid. The 3D transient simulations are taken to capture the liquid film’s dynamic characteristics and spatial distribution. Effects of the inlet dimension, inlet flow rates, surface tension, and surface structuring on the wettability, average velocity, and film thickness are studied systematically. The obtained results show that surface tension is essential for an accurate simulation, while inlet width has no effect on the liquid film parameters in the steady-state flow regime. For low flow rates, wetting area and film thickness both are small, and a suggested range of Reynolds number is chosen to simulate further heat transfer in order to balance the film thickness and dry spots generation. It is shown that a ripple surface structure hinders the liquid film movement, reflected in a lower velocity and a larger film thickness compared to the smooth surface. Lateral movement of a liquid film can also be observed at the structured surface.     

Thermodynamic crisis of boiling

The procedure of calculation of the dynamics of rapid near-wall vaporization near a metallic heater is considered. A physical model of explosive boiling-up on bubbles of fluctuation origin is used. The model is limited to regimes when the motion of bubbles can be ignored and convective flows do not have enough time to develop. It is assumed that the dry spot under a bubble thermally insulates the heater wall. The heat removal is provided by the vicinity of the wetting line (WL). An analytical calculation has become possible on the basis of the well-known exact solution to the problem of the temperature field near the wetting line. Generalization of this solution has led to two new problems. These are the problem of allowance for the difference of the dynamic wetting angle from the right angle and that of allowance for thermocapillary flows. A numerical comparison of the results of calculation of the dry area with the use of different methods of allowing for the “tightness” effect in the dynamics of bubble generation and growth has been made. The evolution of bubbles from their generation to growth limited by the heat supply has been investigated. It has been found that the transition stage in the development of bubbles from Rayleigh to thermal ones is of considerable importance in the processes of explosive boiling-up. The dynamics of change in the dry area and wetting line length prior to the stage of bubble merging into a vapor film has been calculated. A condition of passage of the heat flux through a maximum is found. The applicability of the idea of thermodynamic crisis to calculations of miniature devices is justified. The problem of constructing a model to calculate the interphase surface shape near the wetting line under a developed thermocapillary flow and considerable reactive forces of the vapor flow is formulated. The model is in good agreement with the results of experiments on pulsed superheating of liquids at a rate higher than 1 K/μs.     

Heat transfer and crisis phenomena at boiling in Freon mixture films falling down the structured tube

Investigation results on hydrodynamics, heat transfer, and crisis phenomena in laminar-wave liquid films falling down the surfaces of different geometry are presented in this study. Freon mixture R21/R114 with initial concentration of low-boiling component of 4?C16.6 % was used as the working liquid. The film Reynolds number at the inlet to the experimental section varied from 60 to 700. The heat flux density was changed within 0?C5 W/cm². The images of wave surface of the falling liquid film and formation of dry spots were visualized and recorded by the high-speed video camera. Results of investigation of the wave structure of the film surface, measurements of heat transfer coefficients under the conditions of boiling, and critical heat fluxes in the film flow over the smooth and structured surfaces are presented.     

Vibration characteristics of aluminum surface subjected to ultrasonic waves and their effect on wetting behavior of solder droplets

The vibration characteristics of an aluminum surface subjected to ultrasonic waves were investigated with a combination of numerical simulation and experimental testing. The wetting behavior of solder droplets on the vibrating aluminum surface was also examined. The results show that the vibration pattern of the aluminum surface is inhomogeneous. The amplitude of the aluminum surface exceeds the excitation amplitude in some zones, while the amplitude decreases nearly to zero in other zones. The distribution of the zero-amplitude zones is much less dependent on the strength of the vibration than on the location of the vibration source. The surface of the liquid solder vibrates at an ultrasonic frequency that is higher than the vibration source, and the amplitude of the liquid solder is almost twice that of the aluminum surface. The vibration of the surface of the base metal (liquid solder) correlates with the oxide film removal effect. Significant removal of the oxide film can be achieved within 2 s when the amplitude of the aluminum surface is higher than 5.4 μm or when the amplitude of the liquid solder surface is higher than 10.2 μm.     

三维倾斜平板液氧降膜润湿数值模拟研究

建立三维倾斜平板降膜模型,利用VOF两相流模型计算了液氧降膜的润湿情况,研究了工质物性、倾斜角、液膜入口高度对润湿面积的影响。结果表明:Weber数(We)相同时,液氧和水的润湿比均随Kapitza数(Ka)增大而减小;相同Ka下,液氧和水的润湿比均随We增大而增大,而液氧润湿比一直小于水润湿比,两者的差值也随We增大而增大。拟合得到液氧在液膜入口高度0.4 mm、接触角70°时的界面润湿比经验关联式,拟合值和模拟值相对误差≤±20%;在We=0.76时,液氧的润湿比随倾斜角增大而减小,但降低速率随Ka增大而减小;在倾斜角为90°时,易出现液膜脱离壁面的现象;当We固定时,液氧的润湿比随液膜入口高度增大而增大。     

Interconnection between dynamics of liquid boiling-up and heat transfer crisis for nonstationary heat release

This paper deals with investigation results on crisis phenomena for nonstationary heat release under the conditions of free convection and in falling liquid films. It is shown that the character of the crisis development and the critical heat flux for nonstationary heat release significantly depend on the characteristics determining the temperature head of liquid boiling-up. According to experimental data with an arbitrary temporal function of heating, the character of the dependence between the critical heat flux and the heat release increasing rate is significantly effected by ready nucleation sites. It is found that a change in the nonstationary critical heat flux in the range of high times between impulses for periodic heat release is connected with deactivation of ready nucleation sites on the heat-releasing surface. According to new experimental results, in the studied range of irrigation degree alteration (Re_in = 30–1660), parameters characterizing decay of the falling liquid film with stepped heat release (the distribution of the time of boiling-up expectation along the liquid film, the velocities of movable boundaries in the boiling-up and drying fronts, the intensity of liquid ejection from the heat-releasing surface) complexly depend on the Reynolds number, wave characteristics, and heat flux density. The text was submitted by the authors in English.     

Disintegration of flows of superheated liquid films and jets

This paper represents results on investigating the dynamics of boiling and disintegration of superheated liquid films and jets. The first part deals with experimental study of boiling liquid outflow through short cylindrical and slit channels. Evolution of disintegration of a hot water jet flow is observed both at low and moderate superheating and at high and limit superheating, and also for vaporization mechanisms corresponding to these superheatings. Peculiarities of disintegration of jets through slit and cylindrical channels are noticed. Results on measuring the reactive thrust of the jet through a slit channel under different geometrical conditions behind the channel outlet are represented. The 1/f fluctuations in transient regimes of superheated liquid boiling and in transient regimes of behavior of the jet shape are found. The second part of this article represents results on experimental investigation of nonsteady heat transfer and dynamics of the development of crisis phenomena at boiling of a falling subcooled liquid film in the conditions of stepwise heat release. The experimental data were obtained using synchronized high-speed infrared thermography and video. It is shown that with growth and condensation of vapor bubbles, on the liquid film interface appear large-amplitude waves that lead to considerably increasing local intensity of heat transfer. New data on the boiling incipience temperature in a subcooled liquid film, depending on the heat flux density, are obtained. It is found that the development of boiling crisis is a result of appearance of local dry patches and their subsequent growth by the mechanism of longitudinal thermal conductivity in the heat transfer wall as the equilibrium heat flux density is exceeded.