엔트로피 조건

수치해석 전산유체역학 마지막 수정: 2026-08-27 04:12:41

1. 개요[편집]

엔트로피 조건 Entropy Condition
푸는 문제	약해의 비유일성
역할	물리적으로 옳은 약해 하나를 선별
대표 형태	락스 · 올레이닉 E · 류 · 엔트로피 쌍
유일성 정리	크루즈코프 (1970, 스칼라 다차원)
물리적 근거	점성 소실 극한 · 열역학 제2법칙
수치적 대응물	엔트로피 안정 도식 · 엔트로피 픽스
어겼을 때	팽창 충격(expansion shock)

랭킨-위고니오는 “이 점프가 보존을 깨지 않는가”만 묻는다. 그게 실제로 일어나는 일인지는 묻지 않는다.

엔트로피 조건(entropy condition)은 비선형 쌍곡형 보존법칙의 약해가 무수히 많을 때, 그중 물리적으로 실현되는 해 하나를 골라내기 위해 추가로 부과하는 부등식이다. 보존법칙 $u_t+f(u)_x=0$ 은 해가 불연속을 품기 때문에 미분형으로는 의미가 없고 적분형(약해)으로 정의해야 하는데, 약해의 정의는 유일성을 주지 않는다. 방정식이 답을 하나로 정해 주지 않는다는 뜻이다.

여기서 흔한 오해 하나를 먼저 걷어내자. 랭킨-위고니오 조건은 유일성 조건이 아니다. 그것은 “불연속을 가로질러 플럭스가 상쇄되는가”를 확인하는 필요조건일 뿐이고, 이를 완벽하게 충족하면서도 자연에서는 절대 관측되지 않는 해가 실제로 존재한다. 엔트로피 조건은 그 틈을 메우는 별도의 부등식이며, 미분방정식이 아니라 부등식이라는 점이 본질이다. 등식은 방향을 모르지만 부등식은 시간의 화살을 안다.

2. 비유일성이 얼마나 심각한가[편집]

버거스 방정식 $u_t+(u^2/2)_x=0$ 에 다음 초기조건을 준다.

u(x,0)=\begin{cases}0,& x<0\\ 1,& x>0\end{cases}

두 개의 후보가 있다.

팽창 충격. 불연속을 그대로 유지한 채 속도 $s=(u_L+u_R)/2=1/2$ 로 흘려보낸다. 랭킨-위고니오를 정확히 만족하므로 이것도 약해다. 다만 압축이 아니라 팽창하는 점프라 자연에는 없다.
희박팬. $u=x/t$ 로 부채꼴로 벌어지는 연속해. 이것도 약해다.

둘 다 적분형 보존법칙을 만족하고, 둘 다 초기조건에 부합한다. 방정식만으로는 우열을 가릴 수 없다. 더 나쁜 소식은 이게 두 개짜리 선택 문제가 아니라는 것 — 충격을 여러 조각으로 쪼개 끼워 넣으면 연속체 무한 개의 약해를 만들 수 있다. 그래서 약해 이론은 태생적으로 선별 원리를 요구한다.

물리적 정답은 희박팬이다. 이유는 원래 방정식이 어디서 왔는지를 보면 나온다. 보존법칙은 대개 점성이 있는 방정식의 극한이다.

u_t+f(u)_x=\varepsilon u_{xx},\qquad \varepsilon\to 0^+

$\varepsilon>0$ 이면 해는 매끄럽고 유일하다. 그 유일한 해의 극한을 엔트로피 해로 정의한다. 위 예에서 점성 해는 부채꼴로 벌어지고, 팽창 충격은 어떤 $\varepsilon>0$ 의 극한으로도 나오지 않는다. 이 관점을 완전비선형 1·2계 PDE로 확장한 것이 점성해 이론이고, 이름의 유래도 같다.

3. 락스 조건 — 특성선은 충격으로 들어가야 한다[편집]

가장 직관적인 판정은 락스(Lax, 1957)의 것이다. 진성 비선형(genuinely nonlinear) $k$ -번째 특성장의 충격에 대해

\lambda_k(U_L)>s>\lambda_k(U_R),\qquad \lambda_{k-1}(U_L)<s<\lambda_{k+1}(U_R)

