大尺寸液桥热毛细对流失稳性地面实验研究

开展大尺寸液桥浮力-热毛细对流地面实验, 探究流场转捩的临界条件及临界状态附近的流动情况. 通过粒子图像测速方法(PIV) 获得流体速度场, 研究液桥内部定常和转捩后的流场结构以及流体运动规律;并用红外热像仪测量液桥自由面温度分布, 研究流体流动的时空演化和温度振荡. 实验发现大尺寸半浮区液桥浮力-热毛细对流临界值与几何参数有关, 在大普朗特(Prandtl) 数情况下, 流场存在由稳定态向不稳定态再到混沌的转捩过程, 在临界马兰哥尼(Marangoni) 数附近, 流场内会出现行波现象, 流动模式也会随高径比的变化而发生变化;当继续增大马兰哥尼数, 流动会进入混沌状态.      

矩形液池热毛细对流转捩途径研究

主要研究矩形液池热毛细对流的分岔转捩. 通过测量流体内部温度振荡情况, 详细研究了热毛细对流的转捩过程和转捩途径. 实验发现, 矩形液池热毛细对流的转捩过程依次经历了定常、规则振荡、不规则振荡的阶段. 对于不同普朗特数的硅油在不同长高比情况下, 通向混沌的途径不同. 在转捩过程中, 随着温差的增加, 普朗特数在16 (1cSt) 以下和普朗特数为25 (1.5cSt)、长高比为26 的硅油热毛细对流主要以准周期分岔的转捩方式为主;而普朗特数为25 以上的则以倍周期分岔的转捩方式为主;两种分岔有时还会伴随有切分岔形式的出现.实验中还观察到了表面波动和对流涡胞振荡等现象.      

大普朗特数大液桥浮力-热毛细对流地面实验

通过地面实验研究大尺寸液桥的浮力-热毛细对流. 实验采用2cst硅油(Pr=28.571),研究了不同高径比(A=l/d)和体积比的液桥起振,分析了温度振荡频率及相位变化,探讨了热流体波的问题. 实验液桥的桥柱直径为20mm,由于受重力的限制,建立了3～4.25mm范围内的矮桥. 通过伸入液桥内部不同位置的热电偶的温度信号,发现流场是同时起振的,不同的桥高和体积比有不同的振荡模式,并且随着温差的增加,频率近似以线性增加,各点的振荡相位是一个连续性变化的过程. 不同高径比的液桥转捩到混沌的途径是不一样的.      

双自由面溶质?热毛细液层的不稳定性

溶质?热毛细对流是流体界面的浓度和温度分布不均导致的表面张力梯度驱动的流动, 它主要存在于空间微重力环境、小尺度流动等表面张力占主导的情况中, 例如晶体生长、微流控、合金浇筑凝固、有机薄液膜生长等. 对其流动进行稳定性分析具有重要意义. 本文采用线性稳定性理论研究了双自由面溶质?热毛细液层对流的不稳定性, 得到了两种负毛细力比(η)下的临界Marangoni数与Prandtl数(Pr)的函数关系, 并分析了临界模态的流场和能量机制. 研究发现: 溶质?热毛细对流和纯热毛细对流的临界模态有较大的差别, 前者是同向流向波、逆向流向波、展向稳态模态和逆向斜波, 后者是逆向斜波和逆向流向波. 在Pr较大时, Pr增加会降低流动稳定性; 在其他参数下, Pr增加会增强流动稳定性. 在中低Pr, 溶质毛细力使流动更加不稳定; 在大Pr时, 溶质毛细力的出现可能使流动更加稳定; 在其他参数下, 溶质毛细力会减弱流动稳定性. 流动稳定性不随η单调变化. 在多数情况下, 扰动浓度场与扰动温度场都是相似的. 能量分析表明: 扰动动能的主要能量来源是表面张力做功, 但其中溶质毛细力和热毛细力做功的正负性与参数有关.      

半浮区液桥热毛细振荡流

