• The KSFM Journal of Fluid Machinery
• /
• v.2 no.1 s.2
• /
• pp.103-107
• /
• 1999

Transient flow inside a hollow shaft of a Gas Turbine engine during sudden engine stop may result in non uniform heat transfer coefficients in the shaft due to flow instability similar to steady Taylor vortex, which may decrease the lifetime of the shaft. In the present study, transient Taylor vortex phenomena inside a suddenly stopped hollow shaft are studied analytically. Flow visualization is also performed to study the shape and onset time of Taylor Vortices for various initial rotational speed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.10a
• /
• pp.130-135
• /
• 1995

본 연구에서는 30000rpm으로 회전하는 진공청소기 원심팬의 유동장을 임펠러, 디퓨저, 케이싱을 모두 고려하여 해석하였다. 또한 삼차원으로 배출되는 출구를 간단한 sink 패널로 모델하여 출구의 효과를 충분히 고려하였다. 해석된 유동장 자료를 이용하여 먼 거리에서의 음압을 예측하였다. 예측된 음압자료는 FFT를 이용하여 측정된 값과 주파수 영역에서 비교하였다. 또한 진공청소기 원심팬의 측정자료에서 보이는 광역소음특성이 임펠러에서 흘려지는 후류와류의 교란에 의한 임펠러와 디퓨저 깃의 비정상 힘이 주된 원인임을 밝혔다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.1
• /
• pp.9-18
• /
• 2016

The slender structures such as high rise building or marine riser are highly susceptible to dynamic force exerted by fluid-structure interactions among which vortex-induced vibration(VIV) is the main cause of dynamic unstability of the structural system. If VIV occurs in natural frequency regime of the structure, fatigue failure likely happens by so-called lock-in phenomenon. This study presents the numerical analysis of dynamic behavior of both structure and fluid in the lock-in regimes and investigates the subjacent phenomena to hold the resonance frequency in spite of the change of flow condition. Unsteady and laminar flow was considered for a two-dimensional circular cylinder which was assumed to move freely in 1 degree of freedom in the direction orthogonal to the uniform inflow. Fluid-structure interaction was implemented by solving both unsteady flow and dynamic motion of the structure sequentially in each time step where the fluid domain was remeshed considering the movement of the body. The results show reasonable agreements with previous studies and reveal characteristic features of the lock-in phenomena. Not only the lift force but also drag force are drastically increasing during the lock-in regime, the vertical displacement of the cylinder reaches up to 20% of the diameter of the cylinder. The correlation analysis between lift and vertical displacement clearly show the dramatic change of the phase difference from in-phase to out-of-phase when the cylinder experiences lock-in. From the results, it can be postulated that the change of phase difference and flow condition is responsible for the resonating behavior of the structure during lock-in.

• Journal of Advanced Marine Engineering and Technology
• /
• v.20 no.5
• /
• pp.22-29
• /
• 1996

Flow characteristics of two-dimensional closed square cavities near unsteady critical Reynolds numbers were studied numerically at four Reynolds numbers : $8{\times}10^3,\;8.5{\times}10^3,\;9{\times}10^3\;and\;9.5{\times}10^3.$ A convection conservative difference scheme based upon SOLA to maintain the nearly 2nd-order spatial accuracy is adopted on irregular grid formation. Irregular grid number is $80{\times}80$ and its minimum size is about 1/400 of the cavity height(H) and its maximum is about 1/53 H. The result shows that the critical Reynolds number indicating the emergence of flow wnsteadiness is ranging from Re=$8{\times}10^3\;to\;8.5{\times}10^3$ and their flow patterns reveal periodic fluctuation during transient and fully developed stages. But macroscopic flow behavior in terms of instantaneous and time-mean characteristics represent remarkable difference.

• Aerospace Engineering and Technology
• /
• v.2 no.2
• /
• pp.142-150
• /
• 2003

Turbulent flows and related heat transfer in a square heated duct is investigated by a turbulence model and a large eddy simulation. The cooling channel of calorimeter is modeled to the square duct. The nonlinear k-ε-fμ model of Park et al. [3] is slightly modified and their explicit heat flux model is employed. The Reynolds number is varied in the range 4000≤Reｂ≤20000. The heat transfer is closely linked to the secondary flows which driven by the turbulent motion. Its magnitude is 1~3% of the mean streamwise velocity. The relation of Nu~Re0.8Pr0.34 is validated by comparing with the predicted Nu of k-ε-fμ model. Finally, the coherent structures and thermal fluctuations are scrutinized.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.11a
• /
• pp.305-310
• /
• 2007

For the understanding of the complex flow characteristics in the counter-rotating axial flow fan, it is necessary to investigate the three-dimensional unsteady flow fields in the counter-rotating axial flow fan. This information is also essential for the prediction of the aerodynamic and acoustical characteristics of the counter-rotating axial flow fan. Experimental study on the three-dimensional unsteady flow in the counter-rotating axial flow fan is carried out at the design point(operating condition). Three-dimensional unsteady flow fields in the counter rotating axial flow fan are measured at the cross-sectional planes of the upstream and downstream of each rotor using the $45^{\circ}$ inclined hot-wire. Three-dimensional unsteady flow fields in the counter-rotating axial flow such as the wake, the tip vortex and the tip leakage flow are shown the form of the velocity vectors and the velocity contours.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.10
• /
• pp.859-867
• /
• 2007

