• Fibers and Polymers
• /
• v.6 no.2
• /
• pp.131-138
• /
• 2005

A theoretical model is proposed in order to investigate rheological behavior of bubble suspension with large deformation. Theoretical constitutive equations for dilute bubble suspensions are derived by applying a deformation theory of ellipsoidal droplet [1] to a phenomenological suspension theory [2]. The rate of deformation tensor within the bubble and the time evolution of interface tensor are predicted by applying the proposed constitutive equations, which have two free fitting parameters. The transient and steady rheological properties of dilute bubble suspensions are studied for several capillary numbers (Ca) under simple shear flow and uniaxial elongational flow fields. The retraction force of the bubble caused by the interfacial tension increases as bubbles undergo deformation. The transient and steady relative viscosity decreases as Ca increases. The normal stress difference (NSD) under the simple shear has the largest value when Ca is around 1 and the ratio Of the first NSD to the second NSD has the value of 3/4 for large Ca but 2 for small Ca. In the uniaxial elongational flow, the elongational viscosity is three times as large as the shear viscosity like the Newtonian fluid.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.6
• /
• pp.1518-1528
• /
• 1995

The deformation and departure processes of a bubble attached to the wall are studied experimentally and numerically to understand the phenomena of the nucleate boiling heat transfer enhancement under DC electric field. An air-bubble is injected in a dielectric liquid with different electric fields generated by changing three types of electrode system (Type 1,2 and 3) in the bubble generator. Experimental variables are the electric field strength and the distance and the shape of the electrodes under DC electric field. From experimental results, it is observed that the bubble under Dc electric field is elongated in the same direction as the electric field and the contact angle increases. For the parallel plate electrode which generates a uniform electric field, bubble departure volume doesn't seem to decrease within our experimental range. However, when a needle is raised a few millimeters from the lower electrode to make a nonuniform electric field around the needle, bubble departure volume decreases continuously with the increase of an applied voltage. The reduction effect of bubble departure volume is the most effective under a strong nonuniform electric field generated with Type 3. As the nonuniformity of the electric field due to the shape of a electrode increases, the terminal velocity and the acceleration of a bubble increase largely. For the comparison with visualization results, the deformation of a bubble attached to the electrode is carried out by a numerical method. Numerical results show good agreement qualitatively with experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.6
• /
• pp.2023-2035
• /
• 1996

In the present study, the characteristics of the electrical deformation of a bubble and a drop under a uniform electric field have been investigated to understand EHD heat transfer enhancement by an electric field. The deformation of the bubble and the drop have been studied theoretically using an electric normal stress acting on their interfaces and assured by the numerical analysis and the experiment. From the variation of bubble volume and free energy, it is found that a bubble is compressed in an electric field and free energy had larger value with increasing W and the permittivity of a dielectric fluid. The electric normal stress induced on the interface of the bubble and the drop is different. Because of the surface charge induced at the drop interface, the electric normal stress acting on the drop is much larger than that of the bubble. The drop is, therefore, deformed much more than the bubble. In addition, the experimental and numerical results show that the aspect ratio and the contact angle of the bubble increase with increasing W.

• Korean Chemical Engineering Research
• /
• v.48 no.4
• /
• pp.462-469
• /
• 2010

Gas-liquid 2 phase processes are usually used in chemical, biochemical, environmental engineering and food process. For optimizing these processes, understanding bubble's precise movement and shape are needed. Bubble's movement and shape are effected by liquid's properties-viscosity, surface tension and bubble's properties-size, velocity. This paper deals with experimental data of bubble's movement and shape in high viscous silicone oil. Also, drag coefficient and deformation factor given by other researcher's papers and books are used to predicting and comparing bubble's terminal velocity, drag coefficient, deformation factor and shape with experimental value. Experimental data show that bubble moves faster when it moves in lower viscous silicone oil and it's drag coefficient is bigger when it moves in high viscous silicone oil. Bubble's shape is close to sphere when moving in high viscous silicone. Formulas proposed by Batchelor expect most accurate prediction for bubble's velocity and drag coefficient. Bubble's 2D shape predicted by Batchelor's energy balance, drag coefficient and deformation factor show excellent agreement with experimental bubble's 2D shape.

• Journal of the Korean Society of Visualization
• /
• v.18 no.3
• /
• pp.35-41
• /
• 2020

