• 한국항공우주학회지
• /
• 제38권11호
• /
• pp.1094-1104
• /
• 2010

익형주위 착빙 현상과 관련한 기존 연구는 수치적으로 예측된 결빙형상과 실험 결과를 정성적으로 비교하는데 그쳐, 해석 결과의 정확도에 대한 판단과 외기 조건 변화에 따른 결빙형상의 체계적 분석에 한계가 있었다. 이에 본 논문에서는 실린더와 익형의 결빙 형상을 가용한 실험 및 타 수치해석 결과와 비교, 검증하고 정량적으로 분석하는 연구를 수행하였다. 먼저, 개발된 코드로부터 획득한 결빙형상을 얼음의 최대두께와 얼음의 진행방향, 얼음의 분포, 얼음의 면적을 기준으로 나타낸 후 이를 정량적으로 분석하였다. 정량적 분석을 통해 유사한 결빙 형상을 직관적으로 비교할 수 있었다. 개발한 수치해석 코드는 아랫면에서의 결빙 면적, 얼음의 두께를 작게 예측하였다. 이를 개선하기 위하여 보다 정밀한 유동장 계산을 통해 획득한 입자의 궤적이 요구된다.

• 대한기계학회논문집
• /
• 제11권5호
• /
• pp.755-761
• /
• 1987

본 연구에서는 익형 위에 발생하는 박리기포 주위를 사용한 박리기포 주위 유 동해석에 목적을 두고, 원시변수(primitive variable)를 사용한 부분 포물형 Navier -Stokes 방정식을 사용하여 층류유동에 관한 간단한 기본계산을 통해 비교적 박리기포 가 큰 외부유동(external flow)에도 부분 포물형 방정식이 적용될 수 있음을 보이고자 한다.수치해법은 Galpin 등 이 이차원 관유동(channel flow)에 완전 Navier-Stokes 방정식의 해법으로 사용한 CELS(coupled equation line solver) 방법을 부분 포물형 방정식에 적합하게 수정하여 사용하였다.

• 대한기계학회논문집
• /
• 제12권2호
• /
• pp.286-294
• /
• 1988

본 논문에서는 Bogdanoff가 제시한 열응력 관계식을 이용하여 균일 열유동을 받는 무한체내에 존재하는 균열에 대해서 균열표면에 자유트랙션조건과 열경계 조건을 부여하여 TSIF를 구하고자 한다. 다루고자 하는 커스프 균열의 형태는 하이포사이클로이드형(hypocycloid type), 대칭익형(symmetric airfoil type), 그리고 대칭입술형(symmetric lip type)이다.

• 한국전산유체공학회지
• /
• 제4권1호
• /
• pp.53-61
• /
• 1999

Two-dimensional, unsteady, incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. The compressible code involves a conventional upwind-differenced scheme for the convective terms and LU-SGS scheme for temporal integration. The incompressible code with pseudo-compressibility method also adopts the same schemes as the compressible code. Three two-equation turbulence models are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by predicting the flow around the RAE 2822 transonic airfoil and the NACA 4412 airfoil, respectively. In addition, both the incompressible and compressible code are used to compute the flow over the NLR 7301 airfoil with flap to study the compressible effect near the high-loaded leading edge. The grid systems are efficiently generated using Chimera overlapping grid scheme. Overall, the κ-ω SST model shows closer agreement with experiment results, especially in the prediction of adverse pressure gradient region on the suction surfaces of high-lift airfoils.

• 대한조선학회논문집
• /
• 제36권3호
• /
• pp.50-59
• /
• 1999

본 논문은 비선형 계획법의 하나인 SQP법을 이용한 지면(혹은 해면)효과를 받는 익형(WIG)의 최적단면생성에 관한 연구 결과이다. 유동장의 해석기법으로는 포텐셜패널법을 사용하였고, 지면효과는 경상법(image method)을 사용하여 처리하였다. 수치계산의 검증을 위하여 이미 알려진 날개단면의 압력분포를 주고 그 단면을 역으로 찾아가도록(역해법) 하였다. 본 논문의 목적은 주어진 구속조건(설계조건)중 특히 WIG선에서 중요한 종방향 안정성을 만족하면서 최대양력을 주는 익형의 단면을 생성하는 것이며 꼬리날개의 영향을 고려하였다.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.54-60
• /
• 2012

