양력선 이론

유체역학 전산유체역학 마지막 수정: 2026-08-22 04:51:33

1. 개요[편집]

양력선 이론 Lifting-line theory
정립	Ludwig Prandtl, 괴팅겐 (1918~1919)
별칭	랑체스터-프란틀 날개 이론
모형	스팬 방향 결합와 + 후류 와사슬
핵심 결과	유도항력 계수가 양력계수의 제곱에 비례
최적해	타원 순환 분포 (스팬효율 e = 1)
유효 범위	높은 종횡비 · 비압축 · 부착유동 · 무후퇴각

2차원 익형에는 항력이 없다는 달랑베르의 역설이 있었다. 날개에 끝을 달아 준 순간 항력이 생겼고, 그 항력이 어디서 오는지를 설명한 것이 이 이론이다.

양력선 이론(lifting-line theory)은 유한한 스팬을 가진 날개를 스팬 방향으로 놓인 하나의 결합와선(bound vortex line)과 그로부터 하류로 뻗는 자유와 시트로 치환해, 순환 분포 $\Gamma(y)$ 를 미지수로 하는 적분미분방정식을 세우고 푸는 고전 공기역학 이론이다. 1918~1919년 루트비히 프란틀이 괴팅겐에서 정립했고, 랑체스터가 1907년에 정성적으로 먼저 그림을 그렸다는 이유로 랑체스터-프란틀 이론이라고도 부른다.

이 이론이 준 것은 계산 도구 하나가 아니라 개념 하나다. 3차원 날개의 항력 중 점성과 무관한 성분, 즉 유도항력이 존재하며 그것이 후류에 남긴 운동에너지의 대가라는 사실. 공력해석에서 아무렇지 않게 쓰는 $C_{D,i} = C_L^2/(\pi e\,AR)$ 이라는 한 줄이 통째로 여기서 나온다. 이 문서는 그 한 줄이 어떻게 유도되고, 어디까지만 참인가를 다룬다.

2. 결합와와 말굽와 — 헬름홀츠가 강제한 구조[편집]

출발점은 쿠타-주코프스키 정리다. 2차원 익형의 단위 스팬당 양력은 순환에 비례한다.

L' = \rho_\infty V_\infty \Gamma

3차원 날개에서는 시위와 비틀림이 스팬 방향으로 변하므로 순환도 $\Gamma(y)$ 로 변한다. 그런데 순환이 변하면 그냥 넘어갈 수가 없다. 헬름홀츠의 와정리(1858)가 두 가지를 강제하기 때문이다.

와선(vortex filament)의 세기는 그 선을 따라 일정하다.
와선은 유체 내부에서 끝날 수 없다. 닫힌 고리를 이루거나, 무한히 뻗거나, 경계에서 끝나야 한다.

그러므로 $y$ 에서 $y+dy$ 로 가며 결합와의 세기가 $d\Gamma$ 만큼 줄었다면, 그 차이만큼의 와선이 날개를 떠나 하류로 흘러나가야 한다. 이것이 후류 와사슬(trailing vortex sheet)이고, 시트의 국소 세기는

\gamma(y) = -\frac{d\Gamma}{dy}

이다. 순환이 급격히 변하는 날개끝 근처에서 시트가 가장 강하고, 이 시트는 하류에서 스스로 말려 우리가 아는 날개끝 와류 한 쌍이 된다. 결합와 + 좌우 자유와의 최소 단위를 말굽와(horseshoe vortex)라 부르며, 이륙 순간 뒤에 남겨진 출발와(starting vortex)가 저 멀리서 고리를 닫아 준다. 켈빈의 순환 정리가 요구하는 대로 총 순환은 0으로 유지된다.¹

3. 후류가 만드는 하강기류[편집]

여기서 이론의 핵심 논리가 나온다. 하류로 뻗은 와사슬은 비오-사바르 법칙에 따라 날개 자신의 위치에도 유도속도를 만든다. 반무한 와선 하나가 자기 시작점 평면에서 만드는 유도속도는 무한 와선의 절반이고, 시트 전체를 적분하면 $y_0$ 지점의 하강기류는

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

가 된다( $b$ 는 스팬, 아래 방향을 음으로 잡은 관례). 이 하강기류가 국소 유동 방향을 아래로 꺾으므로, 날개 단면이 실제로 느끼는 받음각은 기하학적 받음각보다 작다. 그 차이가 유도받음각이다.

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

