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

• 한국해안·해양공학회논문집
• /
• 제32권1호
• /
• pp.1-10
• /
• 2020

자유수면을 포함하는 파동장과 같이 단상의 경계가 시간발전에 따라 지속적으로 변화하는 경우나 액상과 기상이 혼합되는 문제에 있어서는 다상유동(multiphase flow) 문제를 적용하는 예가 증가하고 있다. 특히, 파동장과 같은 자유수면의 문제를 취급하는데 있어서는 혼합되지 않는 액상과 기상의 비압축성 뉴턴유체를 고려한 혼상류 모델이 적용되는 경우가 많다. 일반적으로 혼상류 모델은 각상의 경계면에 대한 시간기반 거동추적이 필수적이며, 궁극적으로는 계산의 정도를 좌우한다. 본 연구는 다양한 CFD 수치해석코드에 적용되고 있는 대표적인 VOF-type의 경계면 추적기법들의 이류성능을 평가하였다. 특히, 기존의 전통적인 VOF-type의 경계면 추적기법 및 이류계산에서 발생하는 수치확산을 최소화하기 위해 수치유속(numerical flux)을 제어하는 FCT 법의 효용성을 평가하고, 더불어 CIP 법을 활용한 자유수면 추적성능의 가능성을 고찰하였다. 그 결과, 본 연구에서 적용한 제한된 조건하에서는 수치확산 방지를 위해 수치확산방지 유속을 도입한 FCT-VOF 법이 가장 높은 경계면의 추적성능을 보였다. 본 연구에서 도출되는 결과는 다양한 수치해석코드에 적용되는 자유수면의 추적기법을 선택함에 있어서 중요한 기초자료로 활용될 것으로 기대된다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 춘계학술대회논문집
• /
• pp.215-222
• /
• 2009

In the present study, a numerical analysis on the sloshing in a tank with the harmonic motion was investigated. A VOF method was used for two-phase flows inside the sloshing tank and a source term of the momentum equation was applied for the harmonic motion. This numerical method was verified by comparing its results with the available experimental data. The sloshing in a tank causes the instability of the fluid flows and the fluctuation of the impact pressure on the tank. By these phenomena of the tank sloshing, the sloshing problems such as the failure and the noise of system can be generated. For the reduction of these sloshing problems, the various baffles such as the horizontal/vertical plate baffles and the porous baffles inside the tank are installed. With the installations of these baffles, the characteristics of the liquid behavior in the sloshing tank, the impact pressure on the wall, the amplitude of the free surface near the wall and the sloshing noise were numerically analyzed.

• 한국전산유체공학회지
• /
• 제14권2호
• /
• pp.68-76
• /
• 2009

The sloshing tank causes the instability of the fluid flows and the fluctuation of the impact pressure by the liquid on the tank. These flow characteristics inside the sloshing tank can generate the uncomfortable sloshing noise. In the present study, a numerical analysis for the reduction of a fuel tank sloshing noise was performed. To simulate the flow characteristics in a sloshing tank with partially filled liquid, a VOF method was used for interfacial flows by applying a momentum source term for the sloshing motion in a non-inertial reference frame. This numerical method was verified by comparing its results with the available experimental data. For the reduction of the sloshing noise, the horizontal and vertical baffles and porous media inside a sloshing tank were considered and numerically analyzed in the present study. For various installations of these baffles and porous media, the characteristics of the liquid behavior in the sloshing tank were obtained along with the impact pressure on the wall and the height of the free surface along the wall. These basic results can be used for the design of the actual vehicular fuel tank with the reduced sloshing noise.

• Wind and Structures
• /
• 제14권5호
• /
• pp.435-447
• /
• 2011

A Computational Fluid Dynamics model is presented in this study for the simulation of the complex fluid flows with free surfaces inside the Tuned Liquid Column Dampers in horizontal motion. The characteristics of the fluid model of the TLCD in horizontal motion include the free surface of the multiphase flow and the horizontal moving frame. In this study, the time depend unsteady Standard ${\kappa}-{\varepsilon}$ turbulent model based on Navier-Stokes equations is chosen. The volume of fluid (VOF) method and sliding mesh technique are adopted to track the free surface of water inside the vertical columns of TLCD and treat the moving boundary of the walls of TLCD in horizontal motion. Several model solution parameters comprising different time steps, mesh sizes, convergence criteria and discretization schemes are examined to establish model parametric independency results. The simulation results are compared with the experimental data in the dimensionless amplitude of the water column in four different configured groups of TLCDs with four different orifice areas. The predicted natural frequencies and the head loss coefficient of TLCDs from CFD model are also compared with the experimental data. The predicted numerical results agree well with the available experimental data.

