오일러 방정식 (비점성 유동)

전산유체역학 유체역학 마지막 수정: 2026-07-13 04:28:19

1. 개요[편집]

오일러 방정식 Euler Equations
분야	유체역학 × 전산유체역학 × 기체동역학
가정	비점성(inviscid), 비열전도
방정식 유형	쌍곡형(hyperbolic) 보존 법칙계
보존량	질량 · 운동량 · 에너지
대표 현상	충격파, 팽창파, 접촉 불연속

오일러 방정식(Euler equations)은 점성과 열전도를 무시한 이상 유체의 흐름을 지배하는 질량·운동량·에너지 보존 방정식계이다. 나비에-스토크스 방정식에서 점성 응력항과 열전도항을 통째로 지워버린 형태로, 1757년 레온하르트 오일러가 유도한 유체역학의 원조 방정식이다.¹ 점성이 사라지면서 방정식은 2계에서 1계로 내려앉고, 그 성질도 포물형이 섞인 형태에서 순수 **쌍곡형(hyperbolic)**으로 바뀐다. 이 쌍곡성이야말로 오일러 방정식의 모든 것 — 충격파와 팽창파가 등장하고, 정보가 유한한 속도로 특성선을 따라 전파되며, 리만 솔버라는 전용 도구가 필요해지는 이유다.

“점성을 왜 버리냐”고 물을 수 있다. 답은 간단하다. 고속 압축성 유동에서는 레이놀즈수가 워낙 커서 점성 효과가 얇은 경계층 안에만 갇히고, 바깥의 주 유동은 사실상 비점성이기 때문이다. 그 바깥 영역을 값싸게 푸는 것이 오일러 방정식의 존재 이유다.

2. 보존형 방정식[편집]

오일러 방정식은 보존형(conservation form)으로 쓸 때 가장 아름답다. 1차원의 경우:

\frac{\partial \mathbf{U}}{\partial t} + \frac{\partial \mathbf{F}(\mathbf{U})}{\partial x} = 0

여기서 보존량 벡터 $\mathbf{U}$ 와 플럭스 벡터 $\mathbf{F}$ 는 다음과 같다.

\mathbf{U} = \begin{bmatrix} \rho \\ \rho u \\ \rho E \end{bmatrix}, \qquad \mathbf{F} = \begin{bmatrix} \rho u \\ \rho u^2 + p \\ (\rho E + p)\,u \end{bmatrix}

$\rho$ 는 밀도, $u$ 는 속도, $p$ 는 압력, $E$ 는 단위 질량당 총에너지다. 미지수가 넷( $\rho, u, p, E$ )인데 방정식은 셋이라, 계를 닫으려면 상태방정식(이상기체는 $p = (\gamma-1)\rho e$ )이 하나 더 필요하다. 이 보존형은 유한체적법과 궁합이 완벽하다. 검사체적 경계로 드나드는 플럭스만 따지면 되고, 충격파를 가로질러도 보존량이 정확히 지켜지기 때문이다.²

3. 쌍곡성과 특성선[편집]

오일러 방정식을 준선형 형태로 쓰면 플럭스 자코비안 $\mathbf{A} = \partial \mathbf{F}/\partial \mathbf{U}$ 가 등장한다. 이 자코비안 행렬의 고유값이 바로 정보가 전파되는 속도다.

\lambda_1 = u - a, \qquad \lambda_2 = u, \qquad \lambda_3 = u + a

여기서 $a = \sqrt{\gamma p / \rho}$ 는 음속이다. 세 고유값이 모두 실수이고 고유벡터가 완비되어 있으므로 계는 엄밀히 쌍곡형이다. 이 세 특성 속도는 각각 세 종류의 파동에 대응한다. $u \pm a$ 는 음향파(충격파·팽창파), $u$ 는 접촉 불연속(contact discontinuity)이다. 정보가 유한 속도로 특성선을 따라 흐른다는 사실은 곧 CFL 조건이 시간 스텝을 지배한다는 뜻이기도 하다.³

4. 불연속해와 리만 문제[편집]

오일러 방정식의 가장 극적인 특징은 매끄러운 초기조건에서도 유한 시간에 불연속해(충격파)가 생겨난다는 점이다. 점성이 없으니 이 불연속을 매끄럽게 뭉갤 소산 메커니즘이 방정식 안에 없다. 그래서 고전적 미분 방정식으로는 해가 정의되지 않고, 적분형 보존 법칙을 만족하는 **약해(weak solution)**와 물리적으로 옳은 해를 골라내는 엔트로피 조건이 필요하다.

