• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• 대한토목학회논문집
• /
• 제6권2호
• /
• pp.93-101
• /
• 1986

정지수역(靜止水域)에서 연직상향(鉛直上向)으로 방류(放流)되는 평면부력(平面浮力)?의 거동을 연속방정식(連續方程式), 운동량수송식(運動量輸送式) 및 추적물(追跡物) 수송식(輸送式)에 의하여 수치적(數値的)으로 해석한다. 이 해석은 유함수(流凾數)와 과수송식(過輸送式)을 도입하고 Prandtl 난류모형(亂流模型)을 이용하였다. 방류밀도(放流密度) Froude 수(數)가 4~32인인 부력(浮力)?의 발달된 흐름영역(領域)에서 구한 ?중심선의 속도(速度)와 온도변화(溫度變化), ?주변수역(周邊水域)의 온도분포(溫度分布)와 흐름양상은 기존자료와 잘 맞음을 확인할 수 있다. 적분형해석(積分型解析)에 필요한 퍼짐율과 확산비(擴散比)가 방류밀도(放流密度) Froude 수(數)와 방류구(放流口)로부터의 거리의 함수로 유도된다.

• 대한기계학회논문집
• /
• 제8권3호
• /
• pp.185-194
• /
• 1984

본 연구에서는 2차원 정상 층류 유동을 지배하는 극좌표계로 표시된 Navier- Stokes 방정식을 유한해석법으로 전개하여 극좌표계에서의 그 수치계산적 성질을 검토 하고 유한해석법의 적용 좌표계 범위를 확장하는 한편, 2차원 부채꼴 캐비티 유동에 응용하여 공학적 해석도구로서의 그 가능성을 확인하고 아직 연구가 불충분한 부채꼴 캐비티 유동을 고찰하는데 그 목적이 있다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제7권3호
• /
• pp.595-609
• /
• 2015

The wake characteristics of Contra-Rotating Propeller (CRP) were investigated using numerical simulation and flow measurement. The numerical simulation was carried out with a commercial CFD code based on a Reynolds Averaged Navier-Stokes (RANS) equations solver, and the flow measurement was performed with Stereoscopic Particle Image Velocimetry (SPIV) system. The simulation results were validated through the comparison with the experiment results measured around the leading edge of rudder to investigate the effect of propeller operation under the conditions without propeller, with forward propeller alone, and with both forward and aft propellers. The evolution of CRP wake was analyzed through velocity and vorticity contours on three transverse planes and one longitudinal plane based on CFD results. The trajectories of propeller tip vortex core in the cases with and without aft propeller were also compared, and larger wake contraction with CRP was confirmed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.1-1
• /
• 2009

In this presentation, I talk about various fluid simulation methods that have been developed for computer graphics special effects since 1996. They are all based on CFD but sacrifice physical reality for visual plausability and time. But as the speed of computer increases rapidly and the capability of GPU (graphics processing unit) improves, methods for more physical realism have been tried. In this talk, I will focus on four aspects of fluid simulation methods for computer graphics: (1) particle level-set methods, (2) particle-based simulation, (3) methods for exact satisfaction of incompressibility constraint, and (4) GPU-based simulation. (1) Particle level-set methods evolve the surface of fluid by means of the zero-level set and a band of massless marker particles on both sides of it. The evolution of the zero-level set captures the surface in an approximate manner and the evolution of marker particles captures the fine details of the surface, and the zero-level set is modified based on the particle positions in each step of evolution. (2) Recently the particle-based Lagrangian approach to fluid simulation gains some popularity, because it automatically respects mass conservation and the difficulty of tracking the surface geometry has been somewhat addressed. (3) Until recently fluid simulation algorithm was dominated by approximate fractional step methods. They split the Navier-Stoke equation into two, so that the first one solves the equation without considering the incompressibility constraint and the second finds the pressure which satisfies the constraint. In this approach, the first step introduces error inevitably, producing numerical diffusion in solution. But recently exact fractional step methods without error have been developed by fluid mechanics scholars), and another method was introduced which satisfies the incompressibility constraint by formulating fluid in terms of vorticity field rather than velocity field (by computer graphics scholars). (4) Finally, I want to mention GPU implementation of fluid simulation, which takes advantage of the fact that discrete fluid equations can be solved in parallel.

• International Journal of Aeronautical and Space Sciences
• /
• 제13권3호
• /
• pp.307-316
• /
• 2012

Numerical simulations of 3D aircraft configurations are performed in order to understand the effects of turbulence models on the prediction of aircraft's aerodynamic characteristics. An in-house CFD code that solves 3D RANS equations and two-equation turbulence model equations are used. The code applies Roe's approximated Riemann solver and an AF-ADI scheme. Van Leer's MUSCL extrapolation with van Albada's limiter is also adopted. Various versions of Menter's $k-{\omega}$ SST turbulence models as well as Coakley's $q-{\omega}$ model are incorporated into the CFD code. Menter's $k-{\omega}$ SST models include the standard model, the 2003 model, the model incorporating the vorticity source term, and the model containing controlled decay. Turbulent flows over a wing are simulated in order to validate the turbulence models contained in the CFD code. The results from these simulations are then compared with computational results from the $3^{rd}$ AIAA CFD Drag Prediction Workshop. Numerical simulations of the DLR-F6 wing-body and wing-body-nacelle-pylon configurations are conducted and compared with computational results of the $2^{nd}$ AIAA CFD Drag Prediction Workshop. Aerodynamic characteristics as well as flow features are scrutinized with respect to the turbulence models. The results obtained from each simulation incorporating Menter's $k-{\omega}$ SST turbulence model variations are compared with one another.

