• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.1
• /
• pp.26-32
• /
• 2012

The 3-D unsteady viscous flow around an impulsively started rotating marine propeller is simulated using VIC(Vortex-In-Cell) method which is adequate to analyze the strong vortical flow around complicatedly-shaped body. The computational procedure is governed by the vorticity transport equation in Lagrangian form. In order to solve the equation, a regular grid which is independent to the shape of a body is introduced and each term of the equation is evaluated numerically on the grid by applying immersed boundary concept. In this paper, the overall algorithm including the formulation of governing equations and boundary conditions is described and some computational results are presented with discussing their physical validity.

• Ocean Systems Engineering
• /
• v.3 no.3
• /
• pp.219-236
• /
• 2013

The accurate prediction of motion in waves of a marine vehicle is essential to assess the maximum sea state vs. operational requirements. This is particularly true for small crafts, such as Autonomous Surface Vessels (ASV). Two different numerical methods to predict motions of a SWATH-ASV are considered: an inviscid strip theory initially developed at MIT for catamarans and then adapted for SWATHs and new a hybrid strip theory, based on the numerical solution of the radiation forces by an unsteady viscous, non-linear free surface flow solver. Motion predictions obtained by the viscous flow method are critically discussed against those obtained by potential flow strip theory. Effects of viscosity are analyzed by comparison of sectional added mass and damping calculated at different frequencies and for different sections, RAOs and motions response in irregular waves at zero speed. Some relevant conclusions can be drawn from this study: influence of viscosity is definitely non negligible for SWATH vessels like the one presented: amplitude of the pitch and heave motions predicted at the resonance frequency differ of 20% respectively and 50%; in this respect, the hybrid method with fully non-linear, viscous free surface calculation of the radiation forces turns out to be a very valuable tool to improve the accuracy of traditional strip theories, without the burden of long computational times requested by fully viscous time domain three dimensional simulations.

• Journal of the Korean Institute of Navigation
• /
• v.19 no.2
• /
• pp.13-22
• /
• 1995

This paper presents the numerical analysis results of the viscous flow for a full ship model. The mass and momentum conservation equations are used for governing equations, and the flow field is discretized by the Finite-Volume Method for the numerical calculation. An algebraic grid and elliptic grid generation techniques are adopted for generation of the body-fitted coordinates system, which is suitable to ship's hull forms. Time-marching procedure is used to solve the three-dimensional unsteady problem, where the convection terms are approximated by the QUICK scheme and the 2nd-order central differencing scheme is used for other spatial derivatives. A Sub-Grid Scale turbulence model is used to approximate the turbulence, and the wall function is used at the body surface. Pressure and velocity fields are calculated by the simultaneous iteration method. Numerical calculations were accomplished for the Crude Oil Tanker(DWT 95,000tons, Cb=0.805) model. Calculation results are compared to the experimental results and show good agreements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.80-86
• /
• 2005

This study presents numerical solutions of the two-dimensional Navier-Stokes equations for supersonic unsteady flow in a convergent-divergent nozzle with heat addition. The TVD scheme in generalized coordinates is employed in order to calculate the moving shock waves caused by thermal choking. We discuss on transient characteristics, start and unstart phenomena, fluctuations of specific thrust caused by thermal choking and viscous effects. We prove that the control of separation of boundary layer is the most important key problem to prevent the thermal choking.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1601-1606
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.809-817
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.42-49
• /
• 2008

A wavelet method is presented in order to improve the computational efficiency of two dimensional unsteady flow problems while maintaining the order of accuracy of conventional CFD schemes. First, by using the interpolating wavelet transformation including decomposition and thresholding, an adaptive dataset to a solution is constructed. Then, inviscid and viscous fluxes are calculated only at the points within an adaptive dataset, which enhances the computational efficiency. Second, thresholding step is modified to maintain the spatial and temporal accuracy of conventional CFD schemes automatically by selecting the threshold value between user-defined value and the magnitude of spatial or temporal truncation error. The wavelet method suggested in this study is successfully applied to various unsteady flow problems and it is shown that the computational efficiency is enhanced with maintaining the computational accuracy of CFD schemes.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.42-49
• /
• 2008

A wavelet method is presented in order to improve the computational efficiency of two dimensional unsteady flow problems while maintaining the order of accuracy of conventional CFD schemes. First, by using the interpolating wavelet transformation including decomposition and thresholding, an adaptive dataset to a solution is constructed. Then, inviscid and viscous fluxes are calculated only at the points within an adaptive dataset, which enhances the computational efficiency. Second, thresholding step is modified to maintain the spatial and temporal accuracy of conventional CFD schemes automatically by selecting the threshold value between user-defined value and the magnitude of spatial or temporal truncation error. The wavelet method suggested in this study is successfully applied to various unsteady flow problems and it is shown that the computational efficiency is enhanced with maintaining the computational accuracy of CFD schemes.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.6
• /
• pp.15-20
• /
• 2005

Two-dimensional incompressible Navier-Stokes flow solver is developed using SIMPLER method to study the unsteady viscous flow physics over two-dimensional ellipses. Unsteady viscous flows over various thickness-to-chord ratios of 0.6, 0.8, 1.0, and 1.2 elliptic cylinders are simulated at different Reynolds numbers of 200, 400, and 1,000. This study is focused on the understanding the effects of Reynolds number and elliptic cylinder thickness on the drag and lift forces. The present numerical solutions are compared with available experimental and numerical results and show a good agreement. Through this study, it is observed that the Reynolds number and the cylinder thickness affect not only the frequency of the force oscillations but also the mean values and the amplitudes of the total drag and lift forces significantly.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.189-192
• /
• 2002

In this paper, the submarine model, called DARPA SUBOFF model, has been numerically analyzed to investigate the aerodynamic forces variation in terms of angle of attacks and yaw angles. The SUBOFF model is consisted of the three parts : axisymmetric body, fairwater, and four symmetric stern appendages. Three dimensional unsteady incompressible Wavier-Stokes equation was used on curvilinear multi-block grid system. To validate the present code, the SUBOFF tare hull and an ellipsoid at angle of attacks of $10^{\circ}\;and\;30^{\circ}$ were simulated and a good agreement with experiments was obtained. After the code validation, the flows over SUBOFF model were simulated at three different angle of attacks and yaw angles. The variation of aerodynamic forces in terms of angle of attack and yaw angle were calculated. Also, to understand the flow features around a submarine with variation of yaw and attack angle, the pressure contours and streamlines were plotted.