• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.60-66
• /
• 2000

In order to improve cavitation characteristics for a high-speed propeller, leading edge shape of a 2-D hydrofoil is investigated numerically and experimentally. For flowfield analysis around the leading edge, the incompressible Reynolds Averaged Navier-Stokes(RANS) equation is solved using the standard $k-\varepsilon$ turbulence model and a finite volume method(FVM). The cavitation thickness, which is occurred on hydrofoil surface, is predicted using the panel code. It is shown that the calculation codes predict the experimental trend fairly well. From these results, new hydrofoils are designed

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.1
• /
• pp.9-16
• /
• 2012

This paper presents the numerical results of a simulation of the free surface flow around a blunt bow ship model and focuses on the validation of the proposed method with a brief investigation of the relation between the resistance and free surface behavior. A finite volume method based on the Reynolds Averaged Navier-Stokes (RANS) approach is used to solve the governing flow equations, where the free surface, including wave breaking,is captured by using a two-phase Level-Set (LS) method. For turbulence closure, a two equation k-${\varepsilon}$ model with the standard wall function technique is used. Finally, the numerical results are compared with the available experimental data, showing good agreement.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.3
• /
• pp.183-188
• /
• 2017

This paper provides numerical results for the estimation of the efficiency of KRISO energy saving devices in the design stage. A finite volume method is used to solve Reynolds averaged Navier-Stokes (RANS) equations, where the SST k-$\omega$ model is selected for turbulence closure. The propeller rotating motion is determined using a rigid body motion (RBM) scheme, which is called a sliding mesh technique. The numerical analysis focuses on predicting the power reduction by the designed KRISO devices (K-DUCT) under a self-propulsion condition. The present numerical results show good agreement with the available experimental data. Finally, it is concluded that CFD can be a useful method, along with model tests, for assessing the performance of energy saving devices for propulsion efficiency improvement.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.6
• /
• pp.427-435
• /
• 2015

This paper provides the numerical results of a simulation of the flow around a propeller working beneath the free surface. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes (URANS) equations, where the wave-making problem is solved using a volume-of-fluid (VOF) method. The numerical analysis focuses on the propeller wake structure affected by the free surface, where we consider another free surface boundary condition that treats the free surface as a rigid wall surface. The propeller wake under the effect of these two free surface conditions shows a reduction in the magnitude of the longitudinal and vertical flow velocities, and its vortical structures strongly interact with the free surface. The thrust and torque coefficient under the free surface effect decrease about 3.7% and 3.1%, respectively. Finally, the present numerical results show a reasonable agreement with the available experimental data.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.2
• /
• pp.100-106
• /
• 1997

Viscous flow around an actual ship is calculated by an use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stress is modelled by using k-$\varepsilon$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the complex boundary of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implcit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the disssssssscretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). The subject ship model of actual calculation is 4,410 TEU class container carrier. For 4 geosim models the calculated viscous resistancce values are compared with the model test results and analyzed on their componentss. The resistance performance of an actual ship is predicted very resonably, so this mothod may be utilized as a design tool of hull form.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.2 s.32
• /
• pp.138-146
• /
• 1999

The viscous flow around a ship hull is calculated by the use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stresses are midelled by using the k-${epsilon}$ turbulence model and the law is applied near the body. Body fitted corrdinates are introduced for the treatment of the complex boundary of the ship hull form and the governing equations in the physical domain transformed into ones in the computational domain. The transformed equations are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implicit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the sidcretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). To assure the proprietty of this computing method, HSVA tanker and Dyne hull are calculated ar both model and ship scale Reynolds number. Their reaults of pressure distributions on fore and aft body, axial velocity contours and transverse velocity velocity vectors and viscous resistance coefficients are compared with other's experiments and calculations.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.4
• /
• pp.232-239
• /
• 2013

This paper provides the structure of a Reynolds Averaged Navier-Stokes(RANS) based simulation method and its validation results for the ship motion problem. The motion information of the hull computed from the equations of motion is considered in the momentum equations as the relative fluid motions with respect to a non-inertial coordinates system. A finite volume method is used to solve the governing equations, while the free surface is captured by using a two-phase level-set method and the realizable k-${\varepsilon}$ model is used for turbulence closure. For the validation of the present numerical approach, the numerical results of the resistance and motion tests for DTMB 5415 at two ship speeds are compared against available experimental data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.627-634
• /
• 2008

In a numerical simulation of open channel turbulent flows using RANS (Reynolds averaged Navier-Stokes) equations model equipped with VOF (Volume of Fluid) scheme, the determination of wall roughness for wall function was studied. The roughness constant, based on the law-of-the-wall for flow on rough walls, obtained by experimental works for pipe flows is employed in general wall functions. However, this constant of wall function is the function of Froude number in open channel flows. Thus, the wall roughness should be determined by taking into account the effect of Froude number. In addition, the wall roughness should be corresponding to Manning's roughness coefficient widely used for open channels. In this study, the relation between wall roughness height as an input condition and Manning's roughness coefficient was investigated, and an equation for effective wall roughness height considering the characteristics of numerical models was proposed as a function of Manning's roughness coefficient.

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.