Period-doubling bifurcations of the thermocapillary convection in a floating half zone

This study experimentally explored the fine structures of the successive period-doubling bifurcations of the time-dependent thermocapillary convection in a floating half zone of 10 cSt silicone oil with the diameter d (0)=3.00 mm and the aspect ratio A=l/d (0)=0.72 in terrestrial conditions. The onset of time-dependent thermocapillary convection predominated in this experimental configuration and its subsequent evolution were experimentally detected through the local temperature measurements. The experimental results revealed a sequence of period-doubling bifurcations of the time-dependent thermocapillary convection, similar in some way to one of the routes to chaos for buoyant natural convection. The critical frequencies and the corresponding fractal frequencies were extracted through the real-time analysis of the frequency spectra by Fast-Fourier-Transformation (FFT). The projections of the trajectory onto the reconstructed phase-space were also provided. Furthermore, the experimentally predicted Feigenbaum constants were quite close to the theoretical asymptotic value of 4.669 [Feigenbaum M J. Phys Lett A, 1979, 74: 375-378].

Coupling of thermocapillary convection and evaporation effect in a liquid layer when the evaporating interface is open to air

The coupling mechanism of thermocapillary convection and evaporation effect in evaporating liquids was studied experimentally. The experiments were carried out to study a thin evaporating liquid layer in a rectangular test cell when the upper surface was open to air. By altering the imposed horizontal temperature differences and heights of liquid layers, the average evaporating rate and interfacial temperature profiles were measured. The flow fields were also visualized by PIV method. For comparison, the experiments were repeated by use of another two non-evaporating liquids to study the influence of evaporation effect. The results reveal evidently the role that evaporation effect plays in the coupling with thermocapillary convection.

Tunable ordered droplets induced by convection in phase-separating P2VP/PS blend film

Ordered hexagonal droplets patterns in phase-separating polymeric blend films of polystyrene and poly(2-vinylpyridine) (PS/PVP) formed due to the convection effect by solvent evaporation. The influences of PS molecular weight, solvent evaporation rate, and the weight ratio of PS to PVP on the PVP-rich domains pattern formation and distributions were investigated by atomic force microscope (AFM). Only in an appropriate range of molecular weight of PS, can the ordered pattern form. Too low or too high molecular weight of PS led no ordered pattern due to the viscosity effects. The increase of solvent evaporation rate decreased the mean radius of the PVP-rich domains and the intervals between the centers of the domains due to the enhancement of the viscosity on the top layer of the fluid film. The increase of the weight ratio of PS to PVP decreased mean radius of the PVP-rich domains whereas the intervals between the centers of droplets remained constant. Therefore, the size and the distributions of ordered patterns can be tuned by the polymer molecular weight, the weight ratio of the two components and the solvent evaporation rate.

A HYBRID METHOD OF FRACTIONAL STEPS WITH PREDICTOR-CORRECTOR DIFFERENCE-PSEUDOSPECTRUM FOR NUMERICAL-SOLUTION OF THE CONVECTION-DOMINATED FLOW PROBLEMS

A fractional-step method of predictor-corrector difference-pseudospectrum with unconditional L(2)-stability and exponential convergence is presented. The stability and convergence of this method is strictly proved mathematically for a nonlinear convection-dominated flow. The error estimation is given and the superiority of this method is verified by numerical test.

Effects of buoyancy on Langmuir circulation

Based on the Navier-Stokes equation, an equation describing the Langmuir circulation is derived by a perturbation method when the influences of Coriolis force and buoyancy force are both considered. The approach used in the analysis is similar to the works carried out by Craik and Leibovich [J. Fluid Mech. 73 (1976) 401], Leibovich [J. Fluid Mech. 79 (1977) 715] and Huang [J. Fluid Mech. 91 (1979) 191]. Potential applications of the equation proposed are discussed in the area of Antarctic circumpolar current.

同位网格摄动有限体积格式求解浮力驱动方腔流

利用对流扩散方程的摄动有限体积格式,在Rayleigh数从10~3到10~8的范围内对浮力驱动方腔流动问题作了数值模拟.对流扩散方程的摄动有限体积格式具有一阶迎风格式的简洁形式,使用相同的基点,重构近似精度高,特别是两相邻控制体中心到公共界面的距离相等或不相等,PFV格式公式相同等优点.在数值模拟中,无论均匀网格还是非均匀网格均获得与DSC方法、自适应有限元法、多重网格法等Benchmark解相符较好的数值结果,证明UPFV格式对高Rayleigh数对流传热问题的适用性和有效性.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

DPIV and Interferometry Technique for Velocity Field and Concentration Fields Measurement in Crystal Growth

