• Journal of Ship and Ocean Technology
• /
• v.6 no.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
• /
• v.5 no.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.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

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

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• 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.10a
• /
• 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.

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

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

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