• Journal of the Korean Society of Visualization
• /
• v.21 no.1
• /
• pp.94-102
• /
• 2023

Three design parameters were considered in this study: outlet nozzle angle (30°, 60°, 80°), neck length (1 mm, 3 mm), and flow rate (0.5, 0.6, 0.7, 0.8 lpm). A neck diameter of 0.5 mm induced cavitation flow at a venture nozzle. A secondary transparent chamber was connected after ejection to increase bubble duration and shape visibility. The bubble size was estimated using a Gaussian kernel function to identify bubbles in the acquired images. Data on bubble size were used to obtain Sauter's mean diameter and probability density function to obtain specific bubble state conditions. The degree of bubble generation according to the bubble size was compared for each design variable. The bubble diameter increased as the flow rate increased. The frequency of bubble generation was highest around 20 ㎛. With the same neck length, the smaller the CV number, the larger the average bubble diameter. It is possible to increase the generation frequency of smaller bubbles by the cavitation method by changing the magnification angle and length of the neck. However, if the flow rate is too large, the average bubble diameter tends to increase, so an appropriate flow rate should be selected.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.185-187
• /
• 2009

A density based method with homogeneous cavitation model to investigate cavitation-bubble collapsing behavior is proposed and applied to bubble-shock interaction problems. By applying this method, cylindrical bubbles located in the liquid and incident liquid shock wave are computed. Bubble collapsing behavior, shock-bubble interaction and shock transmission/reflection pattern are investigated.

• Journal of The Korean Astronomical Society
• /
• v.29 no.spc1
• /
• pp.253-254
• /
• 1996

To understand the dynamical structures of stellar wind bubble, one and two-dimensional calculations has been performed. Using FCT Code with cooling effects and assuming constant mass loss rate and ambient medium density, we could divide stellar winds into the regime of slow and fast winds. The slow wind driven bubble shows initially radiative and becomes partially radiative bubble in which shocked stellar wind zone is still adiabatic. In contrast., the fast wind driven bubble shows initially fully adiabatic and becomes adiabatic bubbles with radiative outer shell. We also determine analytically the onset of thin-shell formation time in case of fast wind driven bubble with power-law energy injection and ambient density structure. We solve the line transfer problem with numerical results in order to calculate line profile of [OIII] forbidden line.

• Structural Engineering and Mechanics
• /
• v.15 no.1
• /
• pp.135-149
• /
• 2003

The bubble packing method is implemented for adaptive mesh generation in two and three dimensions. Bubbles on the boundary of a three-dimensional domain are controlled independently of the interior bubbles in the domain, and a modified octree technique is employed to place initial bubbles in the three-dimensional zone. Numerical comparisons are made with other mesh generation techniques to demonstrate the effectiveness of the present bubble packing scheme for two- and three-dimensional domains. It is shown that this bubble packing method provides a high quality of mesh and affordable control of mesh density as well.

• Korean Journal of Mathematics
• /
• v.20 no.3
• /
• pp.323-331
• /
• 2012

We present a high-order potential ow model for the motion of the impulsively accelerated unstable interface of infinite density jump. The Layzer model for the evolution of the interface is extended to high-order. The time-evolution solutions of the bubble and the spike in the interface are obtained from the high-order model. We show that the high-order model gives improvement on the prediction of the evolution of the bubble and the spike.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.17 no.2 s.119
• /
• pp.168-176
• /
• 2007

The mixture sound speed in bubbly fluids is highly dispersive due to differences of the density and compressibility between bubbles and fluids. The dispersion range in bubbly fluids expands to a higher frequency than the resonance frequency of an air bubble. A theoretical model was developed to compute the reduction of radiation noise that is generated by a force applied on an infinite flat plate using a bubble layer as a compliant baffle. For evaluating the effectiveness of a bubble layer in reducing the structure-borne noise of an infinite elastic plate, the noise reduction levels for various parameters such as the thickness of bubble layers, the volume fractions and the distribution types of bubbly fluids are calculated numerically. The noise reduction effect of an air bubble layer on an infinite flat plate is considerable level and similar to the tendency of dispersion of bubbly fluids. It is recommended that the thickness of a bubble layer should be increased with keeping an appropriate volume fraction of an air bubble for the most effective reduction of the radiation noise.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.700-710
• /
• 2008

