Pool Boiling in Microgravity: Recent Results and Perspectives for the Project DEPA-SJ10

Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments have also been performed both in normal gravity and in short-term microgravity in the Drop Tower Beijing. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. In the first case, slight enhancement of heat transfer is observed in microgravity, while diminution is evident for high heat flux in the second one. Lateral motions of bubbles on the heaters are observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drives it to detach from the heaters. The Marangoni effect on the bubble behavior is also discussed. The perspectives for a new project DEPA-SJ10, which has been planned to be flown aboard the Chinese recoverable satellite SJ-10 in the future, are also presented.

Lateral motion and departure of vapor bubbles in nucleate pool boiling on thin wires in microgravity

A space experiment on bubble behavior and heat transfer in subcooled pool boiling phenomenon has been performed utilizing the temperature-controlled pool boiling (TCPB) device both in normal gravity in the laboratory and in microgravity aboard the 22(nd) Chinese recoverable satellite. The fluid is R113 at 0.1 MPa and subcooled by 26 degrees C nominally. A thin platinum wire of 60 mu m in diameter and 30mm in length is simultaneously used as heater and thermometer. Only the lateral motion and the departure of discrete vapor bubbles in nucleate pool boiling are reported and analyzed in the present paper. A scale analysis on the Marangoni convection surrounding a bubble in the process of subcooled nucleate pool boiling leads to formulas of the characteristic velocity of the lateral motion and its observability. The predictions consist with the experimental observations. Considering the Marangoni effect, a new qualitative model is proposed to reveal the mechanism underlying the bubble departure processes and a quantitative agreement can also be acquired.

Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

Two-Phase Flow In Anode Flow Field Of A Small Direct Methanol Fuel Cell In Different Gravities

An in-situ visualization of two-phase flow inside anode flow bed of a small liquid fed direct methanol fuel cells in normal and reduced gravity has been conducted in a drop tower. The anode flow bed consists of 11 parallel straight channels. The length, width and depth of single channel, which had rectangular cross section, are 48.0, 2.5 and 2.0 mm, respectively. The rib width was 2.0 mm. The experimental results indicated that when the fuel cell orientation is vertical, two-phase flow pattern in anode channels can evolve from bubbly flow in normal gravity into slug flow in microgravity. The size of bubbles in the reduced gravity is also bigger. In microgravity, the bubbles rising speed in vertical channels is obviously slower than that in normal gravity. When the fuel cell orientation is horizontal, the slug flow in the reduced gravity has almost the same characteristic with that in normal gravity. It implies that the effect of gravity on two-phase flow is small and the bubbles removal is governed by viscous drag. When the gas slugs or gas columns occupy channels, the performance of liquid fed direct methanol fuel cells is failing rapidly. It infers that in long-term microgravity, flow bed and operating condition should be optimized to avoid concentration polarization of fuel cells.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

Subcooled pool boiling on thin wire in microgravity

A new set of experimental data of subcooled pool boiling on a thin wire in rnicrogravity aboard the 22nd Chinese recoverable satellite is reported in the present paper. The temperature-control led heating method is used. The results of the experiments in normal gravity before and after the flight experiment are also presented, and compared with those in microgravity. The working fluid is degassed R113 at 0.1 MPa and subcooled by 26 degrees C nominally. A thin platinum wire of 60 mu m in diameter and 30 mm in length is simultaneously used as heater and thermometer. It is found that the heat transfer of nucleate pool boiling is slightly enhanced in microgravity comparing with those in normal gravity. It is also found that the correlation of Lienhard and Dhir can predict the CHF with good agreement, although the range of the dimensionless radius is extended by three or more decades above the originally set limit. Three critical bubble diameters are observed in microgravity, which divide the observed vapor bubbles into four regimes with different sizes. Considering the Marangoni effect, a qualitative model is proposed to reveal the mechanism underlying the bubble departure processes, and a quantitative agreement can also be acquired.

微重力准稳态池沸腾中的气泡动力学研究

