• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.5
• /
• pp.400-409
• /
• 2016

In recent, shipping companies have made an enormous effort to improve the operation of vessel in various approaches, due to recession of shipping market and increasing competition among shipping companies. One of important parameters for improving the efficiency of vessel is the resistance performance that consist of friction and residual resistance. Especially, it is recognized that the friction resistance tends to be affected by conditions of vessel’s surface and occupies approximately 70~90% of the total resistance for slow speed ships. In general, the surface of vessel is covered with various type of paint to reduce fouling and corrosion. As time goes by, however, it is so hull roughness would be increased by fouling over the wetted surface that anti-fouling paints, such as CDP(Controlled Depletion Paint), Tin-Free SPC(Self Polishing Co-polymer) or Foul Release, are applied evenly on the hull surface. Nevertheless, these anti-fouling paints could not prevent fouling absolutely. A fundamental study on evaluating ship resistance performance variation due to hull roughness has been performed using a commercial software, Star-CCM+, which solves the continuity and Navier-Stokes equations for incompressible and viscous flow. The results of present simulation for plate are compared with some experimental data available and the effect of surface roughness to ship resistance performance is discussed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.95-104
• /
• 2003

This paper deals with the analytic and FEM analyses of sloshing frequency response of incompressible, invicid and irrotational flow in two dimensional rectangular tank. We use Laplace equation based on potential theory as governing equation. For small amplitude sloshing motion, the linearized free surface condition was applied and the analytic solution as obtained by the separation of variables. To simulate the effect of the energy dissipation due to viscous damping, artificial viscous coefficient is introduced and the divergence of response at resonance frequencies may be avoided by this coefficient. This problem was solved by FEM using 9-node elements in order to predict the maximum amplitude of sloshing response. Numerical results of free surface height, fluid pressure and fluid force show good agreement with those by analytic solution. After verifying the test FEM program, we analyze the frequency response characteristics of sloshing to the fluid height.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.33-42
• /
• 2003

This paper deals with the FEM analysis of nonlinear sloshing of incompressible, invicid and irrotational flow in two dimensional rectangular tank. We use laplace equation based on potential theory as governing equation. For large amplitude sloshing motion, kinematic and dynamic free surface conditions derived from Bernoulli equation are applied. This problem is solved by FEM using 9-node elements. For the time integration and accurate velocity calculation, we introduce predictor-corrector time marching scheme and least square method. Also, numerical stability in tracking of free surface is obtained by direct calculation of free surface location to time variation. Numerical results of sloshing induced by harmonic excitations, while comparing with those of linear theory and references, prove the accuracy and stability. After verification of our program, we analyze sloshing response characteristics to the fluid height and the excitation amplitude.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.1-7
• /
• 2014

