Numerical simulation of bubble growth in a supersaturated solution

In this paper, a Volume of Fluid (VOF) based approach to simulate the growth of a pre-existing bubble in a supersaturated solution is developed and implemented in OpenFOAM^Ⓡ. The model incorporates the Compressive Continuous Species Transfer approach to describe the transport of dissolved gas and surface tension is treated using the Sharp Surface Force method. The driving force for bubble growth is defined using Fick’s 1^st law and a Sherwood number based correlation. The source terms for the governing equations are implemented by extending the work by Hardt and Wondra, J. Comp. Phys. 227 (2008) 5871–5895. The predictions of the proposed solver is compared against theoretical models for bubble growth in supersaturated solutions. The effect of spurious currents, which are generated while modelling surface tension, on bubble growth is also investigated. The proposed approach is used to model the growth of a rising bubble in the supersaturated solution.     

Effect of a viscous fluid flow past a spherical gas bubble on the growth of its radius

A nearly spherical gas bubble expands adiabatically in a viscous incompressible fluid flowing past it. The Rayleigh-Plesset formula for the growth of the bubble radius is modified due to the flow of the viscous fluid.     

The dynamics of the growth of vapour bubbles in a superheated liquid

The boiling of a superheated liquid is investigated assuming that, in the initial state, the liquid and a vapour bubble (or a system of vapour bubbles) are in mechanical and thermal equilibrium. It is shown that the state of a mixture of a liquid and bubbles is unstable due to the action of capillary forces. Linear and non-linear solutions are constructed describing the emergence of the system from the unstable state and also the unbounded growth of a single bubble and the transition into a stable vapour-liquid state when there are bubbles distributed throughout the volume in the initial state.     

Analytic solutions of unstable interfaces for all density ratios in axisymmetric flows

An interface between two fluids subject to an external force is hydrodynamically unstable. We extend the Layzer-type potential flow model for unstable interfaces to the system of finite density ratio in axially symmetric geometry and derive analytic solutions for growth rates of unstable interfaces over all times. The analytic expressions for bubble growth rates at finite times are given for arbitrary Atwood number. Predictions of the analytic solutions for growth rates are in excellent agreements with numerical results.     

Mesoscopic Numerical Study on Flow Boiling Heat Transfer Performance in Channels With Multiple Rectangular Heaters北大核心CSCD

The flow boiling phenomenon in a channel with multiple rectangular heaters under a constant wall temperature was numerically studied with the lattice Boltzmann method. The effects of spacings between heaters, heater lengths and heater surface wettabilities on the bubble morphology, the bubble area and the heat flux on the heater surface, were studied. The results show that, the bubble growth rate increases with the spacing between heaters. The larger the bubble area is, the earlier the nucleated bubbles will leave the heater surface. The corresponding boiling heat transfer performance increases by 12% with the spacing between heaters growing from 250 lattices to 1 000 lattices. On the other hand, the longer the heater length is, the earlier the bubble will nucleate and leave the heater surface, and the better the boiling heat transfer performance will be. The boiling heat transfer performance increases by 13% with the heater length rising from 16 lattices to 22 lattices. In addition, the bubble nucleates later on the hydrophilic surface than on the hydrophobic surface. Compared with the hydrophilic surface, the hydrophobic surface retains residual bubbles after the leaving of bubbles from the heater. The average heat flux and the bubble area of the hydrophilic surface are less than those of the hydrophobic surface. With the contact angle changing from 77° to 120°, the heat transfer performance increases by 26%. Finally, the orthogonal test results indicate that, the wettability of the heat exchanger surface has the greatest influence on the flow boiling heat transfer performance, while the heater length has the least influence. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

气核尺度对空化现象的影响

本文计算了不同尺度的气核对空化现象的影响.结果表明,不同尺度的气核具有不同的空化起始条件及不同的发育和溃灭过程.在计算、模拟空化现象时应当考虑到这一因素.     

水中孤立气泡对声波的耗散作用

本文从气泡壁面线性运动方程出发,对水中孤立气泡的声耗散机理进行了分析。结果表明:气泡对声波具有明显的耗散作用;耗散作用由声散射和声吸收两部份组成;热传导在气泡对声波的耗散作用中起决定性作用;气泡处于谐振时,对声波具有最大的耗散作用。     

Analysis of helium bubble growth in radioactive waste

A discrete kinetic model for the growth of helium bubbles in plutonium is proposed and analyzed. This model captures some relevant qualitative features of the time behavior of the distribution of bubble sizes. Analytic formulae for the solutions are given, which agree reasonably well with the numerical solutions, and a rigorous existence theory is established for three different equivalent formulations.     

Analysis of the cell boundary element methods for convection dominated convection-diffusion equations

