랭킨-위고니오 조건

수치해석 전산유체역학 마지막 수정: 2026-08-26 04:22:51

1. 개요[편집]

랭킨-위고니오 조건 Rankine–Hugoniot condition
형태	s [U] = [F(U)]
유래	랭킨(1870) · 위고니오(1887~89)
성격	약해가 불연속면에서 만족해야 하는 대수 조건
버거스	s = (u_L + u_R) / 2
한계	해를 유일하게 정하지 못함 — 엔트로피 조건 필요
수치적 의미	보존형 도식이 자동으로 만족, 비보존형은 위반

방정식이 미분을 요구하는데 해가 미분 불가능하다면, 미분을 요구하지 않는 형태로 다시 물어봐야 한다.

랭킨-위고니오 조건(Rankine–Hugoniot condition)은 보존형 편미분방정식의 약해에서 불연속면 양쪽의 상태 $U_L,U_R$ 와 그 면의 전파 속도 $s$ 를 잇는 대수 조건으로,

s\,(U_R - U_L) = F(U_R) - F(U_L) \qquad\text{줄여서}\qquad s\,[U] = [F(U)]

의 꼴을 갖는다. 여기서 $[\,\cdot\,]$ 는 면을 가로지르는 점프다.

쌍곡형 보존법칙 $U_t+F(U)_x=0$ 은 매끄러운 초기조건에서도 유한 시간에 불연속을 만든다(버거스 방정식의 특성선 교차가 표준 예). 그 순간부터 고전해는 존재하지 않고, 미분 대신 적분형 보존을 만족하는 약해로 넘어가야 한다. 랭킨-위고니오 조건은 그 약해가 불연속면에서 반드시 지켜야 하는 유일한 국소 조건이며, 뒤집어 말하면 “충격이 얼마나 빨리 가는가”를 결정하는 식이다.

이 문서는 조건 자체와 그 수학적·수치적 함의를 다룬다. 기체동역학적 충격파 현상(경사 충격파, 이탈, 총압 손실)은 충격파가, 셀 경계 불연속을 실제로 푸는 알고리즘은 리만 솔버와 고두노프 도식이 담당한다.¹

2. 약해에서 유도하기[편집]

시험함수 $\phi\in C_c^\infty$ 를 곱해 부분적분하면, 미분을 전부 $\phi$ 로 넘긴 약형식이 정의된다.

\int_0^\infty\!\!\int_{-\infty}^{\infty}\left(U\phi_t + F(U)\phi_x\right)dx\,dt + \int_{-\infty}^{\infty}U_0(x)\phi(x,0)\,dx = 0

이 식은 $U$ 의 미분을 요구하지 않으므로 계단 함수도 후보가 된다. 이제 $U$ 가 곡선 $x=X(t)$ 를 사이에 두고 양쪽에서는 매끈하고 그 위에서만 뛴다고 하자. 곡선을 감싸는 작은 영역을 좌·우로 나눠 각각 발산 정리를 적용하면, 양쪽의 매끈한 부분은 원 방정식을 만족하니 사라지고 곡선 위 선적분만 남는다.

\int \left\{ \dot X\,[U] - [F(U)] \right\}\phi\,d\ell = 0

$\phi$ 가 임의이므로 중괄호가 0이어야 하고, $s=\dot X$ 로 쓰면 $s[U]=[F(U)]$ 다.

물리적으로는 같은 말을 검사체적으로 할 수 있다. 불연속면과 함께 움직이는 좌표계에서는 면이 정지해 있으므로, 면을 통과해 들어온 플럭스와 나간 플럭스가 같아야 한다. 상대속도로 쓰면 $F(U)-sU$ 가 면을 가로질러 연속이라는 뜻이고, 이것이 위 식의 다른 표현이다.

여기서 반드시 짚어야 할 함정이 하나 있다. 어떤 양을 보존시킬지는 물리가 정하는 것이지 미분 규칙이 정하는 게 아니다. 무점성 버거스는 $u_t+(u^2/2)_x=0$ 로도 쓸 수 있지만, 매끈한 해에서는 $(u^2)_t+(2u^3/3)_x=0$ 도 동시에 성립한다. 그런데 두 형태의 랭킨-위고니오 속도는