采用非定常、三维直接数值模拟方法研究大Pr数半浮区液桥热毛细对流从定常流向振荡流的过渡过程．文中详细描述了热毛细振荡流的起振和振荡特征,给出了液桥横截面上振荡流的流场和温度分布．在地面引力场条件下计算的结果与地面实验的结果进行比较,得出液桥水平截面上的流场和温度分布图样以一定的速度旋转,自由表面固定点处流体的环向流速正、负交替变化的一致结论．     

Bingham流体双自由面热毛细液层的稳定性分析

热毛细对流是流体界面温度分布不均导致的表面张力梯度驱动的流动. 它主要存在于空间等微重力环境或小尺度流动等表面张力占主导的情况中. 在很多工业领域, 如晶体生长、聚合物加工、喷墨打印、微流控, 产品质量都与热毛细对流密切相关. 空间3D打印是太空制造的重要技术, 可以支持空间站的在轨长期有人照料的运行和维护, 实现按需制造. 本文以聚合物流体的空间3D打印为应用背景, 采用线性稳定性理论研究了Bingham流体双自由面热毛细液层的稳定性, 得到了在不同Bingham数(B)下的临界Marangoni数(Ma_c)与Prandtl数(Pr)的函数关系,分析了临界模态的流场和能量机制. 研究发现: 该流动的临界模态包括流向波和斜波模态, 与B, Bi和两界面垂直方向上的温差(Q)相关. B和Bi的增加会增强热毛细对流的稳定性. 当Q = 0时, 扰动温度分布分成对称和反对称两种情况. 当Q > 0时, Pr的增加会减弱流动稳定性. 在小Pr情况下, 扰动温度分布在整个流场, 在大Pr情况下, 扰动温度在栓塞区为零. 能量分析表明: 扰动动能的主要能量来源是表面张力做功, 但小Pr数下基本流也有一定贡献.      

大高宽比矩形腔体侧向正弦加热下对流扰动的成长和传热特性

为了研究矩形腔体侧向正弦周期加热条件下对流扰动的成长和传热特性,本研究对流体力学方程组进行了数值模拟。结果表明：随着格拉晓夫数Gr的增加,对流扰动的最大振幅A_max的成长率变大,线性成长阶段的时间变短。对于高宽比A=10,腔体宽度d=2 cm的腔体,普朗特数Pr=6.949的流体,对流扰动的成长率γ_m随着格拉晓夫数Gr变化的经验式为γ_m=9×10^-8 Gr^{1.234 3};对于高宽比A=10、腔体宽度d=6 cm的腔体,普朗特数Pr=0.703的流体,成长率γ_m随着格拉晓夫数Gr的变化关系式为γ_m=8×10^-4 Gr^0.52。当Pr=0.027 2时,热壁面平均努塞尔数■和格拉晓夫数Gr的关联式为■;当Pr=6.949时,■。随着格拉晓夫数Gr数增加,右壁面的传热能力增强。     

液封液桥内热毛细对流的数值模拟

建立了液封液桥(不相溶混的双层同轴液柱)内热毛细对流的物理数学模型，采用涡量-流 函数法、利用有限差分格式对微重力条件下液封液桥内热毛细对流进行了数值模拟，得到了 双层液柱主流区的温度场和流场，证实了液封技术能削弱液桥主流区的热毛细对流，从而提 高浮区晶体生长质量，找出了液封厚度以及内、外层流体物性参数比对液桥内热毛细对流的 影响规律.     

实践十号卫星项目——热毛细对流振荡特征的地面研究

上部开口环形液池在水平径向温度梯度作用下会出现内部温度和自由面的振荡,本文研究了二者发生的临界条件.环形液池内柱加热外壁制冷,以0.5℃/min的速率线性升温得到水平径向温差,T型热电偶测量液层内部单点温度,高精度激光位移传感器测量液层自由面某点形变.随两端温差增加,当超过某一临界温度△T_(cr)时,开始出现振荡.实验结果表明,对同一种硅油,两种振荡的临界条件随液层厚度具有相同的变化趋势.对不同普朗特(Prandtl)数(Pr)的硅油,振荡临界条件临界马兰哥尼数(Ma_(cr))随着邦德数(Bo)的增加而变大.本文工作是中国科学院科学先导专项SJ-10返回式科学实验卫星项目-热毛细对流表面波空间实验研究的地面研究结果,该工作为空间实验提供前期的基础科学研究数据和实验保障.     