결과는 두 가지다. 첫째, 유효 받음각이 깎여 양력이 준다( $\alpha_{\text{eff}} = \alpha - \alpha_i$ ). 둘째, 국소 양력 벡터가 국소 유동에 수직이므로 자유류 기준으로는 뒤로 $\alpha_i$ 만큼 기울고, 그 뒤로 기운 성분이 자유류 방향의 힘 — 즉 항력이 된다.

D_i' = L'\sin\alpha_i \approx L'\,\alpha_i

점성이 전혀 없는데 항력이 나왔다. 달랑베르의 역설이 3차원에서 깨지는 지점이며, 이 항력의 정체는 후류에 영원히 남겨진 회전 운동에너지다. 공짜 양력은 없다.

4. 프란틀의 적분미분방정식[편집]

이제 두 개의 표현을 이어 붙인다. 국소 단면의 양력은 얇은 익형 이론으로 $c_l = a_0(\alpha_{\text{eff}} - \alpha_{L=0})$ 이고( $a_0 \approx 2\pi$ ), 동시에 쿠타-주코프스키로 $c_l = 2\Gamma/(V_\infty c)$ 다. 둘을 같다고 놓고 $\alpha_{\text{eff}} = \alpha - \alpha_i$ 를 대입하면 미지 함수 $\Gamma(y)$ 에 대한 하나의 방정식이 나온다.

\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

프란틀 양력선 이론의 기본 방정식이다. 좌변은 설계자가 정하는 것(스팬 방향 받음각 분포, 즉 기하학적 비틀림), 우변 첫 항은 시위 분포(테이퍼), 셋째 항은 후류가 자기 자신에게 되먹이는 항이다. 미지 함수가 적분 안팎에 동시에 들어 있는 적분미분방정식이고, 적분은 $y = y_0$ 에서 특이점을 가지므로 코시 주치(principal value)로 읽어야 한다.

5. 푸리에 급수 해[편집]

프란틀의 해법은 우아하다. 스팬 좌표를 각으로 바꾸고

y = -\frac{b}{2}\cos\theta, \qquad \theta \in [0, \pi]

순환을 사인 급수로 전개한다.

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

이 형태는 양 끝( $\theta = 0, \pi$ )에서 자동으로 $\Gamma = 0$ 이 되어 날개끝 조건을 공짜로 만족한다. 여기에 글라우어트 적분

\int_0^\pi \frac{\cos n\theta}{\cos\theta - \cos\theta_0}\,d\theta = \frac{\pi \sin n\theta_0}{\sin\theta_0}

를 쓰면 특이적분이 깔끔하게 닫히고, 유도받음각이 급수로 나온다.

\alpha_i(\theta_0) = \sum_{n=1}^{\infty} n A_n \frac{\sin n\theta_0}{\sin\theta_0}

기본 방정식을 $N$ 개의 스팬 위치에서 배치(collocation)하면 $A_1 \dots A_N$ 에 대한 선형 연립방정식이 되고, 이걸 풀면 끝이다. 손계산으로도 $N = 4$ 면 실용적인 정확도가 나왔다는 것이 1920년대에 이 이론이 압도적이었던 이유다.

계수를 알면 성능이 바로 나온다. 양력계수는 첫 항 하나에만 의존한다.

C_L = \pi\,AR\,A_1

여기서 $AR = b^2/S$ 는 종횡비다. 그리고 유도항력계수는 모든 항이 기여한다.

C_{D,i} = \pi\,AR \sum_{n=1}^{\infty} n A_n^2 = \frac{C_L^2}{\pi\,AR}\left[1 + \sum_{n=2}^{\infty} n\left(\frac{A_n}{A_1}\right)^2\right] = \frac{C_L^2}{\pi\,AR}\left(1 + \delta\right)

대괄호 안의 보정항 $\delta = \sum_{n\ge2} n (A_n/A_1)^2$ 은 제곱의 합이라 항상 0 이상이다. 이것을 흡수해 정의한 것이 스팬효율계수다.

e = \frac{1}{1+\delta} \le 1 \qquad \Longrightarrow \qquad C_{D,i} = \frac{C_L^2}{\pi\, e\, AR}

부호 방향을 헷갈리기 쉬운데, 관계는 $e = 1/(1+\delta)$ 이지 $e = 1-\delta$ 가 아니다. $\delta$ 가 커질수록 $e$ 는 작아지고 유도항력은 늘어난다.²

