Thermovibrational Instability Of Rayleigh-Marangoni-Benard Convection In Two-Layer Fluid Systems

The thermovibrational instability of Rayleigh-Marangoni-Benard convection in a two-layer system under the high-frequency vibration has been investigated by linear instability analysis in the present paper. General equations for the description of the convective flow and within this framework, the generalized Boussinesq approximation are formulated. These equations are dealt with using the averaging method. The theoretical analysis results show that the high-frequency thermovibrations can change the Marangoni-Benard convection instabilities as well as the oscillatory gaps of the Rayleigh-Marangoni-Benard convection in two-layer liquid systems. It is found that vertical high-frequency vibrations can delay convective instability of this system, and damp the convective flow down. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Numerical simulation of drop Marangoni migration under microgravity

Thermocapillary motion of a drop in a uniform temperature gradient is investigated numerically. The three-dimensional incompressible Navier-Stokes and energy equations are solved by the finite-element method. The front tracking technique is employed to describe the drop interface. To simplify the calculation, the drop shape is assumed to be a sphere. It has been verified that the assumption is reasonable under the microgravity environment. Some calculations have been performed to deal with the thermocapillary motion for the drops of different sizes. It has been verified that the calculated results are in good agreement with available experimental and numerical results. (C) 2003 Elsevier Ltd. All rights reserved.

Numerical investigation of thermocapillary migration of the drop for large Marangoni numbers

An axisymmetric model is adopted to simulate the problem of unsteady drop thermocapillary motion for large Marangoni numbers. Front tracking methods are used in the investigation. It is found that the non-dimensional drop migration velocity will decrease with increasing Marangoni number. This agrees well with the experimental results obtained from the 4th Shen-Zhou space ship. In the meanwhile, this is also the first time for numerical simulations to verify the experimental phenomenon under large Marangoni numbers.

Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.

微重力条件下悬浮区中的Marangoni对流

本文研究微重力条件下悬浮区中的表面张力梯度驱动对流。采用有限元方法数值研究二维、轴对称悬浮区中的定常流动和温度分布,讨论P_r数及自由面上散热条件的影响。计算表明,P_r数(M_α数)增大时,流区大部分自由表面上温度梯度减小,流区流速较小,对流热输运较扩散效应增强。只要能维持自由面上一定的温度梯度,在较大的P_r数(M_α数)时,仍能得到收敛的定常解。

How to determine critical Marangoni number in half floating zone convection

In the present paper, the coordinated measurements of the temperature profile inside the liquid bridge and the boundary variation of Free surface, in addition to other quantities, were obtained in the same time for the half floating zone convection. The results show that the onset of free surface oscillation is earlier than the one of temperature oscillation during the increasing of applied temperature difference, and the critical Marangoni numbers, defined usually by temperature measurement, are larger than the one defined by free surface measurement, and the difference depends on the volume of liquid bridge. These results induce the question, ''How to determine experimentally the critical Marangoni number?'' Copyright (C) 1996 Elsevier Science Ltd.

Effects of high-frequency vibration on critical marangoni number

The influence of vibration on thermocapillary convection and critical Marangoni number in liquid bridge of half floating zone was discussed for the low frequency range 0.4-1.5 Hz and the intermediate frequency range 2.5-15 Hz in our previous papers. This paper extends the study to high frequency range 15-100Hz. This ground based experiment was completed on the deck of an electromagnetic vibration machine. The results of our experiment shows when the frequency of the applied acceleration is high enough, the amplitude of the time varying part of the temperature response is disappear and the shape of the free surface of the liquid bridge exhibits no fluctuations due to inertia. The critical Marangoni number which is defined to describe the transitions from a peroidical convection in response to vibration to an oscillatory convection due to internal instability is nearly the same as the critical Marangoni number for oscillatory flow in the absence of vibration.

Effects of g-jitter on the critical marangoni number

A half floating zone is fixed on a vibrational deck, which supports a periodical applied acceleration to simulate the effect of g-jitter. This paper deals with the effects of g-jitter on the fluid fields and the critical Marangoni number, which describes the transition from a forced oscillation of thermocapillary convection into an instability oscillatory convection in a liquid bridge of half floating zone with top rod heated. The responses of g-jitter field on the temperature profiles and flow pattern in the liquid bridge were obtained experimentally. The results indicated that the critical Marangoni number decreases with the increasing of g-jitter effect and is slightly smaller for higher frequency of g-jitter with fixed strength of applied gravity.

大Marangoni数下液滴热毛细迁移的数值研究

采用轴对称的非定常模型对大Marangoni数下的液滴热毛细迁移问题进行数值模拟研究，采用界面追踪方法，发现无量纲的液滴迁移速度会随着Marangoni数的增大而减小．数值模拟结果与神舟4号飞船实验中的结果进行了比较，得到一致的变化趋势．同时，这也是第一次在大Marangoni数范围下利用数值模拟的方法证实该实验现象．

