21.19 유한 날개 이론과 유도 항력 (Finite Wing Theory and Induced Drag)

1. 2차원 익형과 3차원 날개의 차이

2차원 익형(airfoil) 이론은 무한 스팬(infinite span)을 가정한 해석으로서, 스팬 방향의 유동 변화가 없는 이상적 조건을 전제한다. 그러나 실제 날개는 유한한 스팬을 가지며, 날개 끝(wing tip)에서 하면의 고압 영역과 상면의 저압 영역 사이의 압력 차이로 인해 날개 끝 주위의 횡방향 유동이 발생한다. 이 횡방향 유동은 날개 끝 와류(wing tip vortex)를 형성하며, 날개 후방으로 연장되는 후예 와류계(trailing vortex system)를 구성한다.

후예 와류계는 날개 주위의 유동장에 하향 유도 속도(induced downwash, $w$ )를 발생시키며, 이는 2차원 해석에서는 나타나지 않는 3차원 고유의 효과이다. 유도 하향류는 국소 자유류 방향을 하방으로 기울여 유효 받음각(effective angle of attack)을 감소시키고, 이에 따라 양력의 감소와 유도 항력(induced drag)의 발생을 초래한다.

2. 프란틀의 양력선 이론

프란틀의 양력선 이론(Prandtl’s lifting line theory, 1918)은 유한 스팬 날개의 공기역학을 체계적으로 해석하는 최초의 이론적 틀을 제공하였으며, 현재까지도 날개 설계의 기본 도구로 활용된다.

양력선 이론은 날개를 스팬 방향으로 배치된 단일 속박 와류(bound vortex)로 모델링한다. 속박 와류의 강도, 즉 순환 $\Gamma(y)$ 는 스팬 방향 위치 $y$ 의 함수이다. 헬름홀츠의 와류 정리(Helmholtz’s vortex theorem)에 의해 와류 강도가 변하는 곳에서는 자유 와류(free vortex)가 방출되어야 하므로, 순환의 스팬 방향 변화율 $d\Gamma/dy$ 에 비례하는 후예 와류가 날개로부터 하류로 방출된다.

후예 와류에 의한 스팬 방향 위치 $y_0$ 에서의 유도 하향류 속도는 비오-사바르 법칙(Biot-Savart law)으로부터 다음과 같이 계산된다:

$w(y_0) = -\frac{1}{4\pi}\int_{-b/2}^{b/2} \frac{d\Gamma/dy}{y_0 - y} \, dy$

여기서 $b$ 는 날개의 스팬이다. 이 적분은 코시 주치(Cauchy principal value)로 해석된다.

21.19.3 프란틀의 기본 적분 방정식

양력선 이론의 핵심은 각 스팬 방향 위치에서 국소 양력 계수가 2차원 익형 이론에 의해 유효 받음각의 함수로 결정된다는 가정이다. 유효 받음각은 기하학적 받음각에서 유도 받음각을 뺀 값이다:

$\alpha_\text{eff}(y_0) = \alpha(y_0) - \alpha_i(y_0)$

유도 받음각은:

$\alpha_i(y_0) = -\frac{w(y_0)}{V_\infty} = \frac{1}{4\pi V_\infty}\int_{-b/2}^{b/2} \frac{d\Gamma/dy}{y_0 - y} \, dy$

쿠타-주코프스키 정리에 의해 국소 양력은 $L'(y) = \rho_\infty V_\infty \Gamma(y)$ 이고, 2차원 양력 곡선으로부터:

$\Gamma(y_0) = \frac{1}{2} V_\infty c(y_0) a_0 [\alpha_\text{eff}(y_0) - \alpha_{L=0}(y_0)]$

이들을 결합하면 프란틀의 기본 적분 방정식이 도출된다:

$\alpha(y_0) = \frac{\Gamma(y_0)}{\pi V_\infty c(y_0)} + \alpha_{L=0}(y_0) + \frac{1}{4\pi V_\infty}\int_{-b/2}^{b/2} \frac{d\Gamma/dy}{y_0 - y} \, dy$

여기서 $a_0 = 2\pi$ 를 가정하였다. 이 적분 방정식을 순환 분포 $\Gamma(y)$ 에 대해 풀면 날개의 양력 분포, 유도 항력, 전체 양력을 결정할 수 있다.

21.19.4 타원 순환 분포와 최소 유도 항력

