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

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.4
• /
• pp.329-335
• /
• 2009

The unsteady panel and time-marching free wake are applied to the rotor aerodynamics and wake behaviour. Numerical results of panel and free wake are compared and validated with experimental data. Using these methods, unsteady rotor aerodynamics in BVI condition is analyzed and discussed in detail.

• International Journal of Aeronautical and Space Sciences
• /
• v.8 no.1
• /
• pp.1-9
• /
• 2007

As a first step toward a complete CFD-CSD coupling for helicopter rotor load analysis, the present study attempts to loosely couple a CFD code with a source-double panel method. The far-field wake effects were calculated by a time-marching free vortex wake method and were implemented into the CFD module via field velocity approach. Unlike the lifting line method, the air loads correction process is not trivial for the source-doublet panel method. The air loads correction process between the source-doublet method and CFD is newly suggested in this work and the computation results are validated against available data for well-known hovering flight conditions.

• Journal of Ocean Engineering and Technology
• /
• v.8 no.2
• /
• pp.31-37
• /
• 1994

3-Dimensional panel method is now developed to the level that one can calculate the lift of a three dimensional body with the same accuracy of wind tunnel test and some current codes can consider the boundary layer effects due to the viscosity and unsteady motion in the calculation of lift. This paper is also aimed to develop these kinds of computing programs, and as a beginning, the authors restricted the problems to the steady potential flow cases. The calculation of 3-Dimensional body, wing and tandem wing carried out, using source panel and vortex ring panel. Finally, the interactions between 3-Dimension symmetric body and a wing are also calculated.

• Journal of Navigation and Port Research
• /
• v.37 no.1
• /
• pp.1-7
• /
• 2013

In this paper the developed prediction technique of wave-making resistance performance for a ship attached with a vertical blade had been verified. Numerical analysis program as a prediction technique had been developed using the Rankine source panel method and the vortex lattice method(VLM). The nonlinearity of the free surface conditions was fully taken into account using the iterative method and the trim and the sinkage of the ship were also considered in the numerical analysis program. Panel cutting method was applied to get hull surface panels. Numerical computations were carried out for a 4000TEU container carrier and the vertical blade was attached 6 different locations astern. To investigate the validity of the numerical analysis program the commercial viscous flow field analysis program FLUENT was used to obtain the viscous flow field around the ship and the model test was performed. The model test results were compared with the numerical analysis results.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.4
• /
• pp.403-410
• /
• 2008

This paper describes the optimization of the rudder section by the genetic algorism based on VLM(Vortex Lattice Method) and panel method. The developed propeller-rudder analysis program has been validated by comparing with experimental data. The research extends to optimize the anti-erosion rudder section of the large container ship. The object function is the amount of pressure at leading edge of rudder which is closely related with erosion phenomena. The optimized rudder has been compared with conventional rudder with NACA 0021 section by analyzing with the developed program. The finally optimized section has low and mild pressure distribution in comparison with the NACA rudder. The experiments is expected to be carried out for the validation of the present optimization and more parametric study of section geometry is also expected to be conducted in the near future.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.3
• /
• pp.189-195
• /
• 2017

A three-dimensional higher order boundary element method based on the B-spline is presented. The method accurately models piecewise continuous bodies and induced velocity potentials using B-spline tensor product representations, and it is capable of obtaining accurate pointwise values for the potential and its derivatives, especially in the trailing edge and tip region of the lift generating body, which may be difficult or impossible to evaluate with constant panel methods. In addition, we implement a wake roll-up and examine the tip vortex formation in the near wake region. The results are compared with existing numerical results and the results of experiments performed out at the cavitation tunnel of Chungnam National University.

• Journal of Institute of Convergence Technology
• /
• v.5 no.1
• /
• pp.37-40
• /
• 2015

Numerical analyses on the aerodynamic characteristics of a counter rotating axial flow fan is carried out using the frequency domain panel method. Front rotor and rear rotor blades of a counter rotating axial fan are designed by using the simplified meridional flow analysis method with the radial equilibrium equation and the free vortex design condition, according to design requirements. Performance characteristics of a counter rotating axial flow fan are estimated for the variation of design parameters such as the hub to tip ratio, the taper ratio and the solidity. Pressure losses were higher at leading edge and hub region of rotor blades. Characteristic curve of the counter rotating fan was overpredicted without consideration of viscous effect.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.10-17
• /
• 2002

A study on the dynamic ground effect on three-dimensional wings is done using an indirect boundary element method(unsteady panel method). An integral equation is obtained by applying Green's theorem on all surfaces of the fluid domain. Constant strength dipole and source panels arc distributed on a wing's surface. The wake sheet is represented by constant strength dipoles. At each time step, a row of wake panels is assumed to be convected from the trailing edge of the wing. The tip vortex behind wings in dynamic ground effect moves outward. The amplitudes of the aerodynamic coefficients for the wings in dynamic ground effect are augmented much more comparing to the case in static ground effect.