22.7 대칭 익형의 양력 해석

1. 대칭 익형의 정의와 기하학적 특성

대칭 익형(symmetric airfoil)은 시위선(chord line)에 대하여 상면과 하면의 형상이 동일한 익형을 의미한다. 즉, 캠버선(camber line)이 시위선과 일치하여 최대 캠버(maximum camber)가 영(zero)인 형상이다. 수학적으로 캠버선의 형상 함수 \bar{y}(x) = 0이 전 시위 영역에서 성립한다.

대표적인 대칭 익형으로는 NACA 0006, NACA 0009, NACA 0012, NACA 0015, NACA 0018 등의 NACA 4자리 계열 익형이 있다. 여기서 앞 두 자리 ’00’은 캠버가 없음을 나타내며, 뒤 두 자리는 최대 두께비(%)를 의미한다. 예를 들어 NACA 0012는 시위 길이 대비 최대 두께가 12%인 대칭 익형이다 (Abbott & Von Doenhoff, 1959).

대칭 익형의 두께 분포(thickness distribution)는 NACA 4자리 계열의 경우 다음과 같이 주어진다.

\frac{y_t}{c} = \frac{t}{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]

여기서 y_t는 시위선으로부터의 반두께(half-thickness), t는 최대 두께비, c는 시위 길이이다.

2. 얇은 익형 이론에 의한 대칭 익형 해석

2.1 경계 조건의 단순화

대칭 익형에서는 캠버선이 시위선과 일치하므로 캠버선의 기울기가 전 영역에서 영이다.

\frac{d\bar{y}}{dx} = 0

이를 얇은 익형 이론의 기본 적분 방정식에 대입하면 다음과 같이 단순화된다.

\frac{1}{2\pi} \int_0^c \frac{\gamma(x_0)}{x - x_0} dx_0 = U_\infty \alpha

여기서 \gamma(x_0)는 시위선 위 위치 x_0에서의 와류 강도(vortex strength), U_\infty는 자유류 속도, \alpha는 받음각이다. 우변이 상수(U_\infty \alpha)가 되어 적분 방정식의 해가 크게 간결해진다.

2.2 글라우어트 변환과 와류 분포

글라우어트 변환(Glauert transformation) x = \frac{c}{2}(1 - \cos\theta)를 적용하면, 와류 강도의 푸리에 급수 전개에서 계수들이 다음과 같이 결정된다.

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

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

d\bar{y}/dx = 0이므로 모든 고차 계수가 영이 된다. 따라서 대칭 익형의 와류 분포는 다음과 같이 단일 항으로 표현된다.

\gamma(\theta) = 2U_\infty \alpha \frac{1 + \cos\theta}{\sin\theta}

이를 x 좌표로 다시 표현하면 다음과 같다.

\gamma(x) = 2U_\infty \alpha \sqrt{\frac{c - x}{x}}

이 와류 분포는 전연(leading edge, x = 0)에서 \gamma \to \infty로 특이성(singularity)을 가지며, 후연(trailing edge, x = c)에서 \gamma = 0으로 쿠타 조건(Kutta condition)을 만족한다.

2.3 전연 특이성의 물리적 의미

전연에서의 와류 강도 특이성은 실제 유동에서 전연 부근의 급격한 유동 가속과 강한 흡입 압력 피크(suction peak)에 대응한다. 실제 익형에서는 유한한 전연 반지름(leading edge radius)이 존재하여 무한대의 속도가 나타나지 않지만, 얇은 익형 이론은 두께를 무시하므로 이러한 특이성이 해에 내재한다. 전연 반지름이 작을수록 실제 유동에서의 흡입 피크도 더 날카로워지며, 이는 전연 실속(leading edge stall) 발생의 주요 인자가 된다 (Anderson, 2017).

3. 양력 계수의 도출

3.1 순환과 양력

대칭 익형의 전체 순환(circulation)은 와류 분포를 시위 전체에 걸쳐 적분하여 구한다.

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

\gamma(\theta) = 2U_\infty\alpha\frac{1+\cos\theta}{\sin\theta}를 대입하면 다음을 얻는다.

\Gamma = U_\infty c \alpha \int_0^\pi (1 + \cos\theta) \, d\theta = \pi c U_\infty \alpha

쿠타-주코프스키 정리(Kutta-Joukowski theorem)에 의하여 단위 스팬당 양력은 다음과 같다.

L' = \rho_\infty U_\infty \Gamma = \pi \rho_\infty U_\infty^2 c \alpha

3.2 양력 계수

양력 계수(lift coefficient)는 다음과 같이 정의된다.

