21.11 익형(Airfoil)의 기하학과 분류 (Airfoil Geometry and Classification)

21.11 익형(Airfoil)의 기하학과 분류 (Airfoil Geometry and Classification)

1. 익형의 정의와 공기역학적 역할

익형(airfoil)은 날개(wing)나 블레이드(blade)의 횡단면(cross-section) 형상을 의미하며, 유동 내에서 양력(lift)을 생성하고 항력(drag)을 최소화하도록 설계된 2차원 형상이다. 익형은 고정익 항공기의 날개, 헬리콥터 및 멀티로터의 로터 블레이드, 프로펠러, 풍력 터빈 블레이드 등 양력을 발생시키는 모든 공력 구조물의 기본 구성 요소이다.

익형의 기하학적 형상은 표면 압력 분포(pressure distribution), 경계층(boundary layer) 거동, 양력 및 항력 특성, 실속(stall) 거동 등 공기역학적 성능 전반을 결정한다. 따라서 익형의 기하학적 매개변수(geometric parameters)를 정확히 정의하고 이해하는 것은 공기역학적 해석과 설계의 출발점이 된다.

2. 익형의 기하학적 매개변수

익형의 형상을 기술하는 주요 기하학적 매개변수는 다음과 같다:

시위선(chord line)과 시위 길이(chord length, c): 시위선은 익형의 앞전(leading edge)과 뒷전(trailing edge)을 연결하는 직선이다. 시위 길이 c는 이 직선의 길이로서, 익형의 기본 참조 길이이며, 모든 기하학적 매개변수의 정규화에 사용된다.

앞전(leading edge): 익형의 전방 끝점으로, 일반적으로 원호(circular arc)에 가까운 곡률을 가진다. 앞전 반경(leading edge radius, r_\text{LE})은 앞전 부근의 곡률 반경으로 정의되며, 이 값은 익형의 실속 특성에 중대한 영향을 미친다. 앞전 반경이 큰 익형은 넓은 받음각 범위에서 점진적 실속(gentle stall)을 나타내는 반면, 앞전 반경이 작은 익형은 급격한 앞전 실속(leading edge stall)을 나타내는 경향이 있다.

뒷전(trailing edge): 익형의 후방 끝점으로, 이상적으로는 첨예한(sharp) 형태를 가진다. 쿠타 조건(Kutta condition)의 만족을 위해 뒷전에서 상면과 하면의 유동이 매끄럽게 합류하여야 하며, 뒷전의 형상은 이 조건의 실현에 영향을 미친다. 실용적으로는 구조적 강도와 제작상의 이유로 유한한 두께의 뒷전(blunt trailing edge)을 갖는 경우가 있다.

캠버선(camber line, 또는 mean line): 익형 상면과 하면의 중간점을 연결하는 곡선이다. 캠버선은 익형의 비대칭성을 나타내며, 캠버의 존재는 영받음각(zero angle of attack)에서도 양력을 생성하게 하는 핵심 요인이다.

최대 캠버(maximum camber)와 최대 캠버 위치: 캠버선과 시위선 사이의 최대 수직 거리를 최대 캠버라 하며, 시위 길이에 대한 백분율로 표현한다. 최대 캠버의 시위 방향 위치는 익형의 압력 분포와 모멘트 특성에 영향을 미친다.

두께 분포(thickness distribution): 시위 방향 각 위치에서 캠버선에 수직한 방향으로 측정한 상면과 하면 사이의 거리이다. 최대 두께(maximum thickness, t_\text{max})는 시위 길이에 대한 백분율(t/c)로 표현되며, 최대 두께의 시위 방향 위치는 익형의 압력 분포 형태를 결정하는 주요 인자이다.

3. 익형 좌표계와 기하학적 구성

익형의 형상은 상면(upper surface)과 하면(lower surface)의 좌표로 기술된다. 시위선을 x축으로, 시위선에 수직한 방향을 y축으로 설정하면, 상면 좌표 y_u(x)와 하면 좌표 y_l(x)로 익형 형상이 완전히 정의된다.

