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

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

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.411-418
• /
• 2004

Using a level set method we develop a shape optimization method applied to energy flow problems in steady state. The boundaries are implicitly represented by the level set function obtainable from the 'Hamilton-Jacobi type' equation with the 'Up-wind scheme.' The developed method defines a Lagrangian function for the constrained optimization. It minimizes a generalized compliance, satisfying the constraint of allowable volume through the variations of implicit boundary. During the optimization, the boundary velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian function. Compared with the established topology optimization method, the developed one has no numerical instability such as checkerboard problems and easy representation of topological shape variations.

• Journal of computational fluids engineering
• /
• v.7 no.4
• /
• pp.28-34
• /
• 2002

A general purpose program for computing 3-D flows has been extended for two-phase flows with topologically complex interfaces. The 3-D interfaces are tracked by employing a coupled level set and volume-of-fluid (CLSVOF) method which not only can calculate an interfacial curvature accurately but also can achieve mass conservation well. The program has been tested through the computations of bubbles rising in a liquid. The numerical results are found to compare well with the results reported in the literature.

• Journal of computational fluids engineering
• /
• v.6 no.3
• /
• pp.19-26
• /
• 2001

An Eulerian-Lagrangian method, so called immersed boundary method, is used for analysing viscous flow around arbitrary bodies, where governing equations are discretized on a regular grid by using a finite volume method. To improve the accuracy of flow near body boundaries, a second-order accurate interpolation scheme is used and a level-set based grid deformation method is presented to construct the adaptive grids around body boundaries. The present scheme is used to simulate steady flow around a semicircular cylinder mounted on the bottom of flow domain and calculated results are validated by results of a body fitted grid method. Finally, present method is applied to a complex flow around multi body and the usefulness is checked by investigating calculated results.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.448-450
• /
• 2010

Free surface liquid jet impingement, which is applicable to cooling of hot plates in a steel-making process, is investigated numerically by solving the conservation equations of mass, momentum and energy in the liquid and gas phases. The free-surface of liquid-gas interface is tracked by an improved level-set method incorporating a sharp-interface technique for accurate imposition of stress and heat flux conditions on the liquid-gas interface. The level-set approach is combined with a non-equilibrium $k-{\omega}$ turbulence model. The computations are made for slit nozzle jets to investigate their flow and cooling characteristics. Also, the effects of jetting angle, velocity and moving velocity of plate on the interfacial motion and the associated flow and temperature fields are quantified.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.340-346
• /
• 2008

A level-set method is developed for computation of drop motions in various engineering applications. Compared with the volume-of-fluid method based on a non-smooth volume-fraction function, the LS method can calculate an interface curvature more accurately by using a smooth distance function. Also, it is straightforward to implement for two-phase flows in complex geometries unlike the VOF method requiring much more complicated geometric calculations. The LS method is applied to simulation of inkjet process, thin film pattering and droplet collisions.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.221-226
• /
• 2008

Numerical simulation is performed for microdroplet deposition on a pre-patterned micro-structure. The level-set method for tracking the liquid-gas interface is extended to treat the immersed (or irregular-shaped) solid surface. The no-slip condition at the fluid-solid interface as well as the matching conditions at the liquid-gas interface is accurately imposed by incorporating the ghost fluid approach based on a sharp-interface representation. The method is further extended to treat the contact angle condition at an immersed solid surface. The present computation of a patterning process using microdroplet ejection demonstrates that the multiphase characteristics between the liquid-gas-solid phases can be used to improve the patterning accuracy.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.221-226
• /
• 2008

Numerical simulation is performed for microdroplet deposition on a pre-patterned micro-structure. The level-set method for tracking the liquid-gas interface is extended to treat the immersed (or irregular-shaped) solid surface. The no-slip condition at the fluid-solid interface as well as the matching conditions at the liquid-gas interface is accurately imposed by incorporating the ghost fluid approach based on a sharp-interface representation. The method is further extended to treat the contact angle condition at an immersed solid surface. The present computation of a patterning process using microdroplet ejection demonstrates that the multiphase characteristics between the liquid-gas-solid phases can be used to improve the patterning accuracy.