C_l = \frac{L'}{\frac{1}{2}\rho_\infty U_\infty^2 c}

이에 대입하면 대칭 익형의 양력 계수는 다음과 같다.

C_l = 2\pi\alpha

이 결과는 얇은 익형 이론의 가장 기본적이고 중요한 결론이다. 양력 계수가 받음각에 대하여 선형적으로 비례하며, 그 기울기(lift curve slope)가 다음과 같음을 나타낸다.

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

3.3 영양력 받음각

대칭 익형에서 양력이 영이 되는 받음각(zero-lift angle of attack)은 다음과 같다.

\alpha_{L=0} = 0

이는 캠버가 없으므로 받음각이 영일 때 상면과 하면의 압력 분포가 대칭이 되어 양력이 발생하지 않기 때문이다.

4. 압력 분포

4.1 압력 계수

대칭 익형 주위의 압력 계수(pressure coefficient) 분포는 와류 강도로부터 유도할 수 있다. 익형 상면과 하면의 압력 차이는 다음과 같다.

\Delta C_p(x) = C_{p,\text{lower}} - C_{p,\text{upper}} = \frac{2\gamma(x)}{U_\infty}

대칭 익형의 경우 이를 대입하면 다음과 같다.

\Delta C_p(x) = 4\alpha\sqrt{\frac{c - x}{x}}

이 분포는 전연(x = 0) 부근에서 매우 큰 값을 가지며, 후연(x = c)을 향하여 단조감소하여 영에 도달한다. 이는 대칭 익형의 양력이 전연 부근에 집중되어 발생함을 의미한다.

4.2 압력 분포와 양력의 관계

압력 차이를 시위 전체에 걸쳐 적분하면 양력 계수를 다시 얻을 수 있다.

C_l = \frac{1}{c}\int_0^c \Delta C_p(x) \, dx = \frac{4\alpha}{c}\int_0^c \sqrt{\frac{c-x}{x}} \, dx = 2\pi\alpha

이는 와류 분포로부터 구한 결과와 일치하며, 해석의 자기 일관성(self-consistency)을 확인시켜 준다.

5. 모멘트 특성

5.1 전연 모멘트 계수

대칭 익형의 전연에 대한 피칭 모멘트 계수(pitching moment coefficient)는 다음과 같다.

C_{m,LE} = -\frac{1}{c^2}\int_0^c x \, \Delta C_p(x) \, dx

와류 분포를 이용하면 다음을 얻는다.

C_{m,LE} = -\frac{\pi}{2}A_0 = -\frac{\pi\alpha}{2} = -\frac{C_l}{4}

5.2 /4 시위점 모멘트 계수

시위의 1/4 지점(x = c/4)에 대한 모멘트 계수는 다음과 같다.

C_{m,c/4} = C_{m,LE} + \frac{C_l}{4} = -\frac{C_l}{4} + \frac{C_l}{4} = 0

대칭 익형에서는 1/4 시위점에 대한 모멘트가 받음각에 관계없이 항상 영이다. 이는 1/4 시위점이 공력 중심(aerodynamic center)인 동시에 압력 중심(center of pressure)이 됨을 의미한다.

5.3 압력 중심

대칭 익형의 압력 중심 위치는 다음과 같이 구해진다.

\frac{x_{cp}}{c} = -\frac{C_{m,LE}}{C_l} = \frac{1}{4}

캠버가 없는 대칭 익형에서는 압력 중심이 받음각에 무관하게 항상 시위의 1/4 지점에 고정된다. 이 특성은 비행 안정성 해석에서 중요한 의미를 가진다.

6. 이론과 실험의 비교

6.1 양력 기울기의 실험적 검증

얇은 익형 이론이 예측하는 양력 기울기 2\pi \; \text{rad}^{-1}은 실제 대칭 익형의 풍동 실험 데이터와 비교할 때 약간의 차이를 보인다. 실험적으로 관측되는 양력 기울기는 일반적으로 이론값보다 다소 낮다.

출처: Abbott & Von Doenhoff (1959)

두께가 증가할수록 실험값이 이론값에서 더 벗어나는 경향이 관찰된다. 이는 두께 효과에 의한 속도 분포의 변화와 점성 효과(경계층 성장)에 기인한다.

6.2 두께 보정

유한한 두께를 고려한 보정된 양력 기울기는 다음과 같이 표현할 수 있다.

a_0 \approx 2\pi\left(1 + 0.77\frac{t}{c}\right)

여기서 t/c는 최대 두께비이다. 이 보정은 두께에 의한 속도 증가 효과를 반영하며, 포텐셜 유동 해석에 기반한 것이다. 그러나 점성 효과까지 포함하면 실제 양력 기울기는 이 보정값보다도 다소 낮아지게 된다 (Katz & Plotkin, 2001).

6.3 점성 효과에 의한 보정

점성 유동에서 경계층의 배제 두께(displacement thickness) \delta^*는 익형의 유효 형상을 변화시킨다. 이에 따라 실효 캠버와 실효 받음각이 변하며, 양력 기울기를 감소시키는 효과를 가져온다. 이를 고려한 양력 기울기의 근사적 보정은 다음과 같다.

a_0^{\text{viscous}} \approx \frac{2\pi}{1 + 2\delta^*/c}

여기서 \delta^*/c는 무차원 배제 두께이다. 일반적인 비행 레이놀즈 수에서 \delta^*/c는 매우 작으므로 점성 보정의 크기도 제한적이지만, 저레이놀즈 수 영역에서는 그 영향이 현저하게 증가한다.

7. 대칭 익형의 공학적 활용

대칭 익형은 받음각이 영일 때 양력이 영이고, 압력 중심이 1/4 시위점에 고정되는 특성으로 인하여 다음과 같은 응용 분야에서 널리 사용된다.

1. 수직 안정판 및 수평 안정판: 항공기의 미익(tail surface)에서 트림 조건의 설정이 용이하다.
2. 헬리콥터 및 멀티로터 블레이드: 블레이드가 회전하면서 양방향 유동을 경험하는 조건에서 대칭 형상이 유리하다.
3. 가동 조종면(control surface): 방향타(rudder), 승강타(elevator) 등의 조종면에서 중립 위치의 힌지 모멘트가 영이 되는 장점이 있다.
4. 풍동 실험 기준 모형: 이론 해석 결과와의 비교가 용이하여 실험 검증용 표준 모형으로 활용된다.

참고 문헌

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

v 0.1.0