• 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 the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.8-17
• /
• 2002

In this paper, the computational analysis of open-water characteristics for three model propellers(P4119, P4842 and 3 podded propeller of KRISO) is done by using a viscous-flow method based on Reynolds-Averaged Navier-Stokes equations. The results are presented for open-water performances, blade-section pressures, and circumferentially-averaged velocity profiles for the all three propeller models. Overall close agreements with available experimental data are shown. However, some discrepancies are also found in the pressure near the leading edge of the propeller blade and the open-water performance of the podded propellers.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.2
• /
• pp.109-123
• /
• 2008

The present paper presents the CFD result for a beam wave test case. An ONR tumblehome ship model with bilge keels is used. The beam wave test is for zero forward speed and roll and heave 2DOF with wave slope $a_k=0.156$ and wavelength ${\lambda}=1.12L_{PP}$, with $L_{PP}$ the ship length. The problems is solved numerically with an unsteady Reynolds averaged Navier-Stokes approach. The free surface flow is computed using a single-phase level-set method and the motions in each time step are integrated using a predictor-corrector iteration approach which uses dynamic overset grids moving with relative ship motion. The predicted CFD results for motions and forces are compared with experimental data, showing a reasonable agreement.

• Journal of Ship and Ocean Technology
• /
• v.7 no.2
• /
• pp.40-47
• /
• 2003

LDV measurements in large cavitation tunnel around a propeller in operation are carried out to provide valuable information for more accurate wake-adapted propeller design and to study hull-propeller interactions. Effective velocities are computed by both the simplified vortex ring method and by RANS solver with the body force representing the propeller load. The former method uses the nominal velocities measured at the propeller plane as an input data of the numerical method and shows a better agreement with experimental data. The latter shows the qualitative agreement and may be used as an alternative design tools in the preliminary design stage.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.28-40
• /
• 2002

Turbulent flows around two-dimensional wing sections in ground effect are analysed by incompressible RANS equations and a finite difference method. The Baldwin-Lomax algebraic turbulence model is used to simulate high Reynolds number flows. The main purpose of this study is to clarify the two-dimensional ground effect and its flow characteristics due to different ground boundary conditions, i.e., moving and fixed bottom boundary. As a first step, to validate the present numerical code, the computational result of Clark-Y(t/C 11.7%) is compared with published numerical results and experimental data. Then, NACA4412 section in ground effect is calculated for various ground clearances with two bottom boundary conditions. According to the computational results, the difference in the lift and moment simulated with the two bottom boundary conditions is negligible, but the drag force simulated by the fixed bottom is to some extent smaller than that by the moving bottom. Therefore, it can be concluded that the drag force measured in a wind tunnel with the fixed bottom could be smaller than that with the moving bottom.

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

• Wind and Structures
• /
• v.18 no.5
• /
• pp.529-548
• /
• 2014

This article presents the unsteady aerodynamic performance of crosswind stability obtained numerically for the ATM train. Results of numerical investigations of airflow past a train under different yawing conditions are summarized. Variations of occurrence flow angle from parallel to normal with respect to the direction of forward train motion resulted in the development of different flow patterns. The numerical simulation addresses the ability to resolve the flow field around the train subjected to relatively large yaw angles with three-dimensional Reynolds-averaged Navier-Stokes equations (RANS). ${\kappa}-{\varepsilon}$ turbulence model solved on a multi-block structured grid using a finite volume method. The massively separated flow for the higher yaw angles on the leeward side of the train justifies the use of RANS, where the results show good agreement with verification results. A method of solution is presented that can predict all aerodynamic coefficients and the wind characteristic curve at variety of angles at different speed.

• 유체기계공업학회:학술대회논문집
• /
• 2004.12a
• /
• pp.437-442
• /
• 2004

The optimization of a blade sweep for a transonic axial compressor rotor (NASA rotor 37) has been performed using a response surface method and a Reynolds-averaged Wavier-Stokes (RANS) flow simulation. Two shape variables of the rotor blade, which are used to define a blade sweep, are introduced to increase an adiabatic efficiency. Data points for response evaluations have been selected by D-optimal design, and linear programming method has been used for an optimization on a response surface. The result shows that the adiabatic efficiency is increased to about 1 percent compared to that of the reference shape of the rotor blade. Relatively high increasement of the adiabatic efficiency is obtained between 20 and 60 percent span. In the present study, backward swept blade is more effective to increase the adiabatic efficiency In the axial compressor rotor.