• Journal of Ship and Ocean Technology
• /
• 제6권2호
• /
• pp.37-50
• /
• 2002

In the present numerical simulation of viscous free surface flow around a ship, two-fluids in-compressible Reynolds-averaged Navier-Stokes equations with the standard $\textsc{k}-\varepsilon$turbulence model are discretized on a regular grid by using a finite volume method. A local level-set method is introduced for capturing the free surface movement and the influence of the viscous layer and dynamic boundary condition of the free surface are implicitly considered. Partial differential equations in the level-set method are discretized with second order ENO scheme and explicit Euler scheme in the space and time integration, respectively. The computational results for the Series-60 model with $C_B=0.6$ show a good agreement with the experimental data, but more validation studies for commercial complicated hull forms are necessary.

• Journal of Ship and Ocean Technology
• /
• 제5권1호
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• 한국수자원학회논문집
• /
• 제52권10호
• /
• pp.697-706
• /
• 2019

다상흐름 모델링 기법과 하이브리드 난류 모델링 기법을 결합한 수치모형을 이용하여 사각형 수로에서의 중력류를 수치모의 하였다. 이 연구에서 적용한 다상흐름 해석기법은 밀도가 큰 중력류 유체, 상대적으로 밀도가 작은 주변류 유체 그리고 자유수면 위에서 흐르는 공기를 3개의 상으로 처리하며, 각 상에 대해서 분리된 흐름 지배방정식을 적용한다. 난류흐름은 벽경계 근처에서는 RANS 모드로 모의하고 벽에서 떨어진 영역에서는 LES 모드로 해석하는 하이브리드 RANS/LES 방법의 일종인 IDDES 기법을 이용하여 해석한다. 이 연구에서 적용한 모델링 기법은 중력류의 머리의 전파속도를 실험값과 일치하게 잘 예측하는 것으로 나타났다. 수치해석 결과는 아울러 낮은 레이놀즈수 난류모형을 이용한 RANS 수치모의에서 이용되는 정도의 격자해상도에서도 큰 규모의 Kelvin-Helmholtz 형식의 경계면 와의 발달과 이들 와가 지속적으로 3차원 형식의 붕괴를 거쳐 작은 난류구조로 분해되면서 난류에너지가 소산되는 현상을 성공적으로 예측함을 보여준다. 적용한 수치모의 기법은 공학적으로 접근 가능한 격자해상도에서 돌출-쪼개짐 흐름 불안정을 동반한 중력류 머리부분의 3차원 거동 특성을 잘 재현하며, 이 결과는 보다 높은 격자해상도에서 구해진 LES 결과에 상응하는 것으로 나타났다. 이 연구결과는 하이브리드 난류모델링 기법과 다상흐름 해석기법을 병합한 수치모형이 자연상태에서 복잡한 중력류의 물리적 거동을 예측하는데 공학적으로 유망한 방법임을 보여준다.

• 대한토목학회논문집
• /
• 제28권6B호
• /
• pp.627-634
• /
• 2008

개수로의 난류흐름을 VOF(Volume of Fluid)기법을 채용한 RANS(Reynolds averaged Navier-Stokes) 방정식 모형을 사용하여 수치모의할 때 벽면함수의 거칠기를 산정하기 위해 고려해야 하는 점들을 연구하였다. 거친 벽면상의 흐름을 위한 벽면함수의 거칠기 상수(roughness constant)는 관수로 흐름의 실험을 통하여 얻어진 값을 사용한다. 그러나 개수로 흐름에서는 이 거칠기 상수가 Froude 수에 따라 변화하므로 이를 고려할 수 있어야 하며, 개수로에서 광범위하게 사용되는 Manning 조도계수에 상응하는 벽면 거칠기 높이(roughness height)를 산정하여 사용할 수 있어야 있다. 본 연구에서는 모형에 입력되는 벽면함수의 거칠기 높이와 Manning 조도계수 사이의 관계를 분석하였다. 이를 바탕으로 수치모형의 특성이 고려되고 Manning 조도계수의 함수로 표현되는 유효 거칠기 높이 산정식을 제안하였다.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2008년도 정기총회 및 학술발표대회
• /
• pp.687-690
• /
• 2008

