摄动有限体积算法及其在两相流中应用

采用一种全新的摄动有限体积(PFV)算法和水平集(Level Set)技术对液液两相系统中液滴坠落进行数值模拟,数值结果表明,PFV新算法具有节点少、精度高,效率高,编程方便等优点,能成功模拟液液两相流动,为两相流动数值模拟提供了一种新的途径.

Investigation on mechanism of critical cavitating flow in liquid jet pumps under operating limits

The critical cavitating flow in liquid jet pumps under operating limits is investigated in this paper. Measurements on the axial pressure distribution along the wall of jet pumps indicate that two-phase critical flow occurs in the throat pipe under operating limits. The entrained flow rate and the distribution of the wall pressure upstream lowest pressure section does not change when the outlet pressure is lower than a critical value. A liquid-vapor mixing shockwave is also observed under operating limits. The wave front moves back and forth in low frequency around the position of the lowest pressure. With the measured axial wall pressures, the Mach number of the two-phase cavitating flow is calculated. It's found that the maximum Mach number is very close to I under operating limits. Further analysis infers a cross-section where Mach number approaches to I near the wave front. Thus, the liquid-vapor mixture velocity should reach the local sound velocity and resulting in the occurrence of operating limits.

A Simple Model For Predicting The Void Fraction Of Gas/Non-Newtonian Fluid Intermittent Flows In Upward Inclined Pipes

In this work, a simple correlation, which incorporates the mixture velocity, drift velocity, and the correction factor of Farooqi and Richardson, was proposed to predict the void fraction of gas/non-Newtonian intermittent flow in upward inclined pipes. The correlation was based on 352 data points covering a wide range of flow rates for different CMC solutions at diverse angles. A good agreement was obtained between the predicted and experimental results. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase intermittent flows.

节流器内液-固两相流固体颗粒冲蚀数值模拟

建立了考虑颗粒碰撞的颗粒冲蚀计算模型,该数学模型包括:在Eulerian坐标系下求解连续相流场;在Lagrangian坐标系下运用离散颗粒硬球模型求解颗粒碰撞;应用半实验关联式求解颗粒冲蚀速率。对水力加砂压裂施工中节流器内液-固两相流的固体颗粒运动和冲蚀特性进行了数值模拟。计算结果表明,固体颗粒密集于节流器入口到出口的一段狭长区域内,冲蚀速率随流体速度呈指数性变化。颗粒直径越大,冲蚀速率也越大。节流器内冲蚀最严重的位置发生在距离节流器出口上边缘10mm以内的局部区域。

Ultrafast dynamics in ZnO thin films irradiated by femtosecond lasers

The damage morphologies, threshold fluences in ZnO films were studied with femtosecond laser pulses. Time-resolved reflectivity and transmissivity have been measured by the pump-probe technique at different pump fluences and wavelengths. The results indicate that two-phase transition is the dominant damage mechanism, which is similar to that in narrow band gap semiconductors. The estimated energy loss rate of conduction electrons is 1.5 eV/ps. (c) 2005 Elsevier Ltd. All rights reserved.

Single-photon source in strong photon-atom interaction regime in quasiperiodic photonic crystals

The antibunching properties of the fluorescence from a two-level ideal system in a 12-fold quasiperiodic photonic crystal are investigated based on the calculated local density of states. We found that the antibunching phenomenon of the fluorescence from two-level ideal systems could be significantly changed by varying their positions, i.e., perfect antibunching and antibunching with damped Rabi oscillation phenomenon occurred in different positions and at different frequencies in photonic crystals as a result of the large differences in the local density of states. This study revealed that the multi-level coherence of fluorescence from a two-level ideal system could be manipulated by controlling the position of the two-level ideal system in photonic crystals and the emission frequency in the photonic band structure. Copyright (C) EPLA, 2008

Uniform mems chip temperatures in the nucleate boiling heat transfer region by selecting suitable, medium boiling number range

