环形浅液池内低Prandtl数流体超临界热毛细对流演化

为研究低Prandtl数（Pr）流体热毛细对流演化过程，对环形浅液池内Pr=0.011的流体热毛细对流进行三维数值模拟.研究发现:当Marangoni数较小时，流动为轴对称稳态流动；当Marangoni数超过某一临界值后，流动失稳并转变为热流体波，其波数随Marangoni数增加而减小，而波动主频增大；随着Marangoni数增加，流动加强，沿周向运动的热流体波演变为沿径向运动的径向波，其波数大大减小；当Marangoni数继续增加时，波动频谱曲线噪声增加，呈广谱特性.因此，在计算范围内热毛细对流的演化过程为:轴对称稳态流动-热流体波-单周期径向波-多周期三维振荡流动.      

微重力环境下蒸发液层热毛细对流的数值模拟

提出了一种新模型来研究由单一物质构成的液层在其纯蒸气中的蒸发.液层置于微重力环境中并且受到水平方向温度梯度的作用,液层的热毛细对流和蒸发耦合在一起,使得气液界面的传热传质规律更加复杂.用理论分析的方法求解了不考虑热毛细效应的纯蒸发模型,得出温度场分布和界面质量流量的解析表达式.对于热毛细对流和蒸发耦合情况,采用有限差分的投影算法同时求解Navier-Stokes方程和能量方程,得到了不同蒸发Biot数和Marangoni数下流场和温度场的稳态数值解.论述了蒸发Biot数和Marangoni数对界面传热传质的影响,提出并解释了蒸发和热毛细对流耦合的三种模式.      

A sounding rocket experiment on the marangoni convection

During a 6-minutes free flight of a sounding rocket an experiment on the thermal Marangoni convection in a right circular cylindric floating zone was carried out to verify the thermal Marangoni convection by flow visualisation for the first time at 10^?4 g. From the motion picture analysis the velocity distributions are obtained. The interdependence between the two velocity maxima and the S-shaped temperature distribution on the free surface is discussed. The influence of the heat transfer caused by the flow of the surrounding air induced by the thermal Marangoni convection itself is considered briefly.     

Microgravity experiments on the transition from laminar to oscillatory thermocapillary convection in floating zones

The onset of oscillatory thermocapillary convection in floating zones has been measured under microgravity in two sounding rocket experiments TEXUS and on board of the satellite SPAS-01. We report on the critical Marangoni numbers measured under μ-g and in the 1-g reference experiments     

Hydrothermal instability of thermocapillary convection in large-Prandtl-number liquid bridges under microgravity

Linear stability analysis was performed to study the mechanism of transition of thermocapillary convection in liquid bridges with liquid volume ratios ranging from 0.4 to 1.2, aspect ratio of 0.75 and Prandtl number of 100. 2-D governing equations were solved to obtain the steady axi-symmetric basic flow and temperature distributions. 3-D perturbation equations were discretized at the collocation grid points using the Chebyshev-collocation method. Eigenvalues and eigenfunctions were obtained by using the Q–R method.The predicted critical Marangoni numbers and critical frequencies were compared with data from space experiments. The disturbance of the temperature distribution on the free surface causes the onset of oscillatory convection. It is shown that the origin of instability is related to the hydrothermal origin for convections in large-Prandtl-number liquid bridges.     

A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges

A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4  Pr  50). We focus on the relationships between the critical Reynolds number Re_c, the azimuthal wave number m, the aspect ratio Γ and the Prandtl number Pr. A detailed Re_c–Pr stability diagram is given for liquid bridges with various Γ. In the region of Pr > 1, which has been less studied previously and where Re_c has been usually believed to decrease with the increase of Pr, we found Re_c exhibits an early increase for liquid bridges with Γ around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocapillary convection at large Prandtl number, and that the stability property of the basic flow is determined by the “effective” part of liquid bridge.     

Thermocapillary migration of bubbles and droplets

The migration of a droplet in a large liquid body possessing a uniform temperature gradient is analyzed in the creeping flow limit for small values of the Marangoni number (ε). The method of matched asymptotic expansions is used for the solution of the field equations. The migration velocity is calculated to 0(ε²) and the result reduces to that of Young, Goldstein and Block [10] in the limit of zero Marangoni number.     

环形液池热毛细对流的线性稳定性研究

对外壁加热的环形液池热毛细对流进行了线性稳定性分析.采用Chebyshev配点法对Pr=6.8、内外径之比为0.5、深宽比A范围为0.25～1.4的数值结果进行分析,发现流动的临界状态均为振荡形式,并且随着A的增大,临界雷诺数减小,相应的临界波数与振荡频率也呈减小趋势.能量分析结果表明,小扰动与基本流相互作用项较小,表面张力在径向做功与周向做功对小扰动的动能变化起主导作用.观察三者与液池深宽比的关系,发现A=0.8时表面张力在径向做功项达到极小值,周向做功项以及小扰动与基本流相互作用项达到极大值.      

