쿠타-주코프스키 정리

1. 개요[편집]

쿠타-주코프스키 정리 Kutta–Joukowski theorem
진술	단위 스팬 양력 L′ = ρ∞ V∞ Γ
힘의 방향	다가오는 유동에 수직 (항력 성분 0)
전제	비점성 · 비압축 · 2차원 · 정상 · 물체 밖 비회전
Γ 를 정하는 것	정리 자신이 아니라 쿠타 조건
짝 정리	달랑베르 역설 (같은 계산의 실수부)
이름	M. W. Kutta (1902) · N. Zhukovsky (1906)

양력은 압력의 적분이다. 그런데 그 적분을 끝까지 하면, 남는 건 순환 하나다.

쿠타-주코프스키 정리(Kutta–Joukowski theorem)는 비점성·비압축·2차원 정상 유동에 놓인 물체가 받는 단위 스팬당 양력이 유입 속도와 물체를 감싼 순환의 곱, 즉 $L' = \rho_\infty V_\infty \Gamma$ 로 정확히 주어지고 그 힘은 유입 방향에 수직이라는 정리다. 형상이 익형이든 원기둥이든 고양이 단면이든 상관없다. 물체의 기하학은 오직 순환 $\Gamma$ 하나에 압축되어 힘에 들어간다.

이 문장이 무서운 이유는 정보량 때문이다. 표면 압력 분포는 $\theta$ 의 함수이고 무한 차원인데, 그중 힘에 기여하는 것은 원방 전개의 딱 한 항이다. 나머지는 전부 상쇄된다. 그래서 20세기 초 항공역학은 “익형을 잘 설계한다”는 문제를 **”$\Gamma$ 를 잘 만든다”**는 문제로 갈아탈 수 있었고, 얇은 익형 이론·양력선 이론·와류격자법·패널법이 전부 그 위에 세워졌다.

다만 정리는 $\Gamma$ 가 얼마인지는 말해 주지 않는다. 그건 쿠타 조건의 일이다. 이 문서는 ”$\Gamma$ 가 주어졌을 때 힘이 무엇인가”만 다루고, ”$\Gamma$ 가 왜 그 값인가”는 그쪽으로 넘긴다.

2. 순환이라는 양[편집]

순환은 닫힌 곡선 $C$ 를 따라 잡은 속도의 선적분이다.

$$\Gamma = \oint_C \mathbf{V}\cdot d\mathbf{s}$$

물체 바깥이 비회전($\nabla\times\mathbf{V}=0$)이면, 스토크스 정리에 의해 물체를 감싸는 어떤 곡선을 잡아도 $\Gamma$ 는 같다. 곡선을 부풀리거나 찌그러뜨려도 값이 안 변한다는 이 성질 덕분에 $\Gamma$ 는 “물체가 들고 있는 고유한 숫자”의 자격을 얻는다. 정리가 형상에 무관해지는 근본 이유가 여기 있다.

부호 관습은 문헌마다 갈려서 사고의 단골 원인이다. 위 정의(반시계 양)를 쓰면 오른쪽으로 흐르는 유동에서 양의 $\Gamma$ 는 아래로 누르는 힘을 준다. 항공 쪽 교재가 $\Gamma$ 를 시계 방향 양으로 뒤집어 정의해 놓고 $L'=\rho V_\infty\Gamma$ 를 그대로 쓰는 이유가 그것이다. 어느 쪽이든 크기는 $\rho_\infty V_\infty|\Gamma|$, 방향은 유입 속도를 $\Gamma$ 의 회전 방향과 반대로 90° 돌린 쪽이라고 외우면 부호 사고를 안 낸다.

3. 유도 1 — 검사체적을 무한히 키운다[편집]

물체를 감싸는 반지름 $R$ 의 큰 원을 검사체적 경계로 잡는다. 정상 유동의 운동량 정리는

$$\mathbf{F} = -\oint_C p\,\mathbf{n}\,ds - \oint_C \rho\,\mathbf{V}(\mathbf{V}\cdot\mathbf{n})\,ds$$

이다. 원방에서 속도를 균일류와 교란으로 쪼개면 $\mathbf{V} = \mathbf{V}_\infty + \mathbf{v}'$ 이고, 물체가 만드는 교란은 크기 순서로 점와류 항 $O(1/R)$ + 더블릿·소스 항 $O(1/R^2)$ 이다. 여기서 세 가지가 동시에 일어난다.

• 베르누이 방정식에서 압력 교란은 $p' \simeq -\rho V_\infty u'$ 로 $O(1/R)$.
• 적분 경로의 길이는 $2\pi R$ 로 $O(R)$.
• 따라서 두 $O(1/R)$ 항의 곱, 즉 $O(1/R^2)$ 짜리는 사라지고, $O(1/R)\times O(R) = O(1)$ 인 항만 살아남는다.

살아남는 항을 모으면 각도 적분에서 $\oint \sin^2$ 꼴만 남아 정확히 $\rho_\infty V_\infty\Gamma$ 가 되고, 유입 방향 성분은 $\oint \sin\theta\cos\theta\,d\theta = 0$ 으로 완벽히 상쇄된다. $R\to\infty$ 극한에서 남는 결론이 곧 정리다.

$$\boxed{\;L' = \rho_\infty V_\infty \Gamma,\qquad D' = 0\;}$$

여기서 두 번째 등식이 바로 달랑베르 역설이다. 양력과 역설은 별개의 두 사실이 아니라 같은 극한 계산의 두 성분이다. 이 점이 중요하다 — 비점성 이론에게 양력을 받아 오면서 항력을 사양할 방법은 없다. 둘은 한 세트로 온다.¹

4. 유도 2 — 복소해석 세 줄[편집]

2차원 비회전 유동은 복소 퍼텐셜 $w(z)$ 로 쓰면 계산이 우스울 만큼 짧아진다. 복소속도 $dw/dz = u - iv$ 의 원방 로랑 전개는

$$\frac{dw}{dz} = V_\infty e^{-i\alpha} - \frac{i\Gamma}{2\pi z} + \frac{c_2}{z^2}+\cdots$$

이고($\alpha$ 는 받음각), 물체에 걸리는 힘은 블라시우스 정리

$$X - iY = \frac{i\rho}{2}\oint_C \left(\frac{dw}{dz}\right)^2 dz$$

로 준다. 제곱을 전개하면 $1/z$ 의 계수는 교차항 $2V_\infty e^{-i\alpha}\cdot(-i\Gamma/2\pi)$ 뿐이고, 유수정리로

$$X - iY = i\rho V_\infty \Gamma e^{-i\alpha}$$

가 나온다. 크기는 $\rho V_\infty\Gamma$, 위상은 유입 방향에서 정확히 90°. 한 줄의 유수(residue)가 양력이고, 그 유수를 만드는 것이 순환이다. 소스 항($z^{-1}$ 의 실수부)이 없다는 사실 — 물체가 유체를 삼키지도 뱉지도 않는다는 사실 — 이 항력 0을 보증한다는 것도 이 전개에서 곧바로 읽힌다.²

5. 원기둥에서 익형으로 — 고전 절차[편집]

정리를 실제 익형에 쓰려면 $\Gamma$ 를 손에 넣어야 하고, 컴퓨터 없던 시절의 답이 등각사상이었다. 절차는 네 단계다.

1. 원기둥 주위 해를 쓴다. 균일류 + 더블릿 + 점와류를 겹치면 반지름 $R$ 원기둥 주위 유동의 완전한 해석해가 나온다. 이때 $\Gamma$ 는 아직 자유 매개변수다.
2. 정체점을 지정한다. 원 위에서 정체점의 위치는 $\Gamma$ 에 따라 미끄러진다. 사상 후 익형의 뒷전이 될 원 위의 그 점에 정체점을 얹으라고 요구하는 것이 쿠타 조건이다. 이 요구가 $\Gamma$ 를 하나로 못 박는다.
3. 주코프스키 변환으로 사상한다. $w = z + a^2/z$ 는 원점 근처 원을 익형으로 편다. 중심을 실축 왼쪽으로 밀면 두께가, 허축 위로 밀면 캠버가 생긴다. 등각사상은 라플라스 방정식을 보존하므로, 사상된 유동장은 여전히 익형 주위의 진짜 해다.
4. 정리를 적용한다. $\Gamma$ 는 사상으로 변하지 않는 양이므로($\oint$ 는 곡선의 상에서 그대로다) 그대로 $L'=\rho V_\infty\Gamma$ 에 넣으면 끝.

2번에서 나오는 순환은 깔끔한 닫힌형이다. 원의 반지름이 $R$, 중심 오프셋이 만드는 캠버각이 $\beta$ 일 때

$$\Gamma = 4\pi R V_\infty \sin(\alpha+\beta)$$

이고, 두께 0·캠버 0 인 평판($R=a$, 시위 $c=4a$)에 넣으면

$$C_l = \frac{L'}{\tfrac12\rho V_\infty^2 c} = 2\pi\sin\alpha \;\approx\; 2\pi\alpha$$

라는, 항공역학에서 가장 유명한 숫자가 떨어진다. 얇은 익형 이론이 적분방정식을 풀어서 얻는 $2\pi$ 를 여기서는 삼각함수 한 줄로 얻는 셈이다. 두께가 있는 주코프스키 익형은 양력기울기가 $2\pi$ 보다 조금 커지며 흔히 $2\pi(1+0.77\,t/c)$ 로 근사한다.

물론 사상이 만드는 익형은 뒷전이 끼인각 0인 첨점이라 실제로 깎을 수 없고, 유한 뒷전각이 필요하면 카르만-트레프츠 변환으로 일반화해야 한다. 사상 자체의 성질은 등각사상이 다루므로 여기서 반복하지 않는다.

6. 전제를 하나씩 깨 보면[편집]

정리가 성립하는 조건은 좁다. 조건 하나가 무너질 때 무슨 일이 벌어지는지 아는 것이 실무 감각이다.

깨는 전제	벌어지는 일	처방
비점성	마찰항력·형상항력이 생기고 실효 순환이 줄어든다	경계층 결합, XFOIL 계열
2차원	순환이 스팬 방향으로 변하고 후류 와류가 생긴다	양력선 이론, 와류격자법
정상	뒷전에서 와도가 계속 방출돼 후류가 되먹임한다	비정상 익형 이론, 바그너·테오도르센
비압축	밀도 변화가 $\Gamma$ 와 압력 적분을 함께 바꾼다	프란틀-글라워트 보정, 초음속은 아예 다른 이론
부착유동	순환이 뒷전과 무관하게 결정된다	난류 모델링 얹은 비정상 CFD

압축성에 대해 한 마디 더. 아음속 선형 이론에서는 정리의 형태가 그대로 살아남고 마하수 효과가 $\Gamma$ 안으로 흡수된다. 하지만 국부적으로 초음속 영역과 충격파가 생기면 원방 전개에 엔트로피를 남기는 파가 끼어들어 조파항력이 나오고, 힘이 유입 방향에 수직이라는 결론이 무너진다. 초음속에서 양력과 항력을 함께 논하려면 압축성 유동의 언어로 갈아타야 한다.

7. 유한 익폭 — 유도항력은 왜 정리 바깥이 아닌가[편집]

흔한 오해가 있다. “쿠타-주코프스키는 항력을 0이라 하는데 실제 날개엔 유도항력이 있으니 정리가 틀렸다.” 틀렸다는 결론은 틀렸다. 정리는 2차원 진술이고, 유한 익폭에서는 그 2차원 진술을 스팬 방향 단면마다 국소적으로 적용한다.

날개 끝에서 압력이 새면서 후류로 와류가 흘러나가고, 그 후류가 날개 위치에 내리흐름(downwash) $w$ 를 만든다. 그러면 단면이 실제로 보는 유입 속도는 $V_\infty$ 가 아니라 그것에 $w$ 를 더한 벡터이고, 방향이 유도각 $\alpha_i \approx w/V_\infty$ 만큼 아래로 기울어져 있다. 정리는 힘이 그 국소 유입 속도에 수직이라고 말하므로, 힘 벡터도 통째로 뒤로 기울어진다. 그 기울어진 힘의 자유류 방향 성분이

$$D_i' = \rho\, w\, \Gamma$$

즉 유도항력이다. 정리를 위반해서 생긴 게 아니라, 정리를 3차원에 성실히 적용하면 자동으로 나오는 항이다. 물체를 감싼 순환의 정의가 스팬 좌표에 의존하게 되는 순간부터 벌어지는 일이고, 그 스팬 분포 $\Gamma(y)$ 를 실제로 푸는 것이 양력선 이론과 와류격자법의 몫이다. 후류를 무한히 하류에서 잘라 보는 트레프츠 평면 계산도 결국 같은 이야기를 에너지 쪽에서 다시 하는 것이다.

3차원 벡터 형태로는 와류 필라멘트 한 토막에 걸리는 단위 길이당 힘이

$$\mathbf{f}' = \rho\,\mathbf{V}\times\boldsymbol{\Gamma}$$

로 쓰인다. 와류격자법 코드에서 각 결합와 요소의 힘을 뽑을 때 쓰는 바로 그 식이며, 2차원은 $\boldsymbol{\Gamma}$ 가 지면에 수직인 특수한 경우다.

8. 수치해석에서의 위치[편집]

CFD 시대에도 이 정리는 죽지 않았다. 쓰임새가 셋이다.

• 저차 코드의 힘 추출기. 와류격자법·패널법·양력선 이론은 압력을 적분하지 않고 순환에서 곧장 힘을 뽑는다. 압력 적분보다 격자 의존성이 훨씬 작다는 실용적 장점이 있다.
• 검증 케이스. 주코프스키 익형은 나비에-스토크스 방정식 솔버가 갖기 어려운 해석해를 준다. 새 솔버의 격자 수렴 시험을 여기서 시작하는 것이 국룰이다.
• 후처리 진단. 익형 주위에 임의의 폐곡선을 잡아 $\Gamma$ 를 적분한 값이 압력 적분 양력과 어긋나면, 그 차이가 곧 점성 손실 + 수치 오차의 크기다. 서로 다른 크기의 곡선에서 $\Gamma$ 가 안 맞으면 격자가 와도를 흘리고 있다는 뜻이다.

마지막 항목은 의외로 자주 사람을 살린다. RANS 결과가 이상한데 이유를 모르겠으면, 익형에서 시위의 2배·5배·10배 거리에 원을 그려 순환을 재 보면 된다. 세 값이 다르면 격자를 의심하고, 세 값이 같은데 양력이 안 맞으면 난류 모델링을 의심한다.³

9. 여담[편집]

이름이 두 사람인 데는 사연이 있다. 쿠타는 1902년 논문에서 곡면 평판의 양력을 순환으로 계산했지만 정리를 일반 형상에 대한 정리로 진술하지는 않았고, 주코프스키가 1906년에 일반적 형태로 발표했다. 랑체스터는 그보다 앞선 1894년경에 순환-양력 그림을 갖고 있었지만 논문 심사에서 거절당해 1907년에야 책으로 냈다.⁴ 그래서 영어권에서 이 정리를 “Kutta–Joukowski”라 부르고, 러시아에서는 “주콥스키 정리”, 영국 일부 문헌은 “Lanchester–Prandtl 계보”를 강조한다. 이 바닥 국룰: 먼저 생각한 사람과 이름이 붙는 사람은 다르다.

같은 정리가 회전하는 공에도 그대로 적용된다(마그누스 효과). 표면 마찰이 공기를 끌고 돌면서 순환을 만들고, $\rho V_\infty\Gamma$ 만큼의 옆힘이 나온다. 축구의 바나나킥과 야구의 커브볼이 이 문서와 같은 식으로 계산된다는 뜻이다. 다만 실제 스포츠공은 표면 거칠기와 경계층 천이가 지배적이라 $\Gamma$ 를 예측하는 쪽이 훨씬 어렵다 — 여기서도 정리는 쉽고 $\Gamma$ 가 어렵다.

10. 관련 문서[편집]

• 쿠타 조건 · 얇은 익형 이론 · 양력선 이론
• 와류격자법 · 패널법 · 공력해석
• 등각사상 · 포텐셜 유동 · 라플라스 방정식
• 와도 · 보텍스 방법 · 베르누이 방정식
• 경계층 · 유동박리 · 압축성 유동
• 달랑베르 역설 · 유도항력 · 마그누스 효과

11. Footnotes[편집]

1. 그래서 “비점성 코드로 항력을 계산했다”는 말은 원리적으로 성립하지 않는다. 실제로 나오는 값은 전부 수치 소산이 만든 가짜 항력이다. 격자를 촘촘히 할수록 그 값이 0으로 내려가는 것을 보고 “수렴한다”고 좋아하는 사람이 있는데, 사실은 코드가 정직해지고 있는 것이다. ↩

2. 블라시우스 정리는 1910년 발표 당시 “복소적분 한 번으로 유체력을 다 준다”는 이유로 꽤나 화제였다. 지금 봐도 계산량 대비 산출이 반칙 수준인데, 대가는 2차원·비점성·정상이라는 좁은 무대다. 반칙적으로 강력한 도구가 반칙적으로 좁은 조건을 다는 것은 수리물리의 국룰이다. ↩

3. 순환을 재는 폐곡선이 후류를 가로지르면 값이 당연히 달라진다. 후류에는 와도가 있어서 그 안쪽은 비회전이 아니기 때문이다. 이걸 모르고 “순환이 안 맞는다”며 격자를 세 번 다시 만든 사례가 실제로 있다. 곡선을 후류 위로 넘겨서 잡아야 한다. ↩

4. 랑체스터의 원고를 반려한 심사평의 요지는 “물리적 근거가 부족하다”였다고 전해진다. 13년 뒤 그 아이디어로 사람이 하늘을 날았다. 학술지 리젝은 예나 지금이나 인류사에 흔적을 남긴다. ↩