캠버선 y_c(x)와 두께 분포 t(x)를 이용하면 상면과 하면의 좌표를 다음과 같이 표현할 수 있다:

x_u = x - \frac{t(x)}{2} \sin\theta, \quad y_u = y_c + \frac{t(x)}{2} \cos\theta

x_l = x + \frac{t(x)}{2} \sin\theta, \quad y_l = y_c - \frac{t(x)}{2} \cos\theta

여기서 \theta = \arctan\left(\frac{dy_c}{dx}\right)는 캠버선의 국소 기울기 각도이다. 캠버가 작은 경우(\theta \approx 0), 이 식은 다음과 같이 간략화된다:

y_u \approx y_c + \frac{t(x)}{2}, \quad y_l \approx y_c - \frac{t(x)}{2}

21.11.4 NACA 4자리 계열 익형

NACA(National Advisory Committee for Aeronautics) 4자리 계열(4-digit series)은 1933년에 개발된 체계적 익형 분류 체계로서, 익형의 기하학적 특성을 네 자리 숫자로 부호화한다 (Jacobs et al., 1933). NACA 4자리 익형의 명명 규칙은 다음과 같다:

  • 첫째 자리: 최대 캠버를 시위 길이의 백분율로 나타낸 값이다.
  • 둘째 자리: 최대 캠버 위치를 시위 길이의 10분의 1 단위로 나타낸 값이다.
  • 셋째·넷째 자리: 최대 두께를 시위 길이의 백분율로 나타낸 값이다.

예를 들어, NACA 2412 익형은 최대 캠버가 시위 길이의 2%, 최대 캠버 위치가 앞전에서 시위 길이의 40% 지점, 최대 두께가 시위 길이의 12%인 익형이다. NACA 0012와 같이 앞의 두 자리가 0인 익형은 캠버가 없는 대칭 익형(symmetric airfoil)이다.

NACA 4자리 계열의 두께 분포는 다음의 다항식으로 정의된다:

\frac{t(x)}{c} = \frac{t_\text{max}/c}{0.20} \left[ 0.2969\sqrt{\frac{x}{c}} - 0.1260\left(\frac{x}{c}\right) - 0.3516\left(\frac{x}{c}\right)^2 + 0.2843\left(\frac{x}{c}\right)^3 - 0.1015\left(\frac{x}{c}\right)^4 \right]

캠버선은 최대 캠버 위치 p 기준으로 전방과 후방에서 서로 다른 포물선으로 정의된다:

y_c = \begin{cases} \frac{m}{p^2}\left[2p\left(\frac{x}{c}\right) - \left(\frac{x}{c}\right)^2\right], & 0 \leq \frac{x}{c} \leq p \\ \frac{m}{(1-p)^2}\left[(1-2p) + 2p\left(\frac{x}{c}\right) - \left(\frac{x}{c}\right)^2\right], & p < \frac{x}{c} \leq 1 \end{cases}

여기서 m은 최대 캠버, p는 최대 캠버의 시위 방향 위치이다.

21.11.5 NACA 5자리 및 수정 계열 익형

NACA 5자리 계열(5-digit series)은 4자리 계열의 캠버선 설계를 확장하여 특정 설계 양력 계수(design lift coefficient)를 구현하도록 개발되었다. 5자리 계열의 명명 규칙에서 첫째 자리는 설계 양력 계수의 3/20에 해당하며, 둘째·셋째 자리는 최대 캠버 위치를, 넷째·다섯째 자리는 최대 두께를 나타낸다. 예를 들어, NACA 23012는 설계 양력 계수 C_{L,\text{design}} = 0.3, 최대 캠버 위치가 시위의 15% 지점, 최대 두께가 12%인 익형이다.

NACA 6자리 계열(6-series)은 층류 유동(laminar flow)의 범위를 최대화하도록 설계된 층류 익형(laminar airfoil)이다. 이 계열은 압력 분포의 역압력 구배(adverse pressure gradient)를 지연시켜 층류 경계층을 시위 방향으로 더 넓은 범위에서 유지하는 것을 설계 목표로 한다. 대표적인 NACA 6자리 익형으로는 NACA 63-215, NACA 64-210 등이 있다 (Abbott & von Doenhoff, 1959).