The property of crystal depends seriously on the solution concentration distribution near the growth surface of a crystal. However, the concentration distributions are affected by the diffusion and convection of the solution. In the present experiment, the two methods of optical measurement are used to obtained velocity field and concentration field of NaClO₃ solution. The convection patterns in sodium chlorate (NaClO₃) crystal growth are measured by Digital Particle image Velocimetry (DPIV) technology. The 2-dimentional velocity distributions in the solution of NaClO3 are obtained from experiments. And concentration field are obtained by a Mach-Zehnder interferometer with a phase shift servo system. Interference patterns were recorded directly by a computer via a CCD camera. The evolution of velocity field and concentration field from dissolution to crystallization are visualized clearly. The structures of velocity fields were compared with that of concentration field.

对流扩散方程的变步长摄动有限差分格式

摄动有限差分(PFD)方法是构造高精度差分格式的一种新方法。变步长摄动有限差分方法是等步长摄动有限差分方法的发展和推广。对需要局部加密网格的计算问题，变步长PFD格式不需要对自变量进行数学变换，且和等步长PFD格式一样，具有如下的共同特点：从变步长一阶迎风格式出发，通过把非微商项(对流系数和源项)作变步长摄动展开，展开幂级数系数通过消去摄动格式修正微分方程的截断误差项求出，由此获得高精度变步长PFD格式。该格式在一、二和三维情况下分别仅使用三、五和七个基点，且具有迎风性。文中利用变步长PFD格式对对流扩散反应模型方程，变系数方程及Burgers方程等进行了数值模拟，并与一阶迎风和二阶中心格式及其问题的精确解作了比较。数值试验表明，与一阶迎风和二阶中心格式相比，变步长PFD格式具有精度高，稳定性与收敛性好的特点。变步长PFD格式与等步长PFD格式相比，变步长PFD解在薄边界层型区域的分辨率得到了明显的提高。

Effects of the orientation of microgravity on the concentration distribution in crystallization from solution

A quasi-steady state growth and dissolution in a 2-D rectangular enclosure is numerically investigated. This paper is an extension to indicate the effects of the orientation of gravity on the concentration field in crystallization from solution under microgravity, especially on the lateral non-uniformity of concentration distribution at the growth surface. The thermal and solute convection are included in this model.

Multiple Axial Cracks in A Coated Hollow Cylinder due to Thermal Shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing. (c) 2005 Elsevier Ltd. All rights reserved.

Tailoring of dendritic microstructure in solidification processing by crucible vibration

We analyzed the effects of both natural convection and forced flows on solid–liquid interface morphology during upward Bridgman solidification of metallic alloys. Experiments were carried out on Al–3.5wt% Ni alloy, for a cylindrical sample. The influence of natural convection induced by radial thermal gradient on solidified microstructure was first analyzed as a function of the pulling rate. Then, the influence of axial vibration on solidification microstructure was experimentally investigated by varying vibration parameters (frequency and amplitude). Experimental results demonstrated that vibrations could be used to either attenuate fluid flow in the melt and obtain a uniform dendritic pattern or to promote a fragmented dendritic microstructure. However, no marked effect was observed for cellular growth. This pointed out the critical role of the mushy zone in the interaction between fluid flow and solidification microstructure.

利用PIV方法测量半浮区液桥热毛细对流的速度场

利用PIV法对硅油液桥热毛细对流的定常速度场进行了实时测量。为了便于测量，液桥上桥端面采取了铜环中嵌透明材料的方法，从液桥的顶部进行观测。当液桥上下桥有温差时，热毛细对流出现；本实验对于不同上下桥的温差，对液桥横剖面内的速度场分别进行了测量，研究外加温差对于流场速度分布的影响；并且在液桥中取了几个典型横截面进行测量，以期对大Pr数液桥的定常速度场有比较全面的定量测量。此外，实验结果也可作为数值模拟计算结果的验证。

Effect of liquid bridge volume on the instability in small-Prandtl-number half zones

A perturbation method is used to examine the linear instability of thermocapillary convection in a liquid bridge of floating half-zone filled with a small Prandtl number fluid. The influence of liquid bridge volume on critical Marangoni number and flow features is analyzed. The neutral modes show that the instability is mainly caused by the bulk flow that is driven by the nonuniform thermocapillary forces acting on the free surface. The hydrodynamic instability is dominant in the case of small Prandtl number fluid and the first instability mode is a stationary bifurcation. The azimuthal wave number for the most dangerous mode depends on the liquid bridge volume, and is not always two as in the case of a cylindrical liquid bridge with aspect ratio near 0.6. Its value may be equal to unity when the liquid bridge is relatively slender.