• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.63-66
• /
• 2002

The present study introduces some useful methods of implementing the Internet numerical analysis program with existing numerical codes for utilizing the Internet environment. The Internet gives developers good environment for development and release. Several methods were suggested, and some of them were implemented with an existing numerical code named SOLA-VOF, a computational fluid dynamics program to solve two-dimensional transient flow problems with free surface. This was reconstructed with Java technologies and compared with the original one. Java technologies have been applied to development of Internet applications for a long time. The objective of this work is to contrive methods of implementing Internet numerical analysis program with existing numerical codes and confirm the possibility of them. Methods using the applet-servlet communication were suggested and implemented. In addition, the Java web services with XML was introduced, which makes possible the cooperation of components. Although the concept has been suggested and developed for business applications, it can also be used for engineering softwares. Therefore, this study will be a preparation for numerical analysis to participate in engineering web services.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.10
• /
• pp.979-988
• /
• 2012

The 3D structure of GDL for fuel cells was measured using high-resolution X-ray tomography in order to study material transport in the GDL. A computational algorithm has been developed to remove noise in the 3D image and construct 3D elements representing carbon fibers of GDL, which were used for both structural and fluid analyses. Changes in the pore structure of GDL under various compression levels were calculated, and the corresponding volume meshes were generated to evaluate the anisotropic permeability of gas within GDL as a function of compression. Furthermore, the transfer of liquid water and reactant gases was simulated by using the volume of fluid (VOF) and pore-network model (PNM) techniques. In addition, the simulation results of liquid water transport in GDL were validated by analogous experiments to visualize the diffusion of fluid in porous media. Through this research, a procedure for simulating the material transport in deformed GDL has been developed; this will help in optimizing the clamping force of fuel cell stacks as well as in determining the design parameters of GDL, such as thickness and porosity.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.604-604
• /
• 2012

국지성 집중호우에 따른 도심지 내수 침수 피해의 주원인으로 하수관거의 설계기준을 초과하는 강우가 침수피해의 주요 원인이며, 도심화로 인해 불투수 면적이 증가함에 따라 유출되는 시간이 짧아 저지대의 피해는 불가피하다. 2010년과 2011년에 100년 이상의 강우사상이 서울시에 연이어 나타나면서 집중호우로 인한 피해지역이 유사하게 나타났으며, 광화문 거리의 연이은 침수는 현재 서울시의 하수관거의 용량과 빗물펌프장 및 저류조 시설로 구성된 기존 수방대책의 한계점을 보이고 있다. 이에 서울시는 광화문 일대의 배수능력을 향상시키기 위하여 효자배수분구 빗물배수터널을 계획하고 있다. 일본, 미국 및 유럽 등지에서는 대심도 지하수로 시설에 대한 수리실험 및 수치 연구를 바탕으로 다양한 지하방수로가 건설되어 국지성 집중 강우에 대해 적절히 대응하고 있으나, 국내의 경우에는 대심도 지하방수로 시설에 대한 연구가 미비하여 지하방수로 설계 지침 및 기술적 자료가 부족한 실정이다. 그러므로 대심도 빗물배수터널 시설에서의 흐름특성 분석에 관한 수리실험 및 수치해석 등의 구체적인 연구가 필요하다고 판단된다. 본 연구에서는 수리모형 실험의 물질적 및 시간적 한계를 극복하기 위하여 일반적으로 3차원 유체거동의 특성분석에 많이 사용되는 Fluent 6.3 모형을 이용하여 대심도 빗물배수터널 시설의 접선식 유입구에 대한 흐름특성을 수치모의 하였다. 접선식 유입구 및 수직갱(drop shaft)에 대한 기하 모형의 격자망은 수치해석의 안정성 확보를 위하여 그림 1과 같이 6면체 격자로 구성하였다. 맨홀 내의 다상유동을 고려하기 위하여 VOF(Volume of Fluid) Scheme을 적용하였으며, 수치해석 방법으로는 비정상류, 1st order implicit method를 사용하였다. Fluent에서의 난류 흐름을 계산하는 방법에는 난류 운동에너지와 난류 에너지 소산율 $\epsilon$의 전달 방정식을 도입한 k-$\epsilon$ 난류 모형을 채택하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2018.05a
• /
• pp.251-251
• /
• 2018

