海洋大气间气泡的气体输运1)

研究球形小气泡在理想流体的波浪场中的气体扩散过程,把小雷诺数下均匀来流绕流球形气泡的气体交换结果与气泡运动方程耦合在一起进行求解.讨论了溶解于水中的气体浓度、波浪、气泡半径、气泡初始深度对单个气泡气体扩散量的影响.由于气泡云对气体的输运,溶解于水中的气体可出现过饱和状态.对10 m/s风速下气泡云的气体输运量进行了计算,得到水中O     

有气体溶解和没有气体溶解的水中石墨烯间的疏水引力

凝聚物理学界发现溶解在水中的气体会在疏水表面吸附,由此认为当疏水物体间距足够小的时候,吸附在疏水表面的气体会相互联通形成纳米气泡桥,纳米气泡桥连结疏水物质形成疏水引力,但是关于纳米气泡桥的形成过程和形态,力学界还没有一个清晰的描述.采用分子动力学方法,研究了在有气体溶解的水中和没有气体溶解的水中两片石墨烯间的引力作用,分析了两种情况下各相密度分布的变化过程、结构相图的变化过程和平均力势的变化过程,详细阐明了纳米气泡桥的形成和消失过程,并定量计算了纳米气泡桥的作用效果和作用距离.模拟结果表明:两片石墨烯在有气体溶解的水中和无气体溶解的水中的疏水引力都是由纳米气泡桥引起的.当石墨烯间距小于0.5nm时,无论水中是否有气体溶解,疏水引力由真空纳米气泡桥引起;当石墨烯间距大于0.5nm时,在没有气体溶解的水中,疏水引力由水蒸气纳米气泡桥引起;而在有气体溶解的水中,疏水引力由所溶气体形成的纳米气泡桥引起.      

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡, 了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提. 基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽, 模型采用S-A IDDES湍流模型进行湍流封闭, 并采用修正的VOF 方法捕捉自由液面, 数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程, 重点分析了结构物附近的水动力和掺气水流运动特性. 结果表明, 建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程, 较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程; 波浪与直立结构物相互作用产生强湍动掺气水流, 其运动过程十分复杂; 掺气流体输运过程中水气界面周围一直伴随着涡的存在, 其中, 气泡的分裂与周围正负涡量剪切作用密切相关, 且其输运轨迹主要受周围流场的影响; 研究揭示了结构物附近湍动能与掺气特性的关系, 发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.      

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡,了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提.基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽,模型采用S-A IDDES湍流模型进行湍流封闭,并采用修正的VOF方法捕捉自由液面,数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程,重点分析了结构物附近的水动力和掺气水流运动特性.结果表明,建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程,较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程;波浪与直立结构物相互作用产生强湍动掺气水流,其运动过程十分复杂;掺气流体输运过程中水气界面周围一直伴随着涡的存在,其中,气泡的分裂与周围正负涡量剪切作用密切相关,且其输运轨迹主要受周围流场的影响;研究揭示了结构物附近湍动能与掺气特性的关系,发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.     

柱形汇聚激波冲击球形重气体界面的实验研究

采用高速纹影法实验研究了柱形汇聚激波与球形重气体界面相互作用的 Richtmyer-Meshkov不稳定性问题. 激波管实验段基于激波动力学理论设计, 将马赫数为1.2 的平面激波转化为柱形汇聚激波, 气体界面由肥皂膜分隔六氟化硫(内)和空气(外)得到. 采用高速摄影机在单次实验中拍摄激波运动的全过程, 对柱形激波的形成进行了实验验证, 并进一步观测了汇聚激波与球形气体界面相互作用过程中的波系发展和气体界面变形以及反射激波同已变形界面二次作用的流场演化. 结果表明: 当柱形汇聚激波穿过气泡界面以后, 气泡左侧界面极点沿激波传播方向保持匀速运动, 气泡右侧界面发展成为射流结构, 气泡主体发展成为涡环结构; 在反射激波的二次作用下, 流场中无序运动显著增强并很快进入湍流混合阶段.     

