• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.415-420
• /
• 2013

천음속 영역에서 비행하는 에어포일 주위의 유동은 아음속 영역과는 달리 충격파(Shockwave)를 동반하고 이에 따라 복잡한 유동 현상이 발생한다. 본 연구는 천음속 영역에서 설계된 에어포일 주위의 유동에서 충격파 발생에 따른 유동변화의 특징들을 관찰하였다. 특히, 충격파에 의해서 발생하는 진동(Buffeting)과 경계층 유동 박리(Shock-induced boundary layer separation)현상에 대한 관찰과 일반적인 에어포일과 초임계 에어포일(Supercritical airfoil)간의 유동 특성의 차이점을 분석하였다. 본 연구를 위해서 EDISON CFD코드가 사용되었다.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.354-357
• /
• 2008

Drag reduction is one of main concerns for commercial aircraft companies than ever because fuel price has been tripled in ten years. In this research, Natural Laminar Flow airfoil is designed and tested to reduce drag at cruise condition, $c_l$=0.3, Re=3.4${\times}$10$^6$ and M=0.6. NLF airfoil is characterized by delayed transition from laminar to turbulent flow, which comes from maintaining favorable pressure gradient to downstream. Transition is predicted by solving Boundary Layer equations in viscous boundary layer and by solving Euler Equation outside the boundary layer. Once boundary layer thickness and momentum thickness are obtained, $e^N$-method is used for transition point prediction. As results, KARI's NLF airfoil is designed and shows better characteristics than NLF-0115. The characteristics are tested and verified at low Reynolds numbers, but at high Reynolds numbers, laminar flow characteristics are not obtainable because of fully turbulent flow over airfoil surfaces. Precious experiences, however, relating NLF airfoil design, subsonic and transonic tests are acquired.

• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2347-2355
• /
• 1996

Purpose of the present study is to develop a computational design program for shape optimization, combining the numerical optimization technique with the flow analysis code. The present methodology is then validated in three cases of aerodynamic shape optimization. In the numerical optimization, a feasible direction optimization algorithm and shape functions are considered. In the flow analysis, the Navier-Stokes equations are discretized by a cell-centered finite volume method, and Roe's flux difference splitting TVD scheme and ADI method are used. The developed design code is applied to a transonic channel flow over a bump, and an external flow over a NACA0012 airfoil to minimize the wave drag induced by shock waves. Also a separated subsonic flow over a NACA0024 airfoil is considered to determine a maximum allowable thickness of the airfoil without separation.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.161-164
• /
• 2002

In the present study, a passive control method, using the porous wall and cavity system, is applied to the shock wave/boundary layer interactions in transonic moist air flow. The two-dimensional, unsteady, compressible Navier-Stokes equations, which are fully coupled with a droplet growth equation, are solved by the third-order MUSCL type TVD finite difference scheme. Baldwind-Lomax turbulence model is employed to close the governing equations. In order to investigate the effectiveness of the present control method, the total pressure losses of the flow and the time-dependent behaviour of shock motions are analyzed in detail. The computed results show that the present passive control method considerably reduces the total pressure losses due to the shock/boundary layer interaction in transonic moist air flow and suppresses the unsteady shock wave motions over the airfoil, as well. It is also found that the location of the porous ventilation significantly influences the control effectiveness.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Journal of Mechanical Science and Technology
• /
• v.15 no.3
• /
• pp.366-373
• /
• 2001

A numerical analysis of shock wave/boundary layer interaction in transonic/supersonic axial flow compressor cascade has been performed by using a characteristics upwind Navier-Stokes method with various turbulence models. Two equation turbulence models were applied to transonic/supersonic flows over a NACA 0012 airfoil. The results are superion to those from an algebraic turbulence model. High order TVD schemes predicted shock wave/boundary layer interactions reasonably well. However, the prediction of SWBLI depends more on turbulence models than high order schemes. In a supersonic axial flow cascade at M=1.59 and exit/inlet static pressure ratio of 2.21, k-$\omega$ and Shear Stress Transport (SST) models were numerically stables. However, the k-$\omega$ model predicted thicker shock waves in the flow passage. Losses due to shock/shock and shock/boundary layer interactions in transonic/supersonic compressor flowfields can be higher losses than viscous losses due to flow separation and viscous dissipation.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.31-34
• /
• 2002

An efficient implicit multigrid method is presented for the Navier-Stokes and $k-{\omega}$ turbulence equations. Freezing and limiting strategies are applied to improve the robustness and convergence of the multigrid method. The eddy viscosity and strongly nonlinear production terms of turbulence are frozen in the coarser grids by passing down the values without update of them. The turbulence equations together with the Navier-Stokes equations, however, are consecutively solved on the coarser grids in a loosely coupled fashion. A simple limit for k is also introduced to circumvent slow-down of convergence. Numerical results for the unseparated and separated transonic airfoil flows show that all computations converge well without any robustness problem and the computing time is reduced to a factor of about 3 by the present multigrid method.

• Journal of computational fluids engineering
• /
• v.4 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.1
• /
• pp.61-66
• /
• 2012

In this study, to find the effect of nonequilibrium condensation on the oscillation of the terminating shock wave in transonic flows, an NACA0014 airfoil flow with nonequilibrium condensation is analyzed using the total variation diminishing (TVD) numerical scheme. Transonic free stream Mach numbers of 0.81-0.87 are tested with variations in the stagnation relative humidity. For the same free stream Mach number and attack angle of ${\alpha}=0^{\circ}$, an increase in the stagnation relative humidity attenuates the strength of the terminating shock and reduces the oscillation of the terminating shock wave. Furthermore, for the same stagnation relative humidity, the larger the free stream Mach number becomes, the shorter the period of the oscillation shock wave is. The excursion distance of the oscillation shock increases with the free stream Mach numbers for the same stagnation relative humidity. Finally, it is found that for the same shock location, the strength of the oscillating shock facing upstream is stronger than that facing downstream.