6. 타원 분포가 최적인 이유[편집]

$\delta \ge 0$ 이고 등호는 $A_2 = A_3 = \cdots = 0$ 일 때만 성립한다. 즉 1항만 남은 순환 분포가 유일한 최적이다. 그 분포는

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

즉 타원 분포다. 증명이라기엔 민망할 만큼 짧다 — $C_L$ 은 $A_1$ 만 보는데 $C_{D,i}$ 는 모든 $A_n^2$ 을 더하니, 같은 양력을 내면서 항력을 줄이려면 $n\ge2$ 를 전부 죽이는 수밖에 없다. 무크(Munk)가 변분법으로 일반화한 결과가 같은 결론에 도달한다.

타원 분포일 때 유도량들이 스팬 방향으로 일정해진다는 것이 물리적 핵심이다.

w = \frac{\Gamma_0}{2b} = \text{const}, \qquad \alpha_i = \frac{C_L}{\pi\,AR} = \text{const}, \qquad C_{D,i} = \frac{C_L^2}{\pi\,AR}

후류에 남기는 운동에너지를 스팬 전체에 균등하게 분산시키는 것이 최소 손실이라는, 제트 추진 효율의 논리와 같은 구조다. 유한 날개의 양력 기울기가 2차원보다 작아지는 것도 여기서 나온다.

a = \frac{dC_L}{d\alpha} = \frac{a_0}{1 + a_0/(\pi\,AR)} \quad (\text{타원 날개})

$a_0 = 2\pi$ 를 넣으면 $AR = 6$ 짜리 날개의 양력 기울기는 2차원 값의 75%로 떨어진다. 글라이더의 날개가 길고 가느다란 이유, 그리고 델타익 전투기가 큰 받음각을 써야 하는 이유가 한 식에 다 들어 있다.

여기서 자주 오해되는 지점이 둘 있다.

타원 평면형이 아니라 타원 “순환 분포”가 최적이다. 무비틀림·일정 익형이면 타원 시위 분포가 타원 순환을 만들지만, 비틀림과 익형 변화로도 얼마든지 만들 수 있다. 현대 여객기가 사다리꼴 평면형에 워시아웃(끝단 비틀림 내림)을 주는 이유가 이것 — 제작이 훨씬 싸면서 $e = 0.95$ 안팎이 나온다.
“타원이 항상 최적”은 제약 조건이 스팬일 때만 참이다. 제약을 날개 뿌리 굽힘 모멘트로 바꾸면 최적 분포가 안쪽에 몰린 종 모양(bell-shaped)으로 바뀐다. 프란틀 본인이 1933년 논문에서 이미 지적한 사실이고, 호르텐 형제의 전익기와 NASA의 Prandtl-D 실험기가 이 분포를 실제로 썼다. 같은 구조 중량이면 종 모양 쪽이 스팬을 더 길게 가져갈 수 있어 총 유도항력이 오히려 작아진다.³

7. 무엇이 성립하지 않는가[편집]

양력선 이론은 가정이 명시적이라 한계도 명시적이다.

높은 종횡비. 결합와를 하나의 선으로 뭉갠 것이므로 시위 방향의 압력 분포는 표현되지 않는다. $AR \gtrsim 4$ 에서 실용적이고, 그 아래로는 오차가 급격히 커진다. 델타익이나 미사일 핀에는 못 쓴다.
무후퇴각. 후퇴각이 붙으면 결합와가 꺾여 각 단면이 자기 자신에게도 유도속도를 만들고, “국소 단면은 2차원 익형처럼 군다”는 전제가 깨진다.
비압축·부착유동. 아음속 압축성은 프란틀-글라우어트 보정으로 얹을 수 있지만 천음속은 불가. 실속 이후는 얇은 익형 이론이 죽으므로 당연히 불가(다만 스팬 방향으로 $c_l$ 분포를 뽑아 어디서 먼저 실속할지 예측하는 데는 여전히 쓴다 — 워시아웃 설계의 근거다).
끝단 특이성. 타원 분포는 날개끝에서 $d\Gamma/dy \to \infty$ 라 유도속도 적분이 병적으로 민감하다. 수치해에서 배치점을 끝단에 너무 붙이면 진동이 생긴다.
평면 후류 가정. 후류 시트가 자유류 방향으로 곧게 흘러간다고 놓았는데, 실제로는 시트가 스스로 말려 올라간다. 고양력 상태에서 무시 못 할 오차가 된다.