21.11.6 현대 익형 설계: 초임계 익형과 저레이놀즈 수 익형

초임계 익형(supercritical airfoil): Whitcomb(1974)에 의해 개발된 초임계 익형은 천음속(transonic) 영역에서의 항력 발산(drag divergence)을 지연시키기 위해 설계되었다. 초임계 익형은 상면의 곡률을 줄여 흡입 피크(suction peak)를 완화하고, 상면에서의 초음속 영역을 약한 충격파로 종결시키며, 하면의 후방에 강한 캠버를 부여하여 양력을 보상한다. 이 설계 철학은 임계 마하 수(critical Mach number)를 높여 순항 효율을 향상시킨다.

저레이놀즈 수 익형(low Reynolds number airfoil): 소형 무인 비행 로봇이 운용되는 저레이놀즈 수(Re < 5 \times 10^5) 영역에서는 층류 박리 거품(laminar separation bubble)의 형성과 층류-난류 천이(laminar-turbulent transition) 현상이 익형 성능을 지배한다. Selig, Eppler 등의 연구자들은 이 영역에 최적화된 다수의 익형을 개발하였다. 대표적으로 Eppler 387, Selig S1223, Wortmann FX 63-137 등이 있으며, 이들은 소형 고정익 UAV의 날개 설계에 널리 활용된다 (Selig et al., 1995).

21.11.7 대칭 익형과 캠버 익형의 공력 특성 비교

대칭 익형(symmetric airfoil): 캠버가 없는 대칭 익형은 시위선에 대해 상하 대칭인 형상을 가진다. 영받음각에서 양력이 영이 되며, 양력 곡선(C_L vs. \alpha)이 원점을 통과한다. 영모멘트 계수(C_{M,\text{ac}})는 영이다. NACA 0012, NACA 0015 등이 대표적인 대칭 익형이다.

캠버 익형(cambered airfoil): 양의 캠버를 가진 익형은 영받음각에서도 양의 양력을 생성한다. 양력이 영이 되는 받음각, 즉 영양력 받음각(zero-lift angle of attack, \alpha_{L=0})은 음의 값을 가진다. 캠버는 양력 곡선을 좌측으로 이동시키며, 영모멘트 계수는 일반적으로 음의 값(기수 하강 모멘트)을 나타낸다.

얇은 익형 이론(thin airfoil theory)에 따르면, 양력 곡선의 기울기 dC_L/d\alpha는 대칭 익형과 캠버 익형 모두에서 2\pi rad^{-1}(\approx 0.1097 deg^{-1})으로 동일하며, 캠버는 양력 곡선의 절편만을 변화시킨다 (Anderson, 2017).

21.11.8 익형 두께비의 공력적 영향

익형의 최대 두께비(t/c)는 공력 성능에 복합적인 영향을 미친다:

  • 양력 특성: 두께비가 증가하면 실제 양력 곡선 기울기가 얇은 익형 이론의 예측값 2\pi보다 다소 커지며, 최대 양력 계수(C_{L,\text{max}})가 증가하는 경향을 보인다. 이는 두꺼운 익형의 큰 앞전 반경이 앞전 부근의 유동 박리를 지연시키기 때문이다.
  • 항력 특성: 두께비 증가는 형상 항력(form drag)의 증가를 초래한다. 두꺼운 익형은 더 큰 역압력 구배를 형성하여 경계층이 두꺼워지고 전단 응력 분포가 변화한다.
  • 모멘트 특성: 두께 분포는 압력 분포에 영향을 미쳐 공력 모멘트에도 간접적 영향을 준다.
  • 구조적 측면: 두께비가 큰 익형은 날개 내부에 더 큰 구조 공간을 제공하므로, 구조 중량을 줄이면서 충분한 강성과 강도를 확보할 수 있다.

