• 한국전산유체공학회지
• /
• 제15권1호
• /
• pp.81-87
• /
• 2010

In this study, we conducted resistance and propulsion performance test of ship composed of the Resistance Test, Propeller Open Water Test and Self Propulsion Test using the CFD(Computational Fluid Dynamics). We used commercial RANS(Reynolds Averaged Navier Stokes equation) solver, as a calculating tool. The unstructured grids were used in a bow and stern of ship, having complex shape, for a convenience of generating grids, and the structured grids were adopted in a central hull and rest of hull having a relatively simple shape which is called hybrid grid method. In addition, The sliding mesh method was adopted to rotate a propeller directly in the Propeller Open Water and Self Propulsion Test. The Resistance Test and Self Propulsion Test were calculated using Volume of Fluid (VOF) model and considering a free surface. And all The three cases were applied realizable k-epsilon model as the turbulence model. The results of calculations were verified for the suitability of calculations by comparing MOERI's EFD results.

• 한국전산유체공학회지
• /
• 제17권2호
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• 한국전산유체공학회지
• /
• 제15권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.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제11권2호
• /
• pp.723-735
• /
• 2019

This paper employs computational tools to investigate the cause of resistance reductions in calm water and waves of the sharp bow form compared to the blunt bow in 66,000 DWT bulk carriers. A more slender shape at the fore-shoulder without a bulbous bow is a prominent feature of the sharp bow. The blunt bow incorporates a bulbous shape. A two-phase unsteady Reynolds averaged Navier-Stokes equations have been solved; and a realizable k-ε model has been applied for the turbulent closure. The free-surface is obtained by solving a VOF equation. The computational results have been validated with model tests carried out at a towing tank. The pressure component of resistance in the sharp bow is reduced by 8.9% in calm water, and 6.4-12.7% in regular head waves. The frictional components of resistance in the sharp and blunt bows are largely the same.

• 한국해양공학회지
• /
• 제22권3호
• /
• pp.24-33
• /
• 2008

In this study, flow and creation and dissipation of vorticity around rectangular floating breakwater is investigated both experimentally and numerically. The PIV system(Particle image velocimetry) is employed to obtain the velocity field in the vorticity of rectangular structure. The numerical model, combined with ${\kappa}-{\varepsilon}$ turbulence model and the VOF method based on RANS equation, is used to analyze the turbulence structure. In the results of this study, the vorticity is found around conner of rectangular structure at all time domain, and creation and dissipation of vorticity are closely related to wave period. Separation points of phase of vortex due to flow separation for longer period waves are faster then for shorter period waves.

• 해양환경안전학회:학술대회논문집
• /
• 해양환경안전학회 2017년도 추계학술발표회
• /
• pp.162-163
• /
• 2017

해양 HNS(Hazardous and Noxious Substances)의 유출 사고 시, 막대한 인명 피해와 환경 훼손을 피하기 위해 유출 사고 조기 예측과 정확한 확산 경로를 예측하는 것이 필수적이다. 본 연구의 최종목적은 전산유체역학을 이용하여 HNS사고가 발생하였을 때 위험구역을 적절히 예측할 수 있는 수치해석기법을 개발하고, 다양한 해양사고조건과 환경영향을 고려하여 근접역에서의 2차원 확산 특성을 고찰하고 확산 현상을 예측하기 위한 모델을 개발하는 것이다. 본 연구에서는 상용코드인 ANSYS FLUENT(V. 17.2)을 사용하여 근접역에서의 2차원 확산특성을 모사하고 분석하였다. 특히, 누출된 HNS의 위치별 농도를 예측하기 위해 종수송방정식(Species Transport Equation)을 이용하였으며 RANS(Reynolds-Averaged Navier-Stokes) 방정식과 표준 $k-{\varepsilon}$ 모델을 이용하여 난류유동을 모사하였다. 해석된 결과는 문헌에서 얻어진 실험데이터와 상호비교하였으며 해수의 유속, HNS의 밀도에 따른 유층 두께, 해수면 HNS 평균 농도 그리고 HNS 전파 속도를 분석하였다. 유층 두께는 해류 유속에 따라 변화하며 변화 경향에 따라 두 구간으로 나눌 수 있다. 해류 전파 속도는 대체로 해류 유속과 선형적 비례관계를 갖는 것으로 나타났다. 해수면 평균 HNS 농도는 해류 유속에 선형적으로 비례하여 감소하며, HNS 밀도가 큰 경우 해수면 평균 HNS 체적 농도는 더 빠르게 감소하게 된다. 이러한 결과는 HNS 확산 특성을 분석하고 관련된 예측모델을 개발하는 데에 기여할 수 있다.

• 한국해양공학회지
• /
• 제37권6호
• /
• pp.225-237
• /
• 2023

In this study, numerical simulations for a single fixed truncated circular cylinder in regular waves were conducted to investigate the nonlinear wave run-up under various dampers and wave period conditions. The present study used the volume of fluid (VOF) technique to capture the air-water interface. The unsteady Reynolds-averaged Navier-Stokes (URANS) equation with the k- 𝜖 turbulence model was solved using the commercial computational fluid dynamics (CFD) software STAR-CCM+. First, a systematic spatial convergence study was conducted to assess the performance and precision of the present numerical wave tank. The numerical scheme was validated by comparing the numerical results of wave run-up on a bare truncated cylinder with the experimental results, and a good agreement was achieved. Then, a series of parametric studies were carried out to examine the wave run-up time series around the truncated cylinder with single and dual dampers in terms of the first- and second-order harmonic and mean set-up components. Additionally, the local wave field and the flow velocity vectors adjacent to the cylinder were evaluated. It was confirmed that under short wave conditions, the high position of the damper led to a noticeable increase in the wave run-ups with significant changes in the first- and second-order harmonic components.

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

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

• 한국수자원학회논문집
• /
• 제43권9호
• /
• pp.801-811
• /
• 2010

급경사에서의 고립파의 처오름을 예측하기 위해 3차원 수치모형에 내부조파기법을 도입하여 수치모형실험을 수행하였다. 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 푸는 LES(large-eddy-simulation) 기반의 수치모형을 사용한다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two-step projection 기법을 사용하여 SANS 방정식을 풀었으며, Poisson 방정식을Bi-CGSTAB 기법을 이용하여 풀고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 고립파를 3차원 공간의 일정 수심상에서 내부조파하여 해석해와 비교한 후 분산오차에 대해 분석하였다. 그리고 고립파를 내부조파하여 급경사에서의 고립파의 처오름 및 처내림 현상을 예측하고 수리모형 실험결과와 비교 및 분석하였다.