柱形装药自由场水中爆炸气泡的射流特性

为研究柱形装药水中爆炸气泡射流特性,采用高速摄影方法在实验水箱中对PETN药柱进行水中爆炸实验,初步获得了有关柱形装药自由场水中爆炸气泡运动的数据资料.研究结果表明:柱形装药自由场水中爆炸形成气泡的运动与球形装药的存在差别,药柱放置方式的变化也会导致水中爆炸气泡产生性质迥异的运动.当药柱水平放置时,气泡表面的运动呈非对...     

波浪场中物质输移离散系数的理论分析

海岸水域波浪引起的物质输移扩散存在着不同于水流引起的物质输移扩散的特征. 通过解析的方法对垂向扩散方程进行求解, 推导出波浪和潮流共同作用下的海岸水域物质输移离散系数的理论公式. 分析中将总的离散系数分为潮流、质量输移流、纯波浪波动作用以及潮流与质量输移流、潮流与波浪波动、质量输移流与波浪波动的相互作用6部分, 通过分析各组成部分特征, 得出海岸水域物质输移中波浪所产生的离散效应主要是质量输移流贡献, 波浪水质点运动及波流相互作用对时间平均离散系数的贡献较小的结论.仅波浪的作用(不包括潮流)的有关结果与他人数值模拟和实验结果的吻合良好, 证明了该理论推导的正确性. 给出了时间平均离散系数及离散系数波动幅值随波浪周期和波高的变化特征,同时将结果应用于垂向环流, 给出了垂向环流对应于质量输移流所产生的离散系数的特征.      

小当量水中爆炸气泡的脉动现象

为研究水中爆炸气泡的脉动现象,设计了用于小当量(10 g TNT)水中爆炸研究的爆炸水箱,在爆炸水箱装置中进行了多次0.125、1.000、3.370、8.000 g TNT当量的PETN球形炸药水中爆炸试验。采用高速摄影系统获得了气泡脉动图片;给出试验范围内气泡脉动过程中气泡直径、速度及加速度随时间变化的拟合曲线。对气泡脉动直径和周期数据进行了分析,这些数据基本符合爆炸相似律。     

静水中气泡上升运动及阻力系数研究

采用模型计算法与实验法结合的方式对静水中气泡上升运动行为进行研究。通过牛顿运动定律,基于不同物理模型,建立气泡在水中运动的微分方程;假设气泡在运动过程中的关键参数取值,推导小气泡在水中浮升过程中的气泡行为预测公式;针对不同流态下的气泡上升关键参数进行适应性分析和算例计算。通过设计气泡上升运动实验,对气泡上升运动公式进行适应性分析,修正关键参数的取值。据此提出一种小气泡上升运动规律的计算方法以及关键参数取值方式及参考区间。     

锌空电池气体扩散电极浸液量对放电性能的影响

锌空电池气体扩散电极在存放和放电过程中,电解液会在毛细力的作用下不断浸入电极。电解液在气体扩散电极中的浸入量与分布情况决定了气体扩散电极中的化学反应活化区,从而影响放电性能。通过实验测量了气体扩散电极开始浸液的4天内浸液量与放电性能的关系,并借助拓扑网络数值模拟电解液浸入多孔介质的过程帮助理解该实验现象。结果表明,随着放电过程的进行,浸液量和分布情况不断变化;气体扩散电极放电性能变化主要分为3个阶段:浸液饱和度为39.4%时放电性能最佳;浸液开始2~24小时进行迅速,浸液饱和度达到81%,放电性能小幅下降;24小时之后浸液增速大幅减缓,放电功率随浸液量增加大幅下降。     

Wave dispersion of a gas in a liquid