A Bubble size distribution model has been developed for the numerical simulation of cryogenic high-speed cavitating flow of the turbo-pumps in the liquid fuel rocket engine. The new model is based on the previous one proposed by the authors, in which the bubble number density was solved as a function of bubble size at each grid point of the calculation domain by means of Eulerian framework with respect to the bubble size coordinate. In the previous model, the growth/decay of bubbles due to pressure difference between bubble and liquid was solved exactly based on Rayleigh-Plesset equation. However, the unsteady heat transfer between liquid and bubble, which controls the evaporation/condensation rate, was approximated by a theoretical solution of unsteady heat conduction under a constant temperature difference. In the present study, the unsteady temperature field in the liquid around a bubble is also solved exactly in order to establish an accurate and efficient numerical simulation code for cavitating flows. The growth/decay of a single bubble and growth of bubbles with nucleation were successfully simulated by the proposed model.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.216-227
• /
• 2004

A new cavitating model by using bubble size distribution based on bubbles-mass has been proposed. Both liquid and vapor phases are treated with Eulerian framework as a mixture containing minute cavitating bubbles. In addition vapor phase consists of various sizes of vapor bubbles, which are distributed to classes based on their mass. The bubble number-density for each class was solved by considering the change of the bubble-mass due to phase change as well as generation of new bubbles due to heterogeneous nucleation. In this method, the bubble-mass is treated as an independent variable, and the other dependent variables are solved in spatial coordinates and bubble-mass coordinate. Firstly, we employed this method to calculate bubble nucleation and growth in stationary super-heated liquid nitrogen, and bubble collapse in stationary sub-cooled one. In the case of bubble growth in super-heated liquid, bubble number-density of the smallest class based on its mass is increased due to the nucleation. These new bubbles grow with time, and the bubbles shift to larger class. Therefore void fraction of each class is increased due to the growth in the whole class. On the other hand, in the case of bubble collapse in sub-cooled liquid, the existing bubbles are contracted, and then they shift to smaller class. It finally becomes extinct at the smallest one. Secondly, the present method is applied to a cavitating flow around NACA00l5 foil. Liquid nitrogen and liquid oxygen are employed as working fluids. Cavitation number, $\sigma$, is fixed at 0.15, inlet velocities are changed at 5, 10, 20 and 50m/s. Inlet temperatures are 90K in case of liquid nitrogen, and 90K and 1l0K in case of liquid oxygen. 110K of oxygen is corresponding to the 90K of nitrogen because of the same relative temperature to the critical one, $T_{r}$=$T/T_c^{+}$. Cavitating flow around the NACA0015 foils was properly analyzed by using bubble size distribution. Finally, the method is applied to a cavitating flow in an inducer of the LE-7A hydrogen turbo-pump. This inducer has 3 spiral foils. However, for simplicity, 2D calculation was carried out in an unrolled channel at 0.9R cross-section. The channel moves against the fluid at a peripheral velocity corresponding to the inducer revolutions. Total inlet pressure, $Pt_{in}$, is set at l00KPa, because cavitation is not generated at a design point, $Pt_{in}$=260KPa. The bubbles occur upstream of the foils and collapse between them. Cavitating flow in the inducer was successfully predicted by using the bubble size distribution.

• 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.

• Journal of the Korean Society of Visualization
• /
• v.15 no.3
• /
• pp.20-26
• /
• 2017

The two-fluid equations are widely used to simulate two-phase flows in a nuclear reactor. For the two-fluid momentum equation, the wall and interfacial drag terms play an important role in predicting a two-phase flow behavior. Since the bubble density is much smaller than the water density, the bubble accelerates faster than the liquid in a nozzle. As a result, the bubble phase becomes faster than the liquid phase in the nozzle. In contrast, the opposite phenomena occur in the diffuser. The purpose of our study is to experimentally show these behaviors in an area-varying channel such as nozzle and diffuser. Experiments were made of turbulent bubbly flows in an area-varying horizontal channel. The velocities of the bubble and liquid phases were measured by the PIV technique. It was shown that the two-phase velocities were no longer close to each other in the area-varying regions. The bubble was faster than the liquid in the nozzle; in contrast, the bubble was slower than the liquid in the diffuser. Code simulations were also performed using the MARS code. By replacing the original wall drag model in the MARS code with Kim (1)'s wall drag partition model, we obtained the simulation results being consistent with experimental observations.