• 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.15 no.4E
• /
• pp.3-8
• /
• 1996

The bubble size distribution is critical information to understand sound propagation and ambient noise in the ocean. To estimate the bubble size distribution in a bubbly water, the sound attenuation has been only in the conventional acoustic bubble sizing method without considering the sound speed variation. However, the effect of the sound speed variation in bubbly water cannot be neglected because of its compressibility variation. The sound attenuation is also affected by the sound speed variation. In this paper, a coherent acoustic bubble sizing inverse technique is introduced as a new bubble sizing technique with considering sound speed variation as well as the sound attenuation. This coherent sizing method is theoretically verified with the bubble distribution functions of single-size, Gaussian, and power-law functions. Its numerical test results with the coherent acoustic bubble sizing method show good agreement with the given bubble distributions.

• Korean Journal of Hazardous Materials
• /
• v.6 no.2
• /
• pp.95-104
• /
• 2018

Gas-liquid bubble column reactors are extensively used in industrial processes. A detailed knowledge of bubble size distribution is needed for determining the mass transfer in gas-liquid film. Experimental data on bubble size distribution and liquid-side mass transfer coefficient($k_L$) were used to calculate the estimated time to saturation in bubble column reactor. Also, the gas flux was evaluated to the liquid-side mass transfer coefficient($k_L$) and solubility data for hydrogen sulfide($H_2S$) and chlorine($Cl_2$) absorption into water. Simulation results show that $H_2S$ absorption time to 50 % of saturation concentrations are 611 sec and 1,329 sec when bubble diameters are 0.5 mm and 4.5 mm, while absorbing 1 % $H_2S$ gas. In case of $Cl_2$, absorption time range 657 to 1,400 sec when bubble size range 0.5 mm to 4.5 mm, while absorbing 1 % $Cl_2$ gas. Calculated simulation results can be used in the design of emergency relief bubble reactors.

• International Journal of Fluid Machinery and Systems
• /
• v.8 no.2
• /
• pp.73-84
• /
• 2015

Aerating hydroturbines have recently been proposed as an effective way to mitigate the problem of low dissolved oxygen in the discharge of hydroelectric power plants. The design of such a hydroturbine requires a precise understanding of the dependence of the generated bubble size distribution upon the operating conditions (viz. liquid velocity, air ventilation rate, hydrofoil configuration, etc.) and the consequent rise in dissolved oxygen in the downstream water. The purpose of the current research is to investigate the effect of location of air injection on the resulting bubble size distribution, thus leading to a quantitative analysis of aeration statistics and capabilities for two turbine blade hydrofoil designs. The two blade designs differed in their location of air injection. Extensive sets of experiments were conducted by varying the liquid velocity, aeration rate and the hydrofoil angle of attack, to characterize the resulting bubble size distribution. Using a shadow imaging technique to capture the bubble images in the wake and an in-house developed image analysis algorithm, it was found that the hydrofoil with leading edge ventilation produced smaller size bubbles as compared to the hydrofoil being ventilated at the trailing edge.

• Journal of the Korean Society of Visualization
• /
• v.17 no.1
• /
• pp.60-68
• /
• 2019

Flow visualizations have been carried out to analyze the characteristics of bubby flows exhausted from a venturi-type bubble generator with an air vent. For various design parameters and operating conditions of the bubble generator, the images of bubbly flows was recorded using a high-speed camera and a microscope. Then the amount and size distribution of bubble was evaluated by an image processing technique. The results show that for increasing the amount of bubble, it is more effective to reduce the venturi throat than to enlarge the air vent diameter. If the water flow rate increases, the bubble generation rate increases but reaches a status of saturation, whose condition depends on Reynolds number at a given air vent diameter. The bubble size increases as the diameter of venturi throat decreases and Reynolds number increases. However, the air vent diameter is not a significant factor on bubble size.

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

• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.509-521
• /
• 2021

In this paper, from a review of the size distribution of the bubbles during pool scrubbing obtained from experiments by EPRI, we apply the bubble size distributions to analyses on the decontamination factors of pool scrubbing via I-COSTA (In-Containment Source Term Analysis). We perform sensitivity studies of the bubble size on the various mechanisms of deposition of aerosol particles in pool scrubbing. We also perform sensitivity studies on the size distributions of the bubbles depending on the diameters at the nozzle exit, the molecular weights of non-condensable gases in the carrier gases, and the steam fractions of the carrier gases. We then perform analyses of LACE-ESPANA experiments and compare the numerical ~ results to those from SPARC-90 and experimental results in order to show the effect of the bubble size distributions.

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

• The Journal of the Acoustical Society of Korea
• /
• v.42 no.4
• /
• pp.305-312
• /
• 2023

Physics-Informed Neural Network (PINN) is used to invert bubble size distributions from attenuation losses. By considering a linear system for the bubble population inversion, Adaptive Learned Iterative Shrinkage Thresholding Algorithm (Ada-LISTA), which has been solved linear systems in image processing, is used as a neural network architecture in PINN. Furthermore, a regularization based on the linear system is added to a loss function of PINN and it makes a PINN have better generalization by a solution satisfying the bubble physics. To evaluate an uncertainty of bubble estimation, deep ensemble is adopted. 20 Ada-LISTAs with different initial values are trained using the same training dataset. During test with attenuation losses different from those in the training dataset, the bubble size distribution and corresponding uncertainty are indicated by average and variance of 20 estimations, respectively. Deep ensemble Ada-LISTA demonstrate superior performance in inverting bubble size distributions than the conventional convex optimization solver of CVX.