Étude expérimentale des transferts de chaleur fluide-paroi induits par convection thermocapillaire : influence du nombre de PrandtlExperimental study of fluid-wall heat transfer induced by thermocapillary convection: influence of the Prandtl number

This article concerns heat transfer due to thermocapillary convection induced by an air bubble introduced under a heated wall into a silicone oil layer. An experimental study is carried out under normal gravity conditions with direct measurements by means of a heat fluxmeter. Effects of the bubble size and the temperature gradient are systematically investigated for two values of viscosity. The nature of the thermocapillary convection state (stationary or oscillatory) is determined by shadowgraphy. Two regimes of heat transfer are made evident. They depend on the Prandtl number but are independent of the state of the flow. To cite this article: C. Reynard et al., C. R. Mecanique 331 (2003).     

Convective self-oscillations near an air-bubble surface in a horizontal rectangular channel

The concentration convection in an isothermal fluid near an air bubble clamped between the vertical walls of a horizontal channel with a rectangular cross-section is studied experimentally and numerically. The channel is filled with an aqueous solution of a surfactant with a nonuniform concentration. As a result of the competition between the gravitational convection in the cavity volume and the Marangoni convection near the bubble surface, an oscillation flow regime is established. This regime is observed experimentally over several hours. In the numerical experiment, the oscillations are obtained in the presence of an initial horizontal surfactant concentration gradient. Against the background of gravitational convection, short bursts of Marangoni convection with ten times greater intensity are observed. The convective flow patterns and the oscillation periods obtained experimentally and numerically are in fairly good agreement.     

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)     

Oscillatory Marangoni convection around the air bubble in a vertical surfactant stratificationConvection de Marangoni oscillante autour d'une bulle d'air dans une stratification verticale surfactante

The solutocapillary Marangoni convection around a gas bubble in the inhomogeneous binary mixture of miscible fluids with a vertical surfactant concentration gradient was studied experimentally. A new phenomenon, the oscillatory instability of the surfactant mass transfer, near the bubble boundary, was detected and investigated. The interpretation of this effect as an interaction between the surfactant adsorption at the bubble free surface and solutocapillary and buoyancy convective mechanisms is proposed. The experimental data on oscillation period in relation to bubble dimensions, time, liquid layer thickness, physico-chemical fluid parameters and concentration gradients are presented and discussed. To cite this article: K. Kostarev et al., C. R. Mecanique 332 (2004).     

Stability of a numerical algorithm for gas bubble modelling

A free-surface-tracking algorithm based on the SOLA-VOF method is analysed for numerical stability when modelling gas bubble evolution in a fluid. It is shown that an instability can arise from the fact that the bubble pressure varies with its volume. A time step stability criterion is introduced which is a function of the natural oscillation period but does not depend on the mesh size. This dependence suggests that the instability is likely to arise in the case of a general motion of a bubble, especially if break-up occurs. The effect is shown using linear Fourier analysis of the discretized equation for radial bubble oscillation and demonstrated numerically using a CFD code FLOW-3D. One- and three-dimensional situations are considered: a bubble in a fluid bounded by two concentric surfaces and a bubble floating in a fluid chamber with and without gravity. In cases where no analytical solution is available, a numerical method for the stability time step limit calculation is suggested based on finding the natural oscillation frequency. The nature of the instability suggests that it can be a feature of any numerical algorithm which models transient fluid flow with a free surface.     

Dissolved gas effects on thermocapillary convection during boiling in reduced gravity environments