半浮区热毛细对流的不稳定性与转捩

近二十年来,微重力流体开展了半浮区液桥热毛细对流的不稳定性与转捩的研究．文中给出了热毛细振荡对流发生的临界参数,分析了液桥几何位形（尺度比,体积比）、物理参数及传热参数对临界Ｍａｘａｎｇｏｎｉ的影响．报导了有关的地面模拟实验,微重力实验以及本问题的线性稳定性分析、能量分析和数值模拟结果,并介绍了定常轴对称热毛细对流通过非定常振荡热毛细对流到湍流的转捩过程和三种热毛细振荡对流的产生机理．     

Three-dimensional numerical simulation of thermocapillary instabilities in floating zones

Instability of thermocapillary convection in liquid bridges of low Prandtl number fluids is investigated by direct three-dimensional and time-dependent simulation of the problem. The field equations are numerically solved explicit in time and with finite difference methods in a staggered cylindrical grid. The numerical results are analyzed and interpreted in the general context of the bifurcation's theory. According to recent stability analyses the computations show that for semiconductor melts the first bifurcation is characterized by the loss of spatial symmetry rather than by the onset of oscillatory flow. When the basic axisymmetric flow field becomes unstable, after a short transient, a three-dimensional supercritical steady state is obtained. It is shown that the flow field organization, depending on the critical wave number, is related to the geometrical aspect ratio of the liquid bridge and that lower is the aspect ratio, higher is the critical wave number and more complex the thermofluid-dynamic field structure.     

Effect of vertical heat transfer on thermocapillary convection in an open shallow rectangular cavity

In order to understand the effect of the vertical heat transfer on thermocapillary convection characteristics in a differentially heated open shallow rectangular cavity, a series of two- and three-dimensional numerical simulations were carried out by means of the finite volume method. The cavity was filled with the 1cSt silicone oil (Prandtl number Pr = 13.9) and the aspect ratio ranged from 12 to 30. Results show that thermocapillary convection is stable at a small Marangoni number. With the increase of the heat flux on the bottom surface, thermocapillary convection transits to the asymmetrical bi-cellular pattern with the opposite rotation direction. The roll near the hot wall shrinks as the Marangoni number increases. At a large Marangoni number, numerical simulations predict two types of the oscillatory thermocapillary flow. One is the hydrothermal wave, which is dominant only in a thin cavity. The other appears in a deeper cavity and is characterized by oscillating multi-cellular flow. The critical Marangoni number for the onset of the oscillatory flow increases first and then decreases with the increase of the vertical heat flux. The three-dimensional numerical simulation can predict the propagating direction of the hydrothermal wave. The velocity and temperature fields obtained by three-dimensional simulation in the meridian plane are very close to those obtained by two-dimensional simulation.     

Influence of rotating magnetic field strength on three-dimensional thermocapillary flow in a floating half-zone model

The effects of rotating magnetic field (RMF) on the three-dimensional thermocapillary flow of semiconductor melt (Pr?=?0.01) in a floating half-zone model under microgravity are investigated numerically by the finite volume method. The results indicate that the thermocapillary flow without magnetic field is a steady three-dimensional convection for Ma?=?40 in a floating half-zone model with As?=?1, and the convection evolves to an oscillatory three-dimensional flow by applying 1–6?mT RMF with 50?Hz rotating frequency. Based on the fast Fourier transform spectrum, the convection is confirmed to be a periodically oscillating flow, the oscillatory main frequency, 1.59?×?10^?3?Hz for 1?mT RMF and 5.84?×?10^?2?Hz for 6?mT RMF, increases with the magnetic strength. However, with increasing the magnetic field strength up to 7?mT, the three-dimensional thermocapillary flow is effectively controlled and the convection turns into a steady axisymmetrical one.     

Onset of unsteady axi-symmetric laminar natural convection in a vertical cylindrical enclosure heated at the wall

