초음속 유동

편집 역사 토론
유체역학 전산유체역학 시뮬레이션 마지막 수정: 2026-08-05 04:57:33

1. 개요[편집]

초음속 유동
Supersonic flow
정의역마하수 M > 1
방정식 성격정상 유동에서 쌍곡형(hyperbolic)
지배 방정식오일러 방정식 / 압축성 나비에-스토크스
주요 수치기법충격 포착 FVM + 리만 솔버 + TVD
대표 현상경사 충격파, 팽창부채, 파동항력

초음속 유동(supersonic flow)은 국소 음속보다 빠른 흐름, 즉 마하수 M=u/a>1M = u/a > 1인 유동이다. 정의 자체는 한 줄이지만 결과는 유체역학 전체의 문법을 바꾼다. 압력 교란이 음속으로 퍼지는데 유체가 그보다 빨리 움직이므로, 어떤 정보도 상류로 거슬러 올라가지 못한다. 아음속 유동에서 물체가 다가온다는 소식을 미리 듣고 비켜서던 공기가, 초음속에서는 물체에 얻어맞고 나서야 반응한다.

이 문서는 초음속 유동을 **“방정식이 쌍곡형이 되었을 때 물리와 수치해석이 어떻게 달라지는가”**의 관점에서 다룬다. 1차원 노즐·디퓨저의 열역학 관계식은 기체동역학이, 랭킨-위고니오 불연속 자체는 충격파가, 리만 문제의 정확해와 근사 솔버 구조는 압축성 유동과 리만 솔버가 각각 자세히 다룬다.

2. 마하 원뿔과 의존 영역[편집]

정지 매질 속을 M>1M>1로 지나가는 점원이 만든 파면들의 포락면은 원뿔이 되고, 그 반각이 마하각이다.

μ=arcsin⁡ ⁣(1M)\mu = \arcsin\!\left(\frac{1}{M}\right)

M→1M \to 1이면 μ→90∘\mu \to 90^\circ(수직 파면), MM이 커질수록 원뿔이 뾰족해진다.1 이 원뿔 바깥은 물체의 존재를 원리적으로 알 수 없는 영역이다. 수치해석 용어로 옮기면 마하 원뿔이 곧 의존 영역(domain of dependence)이고, 원뿔 하류는 영향 영역이다.

정상 초음속 유동에서 오일러 방정식이 쌍곡형이 되는 이유가 이것이다. 아음속에서는 방정식이 타원형이라 도메인 전체를 한꺼번에 풀어야 하지만, 초음속에서는 정보가 특성곡선을 따라 하류로만 흐르므로 상류에서 하류로 행진하며 순차적으로 풀 수 있다. 초음속 노즐 윤곽 설계의 고전 기법인 특성곡선법이 여기서 나왔고, 현대 CFD의 상류 차분(upwind) 철학도 뿌리가 같다. 천음속에서 아음속·초음속 영역이 뒤섞이면 방정식 형태가 공간적으로 바뀌는 혼합형 문제가 되는데, CFD가 유독 천음속을 어려워하는 근본 이유가 이거다.

3. 방향을 바꾸는 두 가지 방법[편집]

초음속 유동이 벽을 따라 방향을 트는 방법은 두 가지뿐이다. 그리고 둘은 전혀 대칭이 아니다.

안쪽으로 꺾이면(압축) 경사 충격파가 선다. 편향각 θ\theta, 충격파각 β\beta, 상류 마하수 M1M_1 사이의 관계가 그 유명한 θ-β-M 관계다.

tan⁡θ=2cot⁡β M12sin⁡2β−1M12(γ+cos⁡2β)+2\tan\theta = 2\cot\beta\,\frac{M_1^{2}\sin^{2}\beta-1}{M_1^{2}\left(\gamma+\cos 2\beta\right)+2}

주어진 M1M_1과 θ\theta에 대해 해가 두 개(약한 해·강한 해) 나오고, 자유 유동에서는 대개 약한 해가 실현된다. 더 중요한 건 θ\theta를 키우다 보면 어느 순간 실수 해가 사라진다는 것이다. 그 임계값 θmax\theta_{max}를 넘으면 충격파가 물체에 붙어 있지 못하고 앞쪽으로 떨어져 나가 이탈 충격파(bow shock)가 된다. 초음속기의 앞전이 칼날처럼 뾰족한 이유이자, 뭉툭한 재진입체 앞에 활 모양 충격파가 서는 이유다.

바깥으로 꺾이면(팽창) 프란틀-마이어 팽창부채가 생긴다. 무한히 많은 마하파가 연속적으로 퍼지는 등엔트로피 과정이며, 편향각은 프란틀-마이어 함수의 차로 계산한다.

