Numerical simulation of drop Marangoni migration under microgravity

Thermocapillary motion of a drop in a uniform temperature gradient is investigated numerically. The three-dimensional incompressible Navier-Stokes and energy equations are solved by the finite-element method. The front tracking technique is employed to describe the drop interface. To simplify the calculation, the drop shape is assumed to be a sphere. It has been verified that the assumption is reasonable under the microgravity environment. Some calculations have been performed to deal with the thermocapillary motion for the drops of different sizes. It has been verified that the calculated results are in good agreement with available experimental and numerical results. (C) 2003 Elsevier Ltd. All rights reserved.

Validity ranges of interfacial wave theories in a two-layer fluid system

In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness d(1), and lower layer thickness d(2), instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehautes plot for free surface waves if water depth ratio r = d(1)/d(2) approaches to infinity and the upper layer water density rho(1) to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of sigma = (rho(2) - rho(1))/rho(2) -> 1.0 and r > 1.0. In the end, several figures of the validity ranges for various interfacial wave theories in the two-layer fluid are given and compared with the results for surface waves.

Drop migration of middle Reynolds number in a vertical temperature gradient

The experimental investigation of the thermocapillary drop migration in a vertical temperature gradient uns performed on ground. Silicon oil and pure soybean oil were used as experimental medium in drops and as continuous phases, respectively, in the present experiment. The drop migration, under the combined effects of buoyancy: and thermocapillarity, was studied for middle Reynolds numbers in order of magnitude O(10(1)). The drop migration velocities depending on drop diameters were obtained. The present experimental results show relatively small migration velocity in comparison with the one suggested by Young et nl. for linear theory of small Reynolds number. An example of flow patterns inside the drop was observed by PIV method.

Numerical investigation on a simulation model of floating zone convection

A simulation model with adiabatic condition at the upper rod and constant temperature at the lower rod is studied numerically in this paper. The temperature distribution in a simulation model is closer to the one in the half part of a floating full zone in comparison with the one in a usual floating half zone model with constant temperature at both rods, because the temperature distribution of a floating full zone is symmetric for the middle plane in a microgravity environment. The results of the simulation model show that the temperature profiles and the how patterns are different from those of the usual floating half zone model. Another type of half zone model, with a special non-uniform temperature distribution at the upper rod and constant temperature at the lower rod, has been suggested by recent experiments. The temperature boundary condition of the upper rod has a maximum value in the center and a lower value near the free surface. This modified simulation model is also simulated numerically in the present paper. Copyright (C)1996 Elsevier Science Ltd.

Numerical simulation of g-jitter effects on half floating zone convection under microgravity environment

The g-jitter effects on the thermocapillary convection in liquid bridge of floating half zone were studied by numerical simulation for unsteady and axi-symmetric model in the cylindrical coordinate system. The g-jitter field was given by a steady microgravity field in addition to an oscillatory low-gravity field, and the effects on the flow field, temperature distribution and free surface deformation were analyzed numerically.

Concentration distribution in crystallization from solution under microgravity

Concentration distribution in crystallization from solution under microgravity is numerically studied. A quasi-steady state growth and dissolution in a 2D rectangular enclosure filled with sodium chlorate (NaClO3) aqueous solution, in which one wall is the growth surface of the crystal and the opposite one is the dissolution surface, is considered. The solute transport process at the growth surface is described by the diffusion-reaction theory with finite interface kinetics coefficient. The results show that the concentration at the growth surface is supersaturated and the supersaturation distribution is of non-uniformity, i.e. the supersaturation in a region facing an incoming flow is high. On the other hand, the non-uniformity of supersaturation at the growth surface is closely related to the gravity level even under microgravity, it exponentially increases as the thermal Rayleigh number on behalf of the gravity level rises.

Hardware development and ground-based studies for drop shaft experiment on oscillatory convection in half floating-zone

Coordinated measurement of temperature, velocity and free surface oscillation were obtained by using the drop shaft facility for microgravity experiments of half floating zone convection. The ground-based studies gave transition from steady to oscillatory convection for multi-quantities measurement.

Effects of g-jitter on the critical marangoni number

A half floating zone is fixed on a vibrational deck, which supports a periodical applied acceleration to simulate the effect of g-jitter. This paper deals with the effects of g-jitter on the fluid fields and the critical Marangoni number, which describes the transition from a forced oscillation of thermocapillary convection into an instability oscillatory convection in a liquid bridge of half floating zone with top rod heated. The responses of g-jitter field on the temperature profiles and flow pattern in the liquid bridge were obtained experimentally. The results indicated that the critical Marangoni number decreases with the increasing of g-jitter effect and is slightly smaller for higher frequency of g-jitter with fixed strength of applied gravity.

Influence of buoyancy on floating-zone convection of small bond number

Gravity may influence the velocity and temperature distributions, bouyancy may induce Rayleigh instability and the instability may be excited due to the change of free surface shape associating with gravity in the thermocapillary convection. These effects have been studied in the present paper. The results show that gravity may have an important effect in thermocapillary oscillatory convection even for the cases of small Bond number experiments either on the ground or in space.

