22.6 얇은 익형 이론(Thin Airfoil Theory)

1. 얇은 익형 이론의 기본 가정과 개요

얇은 익형 이론(Thin Airfoil Theory)은 2차원 비압축성 포텐셜 유동(potential flow) 내에서 두께가 시위(chord)에 비하여 충분히 작은 익형의 공력 특성을 해석적으로 구하는 이론적 체계이다. 이 이론은 Max Munk(1922)와 Hermann Glauert(1926)에 의하여 체계적으로 정립되었으며, 선형 공력학(linear aerodynamics)의 핵심 기반을 이룬다 (Glauert, 1926; Anderson, 2017).

얇은 익형 이론은 다음의 기본 가정에 입각한다.

1. 비압축성 유동(Incompressible Flow): 자유류 마하 수(Mach number)가 약 0.3 이하로서 밀도 변화를 무시할 수 있다.
2. 비점성 유동(Inviscid Flow): 유체의 점성 효과를 무시하며, 포텐셜 유동으로 취급한다.
3. 소교란 가정(Small Perturbation Assumption): 익형에 의한 유동의 교란이 자유류 속도에 비하여 충분히 작다.
4. 얇은 익형 가정: 익형의 최대 두께비(t/c)가 충분히 작아 익형 표면을 시위선(chord line) 또는 캠버선(camber line) 위에 근사할 수 있다.
5. 작은 받음각(Small Angle of Attack): 받음각 \alpha가 충분히 작아 \sin\alpha \approx \alpha, \cos\alpha \approx 1로 근사할 수 있다.

이러한 가정 하에서 익형 주위의 유동 문제는 선형 편미분방정식으로 환원되며, 중첩의 원리(principle of superposition)를 적용할 수 있게 된다.

2. 경계 조건의 선형화

얇은 익형에서 유동이 익형 표면을 따라 흐른다는 불투과 경계 조건(impermeability boundary condition)은 다음과 같이 표현된다. 익형 표면의 형상을 y = f(x)로 나타내면, 표면에서 법선 방향 속도 성분이 영(zero)이어야 한다.

소교란 가정을 적용하면, 전체 속도를 자유류 속도 U_\infty와 교란 속도 성분 (u', v')의 합으로 분해할 수 있다.

u = U_\infty + u', \quad v = v'

여기서 u' \ll U_\infty, v' \ll U_\infty이다. 이때 익형 표면의 경계 조건을 선형화하면 다음과 같다.