ν(M)=γ+1γ−1arctan⁡γ−1γ+1(M2−1)−arctan⁡M2−1\nu(M) = \sqrt{\frac{\gamma+1}{\gamma-1}}\arctan\sqrt{\frac{\gamma-1}{\gamma+1}\left(M^{2}-1\right)} - \arctan\sqrt{M^{2}-1}

γ=1.4\gamma=1.4에서 νmax≈130.45∘\nu_{max} \approx 130.45^\circ인데, 이는 압력 0인 진공으로 팽창시켜도 유동이 그 이상 꺾이지 못한다는 뜻이다. 압축은 불연속으로, 팽창은 매끄럽게 — 이 비대칭은 열역학 제2법칙이 강제한다. 팽창 충격파는 엔트로피를 감소시키므로 존재할 수 없다.2

4. 노즐과 제트[편집]

라발 노즐에서 유동을 초음속으로 가속하는 원리는 면적-마하수 관계에 압축된다.

(AA∗)2=1M2[2γ+1(1+γ−12M2)]γ+1γ−1\left(\frac{A}{A^{*}}\right)^{2} = \frac{1}{M^{2}}\left[\frac{2}{\gamma+1}\left(1+\frac{\gamma-1}{2}M^{2}\right)\right]^{\frac{\gamma+1}{\gamma-1}}

같은 면적비에 아음속 해와 초음속 해가 하나씩 있고, 어느 쪽이 실현되는지는 배압이 정한다. 설계점을 벗어나면 두 가지 병이 난다. 배압이 낮으면 부족팽창(under-expanded) 제트가 노즐 밖에서 팽창부채로 마저 팽창하고, 배압이 높으면 과팽창(over-expanded) 제트가 사각 충격파로 압축된다. 둘 다 충격파와 팽창파가 제트 경계에서 번갈아 반사되며 충격 다이아몬드(shock diamond)를 만든다. 애프터버너 켠 전투기 꽁무니에 보이는 마디 무늬가 바로 그것이며, 예쁘긴 하지만 설계자 입장에선 추력을 흘리고 있다는 증거다.3

과팽창이 심하면 노즐 내부에서 유동이 박리한다. 이때 발생하는 비대칭 측력이 로켓 노즐 구조를 부수는 실제 사고 원인이 되므로, 지상 시험에서 이 영역은 항상 조심스럽게 다룬다.

5. 항력의 새로운 항목[편집]

아음속에서 항력은 마찰과 압력(형상·유도) 항력이었다. 초음속에서는 파동항력(wave drag)이 추가된다. 충격파를 지나며 발생한 엔트로피 생성이 곧 전압력 손실이고, 그 손실이 항력으로 청구된다. 점성이 0이어도 사라지지 않는 항력이라는 점이 특이하다.

천음속에서 국소 초음속 영역과 충격파가 나타나기 시작하면 항력계수가 급등하는데, 이것이 항력 발산(drag divergence)이며 “음속 장벽”의 실체다. 항력 발산 마하수를 올리려면 얇은 익형, 후퇴익, 초임계 익형, 면적 법칙(콜라병 동체)이 동원된다.

점성이 개입하면 충격파-경계층 간섭(SBLI)이 실무의 진짜 난제가 된다. 충격파의 급격한 역압력구배가 경계층을 박리시키고, 박리 거품이 다시 충격파 위치를 바꾸며, 심하면 저주파 자려 진동(버피팅)으로 발전한다. 초음속 흡입구의 언스타트, 천음속 익형의 버피팅이 전부 이 계열이다.

6. 수치해석 — 불연속을 어떻게 다룰 것인가[편집]