s = \frac{u_L+u_R}{2} \qquad\text{vs}\qquad s' = \frac{2}{3}\cdot\frac{u_L^2+u_Lu_R+u_R^2}{u_L+u_R}

로 서로 다르다( $u_L=2,\,u_R=0$ 이면 $1$ 대 $4/3$ ). 둘 다 자기 나름의 약해지만 물리적으로 옳은 것은 하나뿐이고, 어느 쪽인지는 방정식이 어디서 왔는지(무엇이 진짜로 보존되는 양인지)가 결정한다. “미분하면 같은데요”가 통하지 않는 지점이며, 코드가 비보존형으로 짜였을 때 충격 속도가 틀리는 현상의 뿌리도 정확히 여기다.²

3. 스칼라 — 그냥 코드의 기울기[편집]

스칼라 보존법칙 $u_t+f(u)_x=0$ 에서는 조건이 한 줄로 끝난다.

s = \frac{f(u_R)-f(u_L)}{u_R-u_L}

즉 플럭스 함수 그래프 위 두 점을 잇는 코드(chord)의 기울기다. 특성 속도 $f'(u)$ 가 접선의 기울기라는 것과 나란히 놓으면 그림 하나로 정리된다.

버거스 방정식 $f=u^2/2$ → $s=(u_L+u_R)/2$ . 좌우 값의 산술평균이라는 유난히 예쁜 결과가 나오는 이유는 플럭스가 2차식이라 코드 기울기가 중점의 접선 기울기와 같기 때문이다.
교통류(LWR) 모형 $f=\rho v(\rho)$ → $s=(q_R-q_L)/(\rho_R-\rho_L)$ . 정체 충격의 속도가 음수로 나오는 것 — 즉 정체 꼬리가 차량 진행 방향과 반대로 거슬러 올라가는 것 — 이 이 식의 직접적 귀결이다.
얕은 물 방정식 은 스칼라가 아니지만, 질량·운동량 두 점프 조건에서 속도를 소거하면 개수로 유동의 도수 공액수심 관계 $y_2/y_1=\tfrac12(\sqrt{1+8\mathrm{Fr}_1^2}-1)$ 이 그대로 나온다. 도수는 물리적으로 충격파와 같은 대상이라는 말이 수식 수준에서 성립한다는 뜻.

4. 오일러 방정식의 세 점프 조건[편집]

오일러 방정식 (비점성 유동)에서 $U=(\rho,\ \rho u,\ E)^{\mathsf T}$ , $F=(\rho u,\ \rho u^2+p,\ (E+p)u)^{\mathsf T}$ 이므로 조건은 세 개다. 충격과 함께 움직이는 좌표( $v=u-s$ )에서 쓰면 익숙한 꼴이 된다.

\rho_1 v_1 = \rho_2 v_2 \equiv \dot m,\qquad p_1+\rho_1 v_1^2 = p_2+\rho_2 v_2^2,\qquad h_1+\tfrac12 v_1^2 = h_2+\tfrac12 v_2^2

질량·운동량 두 식에서 속도를 소거하면 비체적 $\tau=1/\rho$ 평면에서 레일리 선(Rayleigh line)이 나온다.

\dot m^2 = \frac{p_2-p_1}{\tau_1-\tau_2}

$(\tau,p)$ 평면에서 상태 1과 2를 잇는 직선의 기울기가 $-\dot m^2$ 라는 뜻이라, 기울기는 반드시 음수 — 압력이 오르면 비체적은 줄어야 한다. 여기에 에너지식을 넣어 속도를 완전히 소거하면 상태량만의 관계인 위고니오 곡선(Hugoniot curve, 충격 단열선)이 남는다.

e_2-e_1=\tfrac12(p_1+p_2)\left(\tau_1-\tau_2\right)

주어진 상태 1에서 출발하는 위고니오 곡선과 레일리 선의 교점이 상태 2다. 곡선은 상태방정식이, 직선은 충격 강도가 결정한다.

