• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.389-396
• /
• 2001

Direct numerical simulation of bubble growth and merger process on a single nucleation site during partial nucleate boiling is performed. The equations governing conservation of mass, momentum and energy are solved using a finite difference method combined with a level set method for capturing the vapor-liquid interface. The level set method is modified to include the effects of phase change at the interface and contact angle at the wall. Also, a simplified formulation for predicting the evaporative heat flux in a thin liquid micro-layer is developed and incorporated into the level set formulation. Based on the numerical results, the bubble growth and merger pattern and its effect on the heat transfer are discussed.

• 한국연소학회:학술대회논문집
• /
• 2004.06a
• /
• pp.102-111
• /
• 2004

Partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics to develope a prediction model for the turbulent lift off. The present study is specifically aimed to remedy the problem of the stiff transition of the conditioned partial burning probability across the crossover condition by adopting level-set method which describes propagating or retreating flame front with specified propagation speed. In light of the level-set simulations with two model problems for the propagation speed, the stabilizing conditions for a turbulent lifted flame are suggested. The flame hole dynamics combined with level-set method yields a temporally evolving turbulent extinction process and its partial quenching characteristics is compared with the results of the previous model employing the flame-hole random walk mapping. The probability to encounter reacting' state, conditioned with scalar dissipation rate, demonstrated that the conditional probability has a rather gradual transition across the crossover scalar dissipation rate in contrast to the stiff transition of resulted from the flame-hole random walk mapping and could be attributed to the finite response of the flame edge propagation.

• Journal of the Korean Society of Combustion
• /
• v.9 no.2
• /
• pp.18-29
• /
• 2004

Partial quenching structure of diffusion flames in a turbulent mixing layer has been investigated by the method of flame hole dynamics in oder to develope a prediction model for the phenomenon of turbulent flame lift off. The present study is specifically aimed to remedy the shortcoming of the stiff transition of the conditioned partial burning probability across the crossover condition by employing the level-set method which enables us to include the effect of finite flame edge propagation speed. In light of the level-set simulation results with two models for the edge propagation speed, the stabilizing conditions for turbulent lifted flame are suggested. The flame hole dynamics combined with the level-set method yields a temporally evolving turbulent extinction process and its partial quenching characteristics is compared with the results of the previous model employing the flame-hole random walk mapping based on three critical scalar dissipation rates. The probability to encounter reacting state, conditioned with scalar dissipation rate, demonstrated that the conditional probability has a rather gradual transition across the crossover scalar dissipation rate. Such a smooth transition is attributed to the finite response of the flame edge propagation.

• Proceedings of the KIEE Conference
• /
• 2008.07a
• /
• pp.623-624
• /
• 2008

We present a level set method for numerical shape optimization of electromagnetic systems. The level set method does not only lead to efficient computational schemes, but also is able to handle topological changes such as merging, splitting and even disappearing of connected components. The velocity field on boundaries is obtained by a shape derivative of continuum sensitivity analysis using the material derivative concept and an adjoint variable technique. Two numerical results of dielectric optimization between electrodes showed that the level set method is feasible and effective in solving shape optimization problems of electromagnetic systems.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.4
• /
• pp.36-41
• /
• 2011

In the present paper, a direct forcing/fictitious domain (DF/FD) level set method is proposed to simulate the FSI (fluid-solid interaction) in two-phase flow. The main idea is to combine the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. The DF/FD method is a non-Lagrange-multiplier version of a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. This method does not sacrifice the accuracy and robustness by employing a discrete ${\delta}$ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, easy implementation, and utilization of the original governing equation without modification. Simulations of various water-entry problems have been conducted to validate the capability and accuracy of the present method in solving the FSI in two-phase flow. Consequently, the present results are found to be in good agreement with those of previous studies.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.155-158
• /
• 2009

We developed a numerical method to analyze the contact-line problems, incorporating a dynamic contact angle model. We used level set method to capture free surface. The method is applied to the analysis of dynamic behavior of a droplet in DC electrowetting. The result is compared with an experimental data and result of perturbation method.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.223-227
• /
• 2009

In the present study, a least square weighted residual method and Taylor-Galerkin method were formulated and tested for the discretization of the two hyperbolic type equations of level set method; advection and reinitialization equations. The two approaches were compared by solving a time reversed vortex flow and three-dimensional broken dam flow by employing a four-step splitting finite element method for the solution of the incompressible Navier-Stokes equations. From the numerical experiments, it was shown that the least square method is more accurate and conservative than Taylor-Galerkin method and both methods are approximately first order accurate when both advection and reinitialization phase are involved in the evolution of free surface.

• Journal of Ship and Ocean Technology
• /
• v.6 no.2
• /
• pp.37-50
• /
• 2002

In the present numerical simulation of viscous free surface flow around a ship, two-fluids in-compressible Reynolds-averaged Navier-Stokes equations with the standard $\textsc{k}-\varepsilon$turbulence model are discretized on a regular grid by using a finite volume method. A local level-set method is introduced for capturing the free surface movement and the influence of the viscous layer and dynamic boundary condition of the free surface are implicitly considered. Partial differential equations in the level-set method are discretized with second order ENO scheme and explicit Euler scheme in the space and time integration, respectively. The computational results for the Series-60 model with $C_B=0.6$ show a good agreement with the experimental data, but more validation studies for commercial complicated hull forms are necessary.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.3
• /
• pp.158-165
• /
• 2005

In order to analyze incompressible two-phase flow, Level Set method was applied on pseudo-compressibility formulation. Level Set function is defined as a signed distance function from the phase interface, and gives the information of the each phase location and the geometric data to the flow. In this study, Level Set function transport equation was coupled with flow conservation equations, and owing to pseudo-compressibility technique we could solve the resultant vector equation iteratively. Two-phase flow analysis code was developed on general curvilinear coordinate, and numerical tests of bubble dynamics and surging wave problems demonstrate its capability successfully.