를 요구한다. 첫 줄이 핵심이다. 좌측 특성 속도가 충격보다 빠르고 우측 특성 속도가 충격보다 느리다는 것은, 양쪽에서 특성선이 충격면으로 수렴해 들어간다는 뜻이다. 충격은 정보를 삼키는 구조이지 뱉는 구조가 아니다. 반대로 특성선이 충격에서 뻗어 나오는 배치는 “정보가 무에서 창조되는” 상황이라 물리적으로 말이 안 되고, 그 자리는 희박파가 맡아야 한다.

두 번째 줄은 계(system)에서 필요한 조건으로, 충격을 가로지르는 특성선의 개수를 세면 점프 조건의 개수와 미지수 개수가 맞아떨어진다는 요구다. $n$ -방정식 계에서 $k$ -충격에 들어오는 특성선은 $n+1$ 개여야 하고, 이것이 랭킨-위고니오 $n$ 개 + 속도 $s$ 하나를 결정하는 정보량과 정확히 맞는다.¹

버거스처럼 볼록 플럭스인 스칼라에서는 $\lambda=f'(u)=u$ 이므로 락스 조건은 그냥 $u_L>u_R$ 이다. 왼쪽이 크면 충격, 오른쪽이 크면 희박파. 위 예제는 $u_L=0<u_R=1$ 이므로 충격은 탈락이다.

4. 올레이닉 조건 — 볼록 껍질이 답을 그려 준다[편집]

플럭스가 볼록하지 않으면 락스 조건만으로는 부족하다. 예컨대 $f(u)=u^3/3$ 이나 부호가 바뀌는 곡률을 가진 플럭스에서는, 락스 조건을 만족하는 충격이 여럿 나오거나 충격과 희박파가 붙어 있는 복합파(composite wave)가 정답이 된다. 이때 쓰는 것이 올레이닉(Oleinik)의 E-조건이다.

\frac{f(u)-f(u_L)}{u-u_L}\ \ge\ s\ \ge\ \frac{f(u)-f(u_R)}{u-u_R} \quad\text{(}u_L,\,u_R\text{ 사이의 모든 } u\text{ 에 대해)}

말로 풀면 이렇다. 플럭스 곡선 위의 두 점 $(u_L,f(u_L))$ , $(u_R,f(u_R))$ 을 잇는 현(chord)의 기울기가 곧 충격 속도 $s$ 인데, 그 현이 사이 구간에서 곡선의 어느 쪽에 있어야 하는지를 지정하는 것이다.

기하학적 처방은 더 깔끔하다. 스칼라 리만 문제의 엔트로피 해는 플럭스의 껍질(hull) 로 완전히 결정된다.

$u_L>u_R$ 이면 구간 $[u_R,u_L]$ 위에서 $f$ 의 오목 껍질(upper concave envelope)을 그린다.
$u_L<u_R$ 이면 구간 $[u_L,u_R]$ 위에서 $f$ 의 볼록 껍질(lower convex envelope)을 그린다.

그 껍질에서 직선 구간은 충격(기울기 = 속도), 곡선 구간은 희박파에 대응한다. 볼록 플럭스라면 첫 경우의 오목 껍질이 통째로 현이 되어 단일 충격, 둘째 경우의 볼록 껍질이 $f$ 자신이 되어 순수 희박팬 — 락스 조건과 정확히 같은 답이다. 비볼록 플럭스에서만 껍질이 직선과 곡선을 번갈아 붙여 복합파를 자동으로 생성한다. 자를 대고 그으면 답이 나온다는 점에서, 비볼록 보존법칙을 다룰 때 이만큼 실용적인 도구가 없다.²

계로 넘어가면 대응물이 류(Liu)의 조건이다. 위고니오 궤적을 따라 충격 속도가 감소해야 한다는 요구이고, 비볼록 플럭스나 상전이가 얽힌 계에서 락스보다 강한 선별력을 준다.

4.1. 껍질 구성이 실제로 쓰이는 곳 — 벅클리-레버렛[편집]

교과서 예제로 끝나는 이야기가 아니다. 다공질 매체의 2상 유동(물이 기름을 밀어내는 수공법)은 다르시 법칙에서 출발해 포화도 $S$ 에 대한 스칼라 보존법칙

S_t+\left(\frac{S^2}{S^2+M(1-S)^2}\right)_x=0

으로 정리된다(벅클리-레버렛 방정식, $M$ 은 점도비). 이 분율 유동 함수는 S자 곡선이라 볼록도 오목도 아니다. 락스 조건만 믿고 $S=1$ 에서 $S=0$ 으로 단일 충격을 놓으면 물리적으로 틀린 답이 나온다.