4.1. 위고니오 곡선은 단열선이 아니다[편집]

가장 흔한 오해가 이것이다. 위고니오 곡선은 등엔트로피 곡선(isentrope)과 다르다.

두 곡선은 상태 1에서 만나고, 2차까지 접한다. 그래서 약한 충격은 등엔트로피 압축과 사실상 구별되지 않는다.
3차부터 갈라지고, 그 차이가 곧 충격이 만드는 엔트로피다. 약한 충격에서 $\Delta s \propto (p_2-p_1)^3$ 인 것이 이 접촉 차수의 결과다.
등엔트로피 압축은 무한히 압축할 수 있지만, 위고니오는 압축비에 상한이 있다. 완전기체에서 $\tau_2/\tau_1\to(\gamma-1)/(\gamma+1)$ 이 벽이고, $p_2\to\infty$ 여도 그 벽을 못 넘는다.

완전기체·일정 $\gamma$ 라면 상류 마하수 하나로 전부 닫힌다. 밀도비는

\frac{\rho_2}{\rho_1}=\frac{(\gamma+1)M_1^2}{(\gamma-1)M_1^2+2} \;\xrightarrow[M_1\to\infty]{}\; \frac{\gamma+1}{\gamma-1}

이고 공기( $\gamma=1.4$ )에서는 겨우 6배에서 포화한다. 압력비와 온도비는 $M_1^2$ 에 비례해 무한히 커지는데 밀도는 6배에서 멈추니, 남는 에너지가 전부 내부에너지로 간다. 재진입체가 녹는 이유이자, 극초음속 코드가 실기체 효과(해리·전리)를 넣지 않으면 온도를 수천 K 과대평가하는 이유다. 자세한 유동 현상은 충격파·압축성 유동에 있다.

발열이 있는 경우로 확장하면 위고니오 곡선이 위로 들리는데(반응 위고니오), 레일리 선이 그 곡선에 접하는 특별한 점이 채프먼-주게 조건이고 그것이 폭굉 속도를 결정한다. 자착화·노킹에서 압력파가 자기 유지 폭굉으로 전이하는지 여부가 이 접점 문제로 판정된다.

5. 위고니오 궤적과 리만 문제[편집]

$p$ 개 방정식계에서 미지수는 $U_R$ 의 $p$ 성분과 $s$ 로 $p+1$ 개인데 조건은 $p$ 개다. 따라서 주어진 $U_L$ 에 대해 랭킨-위고니오를 만족하는 상태들은 1-매개변수 곡선족을 이루며, 이를 위고니오 궤적(Hugoniot locus)이라 부른다. 락스의 정리는 각 특성장마다 이런 곡선이 하나씩, 총 $p$ 개가 $U_L$ 을 지나며 다음을 만족한다고 말한다.

$U_R\to U_L$ 극한에서 곡선은 $k$ 번째 우고유벡터 $r_k$ 방향으로 접한다.
같은 극한에서 $s\to\lambda_k(U_L)$ — 약한 충격은 음파 속도로 간다.
같은 점을 지나는 희박파 곡선(적분곡선)과 위고니오 궤적은 2차까지 일치한다.

마지막 항목이 실무적으로 중요하다. 로(Roe) 같은 근사 리만 솔버가 충격과 희박파를 같은 선형화로 뭉개 놓고도 잘 도는 이유가 “약한 파에서는 두 곡선이 어차피 구별되지 않기 때문”이다. 강한 파에서 그 근사가 무너지면 카번클(carbuncle) 현상이나 음의 압력·음의 내부에너지 같은 사고가 난다.

특성장이 선형 퇴화(linearly degenerate)인 경우는 예외적으로 깔끔하다. 오일러의 2번 장이 그렇고, 그 불연속이 접촉 불연속이다. $s=u_1=u_2$ , $[p]=0$ , $[u]=0$ 이고 밀도와 온도만 뛴다. 특성선이 불연속면과 평행하므로 충격처럼 “특성선을 삼키는” 일이 없고, 그래서 수치적으로는 반대 문제 — 아무것도 이 면을 붙잡아 주지 않아 수치 확산으로 계속 번진다 — 가 생긴다. 소드 문제에서 충격은 세 셀 안에 서는데 접촉면만 시간이 갈수록 두꺼워지는 광경이 그것이다.

