• 한국해양공학회지
• /
• 제35권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.

• 대한조선학회논문집
• /
• 제54권3호
• /
• pp.215-226
• /
• 2017

This paper presents numerical results of the performance of a marin propeller in cavitating and non-cavitating flow conditions. The geometry and experimental validation data of the propeller are provided in Potsdam Propeller Test Case(PPTC) in the framework of the second International Symposium on Marine Propulsors 2011(SMP'11) workshop. The PPTC includes open water tests, velocity field measurements and cavitation tests. The present numerical analysis was carried out by using the Reynolds averaged Navier-Stokes(RANS) method on a wall-resolved grid ensuring a y+=1, where the SST k-${\omega}$ model was mainly used for turbulence closure. The influence of the turbulence model was investigated in the prediction of the wake field under a non-cavitating flow condition. The propeller tip vortex flows in both cavitating and non-cavitating conditions were captured through adaptation of additional grids. For the cavitation flows at three operation points, Schnerr-Sauer's cavitation model was used with a Volume-Of Fluid(VOF) approach to capture the two-phase flows. The present numerical results for the propeller wake and cavitation predictions including the open water performance showed a qualitatively reasonable agreement with the model test results.

• 한국해양공학회지
• /
• 제35권1호
• /
• pp.50-58
• /
• 2021

The demand for eco-friendly energy is expected to increase due to the recently strengthened environmental regulations. In particular, the flow inside the pipe used in a cargo handling system (CHS) or fuel gas supply system (FGSS) of hydrogen transport ships and hydrogen-powered ships exhibits a very complex pattern of multiphase-thermal flow, including the boiling phenomenon and high accuracy analysis is required concerning safety. In this study, a feasibility study applying the boiling model was conducted to analyze the multiphase-thermal flow in the pipe considering the phase change. Two types of boiling models were employed and compared to implement the subcooled boiling phenomenon in nucleate boiling numerically. One was the "Rohsenow boiling model", which is the most commonly used one among the VOF (Volume-of-Fluid) boiling models under the Eulerian-Eulerian framework. The other was the "wall boiling model", which is suitable for nucleate boiling among the Eulerian multiphase models. Moreover, a comparative study was conducted by combining the nucleate site density and bubble departure diameter model that could influence the accuracy of the wall boiling model. A comparison of the Rohsenow boiling and the wall boiling models showed that the wall boiling model relatively well represented the process of bubble formation and development, even though more computation time was consumed. Among the combination of models used in the wall boiling model, the simulation results were affected significantly by the bubble departure diameter model, which had a very close relationship with the grid size. The present results are expected to provide useful information for identifying the characteristics of various parameters of the boiling model used in CFD simulations of multiphase-thermalflow, including phase change and selecting the appropriate parameters.

• 한국해양공학회지
• /
• 제25권2호
• /
• pp.7-14
• /
• 2011

Numerical simulations were performed for the evaluation of wave and current loads on a fixed cylindrical substructure model for an ocean wind turbine using the ANSYS-CFX package. The numerical wave tank was actualized by specifying the velocity at the inlet and applying momentum loss as a wave damper at the end of the wave tank. The Volume-Of-Fluid (VOF) scheme was adopted to capture the air-water interface. An accuracy validation of the numerical wave tank with a truncated vertical circular cylinder was accomplished by comparing the CFD results with Morison's formula, experimental results, and potential flow solutions using the higher-order boundary element method (HOBEM). A parametric study was carried out by alternately varying the length and amplitude of the wave. As a meaningful engineering application, in the present study, three kinds of conditions were considered, i.e., cases with current, waves, and a combination of current and progressive waves, passing through a cylindrical substructure model. It was found that the CFD results showed reasonable agreement with the results of the HOBEM and Morison's formula when only progressive waves were considered. However, when a current was included, CFD gave a smaller load than Morison's formula.

• 대한기계학회논문집B
• /
• 제39권11호
• /
• pp.871-884
• /
• 2015

유체이송 기술은 마이크로 유체시스템 개발에서 핵심문제로 인식되고 있다. 최근 명(2014)은 외부동력을 사용하지 않고 액적을 이동시킬 수 있는 새로운 개념을 제안하고, 초기에 반원통형 형상을 가지는 가상의 2차원 액적에 대한 수치해석을 통해 이 개념이 성립함을 보였다. 또한 명과 권(2015)은 친수성/소수성 표면위에서 초기 3차원 반구 형상의 실제 물 액적이송의 메커니즘을 시간에 따른 액적형상과 액적 내부의 운동에너지, 중력에너지, 표면자유에너지 및 압력에너지의 수치해석 결과를 통해 규명하였다. 본 연구는 새로운 개념을 확립시키기 위해 초기 구형액적에 대한 3차원 수치해석을 수행하고, 액적이송의 메커니즘을 모세관력 힘의 불균형 관점에서 액적 형상과 다양한 에너지의 수치해석 결과를 통해 규명하였다.

• 대한토목학회논문집
• /
• 제36권6호
• /
• pp.1037-1048
• /
• 2016

본 연구에서는 Navier-Stokes 방정식과 자유수면 추적에 VOF법을 채용하는 CADMAS-SURF를 이용하여 파랑과 흐름의 비선형 상호간섭현상을 연구하였다. 파-흐름 공존장에서의 유체거동 해석을 위해 CADMAS-SURF를 수정 및 확장하였고, 계산치를 실험치와 비교하여 본 연구의 타당성을 확인할 수 있었다. 본 연구의 수치수로를 파랑-흐름 공존장에서 주어지는 복잡한 물리 현상 규명과 파-흐름-잠제의 상호간섭해석에 적용하여, 유속장, 와도장, 자유수면과 와도와의 관계 등을 논의하였다.

• Coupled systems mechanics
• /
• 제8권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.

• 대한기계학회논문집B
• /
• 제29권6호
• /
• pp.671-686
• /
• 2005

The transient nature and complex geometries of two-phase gas-liquid flows cause fundamental difficulties when measuring flow velocity using an electromagnetic flowmeter. Recently, a current-sensing flowmeter was introduced to obtain measurements with high temporal resolution (Ahn et al.). In this study, current-sensing flowmeter theory was applied to measure the fast velocity transients in slug flows. The velocity fields of axisymmetric gas-liquid slug flow in a vertical pipe were obtained using Volume-of-Fluid (VOF) method, and the virtual potential distributions for the electrodes of finite size were also computed using the finite volume method for simulating slug flow. The output signal prediction for slug flow was carried out from the velocity and virtual potential (or weight function) fields. The flowmeter was numerically calibrated to obtain the cross-sectional liquid mean velocity at an electrode plane from the predicted output signal. Two calibration parameters are proposed for this procedure: a flow pattern coefficient and a localization parameter. The flow pattern coefficient was defined by the ratio of the liquid resistance between the electrodes for two-phase flow with respect to that for single-phase flow, and the localization parameter was introduced to avoid errors in the flowmeter readings caused by liquid acceleration or deceleration around the electrodes. These parameters were also calculated from the computed velocity and virtual potential fields. The results can be used to obtain the liquid mean velocity from the slug flow signal measured by a current-sensing flowmeter.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1951-1956
• /
• 2004

The transient nature and complex flow geometries of two-phase gas-liquid flows cause fundamental difficulties when measuring flow velocity using an electromagnetic flowmeter. Recently, a current-sensing flowmeter was introduced to obtain measurements with high temporal resolution (Ahn et $al.^{(1)}$). In this study, current-sensing flowmeter theory was applied to measure the fast velocity transients in slug flows. To do this, the velocity fields of axisymmetric gas-liquid slug flow in a vertical pipe were obtained using Volume-of-Fluid (VOF) method and the virtual potential distributions for the electrodes of finite size were also computed using the finite volume method for the simulated slug flow. The output signal prediction for slug flow was carried out from the velocity and virtual potential (or weight function) fields. The flowmeter was numerically calibrated to obtain the cross-sectional liquid mean velocity at an electrode plane from the predicted output signal. Two calibration parameters are required for this procedure: a flow pattern coefficient and a localization parameter. The flow pattern coefficient was defined by the ratio of the liquid resistance between the electrodes for two-phase flow with respect to that for single-phase flow, and the localization parameter was introduced to avoid errors in the flowmeter readings caused by liquid acceleration or deceleration around the electrodes. These parameters were also calculated from the computed velocity and virtual potential fields. The results can be used to obtain the liquid mean velocity from the slug flow signal measured by a current-sensing flowmeter.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2015년도 학술발표회
• /
• pp.44-44
• /
• 2015

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