극초음속 유동

유체역학 전산유체역학 계산물리 마지막 수정: 2026-08-06 05:06:24

1. 개요[편집]

극초음속 유동 Hypersonic flow
관례적 기준	$M \gtrsim 5$ (숫자보다 물리가 기준)
특징	얇은 충격층, 엔트로피 층, 점성 상호작용
기체 상태	진동 여기 · 해리 · 전리, $\gamma \neq$ 상수
대표 무대	대기 재진입, 스크램제트, 활공체
연속체 한계	크누센수 $Kn \gtrsim 0.01$에서 붕괴

마하 5는 마하 2보다 두 배 반 빠른 게 아니다. 아예 다른 과목이다.

극초음속 유동(hypersonic flow)은 관례상 마하수 $M \gtrsim 5$ 인 유동을 가리키되, 단순히 더 빠른 초음속이 아니라 얇은 충격층·점성 상호작용·고온 실기체 효과 같은 질적으로 새로운 물리가 지배하기 시작하는 영역이다. 기준이 5인 것에 물리적 필연성은 없다. 실제 기준은 “이상기체 등엔트로피 관계식으로 손계산하던 초음속 유동의 문법이 언제 깨지는가”이고, 그 시점은 물체 형상과 고도에 따라 마하 3에서 시작될 수도, 마하 8까지 미뤄질 수도 있다.¹

무엇이 깨지는지가 이 문서의 본론이다. 극초음속에서는 (1) 공기역학이 화학과 얽히고, (2) 비점성 영역과 경계층의 구분이 흐려지며, (3) 연속체 가정 자체가 고고도에서 무너진다. 압축성 유동의 깔끔한 표들이 하나씩 무효화되는 과정이라고 봐도 된다.

2. 얇은 충격층과 엔트로피 층[편집]

수직 충격파를 지나는 밀도비는 $M \to \infty$ 극한에서 유한한 값에 수렴한다.

\frac{\rho_2}{\rho_1} \to \frac{\gamma+1}{\gamma-1}

$\gamma=1.4$ 면 6이고, 고온에서 해리가 일어나 유효 비열비가 1.2 근처로 떨어지면 10을 넘는다. 밀도가 이렇게 압축되면 같은 질량유량을 통과시키는 데 필요한 단면이 얇아지므로, 이탈 충격파가 물체 표면에 바짝 달라붙는다. 이것이 얇은 충격층(thin shock layer)이며, 구형 앞머리의 충격 이탈거리 $\Delta/R_n$ 은 대략 밀도비의 역수에 비례한다. 재진입체 사진에서 충격파가 기체 표면과 거의 붙어 보이는 건 착시가 아니다.

충격파가 휘어 있으면 유선마다 통과하는 충격파의 세기가 다르고, 따라서 유선마다 엔트로피가 다르다. 정체 유선 근처는 거의 수직 충격을 지나 엔트로피가 가장 높고, 바깥으로 갈수록 약한 경사충격을 지나 낮다. 이 엔트로피 구배 층이 엔트로피 층(entropy layer)이며, 크로코 정리에 따라 여기에는 비점성인데도 와도가 존재한다. 물체를 따라 자라난 경계층은 결국 이 엔트로피 층을 삼키게 되는데(entropy layer swallowing), 그 순간 경계층 바깥 경계조건이 자유류가 아니라 “엔트로피가 높고 밀도가 낮은 층”으로 바뀌면서 표면 마찰과 열유속 예측이 통째로 달라진다. 뭉툭한 코를 가진 비행체의 후방 열유속을 코 반경 없이 계산하면 틀리는 이유다.

3. 점성 상호작용과 공력가열[편집]

극초음속 경계층은 두껍다. 점성 산일이 유동의 운동에너지를 열로 바꿔 경계층 내부 온도를 크게 올리고, 온도가 오르면 밀도가 떨어져 같은 질량유량을 흘리는 데 더 두꺼운 층이 필요하기 때문이다. 층류 경계층 두께는 대략 다음처럼 스케일한다.

\frac{\delta}{x} \sim \frac{M_\infty^2}{\sqrt{Re_x}}

$M^2$ 이 붙어 있다는 게 핵심이다. 마하 10이면 같은 레이놀즈수의 저속 유동보다 경계층이 100배 두꺼워지는 셈이고, 이쯤 되면 경계층이 자기 두께로 바깥 비점성 유동을 밀어내 압력 분포를 바꾼다. 원인과 결과가 서로를 물고 도는 점성 상호작용(viscous interaction)이며, 그 세기를 재는 지표가 상호작용 파라미터다.

\bar{\chi} = \frac{M_\infty^3\sqrt{C_\infty}}{\sqrt{Re_x}}

