两层流体热毛细对流空间实验研究=Space Experiments of Thermo Capillary Convection in Two-liquid Layers

1999年,在我国实践5号卫星上完成了两层流体空间实验,实验研究两层不相混合流体的纯Marangoni对流（温度梯度与界面垂直）与热毛细对流（温度梯度方向与流体界面平行）.前者存在发生Marangoni对流的最小临界温差值△Tc,低于该值流体系统处于静止状态;后者中只要存在沿界面的温度梯度便会产生热毛细对流.空间实验采用石蜡和氟化液两层流体新体系,实现了平整的液-液交界面,并从卫星上传回上万幅数字图像.通过多幅图像叠加处理得到了定量的流速场.数值模拟计算分析了相应工况时对流流动的速度场,两者的流场结构和速度大小基本一致,实验验证了理论模型.

Space experiments of thermocapillary convection in two-liquid layers

In 1999, the space experiments on the Marangoni convection and thermocapillary convection in a system of two immiscible liquid layers in microgravity environment were conducted on board the Chinese scientific satellite SJ-5. A new system of two-layer liquids such as FC-70 liquid and paraffin was used successfully, with the paraffin melted in the space. Two different test-cells are subjected to a temperature gradient perpendicular or parallel to the interface to study the Marangoni convection and thermocapillary convection, respectively. The experimental data obtained in the first Chinese space experiment of fluid are presented. Two-dimensional numerical simulations of thermocapillary convections are carried out using SIMPLEC method A reasonable agreement between the experimental investigation and the numerical results is obtained.

Experimental Study On Liquid Free Surface In Buoyant-Thermocapillary Convection

We investigate the surface deformations of buoyant-thermocapillary convection in a rectangular cavity clue to gravity and temperature gradient between the two sidewalls. The cavity is 52mm x 42mm in horizontal cross section, the thickness of liquid layer h is changed from 2.5mm to 6.5mm. Surface deformations of h = 3.5mm and 6.0mm are discussed and compared. Temperature difference is increased gradually, and the flow in the liquid layer will change from stable convection to unstable convection. Two kinds of optical diagnostic system with image processor are developed for study of the kinetics of buoyant-thermocapillary convection, they give out the information of liquid free surface. The quantitative results are calculated by Fourier transform and correlation analysis, respectively. With the increasing temperature gradient, surface deformations calculated are more declining. It is interesting phenomenon that the inclining directions of the convections in thin and thick liquid layers are different. For a thin layer, the convection is mainly controlled by thermocapillary effect. However, for a thick layer, the convection is mainly controlled by buoyancy effect. The surface deformation theoretically analysed is consistent with our experimental results. The present experiment proves that surface deformation is related to temperature gradient and thickness of the liquid layer. In other words, surface deformation lies on capillary convection and buoyancy convection.

Experimental Study on Bénard-Marangoni Convection in Double Immiscible Liquid Layer

An experimental investigation of Bénard-Marangoni convection has been performed in double immiscible liquid layers of rectangular configuration on the ground. The two kinds of liquid are 10cst silicon oil and FC-70 respectively. The size of rectangular chamber is 100mm×40mm in horizontal cross-section. The evolution processes of convection are observed in the differential thickness ratio of two liquid layers. The critical temperature difference was measured via the detections of fluid convection by a particle image velocimetry (PIV) in the vertical cross-section of the liquid layer. The critical temperature difference or the critical Marangoni number was given. And the influence of the thickness ratio of two liquid layers on the convection instability was discussed. The evolution processes of patterns and temperature distributions on the interface are displayed by using thermal liquid crystal. The velocity distributions on the interface were also obtained. In comparison with the thermocapillary effect, the effect of buoyancy convection will relatively increase when the depth of the liquid layer increases. Because of the coupling of buoyancy and thermocapillary effect, the convection instability is much more complex than that in the microgravity environment. And the critical convection depends on the change of the thickness of liquid layers and also the change of thickness ratio of two liquid layers.

A New Criterion Number for the Boundary Conditions at the Solid/Liquid Interface in Nanoscale