올레이닉 껍질 구성이 곧바로 정답을 준다. $S_R=0$ 점에서 곡선에 접선을 그어 접점 $S^*$ 를 찾으면, 해는 $S=1$ 부터 $S^*$ 까지의 희박파에 이어 $S^*$ 에서 0으로 떨어지는 충격으로 이루어진 복합파가 된다. 석유공학에서 웰지(Welge) 접선 구성이라 부르는 이 작도가 실은 볼록 껍질의 직선 구간을 찾는 작업과 완전히 같은 것이다. 물이 도달하는 선단 시각(breakthrough time)이 이 접점 하나로 결정되므로, 엔트로피 조건이 곧바로 돈이 되는 드문 사례이기도 하다.

5. 크루즈코프 엔트로피 쌍과 유일성[편집]

앞의 조건들은 “불연속 하나”에 대한 국소 판정이다. 이를 해 전체에 대한 약형식 부등식으로 바꾸면 훨씬 강력한 도구가 된다. 엔트로피 쌍(entropy pair) $(\eta,q)$ 는 $\eta$ 가 볼록하고 $q'=\eta'f'$ 를 만족하는 함수 쌍이며, 엔트로피 해는 모든 그런 쌍에 대해

\partial_t\eta(u)+\partial_x q(u)\ \le\ 0

를 분포 의미로 만족한다고 정의한다. 매끄러운 구간에서는 등식이 되고(엔트로피 보존), 충격에서만 진짜 부등식이 되어 엔트로피가 소산된다. 불연속을 가로질러 적분하면 $[q]\le s\,[\eta]$ 라는 점프 형태가 나오는데, 이것이 락스·올레이닉 조건을 포함하는 가장 일반적인 진술이다.

부호가 헷갈리기 쉬운데, 수학적 엔트로피 $\eta$ 는 볼록이고 그 총량이 감소한다. 물리 엔트로피는 오목하고 증가하므로, 압축성 유동에서 쓰는 수학적 엔트로피는 $\eta=-\rho s$ 처럼 부호를 뒤집은 것이다. 하르텐(Harten, 1983)이 오일러 방정식 (비점성 유동)에 대해 $\eta=-\rho h(s)$ 꼴이 언제 볼록인지를 정리했고, 이 볼록성이 뒤에 나올 수치 도식 설계의 전제가 된다.

크루즈코프(Kružkov, 1970)는 스칼라 다차원 보존법칙에서 엔트로피 쌍을 특정 족으로 고정해 유일성을 증명했다. 모든 상수 $k$ 에 대해

\eta_k(u)=|u-k|,\qquad q_k(u)=\operatorname{sgn}(u-k)\,\bigl(f(u)-f(k)\bigr)

이 족만 쓰면 충분하다는 것이 요점이다. 이 부등식을 모든 $k$ 에 대해 만족하는 $L^\infty$ 해는 유일하고, 두 해 사이의 $L^1$ 거리가 시간에 대해 비증가한다( $L^1$ 수축). 증명 기법인 변수 이중화(doubling of variables)는 지금도 보존법칙 이론의 표준 도구다.³

계에서는 사정이 훨씬 나쁘다. 1차원 BV 소데이터에서는 브레상 학파의 결과로 well-posedness가 확립됐지만, 다차원에서는 엔트로피 조건을 만족하는 해가 여전히 무수히 많을 수 있다는 것이 볼록적분(convex integration) 기법으로 구성되어 있다. 즉 “엔트로피 조건을 넣으면 답이 하나”라는 문장은 스칼라에서는 정리이고 다차원 계에서는 희망사항이다.

6. 수치적 대응물 (1) — 엔트로피 픽스[편집]

여기까지가 이론이고, CFD 엔지니어가 이 문서를 검색해서 들어오는 이유는 대개 다음 한 줄 때문이다. 로(Roe) 솔버가 비물리적인 팽창 충격을 뱉었다.

로 솔버는 야코비안을 상수 행렬 $\hat A(U_L,U_R)$ 로 선형화하고, 그 행렬이 $\hat A\,(U_R-U_L)=F(U_R)-F(U_L)$ 를 만족하도록 설계한다(U-property). 이 성질 덕분에 단일 충격을 한 셀 안에 정확히 잡는 대신, 치명적인 부작용이 생긴다. 수치 플럭스의 소산항이 $|\hat\lambda_k|$ 에 비례하는데, 음속점(sonic point, $\lambda_k$ 가 좌우에서 부호를 바꾸는 지점)에서 $\hat\lambda_k=0$ 이 되면 그 특성장의 소산이 완전히 사라진다. 소산이 0이면 도식은 랭킨-위고니오만 보고 판단하게 되고, 정지한 팽창 불연속은 랭킨-위고니오를 만족하므로 그대로 유지된다. 격자를 아무리 줄여도 사라지지 않는다 — 버그가 아니라 도식이 인정하는 정상 해다.

