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

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

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

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

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.10a
• /
• pp.599-604
• /
• 1997

One of important inconsistencies in the six-equation model predictions has been found to be the force experienced by a single bubble placed in a convergent stream of liquid. Various sets of governing equations yield different amount of forces to hold the bubble stationary in a convergent nozzle. By using the first order potential flow theory, it is found that the six-equation model can not be used to estimate the force experienced by a deformed bubble. The theoretical value of the particle stress of a bubble in a convergent nozzle flow has been found to be a function of the Weber number when bubble distortion is allowed. This force has been calculated by using different sets of governing equations and compared with the theoretical value. It is suggested in this study that the bubble size distribution function can be used to remove the presented inconsistency by relating the interfacial variables with different moments of the bubble size distribution function. This study also shows that the inconsistencies in the thermal-hydraulic governing equation can be removed by mechanistic modeling of the phasic interface.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1575-1588
• /
• 2019

Fracture near the U-10Mo/cladding material interface impacts fuel service life. In this work, a mesoscale stress model is developed with the fuel foil considered as a porous medium having gas bubbles and bearing bubble pressure and surface tension. The models for the evolution of bubble volume fraction, size and internal pressure are also obtained. For a U-10Mo/Al monolithic fuel plate under location-dependent irradiation, the finite element simulation of the thermo-mechanical coupling behavior is implemented to obtain the bubble distribution and evolution behavior together with their effects on the mesoscale stresses. The numerical simulation results indicate that higher macroscale tensile stresses appear close to the locations with the maximum increments of fuel foil thickness, which is intensively related to irradiation creep deformations. The maximum mesoscale tensile stress is more than 2 times of the macroscale one on the irradiation time of 98 days, which results from the contributions of considerable volume fraction and internal pressure of bubbles. This study lays a foundation for the fracture mechanism analysis and development of a fracture criterion for U-10Mo monolithic fuels.

• Journal of Korean Society of Water and Wastewater
• /
• v.18 no.2
• /
• pp.235-241
• /
• 2004

The dissolved air flotation (DAF) process has been widely used for removing suspended solids with low density in water. It has been known as measuring the size of microbubbles precisely which move upward rapidly in contact zone is difficult. In this study particle counter monitoring (PCM) method is used to measure the rising microbubble after injection from a nozzle. Size and distribution curve of microbubbles are evaluated at different conditions such as pressure drop at intermediate valve, length of pipeline between saturation tank and nozzle and low pressure. And the efficiency is also checked when it collides with different size floc. The experimental results show the following fact. As the final pressure drop occurred closer to a nozzle, the bubble size became smaller. And small bubble collides with large floc as well as small one because of its physical characteristic. However large bubble collides well with large floc rather than small one since hydrodynamic flow in streamline interferes to collide between two. With performing computational process by mathematical model we have analyzed and verified the size effect between bubble and floc. Collision efficiency is the highest when P/B ratio shows in the range of 0.75 < P/B ratio ($R_{particle/Rbubble}$) < 2.0.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.226-233
• /
• 2022

In the present work, a phase-field model was developed to investigate the influence of recrystallization on bubble evolution during irradiation. Considering the interaction between bubbles and grain boundary (GB), a set of modified Cahn-Hilliard and Allen-Cahn equations, with field variables and order parameters evolving in space and time, was used in this model. Both the kinetics of recrystallization characterized in experiments and point defects generated during cascade were incorporated in the model. The bubble evolution in recrystallized polycrystalline of U-Mo alloy was also investigated. The simulation results showed that GB with a large area fraction generated by recrystallization accelerates the formation and growth of bubbles. With the formation of new grains, gas atoms are swept and collected by GBs. The simulation results of bubble size and distribution are consistent with the experimental results.