热毛细对流模式转换及温度振荡特征

研究了矩形液池中由于两端温差引起的热毛细对流的温度振荡临界条件.在实验室中,设计了一个高分辨率的温度测量系统,用于实时观测并记录流体的温度.该系统主要由热电偶温度传感器、纳伏表和数据采集电脑3部分组成.得到了各种实验条件下温度振荡的临界条件,并且讨论了它与Prandtl数和Bond数之间的关系.利用flow3d软件数值模拟了微重力条件下的热毛细对流,发现了一种由于自由面变形和液层流场相互作用导致的晃荡的现象.     

Attempt to study Marangoni flow of low-Pr-number fluids using a liquid bridge of silver

In order to experimentally investigate the Marangoni flow of low-Prandtl-number fluids in a liquid bridge geometry under the condition of small Marangoni numbers close to the critical Marangoni numbers Ma_c1 and Ma_c2, the formation of a liquid bridge of silver was attempted. The available temperature difference between the upper and lower rods to obtain a small Marangoni number, such as Ma = 50, was calculated for a 5 mm high liquid bridge for several molten metals. For molten silver, the possible temperature difference was estimated to be 16 K, whereas, for molten silicon, this was 0.38 K, which is unrealistic for the purposes of experiments. For silver, a free surface can be obtained in the wide range of oxygen partial pressures, whereas, for molten silicon, the available oxygen partial pressure range is very small; equilibrium oxygen partial pressure for SiO₂ formation is as low as 1.1 × 10⁻¹⁴ Pa. A liquid bridge of molten silver was successfully prepared and temperature oscillation was observed; the estimated Marangoni number was 160 and oscillation frequency was 0.26 Hz.     

上随体Maxwell(UCM)流体热毛细液层弹性失稳的理论分析

通过对上随体Maxwell（UCM）流体热毛细液层线性稳定性的研究分析,发现流场会发生弹性失稳.扰动的增长速率随波数的增加而增加.与牛顿流体不同,UCM流体不存在临界Marangoni数,当波数达到某一临界值时会出现不稳定的弹性扰动波.该临界波数随弹性数和Marangoni数的增加而减小,当弹性数趋近于0时,流体即变为牛顿流体,而相应的临界波数趋于无穷.不同波数及传播方向上,弹性波的波速相同,而其增长速率在特定方向上达到最大.能量分析表明弹性波的扰动能量来自扰动应力做功.      

Manipulation of Colloidal Droplets in Space and the Instability of Thermocapillary Convection in Large Prandtl Number Liquid Bridge in Microgravity

In 2016 and 2017, SJ-10 and TG-2 satellites were launched. In this short paper, we report recent progress on the studies of manipulation of colloidal droplets and instability of thermocapillary convection in large Prandtl number liquid bridge that based on the space experiments boarding SJ-10 and TG-2 satellites, separately. It was shown that the colloidal droplets can be successfully formed and manipulated in microgravity through the patterned substrates. In another aspect, the coffee-ring effect was observed at the first time in space. For the studies of the instability of thermocapillary convection in large Prandtl number liquid bridge in microgravity, our experiments in TG-2 broadened the way of such kind of study and abundant experimental results are emerging.      

相变材料熔化过程中自然对流振荡现象数值计算研究

为了解自然对流对相变材料熔化过程的影响,对相变材料液相区自然对流振荡现象进行了研究,建立了固/液相统一的控制方程,采用焓-孔隙度法和Boussinesq假设,模拟了重力条件下熔化过程,结果表明:当瑞利数(Ra)较小时,液相自然对流具有稳定解,随着Ra数增大,自然对流从稳定流动向周期振荡过渡。     

Marangoni convection in a liquid layer under the simultaneous action of a transverse magnetic field and rotation

Asymptotic stability characteristics with respect to the onset of Marangoni convection in a liquid layer under the simultaneous action of a large transverse magnetic field and a large rate of rotation are investigated. In certain parameter ranges there is a decrease in critical Marangoni number for increasing magnetic field and rotation depending on the coupling between interfacial perturbations and rotation.     

Marangoni and buoyant convection in an open cavity under reduced and under normal gravity