We introduce two kinds of the cell boundary element (CBE) methods for convection dominated convection-diffusion equations: one is the CBE method with the exact bubble function and the other with inexact bubble functions. The main focus of this paper is on inexact bubble CBE methods. For inexact bubble CBE methods we introduce a family of numerical methods depending on two parameters, one for control of interior layers and the other for outflow boundary layers. Stability and convergence analysis are provided and numerical tests for inexact bubble CBEs with various choices of parameters are presented.     

Numerical simulation of multiply connected axisymmetric discontinuous incompressible potential flows

The ascend and evolution of an axisymmetric gas bubble are studied numerically using an inviscid incompressible potential flow model. The volume of the gas bubble varies adiabatically. The transition from a simply connected bubble to a doubly connected toroidal one and its interaction with the free surface are simulated. The change in connectedness is accompanied by a nonzero velocity circulation and a discontinuous velocity potential occurring over an arbitrary toroidal liquid surface enclosing the bubble.     

Numerical modelling of EHD effects on heat transfer and bubble shapes of nucleate boiling

In this article, mathematical and numerical models are developed to study pure electrohydrodynamic (EHD) effects on heat transfer and bubble shapes when an initial bubble attached to a superheated horizontal wall in nucleate boiling. In the modelling of EHD effects on heat transfer, an undeformed bubble is considered; the electric body force and Joule heat are added to the momentum and energy equations; governing equations for heat, fluid flow and electric fields are coupled numerically and solved using a non-orthogonal body-fitted mesh system with necessary interfacial treatments at the gas–liquid boundary. While, to study the pure effect of EHD on the deformation of the bubble, the evaluation of a deformable bubble without heat transfer is simulated by volume of fluid (VOF) method based on an axial symmetric Cartesian coordinate system. The simulations indicate that EHD can effectively enhance heat transfer rate of nucleate boiling by influencing the motion of the ring vortex around the bubble and that bubble can be elongated due to the pull in axial direction and push in the negative radial direction by the electric field force.     

Numerical study of dynamics of single bubbles and bubble swarms

A systematic computational study of the dynamics of gas bubbles rising in a viscous liquid is presented. Two-dimensional simulations are carried out. Both the dynamics of single bubbles and small groups of bubbles (bubble swarms) are considered. This is a continuation of our previous studies on the two-bubble coalescence and vortex shedding [A. Smolianski, H. Haario, P. Luukka, Vortex shedding behind a rising bubble and two-bubble coalescence: a numerical approach, Appl. Math. Model. 29 (2005) 615–632]. The proposed numerical method allows us to simulate a wide range of flow regimes, accurately capturing the shape of the deforming interface of the bubble and the surface tension effect, while maintaining the mass conservation. The computed time-evolution of bubble’s position and rise velocity shows a good agreement with the available experimental data. At the same time, the results on the dynamics of bubble interface area, which are, up to our knowledge, presented for the first time, show how much the overall mass transfer would be affected by the interface deformation in the case of the bubble dissolution. Another set of experiments that are of interest for chemical engineers modelling bubbly flows concerns the bubble swarms and their behavior in different bubble-shape regimes. The ellipsoidal and spherical shape regimes are considered to represent, respectively, the coalescing and non-coalescing bubble swarms. The average rise velocities of the bubble swarms are computed and analyzed for both regimes.     

BEM-based numerical study of three-dimensional compressible bubble dynamics in stokes flow

The dynamics of compressible gas bubbles in a viscous shear flow and an acoustic field at low Reynolds numbers is studied. The numerical approach is based on the boundary element method (BEM), which is effective as applied to the three-dimensional simulation of bubble deformation. However, the application of the conventional BEM to compressible bubble dynamics faces difficulties caused by the degeneration of the resulting algebraic system. Additional relations based on the Lorentz reciprocity principle are used to cope with this problem. Test computations of the dynamics of a single bubble and bubble clusters in acoustic fields and shear flows are presented.     

基于马氏域变的房地产泡沫存在概率研究

借鉴股票市场泡沫的研究方法,采用非齐次马氏域变思想对我国2001年1月到2010年12月的房地产市场泡沫进行了分析.通过Johansen协整检验方法从房价数据中分离出泡沫成份,使用极大似然估计方法估计出房价泡沫在两状态下的方差和转移概率的相关参数,进而求出各时期房价泡沫存在的概率.实证结果表明:第一,样本期间大部分时期内我国房地产价格存在泡沫的概率较小,但2009年4月至12月房价存在泡沫的概率较大;第二,提高贷款利率能够减小从无泡沫到有泡沫的概率,同时增加从有泡沫到无泡沫的概率.     

Numerical simulation of the interaction between vortex ring and bubble plume