A novel method of dispersing a gas in a liquid by pressure pulses generated by waves propagating from a hydrodynamic oscillation generator is proposed. Devices that realize this method, wave dispersers, are created and investigated experimentally. The amplitude-frequency characteristics of the wave processes in the dispersers and the size distribution density functions of the gas bubbles are obtained. It is shown that there are optimal pressure values at the disperser inlet at which a minimum bubble size is achieved. The average diameter of the gas bubbles produced by water wave dispersers in the optimal operation regimes varied on the range from 0.3 to 0.6 mm, depending on the gas flow-rate.     

Bubbles in liquids with phase transition

In the forthcoming second part of this paper a system of balance laws for a multi-phase mixture with many dispersed bubbles in liquid is derived where phase transition is taken into account. The exchange terms for mass, momentum and energy explicitly depend on evolution laws for total mass, radius and temperature of single bubbles. Therefore in the current paper we consider a single bubble of vapor and inert gas surrounded by the corresponding liquid phase. The creation of bubbles, e.g. by nucleation is not taken into account. We study the behavior of this bubble due to condensation and evaporation at the interface. The aim is to find evolution laws for total mass, radius and temperature of the bubble, which should be as simple as possible but consider all relevant physical effects. Special attention is given to the effects of surface tension and heat production on the bubble dynamics as well as the propagation of acoustic elastic waves by including slight compressibility of the liquid phase. Separately we study the influence of the three phenomena heat conduction, elastic waves and phase transition on the evolution of the bubble. We find ordinary differential equations that describe the bubble dynamics. It turns out that the elastic waves in the liquid are of greatest importance to the dynamics of the bubble radius. The phase transition has a strong influence on the evolution of the temperature, in particular at the interface. Furthermore the phase transition leads to a drastic change of the water content in the bubble. It is shown that a rebounding bubble is only possible, if it contains in addition an inert gas. In Part 2 of the current paper the equations derived are sought in order to close the system of equations for multi-phase mixture balance laws for dispersed bubbles in liquids involving phase change.     

On the modeling and simulation of a laser‐induced cavitation bubble

Single cavitation bubbles exhibit severe modeling and simulation difficulties. This is due to the small scales of time and space as well as due to the involvement of different phenomena in the dynamics of the bubble. For example, the compressibility, phase transition, and the existence of a noncondensable gas inside the bubble have strong effects on the dynamics of the bubble. Moreover, the collapse of the bubble involves the occurrence of critical conditions for the pressure and temperature. This adds extra difficulties to the choice of equations of state. Even though several models and simulations have been used to study the dynamics of the cavitation bubbles, many details are still not clearly accounted for. Here, we present a numerical investigation for the collapse and rebound of a laser‐induced cavitation bubble in liquid water. The compressibility of the liquid and vapor are involved. In addition, great focus is devoted to study the effects of phase transition and the existence of a noncondensable gas on the dynamics of the collapsing bubble. If the bubble contains vapor only, we use the six‐equation model for two‐phase flows that was modified in our previous work [A. Zein, M. Hantke, and G. Warnecke, J. Comput. Phys., 229(8):2964‐2998, 2010]. This model is an extension to the six‐equation model with a single velocity of Kapila et al. (Phys. Fluid, 13:3002‐3024, 2001) taking into account the heat and mass transfer. To study the effect of a noncondensable gas inside the bubble, we add a third phase to the original model. In this case, the phase transition is considered only at interfaces that separate the liquid and its vapor. The stiffened gas equations of state are used as closure relations. We use our own method to determine the parameters to obtain reasonable equations of state for a wide range of temperatures and make them suitable for the phase transition effects. We compare our results with experimental ones. Also our results confirm some expected physical phenomena. Copyright © 2013 John Wiley & Sons, Ltd.     

Degassing of water-solved gases in thin pipings during fast pressure drops

 A one-dimensional model is presented, which describes the transient two-phase flow in thin pipes during fast pressure drops and degassing by use of Eulerian and Lagrangian systems. The reduction in dimension is obtained by introduction of a geometry model for bubbly and slug flow regimes. The complete model includes the transient two-phase flow, bubble formation and bubble growth. The flow model predicts rising velocities of bubbles and plugs in arbitrary inclined highly accurate pipes. The mass transfer (diffusion) of the dissolved phase is calculated by the bubble growth model. The quality of the model was examined by simulation of experimental series, whereby water was depressurised from the saturation pressure of the dissolved gas mixture (air), by variation of saturation pressure, pressure gradient and pipe geometry. The results of numerical simulation fit the experimental data well. Received on 17 January 2000     

