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

• Journal of the Korean Mathematical Society
• /
• v.37 no.5
• /
• pp.779-791
• /
• 2000

The nonlinear stage of the evolution of free boundary between a light fluid and a heavy fluid driven by an external force is studied by a potential flow model with a source singlarity. The potential flow model is applied to a bubble and spije evolution for constantly accelerated interface (Rayleigh-Taylor instability) and impulsively accelerated interface (Richtmyer-Meshkow instability). The numerical results of the model show that, in constantly accelerated intergace, bubble grows with constant velocity and the spike falls with gravitational acceleration at later times, while the velocity of the bubble in impulsively accelerated interface decay to zero asymp flow model for the bubble and spike for constantly accelerated interface and impulsively accelerated interface.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.126-134
• /
• 2008

We developed a three dimensional Navier-Stokes code based on the level set method to simulate two phase flows with high density ratio. The Navier-Stokes equations with consideration of the surface tension effects are solved by using SIMPLE algorithm on a non-staggered grid. The present code is validated by simulating two test problems. First one is to simulate a rising bubble inside a cube. The thickness of the interface of the bubble is shown to affect the pressure distribution around the interface. As the thickness decreases, the pressure field around the interface becomes more oscillatory. As the bubble rises, a ring vortex is shown to form around the interface and the bubble eventually develops into an ellipsoidal shape. Merge of two bubbles inside a container is secondly tested to show the robustness of the present code for two phase flow simulation. Numerical results show stable and reliable behavior during the process of merging of two bubbles. The velocity and pressure fields around the interface of bubbles are shown oscillation free during the merging of two bubbles.

• Journal of the Korean Ceramic Society
• /
• v.27 no.4
• /
• pp.543-547
• /
• 1990

Numerous researchers have observed the bubble fromation at SiC/SiO2 interface from 130$0^{\circ}C$ to 1$700^{\circ}C$. According to thermodynamic calculation, the bubble could be formed from the microscopic impurities which result from the chemical vapor deposition of SiC. When C-impurity is present at the interface, it is calculated that the bubble is formed at 1511$^{\circ}C$ and when Si is present, the bubble is formed at 177$0^{\circ}C$. These results are very close to the prior observations, but the calculation can not explain the observation of bubble below 150$0^{\circ}C$.

• Journal of Mechanical Science and Technology
• /
• v.15 no.7
• /
• pp.931-940
• /
• 2001

A numerical method for simulating bubble motion during nucleate boiling is presented. The vapor-liquid interface is captured by a level set method which can easily handle breaking and merging of the interface and can calculate an interfacial curvature more accurately than the VOF method using a step function. The level set method is modified to include the effects of phase change at the interface and contact angle at the wall as well as to achieve mass conservation during the whole calculation procedure. Also, a simplified model to predict the heat flux in a thin liquid microlayer is developed. The method is applied for simulation of a sliding bubble on a vertical surface to further understand the physics of partial boiling. Based on the computed results, the effects of contact angle, wall superheat and phase change on a sliding bubble are quantified.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1891-1896
• /
• 2003

The bubble motion during nucleate boiling in a microchannel is investigated numerically. The liquid-vapor interface is tracked by a level set method which is modified to include the effects of phase change at the interface and contact angle at the wall. The computations are made for various channel sizes, liquid flow rates, and contact angles. Based on the numerical results, the bubble growth pattern and its effect on the flow and heat transfer are discussed.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.69-74
• /
• 2001

It is well known that an air bubble trapped in water emits light at its collapse robustly with a proper forcing amplitude of ultrasound. Instability mechanism which causes deviation from sphericity of bubble wall was investigated theoretically. The rapid change of the bubble wall velocity which is both dependent on the forcing amplitude, was found to be a major factor of instability of the interface. The Rayleigh-Taylor instability which occurs when rapid acceleration is directed from the lighter towards the heavier fluid is found to be not related to the instability of the sonoluminescing gas bubble. A good agreement between the calculation results and experimental data is found.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.2
• /
• pp.265-273
• /
• 2003

Numerical analysis of bubble motion during nucleate boiling is performed by imposing a constant heat flux condition at the base of a heater which occurs in most of boiling experiments. The temporal and spatial variation of a solid surface temperature associated with the bubble growth and departure is investigated by solving a conjugate problem involving conduction in the solid. The vapor-liquid interface is tracked by a level set method which is modified to include the effects of phase change at the interface, contact angle at the wall and evaporative heat flux in a thin liquid micro-layer. Based on the numerical results, the bubble growth pattern and its interaction with the heating solid are discussed. Also, the effect of heating condition on the bubble growth under a micro-gravity condition is investigated.

• 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 B
• /
• v.28 no.8 s.227
• /
• pp.996-1003
• /
• 2004

The bubble motion during nucleate boiling in a microchannel is investigated by numerically solving the equations governing conservation of mass, momentum and energy in the liquid and vapor phases. The liquid-vapor interface is tracked by a level set method which is modified to include the effects of phase change at the interface and contact angle at the wall. Also, the evaporative heat flux from the thin liquid film that forms underneath a growing bubble attached to the wall is incorporated in the analysis. Based on the numerical results, the effects of channel size, contact angle, wall superheat and waiting period on the bubble growth and heat transfer in a microchannel are quantified.