The not only lower but also uniform MEMS chip temperatures can he reached by selecting suitable boiling number range that ensures the nucleate boiling heat transfer. In this article, boiling heat transfer experiments in 10 silicon triangular microchannels with the hydraulic diameter of 55.4 mu m were performed using acetone as the working fluid, having the inlet liquid temperatures of 24-40 degrees C, mass fluxes of 96-360 kg/m(2)s, heat fluxes of 140-420 kW/m(2), and exit vapor mass qualities of 0.28-0.70. The above data range correspond to the boiling number from 1.574 x 10(-3) to 3.219 x 10(-3) and ensure the perfect nucleate boiling heat transfer region, providing a very uniform chip temperature distribution in both streamline and transverse directions. The boiling heat transfer coefficients determined by the infrared radiator image system were found to he dependent on the heat Axes only, not dependent on the mass Axes and the vapor mass qualities covering the above data range. The high-speed flow visualization shows that the periodic flow patterns take place inside the microchannel in the time scale of milliseconds, consisting of liquid refilling stage, bubble nucleation, growth and coalescence stage, and transient liquid film evaporation stage in a full cycle. The paired or triplet bubble nucleation sites can occur in the microchannel corners anywhere along the flow direction, accounting for the nucleate boiling heat transfer mode. The periodic boiling process is similar to a series of bubble nucleation, growth, and departure followed by the liquid refilling in a single cavity for the pool boiling situation. The chip temperature difference across the whole two-phase area is found to he small in a couple of degrees, providing a better thermal management scheme for the high heat flux electronic components. Chen's [11 widely accepted correlation for macrochannels and Bao et al.'s [21 correlation obtained in a copper capillary tube with the inside diameter of 1.95 mm using R11 and HCFC123 as working fluids can predict the present experimental data with accepted accuracy. Other correlations fail to predict the correct heat transfer coefficient trends. New heat transfer correlations are also recommended.

Multi-channel effect of condensation flow in a micro triple-channel condenser

Multi-channel effect is important to understand transport phenomenon in phase change systems with parallel channels. In this paper, visualization studies were performed to study the multi-channel effect in a silicon triple-channel condenser with an aspect ratio of 0.04. Saturated water vapor was pumped into the microcondenser, which was horizontally positioned. The condenser was cooled by the air natural convention heat transfer in the air environment. Flow patterns are either the annular flow at high inlet vapor pressures, or a quasi-stable elongated bubble at the microchannel upstream followed by a detaching or detached miniature bubble at smaller inlet vapor pressures. The downstream miniature bubble was detached from the elongated bubble tip induced by the maximum Weber number there. It is observed that either a single vapor thread or dual vapor threads are at the front of the elongated bubble. A miniature bubble is fully formed by breaking up the vapor thread or threads. The transient vapor thread formation and breakup process is exactly symmetry against the centerline of the center channel. In side channels, the Marangoni effect induced by the small temperature variation over the channel width direction causes the vapor thread formation and breakup process deviating from the side channel centerline and approaching the center channel. The Marangoni effect further forces the detached bubble to rotate and approach the center channel, because the center channel always has higher temperatures, indicating the multi-channel effect. 

Transient flow patterns and bubble slug lengths in parallel microchannels with oxygen gas bubbles produced by catalytic chemical reactions

Transient flow patterns and bubble slug lengths were investigated with oxygen gas (O-2) bubbles produced by catalytic chemical reactions using a high speed camera bonded with a microscope. The microreactor consists of an inlet liquid plenum, nine parallel rectangular microchannels followed by a micronozzle, using the MEMS fabrication technique. The etched surface was deposited by the thin platinum film, which is acted as the catalyst. Experiments were performed with the inlet mass concentration of the hydrogen peroxide from 50% to 90% and the pressure drop across the silicon chip from 2.5 to 20.0 kPa. The silicon chip is directly exposed in the environment thus the heat released via the catalytic chemical reactions is dissipated into the environment and the experiment was performed at the room temperature level. It is found that the two-phase flow with the catalytic chemical reactions display the cyclic behavior. A full cycle consists of a short fresh liquid refilling stage, a liquid decomposition stage followed by the bubble slug flow stage. At the beginning of the bubble slug flow stage, the liquid slug number reaches maximum, while at the end of the bubble slug flow stage the liquid slugs are quickly flushed out of the microchannels. Two or three large bubbles are observed in the inlet liquid plenum, affecting the two-phase distributions in microchannels. The bubble slug lengths, cycle periods as well as the mass flow rates are analyzed with different mass concentrations of hydrogen peroxide and pressure drops. The bubble slug length is helpful for the selection of the future microreactor length ensuring the complete hydrogen peroxide decomposition. Future studies on the temperature effect on the transient two-phase flow with chemical reactions are recommended.

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.

Experimental Investigation of Performance of a Miniature Direct Methanol Fuel Cell in Short-Term Microgravity

Experimental study of a liquid fed direct methanol fuel cell has been conducted in different gravity environments. A small single cell with 5 cm x 5 cm active area has single serpentine channel on the graphite cathode polar plate and 11 parallel straight channels on the graphite anode flow bed. Cell voltage and current have been measured and two-phase flow in anode channels has been in situ visually observed. The experimental results indicate that the effect of gravity on power performance of the direct methanol fuel cell is large when the concentration polarization governs fuel cells operation. Gravitational effect becomes larger at higher current density. Increasing methanol feeding molarity is conducive to weaken the influence of gravity on performance of liquid fed direct methanol fuel cells. Increasing feeding flow rate of methanol solution from 6 to 15 ml/min could reduce the size of carbon dioxide bubbles, while the influence of gravity still exist. Transport phenomena inside direct methanol fuel cells in microgravity is also analyzed and discussed.

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.

微重力气液两相流动与池沸腾传热

综述了近年来中国科学院微重力重点实验室(国家微重力实验室)完成的一系列微重力气液两相流动与池沸腾传热方面的地基实验、飞行实验和理论研究等方面获得的主要成果.在微重力气液两相流动方面, 提出了半理论Weber数模型用于预测微重力条件下气液两相弹一环状流转换, 并采用Monte Carlo方法, 针对气泡初始尺寸对泡一弹状流转换的影响进行数值研究.通过俄罗斯"和平号"空间站与IL-76失重飞机实验, 获得了微重力下的气液两相流型图, 与此同时在地面利用小尺度毛细管模型模拟了微重力气液两相流动特征.实验测量了微重力气液两相流压降, 并基于微重力流动特性建立了一个泡状流压降关联模型.在微重力池沸腾传热方面, 利用我国返回式卫星完成了两次空间实验, 其中, 第22颗返回式卫星搭载铂丝表面R113池沸腾实验采用控制温度的稳态加热方式, 而实践8号育种卫星搭载平面FC-72池沸腾实验则采用控制加热电压的准稳态加热方式.同时, 还进行了地面常重力和落塔短时微重力条件下的对比实验研究.观察到丝状加热表面微重力时轻微的传热强化现象, 而平板加热表面微重力核态池沸腾低热流时传热强化、高热流时传热恶化.微重力实验中观察到气泡脱落前存在横向运动现象, 据此分析了气泡行为与传热之间关系, 并提出了一个预测丝状加热表面气泡脱落直径的半理论模型.旨在对相关领域的进一步发展和空间两相流系统的应用提供数据及理论支持

微小卫星贮箱内两相流界面变形定位及控制管理

液化气推进技术是微小伴随卫星在轨飞行采用的一种新型推进技术. 微小伴星空间调姿、变轨过程需精确的推力控制, 因此必须了解卫星推进系统中推进剂的形位分布. 本文理论分析了空间微重力环境下液化气推进剂气/液界面的形位分布及变化, 并通过落塔实验验证了微重力环境下有效控制液体推进剂的管理方法.

Liquid-liquid extraction of Pa-233(V) with Aliquat 336

The extraction of protactinium with Aliquat 336 (methyl-tri-caprylyl ammonium chloride) in toluene, cyclohexane and chloroform from HCl, HNO3, H2SO4, HClO4, HF and mixed HCl-HF media was investigated by radioactive tracer technique. Distribution ratios of protactinium between the aqueous solution and the organic phase were determined as a function of shaking time, concentrations of acid in aqueous solution phase, extractant concentration and type of diluents in the organic phase. Aliquat 336 can almost quantitatively extract protactinium from strong HCl solution. At the same time, small amounts of HF in HCl solutions have a strong effect on Pa distribution.