초음속 CFD의 역사는 곧 충격파를 계산 격자 위에 어떻게 얹을 것인가의 역사다.

  • 충격 포착 vs 충격 맞춤: 요즘 표준은 충격 포착(shock capturing)이다. 충격파를 격자 두세 칸의 급격한 구배로 자연스럽게 나오게 두는 방식으로, 유한체적법의 보존형 이산화가 필수다. 비보존형으로 풀면 랭킨-위고니오가 안 지켜져 충격파 속도가 틀린 값으로 수렴한다. 격자를 아무리 조여도 안 고쳐지는 종류의 오류라 특히 악질이다.
  • 리만 솔버와 상류 차분: 셀 경계마다 국소 리만 문제를 풀어 플럭스를 만드는 고두노프 계열(Roe, HLLC)이 주력이고, 플럭스 벡터 분할(Steger-Warming, van Leer)과 AUSM 계열(Liou)이 경쟁한다. AUSM은 대류와 압력 플럭스를 분리해 다루며, 강한 충격파에서의 견고함과 저마하수 거동에서 좋은 평판을 얻었다.
  • TVD 제한자와 차수 장벽: 고두노프의 정리에 따라 선형 도식은 2차 이상이면서 단조성을 유지할 수 없다. 그래서 MUSCL 재구성에 비선형 제한자(minmod, van Leer, superbee)를 붙여, 매끄러운 곳에서는 2차·불연속 근처에서는 1차로 스스로 낮아지게 만든다. 자세한 계열 비교는 차분 도식 참고. 제한자가 없으면 충격파 뒤에 진동이 생기고, 밀도나 압력이 음수가 되는 순간 계산은 사망한다.
  • CFL 조건과 안정성: 명시적 시간 적분의 허용 시간 간격은 최대 특성 속도 ∣u∣+a|u|+a가 정한다. 폰 노이만 안정성 해석은 선형화된 방정식에만 엄밀히 적용되므로, 강한 충격파에서는 안정 한계보다 여유 있게(CFL 0.5 안팎) 잡는 것이 현장 국룰이다.
  • carbuncle 현상: 강한 이탈 충격파를 격자선에 정렬시켜 계산하면 정체점 앞에 혹 같은 비물리적 구조가 자라난다. 접촉 불연속을 잘 보존하는 솔버(Roe, HLLC)일수록 잘 걸리는 역설적인 병이며, 원인은 충격파 내부 상태의 횡방향 수치 소산 부족으로 이해된다. 대응책은 충격 근처에서만 소산이 큰 솔버(HLL, Rusanov)로 블렌딩하기, 회전 리만 솔버, 격자를 충격파와 일부러 어긋나게 깔기 등이다.4
  • 엔트로피 고정: Roe 솔버는 음속점 근처에서 소산이 0에 가까워져 물리적으로 불가능한 팽창 충격파를 만들어낸다. 하르텐의 엔트로피 고정으로 고유값 크기에 하한을 두는 것이 표준 처방이다.

여기에 SBLI가 얽히면 난류 모델링까지 문제가 된다. 표준 난류 모델링의 RANS 모델은 충격파 유발 박리의 크기를 체계적으로 틀리게 예측하는 경향이 있어, 정밀한 예측에는 대와류 모사나 하이브리드 RANS-LES가 동원된다. 계산 비용은 당연히 폭등한다.

7. 경계 — 어디까지가 초음속인가[편집]

관례상 M>5M > 5부터를 극초음속 유동으로 따로 부른다. 숫자 자체보다 물리가 바뀌는 것이 기준이다. 충격파 뒤 온도가 수천 K에 이르면 진동 자유도 여기·해리·전리가 일어나 비열비 γ\gamma가 더 이상 상수가 아니고, 충격층이 물체에 바짝 붙으며 점성 간섭과 복사 열전달이 개입한다. 이상기체 관계식으로 돌아가던 초음속의 깔끔한 공식들이 전부 화학 반응 소스항이 붙은 연립계로 바뀌는 순간이다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 마하수는 기준 길이가 필요 없는 몇 안 되는 무차원수다. 그래서 실험실 초음속 풍동에서 잰 M=2M=2 유동장과 실물 비행체의 M=2M=2 유동장은 (레이놀즈수 효과를 뺀 비점성 부분에 한해) 서로 닮은꼴이 된다. 초음속 풍동이 작아도 되는 이유이며, 그럼에도 SBLI 예측이 어려운 것은 정확히 그 “빼놓은 부분” 때문이다. ↩

  2. 예외처럼 보이는 사례가 없지는 않다. 임계점 근처의 조밀 기체(BZT 유체)처럼 기본 기체동역학 미분량이 음수인 매질에서는 팽창 충격파가 열역학적으로 허용된다. 유기 랭킨 사이클 터빈 해석에서 실제로 논의되는 주제이지만, 공기로는 안 된다. ↩

  3. 충격 다이아몬드가 보인다는 건 노즐이 설계점에서 벗어나 있다는 뜻이다. 완벽하게 설계된 고도에서 완벽하게 팽창한 로켓 노즐의 배기는 무늬 없이 밋밋하다. 즉 사진이 멋있을수록 성능은 손해다. 우주 덕후와 추진 엔지니어의 미적 취향이 정확히 반대인 이유. ↩

  4. carbuncle은 1988년 Peery와 Imlay가 보고한 이래 40년 가까이 완전한 해결책이 없다. “정확한 솔버일수록 잘 걸린다”는 성질 때문에 수치 소산을 늘리면 사라지지만, 그러면 애초에 정확한 솔버를 쓴 이유가 사라진다. 정확도와 견고성의 교환이 이보다 노골적인 사례도 드물다. ↩