처방이 엔트로피 픽스(entropy fix)다. 본질은 ” $|\hat\lambda|$ 가 0 근처로 내려가지 못하게 바닥을 깔아 인공 소산을 되살리는 것”이다.

방식	처리	특징
하르텐 $\varepsilon$ -픽스	$\delta$ 이하에서 특성 속도의 절댓값을 포물선으로 매끄럽게 교체	구현 3줄, $\delta$ 를 손으로 고름
하르텐-하이먼	좌우 상태의 특성 속도 차이로 $\delta$ 를 자동 산정	문제 의존성이 적음
HLL 계열 혼합	음속점 부근만 HLL/러스노프 플럭스로 대체	안정하지만 접촉면이 뭉개짐
도식 교체	애초에 엔퀴스트-오셔, HLLC, 러스노프를 씀	픽스 자체가 불필요

주의할 점은 픽스가 공짜가 아니라는 것이다. $\delta$ 를 크게 잡으면 팽창 충격은 확실히 사라지지만 접촉 불연속과 경계층 해상도가 같이 죽는다. 극초음속 유동 코드에서 카벙클 현상을 잡겠다고 픽스를 과하게 걸었다가 벽면 열유속을 통째로 틀리는 사고가 흔한 이유다. 엔트로피 픽스는 필요한 곳에만 최소한으로 거는 게 국룰이다.

반대로 고두노프 도식처럼 정확 리만 솔버를 쓰면 엔트로피 조건이 리만 해 안에 이미 들어 있으므로 픽스가 필요 없다. 러스노프·엔퀴스트-오셔·HLL 같은 소산이 넉넉한 리만 솔버도 마찬가지다. 정확도와 엔트로피 안전성 사이의 거래가 근사 리만 솔버 설계의 절반이다.

7. 수치적 대응물 (2) — 수렴 정리와 엔트로피 안정 도식[편집]

락스-벤드로프 정리는 보존형 도식의 수렴 극한이 약해임을 보장한다. 그런데 지금까지 봤듯이 약해는 여러 개다. 그래서 정리 하나가 더 필요하다. 도식이 이산 엔트로피 부등식

\eta_j^{n+1}\le \eta_j^n-\frac{\Delta t}{\Delta x}\left(Q_{j+1/2}-Q_{j-1/2}\right)

를 정합적인 수치 엔트로피 플럭스 $Q$ 에 대해 만족하면, 그 극한은 엔트로피 해다. 즉 락스-벤드로프(→약해) + 이산 엔트로피 부등식(→선별) 두 장이 있어야 “이 코드가 옳은 답으로 수렴한다”가 완성된다.

이 부등식이 공짜로 따라오는 대표 부류가 단조 도식(monotone scheme)이다. 크랜들-마이다(1980)가 단조 도식의 극한이 엔트로피 해임을 증명했다. 문제는 하르텐-하이먼-락스(1976)가 이미 단조 도식은 최대 1차 정확도임을 보였다는 것 — 안전하지만 뭉개진다. 여기서 벗어나려고 전변분 감소·MUSCL·WENO 도식으로 올라가지만, 고차 도식에서 이산 엔트로피 부등식은 자동으로 따라오지 않는다.

이를 정면으로 설계하는 것이 엔트로피 안정 도식(entropy stable scheme)이다. 타드모어(Tadmor, 1987)의 골격은 두 단계다.

엔트로피 보존 플럭스. 엔트로피 변수 $v=\partial\eta/\partial U$ 로 옮긴 뒤, 수치 플럭스 $F^*$ 가 $(v_{j+1}-v_j)^{\mathsf T}F^*_{j+1/2}=\psi_{j+1}-\psi_j$ 를 만족하도록 잡는다( $\psi=v^{\mathsf T}f-q$ 는 엔트로피 포텐셜). 이러면 이산 엔트로피가 정확히 보존된다.
소산 추가. 여기에 엔트로피 변수 점프 방향으로 음의 준정부호 소산을 얹는다. 그러면 부등식이 확실히 옳은 방향으로 기운다.