本文实验研究了微重力条件下的准稳态核态池沸腾现象中的气泡动力学特征及其对传热特性的影响,发现在低过冷度沸腾中,聚合汽泡表面的强烈振荡促进了整个加热面上的核化过程与核态沸腾传热;而高过冷沸腾中,聚合气泡在表面张力作用下呈球状,难以覆盖整个加热表面,导致核态沸腾向膜态沸腾的过度过程表现为"核态沸腾+局部的干斑扩展"现象,相应的传热曲线没有明显转折.微重力条件下池沸腾临界热流随过冷度和压力的增加而升高,与地面结果定性相符.

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included groundbased tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.

Physical Modeling And Parametric Study On Two-Degree-Of-Freedom Viv Of A Cylinder Near Rigid Wall

Unlike most previous studies on the transverse vortex-induced vibration(VIV) of a cylinder mainly under the wallfree condition (Williamson & Govardhan,2004),this paper experimentally investigates the vortex-induced vibration of a cylinder with two degrees of freedom near a rigid wall exposed to steady flow.The amplitude and frequency responses of the cylinder are discussed.The lee wake flow patterns of the cylinder undergoing VIV were visualized by employing the hydrogen bubble technique.The effects of the gap-to-diameter ratio (e0/D) and the mass ratio on the vibration amplitude and frequency are analyzed.Comparisons of VIV response of the cylinder are made between one degree (only transverse) and two degrees of freedom (streamwise and transverse) and those between the present study and previous ones.The experimental observation indicates that there are two types of streamwise vibration,i.e.the first streamwise vibration (FSV) with small amplitude and the second streamwise vibration (SSV) which coexists with transverse vibration.The vortex shedding pattem for the FSV is approximately symmetric and that for the SSV is alternate.The first streamwise vibration tends to disappear with the decrease of e0/D.For the case of large gap-to-diameter ratios (e.g.e0/D = 0.54～1.58),the maximum amplitudes of the second streamwise vibration and transverse one increase with the increasing gapto-diameter ratio.But for the case of small gap-to-diameter ratios (e.g.e0/D = 0.16,0.23),the vibration amplitude of the cylinder increases slowly at the initial stage (i.e.at small reduced velocity V,),and across the maximum amplitude it decreases quickly at the last stage (i.e.at large Vr).Within the range ofthe examined small mass ratio (m＜4),both streamwise and transverse vibration amplitude of the cylinder decrease with the increase of mass ratio for the fixed value of V,.The vibration range (in terms of Vr ) tends to widen with the decrease of the mass ratio.In the second streamwise vibration region,the vibration frequency of the cylinder with a small mass ratio (e.g.mx = 1.44) undergoes a jump at a certain Vr,.The maximum amplitudes of the transverse vibration for two-degree-of-freedom case is larger than that for one-degree-of-freedom case,but the transverse vibration frequency of the cylinder with two degrees of freedom is lower than that with one degree of freedom (transverse).

Two-Dimensional Hillslope Scale Soil Erosion Model

On a hillslope, overland flow first generates sheet erosion and then, with increasing flux, it causes rill erosion. Sheet erosion (interrill erosion) and rill erosion are commonly observed to coexist on hillslopes. Great differences exist between both the intensities and incidences of rill and interrill erosion. In this paper, a two-dimensional rill and interrill erosion model is developed to simulate the details of the soil erosion process on hillslopes. The hillslope is treated as a combination of a two-dimensional interrill area and a one-dimensional rill. The rill process, the interrill process, and the joint occurrence of rill and interrill areas are modeled, respectively. Thus, the process of sheet flow replenishing rill flow with water and sediment can be simulated in detail, which may possibly render more truthful results for rill erosion. The model was verified with two sets of data and the results seem good. Using this model, the characteristics of soil erosion on hillslopes are investigated. Study results indicate that (1) the proposed model is capable of describing the complex process of interrill and rill erosion on hillslopes; (2) the spatial distribution of erosion is simulated on a simplified two-dimensional hillslope, which shows that the distribution of interrill erosion may contribute to rill development; and (3) the quantity of soil eroded increases rapidly with the slope gradient, then declines, and a critical slope gradient exists, which is about 15-20 degrees for the accumulated erosion amount.

