环形液池内中等Pγ数流体的浮力-热毛细对流

为了了解水平温度梯度作用时环形液池内的浮力-热毛细对流特性,利用有限差分法进行了非稳态三维数值模拟,环形液池外壁被加热,半径为40 mm,内壁被冷却,半径为20 mm,液池深度为3～17mm,液池内流体为0.65cSt的硅油,其Pr数为6.7.模拟结果表明,当水平温度梯度较小时,流动为轴对称稳态流动,随着温度梯度的增加,流动将会失去其稳定性,在浅的液池内(d=3 mm),转化成三维振荡流动,在深的液池内(d≥6mm),转化成三维稳定流动;模拟计算的临界温差及表面温度分布图像与实验结果基本吻合.     

环形液池内中等Pr数流体的浮力-热毛细对流

为了了解水平温度梯度作用时环形液池内的浮力-热毛细对流特性，利用 有限差分法进行了非稳态三维数值模拟，环形液池外壁被加热，半径为40 mm，内壁被冷却， 半径为20 mm，液池深度为3~17 mm，液池内流体为0.65cSt的硅油，其Pr数 为6.7. 模拟结果表明，当水平温度梯度较小时，流动为轴对称稳态流动，随着温度 梯度的增加，流动将会失去其稳定性，在浅的液池内(d=3mm)，转化成三维振荡 流动，在深的液池内(d≥6mm)，转化成三维稳定流动；模拟计算的 临界温差及表面温度分布图像与实验结果基本吻合.     

浮力-热毛细对流表面位型研究

将栅线结构光源系统与图像处理技术相结合,形成一种实时诊断浮力-热毛细对流表面形貌的光学测量系统.研究两端带有温差的矩形池内薄层流体的表面变形.应用图像互相关方法对实验结果进行计算和分析,得到了流体表面变形定量的实验结果.实验结果表明了在浮力-热毛细对流的发展过程中,出现流体表面变形,温度梯度、表面张力、浮力、液固亲润问题、以及壁面温度决定了该表面变形大小和变形方向.     

矩形液池内热毛细对流的流动不稳定性

热毛细对流及其不稳定性是微重力流体科学研究的重要内容. 对该问题的研究不仅有利于人们对微重力环境下流体行为、对流不稳定性和湍流转捩过程等基础物理现象的进一步认识,而且也将促进晶体生长、薄膜制备等空间和地面高新技术的发展. 实验研究了矩形液池中浅液层在水平温度梯度作用下产生的热毛细对流及其稳定性. 实验中,成功地利用PIV (particle image velocimetry) 技术对1mm²/s 硅油液层内的浮力热毛细对流流场结构进行了大量观测. 结果表明,液层中的流场结构经历了多种状态的转变,该过程会受到液层厚度的影响. 当液层厚度较小时,比如当d=2:8mm 时,随着液池两端温差的增大,液层中的流场结构会经历单胞对流到双胞对流再到多胞对流的转变,到达多胞对流状态之后,继续增大温差,对流涡胞的数量会有所减少,而当温差进一步增大到一定程度以后,整个液层转变为三维非定常流动;当液层厚度较大时,比如当d=4:5mm 时,随着温差的增大,流动模式的转变主要体现在水平截面流场截面上面,当温差增大到一定程度以后,在靠近高温端的附近区域会出现具有明显三维效应的"梭形结构",该梭形结构的尺寸随着温差的增大而增长,并在温差超过某个临界值时失去对称性,整个液层转变为三维非定常流动.     

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

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

液桥内热质耦合对流不稳定性及旋转磁场法控制

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

微重力下液封液桥内振荡热毛细对流的数值模拟

对微重力条件下液封液桥(流体采用KF-96硅油和FC-70氟化物,桥高1.4~4毫米,内层直径2或3毫米,外层直径5毫米)和单层液桥(流体采用KF-96硅油,桥高1~1.6毫米,直径2毫米)内的热毛细对流进行了数值模拟,模拟条件与Majima等的实验条件相同。从理论上证实了大温差条件下将出现振荡热毛细对流,确定了发生振荡的临界条件并与相关实验结果进行了对比,同时还计算了振荡频率。     

半浮区液桥热毛细振荡流

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

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

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

三维数值模拟内环加热液池内热毛细对流

为了深入理解热毛细对流不稳定性的特征,数值模拟了内环壁加热、外环壁被冷却的硅油环形浅液池内热毛细对流.发现在内外壁问温差相同的情况下,在内壁加热的液层里,流体表面的温度梯度高于外壁加热情况下10%左右;相应地,发生热流体波的临界Marangoni数Mac=3.06×105(或临界温差△Tc=4.58K)低于外壁加热情况下的Mac=3.34×105(或△Tc=5.0 K).超过临界条件后,内壁加热液池内形成振荡的多涡胞热流体波,它由共存的两组波数相同、频率和振幅相同、沿周向传播方向相反、沿径向凹向高温壁传播的、弯曲螺旋状的热流体波组成,并相互干涉形成点状波纹.     