최적의 두께비는 공력 성능, 구조 요구 사항, 운용 레이놀즈 수 등을 종합적으로 고려하여 결정된다.

21.11.9 최대 두께 위치의 영향

최대 두께의 시위 방향 위치는 익형의 압력 분포 형태를 결정하는 핵심 인자이다:

  • 최대 두께가 전방에 위치한 익형: 앞전 부근에서 급격한 가속과 강한 흡입 피크가 형성되며, 이후 긴 역압력 구배 영역이 뒤따른다. 이러한 압력 분포는 높은 최대 양력 계수를 제공하지만, 경계층의 조기 천이를 유발하여 마찰 항력이 증가할 수 있다.
  • 최대 두께가 후방에 위치한 익형: 상면에서 완만한 순압력 구배(favorable pressure gradient)가 넓은 범위에 걸쳐 유지되어 층류 경계층을 더 오래 보존할 수 있다. NACA 6자리 계열 층류 익형이 이러한 설계 철학을 대표하며, 설계 조건에서 현저히 낮은 마찰 항력을 달성한다.

21.11.10 로봇 공학에서의 익형 선정

비행 로봇의 공기역학적 설계에서 익형의 선정은 임무 요구 사항과 운용 조건에 따라 체계적으로 이루어져야 한다. 소형 무인 비행 로봇의 익형 선정 시 고려하여야 할 주요 인자는 다음과 같다:

레이놀즈 수 영역: 소형 고정익 UAV는 일반적으로 Re = 10^4 \sim 5 \times 10^5 범위에서 운용되며, 이 영역에서 우수한 공력 성능을 나타내는 전용 저레이놀즈 수 익형을 선정하여야 한다. UIUC(University of Illinois at Urbana-Champaign) 저속 익형 데이터베이스는 이 영역의 다양한 익형에 대한 실험 데이터를 제공한다.

설계 양력 계수와 항력 극곡선: 순항 조건에서의 양력 계수 요구값을 만족하면서 최소 항력을 나타내는 익형을 선정한다. 항력 극곡선(drag polar, C_D vs. C_L)의 항력 버킷(drag bucket) 범위가 운용 양력 계수 범위를 포함하는 것이 바람직하다.

실속 특성: 비행 안전성을 위해 급격한 실속보다 점진적 실속 특성을 나타내는 익형이 선호된다. 특히 자율 비행 로봇의 경우, 급격한 실속은 제어 시스템의 안정성을 위협할 수 있다.

로터 블레이드 익형: 멀티로터 시스템의 로터 블레이드에는 넓은 받음각 범위에서 높은 양항비를 나타내며, 낮은 모멘트 계수를 갖는 익형이 적합하다. 블레이드의 반경 방향으로 서로 다른 익형을 배치하는 가변 익형(variable airfoil) 설계도 성능 최적화에 활용된다.

XFOIL(Drela, 1989)과 같은 점성-비점성 연성 해석 코드(viscous-inviscid coupled analysis code)는 저레이놀즈 수 영역에서의 익형 공력 특성을 효율적으로 예측할 수 있어, 비행 로봇의 익형 선정과 설계에 널리 활용되고 있다.

참고 문헌

  • Abbott, I. H., & von Doenhoff, A. E. (1959). Theory of Wing Sections: Including a Summary of Airfoil Data. Dover Publications.
  • Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill Education.
  • Drela, M. (1989). XFOIL: An analysis and design system for low Reynolds number airfoils. In Low Reynolds Number Aerodynamics, Lecture Notes in Engineering, Vol. 54, Springer.
  • Jacobs, E. N., Ward, K. E., & Pinkerton, R. M. (1933). The characteristics of 78 related airfoil sections from tests in the variable-density wind tunnel. NACA Report No. 460.
  • Selig, M. S., Guglielmo, J. J., Broeren, A. P., & Giguère, P. (1995). Summary of Low-Speed Airfoil Data, Vol. 1. SoarTech Publications.
  • Whitcomb, R. T. (1974). Review of NASA supercritical airfoils. ICAS Paper No. 74-10.

버전: 2026-04-16-v1.0