• Journal of Ship and Ocean Technology
• /
• v.7 no.4
• /
• pp.16-31
• /
• 2003

The finite volume based multi-block RANS code, WAVIS developed at KRISO, is used to simulate the turbulent flows around a submarine with the realizable $\textsc{k}-\varepsilon$ turbulence model. RANS methods are verified and validated at the level of validation uncertainty 1.54% of the stagnation pressure coefficient for the solution of the turbulent flows around SUBOFF submarine model without appendages. Another SUBOFF configuration, axisymmetric body with four identical stem appendages, is also computed and validated with the experimental data of the nominal wake and hydrodynamic coefficients. The hydrodynamic forces and moments for SUBOFF model and a practical submarine are predicted at several drift and pitch angles. The computed results are in extremely good agreement with experimental data. Furthermore, it is noteworthy that all the computations at the present study were carried out in a PC and the CPU time required for 2.8 million grids was about 20 hours to get fully converged solution. The current study shows that CFD can be a very useful and cost effective tool for the prediction of the hydrodynamic performance of a submarine in the basic design stage.

• Journal of Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.

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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1287-1292
• /
• 2007

Recently, the frequency of unexpecting heavy rains has been increased due to abnormal climate and extreme rainfall. There was a limit to analyze 1D or 2D stream flow that was applied simple momentum equation and fixed energy conservation. Therefore, hydrodynamics flow analysis in rivers has been needed 3D numerical analysis for correct stream flow interpretation. In this study, CFD model on FLOW-3D was applied to stream flow analysis, which solves three dimenson 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 ${\kappa}-{\varepsilon}$, RNG ${\kappa}-{\varepsilon}$ and LES. Those numerical analysis results have been illustrated by the turbulence energy effects, velocity of flow distributions, water level pressure distributions and eddy flows around the piers at Jangwall bridge in urbarn stream.

• Journal of Ship and Ocean Technology
• /
• v.8 no.1
• /
• pp.10-30
• /
• 2004

High-fidelity RANS simulations are presented for a ducted marine propulsor, including verification ＆ validation (V＆V) using available experimental fluid dynamics (EFD) data. CFDSHIP-IOWA is used with $\textsc{k}-\omega$ turbulence model and extensions for relative rotating coordinate system and Chimera overset grids. The mesh interpolation code PEGASUS is used for the exchange of the flow information between the overset grids. Intervals V＆V for thrust, torque, and profile averaged radial velocity just downstream of rotor tip are reasonable in comparison with previous results. Flow pattern displays interaction and merging of tip-leakage and trailing edge vortices. In interaction region, multiple peaks and vorticity are smaller, whereas in merging region, better agreement with EFD. Tip-leakage vortex core position, size, circulation, and cavitation patterns for $\sigma=5$ also show a good agreement with EFD, although vortex core size is larger and circulation in interaction region is smaller.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.4
• /
• pp.1-10
• /
• 2001

A numerical solution method of the incompressible Reynolds-Averaged Navier-Stokes equations is applied for calculating turbulent flows and performances of a marine propeller in open-water, four-quadrant conditions. Computed propeller flows of the model propeller P4381, for which the experimental data of the open-water performances exist, reveal complex viscous-flow characteristics including three-dimensional flow separations in various off-design conditions and also computed propeller thrusts and torques agree quite well with experimental data except some cases for which severe propeller cavitations occurred in the experiment.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.55-61
• /
• 2008

The Reynolds-Averaged Navier-Stokes(RANS) analysis of the 3-D steady flow around the NREL Phase VI horizontal axis wind turbine(HAWT) rotor was performed. The CFD analysis results were compared with experimental data at several different wind speeds. The present CFD model shows good agreements with the experiments both at low wind speed which formed well-attache flow mostly on the upper surface of the blade, and at high wind speed which blade surface flow completely separated. However, some discrepancy occurs at the relatively high wind speeds where mixed attached and separated flow formed on the suction surface of the blade. It seems that the discrepancy is related to the onset of stall phenomena and consequently separation prediction capability of the current turbulence model. It is also found that strong span-wise flow occurs in stalled area due to the centrifugal force generated by rotation of the turbine rotor and it prevents abrupt reduction of normal force for higher wind speed than the designed value.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.104-111
• /
• 2009

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, the boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Using the newly modified code, numerical simulations were performed and compared with other computational results as well as the experimental data for the supersonic flows over an oblique shock with a bleed region.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.5
• /
• pp.340-345
• /
• 2017

It is important to predict the hydrodynamic maneuvering derivatives, which consist of the forces and moment acting on a hull during a maneuvering motion, when estimating the maneuverability of a ship. The estimation of the maneuverability of a ship with a change in the stern hull form is often performed at the initial design stage. In this situation, a method that can reflect the change in the hull form is necessary in the prediction of the maneuverability of the ship. In particular, the linear hydrodynamics maneuvering derivatives affect the yaw checking motion as the key factors. In the present study, static drift calculations were performed using Computational Fluid Dynamics (CFD) based on Reynolds Average Navier-Stokes (RANS) for a 40-segment hull. A prediction method for the linear hydrodynamic maneuvering derivatives was proposed using the slender body theory from the distribution of the lateral force acting on each segment of the hull. Moreover, the results of a comparison study to the model experiment for KVLCC1 performed by KRISO are presented in order to verify the accuracy of the static drift calculation. Finally, the linear hydrodynamic maneuvering derivatives obtained from both the model test and calculation are compared and presented to verity the usefulness of the method proposed in this study.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.1
• /
• pp.66-76
• /
• 2011

The heat transfer and pressure loss characteristics of a rotating two-pass channel with a guide vane in the turning region have been studied using three-dimensional Reynolds-averaged Navier-Stokes (RANS) analysis, and the shape of the guide vane has been optimized using surrogate modeling optimization technique. For the optimization, thickness, location and angle of the guide vanes have been selected as design variables. The objective function has been defined as a linear combination of the heat transfer and the friction loss related terms with a weighting factor. Latin hypercube sampling has been applied to determine the design points as design of experiments. A weighted-average surrogate model, PBA has been used as the surrogate model. The guide vane in the turning region does not influence the heat transfer in the first passage upstream of the turning region, but enhances largely the heat transfer in the turning region and the second passage. In an example of the optimization, the objective function has been increased by 13.6%.