On-Board Experimental Study Of Bubble Thermocapillary Migration In A Recoverable Satellite

An on-board space experiment of bubble thermocapillary migration was performed in the Chinese 22nd recoverable satellite in 2005. Silicone oil of nominal viscosity 5cSt was used as the continuous phase in the experiment. Air bubbles were injected into the liquid in the same direction as the constant temperature gradient in the liquid. The velocities of bubbles were obtained by recording the paths of the bubbles. The results indicate that the scaled velocity of bubbles decreases with an increase of the Marangoni number extended to 9288, which agrees with the results of previous space experiments and numerical simulation. In addition, the interaction between two bubbles was also observed in the space experiment. The trajectories and the velocities of the bubbles were obtained. The two-bubble experiment results are also consistent with the theoretical analysis.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

A three-dimensional analytical solution of the microheater temperature before bubble nucleation

A three-dimensional analytical solution of the microheater temperature based on heat diffusion equation is developed and compared with experimental results. Dimensionless parameters are introduced to analyze the temperature rise time and the distribution under steady state. To study the microheater temperatures before bubble nucleation, a set of working fluids and microheaters are considered. It is shown that the dimensionless time xi(-)(0) required for the temperature rise from room to 95% of the steady state temperature is about 75, not dependent on working fluids and microheaters. Heat transfer to the surrounding liquid is mainly caused by conduction, not by convection and radiation mechanisms. The microheater length affects the surface temperature uniformity, while its width influences the steady temperatures significantly, yielding the transition from heterogeneous to homogeneous nucleation mechanism from square microheaters to narrow line microheaters. 

Periodic bubble emission and appearance of an ordered bubble sequence (train) during condensation in a single microchannel

Condensation of steam in a single microchannel, silicon test section was investigated visually at low flow rates. The microchannel was rectangular in cross-section with a depth of 30 pm, a width of 800 mu m and a length of 5.0 mm, covered with a Pyrex glass to allow for visualization of the bubble formation process. By varying the cooling rate during condensation of the saturated water vapor, it was possible to control the shape, size and frequency of the bubbles formed. At low cooling rates using only natural air convection from the ambient environment, the flow pattern in the microchannel consisted of a nearly stable elongated bubble attached upstream (near the inlet) that pinched off into a train of elliptical bubbles downstream of the elongated bubble. It was observed that these elliptical bubbles were emitted periodically from the tip of the elongated bubble at a high frequency, with smaller size than the channel width. The shape of the emitted bubbles underwent modifications shortly after their generation until finally becoming a stable vertical ellipse, maintaining its shape and size as it flowed downstream at a constant speed. These periodically emitted elliptical bubbles thus formed an ordered bubble sequence (train). At higher cooling rates using chilled water in a copper heat sink attached to the test section, the bubble formation frequency increased significantly while the bubble size decreased, all the while forming a perfect bubble train flowing downstream of the microchannel. The emitted bubbles in this case immediately formed into a circular shape without any further modification after their separation from the elongated bubble upstream. The present study suggests that a method for controlling the size and generation frequency of microbubbles could be so developed, which may be of interest for microfluidic applications. The breakup of the elongated bubble is caused by the large Weber number at the tip of the elongated bubble induced by the maximum vapor velocity at the centerline of the microchannel inside the elongated bubble and the smaller surface tension force of water at the tip of the elongated bubble.

Study on the kinetics of hydrate formation in a bubble column

Gas hydrate formation experiments were performed using methane in the presence of tetrahydrofuran (THF) in aqueous solution in a transparent bubble column in which a single pipe or a sintered plate was used to produce bubbles. The mole fraction of THF in aqueous solution was fixed at 6%. The hydrate formation kinetic behaviors on the surface of the rising bubble, the mechanical stability of hydrate shell formed on the surface of the bubble, the interactions among the bubbles with hydrate shell were observed and investigated morphologically. The rise velocities of individual bubbles with hydrate shells of different thickness and the consumption rates of methane gas were measured. A kinetic model was developed to correlate the experimentally measured gas consumption rate data. It was found that the hydrate formation rate on the surface of the moving bubble was high, but the formed hydrate shell was not very easy to be broken up. The bubbles with hydrate shells tended to agglomerate rather than merge into bigger bubble. This kind of characteristic of hydrate shell hindered the further formation of hydrate and led to the lower consumption rate of methane. The consumption rate of methane was found to increase with the decrease of temperature or increase of pressure. The increase of gas flux led to a linear increase in consumption rate of methane. It was demonstrated that the developed kinetic model could be used to correlate the consumption rate satisfyingly.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Multiphase technique improves mud-pulse velocity calculations

