• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.138-144
• /
• 1998

The BGK equation, which is the kinetic model equation of Boltzmann equation, is solved using FDDO(finite difference with the discrete-ordinate method) to compute the rarefied flow of monatomic gas. Using reduced velocity distribution and discrete ordinate method, the scalar equation is transformed into a system of hyperbolic equations. High resolution ENO(Essentially Non-Oscillatory) scheme based on Harten-Yee's MFA(Modified Flux Approach) method with Strang-type explicit time integration is applied to solve the system equations. The calculated results are well compared with the experimental density field of NACA0012 airfoil, validating the developed computer code. Next. the computed results of circular cylinder flow for various Knudsen numbers are compared with the DSMC(Direct Simulation Monte Carlo) results by Vogenitz et al. The present scheme is found to be useful and efficient far the analysis of two-dimensional rarefied gas flows, especially in the transitional flow regime, when compared with the DSMC method.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.30-40
• /
• 2007

The implicit discontinuous Galerkin method for the two-dimensional Euler equations was developed on unstructured triangular meshes, which can achieve higher-order accuracy by wing hierachical basis functions based on Legendre polynomials. Numerical tests were conducted to estimate the convergence order of numerical solutions to the Ringleb flow and the supersonic vortex flow for which analytic solutions are available. And, the flows around a circle and a NACA0012 airfoil was also numerically simulated. Numerical results show that the implicit discontinuous Galerkin methods with higher-order representation of curved solid boundaries can be an efficient higher-order method to obtain very accurate numerical solutions on unstructured meshes.

• Proceeding of EDISON Challenge
• /
• 2014.03a
• /
• pp.544-549
• /
• 2014

임계마하수보다 큰 자유흐름 마하수에서는 충격파의 발생으로 인해 급격한 항력증가가 발생하므로, 임계마하수 증가는 고속 공기역학에서 중요한 분야로 다뤄지고 있다. Whitcomb R. T.에 의해 천음속영역에서 순항할 수 있는 초임계익형이 개발되었으나, 충격파 제어 기법들에 대한 실험적인 검증은 형상 제작의 어려움으로 인해 한계를 지닌다. 따라서 본 논문에서는 2D_Comp-2.1_P와 Prandtl-Glauert 압축성 보정식을 이용하여 NACA0012와 RAE2822의 임계마하수를 해석하고, 충격파 제어 장치 중 하나인 Bump를 RAE2822에 설치하여 임계마하수를 향상시키기 위한 연구를 수행하였다. 연구 결과 충격파를 압축파로 분산시켜 충격파의 강도를 약화시키고, 양항비의 4.7% 증가를 확인하였다. 따라서 Bump를 설계한 RAE2822가 기본 익형보다 높은 천음속 조건에서 효율적인 공력특성을 가지는 것을 확인하였다.

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

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.3
• /
• pp.219-225
• /
• 2014

In this study, to find the characteristics of the oscillation of a terminating shock wave in a transonic airfoil flow with non-equilibrium condensation, a NACA00-12,14,15 airfoil flow with non-equilibrium condensation is investigated through numerical analysis of TVD scheme. Transonic free stream Mach number of 0.81-0.90 with the variation of stagnation relative humidity and airfoil thickness is tested. For the free stream Mach number 0.87 and attack angle of ${\alpha}=0^{\circ}$, the increase in stagnation relative humidity attenuates the strength of the terminating shock wave and inactivates the oscillation of the terminating shock wave. For the case of $M_{\infty}=0.87$ and ${\phi}_0=60%$, the decreasing rate in the frequency of the shock oscillation caused by non-equilibrium condensation to that of ${\phi}_0=30%$ amounts to 5%. Also, as the stagnation relative humidity gets larger, the maximum coefficient of drag and the difference between the maximum and minimum in $C_D$ become smaller. On the other hand, as the thickness of the airfoil gets larger, the supersonic bubble size becomes bigger and the oscillation of the shock wave becomes higher.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.77-85
• /
• 2009

A two-dimensional cross-section analysis program based on the finite element method has been developed for composite blades with arbitrary cross-section profiles and material distributions. The modulus weighted approach is used to take into account the non-homogeneous material characteristics of advanced blades. The CLPT (Classical Lamination Plate Theory) is applied to obtain the effective moduli of the composite laminate. The location of shear center for any given cross-sections are determined according to the Trefftz' definition while the torsion constants are obtained using the St. Venant torsion theory. A series of benchmark examples for beams with various cross-sections are illustrated to show the accuracy of the developed cross-section analysis program. The cross section cases include thin-walled C-channel, I-beam, single-cell box, NACA0012 airfoil, and KARI small-scale blades. Overall, a reasonable correlation is obtained in comparison with experiments or finite element analysis results.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.23-33
• /
• 2009

The inherent aeromechanical complexity of a rotor system necessitated the comprehensive analysis code for helicopter rotor system. In the present study, an aerodynamic analysis module has been developed as a part of rotorcraft comprehensive program. Aerodynamic analysis module is largely classified into airload calculation routine and inflow analysis routine. For airload calculation, quasi-steady analysis model is employed based on the blade element method with the correction of unsteady aerodynamic effects. In order to take unsteady effects - body motion effects and dynamic stall - into account, aerodynamic coefficients are corrected by considering Leishman-Beddoes's unsteady model. Various inflow models and vortex wake models are implemented in the aerodynamic module to consider wake induced inflow. Specifically, linear inflow, dynamic inflow, prescribed wake and free wake model are integrated into the present module. The aerodynamic characteristics of each method are compared and validated against available experimental data such as Elliot's induced inflow distribution and sectional normal force coefficients of AH-1G. In order to validate unsteady aerodynamic model, 2-D unsteady model for NACA0012 airfoil is validated against aerodynamic coefficients of McAlister's experimental data.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.166-170
• /
• 2008

The changing process from hovering mode to transition one is of importance to determine a stability of tilt-rotor aircraft, which is utilized in UAV(Unmaned Aerial Vehicle). The analysis on fluid flows and aerodynamic characteristics according to variation of tilting angle of rotor is essential step in development of tilt-rotor. In the present study, the computation domain is divided into the rotating and stationary regions in order to consider the rotating blades. For the convenient realization of various tilting angle as well as application of boundary condition, the whole computation region is constructed into sphere domain. The near farfield boundary condition is adopted. The airfoil used in computation is NACA 0012. The computation results for the hovering mode are validated by comparing with previously conducted experimental results. From the results, the flow fields around rotor blade and the aerodynamic characteristics in transition mode are observed. The computational result will provide the basis for development and performance evaluation of tilt-type aircraft.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.115-119
• /
• 2004

The effect of a continuous blowing or suction on an oscillating 2-D NACA0012 airfoil was investigated numerically for the dynamic stall control. The influence of control parameter variation was also studied in the view point of aerodynamic characteristics. The result showed that the blowing control kept a higher lift drag ratio before stall angle but the dynamic stall angle was not exceed to without control result. As the slot position was closer to leading edge, the positive control effect becomes greater. The stronger jet and the smaller jet angel made more favorable roles on the control performance. In the cases of the suction, the overall control features were similar to those of the blowing, but dynamic stall angle was increased, i.e. suction was more effective to control dynamic stall. It was also founded that the suction control was showed better control effect as the slot position moves to trail edge within thirty percentage of chord length. In the simulation for the jet strength and the jet angle control, the same tendencies were observed to those of blowing cases.