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.

Experimental investigation on bubble coalescence under nonuniform temperature distribution in reduced gravity

Results on bubble coalescences from the space experiment of thermocapillary bubble migration conducted on board the Chinese 22nd recoverable satellite are presented in this paper. Some coalescences of large spherical bubbles under microgravity are observed through bubbles staying at the upper side of the test cell. The data of bubble coalescence time are recorded and compared with theoretical predictions, which is based on a theory to describe the tendency of coalescence connected to chemical potential difference. It is implied that the theory is applicable for the experimental data of bubble coalescence. Moreover, the angle between the line of two bubble centers and temperature gradient falled mostly in the range 20 degrees-40 degrees. (C) 2007 Elsevier Inc. All rights reserved.

The motion of a small bubble in waves

The motion of a single spherical small bubble due to buoyancy in the ideal fluid with waves is investigated theoretically and experimentally in this article. Assuming that the bubble has no effect on the wave field, equations of a bubble motion are obtained and solved. It is found that the nonlinear effect increases with the increase of the bubble radius and the rising time. The rising time and the motion orbit are given by calculations and experiments. When the radius of a bubble is smaller than 0.5mm and the distance from the free surface is greater than the wave height, the results of the present theory are in close agreement with measurements.

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 simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

On the velocity of buckle propagation in a beam on a nonlinear elastic-foundation

The paper revisits a simple beam model used by Chater et al. (1983, Proc. IUTAM Symp. Collapse, Cambridge University Press) to examine the dynamics of propagating buckles on it. It was found that, if a buckle is initiated at a constant pressure higher than the propagation pressure of the model (P-p), the buckle accelerates and gradually reaches a constant velocity which depends upon the pressure, while if it is initiated at P-p, the buckle propagates at a velocity which depends upon the initial imperfection. The causes for the difference are also investigated.

The tangential discontinuous surface of velocity in solar chromosphere

The data of velocity and magnetic fields in the solar photosphere (5324 angstrom) and the chromosphere (4861 angstrom) clearly show the features of tangential discontinuity of velocity in the chromosphere. The velocity fields in and near the solar active region named No. 88029 by the Huairou Station have been analyzed in detail. A lot of magnetohydrodynamic discontinuous surfaces, especially the tangential discontinuities, are shown from the observations. The calculations of the thickness of discontinuous layer and the evolution time of instability agree with the observational results. The variations of the flow field will directly influence the evolutions and changes of the active region as the magnetic field are coupled closely with the plasma motion.

Simultaneous measurement of temperature, density and velocity in gas-flows by modulated photoluminescence

A novel possibility to determine the temperature, density and velocity simultaneously in gas flows by measuring the average value, amplitude of modulation and phase shift of the photoluminescence excited by a temporally or spatially modulated light source is investigated. Time-dependent equations taking the flow, diffusion, excitation and decay into account are solved analytically. Different experimental arrangements are proposed. Measurements of velocity with two components, and temporal and spatial resolutions in the measurements are investigated. Numerical examples are given for N z with biacetyl as the seed gas. Practical considerations for the measurements and the relation between this method and some existing methods of lifetime measurement are discussed.

Coflow jets with large velocity difference - a new aerodynamic principle of flame stabilization and combustion intensification

A new aerodynamic principle of flame stabilization and combustion intensification, the coflow jets with large velocity difference, is described. One or more small high-velocity jets of air or steam, injected off the axis and in the same direction as the low-velocity main fuel-air flow into the combustor, create a large recirculation zone of high turbulence intensity in which the combustibles and high temperature gases are effectively mixed, so that stable and intensive combustion can be maintained even for fuels with poor ignition. A pulverized coal combustor based on the principle mentioned above is shown to be characteristic of excellent combustoom and a simple structure. A number of precombustors of this type are in operation at some power stations and industrial boilers of China. Using such precombustor, successtul startups and part-load operation of the boilers have become available under conditions of unpreheated air and low-grade coal with volatiles as low as 15% and ash content as high as 30%. This principle shows good promise as an attractive new technology of combustion.

Distributions of temperature, velocity and viscosity at symmetrical axis of cylindrical plume in earth's mantle with variable viscosity

The parameters at the symmetrical axis of a cylindrical plume characterize the strength of this plume and provide a boundary condition which must be given to investigate the structure of a plume. For Newtonian fluid with a temperature-and pressure-dependence viscosity, an asymptotical solution of hydrodynamic equations at the symmetrical axis of the plume is found in the present paper. The temperature, upward velocity and viscosity at the symmetrical axis have been obtained as functions of depth, The calculated results have been given for two typical sets of Newtonian rheological parameters. The results obtained show that the temperature distribution along the symmetrical axis is nearly independent of the theological parameters. The upward velocity at the symmetrical axis, however, is strongly dependent on the rheological parameters.

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Bubble-to-Slug Transition for Low-Re Two-Phase Flows

A critical review on the mechanism and models on the bubble-to-slug transition of two-phase gas-liquid flows are presented in the present paper. It is shown that the most possible mechanism controlled the bubble-to-slug transition is the bubble coalescence. Focusing on the bubble-to-slug transition for the low-Re two-phase flow, a simple Monte Carlo method is used to simulate the influence of the initial bubble size on the bubble-to-slug transition. Some secondary factors, such as the liquid viscosity, the surface tension, and the relative slip between the two phases, are ignored in the present study. It is found that the locus of the dimensionless rate of collision is a universal curve. Based on this curve, it is determined that the bubble initial size can affect the phase distribution and flow pattern when its dimensionless value is in the range from 0.03 to 0.4. A simple relationship between the critical void fraction and the initial bubble size is proposed, which agrees very well with the experimental data.

Influence of Evaporation on Marangoni-Benard Instability

Evaporative convection and instability give rise to both scientific and technological interests. Practically, a number of the industrial applications such as thin-film evaporators, boiling technologies and heat pipes concern with the evaporation process of which through the vapor-liquid interface the heat and mass transfer occur. From a physical viewpoint, one of interesting questions is the mechanisms of convection instability in thin-liquid layers induced by the coupling of evaporation phenomenon and Marangoni effect at the mass exchanged interface. Classical theories, including Rayleigh’s and Pearson’s, have only successfully explained convection in a liquid layer heated from below without evaporation. However these theories are unable to explain the convection in an evaporating thin layer, especially liquid layer is cooled from below. In present paper, a new two-sided model is put forward rather than the one-sided model in previous works. In previous works, the vapor is treated as passive gas and dynamics of vapor has been ignored. In this case, the vapor liquid system can be described by one-sided model. In our two-sided model, the dynamics of vapor should be considered. Linear instability analysis of the Marangoni-Bénard convection in the two-layer system with an evaporation interface is performed. We define a new evaporating Biot number which is different from the Biot number in one-sided model and obtain the curves of critical Marangoni number versus wave number. In our theoretical results, the Biot number and the evaporating velocity play a major role in the stability of the vapor-liquid system.

Experimental Study on Subcooled Pool Boiling in Microgravity Utilizing