양력선 이론의 가장 중요한 결과 중 하나는 타원 순환 분포(elliptic circulation distribution)가 주어진 전체 양력에 대해 유도 항력을 최소화하는 최적 분포라는 것이다. 타원 순환 분포는 다음과 같이 표현된다:

$\Gamma(y) = \Gamma_0 \sqrt{1 - \left(\frac{2y}{b}\right)^2}$

여기서 $\Gamma_0$ 는 날개 중앙에서의 최대 순환이다.

타원 순환 분포의 특성:

유도 하향류 속도가 스팬 방향으로 균일하다: $w = -\frac{\Gamma_0}{2b} = \text{const.}$
유도 받음각이 스팬 방향으로 일정하다: $\alpha_i = \frac{C_L}{\pi AR}$
유도 항력 계수가 최소값을 가진다:

$C_{D,i} = \frac{C_L^2}{\pi AR}$

타원 순환 분포를 실현하기 위한 이상적 날개 평면형은 타원 평면형(elliptic planform)이지만, 시위 길이가 스팬 방향으로 변하는 비직사각형 평면형과 적절한 비틀림 분포의 조합에 의해서도 타원 순환 분포에 근사할 수 있다 (Anderson, 2017).

21.19.5 일반 순환 분포와 스팬 효율 인자

타원이 아닌 일반적인 순환 분포에서 유도 항력은 타원 분포보다 항상 크다. 일반적인 순환 분포를 푸리에 사인 급수로 전개하면:

$\Gamma(\theta) = 2bV_\infty \sum_{n=1}^{N} A_n \sin(n\theta)$

여기서 $y = -(b/2)\cos\theta$ 로 변수 치환하였다. 전체 양력 계수와 유도 항력 계수는 다음과 같이 표현된다:

$C_L = \pi AR \cdot A_1$

$C_{D,i} = \pi AR \sum_{n=1}^{N} n A_n^2 = \frac{C_L^2}{\pi AR}(1 + \delta)$

여기서 $\delta = \sum_{n=2}^{N} n(A_n/A_1)^2 \geq 0$ 이다. $\delta$ 는 순환 분포가 타원 분포로부터 벗어나는 정도를 나타내며, 항상 비음(nonnegative)이므로 유도 항력은 타원 분포의 최소값 이상이다.

스팬 효율 인자(span efficiency factor, $e_{\text{span}}$ )를 $e_{\text{span}} = 1/(1+\delta)$ 로 정의하면:

$C_{D,i} = \frac{C_L^2}{\pi e_{\text{span}} AR}$

21.19.6 유도 항력의 물리적 해석

유도 항력의 물리적 본질은 양력 발생에 수반되는 운동 에너지의 후류로의 이송이다. 날개가 양력을 발생시키면 공기에 하방 운동량을 부여하며, 이 하방 운동을 갖는 공기는 운동 에너지를 보유한다. 이 운동 에너지는 날개의 추진 에너지로부터 공급되어야 하므로, 유도 항력으로 나타난다.

후예 와류계의 운동 에너지를 직접 계산하면 유도 항력의 크기를 확인할 수 있다. 트레프츠 평면(Trefftz plane, 날개로부터 충분히 하류의 평면)에서 후예 와류계가 유도하는 유동의 운동 에너지 유속(flux)은 유도 항력 $D_i$ 에 자유류 속도 $V_\infty$ 를 곱한 값, 즉 유도 항력에 의한 동력 손실과 같다:

$D_i V_\infty = \frac{1}{2}\rho_\infty \int\int (v^2 + w^2) V_\infty \, dA$

여기서 $v$ 와 $w$ 는 트레프츠 평면에서의 횡방향 및 수직 방향 유도 속도 성분이다.

3. 유도 항력 저감 설계

유도 항력을 저감하기 위한 주요 설계 전략은 다음과 같다:

높은 종횡비: 유도 항력 계수가 종횡비에 반비례하므로, 종횡비의 증가는 유도 항력 저감에 가장 직접적인 방법이다. 그러나 구조 중량의 증가, 날개 강성 저하, 기동성 감소 등의 trade-off가 존재한다.

타원 양력 분포의 구현: 스팬 방향 양력 분포를 타원에 가깝게 설계하여 스팬 효율 인자를 1에 근접시킨다. 적절한 테이퍼비( $\lambda \approx 0.4 \sim 0.5$ )와 날개 비틀림의 조합으로 타원 분포에 근사할 수 있다.