Air flow over a mountain and the convective boundary-layer

In this paper, the analytical model coupling the convective boundary layer (CBL) with the free atmosphere developed by Qi and Fu (1992) is improved. And by this improved model, the interaction between airflow over a mountain and the CBL is further discussed. The conclusions demonstrate: (1) The perturbation potential temperatures in the free atmosphere can counteract the effect of orographic thermal forcing through entraining and mixing in the CBL. If u(M)BAR > u(F)BAR, the feedback of the perturbation potential temperatures in the free atmosphere is more important than orographic thermal forcing, which promotes the effect of interfacial waves. If u(M)BAR < u(F)BAR, orographic thermal forcing is more important, which makes the interfacial height and the topographic height identical in phase, and the horizontal speeds are a maximum at the top of the mountain. (2) The internal gravity waves propagating vertically in the free atmosphere cause a strong downslope wind to become established above the lee slope in the CBL and result in the hydraulic jump at the top of the CBL. (3) With the CBL deepening, the interfacial gravity waves induced by the potential temperature jump at the top of the CBL cause the airflow in the CBL to be subcritical.

The effect of the surface tension on non-propagating solitary waves

In this paper, the effect of the surface tension is considered carefully in the study of non-propagating solitary waves. The parameter plane of the surface tension and the fluid depth is divided into three regions; in two of them a breather soliton can be produced. In literature the parameters of breather solitons are all in one of the parameter regions. The new region reported here has been confirmed by our experiments. In the third region, the theoretical solution is a kink soliton, but a kind of the non-propagating solitary wave similar to the breather soliton was found in our experiments besides the kink soliton.

A study on stratification of concentration in crystal-growth of alpha-liio3 by holographic-interferometry

A new thermoplastic-photoconductor laser holographic recording system has been used for real-time and in situ observation of alpha-LiIO3 crystal growth. The influence of crystallization-driven convection on the concentration stratification in solution has been studied under gravity field. It is found that the stratification is closely related to the seed orientation of alpha-LiIO3 crystal. When the optical axis of crystal seed C is parallel to the gravity vector g, the velocity of the concentration stratification is two times larger than that in the case of C perpendicular-to g. It needs 40 h for the crystalline system of alpha-LiIO3 to reach stable concentration distribution (expressed as tau) at 47.6-degrees-C. The time tau is not sensitive to the seed orientation. Our results provide valuable data for designing the crystal growth experiments ia space.

Crystal-growth in floating zone with phase-change and thermosolutal convections

Floating zone crystal growth in microgravity environment is investigated numerically by a finite element method for semiconductor growth processing, which involves thermocapillary convection, phase change convection, thermal diffusion and solutal diffusion. The configurations of phase change interfaces and distributions of velocity, temperature and concentration fields are analyzed for typical conditions of pulling rates and segregation coefficients. The influence of phase change convection on the distribution of concentration is studied in detail. The results show that the thermocapillary convection plays an important role in mixing up the melt with dopant. The deformations of phase change interfaces by thermal convection-diffusion and pulling rods make larger variation of concentration field in comparison with the case of plane interfaces.

Excitation mechanism of thermocapillary oscillatory convection

It is suggested that the oscillation of thermocapillary convection may be excited by the buoyancy instability. By means of numerical simulation of the finite-element method, the temperature distributions in the liquid bridge are qualitatively analyzed. The temperature gradient in a certain flow region of liquid bridge may turn to be parallel to the direction of gravity when the temperature difference △T between two boundary rods of liquid bridge is larger than the critical value. The buoyancy instability may be excited, and then the thermocapillary oscillatory convection appears, as the temperature difference increases further. The distribution of the critical Marangoni number in the micro-gravity environment is derived from the data on the ground experiments. The results show that the onset of thermocapillary oscillatory convection is delayed in the case of smaller typical scale of liquid bridge and lower gravity environment.

Approximate theory of natural-convection with small grashof number

To gain some insight into the behaviour of low-gravity flows in the material processing in space, an approximate theory has been developed for the convective motion of fluids with a small Grashof number Gr. The expansion of the variables into a series of Gr reduces the Boussinesq equation to a system of weakly coupled linearly inhomogeneous equations. Moreover, the analogy concept is proposed and utilized in the study of the plate bending problems in solid mechanics. Two examples are investigated in detail, i. e. the 2-dimensional steady flows in either circular or square infinite closed cylinder, which is horizontally imposed at a specified temperature of linear distribution on the boundaries. The results for stream function ψ, velocity u and temperature T are provided. The analysis of the influences of some parameters such as the Grashof number Gr and the Prandtl number Pr, on motions will lead to several interesting conclusions. The theory seems to be useful for seeking for an analytical solutions. At least, it will greatly simplify the complicated problems originally governed by the Navier-Stokes equation including buoyancy. It is our hope that the theory might be applicable to unsteady or 3-dimensional cases in future.