• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.435-449
• /
• 2014

The effects of the gravity field and the free surface on the cavity shape and the drag are investigated through a numerical analysis for the steady supercavitating flow past a simple two-dimensional body underneath the free surface. The continuity and the RANS equations are numerically solved for an incompressible fluid using a $k-{\epsilon}$ turbulence model and a mixture fluid model has been applied for calculating the multiphase flow of air, water and vapor using the method of volume of fluid and the Schnerr-Sauer cavitation model. Numerical solutions have been obtained for the supercavitating flow about a two-dimensional $30^{\circ}$ wedge in wide range of depths of submergence and inflow velocities. The results are presented for the cavity shape, especially the length and the width, and the drag of the wedge in comparison with those of the case for the infinite fluid flow neglecting the gravity and the free surface. The influences of the gravity field and the free surface on the aforementioned quantities are discussed. The length and the width of the supercavity are reduced and the centerline of the cavity rises toward the free surface due to the effects of the gravity field and the free surface. The drag coefficient of the wedge, however, is about the same except for shallow depths of submergence. As the supercavitating wedge is approaching very close to the free surface, it is found the length and the width of a cavity are shorten even though the cavitation number is reduced. Also the present result suggests that, under the influence of the gravity field and the free surface, the length of the supercavity for a certain cavitation number varies and moreover is proportional to the inverse of the submergence depth Froude number.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.5
• /
• pp.355-363
• /
• 2013

Panel methods are essential tools for analyzing a fluid-flow problem around complex three dimensional bodies, but they lack ability to solve viscous effects. On the other hand, CFD methods are considered as powerful tools for analyzing fluid-flow characteristics including viscosity. However, they also have short falls, requiring more computing time and showing different results depending on the selection of turbulence models and grid systems. In this paper a hybrid scheme combining a panel method and a CFD method is suggested. The scheme adopts a panel method for far-field solution where viscous effects are negligible and a CFD method for the solution of RANS equations in near-field where viscous effects are relatively strong. The intermediate region between the far-field and near-field is introduced where the calculated field point velocities by the panel method are given as boundary velocities for the CFD method. To verify the scheme, calculated results, by a panel method, a CFD method and the hybrid scheme, for a two dimensional foil section are compared. The suggested hybrid scheme gives reasonable results, while computation time and memory can be dramatically reduced. By using the hybrid scheme efforts can be concentrated for the local flow near the leading and trailing edges, by providing more dense grid system, where detailed flow characteristics are required.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.2
• /
• pp.147-157
• /
• 2011

A RANS(Reynolds averaged Navier-Stokes) based numerical method is developed for the evaluation of ship resistance and self-propulsion performances. In the usability aspect of CFD for the hull form design, the field grid around practical hull forms is generated by solving a grid Poisson equation based on the hull surface grid generated from station offsets and centerline profile. A body force technique is introduced to model the effects of the propeller in which the propeller loads are obtained from potential flow analysis using an unsteady lifting surface method. The free surface is captured by using a two-phase level-set method and the realizable $k-{\varepsilon}$ model is used for turbulence closure. The hull attitude in vertical plane, i.e., trim and sinkage, is calculated by using a quasi-steady method and then considered in the computation by translating and rotating the grid system according to the values. For the validation of the proposed method, the numerical results of resistance tests for KCS, KLNG, and KVLCC1 and of self-propulsion test for KCS are compared with experimental data.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.5
• /
• pp.60-66
• /
• 2012

One of commonly physical phenomena encountered in pump sump systems in which its significant influence to the hydraulic performance of pump system plays an important role in the field of fluid engineering, is the appearance of free surface and submerged vortices. In this paper, a study of the vortices behavior and their formative mechanism of asymmetry is considered in this paper by using numerical approach. The Reynolds-Averaged Navier-Stokes (RANS) equations and k-omega Shear Stress Transport turbulence model used to describe the properties of turbulent flows, in company with VOF multiphase model, are implemented by Fluent code with multi-block structured grid system. In the numerical simulation, the calculated elevation of air-water interface and vortex core contours are used to classify visually surface vortices as well as submerged vortices. It is shown that the free surface vortex is identified by the concavity of liquid region from the free surface and swirling flow at that own plane. To investigate the distinctive behavior of these vortices corresponding to each given flow rate at the same water level, some numerical testing of them are considered here in such a manner that the flow pattern of surface vortex are obtained similarly to the obtained results from experiment. Furthermore, the influence due to the change of grid refinement and the variation of depth of the concavity are also considered in this paper. From that, these influential factors will be implemented to design a good pump sump with higher performance in the future.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.1-7
• /
• 2014

