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

• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1756-1768
• /
• 2021

The containment filtered venting system (CFVS) filters the atmosphere of the containment building and discharges a part of it to the outside environment to prevent containment overpressure during severe accidents. The Korean CFVS has a tank that filters fission products from the containment atmosphere by pool scrubbing, which is the primary decontamination process; however, prediction of its performance has been done based on researches conducted under mild conditions than those of severe accidents. Bubble behavior in a pool is a key parameter of pool scrubbing. Therefore, the bubble behavior in the pool was analyzed under various injection flow rates observed at the venturi nozzles used in the Korean CFVS using a wire-mesh sensor. Based on the experimental results, void fraction model was modified using the existing correlation, and a new bubble size prediction model was developed. The modified void fraction model agreed well with the obtained experimental data. However, the newly developed bubble size prediction model showed different results to those established in previous studies because the venturi nozzle diameter considered in this study was larger than those in previous studies. Therefore, this is the first model that reflects actual design of a venturi scrubbing nozzle.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.3E
• /
• pp.69-76
• /
• 2008

Measurements of bubble size and distribution in the surface layer of the sea, wind speed, and variation of ocean environments were made continually over a four-day period in an experiment conducted in the South Sea of Korea during 17-20 September 2007. Theoretical background of bubble population model indicates that bubble population is a function of the depth, range and wind speed and bubble effects on sound speed shows that sound speed varies with frequency. Observational evidence exhibited that the middle size bubble population fit the model very well, however, smaller ones can not follow the model probably due to their short lifetime. Meanwhile, there is also a hysteresis effect of void fraction. Observational evidence also indicates that strong changes in sound speed are produced by the presence of swarms of micro bubbles especially from 7 kHz to 50 kHz, and calculation results are consistent with the measured data in the high frequency band, but inconsistent in the low frequency band. Based on the measurements of the sound speed and high frequency transmission configuration in the bubble layer, we present an estimation of underwater acoustic channel capacity in the bubble layer.

• Journal of Energy Engineering
• /
• v.25 no.1
• /
• pp.69-75
• /
• 2016

Bubble size is a key parameter for an accurate prediction of bubble behaviours in the multi-dimensional two-phase flow. In the current STAR CCM+ CFD code, a mechanistic bubble size model $S{\gamma}$ is available for the prediction of bubble size in the flow channel. As another model, Yun model is developed based on DEBORA that is subcooled boiling data in high pressure. In this study, numerical simulation for the gas-liquid two-phase flow was conducted to validate and confirm the performance of $S{\gamma}$ model and Yun model, using the commercial CFD code STAR CCM+ ver. 10.02. For this, local bubble models was evaluated against the air-water data from DEDALE experiments (1995) and Hibiki et al. (2001) in the vertical pipe. All numerical results of $S{\gamma}$ model predicted reasonably the two-phase flow parameters and Yun model is needed to be improved for the prediction of air-water flow under low pressure condition.

• The Journal of the Acoustical Society of Korea
• /
• v.40 no.4
• /
• pp.320-329
• /
• 2021

In this paper, an empirical model of air bubble size to be applied to an air masker for reduction of underwater radiation noise is presented. The proposed model improves the divergence problem under the low-speed flow condition of the existing model derived using Rayleigh's jet instability model and simple continuity condition by introducing a jet flow velocity of air. The jet flow velocity of air is estimated using the bubble size where the liquid is quiescent. In a medium without flow, the size of the bubble is estimated by an empirical method where bubble formation regime is divided into a laminar-flow range, a transition range, and a turbulent-flow range based on the Reynolds number of the injected air. The proposed bubble size model is confirmed to be in good agreement with the Computational Fluid Dynamics (CFD) analysis result and the experimental results of the existing literature. Using the acoustic inversion method, the air bubble population is estimated from the insertion loss measured during the air injection experiment of the air- masker model in a large cavitation tunnel. The results of the experiments and the bubble size model are compared in the paper.

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

• Nuclear Engineering and Technology
• /
• v.56 no.2
• /
• pp.624-635
• /
• 2024

In two-fluid simulations of flow boiling, the modeling of the mean bubble diameter is a key parameter in the closure relations governing the intefacial transfer of mass, momentum, and energy. Monodispersed approach proved to be insufficient to describe the significant variation in bubble size during flow boiling in a heated pipe. A population balance model (PBM) has been employed to address these shortcomings. During nucleate boiling, vapor bubbles of a certain size are formed on the heated wall, detach and migrate into the bulk flow. These bubbles then grow, shrink or disintegrate by evaporation, condensation, breakage and aggregation. In this study, a parametric analysis of the PBM aggregation and breakage models has been performed to investigate their effect on the radial distribution of the mean bubble diameter and vapor volume fraction. The simulation results are compared with the DEBORA experiments (Garnier et al., 2001). In addition, the influence of PBM parameters on the local distribution of individual bubble size groups was also studied. The results have shown that the modeling of aggregation process has the largest influence on the results and is mainly dictated by the collisions due to flow turbulence.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1365-1373
• /
• 1994

A Monte-Carlo simulation model, developed to predict size distribution of air bubbles in turbulent shear flow, is applied to a laboratory-scale problem. Sensitivity to various numerical and physical parameters of the model is analyzed. Practical applicability of the model is explored through comparisons of results with experimental measurements. Bubble size increases with air-water discharge ratio and friction factor. Bubble size decreases with increasing mean flow velocity, but the total bubble surface area in the aeration region remains fairly constant. The effect on bubble size distribution of the longitudinal length increment in the simulation model is negligible. A larger radial length increment yields more small and large bubbles and fewer in between. Bubble size distribution is significantly affected by its initial distribution and the location of air injection. Collision efficiency is introduced to explain the discrepancy between collisions with and without coalescence.

• Journal of Mechanical Science and Technology
• /
• v.15 no.12
• /
• pp.1801-1807
• /
• 2001

In multiphase flne the bubble size and velocity. To achieve this, one of approaches is to utilize laser phase-Doppler anemometry. However, it was found that the second order refraction has great impact on PDA sizing method when the relative refractive index of media is less than one. In this paper, the problem of second order refraction is investigated and a model of phase-size correlation to eliminate the measurement errors is introduced for bubble sizing. As a result, the model relates the assumption of single scattering mechanism in conventional phase-Doppler anemometry. The results of simulations based on this new model by using Generalized Lorenz Mie Theory (GLMT) are compared with those based on the conventional method. An optimization method for accurately sizing air-bubble in water has been suggested.