• 신재생에너지
• /
• 제3권4호
• /
• pp.8-15
• /
• 2007

Liquid water in flow channel is an important factor that limits the steady and transient performance of PEM fuel cells. A computational fluid dynamics study based on the volume-of-fluid [VOF] multi-phase model was conducted to understand the two-phase flow behavior of liquid water in cathode gas channels. The liquid water transport in $180^{\circ}{\Delta}$ bends was investigated, where the effects of surface characteristics (hydrophilic and hydrophobic surfaces], channel geometries (rectangular and chamfered corners], and air velocity in channel were discussed. The two-phase flow behavior of liquid water with hydrophilic channel surface and that with hydrophobic surface was found very different; liquid water preferentially flows along the corners of flow channel in hydrophilic channels while it flows in rather spherical shape in hydrophobic channels. The results showed that liquid water transport was generally enhanced when hydrophobic channel with rounded corners was used. However, the surface characteristics and channel geometries became less important when air velocity was increased over 10m/s. This study is believed to provide a useful guideline for design optimization of flow patterns or channel configurations of PEM fuel cells.

• 해양환경안전학회지
• /
• 제20권3호
• /
• pp.312-317
• /
• 2014

In this paper, the wave run-up height and depression depth around air-water interface-piercing circular cylinder have been numerically studied. The Reynolds Averaged Navier-Stokes equations (RANS) and continuity equations are solved with Reynolds Stress model (RSM) and volume of fluid (VOF) method as turbulence model and free surface modeling, respectively. A commercial Computational Fluid Dynamics (CFD) software "Star-CCM+" has been used for the current simulations. Various Froude numbers ranged from 0.2 to 1.6 are used to investigate the change of air-water interface structures around the cylinder and experimental data and theoretical values by Bernoulli are compared. The present results showed a good agreement with other studies. Kelvin waves behind the cylinder were generated and its wave lengths are longer as Froude numbers increase and they have good agreement with theoretical values. And its angles are smaller with the increase of Froude numbers.

• 한국전산유체공학회지
• /
• 제17권2호
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• 대한조선학회논문집
• /
• 제47권1호
• /
• pp.38-46
• /
• 2010

This study aims at investigating the effect of the baffle height on the liquid sloshing in the three-dimensional (3D) rectangular tank. In order to simulate the 3D incompressible viscous two-phase flow in the 3D tank with partially filled liquid, the present study has adopted the volume of fluid (VOF) method based on the finite-volume method which has been well verified by comparing with the results of the relevant previous researches. The ratio of the baffle height ($h_B$) to filling level (h) has been changed in the range of $0{\leq}h_B/h{\leq}1.2$ to observe the effect on the impact loads on the side wall and free surface behavior. Generally, as baffle height increases, the impact pressure on the wall decreases and the deformation of free surface becomes weaker. However it seemed that a critical ratio of the baffle height existed to reveal the lowest impact pressure on the wall. Consequently, $h_B/h=0.8$ among $h_B/hs$ considered in the study showed the lowest impact pressure.

• 한국수소및신에너지학회논문집
• /
• 제30권6호
• /
• pp.585-592
• /
• 2019

In this study, a flashing boiling phenomenon of pressurized water jet was numerically studied and validated against an experimental data in the literatures. The volume of fluid (VOF) technique was used to consider two-phase behavior of water, while the homogeneous relaxation model (HRM) model was used to provide the velocity of phase change. During the flashing boiling through a nozzle, a mach disk was observed near nozzle exit because of pressure drop resulting from two-phase under-expansion. The flashing jet structure, local distributions of temperature/vapor volume fraction/velocity, and position of the mach disk were examined as nozzle inlet temperature changed.

• 대한기계학회논문집B
• /
• 제30권10호
• /
• pp.1003-1011
• /
• 2006

Droplets are ejected onto a substrate through a nozzle by pushing liquids in flow channels of drop-on-demand devices. The behavior of ejection and formation of droplets is investigated to enhance the physical understanding of the hydrodynamics involved in inkjet printing. The free surface phenomenon of a droplet is described using $CFD-ACE^{TM}$ which employs the volume-of-fluid (VOF) method with the piecewise linear interface construction (PLIC). Droplet formation characteristics are analyzed in various flow regimes with different Ohnesorge numbers. The computational results show that the droplet formations are strongly dependent on the physical properties of working fluids and the inlet flow conditions. In addition, the wetting characteristics of working fluids on a nozzle influence the volume and velocity of a droplet produced in the device. This study may provide an insight into how a liquid droplet is formed and ejected in a piezoelectric inkjet printing device.

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

