기체동역학

전산유체역학 유체역학 마지막 수정: 2026-07-20 04:08:27

1. 개요[편집]

기체동역학 Gas Dynamics
핵심 무차원수	마하수 $Ma$
지배 원리	등엔트로피 관계, 에너지 보존
대표 현상	노즐 유동, 팽창파, 초킹
상위 분야	압축성 유동
응용	로켓 노즐, 초음속기, 풍동

공기도 빨라지면 자기가 스프링인 걸 기억해낸다.

기체동역학(gas dynamics)은 밀도 변화가 유동의 본질을 지배하는 고속 기체 흐름을 열역학과 유체역학의 결합으로 다루는 학문이다. 저속에서는 공기를 밀도 일정한 액체처럼 취급해도 무방하지만, 속도가 음속에 근접하면 압력·밀도·온도가 한 몸으로 얽혀 움직이며 완전히 다른 세계가 열린다. 기체동역학은 바로 그 세계의 문법을 다룬다.

압축성 유동과 상당 부분 겹치는 주제이지만, 관점이 조금 다르다. 압축성 유동 문서가 충격파 불연속과 리만 솔버 같은 수치 기법에 무게를 둔다면, 기체동역학은 노즐·디퓨저·풍동 같은 장치에서 유동이 어떻게 가속·감속되는지를 1차원 열역학 관계식으로 예측하는 고전 공학의 결이 강하다. 로켓 엔진, 초음속 흡입구, 가스터빈 유로 설계의 언어가 여기서 나온다.

2. 마하수와 유동 영역[편집]

기체동역학의 모든 이야기는 마하수 $Ma = u/a$ 에서 시작한다. $a$ 는 국소 음속이고, 이상기체에서는 온도만의 함수다.

a = \sqrt{\gamma R T}

$\gamma$ 는 비열비(공기 1.4), $R$ 은 기체상수, $T$ 는 절대온도다. 마하수 1을 기준으로 유동은 인격이 갈린다. 아음속( $Ma<1$ )에서는 교란이 상류로 거슬러 전파되어 유동이 장애물의 존재를 “미리” 알지만, 초음속( $Ma>1$ )에서는 교란이 하류로만 흘러 유동이 장애물을 만나고 나서야 반응한다. 이 상류 전파 가능 여부의 차이가 노즐 형상 설계부터 충격파 발생까지 모든 것을 결정한다.¹

영역	마하수	성격
아음속	$Ma < 0.8$	교란 상류 전파, 매끄러움
천음속	$0.8 \sim 1.2$	국소 초음속 + 충격파 공존
초음속	$1.2 \sim 5$	경사 충격파, 팽창파 지배
극초음속	$Ma > 5$	실제기체·해리 효과

3. 등엔트로피 유동[편집]

기체동역학의 뼈대는 등엔트로피(isentropic) 유동 가정이다. 충격파가 없고 마찰·열전달을 무시할 수 있는 구간에서는 엔트로피가 보존되며, 정체점(stagnation) 상태량과 국소 상태량이 마하수만의 함수로 깔끔하게 연결된다.

\frac{T_0}{T} = 1 + \frac{\gamma - 1}{2} Ma^2

압력비와 밀도비는 여기에 각각 $\gamma/(\gamma-1)$ , $1/(\gamma-1)$ 지수를 씌우면 나온다. 이 한 줄짜리 관계식들이 실무에서 갖는 위력은 엄청나다. 풍동에서 정체 조건과 마하수만 알면 시험부의 온도·압력이 표 한 장으로 튀어나오고, 로켓 노즐 출구 조건도 손계산으로 잡힌다. 베르누이 방정식의 압축성 확장판이라고 봐도 크게 틀리지 않는다.

4. 노즐 유동과 초킹[편집]

기체동역학의 하이라이트는 노즐이다. 단면적 $A$ 와 속도 $u$ 의 관계는 유명한 면적-속도 관계식으로 요약된다.

\frac{dA}{A} = (Ma^2 - 1)\frac{du}{u}

부호에 주목하자. 아음속( $Ma<1$ )에서는 괄호가 음수라 단면이 좁아져야 유동이 빨라진다. 여기까지는 정원 호스를 손으로 눌러본 사람의 직관과 같다. 그런데 초음속( $Ma>1$ )에서는 괄호가 양수로 뒤집혀 단면이 넓어져야 더 빨라진다. 이 직관 배반이 로켓 노즐이 목(throat) 뒤에서 나팔처럼 벌어지는 이유다.²

