• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.566-574
• /
• 2021

The process of cavitation involves generation, growth, coalescence, and collapse of small bubbles and is tremendously influenced by bubble-bubble interactions. To understand these interactions, a new cavitation model based on the transport equation is proposed herein. The modified Rayleigh-Plesset equation is analyzed to determine the bubble growth rate by assuming equal-sized spherical bubble clouds. The source term in the transport equation is then derived according to the bubble growth rate with the bubble-bubble interaction. The proposed model is validated by various test simulations, including microscopic bubble cloud evolution as well as macroscopical two- and three-dimensional cavitating flows. Compared with previous models, namely the Kunz and Zwart cavitation models, the newly proposed model does not require adjustable parameters and generally results in better predictions both microscopic and macroscopical cases. This model is more physical.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.963-970
• /
• 2007

The objectives of this study are to examine a bubble trap mechanism of the turbidimeter for low turbidity and to acquire it's technology concerned. Reynolds-Averaged Wavier-Stokes equation and Laglangian discrete phase model were applied to analyze a flow field in the bubble trap. 3D hybrid grid system was used to simulate the flow field of bubble trap and numbers of it's node point are about 110,000. From the comparison between the standard $k-{\varepsilon}$ model and the laminar state, it was found that the former estimated less the velocity in the outlet of bubble trap than the latter did, and that the former estimated more the shear stress at the wall of bubble trap than the latter did. And, it was possible to visualize the path of bubbles in the bubble trap and to copy the removal process of bubbles out bubble trap. Also, it was found that nearly most of bubbles in the bubble trap disappeared.

• Korean Chemical Engineering Research
• /
• v.55 no.5
• /
• pp.698-703
• /
• 2017

Bubble flow characteristics in fluidized beds were analyzed by CPFD simulation. A fluidized bed, which had the size of $0.3m-ID{\times}2.4m-high$, was modeled by commercial CPFD $Barracuda^{(R)}$. Properties of bed material were $d_p=150{\mu}m$, ${\rho}_p=2,330kg/m^3$, and $U_{mf}=0.02m/s$. Gas was uniformly distributed and the range of superficial gas velocity was 0.07 to 0.16 m/s. Two other geometries were modeled. The first was a three-dimensional model, and the other was a two-dimensional model of $0.01m{\times}0.3m{\times}2.4m$. Bubble size and rising velocity were simulated by axial and radial position according to superficial gas velocity. In the case of three-dimensional model, simulated bubble rising velocity was different from correlations, because there was zigzag motion in bubble flow, and bubble detection was duplicated. To exclude zigzag motion of bubble flow, bubble rising velocity was simulated in the two-dimensional model and compared to the result from three-dimensional model.

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

• Polymer(Korea)
• /
• v.28 no.1
• /
• pp.51-58
• /
• 2004

The molecular cluster model for the homogeneous bubble nucleation rather than the classical nucleation theory was extended to predict the bubble nucleation events in elastomers(cross-linked polymers), polymers and polymer which are dissolved in the organic solvent. The classical theory assumes the formation of the critical bubble while the molecular cluster model assumes the critical cluster as for the initiation of the bubble nucleation. For the bubble nucleation in elastomers and polymers, the strain energy overcome by a critical bubble was also considered. The calculation results for the number of bubbles nucleated in elastormers and polymer solutions, which are about 10$\^$8/∼10$\^$12/ bubbles/㎤ are in good agreement with observed ones.

• Journal of the Korean Mathematical Society
• /
• v.54 no.2
• /
• pp.663-673
• /
• 2017

We consider the multi-harmonic model for the bubble evolution in the Rayleigh-Taylor instability in two and three dimensions. We extend the multi-harmonic model in two dimensions to a high-order and present a new class of steady-state solutions of the bubble motion. The growth rate of the bubble is expressed by a continuous family of two free parameters. The critical point in the family of solutions is identified as a saddle point and is chosen as the physically significant solution. We also present the multi-harmonic model in the cylindrical geometry and find the steady-state solution of the axisymmetric bubble. Validity and limitation of the model are also discussed.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1228-1233
• /
• 2003

The molecular cluster model for the homogeneous bubble nucleation was extended to predict the bubble nucleation events in elastomers, polymers and polymer solutions. For the bubble nucleation in elastomers and polymers, the strain energy overcome by a critical bubble was also considered. The calculation results for the number of bubbles nucleated are in good agreement with observed ones.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.1
• /
• pp.55-61
• /
• 2011

This paper presents bubble wake measurement results generated by the planing hull. The bubble was generated by SNAME TMB model(No. 4876) with hard chine at the CWC(Circulating Water Channel). ABS(Acoustic Bubble Spectrometer) was used to measure bubble wake measurement. The manufactured model is one meter in length and uniform velocity to generate the bubble at CWC is 3m/s, relatively higher speed than conventional hull form. Measurements were performed successfully and measured results show well the general characteristics of bubble wake generated by planing hull. Furthermore, experimental equations are proposed for the practical use.