계면 양쪽에 서로 다른 두 상태가 주어졌을 때 그 사이에서 벌어지는 파동 구조를 푸는 문제를 리만 문제(Riemann problem)라 한다. 그 해는 충격파, 팽창팬(rarefaction fan), 접촉 불연속의 조합으로 구성되며, 이를 계산하는 것이 리만 솔버다. Godunov, Roe, HLLC 같은 이름들이 여기서 등장하며, 현대 압축성 유동 코드의 심장부를 이룬다.⁴

5. 수치해석과 응용[편집]

오일러 방정식을 푸는 코드는 대개 유한체적법 위에 리만 솔버와 고차 재구성(MUSCL, WENO), TVD 차분 도식을 얹어 충격파 근처의 진동을 억제한다. 실제 활약 무대는 다음과 같다.

초음속·극초음속 공기역학: 항공기·미사일·재진입체 주위의 충격파 구조. 경계층 바깥의 비점성 코어를 오일러로 푼다.
로켓·터빈 내부 유동: 노즐 팽창, 충격 셀 등 압축성 효과가 지배적인 영역.
폭발·천체물리: 충격파 전파, 초신성 잔해 등 강한 불연속이 핵심인 문제.
포텐셜 유동의 상위 호환: 비회전 가정까지 버리고 회전·엔트로피 생성을 허용하므로, 포텐셜 유동으로는 못 잡는 충격 후 유동을 다룬다.

6. 비점성 가정의 한계[편집]

공짜 점심은 없다. 점성을 버린 대가는 명확하다. 첫째, 항력 예측이 반쪽짜리다. 마찰 항력은 물론이고, 경계층 박리로 생기는 압력 항력조차 제대로 못 잡는다. 극단적으로는 무한히 긴 물체 주위 정상 비점성 유동에서 항력이 0으로 나오는 달랑베르의 역설(d’Alembert’s paradox)이 튀어나온다.⁵ 둘째, 경계층·유동박리·난류를 근본적으로 표현하지 못한다. 그래서 실무에서는 오일러 해를 경계층 방정식과 결합(viscous-inviscid interaction)하거나, 애초에 나비에-스토크스 방정식 전체를 푼다. 오일러 방정식은 “빠르고 값싸지만 벽 근처는 못 믿는” 도구로, 개념 설계와 충격파 물리 이해에서 여전히 제 몫을 한다.

7. 관련 문서[편집]

오일러는 이 방정식을 나비에나 스토크스보다 근 한 세기 앞서 유도했다. 점성항을 채워 넣어 나비에-스토크스가 완성된 것이 19세기 중반이니, 역사적으로는 오일러 방정식이 형이고 나비에-스토크스가 동생인 셈이다. ↩
보존형(conservation form)과 비보존형(primitive form)은 매끄러운 해에서는 수학적으로 동치지만, 충격파를 가로지르면 이야기가 달라진다. 비보존형으로 이산화하면 충격파 속도가 틀리게 나온다. “충격파를 풀 거면 무조건 보존형”이 이 바닥 철칙이다. ↩
정보가 유한 속도로 퍼진다는 건 곧 “내 옆 칸만 보면 다음 스텝을 계산할 수 있다”는 뜻이다. 이 국소성이 명시적 시간 적분과 리만 솔버를 가능하게 한다. 반대로 비압축성 유동은 압력이 무한 속도로 전파되는 타원형이라 이 사치를 누리지 못한다. ↩
Godunov(1959)가 리만 문제를 셀 경계마다 푸는 아이디어를 처음 냈을 때만 해도 계산 비용이 살인적이었다. Roe(1981)의 근사 리만 솔버가 나오면서 실용화되었고, 지금은 HLLC가 정확도·강건성의 국룰 자리를 차지하고 있다. ↩
달랑베르의 역설: 점성 없는 비압축·비회전 유동에서 물체가 받는 항력이 0이라는, 경험과 정면으로 배치되는 결과. 이 역설을 해소하려면 아무리 작아도 점성이 만드는 경계층과 그로 인한 박리가 필요하다. 점성을 우습게 보면 항력이 통째로 증발한다는 교훈이다. ↩