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.     

Experimental study on thermocapillary convection

In the present paper, the experimental studies on thermocapillary convection are reviewed. The author‘s interest is mainly focused on the onset of oscillatory thermocapillary convection,the features of oscillatory flow pattern, and the critical Marangoni number related with temperature and free surface oscillation. The coordinated measurement in a microgravity environment of a drops haft is also addressed.     

浮区热毛细对流

概述了浮区中平行于自由面的表面张力梯度驱动热毛细对流领域的研究. 研究兴趣集中于振荡热毛细对流的起振, 或者说从定常流动到振荡流动的转捩. 起振依赖于一系列的临界参数, 临界关系可以表示为这些临界参数的复杂函数. 实验结果表明, 振荡流中速度的变化和平均流动的速度有相同的量级, 而其它量的变化, 比如温度和自由面半径的波动, 相比于它们的平均量而言则要小得多. 因此, 起振应是流体中动力学过程的结果, 该问题是强非线性的. 在过去几十年中, 一些理论模型被引入来研究这个问题, 使用的方法包括理论分析方法、 线性不稳定性分析方法、 能量稳定性分析方法以及非定常的三维直接数值模拟. 其中直接数值模拟被认为是对强非线性过程进行深入分析的最适合方法, 通常能得到和实验较符合的结果. 从振荡热毛细对流向湍流的转捩提供了一个研究混沌行为的新系统, 开创了一个非线性科学的新前沿, 是一个集中了大量近期工作的研究热点. 该文对浮区热毛细对流作了一个回顾, 包括理论模型和分析, 以及实验研究.     

半浮区液桥热毛细振荡流

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

Thermocapillary instability of the equilibrium of a flat layer containing internal heat sources

In the absence of body forces, a factor which has a strong influence on the equilibrium stability of a nonuniformly heated liquid is the dependence of the coefficient of surface tension on the temperature and the thermocapillary effect generated by it. If the equilibrium temperature gradient is sufficiently great, then the presence of the thermocapillary forces on the free surface can lead to the occurrence of convective motion. The monotonie instability of the equilibrium of a flat layer was investigated in [1–3]. Analysis of nonmonotonic disturbances [4] showed that in the case of an undeformable free surface there is no oscillatory instability. In [5] it was found that oscillatory instability is possible if there is a nonlinear dependence of the coefficient of surface tension on the temperature. The present paper is devoted to numerical investigation of the equilibrium stability of a flat layer with respect to arbitrary disturbances. It is shown that for a deformable free boundary there appears an additional neutral curve, which corresponds to monotonie capillary instability. In addition, when the capillary convection mechanism is taken into account, there appears an oscillatory instability, which becomes the most dangerous in the region of small Prandtl and wave numbers.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 27–31, March–April, 1991.I thank V. K. Andreev for a helpful discussion of the work.     

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

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

Thermocapillary convection in two-layer systems in the presence of a surface active agent at the interface

Convective instability in a layered system due to the thermocapillary effect was investigated in [1–5]. In these studies it was shown that the perturbations responsible for equilibrium crisis may build up either monotonically or in an oscillatory fashion. In [6] the stabilizing effect of a surface active agent (SAA) on thermocapillary instability was established for a layer with a free surface. For layers of infinite thickness the effect of SAA on thermocapillary convection was studied in [7–9]. The present investigation is concerned with thermocapillary convection in a system of two layers of finite thickness in the presence of an SAA. Convection due to the lift force is not considered. It is established that the principal result of the action of the SAA is not the stabilizing effect on the monotonic mode but the appearance of a new type of oscillatory instability.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2 pp. 3–8, March–April, 1986.In conclusion the authors wish to thank E. M. Zhukhovitskii for discussind the results.     

Thermocapillary instability of a plane layer with a vertical temperature gradient

The oscillating disturbances in a plane layer with a temperature gradient are analyzed. It is shown that for heating from below taking the deformability of the free surface into account leads to the appearance of short-wave oscillatory instability, which becomes the most dangerous mode. Moreover, the interaction of the capillary and thermocapillary instability mechanisms results in the appearance of oscillating disturbances of a new type, which lead to equilibrium crisis at high Marangoni numbers. It is established that when the free boundary is heated, the onset of convection is possible only with respect to oscillatory disturbances.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.3, pp. 19–23, May–June, 1992.     

