• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.71-76
• /
• 2004

An extension of our recently developed locally linear reconstruction scheme to 2 dimensional incompressible flow solver is presented. The solver is based on a semi-implicit fractional step method in which the convective term is discretized by Adams-Bashforth method and the diffusion term by Crank-Nicolson method. Several numerical examples are tested to demonstrate the mesh type independent accuracy of the solver, which include decaying vortex flow, square cavity flow, and flow around a circular cylinder. The above examples are solved on quadrilateral or hybrid meshes. For all numerical examples, we obtained reasonable results.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.9-13
• /
• 2006

In the present study, a numerical analysis of transient mixer flow is performed considering free surface formation. The flow patterns and free surface shape in a mixers formed by flat paddle and pitched paddle impellers are predicted. In a flat paddle mixer, two flow circulation regions are formed due to strong radial flow, whereas one large circulation is formed in a pitched paddle mixer due to axial downward flow. These differences affect the free surface evolution and shape. It is seen from the results that a flat paddle mixer gives deeper free surface at center region than a pitched paddle mixer. The free surface of 8-blades-flat-paddle mixer is also simulated to compare with the available experimental and simulation results. The present computational results agree reasonably well with the experimental data.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.11-16
• /
• 2002

In this paper, we present two successful results from active controls of flows over a circular cylinder and a sphere for drag reduction. The Reynolds number range considered for the flow over a circular cylinder is 40-3900 based on the free-stream velocity and cylinder diameter, whereas for the flow over a sphere it is $10^{5}$ based on the free-stream velocity and sphere diameter. The successful active control methods are a distributed (spatially periodic) forcing and a high-frequency (time periodic) forcing. With these control methods, the mean drag and lift fluctuations decrease and vortical structures are significantly modified. For example, the time-periodic forcing at a high frequency (larger than 20 times the vortex shedding frequency) produces $50{\%}$ drag reduction for the flow over a sphere at $Re=10^{5}$. The distributed forcing applied to the flow over a circular cylinder results in a significant drag reduction at all the Reynolds numbers investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.751-762
• /
• 1995

The snapshot method is introduced to approximate the coherent structures of planar jet flow. The numerical simulation of instantaneous flow field is analyzed by SIMPLE algorithm. An ensemble of realizations is collected using a sampling condition that corresponds to the passage of a large scale vortex at positions 4 and 6 diameters downstream from the nozzle. With snapshot mothod we could treat the data efficiently and approximate coherent structures inhered in the planar jet flow successfully 94% of total turbulent kinetic energy with 10 terms of Karhunen-Loeve expansions. Finally, In accordance with the recent trend to try to explain and model turbulence phenomena with the existence of coherent structures, in the present study, we express the underlying coherent structures of planar jet flow in the minimum number of modes by calculating Karhunen-Loeve expansions in order to improve to understanding of jet flow and to make the information storage and management in computers easier.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.4
• /
• pp.130-140
• /
• 1997

원형주상체 주위의 유동을 규명하기 위해 회류수조에서 원주방향으로 24개의 위치에 대하여 압력을 계측하였으며, laser sheet을 이용하여 유동을 가시화 하였다. Reynolds수가 4800에서 40000인 범위에 대하여 실험을 수행하였다. 또한, 원형단면체 주위의 비정상 층류유동에 대한 Navier-Stokes방정식의 해를 구하는 수치해석기법을 개발하였다. 효과적인 격자배치를 위하여 H와 O-type의 중첩격자를 사용하였고, 이산화 방법으로는 정규격자시스템에서 유한차분법을 적용하였다. 실험과 수치해석결과에서 뚜렷한 와류박리현상을 볼 수 있었으며, 압력계수 (C$_{p}$ ), 항력계수(C$_{D}$), 스트로얼수(St)를 정량적으로 비교하였을 때, 비교적 잘 일치하는 것을 확인하였다.다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.3917-3925
• /
• 1996

A numerical analysis on three dimensional turbulent incompressible flows through linear cascades of turbine rotor blades with high turning angles has been performed by using a generalized k-.epsilon. model which is a high Reynolds number form and derived by RNG(renormalized group) method to account for the variation of the rate of strain. A second order upwind scheme is used to suppress numerical diffusion in approximating the convective terms. Body-fitted coordinates are adopted to represent the complex blade geometry accurately. For the case without tip clearance, velocity vectors and static pressure contours are shown to be in good agreement with previous experimental results. For the case with tip clearance, the effects of the passage vortex and tip clearance flow on the total pressure loss as well as their interactions are discussed.

• Proceeding of EDISON Challenge
• /
• 2015.03a
• /
• pp.563-568
• /
• 2015

풍력터변의 성능을 높이기 위한 방법으로 최근 해외에서는 쉬라우드를 장착하여 유입 유량을 증진시키는 형상에 대한 새로운 아이디어가 제안되고 시범적으로 적용되고 있다. 본 연구에서는 쉬라우드가 장착된 수평축 풍력터빈에 대해서 EDISON CFD를 이용하여 쉬라우드 내부로 유입되는 질량 유량의 변동을 몇 가지 형상에 대해 수치적으로 비교 분석하였다. 유동장은 비압축성 난류유동으로 가정하였으며, 수치 해석 결과로부터 쉬라우드 주변의 순환의 세기를 형상 변동에 따라 도출하였다. 쉬라우드 형상으로는 캠버를 갖는 goe 417 에어포일을 두 개의 받음각(5도, 10도)에 대해서 수치해석을 수행하였으며, 브림을 갖는 디퓨저 형상(Wind-lens)에 대해서도 유입 유량 변동과 순환 세기에 대해 수치해석을 수행하고 결과를 상호 비교하였다. 본 연구를 통해 쉬라우드가 발생시키는 순환에 의한 유입 유량 증가 현상을 파악할 수 있었으며, 이로써 풍력터빈의 출력을 증대시킬 수 있음을 확인하였다.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2568-2573
• /
• 2008

A drop impact on the liquid film/pool generates several phenomena such as the drop floating, bouncing, formation of vortex ring, jetting, bubble entrapment and splashing. These phenomena depend on the impact velocity, the drop size, the drop properties and the liquid film/pool thickness. These parameters can be summarized by four main dimensionless parameters; Weber number, Ohnesorge number, Froude number and non-dimensional film/pool thickness. In the present study, the phenomena of the splashing and bubble entrapment due to the drop impact on the liquid film/pool were numerically investigated by using a Level Set method for the sharp interface tracking of two distinct phases. After the drop impact, the splashing phenomena with the crown formation and spreading were predicted. Under the specific conditions, the bubble entrapment at the base of the collapsing cavity due to the drop impact was also observed. The numerical results were compared to the available experimental data showing good agreements.

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

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