보존 전하를 포함한 점성 Gubser 흐름의 반정밀 해법 및 유체 시뮬레이션 벤치마크

초록

본 논문은 보존 전하(B, S, Q)를 갖는 점성 Gubser 흐름에 대한 반정밀 해를 제시한다. 온도와 화학 퍼텐셜의 시간·공간 진화를 분석하고, 일정 에너지 밀도 기준의 동결면과 그 법선 벡터를 계산한다. 두 종류의 방정식 상태(EoS1, EoS2)를 비교한 뒤, 새로 개발한 Smoothed Particle Hydrodynamics 기반 코드 CCAKE와의 수치적 일치를 검증한다. 결과는 반정밀 해와 수치 해가 높은 정밀도로 일치함을 보여, 대규모 화학 퍼텐셜을 포함한 상대론적 유체 시뮬레이션의 신뢰성 검증에 유용한 기준을 제공한다.

상세 분석

이 연구는 Gubser 흐름의 대칭성을 활용해 복잡한 점성 유체 방정식을 de Sitter 좌표계로 변환함으로써, 편미분 방정식 집합을 상미분 방정식으로 축소한다는 핵심 아이디어를 갖는다. 특히, 보존 전하가 존재할 경우 에너지‑운동량 텐서와 전하 전류가 직접적으로 결합되지 않지만, 방정식 상태(EoS)를 통해 온도와 화학 퍼텐셜이 간접적으로 연결된다. 저자들은 두 가지 전형적인 방정식 상태를 선택했는데, 첫 번째는 무질량 쿼크‑글루온 플라즈마의 자유 입자 가정에 기반한 식(EoS1)이며, 두 번째는 간단한 선형 결합 형태(EoS2)이다. EoS1은 온도와 화학 퍼텐셜의 4차항까지 포함해 비선형성을 제공하고, EoS2는 T와 μ라는 상수를 도입해 상대적인 스케일을 조정한다.

점성 항은 2계 이완 이론(Israel‑Stewart) 형태의 전이 방정식으로 기술되며, 이완 시간 τ_R은 η/(E+P)와 비례하도록 설정한다. 여기서 η는 η/s = 0.2라는 고정값을 사용해 스케일링을 단순화한다. Gubser 흐름의 대칭성으로 인해 전하 확산이 금지되고, 전하 밀도 n_Y는 단순히 확장률 θ에 의해 감소한다(D n_Y + n_Y θ = 0). 따라서 n_B, n_S, n_Q는 서로 독립적으로 진화한다.

핵심 수식은 온도와 화학 퍼텐셜의 진화 방정식(20a‑b, 22a‑b)이며, 여기서 f(μ̂, T̂) 함수는 EoS1에 따라 복잡한 비선형 형태를 띤다. 이 방정식들은 tanh ρ 항을 통해 de Sitter 좌표계에서의 시간 의존성을 반영한다. 또한, 전단 응력의 축소 변수 (\hat{\bar\pi})는 자체적인 Riccati 형태 방정식(20c)으로 진화한다.

동결면은 일정 에너지 밀도 E_FO를 기준으로 정의되며, F(τ, r) = E_FO − E(τ, r) = 0인 곡면의 법선 벡터 n^μ = ∂^μF 로 구한다. 이때 τ와 r에 대한 편미분을 통해 n^τ와 n^r을 얻고, 정규화된 법선은 유동의 유체 흐름과 직교한다. 이러한 동결면 정의는 화학 퍼텐셜이 큰 경우 온도 감소를 자연스럽게 반영한다는 장점이 있다.

수치 검증에서는 CCAKE 코드가 SPH 방식을 사용해 입자 기반으로 유체를 전개한다. 초기 조건은 반정밀 해와 동일하게 설정하고, 격자 크기와 스무딩 길이를 다양하게 조정해 수렴성을 테스트한다. 결과는 온도·화학 퍼텐셜 프로파일, 전단 응력, 그리고 동결면 위치 모두에서 오차가 1 % 이하로 일치함을 보여준다. 특히, 오프‑대각 전단 성분(π^{τr} 등)에서 약간의 차이가 나타났지만, 전체적인 흐름 구조와 동결면 형태는 거의 동일했다.

이 논문의 주요 기여는 (1) 보존 전하를 포함한 점성 Gubser 흐름의 반정밀 해를 제공함으로써, 복잡한 화학 퍼텐셜 환경에서도 검증 가능한 기준을 마련했다는 점, (2) 두 종류의 방정식 상태가 해에 미치는 영향을 정량적으로 분석했으며, (3) SPH 기반 코드의 정확성을 동결면까지 포함한 전반적인 물리량에 대해 검증함으로써, 향후 고밀도·고전하 물리(예: QCD 임계점 탐색, 중성자별 합성) 시뮬레이션에 바로 적용 가능한 벤치마크를 제시했다는 점이다.