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.

Thermocapillary Migration of Deformable Bubbles at Moderate to Large Marangoni Number in Microgravity

Using the level-set method and the continuum interface model, the axisymmetric thermocapillary migration of gas bubbles in an immiscible bulk liquid with a temperature gradient at moderate to large Marangoni number is simulated numerically. Constant material properties of the two phases are assumed. Steady state of the motion can always be reached. The terminal migration velocity decreases monotonously with the increase of the Marangoni number due to the wrapping of isotherms around the front surface of the bubble. Good agreements with space experimental data and previous theoretical and numerical studies in the literature are evident. Slight deformation of bubble is observed, but no distinct influence on the motion occurs. It is also found that the influence of the convective transport of heat inside bubbles cannot be neglected at finite Marangoni number, while the influence of the convective transport of momentum inside bubbles may be actually negligible.

气泡/液滴Marangoni迁移数值模拟研究

本文利用Level Set方法,数值模拟了微重力情况下气泡/液滴的Marangoni迁移现象,分析了Marangoni数对迁移速度的影响.数值模拟结果表明,随着Marangoni数的增大,非线性热对流效应的影响会逐渐增大,导致沿相界面温度分布趋于均匀,从而降低迁移运动的驱动力,使气泡/液滴的迁移速度随Marangoni数的增加而逐渐减小.

单气泡池沸腾过程中的气泡动力学数值模拟研究

本文以单气泡池沸腾作为核态池沸腾现象的模型体系,以避免汽化核心分布的随机性及气泡间相互作用等引起的复杂性,深入研究核态池沸腾过程中的相变及气泡周围细观流动与传热规律。基于Level Set方法建立了描述单气泡池沸腾问题的数值模型,模拟了气泡的生长与脱落过程,并分析了气泡周围细观流动与传热。

液滴热毛细迁移流场和温度场拓扑结构的演变

热毛细迁移现象是流体颗粒(液滴/气泡)在非均匀温度场中由于界面温度梯度引起的非均匀界面张力驱动的运动。它不仅是流体力学中的经典问题之一,而且在诸如空间材料制备、空间流体和热管理系统等应用中也有着重要的应用。本文利用投影法求解了微重力条件下可变形液滴的轴对称热毛细迁移问题,控制方程组基于Level-set方法和连续界面张力模型。数值计算的终端迁移速度和空间实验结果相一致。计算结果表明,不同Marangoni数(Ma)情形具有相同的流场拓扑结构;但随着Ma数的增加,温度场的拓扑结构变化极大——在较小Ma数情形中,液滴内的最小温度发生在液滴尾部滞止点处。当Ma数超过某个临界值(10～20之间)后,最小温度点跳进液滴内部,并随着Ma数的增加而不断上移;液滴内部冷区最初呈球帽状,但其中部厚度随Ma数增加不断减小,同时向外扩展,形成外缘不断增厚的球壳状冷区;当Ma数超过另一临界值(约100)后,球壳状冷区在液滴轴线处破裂,冷区的拓扑结构转变为环状,最小温度点也随之远离轴线,不断向外、向下移动。温度场拓扑结构的变化反映了对流效应对液滴内部热量传输的影响不断增强,也强烈影响着液滴热毛细迁移终端速度随Ma数的变化特征。

Elliptic Object Detection Based on Shape Preserving and Active Contour

本文通过形状约束方程(组)与一般主动轮廓模型结合,将目标形状与主动轮廓模型融合到统一能量泛函模型中,提出了一种形状保持主动轮廓模型即曲线在演化过程中保持为某一类特定形状。模型通过参数化水平集函数的零水平集控制演化曲线形状,不仅达到了分割即目标的目的,而且能够给出特定目标的定量描述。根据形状保持主动轮廓模型,建立了一个用于椭圆状目标检测的统一能量泛函模型,导出了相应的Euler-Lagrange常微分方程并用水平集方法实现了椭圆状目标检测。此模型可以应用于眼底乳头分割,虹膜检测及相机标定。实验结果表明,此模型不仅能够准确的检测出给定图像中的椭圆状目标,而且有很强的抗噪、抗变形及遮挡性能。

基于形状保持主动轮廓模型长直条的检测

通过形状约束方程(组)与一般主动轮廓模型结合,将目标形状与主动轮廓模型融合到统一能量泛函模型中,提出一种形状保持主动轮廓模型。模型通过参数化水平集函数的零水平集表示某一类特定形状,不仅达到了分割即目标的目的,而且能够给出特定目标的定量描述。根据形状保持主动轮廓模型,建立一个用于长直条状目标检测的统一能量泛函模型,导出相应的Euler-Lagrange常微分方程并用水平集方法实现了长直条状区域的检测。此形状保持模型的一种特殊情况可以用于直线状地平(海天)线提取。实验结果表明,该模型不仅能够准确地检测出给定图像中的长直条状区域而且有很强的抗噪、抗变形及遮挡性能。

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