Transition prediction results are validated with experimental data obtained from a transonic wind tunnel for the KU109C airfoil. A Reynolds-Averaged Navier-Stokes code is simultaneously coupled with the transition transport model of Langtry and Menter and applied to the numerical prediction of aerodynamic performance of the KU109C airfoil. Drag coefficients from the experiment are better correlated to the numerical prediction results using a transition transport model rather than the fully turbulent simulation results. Maximum lift coefficient and drag divergence at the zero-lift condition with Mach number are investigated. Through the present validation procedure, the accuracy and usefulness of both the experiment and the numerical prediction are assessed.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.479-482
• /
• 2002

Riblets with longitudinal grooves along the streamwise direction have been used as an effective flow control technique for drag reduction. A flexible micro-riblet with v-grooves of peak-to-peak spacing of $300{\mu}m$ was made using a MEMS fabrication process of PDMS replica. The flexible micro-riblet was attached on the whole surface of a NACA0012 airfoil with which grooves are aligned with the streamwise direction. The riblet surface reduces drag coefficient about $7.9{\%}\;at\;U_o=3.3m/s$, however, it increases drag about $8{\%}\;at\;U_o=7.0m/s$, compared with the smooth airfoil without riblets. The near wake has been investigated experimentally far the cases of drag reduction ($U_o\;=\;3.3 m/s$) and drag increase ($U_o\;=\;7 m/s$). Five hundred instantaneous velocity fields were measured for each experimental condition using the cross-correlation PIV velocity field measurement technique. The instantaneous velocity fields were ensemble averaged to get spatial distribution of turbulent statistics such as turbulent kinetic energy. The experimental results were compared with those of a smooth airfoil under the same flow condition. The micro-riblet surface influences the near wake flow structure largely, especially in the region near the body surface

• 대한기계학회논문집B
• /
• 제25권6호
• /
• pp.776-785
• /
• 2001

Hot-wire measurements are performed in boundary layers developing on a NACA0012 airfoil over which wakes pass periodically. The Reynolds number based on chord length of the airfoil is 2$\times$10(sup)5 and the wakes are generated by circular cylinders rotating clockwise and counterclockwise around the airfoil. This paper and its companion Part II describe the phenomena of wake-induced transition of the boundary layers on the airfoil using measured data; phase-and time-averaged streamwise mean velocities, turbulent fluctuations, integral parameters and wall skin frictions. This paper describes the background and facility together with results of time-averaged quantities. Due to the passing wake with mean velocity defects and high turbulence intensities, the laminar boundary layer is periodically disturbed at the upstream station and becomes steady-state transitional boundary layer at the downstream station. The velocity defect in the passing wake changes the local pressure at the leading of the airfoil, significantly affects the time-mean pressure distribution on the airfoil and eventually, has influence on the transition process of the boundary layer.

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

The effects of the frequency of upstream gust on the unsteady boundary characteristics on a downstream blade was simulated by using a Navier-Stokes code. The Navier-Stokes code is based on an unstructured finite volume method and uses a low Reynolds k-e turbulence model to close the momentum equations. The MIT flapping foil experiment set-up is used to simulate the interaction between the upstream wake and a blade. The frequency of the upstream wake is simulated by varying rate of pitching motion of the flapping airfoils. Three reduced frequencies. 3.62. 7.24. and 10.86. are simulated. As the frequency increases, the unsteady fluctuation on the surfaces of the downstream hydrofoil is shown to decrease while the upstream flapper wake has larger first harmonics of y-velocity component. The unsteady vortices are shown to interact with each other and. as a result. the upstream wake becomes undiscernible inside the inner layer. The turbulence kinetic energy shows a similar behavior. Limiting streamlines around the trailing edge of the flapper are shown to conform with the unsteady Kutta condition for a round trailing edge. while limiting streamlines around the trailing edge of the hydrofoil conforms with the unsteady Kutta condition for a sharp edge.