라발 노즐

1. 개요[편집]

라발 노즐 de Laval nozzle
형상	수축부 → 목(throat) → 확대부
핵심 관계	면적-속도 관계, 면적-마하수 관계
임계 압력비	$p^*/p_0 = 0.5283$ ($\gamma = 1.4$)
운전 모드	벤투리 · 초킹 · 내부충격 · 과팽창 · 설계점 · 저팽창
이름의 주인	구스타프 드 라발 (1845~1913)

좁히면 빨라진다는 상식이 목에서 뒤집힌다. 그 다음부터는 넓혀야 빨라진다.

라발 노즐(de Laval nozzle)은 수축부와 확대부를 목(throat)으로 이어 붙여, 아음속으로 들어온 기체를 목에서 음속으로 만들고 확대부에서 초음속으로 가속하는 수축-확대 노즐이다. 로켓 엔진, 초음속 풍동, 증기터빈 노즐, 초음속 흡입구의 역방향 설계까지 — 기체를 음속 너머로 밀어붙이는 거의 모든 장치가 이 형상을 쓴다. 다른 방법이 없어서 쓰는 게 아니라, 준1차원 등엔트로피 유동에서는 이것 말고 방법이 없다.

기체동역학이 등엔트로피 관계식 전반을, 초음속 유동이 쌍곡형 유동의 물리와 수치해석을 다룬다면, 이 문서는 노즐이라는 장치 하나에 초점을 맞춘다. 부호가 뒤집히는 이유, 초킹이 무엇을 잠그는지, 배압을 낮출 때 노즐 안에서 무슨 일이 순서대로 벌어지는지가 본론이다.

2. 부호가 뒤집히는 지점[편집]

준1차원 정상 등엔트로피 유동에서 연속·운동량·음속 관계를 정리하면 면적-속도 관계가 나온다.

$$\frac{dA}{A} = (M^2 - 1)\,\frac{du}{u}$$

이 한 줄이 노즐 형상의 전부다. 아음속($M<1$)에서는 괄호가 음수라 $dA<0$일 때 $du>0$ — 좁혀야 빨라진다. 호스 끝을 손가락으로 눌러본 사람의 직관 그대로다. 그런데 $M=1$을 넘는 순간 괄호가 양수로 바뀌어 넓혀야 빨라진다. 밀도 감소율이 속도 증가율을 앞지르기 때문이다. 비압축성 세계에서는 $\rho$가 상수라 $Au$만 지키면 됐지만, 압축성에서는 $\rho A u$가 보존량이고 초음속에서는 $\rho$가 워낙 빠르게 떨어져 $A$까지 키워줘야 수지가 맞는다.¹

그리고 $M=1$은 $dA=0$인 곳, 즉 면적이 극값을 갖는 단면에서만 성립할 수 있다. 목이 아닌 곳에서 음속에 도달하는 매끄러운 해는 존재하지 않는다. 반대로 목에서 반드시 $M=1$이 되는 것도 아니다 — 압력비가 부족하면 목에서 최대 마하수를 찍고 다시 감속하는 벤투리 관 모드로 남는다.

3. 면적-마하수 관계와 두 개의 해[편집]

음속 단면의 면적을 $A^*$로 잡으면, 등엔트로피 유동에서 임의 단면의 면적비는 그 단면의 마하수만으로 결정된다.

$$\frac{A}{A^*} = \frac{1}{M}\left[\frac{2}{\gamma+1}\left(1+\frac{\gamma-1}{2}M^2\right)\right]^{\frac{\gamma+1}{2(\gamma-1)}}$$

이 함수는 $M=1$에서 최솟값 1을 갖고 양쪽으로 단조 증가한다. 즉 주어진 $A/A^*$에 대해 아음속 해와 초음속 해가 하나씩 나온다. 확대비 $A_e/A^* = 4$인 노즐이라면 $\gamma=1.4$에서 $M \approx 0.147$과 $M \approx 2.94$가 둘 다 방정식을 만족한다. 노즐 형상만 보고는 어느 쪽이 실현될지 알 수 없고, 배압(back pressure) $p_a$가 정한다. 형상은 후보를 두 개 제시할 뿐이며 선택은 하류가 한다는 것 — 이게 노즐 해석의 핵심 감각이다.

