• Journal of the Korea Institute of Military Science and Technology
• /
• v.17 no.1
• /
• pp.25-32
• /
• 2014

Recently supercavitating torpedo has been studied because of its high-speed performance as the next generation of underwater weapon systems. In this study, we present a numerical method based on the potential flow. Characteristic features of the shape of supercavities and drag forces are investigated. In addition, we introduce a viscous-potential method to compensate for the effects of viscosity. The results are compared with viscous calculations using a commercial package, FLUENT V13.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.127-138
• /
• 1999

In this study steady state solutions of cavity flows driven by two moving walls are studied. The north and east walls of the cavity are movable where as the remaining two walls are fixed in space. Numerical experiments for three different driving schemes for moving walls are done at two different Reynolds numbers of Re=40 and 400. The first scheme is to accelerate north and east walls simultaneously. In the second one, the north wall is started first and the east wall is accelerated later. In the third one the east wall starts first. It is usually expected that all these three cases yield the same steady state solution after sufficiently long time. However, present numerical experiments show that such a usual belief is valid only when the Reynolds number is low enough (Re=40). At higher Reynolds number (Re=400), the flow develops to three different steady states depending on the history of the boundary condition change.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.3
• /
• pp.41-50
• /
• 1989

A basic formulation of the panel method, which is based on the potential field formulation, is reviewed for the case of two-dimensional hydrofoil problems. Numerical procedures to improve the computational efficiency of the panel method are suggested. By investigating local behavior of the flow around the trailing edge, a wedge type Kutta condition is formulated. By subdividing the trailing edge panels, where dipole strengths of the subdivided panels follow the local behavior of the potential values of the flow outside a wedge, the circulation around a hydrofoil is calculated accurately with a relatively small number of panels. The subdividing technique to improve the accuracy of the numerical Kutta condition is proved to be efficient. A local behavior of the flow around the leading edge is also investigated. By matching the flow around the leading edge with that around a parabola, a very accurate velocity distribution is obtained with relatively small number of panels. An accurate prediction of the stagnation point and the pressure distribution near the leading edge may contribute to improve the accuracy of cavity predictions and boundary layer calculations around hydrofoils.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.35-42
• /
• 2007

This paper describes a recent development on the divergence free basis function based on a hermite stream function and verifies its validity by comparing results with those from a modified residual method known as one of stabilized finite element methods. It can be shown that a proper choice of degrees of freedom at a node with a proper arrangement of the hermite interpolation functions can yield solenoidal or divergent free interpolation functions for the velocities. The well-known cavity problem has been chosen for validity of the present algorithm. The comparisons from numerical results between the present and the modified residual showed the present method yields better results in both the velocity and the pressure within modest Reynolds numbers(Re = 1,000).

• Journal of Ocean Engineering and Technology
• /
• v.25 no.5
• /
• pp.9-14
• /
• 2011

The divergence-free finite elements introduced in this paper are derived from Hermite functions, which interpolate stream functions. Velocity bases are derived from the curl of the Hermite functions. These velocity basis functions constitute a solenoidal function space, and the gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into its solenoidal and irrotational parts, and the decoupled Navier-Stokes equations are then projected onto their corresponding spaces to form appropriate variational formulations. The degrees of the Hermite functions we introduce in this paper are bi-cubis, quartic, and quintic. To verify the accuracy and convergence of the present method, three well-known benchmark problems are chosen. These are lid-driven cavity flow, flow over a backward facing step, and buoyancy-driven flow within a square enclosure. The numerical results show good agreement with the previously published results in all cases.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.69-76
• /
• 1991

A design method based on the surface vorticity distribution is developed to generate a two-dimensional blade section for prescribed velocity distribution in potential flow. The boundary condition used to determine the strength of vorticity distribution requires that the surface of blade section should be a streamline of the resulting flow. In order to obtain the required final geometry of a two-dimensional blade section, an iterative procedure is used. A computer program is developed and several numerical results are presented.

• Computational Structural Engineering
• /
• v.4 no.4
• /
• pp.135-144
• /
• 1991

The finite element technique incorporatating infinite elements is applied to analyzing the general three dimensional wave-structure interaction problems within the limits of linear wave theory. The hydrodynamic forces are assumed to be inertially dominated, and viscous effects are neglected. In order to analyze the corresponding boundary value problems efficiently, two types of elements are developed. One is the infinite element for modeling the radiation condition at infinity, and the other is the fictitious bottom boundary element for the case of deep water. To validate those elements, numerical analyses are performed for several floating structures. Comparisons with the results by using other available solution methods show that the present method incorporating the infinite and the fictitious bottom boundary elements gives good results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.12
• /
• pp.931-939
• /
• 2013

In this paper, we try to find the lifting-line method which is applicable to the conceptual design of aircraft wings, and analyze the accuracy and coverage of the method. Two methods that are extended from the lifting-line theory of Prandtl are selected. One of the methods is Weissinger's method which imposes the velocity boundary condition at the control points located at the quarter chord, and the other is Phillips's method which combines the three-dimensional vortex lifting law. Calculations are performed for an elliptic wing, a swept back wing, and a tapered unswept wing with dihedral angle and geometric twist. The aerodynamic data of the potential flow such as spanwise distributions of circulation and downwash, lift and induced drag are obtained through calculations, and these data are compared with theoretical results and wind tunnel test data. As a result, Weissinger's method showed good accuracy and reliability regardless of wing shapes, but Phillips's method revealed inaccurate results for a swept back wing.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.6 s.150
• /
• pp.638-646
• /
• 2006

In this paper a new hull-panel generation algorithm named 'Panel Cutting Method' was developed to solve the flow phenomena around a ship advancing on the free surface with a constant speed. In this algorithm the non-linearity of the free surface boundary conditions was taken into account using the iterative method and the raised panel was used at each iteration step. Numerical calculations were performed to investigate the validity of the developed algorithm using the series $60(C_B=0.60)$ hull The wave resistance coefficients, the wave patterns and the wave heights were compared between the computed and the experimental results at Fn=0.25 and 0.316. The comparison showed good agreement between computation and experiment.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].