$\bar{\chi}$ 가 커지면 날카로운 평판 앞전에서도 압력이 자유류 값보다 훨씬 높게 나온다. 얇은 익형이 좋다는 초음속의 상식이 극초음속에서 통하지 않는 첫 번째 이유다.

두 번째 이유는 열이다. 정체점 대류 열유속은 널리 쓰이는 상관식에서 $\dot q \propto \sqrt{\rho_\infty / R_n}\,V_\infty^3$ 로 스케일하는데, 코 반경의 제곱근에 반비례한다. 그래서 재진입체는 일부러 뭉툭하게 만든다(앨런-에거스의 뭉툭체 원리). 뭉툭한 코가 세운 강한 이탈 충격파는 운동에너지의 대부분을 물체가 아니라 주변 공기를 데우는 데 써버리고, 물체는 그 잔돈만 받는다. 날카롭게 만들면 항력은 줄지만 앞전이 녹는다. 이 열을 견디는 설계가 열보호시스템 해석이며, 삭마형과 재사용형(타일) 두 계열의 절충 지점이 여기서 정해진다.

경계층의 난류 천이 예측도 판돈이 크다. 난류 열유속은 층류의 서너 배까지 뛰므로 천이 위치 하나가 TPS 두께 전체를 좌우한다. 그런데 극초음속에서는 저속의 토미엔-슐리히팅 파 대신 음향 모드(맥 제2모드)가 지배 불안정이 되어, 저속 경험식이 아예 적용되지 않는다.

4. 고온 실기체 효과[편집]

정체 온도는 $T_0/T = 1+\frac{\gamma-1}{2}M^2$ 로 폭주한다. 고도 30 km 부근(약 227 K)에서 마하 10이면 이상기체 계산상 $T_0$ 가 4700 K를 넘는다. 실제로는 그렇게까지 안 오르는데, 이유가 바로 실기체 효과다 — 에너지의 상당 부분이 온도를 올리는 대신 분자 내부 자유도와 화학결합을 깨는 데 쓰이기 때문이다. 순서는 대략 이렇다.

온도 대역	일어나는 일	결과
800 K 이상	진동 자유도 여기	비열 증가, $\gamma$ 감소
2000~4000 K	산소 해리 $\mathrm{O_2 \to 2O}$	유효 $\gamma$ 급락, 충격층 압축
4000~9000 K	질소 해리 $\mathrm{N_2 \to 2N}$ , NO 생성	다종 혼합기
9000 K 이상	전리, 자유전자 발생	플라즈마 층 형성

이렇게 되면 “ $\gamma$ 는 1.4”라는 전제 위에 세워진 등엔트로피 표·충격파 표가 전부 폐기된다. 해석은 다종 화학종 보존식에 반응 소스항을 붙인 연립계로 넘어가고, 반응 속도는 아레니우스 방정식 형태로 모델링한다(반응속도론).

여기서 결정적인 무차원수가 담쾰러 수 $Da$ — 유동 체류 시간과 화학 반응 시간의 비다. $Da \gg 1$ 이면 화학이 순식간에 따라잡아 국소 평형 유동, $Da \ll 1$ 이면 반응할 틈이 없는 동결 유동이고, 재진입 궤적의 대부분이 하필 $Da \sim 1$ 인 비평형 영역에 걸린다. 평형·동결이라는 두 극한은 각각 계산이 쉬운데, 실제 비행은 둘 사이 어중간한 곳에 있는 셈이다.

열적 비평형도 따로 붙는다. 병진·회전 온도 $T$ 와 진동 온도 $T_v$ 가 서로 다르게 완화되므로, 파크(Park)의 2온도 모델처럼 온도를 두 개 이상 들고 다니면서 반응 속도를 $\sqrt{T T_v}$ 같은 조합으로 평가한다. 전리가 시작되면 자유전자가 만든 플라즈마 층이 전파를 반사·흡수해 통신 두절(blackout) 이 발생한다. 아폴로급 재진입에서 수 분간 교신이 끊겼던 것이 이 현상이며, 전자 밀도가 만드는 차단 주파수 문제라 해석은 플라즈마 시뮬레이션 영역으로 넘어간다.

5. 마하수 독립 원리와 뉴턴 이론[편집]

역설적이게도, 물리가 복잡해지는 만큼 어떤 것들은 오히려 단순해진다. 마하수 독립 원리(Oswatitsch)에 따르면 $M_\infty$ 가 충분히 크면 뭉툭체 주변의 압력계수·항력계수·충격파 형상이 $M_\infty$ 에 거의 무관해진다. 자유류 압력 $p_\infty$ 가 충격 후 압력에 비해 무시할 만큼 작아져, 남는 것은 자유류 운동량뿐이기 때문이다.

그래서 뉴턴이 17세기에 (틀린 전제로) 제안했던 충돌 모형이 이 극한에서 되살아난다. 유체 입자가 표면에 부딪혀 법선 운동량을 전부 잃는다고 가정하면

