Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.

Thermocapillary Flows In Liquid Bridges Of Molten Tin With Small Aspect Ratios

Numerical simulations were conducted to study thermocapillary flows in short half-zone liquid bridges of molten tin with Prandtl number Pr = 0.009, under ramped temperature difference. The spatio-temporal structures in the thermocapillary flows in short half-zone liquid bridges with aspect ratios As = 0.6, 0.8, and 1.0 were investigated. The first critical Marangoni numbers were compared with those predicted by linear stability analyses (LSA). The second critical Marangoni numbers for As = 0.6 and 0.8 were found to be larger than that for As = 1.0. The time evolutions of the thermocapillary flows exhibited unusual features such as a change in the azimuthal wave number during the three-dimensional stationary (non-oscillating) flow regime, a change in the oscillation mode during the three-dimensional oscillatory flow regime, and the decreasing and then increasing of amplitudes in a single oscillation mode. The effects of the ramping rate of the temperature difference on the flow modes and critical conditions were studied as well. In this paper, the experimental observability of the critical conditions was also discussed. (C) 2008 Elsevier Inc. All rights reserved.

Experimental study of thermocapillary convection on a liquid bridge consisting of fluid with low Prandtl number in a floating half-zone

A device of mercury liquid bridge of floating half-zone is designed to experimentally explore thermocapillary convection and its instability of a low Prandtl number liquid. Noncontacted diagnostic techniques were developed to monitor surface flow and surface deformation. The surface flow and the influence of a growing surface film (or skin) on the flow were observed. It is shown that the film is a key factor in changing the behavior associated with the thermocapillary convection. The experiment indicates that the critical Marangoni number should be much higher than that expected by the numerical simulation. The condition and process of surface film growth are discussed. The surface oscillation of the mercury bridge wrapped with ''dirt-film'' was observed, and the characteristics and the frequency associated with this oscillation are given.

Some characteristics of low-speed streaks under sheared air-water interfaces

The characteristics of low-speed fluid streaks occurring under sheared air-water interfaces were examined by means of hydrogen bubble visualization technique. A critical shear condition under which the streaky structure first appears was determined to be u(tau) approximate to 0.19 cm/s. The mean spanwise streak spacing increases with distance from the water surface owing to merging and bursting processes, and a linear relationship describing variation of non-dimensional spacing <(+)over bar> versus y(+) was found essentially independent of shear stress on the interface. Values of <(+)over bar>, however, are remarkably smaller than their counterparts in the near-wall region of turbulent boundary layers. Though low-speed streaks occur randomly in time and space, the streak spacing exhibits a lognormal probability distribution behavior. A tentative explanation concerning the formation of streaky structure is suggested, and the fact that <(+)over bar> takes rather smaller values than that in wall turbulence is briefly discussed.

Experimental Investigation On The Holdup Distribution Of Oil-Water Two-Phase Flow In Horizontal Parallel Tubes

An experimental investigation was conducted to study the holdup distribution of oil and water two-phase flow in two parallel tubes with unequal tube diameter. Tests were performed using white oil (of viscosity 52 mPa s and density 860 kg/m(3)) and tap water as liquid phases at room temperature and atmospheric outlet pressure. Measurements were taken of water flow rates from 0.5 to 12.5 m(3)/h and input oil volume fractions from 3 to 94 %. Results showed that there were different flow pattern maps between the run and bypass tubes when oil-water two-phase flow is found in the parallel tubes. At low input fluid flow rates, a large deviation could be found on the average oil holdup between the bypass and the run tubes. However, with increased input oil fraction at constant water flow rate, the holdup at the bypass tube became close to that at the run tube. Furthermore, experimental data showed that there was no significant variation in flow pattern and holdup between the run and main tubes. In order to calculate the holdup in the form of segregated flow, the drift flux model has been used here.

Computer simulation for effect of Pr(III) on Ca(II) speciation in human interstitial fluid

A multiphase model of metal ion speciation in human interstitial fluid was constructed and the effect of Pr(III) on Ca(II) speciation was studied. Results show that free Ca2+, [Ca(HCO3)], and [Ca(Lac)] are the main species of Ca(II). Because of the competition of Pr(III) for ligands with Ca(II), the percentages of free Ca2+, [Ca(Lac)], and [Ca(His)(Thr)H-3] increase gradually and the percentages of CaHPO4(aq) and [Ca(Cit)(His)H-2] decrease gradually with the increase in the total concentration of Pr(III). However, the percentages of [Ca(HCO3)] and CaCO3(aq) first increase and then begin to decrease when the total concentration of Pr(III) exceeds 6.070 x 10(-4) M.

Computer simulation of insoluble Pr(III) speciation in human interstitial fluid