날개 끝 장치(winglet): 날개 끝에 수직 또는 경사진 소형 날개면을 부착하여 날개 끝 와류의 강도를 약화시키고 유도 항력을 저감한다. Whitcomb(1976)이 제안한 윙렛은 유효 스팬을 증가시키는 효과를 나타내며, 3~7%의 유도 항력 저감이 보고되었다.

비평면 날개(nonplanar wing): 날개를 상반각(dihedral)이나 비평면 형태로 설계하면 후예 와류계의 구조가 변화하여 유도 항력이 감소할 수 있다. 뮌크(Munk)의 스태거 정리(stagger theorem)에 의하면, 유도 항력은 양력면의 전후 방향 배치에 무관하고 전면(Trefftz plane)에서의 투영 형태에만 의존한다.

4. 수치 양력선 방법

프란틀의 적분 방정식을 수치적으로 해석하기 위해 다양한 수치 양력선 방법(numerical lifting line method)이 개발되었다. 가장 기본적인 방법은 순환 분포를 유한 개의 푸리에 사인 급수 항으로 전개하고, 스팬 방향의 이산 제어점에서 프란틀 방정식을 만족시키는 연립 선형 방정식을 구성하여 계수를 결정하는 것이다.

와류 격자법(Vortex Lattice Method, VLM)은 양력선 이론을 3차원으로 확장한 수치적 방법으로서, 날개 표면을 패널(panel)로 분할하고 각 패널에 말발굽 와류(horseshoe vortex)를 배치한다. VLM은 임의 평면형, 후퇴각(sweep), 상반각을 갖는 날개의 양력 분포와 유도 항력을 효율적으로 계산하며, 다중 양력면(multiple lifting surface) 구성에 대한 해석도 가능하다 (Katz & Plotkin, 2001).

5. 양력선 이론의 한계와 확장

프란틀의 양력선 이론은 다음의 가정에 기반하므로 적용에 한계가 존재한다:

날개를 단일 양력선으로 근사하므로, 시위 방향의 압력 분포를 해석할 수 없다.
큰 후퇴각(sweep angle)을 갖는 날개에는 적용이 부적절하다.
낮은 종횡비의 날개에서 정확도가 저하된다.
비선형 공력(큰 받음각, 실속 영역)을 직접 처리하지 못한다.

이러한 한계를 극복하기 위해 양력면 이론(lifting surface theory)과 와류 격자법이 개발되었으며, 더 나아가 패널 방법과 전산 유체역학(CFD)이 임의 형상의 3차원 공력 해석에 활용된다.

비선형 양력선 방법(nonlinear lifting line method)은 각 스팬 위치에서의 국소 양력을 실험적 익형 데이터(실속 포함)로부터 결정하여, 실속 영역의 날개 공력 특성을 근사적으로 예측하는 확장된 방법이다.

6. 로봇 공학에서의 유한 날개 이론 적용

고정익 비행 로봇의 날개 설계에서 유한 날개 이론과 유도 항력의 이해는 필수적이다.

소형 고정익 UAV는 종횡비 $AR = 6 \sim 12$ 범위의 날개를 사용하는 경우가 많으며, 이 범위에서 유도 항력은 전체 항력의 상당 부분(순항 조건에서 30~40%)을 차지한다. 따라서 양력 분포의 최적화와 날개 평면형의 적절한 설계는 비행 효율 향상에 직접적으로 기여한다.

와류 격자법은 소형 UAV의 설계 초기 단계에서 날개 형상의 공력적 평가를 위한 효율적인 도구이며, 날개-꼬리 날개 간의 공기역학적 간섭 해석에도 활용된다. 멀티로터 시스템에서 로터 블레이드의 양력선 이론적 해석은 블레이드 요소-운동량 이론(Blade Element Momentum Theory, BEMT)의 일부를 구성한다.

편대 비행(formation flight)에서 선행 비행체의 날개 끝 와류는 후행 비행체에 상향류(upwash)를 유도하여 후행 비행체의 유도 항력을 감소시킬 수 있다. 이 원리는 다중 UAV의 편대 비행에 의한 에너지 절약 전략의 공기역학적 기초이다 (Lissaman & Shollenberger, 1970).

참고 문헌

Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill Education.
Katz, J., & Plotkin, A. (2001). Low-Speed Aerodynamics (2nd ed.). Cambridge University Press.
Lissaman, P. B. S., & Shollenberger, C. A. (1970). Formation flight of birds. Science, 168(3934), 1003–1005.
Whitcomb, R. T. (1976). A design approach and selected wind-tunnel results at high subsonic speeds for wing-tip mounted winglets. NASA TN D-8260.

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