Currently, Stream flow analysis has been accomplished by one or two dimensional equations and was applied by simple momentum equations and fixed energy conservations which contain many reach uppermost limit. In this study, FLOW-3D using CFD(Computational Fluid Dynamics) was applied to stream flow analysis which can solve three dimensional RANS(Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behavior and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as $k-\varepsilon$, RNG(Renomalized Group Theory) $k-\varepsilon$ and LES(Large Eddy Simulation). Numerical analysis results have been illustrated by the turbulence energy effects, velocity of flow, water level pressure and eddy flows around the side overflow type structures at Jangwall bridge in urban stream.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.137-139
• /
• 2003

As a part of numerical and experimental research works for the prediction and improvement of ship's maneuvering performance, a study on the performance analysis of two different rudders has been carried out. While the planform shape and the aspect ratio of the rudders have been fixed, section shape has been changed. Conventional type of HMRI NP section and special type of dolphin-tail section have been employed. Performances of the rudders have been investigated by using CFD and compared with experimental data obtained in a wind tunnel. A commercial CFD program has been used to solve the RANS equations. Two-equation k-ro model has been applied to close the governing equations. Block-structured grids are used in the numerical calculation. Based upon the calculation results, the rudder with dolphin-tail section has shown a possibility of significantly improving rudder performance if utilized as the section of ship rudders.

• 대한조선학회논문집
• /
• 제50권5호
• /
• pp.282-290
• /
• 2013

The effects of the combination of blade number for forward and after propeller on the propeller shaft forces of a contra-rotating propeller (CRP) system are presented in the paper. The research is performed through the numerical simulations based on the Reynolds-Averaged Navier-Stokes equations (RANS). The simulation results of the present method in open water condition are validated comparing with the experimental data as well as the other numerical simulation results based on the potential method for 4-0-4 CRP (3686+3687A) and 4-0-5 CRP (3686+3849) of DTNSRDC. Two sets of CRP are designed and simulated to study the effect of the combination of blade number in behind-hull condition. One set consists of 3-blade and 4-blade, while the other is 4-blade and 4-blade. A full hull body submerged under the free surface is modeled in the computational domain to simulate directly the wake field of the ship at the propeller plane. From the simulation results, the fluctuations of axial force and moment are dominant in the case of same blade numbers for forward and after propellers, whereas the fluctuations of horizontal and vertical forces and moments are very large in the case of different blade numbers.

• Journal of Ship and Ocean Technology
• /
• 제10권4호
• /
• pp.12-23
• /
• 2006

In this paper, a numerical study is carried out to investigate the turbulent flow around a twin-skeg container ship model with rudders including propeller effects. A commercial CFD code, FLUENT is used with body forces distributed on the propeller disk to simulate the ship stem and wake flows with the propeller in operation. A multi-block, matching, structured grid system has been generated for the container ship hull with twin-skegs in consideration of rudders and body-force propeller disks. The RANS equations for incompressible fluid flows are solved numerically by using a finite volume method. For the turbulence closure, a Reynolds stress model is used in conjunction with a wall function. Computations are carried out for the bare hull as well as the hull with appendages of a twin-skeg container ship model. For the bare hull, the computational results are compared with experimental data and show generally a good agreement. For the hull with appendages, the changes of the stem flow by the rudders and the propellers have been analyzed based on the computed result since there is no experimental data available for comparison. It is found the flow incoming to the rudders has an angle of attack due to the influence of the skegs and thereby the hull surface pressure and the limiting streamlines are changed slightly by the rudders. The axial velocity of the propeller disk is found to be accelerated overall by about 35% due to the propeller operation with the rudders. The area and the magnitude of low pressure on the hull surface enlarge with the flow acceleration caused by the propeller. The propellers are found to have an effect on up to the position where the skeg begins. The propeller slipstream is disturbed strongly by the rudders and the flow is accelerated further and the transverse velocity vectors are weakened due to the flow rectifying effect of the rudder.

• 대한조선학회논문집
• /
• 제53권6호
• /
• pp.494-502
• /
• 2016

The accurate assessment of hull-appendage interaction in the early design stage is important to control the inflow to the propeller plane, which can cause undesirable hydrodynamic effects in terms of cavitation phenomenon. This paper describes a numerical analysis for the flow around a fully appended surface ship model for which KRISO has carried out a model test in the Large Cavitation Tunnel(LCT). This numerical study was performed with the LCT model test in a complementary manner for a good reproduction of the wake distribution of surface ships. A second order accurate finite volume method provided by a commercial computational fluid dynamics(CFD) program was used to solve the governing Reynolds Averaged Navier-Stokes(RANS) equations, where the SST $k-{\omega}$ model was used for turbulence closure. The numerical results were compared to available LCT experimental data for validation. The calculations gave good predictions for the boundary layer profiles on the walls of the empty cavitation tunnel and the wake at the propeller plane of the fully appended hull model in the LCT.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2007년도 춘계학술대회B
• /
• pp.2690-2695
• /
• 2007

This work presents numerical optimization for design of a blade stacking line of a low speed axial flow fan with a fast and elitist Non-Dominated Sorting of Genetic Algorithm (NSGA-II) of multi-objective optimization using three-dimensional Navier-Stokes analysis. Reynolds-averaged Navier-Stokes (RANS) equations with ${\kappa}-{\varepsilon}$ turbulence model are discretized with finite volume approximations and solved on unstructured grids. Regression analysis is performed to get second order polynomial response which is used to generate Pareto optimal front with help of NSGA-II and local search strategy with weighted sum approach to refine the result obtained by NSGA-II to get better Pareto optimal front. Four geometric variables related to spanwise distributions of sweep and lean of blade stacking line are chosen as design variables to find higher performed fan blade. The performance is measured in terms of the objectives; total efficiency, total pressure and torque. Hence the motive of the optimization is to enhance total efficiency and total pressure and to reduce torque.