PARALLEL ADAPTIVE SIMULATION OF A PLUNGING LIQUID JET

This paper is concerned with three-dimensional numerical simulation of a plunging liquid jet. The transient processes of forming an air cavity around the jet, capturing an initially large air bubble, and the break-up of this large toroidal-shaped bubble into smaller bubbles were analyzed. A stabilized finite element method （FEM） was employed under parallel numerical simulations based on adaptive, unstructured grid and coupled with a level-set method to track the interface between air and liquid. These simulations show that the inertia of the liquid jet initially depresses the pool＇s surface, forming an annular air cavity which surrounds the liquid jet. A toroidal liquid eddy which is subse- quently formed in the liquid pool results in air cavity collapse, and in turn entrains air into the liquid pool from the unstable annular air gap region around the liquid jet.     

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.     

Numerical simulations of acoustic emissions due to offshore-pile installation

The use of renewable energy is growing due to economic reasons and a spreading sense of ecology. One possible source for renewable energy is wind energy. It is predicted that the number of offshore wind parks will be increasing within the next ten years. The foundation piles of offshore wind power plants are normally installed by impact driving. This leads to an emission of acoustic noise which possibly influences or even damages the sense of hearing of marine mammals like seals and whales. In this contribution the propagation of acoustic waves in air and water due to pile driving is investigated using the finite element method (FEM) with an explicit time integration scheme. In the analyses a mass representing the pile hammer falls on an elastically modeled offshore monopile. The surrounding media is varied between air and water in order to investigate the acoustic emissions due to the pile driving process with respect to the medium. The numerical results are compared to measurements which can be found in literature to verify the numerical models. Furthermore, as a practical example the numerical model can be used in order to simulate noise reduction techniques such as the bubble curtain. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mathematical modelling of the flotation deinking process

The overall flotation deinking process can be divided into four basic microprocesses:

• 1.(1) collision or capture of an (ink) particle by an air bubble
• 2.(2) adhesion of an (ink) particle to the air bubble by sliding
• 3.(3) development of a three-phase contact at the air bubble/water/particle interface, and
• 4.(4) bubble/particle stability or instability after an aggregate is formed each of these microprocesses have an associated probability that they will occur successfully in a flotation cell.

In this paper, the associated probabilities of each microprocess are employed in the development of a kinetic- or population balance-type model of the overall flotation process. The overall model contains two kinetic constants: the first, k₁ governs the overall probability of a free ink particle successfully intercepting and adhering to an air bubble; the second, k₂ is a measure of the probability that a bubble/particle aggregate pair will become unstable and split to yield a “new” free ink particle.The solution to the kinetic model is presented in terms of k₁ and k₂, which are themselves functions of system parameters such as bubble and particle physical properties (e.g., diameter, density), fluid properties (e.g., viscosity, surface tension), etc. From this solution, a definition of a theoretical flotation efficiency, as well as other system performance parameters are presented.     

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 study of formation of a series of bubbles at a submerged orifice

Bubble formation from a submerged orifice is widely applied in bio-process and chemical reaction systems. In this study, the effects of different orifice diameters and contact angles in Period-I and Period-II regimes are studied systematically on a 2D axisymmetric domain. Simulation results are presented from the formation of the first bubble and explained by means of the surrounding fluid field, bubble interaction, and bubble aspect ratio.The orifice diameter is varied from 0.6 mm to 3 mm. The numerical results show that the detachment time of all bubbles remains constant (in time) for smaller orifice diameters (d_a ≤ 1.5 mm), while the detachment time of the first bubble is different from the rest of the bubbles for larger orifice diameters (d_a ≥ 2 mm), which is due to the different surrounding flow field. Contact angles from 60° to 165° are considered for the gas flow rates in the regime of bubble pairing, and it is observed that the bubble detachment time decreases when the contact angle increases, and it converges to a constant value when the contact angle is larger than 135°. In addition, the transition from period doubling to deterministic chaos (in which there is a variable number of bubbles within each period) is observed.A new scenario of inserting a submerged tube upward into the liquid is considered and compared to the previous cases. It is observed that when the tube is vertically inserted into the liquid, the bubble detachment time is even smaller because of higher influence from the surrounding liquid field, leading to a different phenomenon from the non-inserted tube cases.     

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.     