도로설계기준(2012)에 따르면 노면 배수시설은 측구, 집수정, 배수관, 배수구(빗물받이, 맨홀) 등으로 구성되며 중앙분리대 배수에서 구조가 방호벽일 경우 원칙적으로 집수정과 종배수관, 횡배수관을 설치하여 노면수를 배수한다. 최근 기후변화에 따른 국지성 집중호우로 교량 배수시설의 배수능력 부족으로 교량의 노후화, 노면의 체수와 수막현상을 야기한다. 교량 설계 시, 교량의 외관은 점배수에 의한 배수불량으로 막힘부에서의 식물생육, 부식, 배수구의 오염 등으로 심미성에 악영향을 미치고 있다. 따라서 국토교통부 도로 배수시설 설계 및 관리지침에서는 기존의 집수정으로 유도하여 배수하는 점배수 형태에서 우수의 지체시간을 저감하기 위해 배수구로 즉시 배수되어 배수 효율을 증가시켜주는 선배수시설을 권장하고 있다. 그러나 선배수 구조의 일반적인 내용만을 기술하고 있으므로 보다 상세한 설계기준이 필요하다. 또한 국내에서는 교면 배수를 위한 횡배수구의 배수능력에 관한 설계 개념은 미국의 빗물받이 설계 개념을 그대로 적용하고 있어서 교량 배수시설의 배수능력 증대를 위한 수치 및 수리실험에 관한 기술적인 자료가 부족한 실정이다. 따라서 도로의 흐름해석 및 선배수를 위한 측구 횡배수관의 흐름개선에 관한 구체적인 연구가 필요한 실정이다. 본 연구에서는 배수로 전반의 유향 분석 및 횡유입관 유입부에서의 상세한 흐름해석과 통수능력을 분석하기 위해서는 수리실험을 통한 연구가 필요한 실정이다. 그러나 최적 횡배수관의 형상 및 간격 설정에 대한 수리실험의 물리적 및 시간적 한계를 극복하기 위해 Fluent 모형(Ansys Workbench 13.0)을 활용한 수치모의를 수행하여 모형의 적용성을 검토하였다. 다상유동 해석을 위해 VOF(Volume of Fluid)방법을 적용하였고, 수치해석 방법으로는 비정상류, 난류 모형으로는 standard ${\kappa}-{\varepsilon}$모형을 적용했다. 도로 형상에 따른 우수유출량을 비교 분석하기 위하여 횡경사는 2%로 고정하고 종경사를 2%로 선정하여 수치모의를 통한 배수능력을 분석하고 설계에 직접적인 적용성을 검토하였다. 횡배수관의 간격변화 및 배수공의 위치 변화 등에 따른 차집량을 분석하였으며, 개략적인 수리실험 결과와 수치모의의 차집율을 비교 및 분석하여 Fluent 모형의 적용성을 확인하였다. 또한, 횡배수관의 유입부에서의 유속 변화 및 유출부에서의 유속이 모의가 가능하므로 배수시설에서 보다 정확한 흐름 해석이 가능하여 보다 적정한 배수능력 분석이 가능할 것으로 판단된다.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.4
• /
• pp.313-321
• /
• 2019

The prediction of the hydrodynamic performance of a planing hull vessel is an important and challenging topic for computational fluid dynamic (CFD) applications to naval hydrodynamics. In this paper, the resistance and planing attitude analysis for a Fridsma hull, which is a prismatic planing hull, in still water are numerically studied using OpenFOAM. OpenFOAM is an open source code package based on C++ libraries and the finite volume method (FVM) for the discretization of the RANS equation. The volume of fluid method (VOF) is used to capture the water-air interface and the SST ${\kappa}-{\omega}$ model is used for the turbulence simulation. The overset mesh method is used to capture the large motion of the hull at higher speeds. Before the extensive analysis, uncertainty analyses using various time steps and grid sizes were performed for one ship speed case of Fn = 1.19. The results of the present study are compared with those of a model test, other CFD research, and Savitsky's empirical formula. The results of the present study, following the trend of other CFD results, slightly over predict the resistance and under predict the sinkage and, more significantly, the trim.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.4 no.4
• /
• pp.21-31
• /
• 2001