본 연구에서는 3차원수치파동수로내에 쇄파후의 파랑이 대형연직원주구조물에 작용할 때 작용파력과 구조물에 의한 파랑변형을 수치적으로 해석한다. 수치해석법으로 파랑과 구조물과의 비선형상호간섭에 따른 쇄파현상을 포함하는 복잡한 자유수면의 거동특성을 고정도로 해석할 수 있는 3차원Navier-Stokes운동방정식과 자유수면추적에 3차원VOF(Volume Of Fluid)법을 결합한 강비선형해석법을 적용한다. 3차원파동장내에서 해저는 쇄파를 상대적으로 쉽게 발생시킬 수 있는 경사스텝의 해저(변수심의 경사수역과 일정수심역으로 구성)로 이루어진 경우를 고려하며, 파고의 변화에 따라 쇄파가 경사수역 또는 일정수심역에서 발생하여 일정수심역의 대형연직원주구조물에는 쇄파후의 파랑만이 작용하는 경우로 한정한다. 구조물의 위치 및 입사파랑의 파고변화가 구조물에 작용하는 파력 및 파랑변형에 미치는 특성을 쇄파전후의 파랑을 중심으로 검토하고, 쇄파후 파랑의 전파에 따른 파랑에너지의 변화와 구조물에 작용하는 파력특성과를 연관시켜 논의하여 3차원파동장에서 파랑과 구조물과의 강비선형간섭현상의 특성을 규명한다.

• 한국음향학회지
• /
• 제42권3호
• /
• pp.250-261
• /
• 2023

본 연구에서는 펌프젯 추진기를 대상으로 공동, 비공동 조건에서의 유동 소음원을 규명하기 위하여 추진기의 각 구성품인 덕트와 스테이터, 로터에 의한 소음 기여도를 평가하였으며, 공동과 비공동 조건에서의 소음 수준을 비교하였다. 대형 캐비테이션 터널 내 Suboff 잠수함 선형과 펌프젯 추진기를 대상으로 균일혼상류 가정의 비정상 비압축성 Reynolds averaged Navier-Stokes(RANS) 방정식을 적용하였으며, 이상 유동을 모사하기 위해 Volume of Fluid(VOF) 기법과 Schnerr-Sauer 공동 모델을 적용하였다. 유동해석 결과를 기반으로 수중방사소음을 예측하기 위해 Ffowcs Williams and Hawkings(FW-H) 방정식 기반의 음향상사법을 적용하였으며, 덕트와 스테이터, 로터로 구성된 3개의 비투과성 적분면과 추진기를 감싸는 형태의 2가지 투과성 적분면을 선정하여 소음 기여도를 평가하였다. 소음 예측결과로부터 스테이터는 전체 소음에 대한 직접적인 기여도는 낮으나 덕트와 로터에서의 유동 박리에 의한 소음원 형성에는 영향을 미치는 것을 확인하였으며, 유동이 박리되는 연직상방과 우측방향으로 소음이 크게 방사되었다. 또한 로터에서는 날개의 흡입면과 압력면 간의 압력 섭동에 의해 추진방향으로 소음이 크게 방사되었으며, 투과성적분면을 통해 체적 소음원인 공동의 효과를 반영할 수 있음을 확인하였다.

• 한국수자원학회논문집
• /
• 제43권9호
• /
• pp.801-811
• /
• 2010

급경사에서의 고립파의 처오름을 예측하기 위해 3차원 수치모형에 내부조파기법을 도입하여 수치모형실험을 수행하였다. 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 푸는 LES(large-eddy-simulation) 기반의 수치모형을 사용한다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two-step projection 기법을 사용하여 SANS 방정식을 풀었으며, Poisson 방정식을Bi-CGSTAB 기법을 이용하여 풀고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 고립파를 3차원 공간의 일정 수심상에서 내부조파하여 해석해와 비교한 후 분산오차에 대해 분석하였다. 그리고 고립파를 내부조파하여 급경사에서의 고립파의 처오름 및 처내림 현상을 예측하고 수리모형 실험결과와 비교 및 분석하였다.

• Wind and Structures
• /
• 제38권4호
• /
• pp.261-275
• /
• 2024

Before performing an experimental study on the downburst-generated wave, it is necessary to examine the scale effects and corresponding corrections or compensations. Analysis of similarity is conducted to conclude the non-dimensional force ratios that account for the dynamic similarity in the interaction of downburst with wave between the prototype and the scale model, along with the corresponding scale factors. The fractional volume of fluid (VOF) method in association with the impinging jet model is employed to explore the characteristics of the downburst-generated wave numerically, and the validity of the proposed scaling method is verified. The study shows that the location of the maximum radial wind velocity in a downburst-wave field is a little higher than that identified in a downburst over the land, which might be attributed to the presence of the wave which changes the roughness of the underlying surface of the downburst. The impinging airflow would generate a concavity in the free surface of the water around the stagnation point of the downburst, with a diameter of about two times the jet diameter (D_jet). The maximum wave height appears at the location of 1.5D_jet from the stagnation point. Reynolds number has an insignificant influence on the scale effects, in accordance with the numerical investigation of the 30 scale models with the Reynolds number varying from 3.85 × 10⁴ to 7.30 × 10⁹. The ratio of the inertial force of air to the gravitational force of water, which is denoted by G, is found to be the most significant factor that would affect the interaction of downburst with wave. For the correction or compensation of the scale effects, fitting curves for the measures of the downburst-wave field (e.g., wind profile, significant wave height), along with the corresponding equations, are presented as a function of the parameter G.