Axisymmetric thermocapillary convection in a slowly rotating annular cylindrical container

In a slowly rotating annular cylindrical container the free liquid surface (liquid-gas interface) is subjected to a temperature gradient in radial direction. The temperature dependent surface tension creates a shear stress on the interface which is transmitting a thermocapillary convection in the bulk of the liquid. For constant temperature T ₁ of the inner and T ₂ of the outer wall a steady Marangoni convection takes place, exhibiting a double vortex ring of equal directional flow. For time-oscillatory temperatures of the walls a time-dependent thermocapillary convection appears, which will create on the free liquid surface various wave patterns. They shall, depending on the forcing frequency of the temperature, exhibit resonance peaks. The velocity distribution and the response magnitude inside the container has been determined. Received on 3 September 1997     

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.     

Experimental and numerical investigation of nonlinear thermocapillary oscillations in an annular geometry

The present work investigates thermocapillary flow in a cylindrical configuration using large Prandtl number liquids. The flow is studied using coordinated simulations and experimental optical methods such as PIV and flow visualization. In this way, properties of the oscillatory state can be obtained in great detail. Considerable attention is given to the search for the parameters influencing the onset of the instability. It is found that the onset of oscillations can be correlated using a thermocapillary Reynolds number. The oscillations take the form of a standing wave close to the onset, which is replaced by a travelling wave for stronger forcing. The selection of azimuthal wave number of the oscillatory mode is determined from geometrical parameters, and resembles the wave number selection in vortex ring instabilities. Throughout we obtain good agreement between experiments and simulations using a mathematical model with an undeformed adiabatic free surface.     

Convection in a two-layer system due to the combined action of the Rayleigh and thermocapillary instability mechanisms

In a two-layer system loss of stability may be monotonic or oscillatory in character. Increasing oscillatory perturbations have been detected in the case of both Rayleigh [1, 2] and thermocapillary convection [3–5]; however, for many systems the minimum of the neutral curve corresponds to monotonic perturbations. In [5] an example was given of a system for which oscillatory instability is most dangerous when the thermogravitational and thermocapillary instability mechanisms are simultaneously operative. In this paper the occurrence of convection in a two-layer system due to the combined action of the Rayleigh (volume) and thermocapillary (surface) instability mechanisms is systematically investigated. It is shown that when the Rayleigh mechanism operates primarily in the upper layer of fluid, in the presence of a thermocapillary effect oscillatory instability may be the more dangerous. If thermogravitational convection is excited in the lower layer of fluid, the instability will be monotonic.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 166–170, January–February, 1987.     

Experimental Investigation of Thermocapillary Convection Induced by a Local Temperature Inhomogeneity near the Liquid Surface. 1. Solid Source of Heat

The structure and stability of a thermocapillary flow from a concentrated source of heat located near the free surface of the liquid filling a deep reservoir are experimentally studied. For a certain power of the heat source, oscillatory instability leading to formation of surface waves is observed. Possible mechanisms of the observed instability are discussed.     

Simulations of 2D and 3D thermocapillary flows by a least-squares finite element method

Numerical results for time-dependent 2D and 3D thermocapillary flows are presented in this work. The numerical algorithm is based on the Crank–Nicolson scheme for time integration, Newton's method for linearization, and a least-squares finite element method, together with a matrix-free Jacobi conjugate gradient technique. The main objective in this work is to demonstrate how the least-squares finite element method, together with an iterative procedure, deals with the capillary-traction boundary conditions at the free surface, which involves the coupling of velocity and temperature gradients. Mesh refinement studies were also carried out to validate the numerical results. © 1998 John Wiley & Sons, Ltd.     

Marangoni convection in an annular cylindrical tank

If the free liquid-gas interface of a liquid in an annular container is subjected to a temperature gradient the temperature dependent shear stress on the free liquid surface creates by viscous traction a thermocapillary convection in the bulk of the liquid. For constant temperature T ₁ of the inner and T ₂ of the outer wall a steady Marangoni convection takes place, while for time-oscillatory temperatures of the walls a time-dependent thermocapillary convection appears, which will create on the free liquid surface wave patterns. They shall, depending on the forcing frequency of the temperature, exhibit resonance peaks. The velocity distribution and the response magnitude inside the container has been determined. Received on 11 July 1997     

