直拉法晶体生长过程非稳态流体热流耦合

为了改善复杂对流形态下的晶体生长品质, 提出了一种改进的格子Boltzmann方法研究非稳态熔体流动和传热的耦合性质. 该方法基于不可压缩轴对称D2Q9模型, 构建了包含旋转惯性力和热浮力等外力项的演化关系, 实现了对轴对称旋转流体的速度、温度和旋转角速度的计算与分析. 结果表明, 非稳态熔体中的流、热耦合性质与格拉斯霍夫数和雷诺数的相互作用有关; 通过调节高雷诺数, 可有效抑制熔体中的自然对流, 改善温度分布, 有助于提高单晶的品质. 数值计算结果与实际硅单晶生长试验均证明了所提方法的正确性及有效性.     

Leap behavior of ultrasonic standing waves in the liquids

The generation and behavior of ultrasonic standing waves was modeled using the light cut method for transparent fluid. The oscillations of the fluid surface in initial moment of switching on ultrasound and appearance of standing wave channel were observed. The effect of continuous fluid depth decrease and increase on the behavior of ultrasonic standing wave channel was studied. The ultrasonic standing wave channel floated in the liquid between of the crucible bottom and fluid surface and discretely changed its height by half ultrasonic wavelength with the decrease or increase of the liquid level. This channel had the behavior of a “quasi solid state” and damped of convection.     

DNS of buoyancy-dominated turbulent flows on a bluff body using the immersed boundary method

A novel immersed boundary (IB) method has been developed for simulating multi-material heat transfer problem – a cylinder in a channel heated from below with mixed convection. The method is based on a second-order velocity/scalar reconstruction near the IB. A novel algorithm has been developed for the IB method to handle conjugate heat transfer. The fluid–solid interface is constructed as a collection of disjoint faces of control volumes associated to different material zones. Coupling conditions for the material zones have been developed such that continuity and conservation of the scalar flux are satisfied by a second-order interpolation. Predictions of the local Nusselt number on the cylinder surface show good agreement with the experimental data. The effect of the Boussinesq approximation on this problem was also investigated. Comparison with the variable density formulation suggests that, in spite of a small thermal expansion coefficient of water, the variable density formulation in a transitional flow with mixed convection is preferable.     

MEASUREMENT TECHNIQUES FOR HORIZONTAL FLOW BOILING HEAT TRANSFER WITH PURE AND MIXED REFRIGERANTS

Abstract

A test rig and measurement techniques for horizontal flow boiling of pure and mixed refrigerants are described. Local heat transfer coefficients were measured for R22/R114 and R12/R152a mixtures as well as the corresponding pure components. The test section consists of an S-m-long, 9.1-mm-i.d., electrically heated stainless steel tube and has the distinctive feature of variable heated length. Details of pressure, composition, and fluid and wall temperature measurements are discussed. The composition of subcooled liquid entering test section varied from test to test, and it is recommended that it be measured for each test. For the R22/R114 mixture, local composition measurements in the annular liquid film revealed a composition variation of up to 0.07 mole fraction around the circumference of the heated tube.     

The influence of heat radiation on mixed convection boundary layer flow of a viscoelastic fluid over a circular cylinder with constant surface temperature

The influence of mixed convection boundary layer flow of a viscoelastic fluid over an isothermal horizontal circular cylinder has been analyzed. The boundary layer equations governing the problem are reduced to dimensionless nonlinear partial differential equations and then solved numerically using Keller-box method. Skin friction coefficient and Nusselt number are emphasized specifically. These quantities are displayed against curvature parameter. Effects of mixed convection parameter and radiation-conduction parameter on skin friction coefficient and Nusselt number are illustrated through graphs and table. The boundary layer separation points along the surface of cylinder are also calculated with/without radiation, and a comparison is shown. The presence of radiation helps to reduce the skin friction coefficient in opposing flow case and enhances it for assisting flow case. The increase in value of radiation-conduction parameter helps increase the value of skin friction coefficient and Nusselt number for viscoelastic fluids. The boundary layer separation delays due to thermal radiation.     

Experimental Investigation of Thermocapillary Convection in a Liquid Layer with Evaporating Interface