목에서 초킹이 걸리면 임계 상태량은 정체 상태량만의 함수로 고정된다.

$$\frac{p^*}{p_0} = \left(\frac{2}{\gamma+1}\right)^{\frac{\gamma}{\gamma-1}}, \qquad \frac{T^*}{T_0} = \frac{2}{\gamma+1}$$

$\gamma=1.4$에서 $p^*/p_0 = 0.5283$, $T^*/T_0 = 0.8333$이다. 대기로 뿜는 노즐이라면 정체압이 대기압의 약 1.89배만 넘어도 목은 이미 초킹된다는 뜻이라, 압축공기 배관에서 이 조건은 생각보다 훨씬 자주 걸린다.

4. 초킹 — 하류가 더 이상 상류에 말을 걸지 못한다[편집]

목이 음속에 도달하면 질량유량은 다음 값에 못박힌다.

$$\dot m = \frac{p_0 A^*}{\sqrt{T_0}}\sqrt{\frac{\gamma}{R}}\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}}$$

배압을 더 낮춰도 $\dot m$은 꿈쩍하지 않는다. 이것이 초킹(choking) 이다. 이유는 정보 전파에 있다. 목의 하류 특성속도가 $u - a = 0$ 이하가 되어, 하류에서 만든 압력 신호가 목을 거슬러 올라가지 못한다. 상류는 하류가 얼마나 낮아졌는지 알 방법이 없으니 유량을 바꿀 근거도 없다.² 초킹된 유량을 늘리는 방법은 셋뿐이다 — $p_0$를 올리거나, $T_0$를 낮추거나, $A^*$를 키우거나. 밸브를 더 여는 건 목록에 없다.

이 성질은 계측에서 오히려 축복이다. 임계 유량 벤투리(critical flow venturi)는 초킹된 목의 유량이 하류 조건과 완전히 무관하다는 점을 이용해, 상류 $p_0$·$T_0$만 재고도 유량을 표준급 정확도로 확정한다. 벤투리 관의 아음속 사용법과는 원리가 아예 다르다.

5. 배압을 낮춰가며 보는 운전 모드[편집]

정체 조건을 고정하고 배압 $p_a$만 계속 낮춘다고 하자. 노즐은 다음 순서를 밟는다.

모드	목	확대부	배출
벤투리	$M<1$	감속(디퓨저)	아음속
첫 초킹	$M=1$	감속	아음속
내부 수직충격	$M=1$	초음속 → 충격 → 아음속	아음속
과팽창	$M=1$	전 구간 초음속	경사충격으로 압축
설계점	$M=1$	전 구간 초음속	$p_e = p_a$, 파동 없음
저팽창	$M=1$	전 구간 초음속	팽창 부채꼴로 마저 팽창

주목할 구간은 세 번째다. 목이 초킹된 뒤 배압이 아직 초음속 해에 못 미치면, 확대부 안에서 유동이 초음속으로 가속하다가 수직충격을 만나 아음속으로 떨어지고 나머지 확대부를 디퓨저처럼 쓴다. 충격 위치는 “충격 뒤 아음속 유동이 출구에서 정확히 $p_a$가 되도록” 스스로 자리를 잡는다. 배압을 더 낮추면 충격이 하류로 밀려나고, 결국 출구면에 도달한 뒤 노즐 밖으로 튀어나가 경사충격 계열이 된다.

과팽창이 심하면 확대부 안에서 유동박리가 일어난다. 대략 $p_e/p_a \approx 0.4$ 아래에서 벽면 박리가 시작된다는 경험 기준(서머필드 기준)이 널리 쓰이는데, 이때 발생하는 비대칭 측력은 로켓 노즐 구조를 실제로 부순다. 지상 시험에서 대확대비 노즐이 가장 무서운 순간이 시동·정지 과도기인 이유다.

6. 추력과 고도 보상[편집]

로켓 노즐의 추력은 운동량 항과 압력 항의 합이다.

