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.

Numerical investigation of thermocapillary migration of the drop for large Marangoni numbers

An axisymmetric model is adopted to simulate the problem of unsteady drop thermocapillary motion for large Marangoni numbers. Front tracking methods are used in the investigation. It is found that the non-dimensional drop migration velocity will decrease with increasing Marangoni number. This agrees well with the experimental results obtained from the 4th Shen-Zhou space ship. In the meanwhile, this is also the first time for numerical simulations to verify the experimental phenomenon under large Marangoni numbers.

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.

Internal resonance of an L-shaped beam with a limit stop .1. Free vibration

Real-life structures often possess piecewise stiffness because of clearances or interference between subassemblies. Such an aspect can alter a system's fundamental free vibration response and leads to complex mode interaction. The free vibration behaviour of an L-shaped beam with a limit stop is analyzed by using the frequency response function and the incremental harmonic balance method. The presence of multiple internal resonances, which involve interactions among the first five modes and are extremely complex, have been discovered by including higher harmonics in the analysis. The results show that mode interaction may occur if the higher harmonics of a vibration mode are close to the natural frequency of a higher mode. The conditions for the existence of internal resonance are explored, and it is shown that a prerequisite is the presence of bifurcation points in the form of intersecting backbone curves. A method to compute such intersections by using only one harmonic in the free vibration solution is proposed. (C) 1996 Academic Press Limited

Internal resonance of an L-shaped beam with a limit stop .2. Forced vibration

A limit stop is placed at the elbow of an L-shaped beam whose linear natural frequencies are nearly commensurable. As a result of this hardening device the non-linear system exhibits multiple internal resonances, which involve various degree of coupling between the first five modes of the beam in free vibration. A point load is so placed as to excite several modes and the resulting forced vibration is examined. In the undamped case, three in-phase and two out-of-phase solution branches have been found. The resonance curve is extremely complicated, with multiple branches and interactions between the first four modes. The amplitudes of the higher harmonics are highly influenced by damping, the presence of which can effectively attenuate internal resonances. Consequently parts of the resonance curve may be eliminated, with the resulting response comprising different distinctive branches. (C) 1996 Academic Press Limited

A solution to the upward migration of an intrusion

The evolution of the upward migration of the magma is a nonlinear and unstable problem in mathematics. It is difficult to solve it. And using the numerical method, the solution is relatively tedious and time-consuming. This paper introduces a method of the instantaneous point source to solve the linear and unstable heat conduction equation during the infinite period of time instead of the solution of the nonlinear and unstable heat conduction equation. The results obtained by this method coincide with those by the numerical method, meaning that this method offers a simple way to solve the nonlinear and unstable heat conduction equation.

The effect of variable currents on internal solitary waves

The effect of variable currents on internal solitary waves is described within the context of a variable coefficient Korteweg-de Vries (KdV) equation, and the approximate slowly varying, solitary-wave solution of this equation. The general theory which leads to the variable coefficient KdV equation is described; a derivation for the special case when the solitary wave and the current are aligned in the same direction is given in the Appendix. Using further simplifications and approximations, a number of analytical expressions are obtained for the variation in the solitary wave amplitude resulting from variable shear in the basic current or from when the basic current is a depth-independent flow which is a simple representation of a geostrophic current, tidal flow or inertial wave.

Visualization of turbulence and internal waves in a stratified fluid

A new method for measuring the density, temperature and velocity of N2 gas flow by laser induced biacetyl phosphorescence is proposed. The characteristics of the laser induced phosphorescence of biacetyl mixed with N2 are investigated both in static gas and in one-dimensional flow along a pipe with constant cross section. The theoretical and experimental investigations show that the temperature and density of N2 gas flow could be measured by observing the phosphorescence lifetime and initial intensity of biacetyl triplet (3Au) respectively. The velocity could be measured by observing the time-of-flight of the phosphorescent gas after pulsed laser excitation. The prospect of this method is also discussed.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Two Bubbles Thermocapillary Migration Experimental System

An experimental investigation will be performed on the thermocapillary motion of two bubbles in Chinese return-satellite. The experiment will study the migration process of bubble caused by thermocapillary effect in microgravity environment, and their interaction between two bubbles. The bubble is driven by the thermocapillary stress on the surface on account on the variation of the surface tension with temperature. The interaction between two bubbles becomes significant as the separation distance between them is reduced drastically so that the bubble interaction has to be considered. Recently, the problem has been discussed on the method of successive reflections, and accurate migration velocities of two arbitrarily oriented bubbles were derived for the limit of small Marangoni and Reynolds numbers. Numerical results for the migration of the two bubbles show that the interaction between two bubbles has significant influence on their thermocapillary migration velocities with a bubble approaching another. However, there is a lack of experimental validate for the theoretic results. Now the experimental facility is designed for experimenting time after time. A cone-shaped top cover is used to expel bubble from the cell after experiment. But, the cone-shaped top cover can cause temperature uniformity on horizontal plane in whole cell. Therefore, a metal board with multi-holes is fixed under the top cover. The board is able to let the temperature distribution on the board uniform because of their high heat conductivity, and the bubble can pass through it. In the system two bubbles are injected into the test cell respectively by two sets of cylinder. And the bubbles sizes are controlled by two sets of step-by-step motor. It is very important problem that bubble can be divorced from the injecting mouth in microgravity environment. Thus, other two sets of device for injecting mother liquid were used to push bubble. The working principle of injecting mother liquid is to utilize pressure difference directly between test cell and reservoir

Influence of Wind and Grain Size on Migration of Asymmetric Sand Waves

Large parts of shallow seas are covered by regular seabed patterns and sand wave is one kind of these patterns. The instability of the sedimentary structures may hazard pipelines and the foundations of offshore structures. In the last decade or so, it's a focus for engineers to investigate the movement mechanism of sand waves. Previous theoretical studies of the subject have developed a general model to predict the growth and migration of sand waves, which is based on the two-dimensional vertical shallow water equations and the bed-form deformation equations. Although the relation between wave-current flow and sand bed deformation has been established, the topographic influence has not been considered in the model. In this paper some special patterns, which are asymmetric and close to the reality, are represent as the perturbed seabed and the evolution of sand waves is calculated. The combination of a steady flow induced by wind and a sinusoidal tidal flow is considered as the basic flow. Finally the relations of some parameters (grain size, etc.) and sand waves' growth and migration are discussed, and the growth rate and migration speeds of asymmetric sand waves are carried out.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

The induced flow field by internal solitary wave and its action on cylindrical piles in the stratified ocean

The induced flow fields by internal solitary waves and its actions on cylindrical piles in density stratified ocean with a basic density profile and a basic velocity profile are investigated. Some results, such as the time evolution of flow fields and hydrodynamic forces on the piles are yielded both by theoretical analysis and numerical calculation for general and specific cases. Several kinds of ambient sea conditions of the South China Sea are specified for numerical simulation. Moreover, the effects of relative density difference, depth ratio and wave steepness on maximal total force and total torque are analyzed.

An auto-focusing method for imaging ellipsometry system

An auto-focusing method based on the image brightness gradient sharpness function is presented for imaging ellipsometry system, in which the image plane of the thin-film specimen is not perpendicular to the optical axis. The clear image of a specimen with large area is obtained by moving the imaging sensor in optical axis direction and around its sensitive surface centre successively. The experimental results demonstrate its feasibility.