Thermocapillary convection coupling with the evaporation effect of evaporating liquids is studied experimentally. This study focused on an evaporation liquid layer in a rectangular cavity subjected to a horizontal temperature gradient when the top evaporating surface is open to air, while most previous works only studied pure thermocapillary convection without evaporation. Two liquids with different evaporating rates are used to study the coupling of evaporation and thermocapillary convection, and the interfacial temperature profiles for different temperature gradients are measured. The experimental results indicate evidently the influence of evaporation effect on the thermocapillary convection and interfacial temperature profiles. The steady multicellular flow and the oscillatory multicellular flow in the evaporation liquid layer are observed by using the particle-image-velocimetry method.     

Instabilities of thermocapillary-buoyancy convection in open rectangular liquid layers

This article presents the experimental investigation on instabilities of thermocapillary-buoyancy convection in the transition process in an open rectangular liquid layer subject to a horizontal temperature gradient. In the experimental run,an infrared thermal imaging system was constructed to observe and record the surface wave of the rectangular liquid layer. It was found that there are distinct convection longitudinal rolls in the flow field in the thermocapillary-buoyancy convection transition process. There are different wave characterizations for liquid layers with different thicknesses. For sufficiently thin layers, oblique hydrothermal waves are observed, which was predicted by the linear-stability analysis of Smith Davis in 1983. For thicker layers, the surface flow is distinct and intensified, which is because the buoyancy convection plays a dominant role and bulk fluid flow from hot wall to cold wall in the free surface of liquid layers. In addition, the spatiotemporal evolution analysis has been carried out to conclude the rule of the temperature field destabilization in the transition process.     

Free surface flow between two horizontal concentric cylinders

Results are reported on a combined experimental and numerical investigation of a free surface flow at small Reynolds numbers. The flow is driven by the rotation of the inner of two horizontal concentric cylinders, with an inner to outer radius ratio of 0.43. The outer cylinder is stationary. The annular gap is partially filled, from 0.5 to 0.95 full, with a viscous liquid leaving a free surface. When the fraction of the annular volume filled by liquid is 0.5, a thin liquid film covers the rotating inner cylinder and reenters the liquid pool. For relatively low rotation speeds, the evolution of the film thickness is consistent with the theory for a plate being withdrawn from an infinite liquid pool. The overall liquid flow pattern at this condition consists of two counter-rotating cells: one is around the inner cylinder and the other with weaker circulation rate is in the bottom part of the annulus and nearly symmetric about the vertical axis. With increasing rotation rate, the free surface becomes more deformed, and the dynamics of the stagnation line and the cusp line dividing the cells are tracked as quantitative measures of the interface shape. In addition, the recirculating flow cells lose symmetry and the cusp deforms the free surface severely. A comparison of numerically computed flow which describes the interface by a phase-field method confirms the dynamics of the two cells and the interface deformation. For filling fraction 0.75, the liquid level is slightly above the inner cylinder and a significant decrease in size of the bottom cell with increasing rotation rate is found. For filling fractions approaching unity, the liquid flow consists of one single cell and the surface deformation remains small.     

蒸发液面热对流与温度不连续现象的实验研究

实验研究了存在蒸发界面的水平液层中的热毛细对流对汽/液界面处温度不连续的影响.对于单纯的热毛细对流从理论和实验已有深入研究,但目前国际对蒸发与热毛细对流的耦合尚缺乏研究.特别是近期C.A.Ward等人研究发现了蒸发汽液界面处的温度不连续现象.本文以存在蒸发界面的水平薄液层为研究对象,测量了蒸发界面处的温度分布,研究了普...     

湿消解-原子荧光法测定煤中硒和砷

