• Nuclear Engineering and Technology
• /
• v.44 no.6
• /
• pp.623-638
• /
• 2012

The volume of fluid (VOF) model of FLUENT and the lattice Boltzmann method (LBM) are used to simulate two-phase flows. Both methods are validated for static and dynamic bubble test cases and then compared to experimental results. The VOF method does not reduce the spurious currents of the static droplet test and does not satisfy the Laplace law for small droplets at the acceptable level, as compared with the LBM. For single bubble flows, simulations are executed for various Eotvos numbers, Morton numbers and Reynolds numbers, and the results of both methods agree well with the experiments in the case of low Eotvos numbers. For high Eotvos numbers, the VOF results deviated from the experiments. For multiple bubbles, the bubble flow characteristics are related by the wake of the leading bubble. The coaxial and oblique coalescence of the bubbles are simulated successfully and the subsequent results are presented. In conclusion, the LBM performs better than the VOF method.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.3
• /
• pp.10-17
• /
• 2011

Free energy based lattice Boltzmann method (LBM) has been used to simulate the contact angle and the bubble necking with large density ratio. LBM with the proper contact angle model is able to reduce the spurious currents and eliminate the singularity in the contact lines. The numerical results of the contact angles are satisfied with the Youngs law. For bubble necking flows, simulations are executed for various viscosities and contact angles. The phenomena of the bubble necking are simulated successfully and the subsequent results are presented. The present method is also applicable to the nucleate boiling flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.1
• /
• pp.85-94
• /
• 2002

In this study an adaptive node generation procedure in the Element-free Galerkin (EFG) method using bubble-meshing technique is Proposed. Since we construct the initial configuration of nodes by subdivision of background cell, abrupt changes of inter-nodal distance between higher and lower error regions are unavoidable. This unpreferable nodal spacing induces additional errors. To obtain the smooth nodal configuration the nodal configurations are regenerated by bubble-meshing technique. This bubble meshing technique was originally developed to generate a set of well-shaped triangles and tetrahedra. In odder to evaluate the effect of abrupt changes of nodal spacing, one-dimensional problems with various nodal configurations mere investigated. To demonstrate the performance of proposed scheme, the sequences of making optimal nodal configuration with bubble meshing technique are investigated for several problems.

• Korean Chemical Engineering Research
• /
• v.49 no.1
• /
• pp.83-88
• /
• 2011

Holdup characteristics of small bubbles were investigated in a viscous slurry bubble column. The phase holdup of small bubbles was obtained from the knowledge of total bubble(gas) holdup and large bubble holdup, which were measured by mean of static pressure drop method and dual resistivity probe method, respectively. Effects of gas velocity, viscosity of continuous liquid phase and solid fraction in the slurry phase on the small bubble holdup as well as holdups of total bubble(gas) and large bubble in a viscous slurry bubble column. The small bubble holdup increased with increasing gas velocity but decreased with increasing liquid viscosity or solid fraction in the slurry phase. In addition the fraction of small bubble in the total bubble(gas) holdup increased with increasing gas velocity but decreased with increasing liquid viscosity or solid fraction in the slurry phase. It was revealed that the rising velocity of large bubble did not related to the holdup of small bubble in a viscous slurry bubble column.

• Journal of Mechanical Science and Technology
• /
• v.15 no.7
• /
• pp.931-940
• /
• 2001

A numerical method for simulating bubble motion during nucleate boiling is presented. The vapor-liquid interface is captured by a level set method which can easily handle breaking and merging of the interface and can calculate an interfacial curvature more accurately than the VOF method using a step function. The level set method is modified to include the effects of phase change at the interface and contact angle at the wall as well as to achieve mass conservation during the whole calculation procedure. Also, a simplified model to predict the heat flux in a thin liquid microlayer is developed. The method is applied for simulation of a sliding bubble on a vertical surface to further understand the physics of partial boiling. Based on the computed results, the effects of contact angle, wall superheat and phase change on a sliding bubble are quantified.

• Journal of the Korean Society of Visualization
• /
• v.15 no.2
• /
• pp.48-58
• /
• 2017

To analyze bubbles generated by an ABB (Air Bubble Barrier), we developed image processing procedure and statistical analysis method. Air was discharged from 5 mm nozzle as swarm form at the bottom of 1 m3 water tank. Flow rates of discharged air are ranged from 2 L/min to 20 L/min and these are corresponding to Reynolds number of 1766-17663. Rise velocity of bubble is extracted by using image process pretending intrusive method. Mean equivalent velocity was calculated using void fraction weighting factor. Bubble diameter is obtained and compared with correlations in the literature. Also, we present a correlation according to the result of this study. Mean velocity and mean diameter of bubbles increase with increasing gas Reynolds number. But these parameters show an asymptotic trend when they approach to high Reynolds number.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.48-55
• /
• 2000

A basic research of air bubble locus around the ship hull is performed to know the bubble's behavior which is generated around the hull surface. In this paper. bubble's behavior around the ship hull is calculated and experimented with the variation of bubble size and the location of bubble generation. For the use of the equation for bubble locus, Kawakita's equation which include the effect of bubble buoyancy is adopted. For the numerical simulation. a finite difference method based on the MAC method is used. Also experiments are performed about series 60 hull form at CWC. The results of calculation are compared with those of experiment.

• Nuclear Engineering and Technology
• /
• v.53 no.4
• /
• pp.1119-1126
• /
• 2021

Spacer grid with mixing vane had been widely used in nuclear reactor core. One of the main feather of spacer grid with mixing vane was that strong swirl flow was formed after the spacer grid. The swirl flow not only changed the bubble generation in the near wall field, but also affected the bubble behaviors in the center region of the subchannel. The interaction between bubble and the swirl flow was one of the basic phenomena for the two phase flow modeling in fuel assembly. To obatin better understanding on the bubble behaviors in swirl flow, full three dimension numerical simulations were conducted in the present paper. The swirl flow was assumed in the cylindral calculation domain. The bubble interface was captured by Volume Of Fluid (VOF) method. The properties of saturated water and steam at different pressure were applied in the simulation. The bubble trajectory, motion, shape and force were obtained based on the bubble parameters captured by VOF. The simulation cases in the present study included single bubble with different size, at different angular velocity conditions and at different pressure conditions. The results indicated that bubble migrated to the center in swirl flow with spiral motion type. The lateral migration was mainly related to shear stress magnitude and bubble size. The bubble moved toward the center with high velocity when the swirl magnitude was high. The largest bubble had the highest lateral migration velocity in the present study range. The effect of pressure was small when bubble size was the same. The prelimenery simulation result would be beneficial for better understanding complex two phase flow phenomena in fuel assembly with spacer grid.