Direct numerical simulation of single gas bubbles in pure and contaminated liquids

Disperse gas bubbles play an important role in many industrial applications. Knowing the rising velocity, the interfacial area, or the critical size for break–up or coalescence in different systems can be crucial for the process design. Hence, knowing the fundamental behaviour of a single bubble appears mandatory for the examination of bubble swarms and for the Euler–Lagrange or Euler–Euler modelling of disperse systems. In the present work a level–set–based volume–tracking method is implemented into the CFD–code OpenFOAM to follow the free interface of a single bubble. The volume–tracking method is coupled with a transport model for surfactants on the interface, including adsorption and desorption processes. The dependency of surface tension on the local surfactant concentration on the interface is modelled by a non–linear (Langmuir) equation of state. Marangoni forces, resulting from surface tension gradients, are included. The rise of a single air bubble (i) in pure water and (ii) in the presence of surfactants of different strengths is simulated. The results show good agreement with available correlations from literature. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of wave processes in a vapor-liquid medium

Numerical methods for the simulation of nonlinear wave processes in a vapor-liquid medium with a model two-phase spherical symmetric cell, with a pressure jump at its external boundary are considered. The viscosity and compressibility of the liquid, as well as the space variation of pressure in the vapor, are neglected. The problem is described by the heat equations in the vapor and liquid, and by a system of ODEs for the velocity, pressure, and radius at the bubble boundary. The equations are discretized in space by an implicit finite-volume scheme on a dynamic adaptive grid with grid refinement near the bubble boundary. The total time derivative is approximated by a method of backward characteristics. “Nonlinear” iterations are implemented at each time step to provide a specified high accuracy. The results of numerical experiments are presented and discussed for the critical thermodynamic parameters of water, for some initial values of the bubble radius and pressure jump.     

A model for flow induced by steady air venting and air sparging

A numerical model for the axisymmetric subsurface flow during steady air venting and/or sparging from underground wells along the same vertical axis is presented. A common observation that for steady air injection or extraction, motion of the groundwater is negligible is first justified theoretically. Thus in the two-phase problem, only air is moving, while water remains stationary. Nevertheless there can be significant changes in water saturation because the capillary pressure varies nonlinearly with the air saturation. For pure air venting, the effects of upconing of the water table are studied. For pure air sparging, the role of air compressibility ignored by earlier authors is reassessed, and issues of design interest such as the radius of influence and discharge rates are examined. Comparisons with published measurements are discussed. Physical implications of concurrent venting and sparging are also discussed.     

Update to the MUSIG model in ANSYS CFX for reliable modelling of bubble coalescence and breakup

The MUSIG (Multiple Size Group) model in the commercial CFD code ANSYS CFX is a population balance approach for describing binary bubble coalescence and breakup events. It is widely used in the simulation of poly-dispersed bubbly flows. The purpose of this work is to identify the internal inconsistencies in the discrete method that is applied for the solution of the population balance equation in MUSIG, and to propose an internally consistent one for discretising the source and sink terms that result from bubble coalescence and breakup. The new formulation is superior to the existing ones in preserving both mass and number density of bubbles, allowing arbitrary discretisation schemes and is free of costly numerical integrations. The numerical results on the evolution of bubble size distributions in bubbly flows reveal that the inconsistency in the original MUSIG regarding bubble breakup is non-negligible for both academic and practical cases. The discussion on the effect of internal inconsistency as well as updates to the model presented in this work are necessary and important for calibration of bubble coalescence and breakup models using the MUSIG approach.     

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.     

Numerical path integration method based on bubble grids for nonlinear dynamical systems

A new numerical path integration method based on bubble grids for nonlinear dynamical systems is presented in this paper. The ordinary differential equations for the first and second order moments are derived on the basis of the Gaussian closure method. Then the probability density values on the bubble nodes in the computational domain can be calculated via the obtained method. The good performance of the resulting method is finally shown in the numerical examples by using some specific nonlinear dynamical systems: Duffing oscillator subjected to harmonic and stochastic excitations, and Duffing–Rayleigh oscillator subjected to harmonic and stochastic excitations.     

Vortex shedding behind a rising bubble and two-bubble coalescence: A numerical approach

In this work, we present the computational results on the wake instability in wobbling bubble regime as well as on the coalescence of two bubbles in different shape regimes. This is a continuation of our previous studies on the dynamics of a single gas bubble rising in a viscous liquid (see [A. Smolianski, H. Haario, P. Luukka, Computational Study of Bubble Dynamics, Research Report 86, Lappeenranta University of Technology, Finland]), and we use the same, finite-element/level-set/operator-splitting method that was proposed in [A. Smolianski, Numerical Modeling of Two-Fluid Interfacial Flows, Ph.D. Thesis, University of Jyväskylä, 2001]. The numerical method allows 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 a good mass conservation. Due to the highly unstable and small-scale nature of the considered problems there are very few experimental investigations, but the comparison with available experimental data confirms a good accuracy of our numerical predictions. Our studies show that plausible results can be obtained with two-dimensional numerical simulations, when a single buoyant bubble or a coalescence of two bubbles is considered.     

Small Oscillations of a Soap Bubble

We analyze the spectrum of small oscillations of a soap bubble surrounded by incompressible inviscid air. The analysis is based on the linearization of the exact system for the dynamics of free fluid films. The underlying model takes into account variations in the film's shape and thickness. The resulting dispersion relationship shows the fundamental interplay among these variations. In the limit of vanishing soap bubble mass, the presented dispersion relationship agrees with Rayleigh's classical formula for an oscillating liquid droplet under the influence of surface tension.     

A composite Level Set and Extended-Domain-Eigenfunction Method for simulating 2D Stokes flow involving a free surface

In this paper, the Extended-Domain-Eigenfunction-Method (EDEM) is combined with the Level Set Method in a composite numerical scheme for simulating a moving boundary problem. The liquid velocity is obtained by formulating the problem in terms of the EDEM methodology and solved using a least square approach. The propagation of the free surface is effected by a narrow band Level Set Method. The two methods both pass information to each other via a bridging process, which allows the position of the interface to be updated. The numerical scheme is applied to a series of problems involving a gas bubble submerged in a viscous liquid moving subject to both an externally generated flow and the influence of surface tension.     

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.     

Numerical modelling of bubbly flows with and without heat and mass transfer

In this paper, modelling gas–liquid bubbly flows is achieved by the introduction of a population balance equation combined with the three-dimensional two-fluid model. For gas–liquid bubbly flows without heat and mass transfer, an average bubble number density transport equation has been incorporated in the commercial code CFX5.7 to better describe the temporal and spatial evolution of the geometrical structure of the gas bubbles. The coalescence and breakage effects of the gas bubbles are modelled according to the coalescence by the random collisions driven by turbulence and wake entrainment while for bubble breakage by the impact of turbulent eddies. Local radial distributions of the void fraction, interfacial area concentration, bubble Sauter mean diameter, and gas and liquid velocities, are compared against experimental data in a vertical pipe flow. Satisfactory agreements for the local distributions are achieved between the predictions and measurements. For gas–liquid bubbly flows with heat and mass transfer, boiling flows at subcooled conditions are considered. Based on the formulation of the MUSIG (multiple-size-group) boiling model and a model considering the forces acting on departing bubbles at the heated surface implemented in the computer code CFX4.4, comparison of model predictions against local measurements is made for the void fraction, bubble Sauter mean diameter, interfacial area concentration, and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Good agreement is achieved with the local radial void fraction, bubble Sauter mean diameter, interfacial area concentration and liquid velocity profiles against measurements. However, significant weakness of the model is evidenced in the prediction of the vapour velocity. Work is in progress through the consideration of additional momentum equations or developing an algebraic slip model to account for the effects of bubble separation.     

ON THE STABILITY OF THE RESIDUAL-FREE BUBBLES FOR THE NAVIER-STOKES EQUATIONS

This paper considers the Calerkin finite element method for the incompressible Navier-Stokes equations in two dimensions, where the finite-dimensional spaces employed consist of piecewise polynomials enriched with residual-free bubble （RFB） functions. The stability features of the residual-free bubble functions for the linearized Navier-Stokes equations are analyzed in this work. It is shown that the enrichment of the velocity space by bubble functions stabilizes the numerical method for any value of the viscosity parameter for triangular elements and for values of the viscosity parameter in the vanishing limit case for quadrilateral elements.