리만 문제의 해는 결국 이 곡선들을 이어 붙이는 작업이다. $U_L$ 에서 1번 곡선을 타고 가다가 중간 상태에서 2번, 3번으로 갈아타 $U_R$ 에 도착하는 경로를 찾는 비선형 대수 문제이고, 오일러에서는 압력에 대한 1변수 방정식으로 축약되어 뉴턴-랩슨법으로 풀린다.

6. 랭킨-위고니오만으로는 부족하다[편집]

조건을 만족하는 약해는 유일하지 않다. 버거스에서 $u_L<u_R$ 인 초기 불연속을 생각하면, 그 자리에 그대로 서 있는 불연속도 랭킨-위고니오를 완벽히 만족한다. 하지만 물리적으로 옳은 답은 희박팬이다. 그래서 부등식 하나가 더 필요하고, 이것이 엔트로피 조건이다.

조건	진술	적용 범위
락스(Lax)	$\lambda_k(U_L)>s>\lambda_k(U_R)$ , 그리고 $\lambda_{k-1}(U_L)<s<\lambda_{k+1}(U_R)$	진성 비선형장, 약한 충격
올레이닉 E	코드 기울기가 사이의 모든 $u$ 에서 부등호를 유지	스칼라, 플럭스가 비볼록이어도 성립
류(Liu)	위고니오 궤적을 따라 $s$ 가 감소	계, 비볼록 플럭스
엔트로피 쌍	볼록 $\eta$ 와 짝 플럭스 $q$ 에 대해 $[q] \le s\,[\eta]$	가장 일반적, 약형식 부등식

기하학적 직관은 하나로 통한다. 특성선이 충격으로 들어가야 한다. 충격은 정보를 삼키는 구조여야지 뱉는 구조여서는 안 되고, 뱉는 쪽은 희박파가 담당한다. 이 조건을 명시적으로 넣지 않으면 고두노프 도식의 소박한 변형이나 로 솔버가 음속점에서 정지한 팽창 충격을 만들어 낸다. 엔트로피 픽스란 결국 이 부등식을 손으로 다시 집어넣는 작업이다.

점성 소실 극한으로 정의하는 관점도 같은 답을 준다. $u_t+f(u)_x=\varepsilon u_{xx}$ 의 $\varepsilon\to0^+$ 극한이 유일한 엔트로피 해이고, 버거스에서는 콜-호프 변환이 이 사실을 계산 가능한 적분으로 보여 준다. 계에서는 “충격 내부에 매끄러운 점성 프로파일이 존재하는가”를 묻는 점성 프로파일 판정이 되고, 이것이 락스 조건보다 강한 선별 기준을 주는 경우가 있다(비고전 충격, 상전이 문제).

7. 수치해석에서 왜 중요한가[편집]

락스-벤드로프 정리(1960)가 이 문서 전체의 실무적 결론이다. 보존형 도식

U_j^{n+1}=U_j^n-\frac{\Delta t}{\Delta x}\left(F_{j+1/2}-F_{j-1/2}\right)

이 일관성(consistency)을 갖고, 그 수치해가 유계로 수렴한다면, 극한은 반드시 약해이며 따라서 충격은 랭킨-위고니오 속도로 간다. 도식이 충격 속도 공식을 알 필요가 없다 — 이산 수준에서 플럭스가 상쇄되며 텔레스코핑되기만 하면 조건은 저절로 따라온다. 유한체적법이 충격이 있는 문제의 표준이 된 이유가 이 한 정리다.

반대로 비보존형 이산화는 격자를 아무리 줄여도 틀린 속도로 수렴한다. 잔차는 예쁘게 떨어지고 해는 매끄럽게 보이는데 충격 위치만 계속 틀린, 가장 발견하기 어려운 종류의 버그다.

