프란틀-글라우어트 보정

편집 역사 토론
유체역학 전산유체역학 마지막 수정: 2026-08-23 04:49:05

1. 개요[편집]

프란틀-글라우어트 보정
Prandtl–Glauert Correction
공식$C_p = C_{p0}/\sqrt{1-M_\infty^2}$
적용 범위아음속, 임계 마하수 미만, 얇고 작은 받음각
근거선형화 압축성 퍼텐셜 방정식의 좌표 신축
3차원 확장괴테르트 규칙 (Göthert's rule)
고차 대안카르만-첸 (1941) · 라이톤 (1951)
초음속 대응애커렛 이론, $\sqrt{M_\infty^2-1}$

프란틀-글라우어트 보정(Prandtl–Glauert correction)은 비압축성으로 계산한 압력계수를 1−M∞2\sqrt{1-M_\infty^2} 로 나눠 아음속 압축성 효과를 반영하는 1차 보정이다.

Cp=Cp01−M∞2C_p=\frac{C_{p0}}{\sqrt{1-M_\infty^2}}

Cp0C_{p0} 는 같은 형상·같은 받음각에 대한 비압축성 해이고, M∞M_\infty 는 자유류 마하수다. 분모를 β=1−M∞2\beta=\sqrt{1-M_\infty^2} 로 줄여 쓰는 것이 관례다. 양력계수도 같은 인자로 커지므로 2차원 얇은 익형의 양력곡선 기울기는 2π/β2\pi/\beta 가 된다.

한 줄짜리 나눗셈이 100년 가까이 살아남은 이유는 비압축성 해석 한 번으로 마하수 전 대역을 커버할 수 있다는 경제성 때문이다. 풍동에서 저속으로 잰 데이터를 순항 마하수로 환산하는 것도, 패널법이나 더블릿 격자법 같은 선형 코드에 압축성을 넣는 것도 결국 이 인자 하나다. 대가는 명확하다 — 임계 마하수를 넘으면 그냥 틀린다. 그것도 안전한 쪽이 아니라 위험한 쪽으로 틀린다.

2. 선형화 퍼텐셜 방정식에서의 유도[편집]

비점성·비회전·등엔트로피 유동을 속도 퍼텐셜로 쓰면 전퍼텐셜 방정식이 나오는데, 이것은 비선형이다. 여기에 얇은 물체·작은 받음각 가정을 넣어 자유류에 작은 섭동만 얹힌다고 하자. 즉 Φ=U∞x+φ\Phi=U_\infty x+\varphi 이고 φ\varphi 의 미분이 U∞U_\infty 에 비해 작다고 놓으면, 2차 이상의 항이 떨어져 나가고 선형화 소섭동 퍼텐셜 방정식이 남는다.

(1−M∞2) φxx+φyy+φzz=0(1-M_\infty^2)\,\varphi_{xx}+\varphi_{yy}+\varphi_{zz}=0

M∞<1M_\infty<1 이면 이 방정식은 타원형이다. 그리고 딱 한 항의 계수만 1이 아니므로, 좌표를 늘려서 그 계수를 없앨 수 있다. 이것이 프란틀-글라우어트 변환이다.

ξ=x,η=βy,ζ=βz,β=1−M∞2\xi=x,\qquad \eta=\beta y,\qquad \zeta=\beta z,\qquad \beta=\sqrt{1-M_\infty^2}

이 좌표에서 방정식은 라플라스 방정식이 되고, 압축성 문제는 우리가 이미 풀 줄 아는 비압축성 포텐셜 유동 문제로 바뀐다. 압축성 유동의 어려움을 계산이 아니라 좌표계로 해결한 셈이다.1

압력계수도 같은 가정 아래 선형화되어

Cp≃−2u′U∞=−2U∞∂φ∂xC_p\simeq-\frac{2u'}{U_\infty}=-\frac{2}{U_\infty}\frac{\partial\varphi}{\partial x}

가 되며, 이 관계와 변환을 함께 쓰면 최종 보정식이 나온다.

3. 괴테르트 규칙 — 3차원에서 조심할 것[편집]

여기서 흔히 생략되는 함정이 있다. 좌표를 신축하면 물체 형상도 같이 일그러진다. y,zy,z 방향 치수가 β\beta 배로 줄어들기 때문에, 변환된 공간에서 우리가 푸는 문제는 원래 형상이 아니라 두께비·캠버·받음각·종횡비가 모두 β\beta 배로 줄어든 다른 형상에 대한 비압축성 문제다.

이 대응 관계를 정확히 정리한 것이 괴테르트 규칙(Göthert’s rule)이고, 압력계수 관계는 다음과 같다.

Cp=1β2 CˉpC_p=\frac{1}{\beta^2}\,\bar C_p

Cˉp\bar C_p 는 일그러진 형상에 대한 비압축성 해다. 그렇다면 왜 평소에는 1/β1/\beta 를 쓰는가?

2차원에서는 선형 이론이 형상에 대해 선형이기 때문이다. 얇은 익형 이론에서 Cp0C_{p0} 는 두께비·캠버·받음각에 비례하므로, 일그러진 형상의 해는 원래 형상 해의 β\beta 배다. 대입하면 Cp=(1/β2)(β Cp0)=Cp0/βC_p=(1/\beta^2)(\beta\,C_{p0})=C_{p0}/\beta — 익숙한 그 공식이 나온다.

3차원에서는 이 상쇄가 성립하지 않는다. 종횡비까지 β\beta 배로 줄어드는데, 유한 날개의 양력은 종횡비에 대해 선형이 아니기 때문이다. 그 결과 실제 3차원 날개의 압축성 효과는 1/β1/\beta 보다 약하며, 종횡비가 낮을수록 더 약해진다. 후퇴익 전투기의 낮은 종횡비 날개에 2차원 공식을 그대로 들이대면 압축성 효과를 과대평가한다. 델타익이 천음속에서 상대적으로 얌전한 이유의 일부가 여기 있다.

4. 왜 발산하는가 — 임계 마하수라는 진짜 한계[편집]

M∞→1M_\infty\to1 이면 β→0\beta\to0 이고 Cp→∞C_p\to\infty 다. 이 발산을 프란틀-글라우어트 특이점이라 부른다. 물론 실제 비행기 날개 위의 압력이 무한대가 되지는 않는다. 그 훨씬 전에 이론이 스스로 전제를 배반한다.

이론이 무너지는 지점은 M∞=1M_\infty=1 이 아니라 임계 마하수 McrM_{cr} 이다. 물체 표면 어딘가에서 국소 유속이 처음으로 음속에 도달하는 자유류 마하수이며, 그 위에서는 국소 초음속 영역과 충격파가 생겨 유동이 더 이상 등엔트로피·비회전이 아니게 된다. 국소 마하수가 1이 되는 지점의 압력계수는 등엔트로피 관계에서 정확히 구할 수 있다.

Cp,crit=2γMcr2[(2+(γ−1)Mcr2γ+1)γ/(γ−1)−1]C_{p,\mathrm{crit}}=\frac{2}{\gamma M_{cr}^2} \left[\left(\frac{2+(\gamma-1)M_{cr}^2}{\gamma+1}\right)^{\gamma/(\gamma-1)}-1\right]

여기에 프란틀-글라우어트로 늘린 최소 압력계수 곡선 Cp,min⁡(M∞)=Cp0,min⁡/βC_{p,\min}(M_\infty)=C_{p0,\min}/\beta 를 겹쳐 그려 두 곡선의 교점을 읽으면 그 익형의 임계 마하수 추정치가 나온다. 이것이 교과서의 고전적 McrM_{cr} 추정법이며, 보정식의 유효 범위를 보정식 자신으로 정하는 재미있는 절차다.2

McrM_{cr} 을 넘으면 조종면 효율이 떨어지고 항력이 급증하는 이른바 항력 발산이 시작되므로, 결국 프란틀-글라우어트는 자신이 쓸모없어지는 지점을 알려주는 데까지가 자기 몫이다.

5. 고차 보정 — 카르만-첸과 라이톤[편집]

프란틀-글라우어트의 근본 약점은 국소 압력에 무관하게 자유류 마하수만으로 보정한다는 것이다. 흡입 피크가 큰 곳일수록 국소 마하수가 높아 압축성이 더 세게 작용해야 하는데, 그 정보가 식에 들어 있지 않다.

카르만-첸 규칙(Kármán–Tsien)은 압력-밀도 관계를 접선 직선으로 근사하는 이른바 접선기체(tangent gas) 가정에서 유도되며, Cp0C_{p0} 자체가 분모에 들어가 국소 압력 정보를 반영한다.

Cp=Cp0β+[M∞21+β]Cp02C_p=\frac{C_{p0}}{\displaystyle\beta+\left[\frac{M_\infty^2}{1+\beta}\right]\frac{C_{p0}}{2}}

라이톤 보정(Laitone)은 프란틀-글라우어트를 자유류가 아닌 국소 마하수로 평가하는 발상에서 나오며, 비열비 γ\gamma 가 명시적으로 등장한다.

Cp=Cp0β+[M∞2(1+γ−12M∞2)2β]Cp0C_p=\frac{C_{p0}}{\displaystyle\beta+\left[\frac{M_\infty^2\left(1+\frac{\gamma-1}{2}M_\infty^2\right)}{2\beta}\right]C_{p0}}

세 식은 Cp0→0C_{p0}\to0 이나 M∞→0M_\infty\to0 극한에서 서로 일치하고, 흡입 피크가 깊어질수록 갈라진다. 흡입면(Cp0<0C_{p0}<0)에서 분모를 비교하면 순서가 대수적으로 정해지는데, 프란틀-글라우어트가 가장 낙관적(압축성을 가장 약하게 봄), 라이톤이 가장 보수적, 카르만-첸이 그 사이다. 어느 쪽이 실측에 더 가까운지는 익형과 조건에 따라 갈리지만, 산업에서 가장 널리 채택된 것은 카르만-첸이고 익형 해석 코드 XFOIL 도 이것을 쓴다. 다만 어느 것도 충격파를 모른다는 사실은 변하지 않는다. 임계 마하수를 넘으면 셋 다 틀린다.

6. 초음속 — 애커렛 이론[편집]

M∞>1M_\infty>1 이면 선형화 방정식의 첫 항 부호가 뒤집힌다.

(M∞2−1) φxx−φyy=0(M_\infty^2-1)\,\varphi_{xx}-\varphi_{yy}=0

타원형이던 방정식이 쌍곡형으로 바뀌고, 이것은 사실상 파동방정식이다. 정보가 사방으로 퍼지지 않고 마하선을 따라서만 전파되며, 해는 f(x−λy)f(x-\lambda y) 꼴로 아래로 흐르는 파동이 된다(λ=M∞2−1\lambda=\sqrt{M_\infty^2-1}). 얇은 물체에 대한 결과가 애커렛(Ackeret) 이론이다.

Cp=2θM∞2−1C_p=\frac{2\theta}{\sqrt{M_\infty^2-1}}

θ\theta 는 표면의 국소 기울기(자유류에 대한 각도, 라디안)다. 아음속과 비교하면 세 가지가 인상적으로 다르다.

  • 압력이 국소 기울기만으로 결정된다. 상류의 다른 부분이 무엇을 하든 상관없다. 아음속에서는 형상 전체가 모든 점의 압력에 관여했다.
  • 분모가 M2−1\sqrt{M^2-1} 로 바뀐다. 아음속의 1−M2\sqrt{1-M^2} 와 형태가 짝을 이루며, 마하 1 양쪽에서 모두 발산한다.
  • 파항력이 생긴다. 비점성인데도 항력이 0이 아니다. 달랑베르의 역설이 초음속에서는 성립하지 않는다는 뜻으로, 초음속 유동의 근본적 특징이다.

정리하면 M=1M=1 근방은 아음속 이론과 초음속 이론이 양쪽에서 동시에 발산하며 만나는 구멍이다. 선형 이론이 원리적으로 다룰 수 없는 영역이 딱 거기다.

7. 천음속에서 왜 못 쓰나[편집]

수학적으로 답이 명확하다. 소섭동 전개에서 우리가 버린 항 중 가장 큰 것은 대략 M∞2(γ+1)(φx/U∞)φxxM_\infty^2(\gamma+1)(\varphi_x/U_\infty)\varphi_{xx} 인데, 이것을 버릴 수 있었던 근거는 계수 (1−M∞2)(1-M_\infty^2) 가 그보다 훨씬 크다는 것이었다. M∞→1M_\infty\to1 이면 그 근거가 사라진다. 두 항이 같은 크기가 되므로 비선형 항을 살려야 하고, 그 결과가 천음속 소섭동 방정식(TSD)이다.

[1−M∞2−M∞2(γ+1)φxU∞]φxx+φyy+φzz=0\left[1-M_\infty^2-M_\infty^2(\gamma+1)\frac{\varphi_x}{U_\infty}\right]\varphi_{xx} +\varphi_{yy}+\varphi_{zz}=0

대괄호 안의 부호가 해 자신에 의해 결정되므로, 이 방정식은 한 유동장 안에서 타원형 영역(아음속)과 쌍곡형 영역(국소 초음속)이 공존하는 혼합형이다. 그 경계에 충격파가 앉는다. 선형 이론의 세계관 — 해를 겹쳐 더할 수 있고 형상에 비례한다 — 이 통째로 무너지는 지점이다.

이 실패에는 값비싼 실무적 대가가 따른다. 천음속 플러터의 경계가 뚝 떨어지는 이른바 천음속 딥(transonic dip)은 익면 충격파가 구조 진동에 따라 앞뒤로 움직이며 공기력 위상을 바꾸는 데서 오는데, 프란틀-글라우어트로 압축성을 넣은 선형 공기력 모델은 이 골짜기를 아예 예측하지 못한다. 그리고 하필 그 골짜기가 여객기 순항 마하수 근처다. 천음속 영역에서 전산유체역학이 대안이 아니라 필수인 이유다.

8. 그런데도 여전히 쓴다[편집]

이렇게 한계가 뚜렷한데도 프란틀-글라우어트는 죽지 않았다.

  • 선형 공력 코드의 표준 부속. 와류 격자법·패널법·더블릿 격자법은 근본적으로 라플라스 방정식 풀이라, 압축성을 넣는 유일하게 값싼 길이 괴테르트 신축이다. 항공탄성학 실무에서 마하수별 공기력 영향계수 행렬을 만드는 방식이 정확히 이것이다.
  • 풍동 데이터 환산. 저속 풍동에서 얻은 데이터를 실제 비행 마하수로 옮길 때 1차 보정으로 쓴다.
  • 개념 설계 단계. 형상 후보 수백 개를 훑는 단계에서 CFD를 돌릴 수는 없다. 나눗셈 하나로 마하수 효과의 경향을 얻는 것은 충분히 남는 장사다.
  • 게임과 실시간 시뮬레이션. 프레임당 수십 마이크로초 안에 공력을 계산해야 하는 물리 엔진이나 비행 게임에서는 마하수별 공력 테이블을 다 들고 있기 부담스러울 때가 많다. 저속 계수 하나에 1/β1/\beta 를 곱하는 것으로 “빠를수록 양력이 커진다”는 정성적 거동을 공짜로 얻는다. 다만 β\beta 가 0으로 가는 것을 막기 위해 마하 0.7~0.8 근처에서 클램프하는 것이 필수다. 안 그러면 플레이어가 급강하할 때 양력이 발산해 기체가 우주로 사라진다.3

요컨대 프란틀-글라우어트는 정확한 답을 주는 도구가 아니라, 압축성이 개입한다는 사실을 가장 싸게 반영하는 도구다. 어디까지 믿을지를 아는 사람이 쓰면 유용하고, 모르고 쓰면 임계 마하수 위에서 조용히 틀린 숫자를 준다.4

9. 관련 문서[편집]

10. Footnotes[편집]

  1. 프란틀이 1920년대 초 강의에서 이 착상을 내놓았고, 글라우어트가 1928년 논문으로 정식화하면서 두 사람의 이름이 붙었다. 초음속 쪽의 대응물은 그보다 앞선 1925년 애커렛의 결과이고, 3차원 형상 대응 관계를 바로잡은 괴테르트의 정리는 1940년경이다. 즉 “압축성 보정”이라는 한 덩어리로 배우는 내용이 실제로는 20년에 걸쳐 서로 다른 사람들이 채워 넣은 조각들이다. ↩

  2. 자기 유효 범위를 자기가 알려주는 근사식은 흔치 않다. 대개는 “언제부터 틀리는지”를 실험이나 상위 이론이 따로 알려줘야 하는데, 이 경우는 발산하는 근사식과 정확한 등엔트로피 관계가 만나는 지점이 곧 유통기한이 된다. 우아하다고 해야 할지, 뻔뻔하다고 해야 할지. ↩

  3. 실제 비행 시뮬레이터는 대개 마하수별 계수 테이블을 보간해서 쓴다. 프란틀-글라우어트는 테이블을 만들 예산이 없는 프로젝트, 또는 저속 데이터밖에 없는 기체를 급조로 모델링할 때의 임시방편에 가깝다. 클램프를 잊은 코드가 마하 0.99에서 양력을 7배로 만드는 사고는 아케이드 비행 게임의 오랜 전통이다. ↩

  4. 천음속 비행기 주위에 생기는 원뿔 모양 수증기 구름을 “프란틀-글라우어트 특이점”이라 부르는 관행이 인터넷에 널리 퍼져 있는데, 이 이름은 부정확하다. 그 구름은 팽창으로 국소 온도가 이슬점 아래로 내려가 수증기가 응결한 것이고, 충격파의 위치와 반드시 일치하지도 않는다. 발산하는 수식 이름이 붙었다는 이유만으로 사진이 더 멋있어 보이는 효과는 인정한다. ↩