基于POD-Galerkin降维方法的热毛细对流分岔分析

作为流动与传热相互耦合的非线性过程, 热毛细对流有着复杂的转捩过程, 探究流场和温度场随参数变化而发生的分岔现象, 是热毛细对流研究的一个重要课题. 基于本征正交分解的POD-Galerkin降维方法可以通过提取特征模态, 构建低维模型, 实现流场的快速计算. 数值分岔方法可以通过求解含参数动力系统的分岔方程, 直接计算稳定解和分岔点. 探究了将直接数值模拟方法、POD-Galerkin降维方法、数值分岔方法的优势结合, 以提高热毛细对流转捩过程分析效率的可行性. 利用直接数值模拟得到的流场和温度场数据, 构建了不同体积比下, 二维有限长液层热毛细对流的POD-Galerkin低维模型, 在低维模型上采用数值积分及数值分岔方法计算了分岔点, 得到了低维方程的分岔图. 在一定参数范围内, 在低维模型上模拟热毛细对流, 对雷诺数和体积比进行参数外推, 通过与直接数值模拟的结果对比, 验证了低维模型的准确性与鲁棒性. 说明了低维方程可以定性反映原高维系统的流动特性, 而定量方面, 由低维模型和直接数值模拟计算得到的周期解频率的相对误差大约为5%. 验证了利用POD-Galerkin降维方法研究热毛细对流的可行性.      

A numerical study of the longitudinal thermoconvective rolls in a mixed convection flow in a horizontal channel with a free surface

This paper presents a numerical study of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The coupled Navier–Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects (Marangoni effect). When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable (Ri ⋍ O(1) and Bd = Ra/Ma ⋍ O(1)), the results show that the development of instabilities in the form of steady longitudinal convective rolls is similar to those encountered in the Poiseuille–Rayleigh–Bénard flow. The number and spatial distribution of these rolls along the channel depend on the flow conditions. The objective of this work is to study the influence of parameters, such as the Reynolds, Rayleigh and Biot numbers, on the flow patterns and heat transfer characteristics. The effects of variations in the surface tension with temperature gradients (Marangoni effect) are also considered.     

热毛细对流的液桥几何位形及浮力效应

本文讨论重力对不同高度、直径比液桥的热毛细对流的影响。当液桥高度、直径比增大时,液桥中的等流函数线呈双涡结构,这种流动图样并不必然与热毛细振荡流相联系。在地面热毛细对流实验中模拟空间微重力情况,液桥高度需小于1.5mm。在微重力环境中,液桥内的流场和温度分布介于地面相同参数液桥的上部加热和下部加热两种结果之间。因此,可以用地面实验结果估计空间液桥的对流和热输运情况。     

Extensional and shear rheometry of oriented triblock copolymers

 The effects of extensional flow orientation on the rheological properties of two poly(styrene)-poly(ethylene-co-butylene)-poly (styrene) (PS-PEB-PS) triblock copolymers containing either spherical or cylindrical PS microdomains were studied by oscillatory shear and oscillatory extensional experiments. Extensional measurements revealed that below the PS block glass transition temperature pre-oriented triblocks display highly anisotropic mechanical properties. For both polymers, the storage modulus E ′ is higher along the flow direction. Above the PS glass transition temperature the materials are no longer anisotropic and the same storage moduli are obtained along the flow direction and perpendicular to it. Above the PS glass transition temperature the rheological behaviour parallel and perpendicular to the flow direction was also probed in pre-oriented and non-oriented samples by oscillatory shear rheometry. At high frequencies, the mechanical response of the triblocks was found to be independent of the orientation for both copolymers while at low frequencies a strong effect of the flow orientation could be observed. For both polymers the value of the storage modulus was found to be lower along the flow direction that perpendicular to it. This was explained by the ability of PS blocks to relax more easily along the flow direction. Received: 10 September 1999/Accepted: 1 October 1999