• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.21-25
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Fire Science and Engineering
• /
• v.30 no.6
• /
• pp.111-117
• /
• 2016

The present study examined the free surface flow of a liquid sheet near a sprinkler head using a Computational Fluid Dynamics (CFD) model and considered the feasibility of the empirical model for predicting the initial spray characteristics of the sprinkler head through a comparison of the CFD results. The CFD calculation for a simplified sprinkler geometry considering the nozzle and deflector were performed using the commercially available CFD package, CFX 14.0 with the standard $k-{\varepsilon}$ turbulence model and theVolume of Fluid (VOF) method. The predicted velocity of the empirical model at the edge of deflector were in good agreement with that of the CFD model for the flat plate region but there was a certain discrepancy between the two models for the complex geometry region. The mean droplet diameter predicted by the empirical model differed significantly from the measured value of the real sprinkler head. On the other hand, the empirical model can be used to understand the mechanism of droplet formation near the sprinkler head and predict the initial spray characteristics for cases without experimental data.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.6
• /
• pp.750-757
• /
• 2016

This study presents the wave run-up height around single and multiple surface-piercing cylinders according to wave period and steepness. In order to simulate 3D incompressible viscous two-phase turbulent flow, the present study employed a volume of fluid (VOF) method with realizable $k-{\varepsilon}$ turbulence model based on commercial Computational Fluid Dynamics (CFD) software, "STAR-CCM". The wave periods at model scale were 1.269s and 1.692s for a single cylinder and 1.716s for multiple cylinders. In each case, wave steepness of has 1/30 and 1/16 were used, respectively. Consequently, the results for wave run-up height with regard to wave steepness and period were compared with those of relevant previous experimental studies. The numerical simulation results showed a good qualitative agreement with experiments.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.41-49
• /
• 2006

A modified Volume-of-Fluid (VOF) numerical method is used to predict the dynamics of a liquid drop of a low viscosity dilute polymer solution, forming in air from a circular nozzle. Viscoelastic effects are rep-resented using an Oldroyd-B model. Predicted drop shapes are compared with experimental observations. The main features, including the timing of the shape evolution and the 'bead-on-a-string' effect, are well reproduced by the simulations. The results confirm published conclusions of the third author, that the deformation is effectively Newtonian until near the time of Newtonian pinch-off and that the elastic stress becomes large in the pinch region due to the higher extensional flow there.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.704-708
• /
• 2017

A spray characteristics is carried out in a numerical simulation of swirl coaxial injector. The water and nitrogen are the oxidizer and fuel is used in cold flow condition. The simulation is carried out in 3d model with varying recess length. Reynolds stress turbulence and volume of fluid model were chosen to perform the simulation. The spray characteristics have been investigated as well as the influence of the inlet swirl strength of the internal flow. Effect of recess length is studied for the axial and radial velocity decreased with a reduced length of inner injector due to the decline vortex intensity.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.341-346
• /
• 2004

Motion of a bubble inside narrow tube is numerically studied. The numerical code assumes axi-symmetric incompressible flow field. The surface of the bubble is captured by VOF (Volume Of Fluid) method, and it is advected by MARS (Multiphase Advection and Reconstruction Scheme). Air bubble inside water is first studied, and it was found that a strong vortex, which is induced by the pressure difference caused by the surface tension, is formed at the rear part of the bubble. Then flow parameters are parametrically varied to understand the correlation between the bubble shape, the bubble velocity, and the flow parameters. The parametric study revealed that the aspect ratio of the bubble mainly depends on We number, and the oscillation of the bubble speeds is dependent on Re number.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.1-10
• /
• 1998

서로 섞이지 않는 두 비압축성 유체의 유동을 해석하기 위하여 VOF 방법에 기초한 수치 기법을 개발하였다. 유체간의 계면형상의 거동은 유동장내의 유체의 점유체적비의 변화에 의해 묘사되는데 이를 지배하는 이동방정식을 풀기 위한 새로운 대류항 계산법을 고안하였다. 대류항은 유체계면의 방향에 따라 풍상법과 역풍상법의 적절한 조합을 취하여 계산하는데 여기에 대각방향의 상류효과를 포함시켜 시간에 대한 2차 정확도를 갖도록 하였다. 또한 이 방법을 유량보정수송(FCT)법과 결합시켜 해의 단조성을 보장하였다. 몇 가지 단순 문제에 대한 시험 결과 이 기법이 수치오차에 의한 계면형상의 변형과 파손을 감소시킴을 확인하였다.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.10a
• /
• pp.269-272
• /
• 2005

For the three decades, the mold-filling of injection molding process was modeled as Hele-Shaw model. However, this model can not consider the 3D effect. In this paper, numerical simulations of three dimensional mold-filling during the filling phase were performed. The governing equations were discretized by segregated finite element method, which used equal order interpolation for pressure and velocity fields. The iterative linear equation solver (JCG, SOR) was employed for the solution of the momentum and pressure equations. Volume of Fluid (VOF) was employed for the melt front advancement. To check the validity of the numerical results, the results were compared with the experimental ones. The agreements between the experiment and the numerical results were found to be satisfactory.

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