煤中微量元素砷、硒的含量与环境安全和居民健康密切相关,快速准确测定煤中砷、硒含量具有重要的理论和实践意义。现有中国国标(GB/T 16415—2008, GB/T 3058—1996)测定煤中总砷、总硒含量为干法分析法,该方法具有操作步骤复杂,误差大的缺陷。为了准确、快速测定煤中总砷、总硒的含量,克服干法分析方法的缺陷,在长期摸索的基础上,借鉴土壤中砷、硒的测试方法,建立了电热板混酸消解-氢化物发生-原子荧光测定煤中总砷、总硒的方法,即湿法分析方法。方法：准确称取0.05～0.10 g煤样(100目),转移入100 mL玻璃烧杯中,加入10 mL硝酸、2 mL高氯酸,加盖玻璃表面皿,于通风橱中静置过夜。次日,在电热板(180 ℃)上加热消解至白色浓烟冒出,液体透明,残余样品呈白色或者灰白色。取下冷却后,加入6 mol·L-1盐酸溶液3 mL,再次放在电热板上加热消解。当再次有白色浓烟开始冒出时,取下冷却后,加入1 mL浓盐酸。少量亚沸水冲洗烧杯和表面皿内壁,将消解液转移入25 mL试管,亚沸水定容至刻度线,摇匀待测硒;取3 mL消解液,转移入15 mL试管,加入1 mL浓盐酸,1 mL硫脲+抗坏血酸还原剂,亚沸水定容至刻度线,摇匀待测砷;用原子荧光法测定消解液中砷、硒含量。用该方法对煤成分分析标准物质(GBW11115, GBW11117)中砷、硒含量进行测定,测得砷回收率在99.7%～100.3%之间,相对标准偏差范围为5.6%～6.0%;硒相对标准偏差范围为11.1%～13.5%;对应砷和硒的检测限分别为0.05和0.01 μg·L-1。该方法可以准确测定煤中砷、硒含量,且较国标干法分析方法具有操作步骤少、准确性高等优点。     

Influence of the meridional flow and thermomagnetic convection on characteristics of magnetic fluid seal

For the low-speed magnetic fluid seals, the influence of the meridional flow, induced by the shaft rotation, on the distribution of magnetic particles concentration, is studied. Influence of the thermomagnetic convection on the structure of this flow and on the temperature distribution in high-speed magnetic fluid seals is investigated also. The problems were examined by numerical methods. It is discovered that even very slow rotation of the shaft homogenises distribution of the magnetic particles concentration in the seal and thereby enlarges its operation life. For high-speed seals thermomagnetic convection provides the penetration of the fluid flow in the region of the narrow gap and levels off the temperature distribution decreasing its maximum value and thereby enlarges its operation life too. It is found also that the influence of thermomagnetic convection grows with the viscosity increasing.     

Natural convection with damped thermal flux in a vertical circular cylinder

Natural convection flows of an incompressible Newtonian fluid inside a circular cylinder are studied. The heat transfer process is described by a generalized fractional constitutive equation for the thermal flux-temperature gradient. Caputo time-fractional derivative operator, which provides the damping of thermal flux, is considered into the studied model.Analytical solutions to the fluid temperature, thermal flux, fluid velocity and volume flow rate are obtained with the integral transforms method (Laplace transform and finite Hankel transform).Temperature behaviors for small and large values of the time t, as well as the post-transient and transient velocity components are determined. The influence of the memory parameter (the order of the time-fractional derivative) on the temperature, thermal flux, velocity and the volume flow rate is numerically and graphically studied.     

磁场对水平温度梯度下双层流体热对流的影响

在水平温度梯度下,双层流体交界面的表面张力会出现梯度,驱动热毛细对流运动,造成热剪切层内的扰动.本文数值模拟了不同重力条件下,双层流体内的对流现象,得出了在微重力时,对流运动将引起热剪切层内强烈的扰动.为了减弱这种扰动,我们利用磁场对流体的运动进行控制.为此,又对微重力条件下,不同方向应用磁场下的热剪切层内扰动行为进行了数值研究,结果显示,磁场对热剪切层稳定性有促进作用,加入法向的应用磁场最为有效.     

Experimental Study on Liquid Free Surface in Buoyant-Thermocapillary Convection