Dissipative particle dynamics simulation of wettability alternation phenomena in the chemical flooding process

Wettability alternation phenomena is considered one of the most important enhanced oil recovery (EOR) mechanisms in the chemical flooding process and induced by the adsorption of surfactant on the rock surface. These phenomena are studied by a mesoscopic method named as dissipative particle dynamics (DPD). Both the alteration phenomena of water-wet to oil-wet and that of oil-wet to water-wet are simulated based on reasonable definition of interaction parameters between beads. The wetting hysteresis phenomenon and the process of oil-drops detachment from rock surfaces with different wettability are simulated by adding long-range external forces on the fluid particles. The simulation results show that, the oil drop is liable to spread on the oil-wetting surface and move in the form of liquid film flow, whereas it is likely to move as a whole on the water-wetting surface. There are the same phenomena occuring in wettability-alternated cases. The results also show that DPD method provides a feasible approach to the problems of seepage flow with physicochemical phenomena and can be used to study the mechanism of EOR of chemical flooding.

Laser-driven coherent betatron oscillation in a laser-wakefield cavity

The origin of beam disparity in emittance and betatron oscillation orbits, in and out of the polarization plane of the drive laser of laser-plasma accelerators, is explained in terms of betatron oscillations driven by the laser field. As trapped electrons accelerate, they move forward and interact with the laser pulse. For the bubble regime, a simple model is presented to describe this interaction in terms of a harmonic oscillator with a driving force from the laser and a restoring force from the plasma wake field. The resulting beam oscillations in the polarization plane, with period approximately the wavelength of the driving laser, increase emittance in that plane and cause microbunching of the beam. These effects are observed directly in 3D particle-in-cell simulations.

Pondermotive shift of photoelectron spectra in intense laser field

In a recent experimental work on the excess photon detachment (EPD) of H- ions [Phys. Rev. Lett. 87 (2001) 243001] it has been found that the ponderomotive shift of each EPD peak increases with the order of the EPD channel. By using a nonperturbative quantum scattering theory, we obtain the kinetic energy spectra for the differential detachment rate along the laser polarization for several laser intensities. It is demonstrated that higher order EPD peaks are produced mainly at relatively higher laser intensities. By calculating the overall EPD spectra with varying laser intensities, it is found that the ponderomotive shift of each EPD peak increases with the order of the EPD channel. Our calculations are in good agreement with the experimental observation. It is found that different EPD channels occur mainly when the laser field reaches some values, thus the intensity distribution of the laser field is responsible for the varying ponderomotive shifts.

Structure-change-dependent transmittance of PtOx thin film in super-resolution near-field structure: Diffraction effect analysis

To study working mechanism of super-resolution near-field structure (super-RENS) optical disk from a far-field optics view is very necessary because of the actual far-field writing/readout process in the optical disk system. A Gaussian diffraction model based on Fresnel-Kirchhoff diffraction theory of PtOx-type super-RENS has been set up in this Letter. The relationship between micro-structural deformation (change of bubble structure and refractive index profile) with far-field optical response of PtOx thin film has been studied with it in detail. The simulation results are in good agreement with the experimental results reported in literatures with a designed configuration. These results may provide more quantitative information for better understanding of the working mechanism of metal-oxide-type super-RENS. (c) 2007 Elsevier B.V. All rights reserved.

Study on readout durability of super-RENS disk

Characteristics essential for the readout durability of a superresolution near-field structure (super-RENS) disk are studied experimentally by using a home-built optical measuring setup and atomic force microscope, based on a simplified PtOx super-RENS disk. The experimental results show that for a super-RENS disk with constant structure and materials, readout signals including transmittance and reflectance vary with changes in bubble shape and size, indicating that the readout durability of the disk has a strong dependence on bubble stability, which is closely related to the thickness of the cover layer, the recording power and readout power, and the mechanical properties of the dielectric layer. Based on our experimental results, the main direction for improving readout durability is also proposed. (c) 2008 Optical Society of America.