We conduct numerical simulations of the interaction of a deformable structure with two-phase compressible flow. The finite volume method (FVM) is used to simulate fluid phenomena including a shock wave, a gas bubble, and the deformation of free surface. The deformation of a floating structure is computed with the finite element method (FEM). The compressible two-phase volume of fluid (VOF) method is used for the generation and development of a cavitation bubble, and the immersed boundary method (IBM) is used to impose the effect of the structure on the fluid domain. The result of the simulation shows the generation of a shock wave, and the expansion of the bubble. Also, the deformation of the structure due to the hydrodynamic loading by the explosion is identified.

• Proceedings of the Korean Fiber Society Conference
• /
• 2003.10a
• /
• pp.47-48
• /
• 2003

Dilute bubble suspensions are prepared by introducing carbon dioxide bubbles into polyol resin. The apparent shear viscosity is measured with a wide gap parallel plate rheometer. A numerical simulation for deformation of a single bubble suspended in a Newtonian fluid is conducted by using a finite volume method (FVM) where multigrid algorithms are incorporated. Transient and steady results of bubble deformation were obtained and were in good agreement with experimental results. At high capillary number, viscosity of the suspension increases as the volume fraction increases, while at low capillary number, the viscosity decreases as the volume fraction increases.

• Korea-Australia Rheology Journal
• /
• v.16 no.1
• /
• pp.47-54
• /
• 2004

Low Reynolds number, dilute, and surfactant-free bubble suspensions are prepared by mechanical mixing after introducing carbon dioxide bubbles into a Newtonian liquid, polyol. The apparent shear viscosity is measured with a wide-gap parallel plate rheometer by imposing a simple shear flow of capillary numbers(Ca) of the order of $10^{-2}$ ~ $10^{-1}$ and for various gas volume fractions ($\phi$). Effects of capillary numbers and gas volume fractions on the viscosity of polyol foam are investigated. At high capillary number, viscosity of the suspension increases as the gas volume fraction increases, while at low capillary number, the viscosity decreases as the gas volume fraction increases. An empirical constitutive equation that is similar to the Frankel and Acrivos equation is proposed by fitting experimental data. A numerical simulation for deformation of a single bubble suspended in a Newtonian fluid is conducted by using a newly developed two-dimensional numerical code using a finite volume method (FVM). Although the bubble is treated by a circular cylinder in the two dimensional analysis, numerical results are in good agreement with experimental results.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.2
• /
• pp.165-171
• /
• 2009

We developed a 'multi-point vibration acceleration method' for accurately predicting the cavitation intensity in pumps. Pressure wave generated by cavitation bubble collapse propagates and causes pump vibration. We measured vibration accelerations at several points on a casing, suction and discharge pipes of centrifugal and mixed-flow pumps. The measured vibration accelerations scattered because the pressure wave damped differently between the bubble collapse location and each sensor. In a conventional method, experimental constants are proposed without evaluating pressure propagation paths, then, the scattered vibration accelerations cause the inaccurate cavitation intensity. In our method, we formulated damping rate, transmittance of the pressure wave, and energy conversion from the pressure wave to the vibration along assumed pressure propagation paths. In the formulation, we theoretically defined a 'pressure propagation coefficient,' which is a correlation coefficient between the vibration acceleration and the bubble collapse pressure. With the pressure propagation coefficient, we can predict the cavitation intensity without experimental constants as proposed in a conventional method. The prediction accuracy of cavitation intensity is improved based on a statistical analysis of the multi-point vibration accelerations. The predicted cavitation intensity was verified with the plastic deformation rate of an aluminum sheet in the cavitation erosion area of the impeller blade. The cavitation intensities were proportional to the measured plastic deformation rates for three kinds of pumps. This suggests that our method is effective for estimating the cavitation intensity in pumps. We can make a cavitation intensity map by conducting this method and varying the flow rate and the net positive suction head (NPSH). The map is useful for avoiding the operating conditions having high risk of cavitation erosion.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.1926-1932
• /
• 2001

In this second part of the paper, the automatic mesh generation and remeshing algorithm using bubble packing method is applied to the nonlinear problem. The remeshing/refinement procedure is necessary in the large deformation process especially because the mesh distortion deteriorates the convergence and accuracy. To perform the nonliear analysis, the transfer of state variables such as displacement and strain is added to the algorithm of Part 1. The equilibrium equation based on total Lagrangian formulation and elasto-viscoplastic model is used. For the numerical experiment, the upsetting process including the contact constraint condition is analyzed by two refinement criteria. And from the result, it is addressed that the present algorithm can generate the refined meshes easily at the largely deformed area with high error.