또 하나, 유도항력은 점성 항력이 아니다. 실제 항력은 여기에 마찰항력과 형상항력이 더해지며, 게다가 받음각이 커지면 점성 항력도 $C_L^2$ 에 비례하는 성분을 갖는다. 그래서 항공기 전체 극곡선 $C_D = C_{D,0} + C_L^2/(\pi e_0 AR)$ 에 등장하는 오스왈드 효율계수 $e_0$ 는 양력선 이론의 스팬효율 $e$ 와 다른 물건이다. 잘 만든 날개의 $e$ 는 0.9~~0.99인데 전기체 $e_0$ 는 흔히 0.7~~0.85다. 그 차이는 동체 간섭과 양력 의존 점성 항력이 먹은 것이다.

8. 후속 — 양력면에서 CFD까지[편집]

양력선의 한계를 하나씩 걷어 낸 계보가 그대로 전산공력해석의 역사다.

양력면 이론(lifting-surface theory): 결합와를 선이 아니라 시위 방향으로도 퍼진 와도 시트로 놓는다. 저종횡비와 후퇴익을 다룰 수 있게 된다. 바이싱거의 준양력면법이 중간 단계.
와류격자법(vortex lattice method, VLM): 양력면을 말굽와 격자로 이산화하고 각 패널 제어점에서 유동 접선 조건을 부과해 선형 연립방정식을 푼다. 사실상 양력선 이론의 다차원 이산 버전이며, 오늘날에도 초기 설계 단계의 주력이다. 수백 개 패널이면 노트북에서 순식간에 돌아가는데 정확도는 부착유동 한정으로 CFD와 몇 % 안에서 붙는다.
패널법: 물체 표면 전체에 소스·더블릿 특이점을 분포시켜 라플라스 방정식을 경계요소법 방식으로 푼다. 두께 효과와 동체 간섭까지 들어온다. 포텐셜 유동 가정은 여전하다.
오일러/RANS CFD: 점성과 압축성, 박리까지. 대신 격자 하나 만드는 데 걸리는 시간이 VLM 해석 천 번보다 길다.

계보 전체가 상위 기법으로 대체된 것처럼 보이지만 현실은 다르다. 개념 설계 단계에서 스팬과 비틀림을 정하는 일은 여전히 양력선·VLM의 몫이고, 최적화 루프를 수천 번 돌려야 하는 다분야 설계 최적화에서는 값싼 저차 모형이 없으면 아예 시작이 안 된다. CFD는 결정을 검증하고, 양력선 이론은 결정을 만든다.

9. 관련 문서[편집]

그래서 원리적으로는 비행기가 이륙한 활주로 끝에 출발와가 아직도 떠 있다. 물론 점성이 몇 분 안에 잡아먹지만, 대형기의 날개끝 와류가 후속기에 위험할 만큼 오래 남는 것은 같은 이유다. 관제에서 항공기 사이 간격을 기종 무게별로 벌리는 규정이 전부 이 이론의 실무 버전이다. ↩
시험에서 가장 흔한 실수가 $e$ 와 $\delta$ 를 뒤집는 것이다. 외우기 좋은 방향은 이렇다 — $\delta$ 는 “타원에서 벗어난 정도”라 작을수록 좋고, $e$ 는 “효율”이라 클수록 좋다. 둘은 역수 관계로 묶여 있고 타원에서 각각 0과 1이다. ↩
새가 V자 편대로 나는 이유도 이 그림 안에 있다. 앞 새의 날개끝 와류는 바깥쪽에서 상승기류를 만들므로, 뒤 새가 그 자리에 날개끝을 놓으면 유도받음각이 줄어 유도항력이 덜 든다. 2014년 붉은볼따오기 편대를 GPS와 가속도계로 측정한 연구에서, 새들이 상승기류 위치를 정확히 찾을 뿐 아니라 날갯짓 위상까지 맞춘다는 것이 확인됐다. 윙렛도 같은 계산의 다른 답인데, 흔한 오해와 달리 “와류를 막는 판때기”가 아니라 와류가 만든 유동장 속에서 안쪽으로 양력을 내고 그 힘의 전방 성분으로 항력을 상쇄하는 작은 날개다. ↩