大Reynolds数液滴Marangoni迁移的空间实验装置

讨论了液滴热毛细迁移的空间实验装置，该实验是中国神舟号飞船中的一个项目。与国际上同类装置比较，它具有不需要宇航员直接操作，领先程序自动控制及遥科学手段完成全部科学实验的能力。它不仅能完成大Reynolds数的液滴迁移实验，而且在这类空间装置上使用了等厚型光学干涉系统。

过冷核态池沸腾中的Marangoni效应

分析了过冷核态池沸腾过程中气泡横向运动现象,指出其成因在于相邻气泡界面温度差引起的Mara- ngoni对流对周围液体的吸引.通过对该流动的尺度分析,得到了气泡横向运动特征速度及其可观测度的估算公式,其预测结果与实验观测相一致.特别是对极细小的初始核化气泡,该公式预测了强烈的横向Marangoni对流会导致气泡顶端微射流的形成.该效应在加热面水平向下或微重力沸腾等气泡脱落频率较低的情形中尤为重要.

纵向非均匀多孔介质Marangoni对流不稳定性

理论研究了纵向非均匀多孔介质中流体表面张力驱动的对流不稳定性、充满液体的多孔介质层从下方加热，上方自由表面冷却，形成可引起多孔介质液层Marangoni—Benard对流流动的纵向温度梯度．采用线性化的Brinkman．Forchheimier方程作为控制方程组，对孔隙率分别为线性函数、正弦三角函数分布的非均匀多孔介质液层的Marangoni—Benard问题进行了线性稳定性分析、通过采用Chebyshev-Tau谱方法求解广义特征值问题，得到了系统临界Marangoni数随无量纲波数变化的中性稳定性曲线，分析和比较了孔隙率的变化对液层对流稳定性和流场结构的影响，获得了纵向非均匀多孔介质液层不稳定性现象的新特征．

Oscillatory Instability of Rayleigh-Marangoni-Benard Convection in Two-Layer System

The convective instabilities in two or more superposed layers heated from below were studied extensively by many scientists due to several interfacial phenomena in nature and crystal growth application. Most works of them were performed mainly on the instability behaviors induced only by buoyancy force, especially on the oscillatory behavior at onset of convection (see Gershuni et. Al.(1982), Renardy et. Al. (1985,2000), Rasenat et. Al. (1989), and Colinet et. Al.(1994)) . But the unstable situations of multi-layer liquid convection will become more complicated and interesting while considering at the same time the buoyancy effect combined with thermocapillary effect. This is the case in the gravity reduced field or thin liquid layer where the thermocapillary effect is as important as buoyancy effect. The objective of this study was to investigate theoretically the interaction between Rayleigh-Bénard instability and pure Marangoni instability in a two-layer system, and more attention focus on the oscillatory instability both at the onset of convection and with increasing supercriticality. Oscillatory behavious of Rayleigh-Marangoni-Bénard convective instability (R-M-B instability) and flow patterns are presented in the two-layer system of Silicon Oil (10cSt) over Fluorinert (FC70) for a larger various range of two-layer depth ratios (Hr=Hupper/Hdown) from 0.2 to 5.0. Both linear instability analysis and 2D numerical simulation (A=L/H=10) show that the instability of the system depends strongly on the depth ratio of two-layer liquids. The oscillatory instability regime at the onset of R-M-B convection are found theoretically in different regions of layer thickness ratio for different two-layer depth H=12,6,4,3mm. The neutral stability curve of the system displaces to right while we consider the Marangoni effect at the interface in comparison with the Rayleigh-Bénard instability of the system without the Marangoni effect (Ma=0). The numerical results show different regimes of the developing of convection in the two-layer system for different thickness ratios and some differences at the onset of pure Marangoni convection and the onset of Rayleigh-Bénard convections in two-layer liquids. Both traveling wave and standing wave were detected in the oscillatory instability regime due to the competition between Rayleigh-Bénard instability and Marangoni effect. The mechanism of the standing wave formation in the system is presented numerically in this paper. The oscillating standing wave results in the competition of the intermediate Marangoni cell and the Rayleigh convective rolls. In the two-layer system of 47v2 silicone oil over water, a transition form the steady instability to the oscillatory instability of the Rayleigh-Marangoni-Bénard Convection was found numerically above the onset of convection for ε=0.9 and Hr=0.5. We propose that this oscillatory mechanism is possible to explain the experimental observation of Degen et. Al.(1998). Experimental work in comparison with our theoretical findings on the two-layer Rayleigh-Marangoni-Bénard convection with thinner depth for H<6mm will be carried out in the near future, and more attention will be paid to new oscillatory instability regimes possible in the influence of thermocapillary effects on the competition of two-layer liquids