• Journal of Mechanical Science and Technology
• /
• 제20권6호
• /
• pp.840-848
• /
• 2006

We carried out numerical studies to investigate the single- and two-phase flow characteristics in the single- and multi-channels. We used the finite volume method to solve the mass and momentum conservation equations. The volume of fluid model is used to predict the two-phase flow in the channel. We obtained the distribution of velocity fields, pressure drop and air volume fraction for different water mass flow rates. We also calculated the distribution of mass flow rates in the multi-channels to understand how the flow is distributed in the channels. The calculated results for the single- and two-phase flow are partly compared with the present experimental data both qualitatively and quantitatively, showing relatively good agreement between them. The numerical scheme used in this study predicts well the characteristics of single-and two-phase flow in a multi-channel.

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

• Nuclear Engineering and Technology
• /
• 제53권4호
• /
• pp.1119-1126
• /
• 2021

Spacer grid with mixing vane had been widely used in nuclear reactor core. One of the main feather of spacer grid with mixing vane was that strong swirl flow was formed after the spacer grid. The swirl flow not only changed the bubble generation in the near wall field, but also affected the bubble behaviors in the center region of the subchannel. The interaction between bubble and the swirl flow was one of the basic phenomena for the two phase flow modeling in fuel assembly. To obatin better understanding on the bubble behaviors in swirl flow, full three dimension numerical simulations were conducted in the present paper. The swirl flow was assumed in the cylindral calculation domain. The bubble interface was captured by Volume Of Fluid (VOF) method. The properties of saturated water and steam at different pressure were applied in the simulation. The bubble trajectory, motion, shape and force were obtained based on the bubble parameters captured by VOF. The simulation cases in the present study included single bubble with different size, at different angular velocity conditions and at different pressure conditions. The results indicated that bubble migrated to the center in swirl flow with spiral motion type. The lateral migration was mainly related to shear stress magnitude and bubble size. The bubble moved toward the center with high velocity when the swirl magnitude was high. The largest bubble had the highest lateral migration velocity in the present study range. The effect of pressure was small when bubble size was the same. The prelimenery simulation result would be beneficial for better understanding complex two phase flow phenomena in fuel assembly with spacer grid.

• 한국해안·해양공학회논문집
• /
• 제20권1호
• /
• pp.1-13
• /
• 2008

본 연구에서는 복수의 대형 연직 주상구조물을 중심으로 구조물에 작용하는 비선형파력과 구조물에 의한 비선형파랑변형을 검토하며, 수치해석의 이론으로 쇄파나 강비선형 파랑을 고정도로 해석할 수 있는 Navier-Stokes방정식에 근거한 3차원 VOF(Volume Of Fluid)법을 적용하였다. 본 수치해석의 타당성을 검증하기 위해 파랑변형과 파력에 대한 기존의 수리실험결과와 본 수치해석결과를 비교 검토하였으며, 이로부터 2기로 구성된 임의형상의 연직 주상구조물에 대해 구조물 형상, 구조물 배치간격, 파의 입사각도 변화에 따른 파력과 파랑변형의 특성을 논하였다.

• 한국항해항만학회지
• /
• 제34권6호
• /
• pp.417-422
• /
• 2010

정수중 부선 주위의 유동 특성을 조사 및 분석하기 위하여 부선의 저항성능을 모형시험과 CFD를 이용한 수치계산에 의하여 해석하였다. 모형실험은 무한수심 조건에서의 회류수조에서 저항과 선측파형 등을 관측하였으며, 유동장 해석을 위하여 VOF(Volume of Fluid)법을 기초로 한 수치 시뮬레이션을 이행하였다. 부선의 형태에 따르는 유동특성을 파악하기 위하여 서로 다른 모습의 부선 모형을 제작하였으며, 모형실험은 설계속도를 7kts로 정하고 역조와 순조의 조류 영향을 고려하여 5kts~10kts의 구간에서 실험을 실시하여 각 속도에 따르는 저항 특성을 비교 검토 하였다. 수치 시뮬레이션은 일정 속도에서 만재상태를 기준으로 하여 계산을 수행하였으며 수조 실험 결과와 상호 비교함으로서 본 수치계산법의 타당성을 검증할 수 있었으며, 부선의 저항성능 및 선체 주위 유동장의 특성을 분석하였다. 저속일 때 25% 정도의 저항의 차이가 발생하고, 7kts부터 유효마력이 급격히 증가하는 현상을 파악할 수 있었다.