C_p = 2\sin^2\theta

가 나온다( $\theta$ 는 자유류와 표면이 이루는 각). 아음속에서는 완전히 엉터리인 이 식이 극초음속 뭉툭체에서는 놀랍도록 잘 맞는다. 실무에서는 정체점 값을 수직충격 뒤 정체압으로 교정한 수정 뉴턴 이론 $C_p = C_{p,\max}\sin^2\theta$ 를 쓰며, 형상 최적화 초기 단계에서 CFD 없이 압력 분포를 훑는 도구로 지금도 쓰인다. 여기에 극초음속 상사 파라미터 $K = M_\infty\theta$ 까지 더하면, 마하수와 두께비가 다른 형상들을 하나의 곡선으로 묶을 수 있다.

6. 수치해석 — 무엇이 부러지는가[편집]

극초음속 CFD는 초음속 유동 코드에 화학을 얹으면 되는 일이 아니다. 다음 네 가지가 순서대로 발목을 잡는다.

충격 포착 도식의 병리. 강한 이탈 충격파를 격자선에 정렬시켜 계산하면 정체선 앞에 혹 같은 비물리적 구조가 자라난다(carbuncle). 접촉 불연속을 잘 보존하는 정확한 리만 솔버일수록 잘 걸리는 역설적인 병이라, 충격 근처에서만 소산이 큰 솔버로 블렌딩하거나 격자를 충격면과 일부러 어긋나게 까는 식으로 회피한다. 정체점 열유속은 격자 품질에 극도로 민감해서, 벽면 셀 레이놀즈수를 충분히 낮추지 않으면 열유속이 수십 % 단위로 틀린다.
화학 소스항의 강성. 반응 시간 척도가 유동 시간 척도보다 수 자릿수 짧아 방정식이 강성 방정식이 된다. 명시적 시간적분은 화학이 정한 시간 간격에 묶여 사실상 불가능하므로, 소스항만 암시적으로 처리하는 점-암시(point-implicit) 기법이나 연산자 분리가 표준이다.
복사 열전달. 진입 속도가 10 km/s를 넘으면 충격층 기체 자체가 빛나며 복사로 물체를 데운다. 복사 전달 방정식을 유동과 연성해 풀어야 하고, 스펙트럼 선을 다 담으면 비용이 폭발해 밴드 모델로 압축한다.
연속체 가정의 붕괴. 고고도로 갈수록 평균자유행로가 길어져 크누센수가 커진다. $Kn \lesssim 0.01$ 까지가 연속체, 0.01~~0.1은 슬립, 0.1~~10은 천이, 그 위는 자유분자 영역이다. 나비에-스토크스가 유효하지 않은 영역에서는 볼츠만 방정식 기반의 직접 시뮬레이션 몬테카를로로 갈아타야 하며, 하나의 재진입 궤적이 자유분자 → 천이 → 연속체를 전부 통과하기 때문에 코드 두 벌을 이어 붙이는 하이브리드 해석이 흔하다.²

여기에 지상 시험의 한계가 겹친다. 극초음속 풍동은 마하수와 레이놀즈수는 맞춰도 정체 엔탈피(즉 실제 화학 상태) 를 동시에 맞추기가 극히 어렵고, 팽창관·충격관 시설은 시험 시간이 밀리초 단위다. 결국 검증 및 확인의 실험 근거가 만성적으로 부족한 분야이며, 비행 시험 데이터 한 세트의 가치가 유난히 높은 이유이기도 하다. 스크램제트 개발이 유독 더디게 진행되는 배경에도 이 검증 공백이 있다.³

7. 관련 문서[편집]

앤더슨의 교과서가 못 박듯이 “극초음속은 마하수가 아니라 현상으로 정의된다”. 실제로 뭉툭체에서는 마하 3~4부터 실기체 효과가 나타나기 시작하고, 아주 날카로운 쐐기라면 마하 8에서도 이상기체 표가 그럭저럭 맞는다. 5는 교과서 목차를 나누기 좋은 숫자일 뿐이다. ↩
그래서 재진입 해석은 코드 한 벌로 끝나는 법이 없다. DSMC 코드, 화학 비평형 CFD 코드, 복사 코드, 열구조 코드가 순서대로 붙고, 각 접합면에서 물리량을 넘겨주다 보면 어디서 오차가 났는지 아무도 모르게 된다. “일단 돌려”의 난이도 최상급 버전. ↩
X-43A가 2004년에 마하 9.6대, X-51A가 마하 5대에서 스크램제트 연소 비행을 실증했지만, 스크램제트 연소 지속 시간은 각각 10초 남짓·200초 남짓이었다. 실험실에서 밀리초, 비행에서 수 분 — 이 사이의 간극을 메우는 것이 이 분야의 영원한 현생이다. ↩