We investigate the surface deformations of buoyant-thermocapillary convection in a rectangular cavity due to gravity and temperature gradient between the two sidewalls. The cavity is 52mm×42 mm in horizontal cross section, the thickness of liquid layer h is changed from 2.5 mm to 6.5 mm. Surface deformations of h = 3.5 mm and 6.0mm are discussed and compared. Temperature difference is increased gradually, and the flow in the liquid layer will change from stable convection to unstable convection. Two kinds of optical diagnostic system with image processor are developed for study of the kinetics of buoyant-thermocapillary convection, they give out the information of liquid free surface. The quantitative results are calculated by Fourier transform and correlation analysis, respectively. With the increasing temperature gradient, surface deformations calculated are more declining. It is interesting phenomenon that the inclining directions of the convections in thin and thick liquid layers are different. For a thin layer, the convection is mainly controlled by thermocapillary effect. However, for a thick layer, the convection is mainly controlled by buoyancy effect. The surface deformation theoretically analysed is consistent with our experimental results. The present experiment proves that surface deformation is related to temperature gradient and thickness of the liquid layer. In other words, surface deformation lies on capillary convection and buoyancy convection.     

微重力下凝结和沸腾着的汽泡周围流场

1引言尽管对沸腾换热已有大量的研究，但对沸腾过程中驱动汽泡外流动的机理仍没有透彻理解，通常认为，汽泡外流动是由液体自然对流引起的。但是，近来的微重力实验表明［‘］；在没有自然对流的微重力场中，沸腾换热同样很剧烈。为了弄清正常重力和微重力环境中驱动换热的机理，很有必要作进一步的研究。本文数值分析了微重力环境下驱动壁面汽泡周围液体流动的基本机理，深入分析了Marangoni效应的影响，同时考虑表面凝结和蒸发过程的作用。2理论模型汽泡外流动可由二维层流N－S方程及烙方程来描述，流动是由Marangoni边界条件驱动，该条…     

Microfluidic mixing via acoustically driven chaotic advection

Mixing presents a notoriously difficult problem in small amounts of fluids. Herein, surface acoustic waves provide a convenient technique to generate time-dependent flow patterns. These flow patterns can be optimized in such a way that advected particles are mixed most efficiently in the fluid within a short time compared to the time pure diffusion would take. Investigations are presented for the mixing efficiency of a flat cylinder that is driven by two surface acoustic waves. The experimental results favorably agree with model calculations of the flow patterns and the advective transport.     

流向对管束外混合对流影响的数值模拟

用格子玻尔兹曼方法对微圆管管束外混合对流的换热情况进行了数值模拟研究。在充分验证程序正确性的基础上,对混合流动绕流时不同壳侧进液角度,不同雷诺数Re以及瑞利数Ra等方面对传热的影响进行模拟研究,对衡量传热效果的Nu数及温度场进行了对比,总结出了换热规律。为高效微管换热器的开发提供了理论指导。     

MHD mixed convection flow near the stagnation-point on a vertical permeable surface

The steady magnetohydrodynamic (MHD) mixed convection boundary layer flow of a viscous and electrically conducting fluid near the stagnation-point on a vertical permeable surface is investigated in this study. The velocity of the external flow and the temperature of the plate surface are assumed to vary linearly with the distance from the stagnation-point. The governing partial differential equations are first transformed into ordinary differential equations, before being solved numerically by a finite-difference method. The features of the flow and heat transfer characteristics for different values of the governing parameters are analyzed and discussed. Both assisting and opposing flows are considered. It is found that dual solutions exist for both cases, and the range of the mixed convection parameter for which the solution exists increases with suction.     

Mathematical modeling of convective-conductive heat transfer in a rectangular domain in a conjugate statement

The results of mathematical modeling of convection of a viscous incompressible liquid in a rectangular domain with sources of mass input and output are presented. A conjugate statement within the framework of the Boussinesq approximation is used. The regimes of forced and mixed convection in a domain have been investigated. The domain has two vertical walls and one horizontal wall of finite thickness, two zones of liquid input and output, and a free surface. A plane nonstationary problem within the framework of the Navier-Stokes model for the liquid phase and the heat conduction equation for the solid phase are considered. The distributions of the hydrodynamic parameters and temperatures characterizing the main regularities of the processes under investigation have been obtained. Circulation flows have been identified. The vortex formation mechanism and the temperature distribution in the solution domain under the regimes of forced and mixed convection and different locations of mass input and output zones have been analyzed. It has been found that natural convection should be taken into account when modeling convective heat transfer with Gr number values from 10⁵ and higher.