• 한국전산유체공학회지
• /
• 제12권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.

• 한국화재소방학회논문지
• /
• 제30권6호
• /
• pp.111-117
• /
• 2016

본 연구는 스프링클러 헤드 근처에서 형성되는 액막의 자유표면 유동에 대해 CFD 모델을 적용하여 해석하고 스프링클러의 초기분무 특성 예측을 위한 기존 이론식의 결과와 비교를 통해 이론 모델의 타당성을 검토하였다. CFD 해석은 상용 해석프로그램인 CFX 14.0을 이용하였으며 노즐과 디플렉터로 이루어진 단순형상에 대해 표준난류모델과 VOF법을 적용하여 해석을 수행하였다. 평판부의 디플렉터 끝단에서 속도분포는 CFD 해석과 경험식이 매우 잘 일치된 결과를 보였으나 기하학적 형상이 복잡한 부분에서는 속도분포의 차이를 보였다. 이론모델에서 예측된 평균액적크기는 실제 스프링클러 헤드에서 측정된 평균액적크기에 대한 이전 연구결과와 큰 차이를 보였다. 그러나 이론 모델은 스프링클러 헤드의 초기 액적형성과정의 메커니즘을 이해하고 실험적 접근이 용이하지 않은 상황에서 분무 액적의 특성을 예측하는데 유동한 도구로 활용될 수 있다.

• 해양환경안전학회지
• /
• 제22권6호
• /
• pp.750-757
• /
• 2016

본 연구에서는 파 주기와 기울기에 따른 단일 실린더와 다중 실린더의 wave run-up을 추정하였다. 3차원 비압축성 점성 난류 유동이며 서로 다른 상을 가지는 이상유체에 대한 수치해석을 수행하기 위하여, 본 연구에서는 전산 유체 역학 상용코드인 "STAR-CCM+"을 기반으로 하여 VOF기법과 realizable $k-{\varepsilon}$ 난류 모델을 사용하였다. 모델스케일에 대한 파 주기는 단일 실린더의 경우 1.269초와 1.692초이고, 다중 실린더의 경우 1.716초이다. 각 케이스 별로 1/30, 1/16의 파 기울기를 가진다. 최종적으로 파 기울기와 주기에 대한 wave run-up 추정 결과는 관련 실험 결과와 비교하였다. 수치해석 결과는 실험과 비교하여 정성적으로 유사함을 확인하였다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제6권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
• /
• 제18권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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2017년도 제48회 춘계학술대회논문집
• /
• 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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제3권1호
• /
• pp.1-10
• /
• 1998

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

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2005년도 추계학술대회 논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제7권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.