고압축성 자기유체역학을 위한 강인 수치 스키마

본 논문은 이상 자기유체역학(MHD) 방정식의 수치 해석에 있어 가장 어려운 고압축성·고마하 흐름을 안정적으로 풀기 위한 새로운 스키마를 제시한다. 저자들은 FLASH 코드의 기존 MHD 모듈을 개조하여, 엔트로피 안정성을 보장하는 HLL‑type 근사 리만 솔버(HLL3R, HLL5R)와 양성 보존을 위한 특수 제한자를 적용한 MUSCL‑Hancock 재구성 방식을 도입하였다. 핵심 아이디어는 다음과 같다. 1. **엔트로피 안정성**: HLL3R/HLL5R는 이완 모델을 기반으로 하며, Chapman‑Enskog 전개를 통해 엔트로피 불등식(ρs)_t+(ρus)_x≥0을 만족한다. 이는 전통적인 Roe 솔버가 음의 밀도·내부 에너지를 생성할 위험을 근본적으로 차단한다. 2. **양성 보존**: 원시 변수(ρ, u, B, p)를 이용한 재구성에 MC 제한자를 적용하고,