근사 리만 솔버 설계에도 조건이 직접 들어간다. 로 평균 행렬 $\hat A(U_L,U_R)$ 가 만족해야 하는 세 요건 중 하나가

\hat A(U_L,U_R)\,(U_R-U_L)=F(U_R)-F(U_L)

인데, 이것은 이산화된 랭킨-위고니오 조건 그 자체다(U-property). 이 성질 덕분에 로 솔버는 단일 충격을 정확히 한 셀 안에 해상한다. 오일러에서 이 요건을 푸는 해가 그 유명한 $\sqrt{\rho}$ 가중 평균이고, 로가 논문에서 “우연히 찾았다”고 쓴 그 평균이다.³

몇 가지 남는 주제.

비보존형 계. 다상유동의 7-방정식 모형이나 층별 얕은 물처럼 $U_t+A(U)U_x=0$ 로만 쓰이는 계에서는 $A(U)U_x$ 가 불연속에서 초함수 곱이 되어 랭킨-위고니오 자체가 정의되지 않는다. 달 마소-르플로슈-뮈라(DLM) 이론은 상태공간에 경로족을 지정해 점프를 정의하는데, 경로를 바꾸면 충격 속도가 바뀐다. 즉 답이 모형화 선택에 의존한다.
충격 적합 vs 충격 포착. 랭킨-위고니오를 경계조건으로 명시적으로 걸고 불연속면을 격자면으로 추적하는 것이 충격 적합(shock fitting)이다. 정확도는 압도적이지만 충격이 갈라지거나 생겨날 때 위상 관리가 지옥이라, 현대 코드는 WENO 도식·MUSCL 기반 포착으로 통일됐다.
검증 시 확인할 것. 충격 속도가 맞는지는 검증 및 확인에서 가장 값싸게 잡히는 항목이다. 소드 문제를 돌린 뒤 충격 위치를 해석해와 비교해서 어긋나면, 대개 플럭스가 보존형이 아니거나 소스항 처리가 보존을 깨고 있다.

8. 관련 문서[편집]

이름의 지분은 좀 복잡하다. 스토크스가 1848년에 질량·운동량 점프 조건을 이미 적었지만 에너지 보존 문제로 비판받고 스스로 철회했고, 리만은 1860년 논문에서 충격을 가로질러 엔트로피가 보존된다고 가정해 틀린 점프 조건을 얻었다. 에너지 보존을 제대로 넣어 지금의 형태를 세운 것이 랭킨(1870)과 위고니오(1887~89)다. 그러면서도 초기 불연속을 두 상태로 놓고 푸는 그 문제 설정에는 여전히 리만의 이름이 붙어 있다는 점이 이 바닥 명명법의 정취를 잘 보여 준다. 틀린 답을 냈어도 옳은 질문을 먼저 던지면 이름이 남는다. ↩
이 애매함은 이론적 호기심이 아니라 실제 코드 설계 문제다. 압축성 코드에서 총에너지 대신 엔트로피나 압력을 보존 변수로 삼는 정식화(이른바 압력 기반·엔탈피 기반 형태)가 종종 유혹적인데, 매끈한 해에서는 동등해도 충격을 지나는 순간 답이 달라진다. 저마하 유동에서는 잘 돌던 코드가 마하 2에서 충격 위치를 몇 셀씩 틀리기 시작하면, 격자보다 먼저 “무엇을 보존시키고 있는지”를 확인해야 한다. ↩
로(P. L. Roe)의 1981년 논문은 이 평균을 “매개변수 벡터 $(\sqrt\rho,\sqrt\rho u,\sqrt\rho H)$ 를 도입하면 플럭스 차이가 그 벡터에 대해 2차형식이 되어 자동으로 풀린다”는 방식으로 얻었다. 나중에 임의의 상태방정식으로 확장하려는 시도가 수십 편 나왔는데, 대부분 “적분 평균을 수치적으로 구해야 한다”는 결론에 도달했다. 완전기체에서만 공짜였던 것이다. ↩