We give the first full analysis of the microgravity experiment MKB ( arangoni onvektion im offenen oot) conducted on board of the space shuttle during the D1-mission in 1985, together with that of some ground based experiments on thermocapillary flow in a rectangular liquid volume with free upper surface. Streamlines and velocity profiles under μ-g and l-g are presented and discussed. Strong thermocapillary flow was observed under μ-g after the rupture of a dirt-film suppressing surface tension forces. The Nusselt-number under microgravity due to thermocapillary flow for Ma = 4.6 · 10⁵ is at least Nu = 6 ± 1. Under l-g a separation of the surface tension driven convection roll of hot fluid circulating on top of colder fluid was found for large Ma. The coupling and decoupling of buoyant and thermocapillary forces has been further studied in an experiment varying the temperature gradient along the free surface independently from that in the bulk fluid.     

EFFECTS OF CONVECTIVE FLUID MOTION UPON OXIDE CRYSTAL GROWTH IN HIGH TEMPERATURE SOLUTION

For understanding of the influence of convective flow on crystal growth, space high temperature in situ observation instrument (SHITISOI) is dedicated to visualize and record the whole growth process of oxide crystals in high temperature up to 1000°C. Model experiments using transparent liquids such as KNbO3 and a mix ture of Li2B4O7+KNbO3 were chosen to investigate effects on ground and in space.On the earth, an investigation of growth kinetics of KNbO3 crystal related to two different states of convection: diffusive-advective flow and diffusive-convective flow,has been performed. The per unit length of a step e is calculated from the exper imental data for two different states of convection. Analyses of these data show the effect of buoyancy convection is to enhance the sharpness of the interface. The growth of KNbO3 crystals from solution of KNbO3+Li2B4O7 was investigated in space. The streamlines of the steady thermocapillary convection in Li2B4O7 solvent was observed. Due to thermocapillary convection, KNbO3 crystal grains grew and filled the whole solution homogeneously. Earth-based quenching experiments are de signed in order to study polyhedral instability of KNbO3 crystal, which is controlled by diffusion mechanism limitation. In all cases, when the crystal was nucleated near air/solution surface, it lost its polyhedral stability and varied from polyhedrons to dedrites. The thickness of diffusion mechanism limitation layer is about 60μm.     

Marangoni convection when the surface tension increases with the temperature in normal and low gravity conditions

Thermocapillary convection has been studied in n-heptanol aqueous solutions whose surface tension is increasing with the temperature.

The fluid was confined in a parallelipipedic enclosure and a thermal gradient was imposed parallel to the free liquid/gas interface. The motions induced by the thermocapillary forces have been studied under low gravity conditions during Texus sounding rocket experiments and during the Spacelab D1 mission.

The combined thermocapillary and buoyancy convection have been extensively studied during ground based experiments.

The influence of the aspect ratio and of the alcohol concentration were investigated. Detailed velocity field was determined at the steady state under normal gravity conditions by Laser Doppler Anemometry (LDA).

The aspect ratio was of 0.3 and the temperatures imposed to the lateral sides of the cell were respectively 45°C and 60°C. With these experimental conditions, two superposed contrarotative cells were observed with rapid motions in the surface from the cold to the hot side. This convective pattern was also observed during the normal and high gravity periods of parabolic flights but during the low gravity period of the parabola the motions stopped everywhere in the fluid even in the surface and reappeared immediatly at the pull out of the parabola.     

Nusselt number for convection driven by tidal and radiogenic heating

The ‘traditional’ formulas giving the Nusselt number Nu as a function of Rayleigh number Ra cannot be used for low and moderate values of Ra. Moreover, the recent progress in 3D numerical modeling makes possible to determine the Nusselt number Nu as a function of Rayleigh number Ra for convection driven by radiogenic and tidal heating. We found that for low and moderate Ra: Nu(Ra) = ε(Ra+ξ)^λ where λ depends on rheology and boundary conditions, ε depends only on the mode of heating, and ξ = ε^−1/λ · Nu(Ra) makes possible to develop a parameterized theory of convection in medium size icy satellites. We also indicate some differences between convection driven by tidal and radiogenic heating and convection driven exclusively by radiogenic heating.     

Numerical study of forced convective heat transfer around a spherical aerostat

Forced convective heat transfer is one of the major factors that dominate the thermal behaviors of aerostats. Due to the large physical size, the convection around an aerostat has high Reynolds numbers. The existing forced convective heat transfer correlations are limited to the Reynolds number lower than 10⁵, which are not appropriate for aerostat applications. Therefore, it is necessary to obtain a convective heat transfer correlation applicable to spherical aerostats at high Reynolds numbers. In this paper, steady convective heat transfer from an isothermal spherical aerostat is numerically investigated. The numerical simulation is carried out by commercial computational fluid dynamic software with the Reynolds number from 20 to 10⁸. The average Nusselt numbers are obtained and compared with those of available in literature. Based on regression and optimization with software, a new piecewise correlation of Nusselt number is proposed. The verification shows that the new correlation is reliable.     

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

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