Bifurcation of thermocapillary convection in a shallow annular pool of silicon melt

In order to understand the nature of surface patterns on silicon melts in industrial Czochralski furnaces, we conducted a series of unsteady three-dimensional numerical simulations of thermocapillary convections in thin silicon melt pools in an annular container. The pool is heated from the outer cylindrical wall and cooled at the inner wall. Bottom and top surfaces are adiabatic. The results show that the flow is steady and axisymmetric at small temperature difference in the radial direction. When the temperature difference exceeds a certain threshold value, hydrothermal waves appear and bifurcation occurs. In this case, the flow is unsteady and there are two possible groups of hydrothermal waves with different number of waves, which are characterized by spoke patterns traveling in the clockwise and counter-clockwise directions. Details of the flow and temperature disturbances are discussed and number of waves and traveling velocity of the hydrothermal wave are determined. The project supported by the National Natural Science Foundation of China (50476042) and the Scientific Research Foundation for the Returned Overseas Chinese Scholars, Ministry of Education of China. The English text was polished by Yunming Chen.     

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.     

Onset of Rayleigh–Bénard convection in cold water near its density maximum in vertical annular containers

A set of three-dimensional numerical simulations of Rayleigh–Bénard convection in cold water near its density maximum in vertical annular containers is performed with the aim of determining the critical Rayleigh number at the onset of convection and the primary flow patterns for different geometric dimensions and density inversion parameters. The Prandtl number of cold water is about 11.57. The annular container is heated from below and cooled from above. The inner and outer sidewalls are considered to be perfectly adiabatic. The results obtained show that the critical Rayleigh number at the onset of convection increases with increase in the density inversion parameter and the radius ratio and with decrease in the aspect ratio. When the radius ratio is small, the flow patterns in vertical annular containers are similar to those in cylindrical containers. At large radius ratios the flow pattern is relatively simple, with several convective rolls observable along the azimuthal direction and similar with those characteristic of Rayleigh–Bénard convection in the Boussinesq fluid. The stratified flow phenomenon is found to exist at moderate values of the density inversion parameter. The results are compared with those obtained in the Boussinesq fluid to reveal the effect of the density inversion parameter.     

Transient natural convection between concentric and vertically eccentric spheres with mixed boundary conditions

Transient analysis has been investigated numerically to determine heat transfer by natural convection between concentric and vertically eccentric spheres with constant heat flux on the inner wall and a specified isothermal temperature on the outer wall. The governing equations, in terms of vorticity, stream function and temperature are expressed in a spherical polar coordinate system. The alternating direction implicit method and the successive over-relaxation techniques are applied to solve the finite difference form of governing equations. A physical model is introduced which accounts for the effects of fluid buoyancy as well as eccentricity of the outer sphere. Transient solutions of the entire flow field are obtained for a range of modified Rayleigh number (10³<Ra?<5×10⁵), for a Prandtl number of 0.7 and a radius ratio of 2.0, with the outer sphere near the top and bottom of the inner sphere (ε=±0.625). Results of the parametric study conducted further reveal that the heat and flow fields are primarily dependent on the modified Rayleigh number and the eccentricity of the spherical annulus. The results of average Nusselt numbers are also compared with the results obtained for flow between two isothermal spheres.     

A Modified Nested-Gridding for Upscaling–Downscaling in Reservoir Simulation

This article reports a numerical study of double-diffusive convection in a fluid-saturated vertical porous annulus subjected to discrete heat and mass fluxes from a portion of the inner wall. The outer wall is maintained at uniform temperature and concentration, while the top and bottom walls are adiabatic and impermeable to mass transfer. The physical model for the momentum equation is formulated using the Darcy law, and the resulting governing equations are solved using an implicit finite difference technique. The influence of physical and geometrical parameters on the streamlines, isotherms, isoconcentrations, average Nusselt and Sherwood numbers has been numerically investigated in detail. The location of heat and solute source has a profound influence on the flow pattern, heat and mass transfer rates in the porous annulus. For the segment located at the bottom portion of inner wall, the flow rate is found to be higher, whereas the heat and mass transfer rates are higher when the source is placed near the middle of the inner wall. Further, the average Sherwood number increases with Lewis number, while for the average Nusselt number the effect is opposite. The average Nusselt number increases with radius ratio (λ); however, the average Sherwood number increases with radius ratio only up to λ = 5, and for λ > 5 , the average Sherwood number does not increase significantly.     

环形浅液池内热流体波的本质特征