The present study numerically investigated the effect of the advance ratio on the wake characteristics of the marine propeller in the propeller open water test. Therefore, a wide range of the advance ratio(0.2${\kappa}-{\omega}$SST Model are considered. The three-dimensional vortical structures of tip vortices are visualized by the swirl strength, resulting in fast decay of the tip vortices with increasing the advance ratio. Furthermore, to better understanding of the wake evolution, the contraction ratio of the slip stream for different advance ratios is extracted from the velocity fields. Consequently, the slip stream contraction ratio decreases with increasing the advance ratio and successively the difference of the slip stream contraction ratio between J=0.2 and J=0.8 is about 0.1R.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.245-250
• /
• 2002

This report describes the effects of following flaws due to ship's propeller on the fish farm structure when the ship's propeller is operated in full power. This study is applied an incompressible newtonian fluid theory, which is governed the Navier-Stokes equation. For the numerical solution, Neumann equation are applied as the boundary conditions. The result shows that the flow velocity near the fish farm is 1.0 m/sec. The actual measurement carries out by using propeller type velocimeter in order to measure the velocity of following flows and currents around the fish farm area. The result shows that the maximum velocity near the fish farm structure is 1.2 m/sec in depth of 1.5 m. This velocity is used for calculation of external force on the fish farm structure. The results of structural strength of the fish farm structures show that the actual maximum bending moment and bending stress are less than the damage strength of material. So the fish farm structure is not affected by the following flows and currents of ship's propeller.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.2
• /
• pp.151-161
• /
• 2012

The present study numerically investigates the aerodynamic characteristics of two-dimensional wings in the vicinity of the ground by solving two-dimensional steady incompressible Navier-Stokes equations with the turbulence closure model of the realizable k-${\varepsilon}$ model. Numerical simulations are performed at a wide range of the normalized ground clearance by the chord length ($0.1{\leq}h/C{\leq}1.25$) for the angles of attack ($0^{\circ}{\leq}{\alpha}{\leq}10^{\circ}$) in the prestall regime at a Reynolds number (Re) of $2{\times}10^6$ based on free stream velocity $U_{\infty}$ and the chord length. As the physical model of this study, a cambered airfoil of NACA 4406 has been selected by a performance test for various airfoils. The maximum lift-to-drag ratio is achieved at ${\alpha}=4^{\circ}$ and h / C = 0.1. Under the conditions of ${\alpha}=4^{\circ}$ and h / C = 0.1, the effect of the Reynolds number on the aerodynamic characteristics of NACA 4406 is investigated in the range of $2{\times}10^5{\leq}Re{\leq}2{\times}10^9$. As Re increases, $C_l$ and $C_d$ augments and decreases, respectively, and the lift-to-drag ratio increases linearly.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.127-136
• /
• 2022

In this work a multi-fidelity non-intrusive polynomial chaos (MF-NIPC) has been applied to a structural wind engineering problem in architectural design for the first time. In architectural design it is important to design structures that are safe in a range of wind directions and speeds. For this reason, the computational models used to design buildings and bridges must account for the uncertainties associated with the interaction between the structure and wind. In order to use the numerical simulations for the design, the numerical models must be validated by experi-mental data, and uncertainties contained in the experiments should also be taken into account. Uncertainty Quantifi-cation has been increasingly used for CFD simulations to consider such uncertainties. Typically, CFD simulations are computationally expensive, motivating the increased interest in multi-fidelity methods due to their ability to lev-erage limited data sets of high-fidelity data with evaluations of more computationally inexpensive models. Previous-ly, the multi-fidelity framework has been applied to CFD simulations for the purposes of optimization, rather than for the statistical assessment of candidate design. In this paper MF-NIPC method is applied to flow around a rectan-gular 5:1 cylinder, which has been thoroughly investigated for architectural design. The purpose of UQ is validation of numerical simulation results with experimental data, therefore the radius of curvature of the rectangular cylinder corners and the angle of attack are considered to be random variables, which are known to contain uncertainties when wind tunnel tests are carried out. Computational Fluid Dynamics (CFD) simulations are solved by a solver that employs the Finite Element Method (FEM) for two turbulence modeling approaches of the incompressible Navier-Stokes equations: Unsteady Reynolds Averaged Navier Stokes (URANS) and the Large Eddy simulation (LES). The results of the uncertainty analysis with CFD are compared to experimental data in terms of time-averaged pressure coefficients and bulk parameters. In addition, the accuracy and efficiency of the multi-fidelity framework is demonstrated through a comparison with the results of the high-fidelity model.

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.192-198
• /
• 2004

A numerical method is developed to solve a two-dimensional diffraction problem for a body located in a sediment pocket where a heavier muddy water is trapped. In the present study, the wave exciting forces acting on a submerged body in the water-sediment interface by an incident wave is investigate. It is assumed that the heavier mud is trapped locally in a sediment pocket. A mathematical formulation is made in the scope of the potential theory. The fluid is assumed to be inviscid, incompressible and its motion irrotational. The boundary conditions on the unknown free surface and interface are linearized. As a method of solution, the localized finite-element method is adopted. In the method, the computation domain is reduced by utilizing the complete set of analytic solutions known in the infinite subdomain to be truncated by introduction of an appropriate juncture conditions. The main advantage of this method is that any complex geometry of the boundaries can be easily accommodated. Computations are carried out for mono-chromatic plane progressive surface waves normally incident on the domain. Numerical results are compared with those obtained by Lassiter based on Schwingers variational method. Good Agreements are obtained in general. Another numerical computations are made for the cases with and without a body in the sediment pocket.