\frac{v'(x, 0)}{U_\infty} = \frac{df}{dx} = \frac{d\bar{y}}{dx} - \alpha

여기서 \bar{y}(x)는 캠버선의 형상 함수이고, \alpha는 받음각이다. 이 조건은 익형 표면이 아닌 시위선(y = 0) 위에서 적용되며, 이것이 얇은 익형 이론의 핵심적인 선형화이다.

3. 와류 분포를 이용한 모형화

얇은 익형 이론에서는 익형을 시위선 위에 분포된 와류 시트(vortex sheet)로 대체한다. 시위선 위의 위치 x에서 단위 길이당 와류 강도를 \gamma(x)로 정의하면, 위치 x_0에서의 와류에 의하여 위치 x에서 유도되는 수직 방향 속도(induced velocity)는 비오-사바르 법칙(Biot-Savart law)의 2차원 형태에 의하여 다음과 같다.

dv(x) = -\frac{\gamma(x_0)}{2\pi(x - x_0)} dx_0

전체 와류 시트에 의한 x 위치에서의 유도 수직 속도는 다음의 적분으로 주어진다.

v'(x, 0) = -\frac{1}{2\pi} \int_0^c \frac{\gamma(x_0)}{x - x_0} dx_0

여기서 c는 시위 길이이며, 적분은 코시 주치(Cauchy principal value)의 의미에서 수행된다.

4. 기본 적분 방정식

경계 조건과 와류 유도 속도를 결합하면 얇은 익형 이론의 기본 적분 방정식을 얻는다.

\frac{1}{2\pi} \int_0^c \frac{\gamma(x_0)}{x - x_0} dx_0 = U_\infty \left( \alpha - \frac{d\bar{y}}{dx} \right)

이 식은 미지 함수 \gamma(x_0)에 대한 특이 적분 방정식(singular integral equation)이다. 이를 풀기 위하여 쿠타 조건(Kutta condition)이 추가적인 경계 조건으로 부과된다. 쿠타 조건은 후연(trailing edge)에서 와류 강도가 영이 되어야 한다는 조건이다.

\gamma(c) = 0

5. 글라우어트 변환과 급수 해

이 적분 방정식의 해를 구하기 위하여 글라우어트 변환(Glauert transformation)을 도입한다. 변수 x를 각도 변수 \theta로 치환한다.

x = \frac{c}{2}(1 - \cos\theta), \quad 0 \leq \theta \leq \pi

이 변환에 의하여 x = 0은 \theta = 0(전연, leading edge)에, x = c는 \theta = \pi(후연, trailing edge)에 대응한다. 와류 강도를 다음과 같은 푸리에 급수(Fourier series)로 전개한다.

\gamma(\theta) = 2U_\infty \left[ A_0 \frac{1 + \cos\theta}{\sin\theta} + \sum_{n=1}^{\infty} A_n \sin(n\theta) \right]

첫 번째 항은 전연에서의 특이성(singularity)을 포함하며 쿠타 조건(\theta = \pi에서 \gamma = 0)을 자동으로 만족한다. 두 번째 항의 \sin(n\theta) 급수는 \theta = 0과 \theta = \pi 모두에서 영이 되어 정규적(regular)이다.

이 전개를 기본 적분 방정식에 대입하고, 글라우어트의 적분 공식(Glauert’s integral formula)을 적용하면 다음을 얻는다.

\frac{1}{\pi} \int_0^\pi \frac{\cos(n\theta_0)}{\cos\theta - \cos\theta_0} d\theta_0 = \frac{\sin(n\theta)}{\sin\theta}

이를 이용하면 경계 조건으로부터 각 계수를 결정할 수 있다.

A_0 = \alpha - \frac{1}{\pi} \int_0^\pi \frac{d\bar{y}}{dx} d\theta

A_n = \frac{2}{\pi} \int_0^\pi \frac{d\bar{y}}{dx} \cos(n\theta) \, d\theta, \quad n = 1, 2, 3, \ldots

여기서 d\bar{y}/dx는 \theta의 함수로 표현된 캠버선의 기울기이다.

6. 공력 계수의 도출

6.1 양력 계수

와류 분포가 결정되면, 쿠타-주코프스키 정리(Kutta-Joukowski theorem)에 의하여 단위 스팬당 양력은 전체 순환(circulation)에 비례한다.

L' = \rho_\infty U_\infty \Gamma

전체 순환은 다음과 같다.

\Gamma = \int_0^c \gamma(x) \, dx = \frac{c}{2} \int_0^\pi \gamma(\theta) \sin\theta \, d\theta

이를 계산하면 다음을 얻는다.

\Gamma = \pi c U_\infty \left( A_0 + \frac{A_1}{2} \right)

따라서 양력 계수(lift coefficient)는 다음과 같다.

C_l = \frac{L'}{\frac{1}{2}\rho_\infty U_\infty^2 c} = 2\pi \left( A_0 + \frac{A_1}{2} \right)

이를 전개하면 다음과 같이 표현할 수 있다.

C_l = 2\pi \left[ \alpha - \frac{1}{\pi}\int_0^\pi \frac{d\bar{y}}{dx} d\theta + \frac{1}{\pi}\int_0^\pi \frac{d\bar{y}}{dx}\cos\theta \, d\theta \right]

6.2 양력 기울기

대칭 익형(d\bar{y}/dx = 0)의 경우, A_0 = \alpha이고 A_n = 0이므로 양력 계수는 다음과 같이 단순화된다.

C_l = 2\pi\alpha

이로부터 양력 기울기(lift curve slope)는 다음과 같다.

\frac{dC_l}{d\alpha} = 2\pi \approx 6.283 \; \text{rad}^{-1}

이 값은 얇은 익형 이론의 가장 중요한 결과 중 하나로서, 2차원 비압축성 유동에서 익형 형상에 무관하게 양력 기울기가 2\pi임을 의미한다.

6.3 영양력 받음각

캠버가 있는 익형의 경우, 양력이 영이 되는 받음각(zero-lift angle of attack) \alpha_{L=0}는 C_l = 0 조건으로부터 다음과 같이 구해진다.

\alpha_{L=0} = \frac{1}{\pi}\int_0^\pi \frac{d\bar{y}}{dx}(1 - \cos\theta) \, d\theta \cdot (-1)

보다 명시적으로 표기하면 다음과 같다.

\alpha_{L=0} = -\frac{1}{\pi}\int_0^\pi \frac{d\bar{y}}{dx}(\cos\theta - 1) \, d\theta

캠버가 양(positive)인 익형에서는 \alpha_{L=0} < 0이 되어, 받음각이 영인 상태에서도 양력이 발생한다.

7. 모멘트 계수

7.1 전연 모멘트 계수

전연(leading edge)에 대한 피칭 모멘트 계수(pitching moment coefficient)는 와류 분포로부터 다음과 같이 계산된다.

C_{m,LE} = -\frac{\pi}{2}\left( A_0 + A_1 - \frac{A_2}{2} \right)

이를 양력 계수와 관련지으면 다음과 같다.

C_{m,LE} = -\frac{C_l}{4} - \frac{\pi}{4}(A_1 - A_2)

7.2 공력 중심과 1/4 시위점 모멘트

얇은 익형 이론에서 공력 중심(aerodynamic center), 즉 모멘트 계수가 받음각에 대하여 변하지 않는 점은 시위의 1/4 지점(x = c/4)에 위치한다. 이는 다음과 같이 확인할 수 있다. 1/4 시위점에 대한 모멘트 계수는 다음과 같다.

C_{m,c/4} = C_{m,LE} + \frac{C_l}{4} = -\frac{\pi}{4}(A_1 - A_2)

A_1과 A_2는 캠버선의 형상에만 의존하고 받음각 \alpha에 무관하므로, C_{m,c/4}는 받음각에 대하여 일정하다. 이것이 1/4 시위점이 공력 중심이 되는 이유이다.

대칭 익형의 경우 A_1 = A_2 = 0이므로 C_{m,c/4} = 0이다.

7.3 압력 중심

압력 중심(center of pressure) x_{cp}는 전체 모멘트가 영이 되는 점으로 다음과 같이 정의된다.

\frac{x_{cp}}{c} = -\frac{C_{m,LE}}{C_l} = \frac{1}{4} + \frac{\pi}{4C_l}(A_1 - A_2)

대칭 익형에서는 x_{cp}/c = 1/4로 일정하지만, 캠버 익형에서는 받음각에 따라 변한다.

8. 이론의 적용 범위와 한계

얇은 익형 이론은 다음의 조건에서 실제 공력 데이터와 양호한 일치를 보인다.

• 받음각이 약 \pm 8° 이내의 선형 영역
• 두께비 t/c가 약 12% 이하
• 레이놀즈 수(Reynolds number)가 충분히 높아 층류 박리(laminar separation)가 발생하지 않는 조건

반면 이 이론은 다음의 현상을 예측하지 못한다.

실제 익형의 양력 기울기는 점성 효과와 두께 효과에 의하여 이론값 2\pi보다 다소 낮은 값을 보이며, 일반적으로 실험에서는 약 5.7 \sim 6.0 \; \text{rad}^{-1} 범위가 관측된다 (Abbott & Von Doenhoff, 1959).

9. 드론 및 소형 항공기 설계에의 의의

얇은 익형 이론은 단순한 해석적 형태로 양력, 모멘트, 공력 중심 등 핵심 공력 특성을 제공하므로, 드론 및 소형 무인 항공기의 초기 설계 단계에서 익형 선정과 공력 성능 추정에 광범위하게 활용된다. 특히 양력 기울기 2\pi는 비행 제어 시스템의 공력 모델에서 기본 매개변수로 사용되며, 캠버에 의한 영양력 받음각의 변화는 익형의 설계 트림(trim) 조건 결정에 핵심적인 정보를 제공한다.

또한 이 이론은 패널법(panel method)이나 전산유체역학(CFD) 등 보다 정밀한 수치 해석 방법의 검증(validation)에도 기준 해(benchmark solution)로서 활용된다.

참고 문헌

• Abbott, I. H., & Von Doenhoff, A. E. (1959). Theory of Wing Sections. Dover Publications.
• Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill Education.
• Glauert, H. (1926). The Elements of Aerofoil and Airscrew Theory. Cambridge University Press.
• Katz, J., & Plotkin, A. (2001). Low-Speed Aerodynamics (2nd ed.). Cambridge University Press.
• Munk, M. M. (1922). General Theory of Thin Wing Sections. NACA Report, No. 142.

v 0.1.0