Acoustics of a Bubbly Liquid with a Weak Chemical Reaction in the Gaseous Phase

The influence of the composition and thermophysical properties of gas-liquid bubbly systems with a dissociating component in the gaseous phase on the laws of small-disturbance propagation and attenuation is investigated. It is found that the reacting gas component in the bubbles significantly affects the sonic-wave attenuation coefficient in the bubbly liquid. This follows from the fact that when a gas bubble is compressed isothermally, a recombination reaction occurs which prevents pressure growth in the bubble.Small-disturbance propagation in bubbly liquids was investigated in a number of publications discussed in review [1]. The acoustics of a bubbly liquid with a gas phase containing active admixtures are of both methodical and practical interest. The dynamics of such multicomponent bubbles were investigated in [2].     

Dynamics of a “collective” bubble in a magma melt flow behind the decompression wave front

Specific features of the dynamics of the wave field structure and growth of a “collective” bubble behind the decompression wave front in the “Lagrangian” section of the formed cavitation zone are numerically analyzed. Two cases are considered: with no diffusion of the dissolved gas from the melt to cavitation nuclei and with the diffusion flux providing an increase in the gas mass in the bubbles. In the first case, it is shown that an almost smooth decompression wave front approximately 100 m wide is formed, with minor perturbations that appear when the front of saturation of the cavitation zone with nuclei is passed. In the case of the diffusion process, the melt state behind the saturation front is principally different: jumps in mass velocity and viscosity are observed in the vicinity of the free surface, and the pressure in the “collective” cavitation bubble remains unchanged for a sufficiently long time interval, despite the bubble growth and intense diffusion of the gas from the melt. It is assumed that the diffusion process (and, therefore, viscosity) actually become factors determining the dynamics of growth of cavitation bubbles beginning from this time interval. A pressure jump is demonstrated to form near the free surface.     

Dissipation of the energy of sound waves in small gas bubbles

One of the main factors affecting the dynamics of homogeneous solution type pulse reactors is the formation of gas bubbles on the fission-fragment tracks [1, 2]. The behavior of the reactor depends very considerably on the size (10^?5 cm) and growth rate of these bubbles [2], and it is, accordingly, a very important matter to study these properties. One convenient means of doing this lies in the acoustic method. The behavior of gas bubbles in the field of a sound wave has been studied in a large number of papers and reviews [3, 4]. In this paper we shall see the approximation of a sound wave of small amplitude to consider the dissipation of sound-wave energy in a gas bubble, at the same time allowing for inertia, surface tension, viscosity, heat transfer, and the diffusion of gas through the surface of the bubble.     

Disengagement of the gas phase in bubble columns

The technique of gas disengagement is popularly used to assess the bubble size distribution in bubble columns. The technique involves the dynamic measurement of dispersion height when the gas supply is stopped. In this paper a mathematical model has been proposed for the process of dynamic gas disengagement. It has been shown that the initial faster disengagement is due to the presence of internal liquid circulation and not due to the presence of very large bubbles. Further, slower disengagement has been attributed to the transition from heterogeneous dispersion to homogeneous dispersion. The new model also explains the effects of superficial gas velocity, column diameter, column height and liquid phase physical properties on the gas disengagement.     

Numerical simulation of a single bubble by compressible two‐phase fluids

The present work deals with the numerical investigation of a collapsing bubble in a liquid–gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non‐conservative transport equations of the pressure law coefficients. The original first‐order discretization is extended to higher order applying second‐order ENO reconstruction to the primitive variables. The derivation of the non‐conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale‐based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2‐D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser‐induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd.