얇은 익형 이론

1. 개요[편집]

얇은 익형 이론 Thin airfoil theory
정립	Max Munk (1922) · Hermann Glauert (1926)
모형	캠버선 위의 와류시트 γ(x)
지배식	코시 특이핵을 가진 1종 적분방정식
해법	글라우어트 변환 + 푸리에 사인 급수
대표 결과	cl = 2π(α − αL=0) · 공력중심 c/4
유효 범위	얇음 · 작은 캠버 · 작은 받음각 · 부착유동

두께를 0으로, 받음각을 작게, 캠버를 얕게. 세 번 양보하고 나면 익형 문제가 손으로 풀린다.

얇은 익형 이론(thin airfoil theory)은 익형을 두께 없는 캠버선으로 치환하고 그 위에 와류시트 $\gamma(x)$ 를 분포시킨 뒤, 유동이 캠버선에 접해야 한다는 조건을 특이적분방정식으로 세워 해석적으로 푸는 2차원 공력 이론이다. 1920년대에 뭉크와 글라우어트가 완성했고, 100년이 지난 지금도 양력계수가 받음각에 선형이고 기울기가 $2\pi$ 이며 공력중심이 1/4 시위점에 있다는 세 문장을 세상에 공급하는 유일한 출처다.

이 이론의 지위는 묘하다. 정확도로만 보면 패널법이나 CFD에 상대가 안 된다. 그런데 양력선 이론이 각 단면의 양력 기울기 $a_0$ 를 요구할 때, 와류격자법이 격자 배치의 정당성을 요구할 때, 실험자가 “이 익형은 왜 $\alpha_{L=0}$ 이 −2°인가”를 물을 때 — 답을 주는 것은 늘 이쪽이다. 수치해가 참고할 정답이 있는 유일한 익형 이론이라는 게 진짜 값어치다.¹

2. 세 번의 양보[편집]

가정은 명시적이고, 세 가지다.

• 두께 $\to$ 0. 두께 분포는 대칭이라 양력에 1차 기여를 하지 않는다는 것이 근거다. 두께는 소스 분포로, 캠버·받음각은 와류 분포로 갈라지며 선형화된 문제에서 둘이 분리된다.
• 캠버 $z(x)$ 와 그 기울기가 작다. $|dz/dx| \ll 1$.
• 받음각이 작다. $\sin\alpha \approx \alpha$, $\cos\alpha \approx 1$.

이 셋을 받아들이면 결정적인 단순화가 따라온다. 경계조건을 캠버선이 아니라 시위선(chord line) 위에서 부과해도 1차 정확도로는 같다. 미지 함수가 놓이는 자리가 곡선에서 직선 $0 \le x \le c$ 로 내려앉는 것이며, 이 한 번의 평탄화가 문제를 손으로 풀 수 있게 만든다. 와류격자법이 3차원에서 캠버면 대신 평면에 격자를 까는 것도 같은 논리의 연장이다.

3. 와류시트와 접선 조건[편집]

캠버선을 세기 $\gamma(\xi)$ 인 와류시트로 대체한다. 시트는 그 자체로 라플라스 방정식을 만족하므로 남은 일은 경계조건뿐이다. 세기 $\gamma(\xi)d\xi$ 인 미소 와류가 위치 $x$ 에 만드는 수직 유도속도는

$$dw(x) = -\frac{\gamma(\xi)\,d\xi}{2\pi (x-\xi)}$$

이고, 시트 전체를 적분한 $w(x)$ 가 자유류의 캠버선 수직 성분 $V_\infty(\alpha - dz/dx)$ 를 정확히 상쇄해야 유동이 캠버선에 접한다. 그 결과가 얇은 익형 이론의 기본 방정식이다.

$$\frac{1}{2\pi}\int_0^c \frac{\gamma(\xi)}{x-\xi}\,d\xi = V_\infty\left(\alpha - \frac{dz}{dx}\right)$$

핵이 $1/(x-\xi)$ 라 적분은 $\xi = x$ 에서 특이하며, 코시 주치로 읽어야 한다. 형태로는 코시 핵을 가진 제1종 특이적분방정식(에어포일 방정식이라 불리기도 한다)이고, 이런 방정식은 해가 유일하지 않다. 부족한 한 조건을 채우는 것이 쿠타 조건이다.

$$\gamma(c) = 0$$

즉 뒷전에서 와류시트 세기가 사라져야 한다. 이 한 줄을 넣기 전까지는 이론이 아직 양력을 모른다.

4. 글라우어트 변환[편집]

특이적분을 정면으로 다루는 대신, 글라우어트는 좌표를 바꾼다.

$$\xi = \frac{c}{2}(1-\cos\theta), \qquad x = \frac{c}{2}(1-\cos\theta_0), \qquad \theta,\theta_0 \in [0,\pi]$$

앞전이 $\theta=0$, 뒷전이 $\theta=\pi$ 에 대응한다. 그리고 미지 함수를 다음 급수로 전개한다.

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

이 형태는 우연이 아니라 답을 알고 고른 것이다. 두 가지가 미리 심어져 있다.

• $\theta \to \pi$(뒷전)에서 $(1+\cos\theta)/\sin\theta = \cot(\theta/2) \to 0$ 이고 $\sin n\pi = 0$ 이므로 $\gamma(c)=0$ — 쿠타 조건이 자동 만족된다.
• $\theta \to 0$(앞전)에서 $\cot(\theta/2)\to\infty$, 즉 $\gamma \sim 1/\sqrt{x}$ 로 발산한다. 이건 결함이 아니라 물리적으로 의미 있는 특이성이다 — 두께 0인 평판의 앞전은 곡률 반경이 0이므로 유동이 무한 속도로 돌아가야 하며, 그 무한 속도가 만드는 유한한 흡입력이 앞전 흡입력(leading-edge suction)이다. 이 항이 없으면 평판의 양력이 자유류에 수직이 아니게 나온다.

이제 저 유명한 글라우어트 적분

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

를 쓰면 특이적분이 깨끗하게 닫히고, 기본 방정식이 항별로 분해되어 계수가 캠버선 기울기의 푸리에 계수로 곧장 나온다.

$$A_0 = \alpha - \frac{1}{\pi}\int_0^{\pi} \frac{dz}{dx}\,d\theta_0, \qquad A_n = \frac{2}{\pi}\int_0^{\pi} \frac{dz}{dx}\cos n\theta_0 \,d\theta_0 \quad (n \ge 1)$$

받음각은 $A_0$ 에만 들어가고, $A_1, A_2, \dots$ 는 캠버선 형상만으로 결정된다는 점이 이후 모든 결과의 열쇠다.

5. 결과 — 양력과 모멘트[편집]

순환은 시트를 적분하면 나온다. $\Gamma = \int_0^c \gamma\,d\xi = cV_\infty(\pi A_0 + \tfrac{\pi}{2}A_1)$ 이고, 쿠타-주코프스키 정리 $L' = \rho V_\infty \Gamma$ 를 거쳐

$$c_l = \pi\left(2A_0 + A_1\right)$$

를 얻는다. $A_0$ 의 정의를 대입해 정리하면 교과서의 그 식이 된다.

$$c_l = 2\pi\left(\alpha - \alpha_{L=0}\right), \qquad \alpha_{L=0} = -\frac{1}{\pi}\int_0^{\pi}\frac{dz}{dx}\left(\cos\theta_0 - 1\right)d\theta_0$$

여기서 두 가지가 동시에 나온다.

• 양력 기울기는 캠버와 무관하게 항상 $2\pi$ /rad 이다. 캠버가 하는 일은 곡선을 위아래로 통째로 옮기는 것뿐이다.
• 영양력 받음각 $\alpha_{L=0}$ 은 캠버선만으로 결정된다. 대칭 익형이면 $dz/dx \equiv 0$ 이라 $\alpha_{L=0}=0$ 이고, $c_l = 2\pi\alpha$ 로 끝난다.

모멘트도 같은 계수들로 떨어진다. 앞전 기준 피칭 모멘트 계수는

$$c_{m,\text{le}} = -\frac{\pi}{2}\left(A_0 + A_1 - \frac{A_2}{2}\right)$$

이고, 이를 1/4 시위점으로 옮기면($c_{m,c/4} = c_{m,\text{le}} + c_l/4$) $A_0$ 항이 정확히 소거된다.

$$\boxed{\;c_{m,c/4} = \frac{\pi}{4}\left(A_2 - A_1\right)\;}$$

우변에 $\alpha$ 가 없다. 받음각을 어떻게 바꿔도 1/4 시위점 둘레의 모멘트가 변하지 않는다 — 이것이 공력중심(aerodynamic center)의 정의이므로, 얇은 익형의 공력중심은 $x_{ac} = c/4$ 다. 대칭 익형이면 $A_1 = A_2 = 0$ 이라 $c_{m,c/4}=0$ 이고, 압력중심마저 1/4 시위점에 붙는다. 캠버가 있으면 압력중심은

$$\frac{x_{cp}}{c} = \frac{1}{4}\left[1 + \frac{\pi}{c_l}(A_1 - A_2)\right]$$

로 $c_l$ 에 따라 움직이지만, 공력중심은 꿈쩍도 하지 않는다. 비행기의 세로 안정성을 논할 때 압력중심이 아니라 공력중심을 쓰는 이유가 정확히 이 성질이다.²

6. 예시로 감 잡기[편집]

가장 흔한 두 캠버선을 넣어 보면 이론이 무슨 일을 하는지 보인다.

캠버선	$\alpha_{L=0}$	$c_{m,c/4}$	비고
평판·대칭 익형	$0^\circ$	$0$	$c_l = 2\pi\alpha$
NACA 4자리 2% 캠버(2412 계열)	약 $-2^\circ$	약 $-0.05$	실측과 잘 맞음
뒤로 젖힌 리플렉스 캠버	음(작음)	$\approx 0$ 또는 양	무미익기가 쓰는 이유

세 번째 줄이 실용적으로 재미있다. $c_{m,c/4} = (\pi/4)(A_2-A_1)$ 이므로 $A_2$ 를 키우면 모멘트를 0이나 양으로 만들 수 있고, 이는 캠버선 뒷부분을 위로 젖히는(reflex) 형상에 해당한다. 꼬리날개 없이 세로 트림을 맞춰야 하는 전익기와 행글라이더가 이 계산의 직접 수혜자다.

7. 실험과 얼마나 맞나[편집]

$2\pi$ /rad 를 도 단위로 바꾸면 $2\pi/57.30 = 0.1097\ /^\circ$ 다. 그리고 실제 익형의 측정 양력 기울기는 레이놀즈수와 표면 상태에 따라 대략 $0.10 \sim 0.11\ /^\circ$ 에 들어온다. 두께도 점성도 무시한 이론치고는 어이없을 만큼 잘 맞는데, 이건 두 오차가 서로 반대 방향이라 상당 부분 상쇄되기 때문이다.

• 두께는 기울기를 올린다. 두께가 있으면 표면 속도가 캠버선 위보다 빨라져 순환이 커진다. 글라우어트의 근사식 $dc_l/d\alpha \approx 2\pi(1 + 0.77\,t/c)$ 에 따르면 12% 두께 익형은 이론상 $2\pi$ 보다 9%쯤 크다. 실제로 패널법으로 NACA 0012 를 풀면 $2\pi\alpha$ 보다 10% 안팎 큰 값이 나온다.
• 점성은 기울기를 내린다. 뒷전 근처 경계층의 배제두께가 실효 캠버를 깎고(decambering) 실효 뒷전을 두껍게 만들어 순환을 줄인다. 레이놀즈수가 낮을수록 이 효과가 크다.

$\alpha_{L=0}$ 예측은 기울기보다도 잘 맞는다. 캠버선의 1차 정보만 쓰는 양이라 두께 효과가 거의 안 들어오기 때문이다. 공력중심 역시 실측치가 대개 $0.24c \sim 0.27c$ 로 이론값 $0.25c$ 를 감싸고 있다.

8. 어디서 깨지는가[편집]

• 실속. 이론에는 점성이 없으니 $\alpha$ 를 아무리 키워도 $c_l$ 이 선형으로 자란다. 실제 익형은 $\alpha \approx 12\sim18^\circ$ 에서 유동박리로 무너진다. 얇은 익형 이론은 실속각을 예측하지 못할 뿐 아니라 실속의 존재조차 모른다.
• 앞전 특이성. $\gamma \sim 1/\sqrt{x}$ 는 두께 0의 산물이다. 실제 익형의 앞전 반경이 유한하므로 그 근처의 압력 분포는 완전히 다르며, 얇은 익형 이론으로는 앞전 흡입 피크를 계산할 수 없다. 얇고 앞전이 날카로운 익형이 저받음각에서도 앞전 박리 버블을 만드는 현상 역시 이 이론 바깥이다.
• 두꺼운 익형·큰 캠버. 선형화가 근거를 잃는다. 두께비 12% 안팎이 흔히 언급되는 실용 한계이며, 그 위로는 두께 보정 없이 쓰면 안 된다.
• 압축성. 비압축 가정이므로 아음속에서는 프란틀-글라우어트 보정으로 $c_l$ 을 $1/\sqrt{1-M^2}$ 배 해 얹지만, 천음속에서 충격파가 서면 그 보정마저 무너진다.
• 3차원. 이 이론은 무한 스팬 2차원 익형의 것이다. 유한 날개에서는 후류가 유효 받음각을 깎아 기울기가 $2\pi$ 보다 작아지며, 그 이야기는 양력선 이론과 와류격자법의 몫이다.

9. 남긴 것[편집]

얇은 익형 이론이 후대에 남긴 진짜 유산은 숫자 $2\pi$ 가 아니라 구조다. “미지의 와도 분포 + 유동 접선 조건 = 특이적분방정식”이라는 틀은 그대로 3차원으로 확장되어 양력선 이론의 적분미분방정식이 되고, 그것을 이산화하면 와류격자법이 되며, 물체 표면으로 옮겨 두께까지 실으면 패널법이 된다. 100년치 저차 공력해석이 전부 이 한 장의 확장판인 셈이다.

역설적으로 이론이 가장 자주 쓰이는 곳은 익형 설계가 아니라 검증이다. 새로 짠 와류격자법 코드를 종횡비 100짜리 평판 날개에 돌려 $c_l/\alpha$ 가 $2\pi$ 근처로 오는지 보는 것이 이 바닥의 첫 번째 단위 테스트고, 여기서 어긋나면 그 코드는 어디선가 틀린 것이다. 손으로 푼 정답이 하나 있다는 것은 수치해석에서 생각보다 귀한 자산이다.³

10. 관련 문서[편집]

• 쿠타 조건 · 쿠타-주코프스키 정리
• 와류격자법 · 패널법
• 양력선 이론 · 공력해석
• 포텐셜 유동 · 라플라스 방정식 · 와도
• 경계층 · 유동박리 · 레이놀즈수
• 프란틀-글라우어트 보정 · 압축성 유동
• NACA 익형 · 등각사상 · 복소해석

11. Footnotes[편집]

1. 이론이 나오기 전 익형 해석의 정석은 등각사상이었다. 원기둥의 해를 알고 있으니 사상하면 된다는 발상인데, 사상 함수가 존재하는 형상만 풀 수 있다는 치명적 제약이 있었다. 얇은 익형 이론은 “임의의 캠버선”을 받아들인다는 점에서 처음으로 설계자가 원하는 형상을 계산할 수 있게 해 준 도구였다. ↩

2. 학부 비행역학 시험에서 압력중심과 공력중심을 뒤집는 사고가 매년 반복된다. 외우기 좋은 구분은 이렇다 — 압력중심은 실제로 힘이 걸리는 자리라 받음각 따라 돌아다니고(양력이 0에 가까워지면 무한대로 도망간다), 공력중심은 모멘트가 안 변하는 자리라 얌전히 c/4에 앉아 있다. 안정성 계산에서 도망가는 점을 기준으로 삼으면 방정식이 지옥이 된다. ↩

3. 반대로 이 테스트만 통과하고 안심하는 것도 위험하다. $2\pi$ 는 2차원 극한을 맞히는 것뿐이라, 스팬 방향 이산화가 엉망이어도 종횡비를 충분히 키우면 통과한다. 유도항력이 $C_L^2/(\pi\,e\,AR)$ 로 오는지, 타원 날개에서 $e$ 가 1에 붙는지까지 봐야 코드가 3차원을 제대로 푼다고 말할 수 있다. ↩