수축-확대 노즐(라발 노즐)에서 목의 마하수가 1에 도달하면 초킹(choking) 이 발생한다. 하류 압력을 아무리 낮춰도 질량유량이 더는 늘지 않는 상태로, 목이 유량의 병목을 물리적으로 잠가버린다. 제트·로켓 엔진 설계의 기본 상식이자, 초심자가 “왜 밸브를 더 열었는데 유량이 그대로죠?”라며 처음 부딪히는 벽이기도 하다.

5. 초음속 유동[편집]

초음속 유동(supersonic flow)은 $Ma>1$ , 즉 유속이 국소 음속을 넘는 흐름으로, 기체동역학이 가장 극적으로 드러나는 무대다. 초음속의 핵심 성질은 하나로 요약된다 — 교란이 상류로 전파되지 못한다. 그래서 초음속 물체는 자신이 다가온다는 신호를 앞쪽 공기에 보내지 못하고, 공기는 물체가 코앞에 닿아서야 비명을 지른다. 그 비명이 바로 충격파이며, 물체를 지나 지면에 닿으면 소닉붐으로 들린다.

초음속 유동에서 방향을 바꾸는 방식은 두 가지뿐이다. 흐름이 자기 안쪽으로 꺾이면(압축) 경사 충격파(oblique shock)가 서면서 압력이 계단처럼 뛰고 마하수가 떨어진다. 반대로 흐름이 바깥쪽으로 꺾이면(팽창) 프란틀-마이어 팽창파(Prandtl–Meyer expansion fan)라는 매끄러운 부채꼴 파동을 통해 연속적으로 가속하며 압력이 떨어진다. 압축은 계단으로, 팽창은 부드럽게 — 이 비대칭은 열역학 제2법칙이 강제하는 자연의 규칙이다.³ 충격파의 상세한 랭킨-위고니오 관계와 수치 처리는 압축성 유동 문서가 자세히 다룬다.

초음속 CFD에서는 정보가 특성선(characteristic line)을 따라 유한한 속도로 전파된다는 쌍곡형 성질 덕분에, 순수 초음속 영역은 하류로 진행하며 순차적으로 풀 수 있다. 이 아이디어를 노즐·초음속기 설계에 응용한 고전 기법이 특성곡선법이다.

6. 수치해석과 실제기체 효과[편집]

기체동역학을 전산유체역학으로 옮기면 몇 가지 대가를 치른다. 압축성 명시적 솔버는 CFL 조건이 $u+a$ 기준이라 시간 간격이 잘게 쪼개지고, 충격파 근처에서는 매끄러움을 전제한 고차 도식이 진동을 일으켜 TVD·WENO 같은 억제 장치가 필요하다. 저마하 영역에서는 음속이 유속보다 훨씬 커서 방정식이 강성(stiff)해지는 아이러니도 있다.

극초음속( $Ma>5$ )에 가면 이상기체 가정마저 무너진다. 충격파 뒤 온도가 수천 K에 이르러 산소·질소가 해리·전리되고, 화학 반응과 복사 열전달이 개입하며 $\gamma$ 가 더 이상 상수가 아니게 된다. 재진입체와 극초음속 비행체 해석이 유독 어려운 이유가 여기 있다. 이쯤 되면 “일단 돌려”라는 CFD 국룰도 화학 반응 모델 앞에서 잠시 겸손해진다.

7. 관련 문서[편집]

이 “상류로 소식을 전할 수 있느냐”가 아음속과 초음속을 가르는 진짜 경계다. 아음속 새는 다가오는 차 소리를 미리 듣고 피하지만, 초음속 총알은 자기 소리보다 먼저 도착한다. 자연이 정보에 속도 제한을 걸어둔 셈. ↩
노즐이 목 뒤에서 벌어지는 걸 처음 본 사람은 대개 “좁혀야 빨라지는 거 아냐?”라고 되묻는다. 맞다, 아음속에서는. 초음속에서는 정반대다. 스웨덴 엔지니어 구스타프 드 라발이 증기터빈 노즐에서 이 원리를 실용화해 라발 노즐이라는 이름이 붙었다. ↩
팽창 충격파는 엔트로피를 감소시켜 열역학 제2법칙을 위반하므로 자연에 존재하지 않는다. 수학적으로는 방정식이 두 해를 다 뱉지만, 우주는 엔트로피가 증가하는 쪽 하나만 허락한다. 대학원생이 코드에서 팽창 충격파를 발견하면 그건 버그다. ↩