The objective of this study is to develop and validate a numerical method which can handle the multi-rotor aerodynamic characteristics. For the purpose of power estimation, table look-up method is implemented to the existing unsteady panel code that is coupled with a time-marching free wake model. Also, the Reynolds number scaling is implemented for the application to various regions of Reynolds number. The computed results are validated against the available experimental data for coaxial and tandem rotors. In the validation case for the coaxial rotor, more accurate result is acquired when the thickness effect is considered. The wake instability problem occurs at a particular separation distance between the rotors for tandem rotors. The wake instability is avoided by setting the single-rotor wake geometry as the initial wake geometry for the multi-rotor analysis. The estimated result for rotor separation effect is compared with the result of the momentum theory.

• Journal of Advanced Marine Engineering and Technology
• /
• v.31 no.2
• /
• pp.152-158
• /
• 2007

Most of experimental visualizations and numerical results on the flow field separated form a leading edge around an unsteady foil show a continuous streakline from the leading edge and large reverse flow between the streakline and the suction surface. However, they have not exactly clarified yet the dynamic behavior of vortices separated from the leading edge because separation around an unsteady foil is very complicated phenomenon due to many parameters. In the present study the flow fields around pitching foils have been visualized by using a Schlieren method with a high speed camera in a wind tunnel at low Reynolds number regions. It has been observed that small vortices are shed discretely from the leading and trailing edge and that they stand in line on the integrated streakline of separation shear layer. By counting vortices in the VTR frames it was clarified that the number of vortex shedding from the leading and trailing edge during one pitching cycle strongly depends on the non-dimensional pitching rate. Futhermore the vortices moving up to the leading edge on the suction surface of the pitching foil are visualized. They play an important role to balance the number of vortex shedding from both edges.

• Korean Chemical Engineering Research
• /
• v.58 no.2
• /
• pp.319-324
• /
• 2020

Electroconvective instability is a non-linear transport phenomenon which can be found in ion-selective transport system such as electrodialysis, Galvanic cell and electrolytic cell. The instability is triggered by the fluctuation of space charge layer in adjacent of ion-selective surface, leading to increase of mass transport rate. Thus, in the aspect of mass transport, the instability has an important meaning. Although recent experimental techniques have opened up an avenue to direct visualize the instability, fundamental investigations have been conducted in limited area due to several experimental limitations. In this work, the electroconvective instability under potentiostatic mode was solved by numerical method in order to demonstrate correlation between current-time curve and the instability behavior. By rigorous time-resolved analysis, the transition behaviors can be divided into three stages; formation of space charge layer - growth of electroconvective instability - steady state. Furthermore, scaling laws of transition time were numerically obtained according to applied voltage as well.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.75-75
• /
• 2011

본 논문에서는 진동중인 교량 거더에 작용하는 풍하중을 산정하고 그에 따른 플러터 발생풍속을 예측하기 위하여 분산형 전산환경을 활용한 수치해석 연구를 수행하였다. 분산형 전산환경은 웹 포탈을 기반으로 수치해석 환경을 제공하는 수치풍동 시스템으로서, 전산유체역학(CFD : Computational Fluid Dynamics)에 대한 전문지식이 부족한 사용자들도 격자생성, 수치해석자를 이용한 계산, 가시화 등의 전 과정을 편리하게 수행할 수 있는 차세대 토목분야 연구 환경이다. 본 시스템은 그리드스피어(GfidSphere)를 기반으로 구성되었으며, 기본적으로 사용자 관리, 세션 관리, 그룹 관리, 레이아웃 관리 등을 제공하여 사용자가 포탈을 통해서 다양한 서비스를 쉽게 사용할 수 있는 환경을 구축하도록 도와준다. 수치해석을 위한 유체 지배방정식은 2차원 비정상 비압축성 RANS(Reynolds-Averaged Navier-Stokes) 방정식이며, pseudo compressibility 방법을 적용하였다. 비정상 유동장을 해석하기 위하여 이중시간 전진법(dual time stepping)을 사용하였으며, 수렴가속화를 위해 Multi-grid 기법을 적용하였다. 또한 난류 유동장 해석을 위해서 $k-{\omega}$ SST 난류 모델을 사용하였으며, 난류 천이 과정에서의 유동을 모사하기 위하여 Total stress limitation 방법을 적용하였다. 교량 거더의 연직과 회전방향의 2자유도 움직임을 모사하기 위하여 동적격자 기법을 도입하였다. 교량 거더 주변의 비정상 유동해석 결과를 통해, 거더 표면에서 떨어져나가는 크고 작은 와류의 영향으로 양력 및 모멘트 계수 그래프가 중첩된 진폭과 주기를 갖고 주기적으로 나타나는 것을 확인할 수 있었다. 또한 계산된 비정상 공기력을 적용한 2자유도 플러터 방정식을 통하여 플러터 발생풍속을 산정하였다. 최종적으로 본 연구에서 계산된 결과의 타당성을 검증하기 위하여 수치적으로 구한 플러터 발생풍속과 기존의 실험 및 수치해석 결과를 비교하였으며, 결과는 잘 일치하였다.