• 한국가시화정보학회:학술대회논문집
• /
• 한국가시화정보학회 2003년도 추계학술대회 논문집
• /
• pp.39-42
• /
• 2003

Leading edge extension(LEX) in a highly swept shape applied to a delta wing features the modern air-fighters. The LEX vortices generated upon the upper surface of the wing at high angle of attack enhance the lift force of the delta wing by way of increased negative suction pressure over the surfaces. The present 3-D stereo PIV includes the Identification of 2-D cross-correlation equation, stereo matching of 2-D velocity vectors of two cameras, accurate calculation of 3-D velocity vectors by homogeneous coordinate system, removal of error vectors by a statistical method followed by a continuity equation criterion and so on. A delta wing model with or without LEX was immersed in a circulating water channel. Two high-resolution, high-speed digital cameras$(1280pixel\times1024pixel)$ were used to allow the time-resolved animation work. The present dynamic stereo PIV represents the complicated vortex behavior, especially, in terms of time-dependent characteristics of the vortices at given measuring sections. Quantities such as three velocity vector components, vorticity and other flow information can be easily visualized via the 3D time-resolved post-processing to make the easy understanding of the LEX effect or vortex emerging and collapse which are important phenomena occurring in the field of delta wing aerodynamics.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1672-1677
• /
• 2004

Leading edge extension(LEX) in a highly swept shape applied to a delta wing features the modem air-fighters. The LEX vortices generated upon the upper surface of the wing at high angle of attack enhance the lift force of the delta wing by way of increased negative suction pressure over the surfaces. The present 3-D stereo PIV includes the Identification of 2-D cross-correlation equation, stereo matching of 2-D velocity vectors of two cameras, accurate calculation of 3-D velocity vectors by homogeneous coordinate system, removal of error vectors by a statistical method followed by a continuity equation criterion and so on. A delta wing model with or without LEX was immersed in a circulating water channel. Two high-resolution, high-speed digital cameras($1280pixel{\times}1024pixel$) were used to allow the time-resolved animation work. The present dynamic stereo PIV represents the complicated vortex behavior, especially, in terms of time-dependent characteristics of the vortices at given measuring sections. Quantities such as three velocity vector components, vorticity and other flow information can be easily visualized via the 3D time-resolved post-processing to make the easy understanding of the LEX effect or vortex emerging and collapse which are important phenomena occurring in the field of delta wing aerodynamics.

• Journal of Advanced Marine Engineering and Technology
• /
• 제36권8호
• /
• pp.1091-1096
• /
• 2012

유한체적법을 기반으로 나비에 스톡스 방정식을 비구조격자로 풀어 실린더 주위의 공력특성을 규명하였다. 보텍스, 속도, 압력, 잔차, 항력계수 등의 데이터를 가지고 분석하였고 레이놀즈 수는 50, 100이다. 유동특성은 Re>50에서 주기적으로 진동하는 소용돌이를 후류에 형성하며 이 현상은 이동하는 실린더에서도 유사한 현상을 보여 주었다. 지면효과는 실린더 위쪽에서 형성된 소용돌이가 벽면에 근접할수록 실린더의 후방으로 길게 늘어나는 형상을 보이고, 실린더와 근접벽면 사이의 유속이 정체되어 실린더와 벽면 사이의 간격이 0.6 에서는 근접 평판과 실린더 사이의 유동이 거의 끊겨짐을 알 수 있었다. 본 수치계산의 검증을 위하여 항력계수를 타 연구결과와 비교하였다.

• 한국음향학회지
• /
• 제21권2호
• /
• pp.134-141
• /
• 2002

본 연구에서 곡선 와동 요소와 CVC (Constant Vorticity Contour) 후류 모델이 수평축 풍력 터빈의 공력 성능 및 소음을 예측하기 위해 사용되었다. 또한 2차 회귀에 근거하여 회전수의 예측을 하였다. 광대역 소음을 예측은 경험식에 근거한 방법을 사용하였다. 직선 와동요소 대신에 BCVE (Basic Curved Vortex Element)와 SIVE(Self Induction Vortex Element)를 사용하는 곡선 와동 요소를 사용하였으며 CVC 후류 모델에서 와동의 강도는 블레이드의 스팬방향 및 와동을 따라서 일정하다고 가정하였다. 이렇게 만들어진 자유 후류는 와동 격자를 대치한다. 본 방법은 여타의 방법에 비해 휠씬 적은 계산 시간을 요구하며 후류의 정확한 구조를 모사할 수 있었다. 검증을 위해서는 김준모의 실험과 Zond사의 Z-40FS의 모델을 성능 예측 결과와 WTS-4와 USWP를 소음 예측 결과와 비교하였다. 계산 결과는 실제의 실험치와 잘 일치하는 것으로 나타났다.