A simple, but important three-atom model was proposed at the solid/liquid interface, leading to a new criterion number, lambda, governing the boundary conditions (BCs) in nanoscale. The solid wall is considered as the face-centered-cubic (fcc) structure. The fluid is the liquid argon with the well-known LJ potential. Based on the concept, the two micro-systems have the same BCs if they have The same criterion number. The degree of the locking BCs is enhanced when lambda equals to 0.757. Such critical criterion number results in the substantial epitaxial ordering and one, two, or even three liquid layers are locked by the solid wall, depending on the coupling energy scale ratio of the solid and liquid atoms. With deviation from the critical criterion number, the flow approaches the slip BCs and there are little ordering structures within the liquid. Always at the same criterion number, the degree of the slip is decreased or the locking is enhanced with increasing the coupling energy scale ratio of the solid and liquid atoms. The above analysis is well confirmed by the molecular dynamics (MD) simulation. The slip length is well correlated in terms of the new criterion number. The future work is suggested to extend the present theory for other microstructures of the solid wall atoms and quasi-LJ potentials.

Coupling of thermocapillary convection and evaporation effect in a liquid layer when the evaporating interface is open to air

The coupling mechanism of thermocapillary convection and evaporation effect in evaporating liquids was studied experimentally. The experiments were carried out to study a thin evaporating liquid layer in a rectangular test cell when the upper surface was open to air. By altering the imposed horizontal temperature differences and heights of liquid layers, the average evaporating rate and interfacial temperature profiles were measured. The flow fields were also visualized by PIV method. For comparison, the experiments were repeated by use of another two non-evaporating liquids to study the influence of evaporation effect. The results reveal evidently the role that evaporation effect plays in the coupling with thermocapillary convection.

Experimental Study on Benard-Marangoni Convection by PIV and TCL

An experimental investigation of Benard-Marangoni convection has been performed in double immiscible liquid layers of rectangular configuration. The two kinds of liquid are 10cst silicon oil and FC-70 respectively. The velocity fields in the vertical cross-section are obtained by PIV. Flow patterns and/or temperature distributions on the horizontal interface are displayed by using thermal color liquid crystal (TLC), and the velocity distributions on the interface were also obtained with the help of the serial particle image of TCL. The evolution processes of convection are observed in the differential thickness ratio of two liquid layers, and the convection styles are discussed.

The Hydrothermal Wave Of Large-Prandtl-Number Fluid In A Shallow Cavity

The hydrothermal wave was investigated numerically for large-Prandtl-number fluid (Pr = 105.6) in a shallow cavity with different heated sidewalls. The traveling wave appears and propagates in the direction opposite to the surface flow (upstream) in the case of zero gravity when the applied temperature difference grows and over the critical value. The phase relationships of the disturbed velocity, temperature and pressure demonstrate that the traveling wave is driven by the disturbed temperature, which is named hydrothermal wave. The hydrothermal wave is so weak that the oscillatory flow field and temperature distribution can hardly be observed in the liquid layer. The exciting mechanism of hydrothermal wave is analyzed and discussed in the present paper.

Fractal features of oscillatory convection in the half-floating zone

Unsteady and two-dimensional numerical simulation is applied to study the transition process from steady convection to turbulence via subharmonic bifurcation in thermocapillary convection of a liquid bridge in the half-floating zone. The results of numerical tests show clearly the fractal structure of period-doubling bifurcations, and frequency-locking at f/4, f/8, f/16 with basic frequency f is observed with increasing temperature difference. The Feigenbaum universal constant is given by the present paper as delta(4) = 4.853, which can be compared with the theoretical value 4.6642016.

Transition from steady to oscillatory convection with chaotic feature in thermocapillary convection

The transition process from steady to turbulent convection via subharmonic bifurcation in thermocapillary convection of half floating zone was studied by numerical simulation and experimental test. Both approaches gave structure of period doubling bifurcations in the present paper, and the Feigenbaum universal law was checked for the system of thermocapillary convection.

G-jitter effects on half floating-zone convection in intermediate-frequency range