$$F = \dot m\,u_e + (p_e - p_a)A_e$$

$p_e = p_a$인 설계점에서 두 번째 항이 사라지고, 이때 주어진 $A_e$에 대해 추력이 최대가 된다. 문제는 $p_a$가 고도에 따라 변한다는 것이다. 지상 최적으로 설계하면 진공에서 저팽창이라 손해고, 진공 최적으로 설계하면 지상에서 과팽창이라 박리 위험까지 안는다. 1단 엔진은 대개 확대비를 낮게(예: 20 내외), 상단 엔진은 높게(100 이상) 잡는 절충으로 끝내며, 그래서 같은 엔진 계열도 지상용과 진공용 노즐 확대비가 다르다. 노즐 하나로 전 고도를 커버하려는 시도가 에어로스파이크 노즐이지만, 냉각과 무게 문제로 실용화는 여전히 더디다. 성능 지표로는 $F$ 자체보다 비추력이 더 자주 쓰인다.

7. 윤곽 설계와 수치해석[편집]

준1차원 이론은 면적 분포만 말해줄 뿐 벽면 곡선은 정해주지 않는다. 단순 원뿔 노즐은 출구 유동이 축에 대해 반각만큼 벌어져 있어 추력에 $\lambda = (1+\cos\alpha)/2$ 만큼의 발산 손실이 생긴다(반각 $15^\circ$면 약 1.7% 손실). 이 손실을 줄이려고 확대부를 종 모양으로 휘게 만든 것이 벨 노즐이고, 그 윤곽은 특성곡선법으로 설계한다. 목 직후 팽창부에서 프란틀-마이어 팽창으로 유동을 벌린 뒤, 벽을 다시 안쪽으로 꺾어 마하파를 상쇄시켜 출구에서 균일 평행류를 만드는 것이 표준 절차다.

CFD로 넘어가면 오일러 방정식 (비점성 유동)을 유한체적법으로 푸는 것이 기본이고, 내부 충격을 잡아야 하므로 보존형 이산화와 리만 솔버가 필수다. 준1차원 노즐 문제는 CFD 코드 검증의 고전 벤치마크이기도 하다 — 등엔트로피 해와 수직충격 해를 해석적으로 알기 때문에 격자 수렴과 충격 위치 오차를 정량 비교할 수 있어서, 검증 및 확인 절차에서 초급 문제로 자주 등장한다.³ 실제 엔진 해석에서는 여기에 경계층 배제두께 보정, 다종 화학 평형(연소 시뮬레이션), 벽면 재생냉각 열유속이 얹히면서 문제가 통째로 커진다.

8. 관련 문서[편집]

• 기체동역학 · 압축성 유동
• 초음속 유동 · 극초음속 유동
• 충격파 · 프란틀-마이어 팽창
• 특성곡선법 · 리만 솔버
• 오일러 방정식 (비점성 유동) · 유한체적법
• 유동박리 · 터보기계 해석
• 무차원수 · 검증 및 확인

9. Footnotes[편집]

1. 스웨덴 엔지니어 구스타프 드 라발이 1888년 증기터빈 노즐에 이 형상을 실용화하면서 이름이 붙었다. 정작 그는 초음속이라는 개념을 명시적으로 다루지 않았고, “이렇게 만드니 증기가 더 빨리 나가더라”는 실험적 발견에 가까웠다. 이론적 정리는 프란틀 학파가 20세기 초에 마무리했다. ↩

2. 그래서 초킹된 노즐 상류에 앉아 있는 압력계는 하류에서 벌어지는 참사를 끝내 모른다. 밸브를 아무리 열어도 유량이 안 늘어난다며 계측 장비를 의심하는 신입은 매년 재생산된다. 장비는 멀쩡하고 물리가 잠근 것이다. ↩

3. 이 문제의 잔인한 점은 답이 이미 종이에 있다는 것이다. 코드가 틀리면 변명할 여지가 없다. “격자가 다 했다”고 우기려 해도 격자를 조이면 해석해로 수렴해야 정상이라, 안 수렴하면 그건 격자가 아니라 플럭스 함수 문제다. ↩