Tuning the Energy Level of Organic Sensitizers for High-Performance Dye-Sensitized Solar Cells

Cost-effective organic sensitizers will play a pivotal role in the future large-scale production and application of dye-sensitized solar cells. Here we report two new organic D-pi-A dyes featuring electron-rich 3,4-ethylenedioxythiophene- and 2,2'-bis(3,4-ethylenedioxythiophene)-conjugated linkers, showing a remarkable red-shifting of photocurrent action spectra compared with their thiophene and bithiophene counterparts. On the basis of the 3-f{5'-[N,N-bis(9,9-dimethylfluorene-2-yl)phenyl]-2,2'-bis(3,4-ethylenedioxythiophene)-5-yl}2-cyanoacrylic acid dye, we have set a new efficiency record of 7.6% for solvent-free dye-sensitized solar cells based on metal-free organic sensitizers. Importantly, the cell exhibits an excellent stability, keeping over 92% of its initial efficiency after 1000 h accelerated tests under full sunlight soaking at 60 degrees C. This achievement will considerably encourage further design and exploration of metal-free organic dyes for higher performance dye-sensitized solar cells.

Analysis of the sea-level variability along the Chinese coast and estimation of the impact of a CO2-perturbed atmospheric circulation

The relationship between monthly sea-level data measured at stations located along the Chinese coast and concurrent large-scale atmospheric forcing in the period 1960-1990 is examined. It is found that sea-level varies quite coherently along the whole coast, despite the geographical extension of the station set. A canonical correlation analysis between sea-level and sea-level pressure (SLP) indicates that a great part of the sea-level variability can be explained by the action of the wind stress on the ocean surface. The relationship between sea-level and sea-level pressure is analyzed separately for the summer and winter half-years. In winter, one factor affecting sea-level variability at all stations is the SLP contrast between the continent and the Pacific Ocean, hence the intensity of the winter Monsoon circulation. Another factor that affects coherently all stations is the intensity of the zonal circulation at mid-latitudes. In the summer half year, on the other hand, the influence of SLP on sea-level is spatially less coherent: the stations in the Yellow Sea are affected by a more localized circulation anomaly pattern, whereas the rest of the stations is more directly connected to the intensity of the zonal circulation. Based on this analysis, statistical models (different for summer and winter) to hindcast coastal sealevel anomalies from the large-scale SLP field are formulated. These models have been tested by fitting their internal parameters in a test period and reproducing reasonably the sea-level evolution in an independent period. These statistical models are also used to estimate the contribution of the changes of the atmospheric circulation on sea-level along the Chinese coast in an altered climate. For this purpose the ouput of 150 year-long experiment with the coupled ocean-atmosphere model ECHAM1-LSG has been analyzed, in which the atmospheric concentration of greenhouse gases was continuously increased from 1940 until 2090, according to the Scenario A projection of the Intergovermental Panel on Climate Change. In this experiment the meridional (zonal) circulation relevant for sea-level tends to become weaker (stronger) in the winter half year and stronger (weaker) in summer. The estimated contribution of this atmospheric circulation changes to coastal sea-level is of the order of a few centimeters at the end of the integration, being in winter negative in the Yellow Sea and positive in the China Sea with opposite signs in the summer half-year.

COASTAL SEA-LEVEL AND THE LARGE-SCALE CLIMATE STATE - A DOWNSCALING EXERCISE FOR THE JAPANESE ISLANDS

A major problem which is envisaged in the course of man-made climate change is sea-level rise. The global aspect of the thermal expansion of the sea water likely is reasonably well simulated by present day climate models; the variation of sea level, due to variations of the regional atmospheric forcing and of the large-scale oceanic circulation, is not adequately simulated by a global climate model because of insufficient spatial resolution. A method to infer the coastal aspects of sea level change is to use a statistical ''downscaling'' strategy: a linear statistical model is built upon a multi-year data set of local sea level data and of large-scale oceanic and/or atmospheric data such as sea-surface temperature or sea-level air-pressure. We apply this idea to sea level along the Japanese coast. The sea level is related to regional and North Pacific sea-surface temperature and sea-level air pressure. Two relevant processes are identified. One process is the local wind set-up of water due to regional low-frequency wind anomalies; the other is a planetary scale atmosphere-ocean interaction which takes place in the eastern North Pacific.