The mechanisms by which thermocapillary convection arises during boiling of nominally pure fluids in low-g environments are currently not known. It has recently been suggested that small amounts of dissolved gas within the bulk liquid can accumulate within the vapor bubble, forming localized concentration gradients that results in a temperature gradient to form along the liquid–vapor interface that drives thermocapillary convection. This hypothesis was tested by boiling > 99.3% pure n-perfluorohexane with and without noncondensible gas in a low-g environment using a 7.0 × 7.0 mm² microheater array to measure time and space resolved heat transfer at various wall superheats. The thermocapillary convection around the primary bubble that formed in the gassy fluid was found to be much weaker than in the degassed fluid, and the primary bubble diameter was much larger in the gassy fluid due to the accumulation of noncondensible gas within the bubble. The results suggest that the accumulation of noncondensible gas in the bubble can result in temperature variations along the interface but due to the increased vapor/gas bubble size, the driving thermocapillary temperature gradient along the interface is significantly reduced and result in much weaker thermocapillary flow. The highest CHF values in a reduced gravity environment (19 W/cm²) occurred when the fluid was highly subcooled and degassed.     

水中爆炸气泡脉动现象的实验研究

在实验水箱中采用高速摄影技术获得了炸药水中爆炸气泡脉动过程图像。实验结果表明:短长径比（1.05）PETN炸药水中爆炸产生的气泡形状近似球形,气泡脉动周期和半径随炸药质量的增加而增大。在膨胀阶段,气泡中心位置保持不变,随着体积的不断缩小,气泡上浮越来越明显。受重力影响,气泡下表面收缩速度高于上表面,在收缩至最小时下表面向上冲顶,气泡溃灭形成水射流。     

Numerical and experimental study of the heat transfer and fluid flow by thermocapillary convection around gas bubbles

 At liquid–gas or liquid–liquid interfaces thermocapillary or Marangoni convection develops in the presence of a temperature or concentration gradient along the interface. This convection was not paid much attention up to now, because under terrestrial conditions it is superimposed by the strong buoyancy convection. In a microgravity environment, however, it is the remaining mode of natural convection. During boiling in microgravity it was observed at subcooled conditions. Therefore the question arises about its contribution to the heat transfer. Thus the thermocapillary convection was intensively studied at single gas bubbles in various liquids both experimentally and numerically. Inside a temperature gradient chamber, the overall heat transfer around single bubbles of different volume was measured with calorimetry and the liquid flow with PIV and LDV. In parallel to the experiment, a 2-dimensional mathematical model was worked out and the coupled heat transfer and fluid flow was simulated with a CV-FEM method both under earth gravity level and under microgravity. The results are described in terms of the dimensionless Nusselt-, Peclet-, Marangoni-, Bond- and Prandtl-number. Received on 23 August 1999     

Note on buble motion in a rotating liquid under residual gravity

The motion of a small gas bubble, presumed to retain its geometrical shape and contained in a rotating liquid, has been investigated. The fluid system liquid-gas is subject to the influence of a reduced gravitational field. It is demonstrated that under certain conditions (spin axis and direction of gravity are perpendicular to each other) the bubble travels on a circular path about the axis of rotation, as seen from an observer moving with the bulk of the liquid.     

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

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

Marangoni效应下气泡上升过程的FTM模拟

采用界面跟踪法FTM(Front-Tracking Method),并结合CSF(continuum surface force)模型,研究了在垂直方向上温度分布不均匀的对称流场中由Marangoni效应引起的气泡上升运动问题。模拟了在不同的M a数和Pr数下单气泡及同轴双气泡的运动。研究结果表明,在不同的M a数下气泡的运动差异较大,M a数越大,气泡运动至稳态时的速度越大,且气泡运动的最大速度值与M a数呈正相关关系;增大Pr数所造成的粘度增大或热扩散率减小将削弱气泡的迁移运动;Marangoni对流中双气泡的局部运动证实了温度梯度和气泡运动速度紧密相关。     

Steady motion of a rigid disk of finite thickness on a horizontal plane

The article discusses the steady motion of a rigid disk of finite thickness rolling on its edge on a horizontal plane under the influence of gravity. The governing equations are presented and two cases allowing for a steady-state solution are considered: rolling on consistently rough ground and rolling on perfectly smooth ground. The conditions of steady motion are derived for both kinds of ground and it is shown that the possible steady motion of a disk is either on a straight line or in a circle. Oscillations about steady state are discussed and conditions for stable motion established. The bifurcations of steady motions on a smooth surface are also considered.     

