• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.186-191
• /
• 1999

An Euler solver is developed to predict accurate aerodynamic data such as lift coefficient, drag coefficient, and moment coefficient. The conservation law form of the compressible Euler equations are used in the generalized curvilinear coordinates system. The Euler solver uses a finite volume method and the second order Roe's flux difference splitting scheme with min-mod flux limiter to calculate the fluxes accurately. An implicit scheme which includes the boundary conditions is implemented to accelerate the convergence rate. The multi-block grid is integrated into the flow solver for complex geometry. The flowfields are analyzed around NACA 0012 airfoil in the cases of $M_{\infty}=0.75,\;\alpha=2.0\;and\;M_{\infty}=0.80,\;\alpha=1.25$. The numerical results are compared with other numerical results from the literature. The final goal of this research is to prepare a robust and an efficient Navier-Stokes solver eventually.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.90-97
• /
• 1999

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.171-176
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler , thin layer Navier- Stokes and full Navier-Stokes ones, are solved using implicit LU-ADI decomposition scheme. The feasible direction algorithm with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• Proceeding of EDISON Challenge
• /
• 2016.11a
• /
• pp.35-43
• /
• 2016

본 연구에서는 EDISON의 해석자를 활용하여 선박의 최적 운항 자세 도출과 같은 실용적 문제 해석뿐 만 아니라 자유 수면 및 대칭 경계 조건이 수직한 몰수 익형 주위 유동에 미치는 영향을 분석하는 학술적 연구를 수행하였다. 선박의 자세에 따른 저항 변화를 분석한 결과 0.5m 선수 트림에서 선체 저항이 가장 작은 것으로 나타났으며, 이는 유동 가속에 의한 선수 어깨부의 낮은 압력 및 선미부에서의 압력 회복에 의한 것이다. 반면, 1.0m 선미 트림에서 선체 저항이 가장 큰 것으로 나타났는데, 평형 상태보다 선미부의 압력 회복이 약하기 때문이다. 또한 자유 수면과 대칭 경계 조건이 날개 성능에 미치는 영향을 분석한 결과, 비현실적 대칭경계 조건으로 인해 날개 양력이 13%~16% 크게 나타났으므로 대칭경계조건을 사용할 경우에는 이러한 오차를 감안해야 한다.

• Proceeding of EDISON Challenge
• /
• 2016.03a
• /
• pp.638-642
• /
• 2016

현재 양력발생의 원리는 긴 경로 이론과 동시 통과 이론 등으로 알려져 있으며, 이 이론들은 간단히 설명하려는 목적으로 인해 양력발생의 원리가 와전되어왔다. 본 연구에서는 여러 양력발생의 원리 중 Kutta-Joukowski의 이론을 이용하여 순환을 통한 양력발생의 관계에 대해 연구를 수행하였다. 해석 모델은 대칭형 에어포일인 NACA0012를 이용하였다. 해석 프로그램으로는 EDISON을 이용하였으며 시간에 따른 와도를 확인하기 위해 비정상상태 해석을 진행하였다. 그 결과 유동이 시간의 경과에 따라 에어포일 뒷전으로 밀려나면서 Kutta-Joukowski의 이론에 중요한 논점인 출발 와류가 생성됨을 확인할 수 있었다. 이로 인해 순환이 형성되어 양력이 발생함을 추측할 수 있다. 추가적으로 받음각과 출발 와류 사이의 영향성을 연구하였다. 받음각이 커짐에 따라 출발 와류가 에어포일 뒷전에서 더 크게 발생하였고, 받음각이 낮아질수록 출발 와류의 진행이 가속화됨을 확인할 수 있었다.

• Proceeding of EDISON Challenge
• /
• 2016.03a
• /
• pp.593-597
• /
• 2016