In the present study laminar transition to oscillatory convection of fluids having different Prandtl numbers in a laterally heated vertical cylindrical enclosure for different aspect ratios (melt height to crucible radius) of 2–4 is investigated numerically for 0.01 ≤ Pr ≤ 10. Numerical solution to two-dimensional axisymmetric transient Navier Stokes equations and energy equation were solved by finite volume method using SIMPLE algorithm. Numerical results illustrate that there exists a critical Rayleigh number for each Prandtl number beyond which sustained laminar oscillatory flow sets in. The oscillatory regime was characterised by the oscillation of the average kinetic energy and average thermal energy of the melt. For a given aspect ratio, critical Rayleigh number increases with Pr upto 1 and then flattens. It was observed that for low Prandtl number fluids, Pr < 1.0, critical Rayleigh number is found to increase with increase in aspect ratio while for high Prandtl number fluids, Pr ≥ 1.0, it is found to decrease with increase in aspect ratio. The influence of aspect ratio on the transient behaviour of the melt volume below and above the critical Rayleigh number was studied.     

Two bifurcation transition processes in floating half zone convection of larger Prandtl number fluid

Processes of the onset oscillation in the thermocapillary convection under the Earth's gravity are investigated by the numerical simulation and experiments in a floating half zone of large Prandtl number with different volume ratio. Both computational and experimental results show that the steady and axisymmetric convection turns to the oscillatory convection ofm=1 for the slender liquid bridge, and to the oscillatory convection before a steady and 3D asymmetric state for the case of a fat liquid bridge. It implies that, there are two critical Marangoni numbers related, respectively, to these two bifurcation transitions for the fat liquid bridge. The computational results agree with the results of ground-based experiments. The project supported by the National Natural Science Foundation of China (19789201) and the Ministry of Science and Technology of China (95-yu-34)     

Effects of a.c. Electric Field and Rotation on Bénard–Marangoni Convection

The coupled buoyancy and thermocapillary instability, the Bénard–Marangoniproblem, in an electrically conducting fluid layer whose upper surface is deformed and subject to a temperature gradient is studied. Both influences of an a.c. electric field and rotation are investigated. Special attention is directed at the occurrence of convection both in the form of stationary motion and oscillatory convection. The linear stability problem is solved for different values of the relevant dimensionless numbers, namely the a.c. electric Rayleigh number, the Taylor, Rayleigh, Biot, Crispation and Prandtl numbers. For steady convection, it is found that by increasing the angular velocity, one reinforces the stability of the fluid layer whatever the values of the surface deformation and the applied a.c. electric field. We have also determined the regions of oscillatory instability and discussed the competition between both stationary and oscillatory convections. This revised version was published online in July 2006 with corrections to the Cover Date.     

INSTABILITY OF COUPLED THERMO-SOLUTE CAPILLARY CONVECTION IN LIQUID BRIDGE AND CONTROL BY ROTATING MAGNETIC FIELD

浮区法因具有无坩埚接触污染的生长优点而成为生长高完整性和高均匀性单晶材料的重要技术.但熔体中存在的毛细对流会给浮区法晶体生长带来极大挑战,这是由于对流的不稳定会导致晶体微观瑕疵的产生和宏观条纹等缺陷的形成.为了提高浮区法生长单晶材料的品质,研究浮区法晶体生长中毛细对流特性及如何控制其不稳定性显得尤为重要.本文采用数值模拟的方法对半浮区液桥内Si_xGe_1-x体系中存在的热质毛细对流展开研究并施加旋转磁场对其进行控制.结果表明：纯溶质毛细对流表现为二维轴对称模式,温度场主要由热扩散作用决定,而浓度场则由对流和溶质扩散共同支配;纯热毛细对流呈现三维稳态非轴对称流动,浓度分布与熔体内热毛细对流的流向密切相关,等温线在对流较大的区域发生弯曲;耦合溶质与热毛细对流则为三维周期性旋转振荡流.施加旋转磁场后,熔体周向速度沿径向向外增大,熔体内浓度场和流场均呈现二维轴对称分布.