Measurement while drilling (MWD) has become a popular survey technology to monitor directional data, drilling data, formation evaluation data and safety data in the world. And closed loop drilling shows promise in recent years. Obviously, the method of tr

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, "turbulent spots" appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed neat sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces.

Detection and Characterization of Long-Pulse Low-Velocity Impact Damage in Plastic Bonded Explosives

Damage not only degrades the mechanical properties of explosives, but also influences the shock sensitivity, combustion and even detonation behavior of explosives. The study of impact damage is crucial in the vulnerability evaluation of explosives. A long-pulse low-velocity gas gun with a gas buffer was developed and used to induce impact damage in a hot pressed plastic bonded explosive. Various methods were used to detect and characterize the impact damage of the explosive. The microstructure was examined by use of polarized light microscopy. Fractal analysis of the micrographs was conducted by use of box counting method. The correlation between the fractal dimensions and microstructures was analyzed. Ultrasonic testing was conducted using a pulse through-transmission method to obtain the ultrasonic velocity and ultrasonic attenuation. Spectra analyses were carried out for recorded ultrasonic signals using fast Fourier transform. The correlations between the impact damage and ultrasonic parameters including ultrasonic velocities and attenuation coefficients were also analyzed. To quantitatively assess the impact induced explosive crystal fractures, particle size distribution analyses of explosive crystals were conducted by using a thorough etching technique, in which the explosives samples were soaked in a solution for enough time that the binder was totally removed. Impact induces a large extent of explosive crystal fractures and a large number of microcracks. The ultrasonic velocity decreases and attenuation coefficients increase with the presence of impact damage. Both ultrasonic parameters and fractal dimension can be used to quantitatively assess the impact damage of plastic bonded explosives.

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.

A modified thermal viscoplastic constitutive law involving the effect of temperature rise rate

At high temperature rise rate, the mechanical properties of 10 # steel were determined experimentally in a very wide range of temperature and strain rates. A new constitutive relationship was put forward, which can fit with the experimental results and describe various phenomena observed in our experiments. Meanwhile, some interesting characteristics about the temperature rise rate, strain and strain rate hardening and thermal softening are also shown in this paper. Finally, the reliability of the constitutive law and the correctness of the constitutive parameters were verified by comparing the calculation results with the experimental data.

Fourth order accurate compact scheme with group velocity control (GVC)

For solving complex flow field with multi-scale structure higher order accurate schemes are preferred. Among high order schemes the compact schemes have higher resolving efficiency. When the compact and upwind compact schemes are used to solve aerodynamic problems there are numerical oscillations near the shocks. The reason of oscillation production is because of non-uniform group velocity of wave packets in numerical solutions. For improvement of resolution of the shock a parameter function is introduced in compact scheme to control the group velocity. The newly developed method is simple. It has higher accuracy and less stencil of grid points.

Quantifying the Effects of Contact Duration, Loading Rate, and Approach Velocity on P-Selectin-PSGL-1 Interactions Using AFM

Kinetics and its regulation by extrinsic physical factors govern selectin-ligand interactions that mediate tethering and rolling of circulating cells on the vessel wall under hemodynamic forces. While the force regulation of off-rate for dissociation of selectin-ligand bonds has been extensively studied, much less is known about how transport impacts the on-rate for association of these bonds and their stability. We used atomic force microscopy (AFM) to quantify how the contact duration, loading rate, and approach velocity affected kinetic rates and strength of bonds of P-selectin interacting with P-selectin glycoprotein ligand I (PSGL-1). We found a saturable relationship between the contact time and the rupture force, a biphasic relationship between the adhesion probability and the retraction velocity, a piece-wise linear relationship between the rupture force and the logarithm of the loading rate, and a threshold relationship between the approach velocity and the rupture force. These results provide new insights into how physical factors regulate receptor-ligand interactions.

Turbulence in laminar premixed V-flames

Strong velocity fluctuations had been found in the laminar premixed V-flames. These velocity fluctuations are closely related to the chemical reaction. But the effects of the upstream combustible mixture velocity on the velocity fluctuations inside the flame are quite weak. The probability distribution function (PDF) of the velocity in the centre region of the flame appears "flat top" shaped. By analyzing the experiment results the flame-flow interactions are found to affect the flame not only at large scale in the flow field but also at small scale inside the flame. These effects will give rise to flame generated small scale turbulences.

Scheme Construction with Numerical Flux Residual Correction (NFRC) and Group Velocity Control (GVC)

For simulating multi-scale complex flow fields like turbulent flows, the high order accurate schemes are preferred. In this paper, a scheme construction with numerical flux residual correction (NFRC) is presented. Any order accurate difference approximation can be obtained with the NFRC. To improve the resolution of the shock, the constructed schemes are modified with group velocity control (GVC) and weighted group velocity control (WGVC). The method of scheme construction is simple, and it is used to solve practical problems.