When an oil-spilling accident occurs at sea, it is of the primary importance to predict the amount of oil leakage for the swift response and decision-making. The simplest method of oil-leakage estimation is based on the hydrostatic pressure balance between oil inside the tank and seawater outside of leakage hole, that is the so-called Torricelli equilibrium relation. However, there exists discrepancy between the reality and the Torricelli relation, since the latter is obtained from the quasi-steady treatment of Bernoulli equation ignoring viscous friction. A preliminary experiment has been performed to find out the oil-leaking speed and shape. Soy-bean oil inside the inner tank was ejected into water of the outer tank through four different leakage holes to record the amount of oil leakage. Furthermore, a CFD (Computational Fluid Dynamics) method was utilized to simulate the experimental situation. The Wavier-Stokes equations were solved for two-density flow of oil and water. VOF method was employed to capture the shape of their interface. It is found that the oil-leaking speed varies due to the frictional resistance of the leakage hole passage dependent on its aspect ratio. The Torricelli factor relating the speed predicted by using the hydrostatic balance and the real leakage speed is assessed. For the present experimental setup, Torricelli factors were in the range of 35%~55% depending on the aspect ratio of leakage holes. On the other hand, CFD results predicted that Torricelli factor could be 52% regardless of the aspect ratio of the leakage holes, when the frictional resistance of leakage hole passage was neglected.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.245-245
• /
• 2022

최근 국내 기후변화에 따른 국지성 집중호우로 인한 시간당 강우량의 증가로 도로부 유출량의 증가와 배수관거에서의 내수배제 불량에 따른 도심지 내수침수 피해가 증가함에 따라 이를 해결하기 위한 우수유출저감시설이 설치되고 있다. 그러나 대단위의 지하 저류시설의 지속적인 설치는 과밀화된 도심지에서 설치 지하공간의 구조적인 한계 및 적정 설치 위치의 미확보 등의 다양한 문제가 발생하여 저류시설의 침수저감 효과에 대한 추가적이고 새로운 저류시설에 대한 연구가 필요한 실정이다. 이에 내수 침수 저감 및 배수 능력 향상을 위한 도로 배수시설과 연계된 도로 측구부 저류시스템 구축이 필요하다. 이를 위해 역류 방지 및 노면수 저류 빗물받이에 적용되는 부력식 역류차단장치를 개발하였으며, 역류차단장치의 최적 형상 개발을 위해서 기존 빗물받이 연결관과의 통수능 비교 및 분석이 필요한 실정이다. 따라서 본 연구에서는 기존 빗물받이 연결관 및 연결관 내에 역류차단장치가 적용된 역류차단 빗물받이의 흐름분석을 위해 Fluent 모형을 이용하여 3차원 수치모의를 수행하였다. 수치모의 구성으로는 전체 형상을 40×50cm의 빗물받이 유입부와 50×50cm의 빗물받이로 결정하고 격자는 빗물받이 내부의 복잡한 3차원 흐름을 모의하기 위해 1.2~2mm 크기로 생성하였다. 다상유동 해석을 위해 VOF(Volume of Fluid)방법을 적용하였고, 수치해석 방법으로는 비정상류, 난류 모형으로는 SST k-𝜔모형을 적용하였다. 해석조건으로는 김정수(2021) 등이 제시한 4차선 기준 설계빈도별(5~30년) 빗물받이 유입유량을 산정하여 빗물받이 유입조건으로 선정하였으며, 빗물받이와 연결관에서의 통수능력 분석 조건으로는 빗물받이에 기존 연결관이 부착된 조건과 연결관 내에 역류차단장치가 설치되어 역류차단장치가 개방된 조건에서의 통수능을 비교하였으며, 역류상황을 가정한 연결관에서의 통수능을 비교하기 위하여 역류차단장치의 개폐정도를 15도(통수단면 33%감소) 닫힌 상태 및 30도(통수단면 67% 감소) 닫힌 상태 조건을 대상으로 빗물받이와 연결관에서의 흐름을 모의하였다. 수치모의 결과 역류차단장치의 계폐조건에 상관없이 5년 빈도유입량 조건에서는 완전 배수가 되었으며, 개폐조건 15도에서는 10년 빈도의 유입량에서는 완전 배수가 되었으나 20년 빈도 이상의 유입량 조건에서 빗물받이 유입부로의 역류가 발생하였으며, 개폐조건 30도에서는 5년 빈도 이상 유입량 조건에서 빗물받이 유입부로 역류가 발생하는 것으로 나타났다. 특히, 30년 빈도 이상의 유입량부터는 빗물받이 연결관 내에 역류차단장치 개페조건과 관계없이 빗물받이 유입부로의 역류로 인한 도로 침수가 발생하기 때문에 유휴공간인 도로 측구부를 저류공간으로 활용할 수 있는 도로 측구부 저류시스템의 구축은 필수적이라고 판단되며, 유량 조건에 따른 빗물받이 내부 와 흐름과 유출부에서의 유속 변화 특성을 확인하였다. 그러므로 측구 저류조 개발 형상과 연결한 3차원 흐름의 구현 및 분석에 Fluent 모형의 적용이 가능하다고 판단된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.234-234
• /
• 2022

