• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.1
• /
• pp.69-84
• /
• 2018

In order to analyze the characteristics of propeller exciting force, the hybrid grid is adopted and the numerical prediction of KCS ship model is performed for hull-propeller-rudder system by Reynolds-Averaged Navier Stokes (RANS) method and volume of fluid (VOF) model. Firstly, the numerical simulation of hydrodynamics for bare hull at oblique state is carried out. The results show that with the increasing of the drift angle, the coefficients of resistance, side force and yaw moment are constantly increasing, and the bigger the drift angle, the worse the overall uniformity of propeller disk. Then, propeller bearing force for hull-propeller-rudder system in oblique flow is calculated. It is found that the propeller thrust and torque fluctuation coefficient peak in drift angle are greater than that in straight line navigation, and the negative drift angle is greater than the positive. The fluctuation peak variation law of coefficient of side force and bending moment are different due to various causes.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• 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.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• 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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3B
• /
• pp.323-334
• /
• 2008

Present study examined the functionality of the solitary wave (tsunami) control of the two-rowed porous submerged breakwater by numerical experiments, using a numerical wave tank which is based on the Navier-Stokes equation to explain fluid fields and uses a Volume of Fluid (VOF) method to capture the free water surface. Solitary wave was generated by the internal wave source installed within the computational zone in the numerical wave tank and its wave transformations by structure were compared with those in the previous study. Comparisons with the precious numerical results showed a good agreement. Based on these results, several tow-dimensional numerical modeling investigations of the water fields, including wave transformations, reflection, transmission and energy flux, by the one- and two-rowed permeable submerged breakwater under solitary waves were performed. Even if, it is a research of the limited scope, in case of two-rowed permeable submerged breakwater with $h_0/h=0.925$ ($h_0$ is height of submerged breakwater and h is water depth), the wave height damping in range of $l/L_{eff}>0.4$($L_{eff}$ is effective distance of solitary wave) can reach nearly 60% of the incident wave height. In addition, it is found that reflection coefficient increases nearly 47% and transmission coefficient decreases nearly 18% than one-rowed one. The numerical results revealed that the tow-rowed submerged breakwater can control the incident solitary wave economically and more efficiently than the one-rowed one.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.289-301
• /
• 2009

In this work, wave run-up heights and resultant wave forces on a vertical revetment due to tsunami (solitary wave) are investigated numerically using a numerical wave tank model called CADMAS-SURF (CDIT, 2001. Research and Development of Numerical Wave Channel (CADMAS-SURF). CDIT library, No. 12, Japan.), which is based on a 2-D Navier-Stokes solver, coupled to a volume of fluid (VOF) method. The third order approximate solution (Fenton, 1972. A ninth-order solution for the solitary wave. J. of Fluid Mech., Vol. 53, No.2, pp.257-271) is used to generate solitary waves and implemented in original CADMAS-SURF code. Numerical results of the wave profiles and forces are in good agreements with available experimental data. Using the numerical results, the regression curves determined from the least-square analysis are proposed, which can be used to determine the maximum wave run-up height and force on a vertical revetment due to tsunami. In addition, the capability of CADMAS-SURF is demonstrated for tsunami wave forces acting on an onshore structure using various configuration computations including the variations of the crown heights of the vertical wall and the position of the onshore structure. Based on the numerical results such as water level, velocity field and wave force, the direct effects of tsunami on an onshore structure are discussed.

• 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 the Society of Naval Architects of Korea
• /
• v.54 no.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.

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

• 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의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.