• Journal of Ocean Engineering and Technology
• /
• v.21 no.1 s.74
• /
• pp.31-39
• /
• 2007

Recently, the nonlinear dynamic responses among waves, submarine pipeline and seabed have become a target of analyses for marine geotechnical and coastal engineers. Specifically, the velocity field around the submarine pipeline and the wave-induced responses of soil, such as stress and strain inside seabed, have been recognized as dominant factors in discussing the stability of submarine pipeline. The aim of this paper is to investigate nonlinear dynamic responses of soil in seabed, around submarine pipeline, under wave loading. In order to examine wave-induced soil responses, first, the calculation is conducted in the whole domain, including wave field and the seabed, using the VOF-FDM method. Then, velocities and pressures, which are obtained on the boundary between the wave field and the seabed, are used as the boundary condition to compute the wave-induced stress and strain inside seabed, using the poro-elastic FEM model, which is based on the approximation of the Biot's equations. Based on the numerical results, the characteristics of wave-induced soil responses around submarine pipeline are investigated, in detail, inrelation to relative separate distance of the submarine pipeline from seabed. Also, the velocity field around the submarine pipeline is discussed.

• 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
• /
• 2015.05a
• /
• pp.44-44
• /
• 2015

저수지나 하천 사면에서 발생하는 산사태와 토석류는 저수지와 하천 수체에 충격을 가한다. 이로 인해 발생하는 수면 충격파는 전파되어 반대편 제방으로 파의 처오름 또는 댐 제체위로의 물넘이로 큰 피해를 줄 수 있다. 최근 외국에서는 2차원 충격파 생성 및 전파의 기본 과정을 구명하기 위한 실험적 연구가 이뤄지고 있으며, 이들 연구들은 충격파의 발생과 전파, 사면활동 물질과 수체의 상호작용 그리고 자유 수면과 유속분표의 발달에 대한 자세한 관측 자료를 제시하고 있다. 아울러 충격파에 영향을 주는 지배 매개변수를 제시하고 있다. 하지만, 이러한 실험적 연구의 최근 진보에도 불구하고, 이들 지배 매개변수를 고려한 충격파 지배공식들은 대상 지역의 복잡한 바닥 지형이나, 평면적 지형 변화를 단순한 추정치로만 고려하게 된다. 따라서 복잡한 지형조건에서 토석류와 수체의 상호작용과 수면 충격파의 전파를 합리적으로 해석하는 데는 한계가 있다. 이 경우 수치모델링 기법을 대안으로 적용할 수 있으나, 수치모델링은 수면에서 충격파의 전파와 수중에서 토석류의 전파를 동시에 모의해야 하고, 뉴턴 유체와 비뉴턴 유체의 특성을 동시에 고려해야하므로 수치해석 연구자들에게는 하나의 큰 도전사항이다. 이 연구는 경계면 포착기법을 이용한 계산유체동력학 기법을 이용하여 사면활동과 이로 인한 정지 수역에서의 충격파의 발생 및 전파를 재현하기 위한 수치 모델링 기법을 개발하는 것이 목적이다. 사면활동과 수면의 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축 기법을 이용한다. 지배 방정식은 비압축성(incompressible) 질량 보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 큰와 모의 계열의 난류 모델링 기법을 적용하여 충격파의 전파와 붕괴에 대한 난류의 영향을 고려하였다. 토석류는 비뉴턴 흐름저항 관계식을 적용하여 그 흐름특성을 재현하였다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. 외국의 연구자들이 관측하여 제시한 길이 11 m 그리고 폭 0.5 m의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.

• Coupled systems mechanics
• /
• v.8 no.4
• /
• pp.339-350
• /
• 2019

A numerical simulation of the incompressible multiphase hydraulic jump flow was performed to compare the interface prediction through the use of the three RANS turbulence models: $k-{\varepsilon}$, $RNGk-{\varepsilon}$ and SST $k-{\omega}$. A three dimensional no submerged hydraulic jump and a two dimensional submerged hydraulic jump were modeled. Both the geometry and the mesh were created using the open source Gmsh code. The project's geometry consists of a rectangular channel with length and height differences between the two dimensional and three dimensional simulations. Uniform hexahedral cells were used for the mesh. Three refining meshes were constructed to allow to verify simulation convergence. The Volume of Fluid (abbr. VOF) method was used for treatment of the air-water surface. The turbulence models were evaluated in three distinct set up configurations to provide a greater accuracy in the flow representation. In the two-dimensional analysis of a submerged hydraulic jump simulation, the turbulence model RNG RNG $k-{\varepsilon}$ provided a better interface adjust with the experimental results than the model $k-{\varepsilon}$ and SST $k-{\omega}$. In the three-dimensional simulation of a no-submerged hydraulic jump the k-# showed better results than the SST $k-{\omega}$ and RNG $k-{\varepsilon}$ capturing the height and length of the ledge with a better fit with the experimental results.

• Journal of Ocean Engineering and Technology
• /
• v.35 no.1
• /
• pp.38-49
• /
• 2021

Slug flow is the most common multi-phase flow encountered in oil and gas industry. In this study, the hydrodynamic features of flow in pipes investigated numerically using computational fluid dynamic (CFD) simulations for the effect of slug flow on the vertical and bent pipeline. The compressible Reynold averaged Navier-Stokes (RANS) equation was used as the governing equation, with the volume of fluid (VOF) method to capture the outline of the bubble in a pipeline. The simulations were tested for the grid and time step convergence, and validated with the experimental and theoretical results for the main hydrodynamic characteristics of the Taylor bubble, i.e., bubble shape, terminal velocity of bubble, and the liquid film velocity. The slug flow was simulated with various air and water injection velocities in the pipeline. The simulations revealed the effect of slug flow as the pressure occurring in the wall of the pipeline. The peak pressure and pressure oscillations were observed, and those magnitudes and trends were compared with the change in air and water injection velocities. The mechanism of the peak pressures was studied in relation with the change in bubble length, and the maximum peak pressures were investigated for the different positions and velocities of the air and water in the pipeline. The pressure oscillations were investigated in comparison with the bubble length in the pipe and the oscillation was provided with the application of damping. The pressures were compared with the case of a bent pipe, and a 1.5 times higher pressures was observed due to the compression of the bubbles at the corner of the bent. These findings can be used as a basic data for further studies and designs on pipeline systems with multi-phase flow.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1037-1048
• /
• 2023

Fishways are constructed to ensure the fish migration because river-crossing structures such as dams and weirs cut off the stream longitudinal connectivity and influence on aquatic ecosystems. However, the passage efficiency of fishes varies depending on flow characteristics in the fishway and fish species. In this study, three-dimensional numerical simulations are carried out using a RANS model and the volume of fluid method for resolving free surface fluctuations to calculate the turbulent flow in the pool and weir type fishway. The Flow velocity and turbulent kinetic energy in the pool of fishway are analyzed according to variation of the upstream water level and the length of pool. The present numerical simulations reasonably well reproduce the stream flow and plunging flow characteristics in the pool. The simulation results show that the stream flow changes to the plunging flow as the length of the pool increases. When the upstream level increases, the stream flow becomes more evident. Key parameters related to the fish migration within the fishway such as the flow velocity and the turbulent kinetic energy are examined to assess the ascending possibility of fishes.

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

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