The g-jitter influence on thermocapillary convection and critical Marangoni number in a liquid bridge of half-floating rone was discussed in the low frequency range of 0.4 to 1.5 Hz in a previous paper. This paper extended the experiments to the intermediate frequency range of 2 to 18 Hz, which htrs often been recorded as vibration environment of spacecrafts. The experiment was completed on the deck of a vibration machine, which gave a periodical applied acceleration to simulate the effects of g-jitter. The experimental results in the intermediate frequency range are different from that in the low frequency range. The velocity field and the shape of the free surface have periodical fluctuations in response to g-jitter. The amplitude of the periodical varying part of the temperature response decreases obviously with increasing frequency of g-jitter and vanishes almost when the frequency of g-jitter is high enough. The critical Marangoni number is defined to describe the transition from a periodical convection in response to g-jitter to an oscillatory convection due to internal instability, and will increase with increasing g-jitter frequency. According to the spectral analysis, it can be found that the oscillatory part of temperature is a superposition of two harmonic waves if the Marangoni number is larger than a critical value.

Influence of buoyancy on floating-zone convection of small bond number

Gravity may influence the velocity and temperature distributions, bouyancy may induce Rayleigh instability and the instability may be excited due to the change of free surface shape associating with gravity in the thermocapillary convection. These effects have been studied in the present paper. The results show that gravity may have an important effect in thermocapillary oscillatory convection even for the cases of small Bond number experiments either on the ground or in space.

Experimental study on oscillatory thermocapillary convection

The transition processes from steady flow into oscillatory flow in a liquid bridge of the half floating zone are studied experimentally. Two methods of noncontacted diagnoses are developed to measure the distribution of critical Marangoni numbers described by the onset of the oscillation st the free surface of the liquid bridge.The experimental results obtained for both cases of the upper rod heated and the lower rod heated agree with the prediction by Rayleigh's instability theory.The sensitive relations between the relatively fat or slender liquid bridge and the onset of oscillatory convection are also discussed to reveal the insight of the pressure distribution near the free surface. The experiments have been performed in a small liquid bridge, where the Bond number is much smaller than 1, and the results can be used to simulate the experiment in the microgravity environment.

蒸发液层对流稳定性理论与实验研究进展

综述了蒸发液层稳定性理论分析和蒸发对流实验研究方面的最新进展．首先回顾历史上经典的单层流Marangoni不稳定性分析．利用经典的单层流模型可以很好解释无蒸发液层的稳定性特性，但是由于经典的单层流模型没有考虑蒸汽层与液层之间的动力学耦合关系，所以不能完全解释蒸发液层的Marangoni不稳定性特性．有的学者建立了考虑蒸汽层与液层的热耦合与动力学耦合关系的两层流模型，并采用了界面温度连续这样的假设．而在实验的观测中，蒸发界面处的温度是不连续的，特别是在蒸发量比较大的情况下，汽／液界面处温度跳跃很明显．由于界面温度连续假设在处于非平衡状态的系统中是不成立的，所以这些模型虽然能给出一些新的有关系统稳定性的特性，但还是不能完全解释蒸发液层的Marangoni不稳定性的特性，特别是为什么从底部冷却液层的时候，在实验中仍然能够观察到Marangoni对流涡胞的出现的原因．本文总结了前人的研究成果，同时给出了蒸发系统的动力学建模过程和实验研究方法，并对各种模型的稳定性特性进行了总结．最后，指出了现有理论中存在的问题和有待进一步研究的问题．

Influence of Evaporation on Marangoni-Benard Instability

Evaporative convection and instability give rise to both scientific and technological interests. Practically, a number of the industrial applications such as thin-film evaporators, boiling technologies and heat pipes concern with the evaporation process of which through the vapor-liquid interface the heat and mass transfer occur. From a physical viewpoint, one of interesting questions is the mechanisms of convection instability in thin-liquid layers induced by the coupling of evaporation phenomenon and Marangoni effect at the mass exchanged interface. Classical theories, including Rayleigh’s and Pearson’s, have only successfully explained convection in a liquid layer heated from below without evaporation. However these theories are unable to explain the convection in an evaporating thin layer, especially liquid layer is cooled from below. In present paper, a new two-sided model is put forward rather than the one-sided model in previous works. In previous works, the vapor is treated as passive gas and dynamics of vapor has been ignored. In this case, the vapor liquid system can be described by one-sided model. In our two-sided model, the dynamics of vapor should be considered. Linear instability analysis of the Marangoni-Bénard convection in the two-layer system with an evaporation interface is performed. We define a new evaporating Biot number which is different from the Biot number in one-sided model and obtain the curves of critical Marangoni number versus wave number. In our theoretical results, the Biot number and the evaporating velocity play a major role in the stability of the vapor-liquid system.