GLM 기반 CGL 방정식의 엔트로피 안정 스키마와 발산 억제

본 논문은 플라즈마 흐름을 기술하는 Chew‑Goldberger‑Low(CGL) 방정식에 대한 새로운 수치 해법을 제시한다. CGL 방정식은 비보존형 항을 포함하는 하이퍼볼릭 PDE 시스템으로, 압력 텐서가 자기장 방향에 따라 회전하는 특성을 갖는다. 이러한 비등방향 압력은 전통적인 MHD 모델이 가정하는 등방향 압력과 달리, 플라즈마가 충돌이 거의 없고 국부 열평형이 성립하지 않을 때 필수적이다. 그러나 비보존 항은 수치 해석 시 경로 의존성을 초래하고, 특히 ∇·B=0 조건을 만족시키지 못하면 물리적 비현실성을 야기한다. 1. **GLM‑CGL 모델 구축** 저자들은 Generalized Lagrange Multiplier(GLM) 기법을 CGL 방정식에 적용해 GLM‑CGL 시스템을 도출한다. 스칼라 라그랑주 승수 Ψ와 청소 속도 c_h를 도입해 자기장 발산을 억제하는 추가 방정식(∂_tΨ + c_h^2∇·B + v·∇Ψ = 0)을 포함한다. 기존 GLM‑MHD와 달리, GLM‑CGL에서는 Ψ가 0에 수렴함에 따라 원래 CGL 방정식으로 복원된다. 또한, 총 에너지 정의에 Ψ^2/2 항을 추가해 에너지 보존을 유지한다. 2. **엔트로피 분석 및 비보존 항 재구성** 엔트로피 함수 S=−ρs와 엔트로피 플럭스 q_i를 정의하고, 시스템이 엔트로피 불안정 없이 진화함을 보이기 위해 비보존 항을 선택적으로 재배치한다. 구체적으로, 보존 항 중 일부를 비보존 형태로 전환해 V·(비보존 항)=0이 되도록 설계하였다. Lemma 3.1을 통해, Ψ와 B에 대한 GLM 항이 엔트로피 생산에 기여하지 않으며, ∇·B=0일 경우 엔트로피 방정식이 보존 형태임을 증명한다. 3. **시스템 대칭화와 고유값·고유벡터** 비보존 항을 포함한 전체 시스템을 Φ′(V)·(∇·B) 형태로 재작성함으로써, 엔트로피 변수와 결합된 비보존 항이 시스템을 대칭화 가능하게 만든다. 이를 통해 8개의 파동 속도(알파, c_h, 알파±c_h 등)를 도출하고, 엔트로피 스케일링된 오른쪽 고유벡터를 구해 고차 정확도 스키마에 활용한다. 4. **엔트로피 안정 수치 스키마 설계** 공간 이산화는 균일 격자에 대해 반정밀 차분 형태로 전개한다. 보존 항은 엔트로피 보존 수치 플럭스(중심 플럭스 + 엔트로피 스케일링된 확산 연산자)를 사용하고, 비보존 항은 적절한 차수의 중앙 차분으로 처리한다. 엔트로피 보존 조건을 만족하도록 수치 플럭스 f̂_i는 V·(f̂_i−f_i)=ΔΨ_i 관계를 만족한다. 고차 재구성을 위해 sign‑preserving reconstruction과 entropy‑scaled eigenvectors를 이용한 고차 확산 연산자를 설계하였다. 시간 적분은 SSP‑Runge‑Kutta 방법을 적용해 전체 스키마의 강인성을 확보한다. 5. **수치 실험** 1차원 및 2차원 테스트 케이스(오리엔테이션 전파, 충격 파동, 플라즈마 불안정, 복합 파동 상호작용)를 수행하였다. 모든 2차원 실험에서 GLM‑CGL 모델은 ∇·B 오류를 현저히 감소시켰으며, 엔트로피 감소 법칙을 만족하였다. 특히, GLM 없이 CGL 방정식을 직접 해석했을 때 발생하는 자기장 발산과 엔트로피 비보존 현상이 GLM‑CGL에서는 거의 사라졌다. 6. **결론** GLM 기법을 CGL 방정식에 적용하고, 비보존 항을 엔트로피 변수와의 내적이 0이 되도록 재구성함으로써, 엔트로피 안정성과 자기장 발산 억제를 동시에 만족하는 수치 스키마를 성공적으로 개발하였다. 제안된 방법은 고차 정확도와 강인성을 유지하면서 플라즈마 물리학에서 비등방향 압력 효과를 정확히 포착할 수 있다. 향후 연구에서는 비정형 격자와 복합 경계 조건에 대한 확장 및 실험적 검증을 진행할 계획이다.