• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• 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).

• Ocean Systems Engineering
• /
• 제8권2호
• /
• pp.167-182
• /
• 2018

RTI (Rayleigh-Taylor instability) is investigated by a multi-liquid MPS (Moving Particle Semi-implicit) method for both viscous and inviscid flows for various density differences, initial-disturbance amplitudes, viscosities, and surface tensions. The MPS simulation can be continued up to the late stage of high nonlinearity with complicated patterns and its initial developments agree well with the linear theoretical results. According to the relevant linear theory, the difference between inviscid and viscous fluids is the rising velocity at which upward-mushroom-like RTI flow with vortex formation is generated. However, with the developed MPS program, significant differences in both growing patters and developing speeds are observed. Also, more dispersion can be observed in the inviscid case. With larger Atwood (AT) number, stronger RTI flows are developed earlier, as expected, with higher potential-energy differences. With larger initial disturbances, quite different patterns of RTI-development are observed compared to the small-initial-disturbance case. If AT number is small, the surface tension tends to delay and suppress the RTI development when it is sufficiently large. Interestingly, at high AT number, the RTI-suppressions by increased surface tension become less effective.

• 한국항공우주학회지
• /
• 제48권1호
• /
• pp.1-11
• /
• 2020

항공기에서 발생한 후류의 거동은 항공기의 성능과 비행안정성에 큰 영향을 미친다. 본 연구에서는 푸리에-스펙트럴법을 사용하여 근접 편대비행을 하고 있는 항공기 날개 후류 거동에 관하여 연구하였다. 초기와의 순환강도, 상대위치 등에 따라 와들의 거동이 복잡한 양상으로 나타났다. 와의 순환강도 값이 큰 경우 와들의 이동이 크게 나타났다. 초기에는 가까이 위치한 와의 영향을 받아 이동하지만 이동하는 과정에서 다른 와와 가깝게 되면 새로운 경로를 형성하였다. 점성이 클수록 와의 반경이 증가하여 와 반경 근처로 새롭게 진입하는 와와 상호작용이 강해진다. 향후 항공기가 이착륙 시에 발생하는 지면 효과를 고려한 후류 거동 해석 연구를 수행하고자 한다.

• 대한기계학회논문집B
• /
• 제38권9호
• /
• pp.731-737
• /
• 2014

와류실은 압축된 공기를 이용하여 고온과 저온 가스로 분리할 수 있는 단순한 장치로, 차세대 새로운 열교환기로 각광받고 있으나, 와류실 내부에서 발행하는 물리적 유동특성에 대해 아직까지 많이 알려지지 않았다. 본 연구에서는 온도 분리 현상을 조사하기 위해 실험 및 수치해석을 수행하였다. 공급 압력에 따른 온도 변화를 측정하기 위하여 다수의 압력 및 온도 센서를 사용하였으며, CFD 기법을 적용하여 3차원 비정상 압축성 유동장을 조사하였다. 연구를 통해 온도 분리 현상은 점성일과 밀접한 관계가 있는 공급 압력과 와류실의 직경에 영향을 받았으며, 와류실에서 발생하는 온도분리 현상은 압력구배파의 개념으로 확증할 수 없었다.

• 산업기술연구
• /
• 제15권
• /
• pp.153-168
• /
• 1995

The purpose of this study is to develop a method for predicting the aerodynamic performance of the subsonic airfoils in the 2-dimensional, steady and viscous flow. For this study, the airfoil geometry is specified by adopting the longest chord line system and by considering local surface curvature. In case of the inviscid-incompressible flow, the analysis is accomplished by the linearly varying strength vortex panel method and the Karman-Tsien correction law is applied for the inviscid-compressible flow analysis. The Goradia's integral method and the Truckenbrodt integral method are adopted for the boundary layer analysis of the laminar flow and the turbulent flow respectively. Viscous and inviscid solutions are converged by the Lockheed iterative calculating method using the equivalent airfoil geometry. And the analysis of the seperated flow is performed using the Dvorak and Maskew's method as the basic method. The wake effect is also considered and its geometry expressed by the formula of Summey & Smith when no seperation occurs. A computational efficiency is verified by the comparison of the computational results with experimental data and by the shorter execution time.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2002년도 유체기계 연구개발 발표회 논문집
• /
• pp.244-249
• /
• 2002

It is known that tip clearance flows reduce the pressure rin, flow range and efficiency of the turbomachinery. So, the clear understanding about flow fields in the tip region is needed to efficiently design the turbomachinery. The Navier-Stokes code with the proper treatment of the boundary conditions has been developed to analyze the three-dimensional steady viscous flow fields in the transonic rotating blades and a numerical study has been conducted to investigate the detail flow physics in the tip region of transonic rotor, NASA Rotor 67. The computational results in the tip region of transonic rotors show the leakage vortices, leakage flow from pressure side to suction side and their interaction with a shock Depending on the operating conditions, the position of shock-wave on the blade surface are v8y different close to the blade tip of the transonic compressor rotor. The shock-wave position dose to the blade tip had the dose relationship with the starting position of leakage vortex and the direction of leakage flow.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2007년도 추계 학술대회논문집
• /
• pp.171-176
• /
• 2007

Supersonic viscous flow over a 5 degree half angle cone studied computationally with three-dimensional Navier-Stokes equations. Steady asymmetric solutions of 5-deg half angle cone show that the asymmetric flow separation is caused by convective instability. The angle of attack, Reynolds number, and Mach number affected the side force variation that is caused by asymmetric vortical flow.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 춘계학술대회논문집
• /
• pp.25-28
• /
• 2008

An aerodynamic calculation in hover of KUH main rotor blade is performed using a three-dimensional unstructured hybrid mesh viscous flow solver. The flow solver utilizes a vertex-centered finite-volume scheme that is based on the Roe's flux-difference splitting with an implicit Jacobi/Gauss-Seidel time integration. The eddy viscosity are estimated by the Spalart-Allmaras one-equation turbulence model. A solution-adaptive mesh refinement technique is used for efficient capturing of the tip vortex. Calculations are performed at several operating conditions with varying collective pitch setting for KUH main rotor blade in hover. Good agreements are obtained between the present and other results using HOST and CAMRAD II in overall rotor performance. It is demonstrated that the present vertex-centered flow solver is an efficient and accurate tool for the assessment of rotor performance in hover.

• 대한기계학회논문집B
• /
• 제28권12호
• /
• pp.1600-1607
• /
• 2004

Performance of a regenerative pump is evaluated based on the calculated through flows using the CFX-TASCflow code. Flow calculations are performed in one vane to vane space of the impeller and side channel. The flow is very complex three dimensional with a strong radial vortex due to centrifugal force and an axial vortex due to re-circulating flow between the impeller and the side channel. Momentum exchange on the plane between the impeller and the side channel are evaluated to estimate design parameters and viscous losses in the pump. The present study contributes to showing the capability of flow simulation of complex flow in the regenerative pump by comparing the calculated performance with the measured value.

• 동력기계공학회지
• /
• 제9권1호
• /
• pp.23-29
• /
• 2005

Numerical analysis of three-dimensional viscous flow-fields in the turbine rotor passages was carried out to investigate flow physics including the interaction between secondary vortices, tip leakage vortex, and the rotor wake. The blade tip geometry was accurately modeled adopting the embedded H grid system. An explicit four-stage Runge-Kutta scheme was used for the time integration of both the mean flow and turbulence equations. The computational results for the entire turbine rotor flows, particularly the tip clearance flow and the secondary flows, were interpreted and compared with the experimental data from the Penn State turbine stage. The predictions for major features of the flow field have been found to be in good agreement with the experimental data.