In present paper a three-dimensional Vortex-In-Cell method with two-way coupling effect was developed to study the bubble plume entrainment by a vortex ring. In this method the continuous flow was calculated by the three-dimensional Vortex-In-Cell method and the bubbles are tracked through bubble motion equation. Two-way coupling effect between continuous flow and dispersed bubbles is considered by introducing a vorticity source term, which is induced by the change of void fraction gradient in each computational cell. After validated by the comparison between experimental measurements and simulation results for the motion of vortex rings and the rising velocity of bubble plume, present method is implemented to simulate the interaction between an evolving vortex ring and a rising bubble plume. It was found that there is little effect of the bubble entrainment to the total circulation of vortex ring while the effect of bubble entrainment to the vortex ring structure is quite obvious. The bubble entrainment by the vortex ring not only changed the vorticity distribution in the vortex structure, but also displaced the positions of the vortex cores. The vorticity in the lower vortex core of the vortex ring decreases more than that in the upper vortex core of the vortex ring while the vortex core in the upper part of the vortex ring is displaced to the center of vortex ring by the entrained bubbles. Smaller bubbles are easier to be entrained by the large scale vortex structure and the transportation distance is in inverse proportion to bubble diameter.     

Modelling bubble rise and interaction with a glass surface

A theoretical model has been developed to analyse bubble rise in water and subsequent impact and bounce against a horizontal glass plate. The multiscale nature of the problem, where the bubble size is on the millimetre range and the film drainage process happens on the micrometre to nanometre scale requires the combined use of different modelling techniques. On the macro scale we solve the full Navier–Stokes equations in cylindrical coordinates to model bubble rise whereas modelling film drainage on the micro scale is based on lubrication theory because the film Reynolds number becomes much smaller than unity. Quantitative predictions of this model are compared with experimental data obtained using synchronised high-speed cameras. Video recording of bubble rise and bounce trajectories are combined with interferometry data to deduce the position and time-dependent thickness of the thin water film trapped between the deformed bubble and the glass plate. Bubble rise velocity indicated that the boundary condition at the bubble surface was tangentially immobile. Quantitative comparisons are presented for bubbles of different size to quantify similarities and differences.     

The effect of surfactant on bubble and thread dynamics

Two examples from the dynamics of surfactant-laden interfacial flow are considered. In the first example, a bubble is rapidly stretched by an imposed flow to form a dumb-bell shape, then the imposed flow is relaxed, so that the bubble evolves solely under the action of surfactant-modified surface tension. In the second example, a surfactant-coated bubble is continually stretched by a steady extensional flow. At sufficiently small strain rates, steady bubble shapes are found, whereas at larger strain rates a long-wave model of the dynamics predicts behavior that is reminiscent of experimentally observed tip-streaming. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Prediction of initial motion of a gas bubble in liquids

Numerical predictions are presented for the motion and distortion of a single gas bubble rising through the liquid. The computations were made with an implicit finite-difference procedure which solves the transient equations of motion throughout the bubble and the liquid, such that the free surface between the gas bubble and the liquid is not a boundary of the computational domain.The predictions compare well with the experimental results of others. Computations are presented for bubble sizes from 0.02 to 0.05 m radius and for bubbles of different gas densities rising in liquids of different densities. Surface tension effects are neglected.     

狭窄流域内气泡破裂现象数值模拟

基于边界积分法建立狭窄流域内气泡破裂数值模型,开发相应的计算程序,分别模拟对称破裂与非对称破裂两类典型工况并与相应实验结果进行对比,计算值与实验值吻合很好,表明三维数值模型的有效性．从狭窄流域内气泡运动的基本现象入手,基于开发的程序系统地研究气泡的对称破裂与非对称破裂,在已有数值研究成果和实验数据基础上,提出气泡破裂的可行性准则,研究分裂后子气泡的动力学特性,并分析距离参数对气泡破裂特性及子气泡动力学行为的影响,总结相应规律,旨在为相关气泡破裂特性研究提供参考．     

Weak oscillations of a gas bubble in a spherical volume of compressible liquid

The following spherically symmetric problem is considered: a single gas bubble at the centre of a spherical flask filled with a compressible liquid is oscillating in response to forced radial excitation of the flask walls. In the long-wave approximation at low Mach numbers, one obtains a system of differential-difference equations generalizing the Rayleigh-Lamb-Plesseth equation. This system takes into account the compressibility of the liquid and is suitable for describing both free and forced oscillations of the bubble. It includes an ordinary differential equation analogous to the Herring-Flinn-Gilmore equation describing the evolution of the bubble radius, and a delay equation relating the pressure at the flask walls to the variation of the bubble radius. The solutions of this system of differential-difference equations are analysed in the linear approximation and numerical analysis is used to study various modes of weak but non-linear oscillations of the bubble, for different laws governing the variation of the pressure or velocity of the liquid at the flask wall. These solutions are compared with numerical solutions of the complete system of partial differential equations for the radial motion of the compressible liquid around the bubble.