Three-dimensional convective interactions in an inclined longitudinal air layer

The results of a three-dimensional numerical simulation of thermal-gravitational convection for an arbitrary direction of the gravity force are presented. The calculations were carried out for a rectangular domain with the solid-wall aspect ratio 4:0.5:1. The cavity filled with air was heated from below. The convection regime was laminar with Rayleigh numbers Ra = 10³–10⁵. The cavity was gradually turned about the short side in different directions in the range 0–90° to study how the initial conditions affect the convection. The average heat transfer characteristics and the spatial flow patterns are presented. A comparison with two-dimensional calculations is given. The conditions of the occurrence of a hysteresis in the steady state of convection are found, a regime map is constructed, and the specific features of the interaction of heat fluxes are described.     

Experimental study of the boundary layer over an airfoil in plunging motion

This is an experimental study on the boundary layer over an airfoil under steady and unsteady conditions.It specifically deals with the effect of plunging oscillation on the laminar/turbulent characteristics of the boundary layer.The wind tunnel measurements involved surfacemounted hot-film sensors and boundary-layer rake.The experiments were conducted at Reynolds numbers of 0.42×10 6 to 0.84 × 10 6 and the reduced frequency was varied from 0.01 to 0.11.The results of the quasi-wall-shear stress as well as the boundary layer velocity profiles provided important information about the state of the boundary layer over the suction surface of the airfoil in both static and dynamic cases.For the static tests,boundary layer transition occurred through a laminar separation bubble.By increasing the angle of attack,disturbances and the transition location moved toward the leading edge.For the dynamic tests,earlier transition occurred with increasing rather than decreasing effective angle of attack.The mean angle of attack and the oscillating parameters significantly affected the state of the boundary layer.By increasing the reduced frequency,the boundary layer transition was promoted to the upstroke portion of the equivalent angle of attack,but the quasi skin friction coefficient was decreased.     

电场作用下气泡上升行为特性的数值计算研究

利用电场控制气泡形态及运动,强化气液相间传热传质是电流体动力学的重要研究内容之一. 然而目前多数研究集中在非电场下的气泡动力学上,对于电场下的气泡行为特性及电场的作用机制仍需开展深入研究. 本研究对电场作用下单个气泡在流体中上升过程的动力学行为进行了数值模拟研究. 在建立二维模型的基础上求解电场方程与Navier-Stokes方程,并采用水平集方法捕捉了上升气泡的位置及形状. 模拟结果的准确性与有效性通过与前人实验和数值结果进行对比得到了验证. 通过改变雷诺数、邦德数和电邦德数等不同参数研究了电场下液体黏度、表面张力和电场力对气泡运动变形的影响. 计算结果表明,电场对气泡的动态特性有显著影响. 非电场情况下液体黏度和表面张力较大时气泡基本维持球状,反之气泡发生变形并逐步达到稳定状态. 此外,电场作用使气泡在初始上升阶段发生剧烈形变,随着不断上升,气泡形变程度不断减小,且气泡的上升速度和长径比均出现振荡. 垂直电场使气泡的上升速度有较大的提高,且随着电邦德数的增大,难以达到相对稳定的状态.      

热毛细对流温度场ESPI检测研究

在地面上用二维液体盒研究了热壁下气泡周围液体中的热毛细对流现象,并应用电子散斑干涉技术（ESP）对热毛细对流温度场进行了实时检测研究,给出了部分典型的实验结果。     

Conditions of Rayleigh-Bénard convection onset and heat transfer in a near-critical medium

Thermal gravity convection in a horizontal layer of compressible perfect gas heated from below and a van der Waals gas near the critical state is investigated. The characteristics of the isentropic equilibrium of a compressible medium with a van der Waals equation of state are considered. The known conditions of convection onset in the perfect and van der Waals gases are checked on the basis of a solution of the complete and linearized equations. The restrictions imposed in deriving the known formulas for the adiabatic temperature gradients used in the conditions of absence and onset of convection are discussed. The characteristics of the convective heat transfer are examined, including the causes of the heat-transfer deterioration in the near-critical medium above the hydrostatic equilibrium threshold.     

