압축성 다상유동 수치해석이 단상 유동과 근본적으로 다른 이유, 그리고 이 분야가 왜 아직도 활발한 연구 주제인지를 정리합니다.

압축성 다상유동이란?#

유체역학에서 **다상유동(multiphase flow)**은 두 가지 이상의 상(phase)이 공존하는 흐름을 말합니다. 물과 공기, 연료 액적과 고온 가스, 수중 폭발에서의 기체 버블 등이 대표적입니다.

여기에 **압축성(compressibility)**이 더해지면 문제가 급격히 복잡해집니다. 각 상이 서로 다른 상태방정식(EOS)을 가지며, 계면(interface)을 사이에 두고 충격파(shock wave)가 발생하고 반사, 투과, 굴절됩니다.

단상 압축성 유동과의 차이#

단상 Euler 방정식은 다음과 같습니다:

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

여기서 보존 변수 $\mathbf{U} = [\rho, \rho u, E]^T$ 이고, 단일 EOS (예: ideal gas $p = (\gamma - 1)\rho e$ )로 시스템이 닫힙니다.

다상유동에서는 계면을 경계로 EOS가 바뀝니다. 기체 영역에서는 $\gamma = 1.4$ 인 ideal gas를, 액체 영역에서는 stiffened gas EOS를 쓴다고 하면:

p = (\gamma - 1)\rho e - \gamma p_\infty

문제는 수치적으로 계면 근처에서 두 EOS가 **혼합(mixing)**될 때 발생합니다. 보존 변수 기반으로 단순히 평균을 내면 비물리적인 압력 진동(pressure oscillation)이 생깁니다.

가장 유명한 문제입니다. 1994년 Abgrall이 지적한 이래, 이 문제를 해결하는 것이 다상유동 수치 기법의 출발점이 되었습니다.

"압축성 다상유동 수치해석의 역사는 곧 압력 진동과의 싸움의 역사이다."

단순한 보존형 유한체적법(conservative finite volume)은 계면에서 $\gamma$ 가 중간값을 갖게 되면서 spurious oscillation을 만들어냅니다.

실제 문제에서는 ideal gas만으로 부족합니다. 물은 Tait EOS 또는 stiffened gas EOS, 고압 상태에서는 Jones-Wilkins-Lee (JWL) EOS, 반응성 유동에서는 화학 반응과 결합된 EOS가 필요합니다.

각 EOS에 대해 Riemann solver를 일반화해야 하므로 구현 복잡도가 크게 올라갑니다.

계면을 어떻게 다룰 것인가는 크게 두 가지 접근이 있습니다:

접근법	방법	장점	단점
Sharp interface	Ghost Fluid, Front Tracking	계면이 선명	위상 변화 처리 어려움
Diffuse interface	5-equation, 6-equation model	위상 변화 자연스럽게 처리	계면 스미어링

압축성 다상유동 시뮬레이션은 다양한 공학 분야에서 핵심입니다:

이 분야의 핵심 도구인 Riemann 문제와 Godunov 기법을 다룹니다. 단상 Euler 방정식의 Riemann 문제를 먼저 이해한 뒤, 이것이 다상유동으로 어떻게 확장되는지 살펴보겠습니다.