为了了解环形浅液池内热毛细对流的特征,利用有限容积法对外壁受热、内壁冷却、厚度为1 mm的环形浅液池内硅油的热毛细对流进行了三维数值模拟.结果表明,当Marangoni(Ma)数小于临界值时,随着Ma数的增加,内、外壁附近的温度梯庹上升,稳定的二维轴对称流动的径向速率增加;超过临界Ma数后,漉动转化为三维振荡对流,形成热流体波.沿径向运动的同时,伴随着热流体波的传播流体质点成对地绕顺时针和逆时针方向旋转.热流体波的周向传播速度较快,而流体质点的周向速度很小.分析发现,热流体波为对数螺线形波纹,其传播角为常数;随着Ma数的增大,传播角增大.     

Numerical investigation of a spatially developing turbulent natural convection boundary layer along a vertical heated plate

Large-eddy simulation (LES) on a spatially developing natural convection boundary layer along a vertical heated plate was conducted. The heat transfer rate, friction velocity, mean velocity and temperature, and second-order turbulent properties both in the wall-normal and the stream-wise direction showed reasonable agreement with the findings of past experiments. The spectrum of velocity and temperature fluctuation showed a -2/3-power decay slope and -2-power decay slope respectively. Quadrant analysis revealed the inclination on Q1 and Q3 in the Reynolds stress and turbulent heat flux, changing their contribution along the distance from the plate surface. Following the convention, we defined the threshold region where the stream-wise mean velocity takes local maximum, the inner layer which is closer to the plate than the threshold region, the outer layer which is farther to the plate than the threshold region. The space correlation of stream-wise velocity tilted the head toward the wall in the propagating direction in the outer layer; on the other hand, the correlated motion had little inclination in the threshold region. The time history of the second invariant of gradient tensor Q revealed that the vortex strength oscillates both in the inner and the outer layers in between the laminar and the transition region. In the turbulent region, the vortex was often dominant in the outer layer. Instantaneous three-dimensional visualization of Q revealed the existence of high-speed fluid parcels associated with arch-shape vortices. These results were considered as an intrinsic structure in the outer layer, which is symmetrical to the structure of canonical smooth/rough wall bounded layer flow in forced convection.     

Experimental investigation of micro-scale temperature transients in sub-cooled flow boiling on a horizontal heater

Surface temperature fluctuations that occur locally underneath departing bubbles in pool boiling are shown to result in local heat transfer coefficients ranging from 1 to 10 kW/cm². These estimates were reported in the literature involved both numerical and experimental approaches. Significantly higher heat fluxes are associated with flow boiling than pool boiling under similar conditions of wall superheat and liquid subcooling (e.g. at boiling inception and at critical heat flux). These enhancements are primarily caused by the convective transport, acceleration/distortion of the bubble departure process as well as the resultant potential enhancement of the local surface temperature fluctuations.In this study we measure the surface temperature fluctuations using temperature micro/nano-sensors fabricated on a silicon wafer during flow boiling on the silicon wafer which is heated from below. The silicon wafer is clamped on a constant heat flux type calorimeter consisting of a vertical copper cylinder with embedded cartridge heaters and K-type thermocouples. Micro/nano-thermocouples (thin film thermocouples or “TFT”) are fabricated on the surface of the silicon wafer. High speed data acquisition apparatus is used to record temperature data from the TFT at 1 kHz. A fluorinert was used as the test fluid (PF-5060, manufacturer: 3M Co.). The calorimeter and surface temperature measurement apparatus is housed in a test section with glass walls for visual observation. The liquid is pumped from a constant temperature bath to maintain a fixed subcooling during the experiments under steady state conditions. The transient temperature data from the FFT array during flow boiling on the silicon wafer is analyzed using fast Fourier transform (FFT). The FFT data is analyzed as a function of the wall heat flux and wall superheat. The number of temperature peaks in the FFT data is observed to increase with increase in wall heat flux and the peaks are found to cover a wider spectrum with peaks at higher frequencies with enhancement of heat flux. The surface temperature fluctuations, especially at small length and time scales, are perturbed potentially by the coupled hydrodynamic and thermal transport processes, resulting in enhanced local and global heat flux values. Boiling incipience condition and the flow boiling data are compared with correlations reported in the literature.     

Flow pattern and boiling heat transfer of CO2 in horizontal small-bore tubes

Increasing attention has been focused on carbon dioxide (CO₂) heat pump system where the temperature level is rather low, while the operating pressure is rather high. In this system, the density difference between vapor and liquid becomes rather small, which significantly affects flow patterns. Low surface tension and latent heat also have significant influence on two-phase flow patterns and heat transfer. This paper describes experimental and numerical investigation on flow patterns and heat transfer characteristics of boiling flow CO₂ at high pressure in horizontal small-bore tubes ranging from 1.0 mm to 3.0 mm I.D. Even though the density difference is rather small at high pressure, phase stratification takes place, which leads to the intermittent dryout at the upper wall. So far developed discrete bubble model by the authors for vertical flows is modified so as to include horizontal flow mechanisms. The predicted flow patterns with this new model agree on the whole with the experimental observation.