A multi-phase model of Pr(III) speciation in human interstitial fluid was constructed and insoluble Pr(III) speciation was studied. When the total concentration of Pr(III) is below 8.401E-10 mol/L, soluble Pr(III) species are main species. With rising the total concentration of Pr(III), Pr(III) is firstly bound to phosphate to form precipitate of PrPO4, then bound to carbonate and another precipitate of Pr-2(CO3)(3) was obtained. When the total concentration is between 1.583E-9 mol/L and 4.000E-3 mol/L, the insoluble species are predominant Pr(III) species.

Damping of Vertically Excited Surface Wave in Weakly Viscous Fluid

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.

Two bifurcation transitions of the floating half zone convection in a fat liquid bridge of larger Pr

The transient process of the thermocapillary convection was obtained for the large Pu floating half zone by using the method of three-dimensional and unsteady numerical simulation. The convection transits directly from steady and axisymmetric state to oscillatory flow for slender liquid bridge, and transits first from steady and axisymmetric convection to the steady and non-axisymmetric convection, then, secondly to the oscillatory convection for the fatter liquid bridge. This result implies that the volume of liquid bridge is not only a sensitive critical parameter for the onset of oscillation, but also relates to the new mechanism for the onset of instability in the floating half zone convection even in case of large Prandtl number fluid.

The Hydrothermal Wave Of Large-Prandtl-Number Fluid In A Shallow Cavity

The hydrothermal wave was investigated numerically for large-Prandtl-number fluid (Pr = 105.6) in a shallow cavity with different heated sidewalls. The traveling wave appears and propagates in the direction opposite to the surface flow (upstream) in the case of zero gravity when the applied temperature difference grows and over the critical value. The phase relationships of the disturbed velocity, temperature and pressure demonstrate that the traveling wave is driven by the disturbed temperature, which is named hydrothermal wave. The hydrothermal wave is so weak that the oscillatory flow field and temperature distribution can hardly be observed in the liquid layer. The exciting mechanism of hydrothermal wave is analyzed and discussed in the present paper.

Progress In Modeling Of Fluid Flows In Crystal Growth Processes

Modeling of fluid flows in crystal growth processes has become an important research area in theoretical and applied mechanics. Most crystal growth processes involve fluid flows, such as flows in the melt, solution or vapor. Theoretical modeling has played an important role in developing technologies used for growing semiconductor crystals for high performance electronic and optoelectronic devices. The application of devices requires large diameter crystals with a high degree of crystallographic perfection, low defect density and uniform dopant distribution. In this article, the flow models developed in modeling of the crystal growth processes such as Czochralski, ammonothermal and physical vapor transport methods are reviewed. In the Czochralski growth modeling, the flow models for thermocapillary flow, turbulent flow and MHD flow have been developed. In the ammonothermal growth modeling, the buoyancy and porous media flow models have been developed based on a single-domain and continuum approach for the composite fluid-porous layer systems. In the physical vapor transport growth modeling, the Stefan flow model has been proposed based on the flow-kinetics theory for the vapor growth. In addition, perspectives for future studies on crystal growth modeling are proposed. (c) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

An Approximate Method For Calculation Of Fluid Force And Response Of A Circular Cylinder At Lock-In

In this paper, equations calculating lift force of a rigid circular cyclinder at lock-in uniform flow are deduced in detail. Besides, equations calculating the lift force on a long flexible circular cyclinder at lock-in are deduced based on mode analysis of a multi-degree freedom system. The simplified forms of these equations are also given. Furthermore, an approximate method to predict the forces and response of rigid circular cyclinders and long flexible circular cyclinders at lock-in is introduced in the case of low mass-damping ratio. A method to eliminate one deficiency of these equations is introduced. Comparison with experimental results show the effectiveness of this approximate method.

The Head-On Colliding Process Of Binary Liquid Droplets At Low Velocity: High-Speed Photography Experiments And Modeling

The experimental and theoretical studies are reported in this paper for the head-on collisions of a liquid droplet with another of the same fluid resting on a solid substrate. The droplet on the hydrophobic polydimethylsiloxane (PDMS) substrate remains in a shape of an approximately spherical segment and is isometric to an incoming droplet. The colliding process of the binary droplets was recorded with high-speed photography. Head-on collisions saw four different types of response in our experiments: complete rebound, coalescence, partial rebound With conglutination, and coalescence accompanied by conglutination. For a complete rebound, both droplets exhibited remarkable elasticity and the contact time of the two colliding droplets was found to be in the range of 10-20 ms. With both droplets approximately considered as elastic bodies, Hertz contact theory was introduced to estimate the contact time for the complete rebound case. The estimated result Was found to be on the same order of magnitude as the experimental data, which indicates that the present model is reasonable. (C) 2008 Elsevier Inc. All rights reserved.