High-fidelity numerical simulation of the flow field around a NACA-0012 aerofoil from the laminar separation bubble to a full stall

In the present work, large eddy simulations of the flow field around a NACA-0012 aerofoil near stall conditions are performed at a Reynolds number of 5 × 10⁴, Mach number of 0.4, and at various angles of attack. The results show the following: at relatively low angles of attack, the bubble is present and intact; at moderate angles of attack, the laminar separation bubble bursts and generates a global low-frequency flow oscillation; and at relatively high angles of attack, the laminar separation bubble becomes an open bubble that leads the aerofoil into a full stall. Time histories of the aerodynamic coefficients showed that the low-frequency oscillation phenomenon and its associated physics are indeed captured in the simulations. The aerodynamic coefficients compared to previous and recent experimental data with acceptable accuracy. Spectral analysis identified a dominant low-frequency mode featuring the periodic separation and reattachment of the flow field. At angles of attack α ≤ 9.3°, the low-frequency mode featured bubble shedding rather than bubble bursting and reformation. The underlying mechanism behind the quasi-periodic self-sustained low-frequency flow oscillation is discussed in detail.     

Nonlinear dynamics of hydrostatic internal gravity waves

Stratified hydrostatic fluids have linear internal gravity waves with different phase speeds and vertical profiles. Here a simplified set of partial differential equations (PDE) is derived to represent the nonlinear dynamics of waves with different vertical profiles. The equations are derived by projecting the full nonlinear equations onto the vertical modes of two gravity waves, and the resulting equations are thus referred to here as the two-mode shallow water equations (2MSWE). A key aspect of the nonlinearities of the 2MSWE is that they allow for interactions between a background wind shear and propagating waves. This is important in the tropical atmosphere where horizontally propagating gravity waves interact together with wind shear and have source terms due to convection. It is shown here that the 2MSWE have nonlinear internal bore solutions, and the behavior of the nonlinear waves is investigated for different background wind shears. When a background shear is included, there is an asymmetry between the east- and westward propagating waves. This could be an important effect for the large-scale organization of tropical convection, since the convection is often not isotropic but organized on large scales by waves. An idealized illustration of this asymmetry is given for a background shear from the westerly wind burst phase of the Madden–Julian oscillation; the potential for organized convection is increased to the west of the existing convection by the propagating nonlinear gravity waves, which agrees qualitatively with actual observations. The ideas here should be useful for other physical applications as well. Moreover, the 2MSWE have several interesting mathematical properties: they are a system of nonconservative PDE with a conserved energy, they are conditionally hyperbolic, and they are neither genuinely nonlinear nor linearly degenerate over all of state space. Theory and numerics are developed to illustrate these features, and these features are important in designing the numerical scheme. A numerical method is designed with simplicity and minimal computational cost as the main design principles. Numerical tests demonstrate that no catastrophic effects are introduced when hyperbolicity is lost, and the scheme can represent propagating discontinuities without introducing spurious oscillations.      

Oscillations of a gas bubble in a non-Newtonian liquid under the action of an acoustic field

The evolution of the radius of a spherical cavitation bubble in an incompressible non-Newtonian liquid under the action of an external acoustic field is investigated. Non-Newtonian liquids having relaxation properties and also pseudoplastic and dilatant liquids with powerlaw equation of state are studied. The equations for the oscillation of the gas bubble are derived, the stability of its radial oscillation and its spherical form are investigated, and formulas are given for the characteristic frequency of oscillations of the cavitation hollow in a relaxing liquid. The equations are integrated numerically. It is shown that in a relaxing non-Newtonian liquid the viscosity may lead to the instability of the radial oscillations and the spherical form of the bubble. The results obtained here are compared with the behavior of a gas bubble in a Newtonian liquid.