도로에서의 우수를 원활하게 처리하기 위해서 빗물받이 및 연결관 등의 노면 배수시설이 설치되고 있으며, 노면 배수는 측구부를 통해 흘러 빗물받이 유입부로 차집되고 연결관을 통해 하수관거로 배수된다. 그러나 최근 국내 기상패턴의 변화로 국지성 집중호우와 같이 시간당 강우량 증가로 도로부와 저지대에서 배수시설의 배수불량에 따른 도심지 내수침수 피해가 발생하고 있다. 이에 정부에서는 다양한 우수관거 개선사업, 빗물펌프장, 지하저류조와 같은 방재시설을 설치하고 있으나 우수유출저감시설은 대규모 예산이 소요되고 실제 침수지역에 피해 저감효과에 대한 효용성 문제에 대한 제기뿐만 아니라 과밀화된 도심지에서는 지하공간 활용에 한계가 있는 실정이므로 도심지의 다양한 공간을 활용한 도시 배수 및 저류시설에 대한 연구가 필요하다. 따라서 본 연구에서는 유휴 공간인 도로 측구부 공간을 활용하여 도로 노면수를 저류 및 지체할 수 있는 노면수 측구 저류시설의 개념을 제시하고 측구저류조의 활용성을 판단하기 위하여 빗물받이 유입구, 빗물받이, 측구 저류조 및 빗물받이와 측구저류조 연결부에서의 노면수 유입, 유출 및 저류 등의 다양한 흐름 변화를 확인하기 위하여 Fluent 모형의 적용성을 분석하였다. 수치모의 전체 형상은 50x50cm 크기의 빗물받이를 기준으로 양쪽에 2m 길이의 측구 저류조를 원형관으로 연결하여 1/5 축소모형으로 구성하고 격자는 빗물받이 유입부, 빗물받이 및 측구 저류조 내부의 복잡한 3차원 흐름을 모의하기 위해 사면체와 육면체로 조밀하게 생성하였다. 다상유동해석을 위해 VOF(Volume of Fluid)방법을 적용하였고, 수치해석 방법으로는 비정상류, 난류 모형으로는 SST k-ω모형을 적용하였다. 수치모의 조건으로는 설계빈도별(5~30년) 우수유출량을 산정하여 유입 유량별 기존 빗물받이 유입부에서의 유입흐름, 빗물받이 내부에서의 와 발생흐름, 측구 저류조 및 연결관에서의 흐름을 구현하여 분석하였다. 수치모의 결과 빗물받이 유입부에서 연결관을 통한 측구 저류조로 유입되는 유입흐름과 빗물받이 하단부의 배수관을 통해 유출되는 흐름을 정상적으로 구현하였으며, 빗물받이 유입부 및 측구 저류조 연결관에서의 유속변화도 확인할 수 있었다. 또한 빗물받이와 측구 저류조에서 다양한 흐름을 구현하기 위한 Flunet 모형의 적용성을 검토하였으며, 향후 수리실험을 통하여 실제 흐름과의 매개변수 최적화 및 다양한 도로 조건의 변화를 고려한 수치모의 분석을 통하여 지속적인 모형의 검증이 가능할 것으로 판단된다.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.17 no.3
• /
• pp.182-188
• /
• 2014

In an oscillating water column (OWC)-type wave energy conversion system, the performance of the OWC chamber depends on the chamber shape, as well as the incident wave direction and pressure drop produced by the turbine. Although the previous studies on OWC chambers have focused on wave absorbing performance in ideal operating conditions, incident waves do not always arrive normally to the OWC chamber in real sea conditions, especially in fixed devices. The present study deals with experiments and numerical calculations to investigate the effects of wave direction on the performance of the OWC chamber. The experiments were carried out in a three-dimensional wave basin for five different wave directions, including the effect of turbine using the corresponding orifice. The wave elevation inside the chamber was measured at the center point under various incident wave conditions. The numerical study was conducted by using a numerical wave tank-based volume-of-fluid model to compare the results with experimental data and to reveal the detailed flows around the chamber.