항공기의 비행 안정성을 판별하기 위해서는 정적(static)안정성과 동적(dynamic)안정성을 모두 고려해야 한다. 본 연구는 EDISON_CFD의 KFLOW_EDISON_2D3DOF, 2D_Com_P해석자를 이용해 2차원익형의 정적(static)계수와 동적(dynamic)계수를 구해 항공기 전체형상의 정적 안정성 및 동적 안정성을 분석하였다. 그 중 2015년도 인간 동력 항공기대회 건국대학교 참가 기체를 예시로 연구를 진행하였다. 주익 익형 SG6043과 미익 익형 NACA0012을 EDISON_CFD를 이용하여 각 익형의 모멘트계수를 구하고 그 결과값을 이용하여 기체의 무게중심에 대한 각 익형의 감쇠계수를 계산하였다. 그리고 주익과 미익의 면적비율과 코드길이의 비율을 고려해 전체 항공기형상의 안정성과 비행 형태를 분석하였다.

• Journal of Wind Energy
• /
• v.5 no.1
• /
• pp.33-42
• /
• 2014

The new aerofoil, KA2 was designed to apply to the wind turbine blade. For the aerofoil, numerical analysis was performed to review aerodynamic characteristics like lift and drag coefficient. And they are verified with test data using the digital wind tunnel and test samples from 3D printer. The digital wind tunnel was developed to test wing in the small laboratory, and verified with test of NACA0012 airfoil. KA2 aerofoil is asymmetric, and has the thickness ratio of 14%, and 12 degree of AOA at the maximum lift coefficient of 1.3. In this paper, aerodynamic characteristics from numerical and test approaches will be proposed with AOA in detail. Therefore, this aerofoil will be used for the design of wind turbine blade.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.12
• /
• pp.785-792
• /
• 2016

In the present study, the effects of non-equilibrium condensation on the drag divergence Mach number with the angle of attack in a transonic 2D moist air flow of NACA0012 are investigated using the TVD finite difference scheme. For the same ${\alpha}$, the maximum upstream Mach number of the shock wave, Mmax, and the size of supersonic bubble decrease with the increase in ${\Phi}_0$. For the same $M_{\infty}$, ${\Phi}_0$, and $T_0$, the length of the non-equilibrium condensation zone ${\Delta}_z$ decreases with increasing ${\Phi}_0$. On the other hand, because of the attenuating effect of non-equilibrium condensation on wave drag, which is related to the interaction between the shock wave and the boundary layer, the drag coefficient $C_D$ decreases with an increase in ${\Phi}_0$ for the same $M_{\infty}$ and ${\alpha}$. For the same ${\alpha}$, $M_D$ increases with increasing ${\Phi}_0$, while $M_D$ decreases with an increase in ${\alpha}$.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.12
• /
• pp.1009-1015
• /
• 2014

The present study investigated the effects of non-equilibrium condensation with the angle of attack on the coefficients of pressure, lift, and drag in the transonic 2-D flow of NACA0012 by numerical analysis of the total variation diminishing (TVD) scheme. At $T_0=298k$ and ${\alpha}=3^{\circ}$, the lift coefficients for $M_{\infty}=0.78$ and 0.81 decreased monotonically with increasing ${\Phi}_0$. In contrast, for $M_{\infty}$ corresponding to the Mach number of the force break, $C_L$ increased with ${\Phi}_0$. For ${\alpha}=3^{\circ}$ and ${\Phi}_0=0%$, $C_D$ increased markedly as $M_{\infty}$ increased. However, at ${\Phi}_0=60%$ and ${\alpha}=3^{\circ}$, which corresponded to the case of the condensation having a large influence, $C_D$ increased slightly as $M_{\infty}$ increased. The decrease in profile drag by non-equilibrium condensation grew as the angle of attack and stagnation relative humidity increased for the same free stream transonic Mach number. At ${\Phi}_0=0%$, the coefficient of the wave drag increased with the attack angle and free stream Mach number. When ${\Phi}_0$ > 50%, the coefficient of the wave drag decreased as ${\alpha}$ and $M_{\infty}$ increased. Lowering ${\Phi}_0$ and increasing $M_{\infty}$ increased the maximum Mach number.

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.