The present study numerically investigated the effect of the advance ratio on the wake characteristics of the marine propeller in the propeller open water test. Therefore, a wide range of the advance ratio(0.2${\kappa}-{\omega}$SST Model are considered. The three-dimensional vortical structures of tip vortices are visualized by the swirl strength, resulting in fast decay of the tip vortices with increasing the advance ratio. Furthermore, to better understanding of the wake evolution, the contraction ratio of the slip stream for different advance ratios is extracted from the velocity fields. Consequently, the slip stream contraction ratio decreases with increasing the advance ratio and successively the difference of the slip stream contraction ratio between J=0.2 and J=0.8 is about 0.1R.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.281-281
• /
• 2017

하천에서 하도불안정(stream instability)으로 인하여 하상의 형태가 변화하고 하상파(sand wave)가 발생한다. 사련(ripple), 사구(dune) 등과 같은 하상파는 흐름저항을 유발하여 홍수시 수위를 증가시킨다. 수리실험 및 수치모형을 이용하여 사련 및 사구의 발달과정 그리고 이를 지나는 난류흐름에 대한 연구가 국외에서는 이루어지고 있지만 국내의 경우 거의 찾아보기 힘들다. 수치모형을 활용한 연구는 주로 횡방향으로 하상파가 일정하다는 가정하에 연직 2차원 수치모형을 적용하였으나 최근 컴퓨터 기술 및 수치기법의 고도화로 3차원 RANS(Reynolds averaged Navier-Stokes) 또는 LES(Large eddy simulation)를 이용한 수치모형이 개발되고 있다. 본 연구에서는 LES에 유사이송 및 하상변동 모형과 결합하여 사구발달에 대한 수치모의를 수행하였다. LES와 유사이송 및 하상변동 모형의 결합은 순간유속성분을 하상변동모형에 직접 적용되기 때문에 난류영향을 고려할 수 있는 것이 장점이다. 특히 사구의 발달에 따라 복잡한 흐름이 발생하며 3차원 와구조가 발생하므로 난류특성의 고려는 필수적이다. 수치모의는 Delft Hydraulics (Bakker et al., 1986)에서 수행한 수리실험 T39를 활용하였다. 수리실험은 길이 100 m, 폭 0.5 m 개수로에서 수행되었으며 평균유속은 0.611 m/s, 수심은 0.436 m이다. 하상파 실험에 사용된 유사입경은 0.78 mm 균일사를 사용하였다. 수치모의 조건은 수리실험과 동일하게 하였으나 계산시간의 효율을 고려하여 흐름방향의 계산영역은 4.0 m로 하고 주기경계조건(periodic boundary condition)을 부여하여 계산을 수행하였다. 수치모의 계산은 사구의 길이 및 파고가 평형상태에 이를 때까지 수행되었다. 수치모의 통해 사구발달에 따른 흐름 및 하상변동 특성을 분석하였다.

• Journal of Ship and Ocean Technology
• /
• v.12 no.3
• /
• pp.36-54
• /
• 2008

A jet stream applied tangential to a curved surface in fluid increases lift force by strengthening circulation around the surface and this phenomenon is known as the Coanda effect. Many experimental and numerical studies have been performed on the Coanda effect and the results found to be useful in various fields of aerodynamics. Recently, preliminary studies on Coanda control surface are in progress to look for practical application in marine hydrodynamics since various control surfaces are used to control behaviors of ships and offshore structures. In the present study, the performance of a Coanda control surface with different geometries of the jet injection nozzle was surveyed to assess applicability to ship rudders. A numerical simulation was carried out to study flow characteristics around a section of a horn type rudder subjected to a tangential jet stream. The RANS equations, discretized by a cell-centered finite volume method were used for this computation after verification by comparing to the experimental data available. Special attentions have been given to the sensitivity of the lift performance of a Coanda rudder to the location of the slit (outlet) and intake of the gap between the horn and rudder surface at the various angles of attack. It is found that the location of the water intake is important in enhancing the lift because the gap functions as a conduit of nozzle generating a jet sheet on the rudder surface.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.4
• /
• pp.209-220
• /
• 2019

Even though the scale of hard structures for beach stabilization should carefully be determined such that these structures do not interrupt the great yearly circulation process of beach sediment in which the self-healing ability of natural beach takes places, massive hard structures such as the submerged breakwater of wide-width are frequently deployed as the beach stabilization measures. On this rationale, asymmetric ripple mat by Irie et al. (1994) can be the alternatives for beach stabilization due to its small scale to replace the preferred submerged breaker of wide-width. The effectiveness of asymmetric ripple mat is determined by how effectively the vortices enforced at the contraction part of flow area over the mat traps the sediment moving toward the offshore by the run-down. In order to verify this hypothesis, we carry out the numerical simulations based on the Navier-Stokes equation and the physically-based morphology model. Numerical results show that the asymmetric ripple mat effectively capture the sediment by forced vortex enforced at the apex of asymmetric ripple mat, and bring these trapped sediments back to the beach, which has been regarded to be the driving mechanism of beach stabilization effect of asymmetric ripple mat.