• Nuclear Engineering and Technology
• /
• v.56 no.6
• /
• pp.2029-2038
• /
• 2024

A deep understanding of the characteristics and mechanism of geyser boiling and capillary pumping is necessary to optimize a high-temperature sodium heat pipe. In this work, the Volume of Fluid (VOF) two-phase model and the capillary force model in the mesh wick were used to model the complex phase change and fluid flow in the heat pipe. Computational Fluid Dynamics (CFD) simulations successfully predicted the process of bubble nucleation, growth, aggregation, and detachment from the wall in the liquid pool of the evaporation section of the heat pipe in horizontal and tilted states, as well as the reflux phenomenon of capillary suction within the wick. The accuracy and stability of the capillary force model within the wick were verified. In addition, the causes of geyser boiling in heat pipes were analyzed by extracting the oscillation distribution of heat pipe wall temperature. The results show that adding the capillary force model within the wick structure can reasonably simulate the liquid backflow phenomenon at the condensation; Under the horizontal and inclined operating conditions of the heat pipe, the phenomenon of local dry-out will occur, resulting in a sharp increase in local temperature. The speed of bubble detachment and the timely reflux of liquid sodium (condensate) replenishment in the wick play a vital role in the geyser temperature oscillation of the tube wall. The numerical simulation method and the results of this study are anticipated to provide a good reference for the investigation of geyser boiling in high-temperature heat pipes.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.3
• /
• pp.129-145
• /
• 2019

For the design of composite breakwater as representative one of the coastal and harbor structures, it has been widely discussed by the researchers about the relation between the behavior of excess-pore-water pressure inside the rubble mound and seabed caused by the wave load and its structural failure. Recently, the researchers have tried to verify its relation through the numerical simulation technique. The above researches through numerical simulation have been mostly applied by the linear and nonlinear analytic methods, but there have been no researches through the numerical simulation by the strongly nonlinear mutiphase flow analytical method considering wave-breaking phenomena by VOF method and turbulence model by LES method yet. In the preceding research of this study, olaFlow model based on the mutiphase flow analytical method was applied to the nonlinear interaction analysis of regular wave-composite breakwater-seabed. Also, the same numerical techniques as preceding research are utilized for the analysis of irregular wave-composite breakwater-seabed in this study. Through this paper, it is investigated about the horizontal wave pressures, the time variations of excess-pore-water pressure and their frequency spectra, mean flow velocities, mean vorticities, mean turbulent kinetic energies and etc. around the caisson, rubble mound of the composite breakwater and seabed according to the changes of significant wave height and period. From these results, it was found that maximum nondimensional excess-pore water pressure, mean turbulent kinetic energy and mean vorticity come to be large equally on the horizontal plane in front of rubble mound, circulation of inflow around still water level and outflow around seabed is formed in front of rubble caisson.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.3
• /
• pp.575-583
• /
• 2017

Various types of flow patterns around the stepped weir and spillway, such as the skimming flow over such structures and the wave-type flow with a standing undular hydraulic jump and roller downstream of the structures, are developed in open channels. Unsteady three-dimensional numerical simulations are carried out using a hybrid RANS-LES turbulence modeling approach and the volume of fluid method for resolving free surface fluctuations to represent the turbulent flow including the skimming flow and wave-type flow over a stepped weir installed in a rectangular channel. The comparison of numerical results with an existing experimental measurement reveals that the present numerical simulations reasonably well reproduce the turbulent flow passing the stepped weir, in terms of time-averaged velocity profiles at selected locations downstream of the weir, flow topology characterized by the wave-type and skimming flows, the maximum height and length of the standing wave and the length of reattachment of recirculating zone. The numerical result further elucidates the distinct flow behaviors of the wave-type and skimming flow by presenting instantaneous intense variations of free surface and velocity vectors, the distributions of Reynolds shear stress and turbulent kinetic energy and three-dimensional complex features of coherent structures and total pressure distribution.

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

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

• Journal of the Korea Society of Computer and Information
• /
• v.16 no.11
• /
• pp.67-76
• /
• 2011

This paper describes the development of interactive visualization techniques and a program that allow us to visualize the structure of the volume data and interactively select and visualize the isosurface components using contour tree. The main characteristic of this technique is to provide an algorithm that draws the contour tree in 2D plane in a way that users easily understand the tree, and to provide an algorithm that can efficiently extract an isosurface component utilizing GPU's parallel architecture. The main characteristic of the program we developed through implementing the algorithms is to provide us with an interactive user interface based on the contour tree for extracting an isosurface component and visualization that integrates with previous isosurface and volume rendering techniques. To show the excelland vof our methods, we applied 3D biomedical volume data to our algorithms. The results show that we could interactively select the isosurface components that represent a polypeptide chain, a ventricle and a femur respectively using the user interface based on our contour tree layout method, and extract the isosurface components with 3x-4x higher speed compared to previous methods.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.7
• /
• pp.930-938
• /
• 2018

Underwater noise radiated from submarines is directly related to the probability of being detected by the sonar of an enemy vessel. Therefore, minimizing the noise of a submarine is essential for improving survival outcomes. For modern submarines, as the speed and size of a submarine increase and noise reduction technology is developed, interest in flow noise around the hull has been increasing. In this study, a noise analysis technique was developed to predict flow noise generated around a submarine shape considering the free surface effect. When a submarine is operated near a free surface, turbulence-induced noise due to the turbulence of the flow and bubble noise from breaking waves arise. First, to analyze the flow around a submarine, VOF-based incompressible two-phase flow analysis was performed to derive flow field data and the shape of the free surface around the submarine. Turbulence-induced noise was analyzed by applying permeable FW-H, which is an acoustic analogy technique. Bubble noise was derived through a noise model for breaking waves based on the turbulent kinetic energy distribution results obtained from the CFD results. The analysis method developed was verified by comparison with experimental results for a submarine model measured in a Large Cavitation Tunnel (LCT).

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5B
• /
• pp.473-482
• /
• 2009

This study numerically discusses wave forces acting on a vertical wall such as breakwaters or revetments, subjected to incident undular or turbulent bores. Due to the complex hydrodynamics of bore, its wave forces have been predicted, mainly through laboratory experiments. Numerical simulations in this paper were carried out by CADMAS-SURF(CDIT, 2001), which is based on Navier-Stokes momentum equations and VOF method (Hirt and Nichols, 1981) for tracking free water surface. Its original source code was also partly revised to generate bore in the numerical water channel. Numerical raw data computed by CADMAS-SURF included great strong spike phenomena that show the abrupt jumps of wave loads. To resolve this undesired noise of raw data, the band-pass filter with the frequency of 5Hz was utilized. The filtered results showed reasonable agreements with the experimental results performed by Matsutomi (1991) and Ramsden (1996). It was confirmed that CADMASSURF can be applied to the design of coastal structures against tsunami bores. In addition, the transformation process and propagation speed of bores in the same 2-d water channel were discussed by the variations of water level for time and space. The numerical results indicated that the propagation speed of bore was changed due to the nonlinear interactions between negative and reflected waves.

• 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 모형의 적용성을 검토하였으며, 향후 수리실험을 통하여 실제 흐름과의 매개변수 최적화 및 다양한 도로 조건의 변화를 고려한 수치모의 분석을 통하여 지속적인 모형의 검증이 가능할 것으로 판단된다.