엔트로피 변수가 계를 대칭화한다는 고두노프-모크-하르텐의 구조가 이 설계의 수학적 기반이다. 이 계열은 고차 정확도를 유지하면서도 엔트로피 부등식을 이산 수준에서 증명할 수 있어서, 2010년대 이후 불연속 갤러킨법이나 고차 유한차분에서 표준적인 안정화 전략으로 자리 잡았다. 물론 엔트로피 안정성이 양의 밀도·압력을 보장하지는 않는다 — positivity preserving은 별개의 장치가 필요하다.

8. 조건이 하지 못하는 일[편집]

엔트로피 조건을 만병통치약으로 오해하는 것도 흔한 함정이라 경계를 명확히 해 둘 필요가 있다.

양의 밀도·압력을 보장하지 않는다. 엔트로피 안정 도식이 강한 팽창이나 진공 근처에서 음의 압력을 뱉는 일은 얼마든지 있다. positivity preserving은 별도의 제한자로 강제해야 하고, 둘은 서로를 함의하지 않는다.
다차원 계의 유일성을 주지 않는다. 앞서 말한 대로 볼록적분 구성으로 엔트로피 조건을 만족하는 해를 무수히 만들 수 있다. “물리적으로 옳은 해”를 고르는 데는 조건이 하나 더 필요하다는 뜻이고, 그 조건이 무엇인지는 미해결이다.
모형이 틀린 것을 고쳐 주지 않는다. 상변화가 있는 계, 비이상 상태방정식, 비고전(nonclassical) 충격이 나오는 계에서는 고전적 엔트로피 조건이 실제 물리와 다른 해를 고르는 사례가 있다. 이때는 충격 내부의 점성 프로파일이 실제로 존재하는지를 묻는 점성 프로파일 판정 같은, 원래 물리에 더 가까운 기준으로 갈아타야 한다.
수치 소산과 물리 소산은 다르다. 엔트로피 부등식은 소산의 부호만 강제하고 크기는 강제하지 않는다. 부등식을 만족하면서도 충격을 열 셀에 걸쳐 뭉개는 도식은 얼마든지 있고, 그건 “안정하지만 틀린” 결과다.

9. 실무 체크리스트[편집]

팽창 충격 진단. 소드 문제나 준정지 충격 문제를 돌렸을 때 팬 안에 계단이 하나 서 있으면 십중팔구 음속점 문제다. 격자를 두 배로 줄여도 계단 폭이 그대로면 확정.
엔트로피 픽스는 최소한으로. $\delta$ 를 물리 문제의 특성 속도 스케일에 묶어라. 절대값 상수로 박아 넣으면 무차원화가 바뀌는 순간 조용히 틀린다.
비볼록 플럭스를 의심하라. 다상유동, 상전이, 교통류, 비이상 기체 상태방정식에서는 플럭스가 볼록이 아니다. 락스 조건만 코딩해 두면 복합파를 놓치고, 이때는 올레이닉 껍질 구성으로 참값을 직접 그려 검증해야 한다.
엔트로피 잔차를 모니터링하라. 셀별 수치 엔트로피 생산량을 후처리로 뽑으면 충격 위치와 과도한 소산 영역이 한눈에 보인다. 실질적으로 공짜인 오차 지시자라 적응 격자 세분화의 지표로도 쓴다.

10. 관련 문서[편집]

이 계수 맞추기는 락스 조건이 왜 “그럴듯한 물리적 서사”가 아니라 잘 정의된 문제의 조건인지를 보여 준다. 특성선이 하나 모자라면 미결정, 하나 많으면 과결정이다. 선형대수 시험 문제와 똑같은 논리인데 무대만 시공간일 뿐. ↩
이 구성 때문에 스칼라 비볼록 리만 문제는 “종이와 자만 있으면 정확해가 나오는” 몇 안 되는 비선형 PDE 문제가 됐다. 학부생에게 볼록 껍질 그리기를 시켜 보면 껍질을 잘못 잡아 충격 하나를 통째로 빠뜨리는데, 그게 정확히 코드가 하는 실수이기도 하다. ↩
변수 이중화는 이름 그대로 $u(x,t)$ 와 $v(y,s)$ 를 서로 다른 변수로 놓고 두 부등식을 더한 뒤 델타족으로 대각선에 밀어붙이는 기법이다. 처음 보면 마술 같고 두 번째 봐도 마술 같다. 1970년 논문 한 편으로 스칼라 보존법칙 이론이 사실상 종결됐다는 평가를 받는다. ↩