Centrifugal instability of compressible flows and the hydrodynamic stability of accretion disks

기술 요약: 압축성 흐름의 원심 불안정성과 강착 원반의 유체역학적 안정성

문제 제기
회전하는 흐름에서의 원심 불안정성은 전통적으로 레이leigh(Rayleigh) 기준에 의해 지배되며, 이는 특정 각운동량이 반경에 따라 감소할 때 불안정성이 발생한다는 원리이다. 이 기준은 비점성, 비압축성 유체를 위해 도출되었다. 그러나 최근 압력에 의해 지지되는 압축성 상대론적 회전(Gourgouliatos & Komissarov 2018b)에 대한 분석은 이 기준이 흐름의 마하수($M$)에 명시적으로 의존하는 일반화된 형태를 가질 수 있음을 시사했다. 뉴턴 역학적 한계에서 이 "KG 기준"은 $\Psi = \rho (v_\phi R)^2$일 때, $d \ln \Psi / d \ln R < M^2$이면 흐름이 불안정하다고 가정한다. 이는 고전적인 레이leigh 기준(즉, $d(v_\phi R)^2/dR > 0$) 하에서는 안정적인 흐름이라 할지라도, 마하수가 충분히 높다면 불안정해질 수 있음을 의미한다.

이러한 이론적 예측은 천체 강착 원반에 대한 기존의 이해에 도전한다. 얇은 케플러 원반는 초음속 속도로 회전함에도 불구하고 레이leigh 기준을 만족하기 때문에 유체역학적으로 안정적인 것으로 널리 간주된다. 저자들은 이러한 계의 높은 마하수와 압축성이 결래, 이들이 실제로 원심 불안정성에 취약해질 수 있다는 점을 제기하며, 이는 비압축성을 가정하거나 마하수에 대한 명시적 의존성을 무시하는 표준 안정성 분석에서 흔히 간과되는 요소이다.

방법론
저자들은 뉴턴 버전의 KG 기준을 검증하기 위해, 비상대론적, 천음속 회전 흐름에 대한 축대칭 수치 시뮬레이션을 수행하였다. 시뮬레이션은 이상 기체 상태 방정식($\gamma = 5/3$)을 사용하는 AMR-VAC 코드를 사용하여 이상 압축성 유체 역학 방정식(연속 방정식, 운동량 방정식, 에너지 방정식)을 풀었다.

시뮬레이션 설정은 다음과 같다:

• 초기 평형: 균일 밀도($\rho = \rho_0$)와 중심 핵에서 강체 회전에서 외부로의 거듭제곱 법칙 붕괴($v_\phi \propto R^{-n}$)로 전이되는 회전 프로파일을 가진 구성.
• 압력 균형: 압력 분포는 반경 방향 힘의 균형 방정식 $dp/dR = \rho v_\phi^2 / R$을 만족하도록 도출되었다.
• 섭동: 잠재적 불안정성을 유발하기 위해 작은 방위각 속도 섭동($\delta v_\phi = \sin(kz)$)으로 평형 상태를 교란하였다.
• 매개변수 공간: 저자들은 거듭제곱 지수 $n$(0.2에서 1.0 사이)과 중심 압력 $p_0$를 변화시켜 중심 경계에서의 마하수($M_c$)가 다른 모델들을 생성하였다. 이를 통해 레이leigh 기준에 따르면 안정적이지만 KG 기준에 따르면 불안정한 구성을 테스트할 수 있었다.

주요 기여 및 결과

1. KG 기준의 검증: 수치적 결과는 이론적인 KG 기준과 완벽하게 일치함을 보여주었다. 불안정할 것으로 예측된 모델들(여기서 $2(1-n) < M^2$)은 섭동의 급격한 성장을 보인 반면, 안정적이라고 예측된 모델들은 안정성을 유지하였다. 이는 압축성 흐름에서의 원심 불안정성이 마하수에 강하게 의존하며, 이는 비압축성 한계에서는 나타나지 않는 특징임을 확인시켜 준다.
2. 비선형 진화: 시뮬레이션은 이 불안정성의 비선형 발달을 드러냈다. 섭동은 불안정한 영역의 안쪽 경계에서 가장 빠르게 성장하였다. 불안정성이 진화함에 따라, 흐름은 중심 핵으로부터 멀어지는 반경 방향 스트림(손가락 모양)을 형성하였다. 이 스트림들은 주변 유체와 상호작용하여 켈빈-헬름홀츠 불안정성의 특징인 버섯 모양의 구조와 코쿤(cocoon)을 형성하였다. 이 과정은 넓은 혼합을 보이는 난류 층으로 이어졌으며, 이는 초기 선형 불안정 영역을 크게 벗어나 확장되었다.
3. 강착 원반 안정성의 재평가: 저자들은 압력 힘과 중심의 컴팩트 천체 중력을 모두 포함하도록 휴리스틱 안정성 분석을 확장하였다. 그들은 마하수에 대한 명시적 의존성을 유지하는 일반적인 불안정성 기준을 도출하였다.
   • 중력이 사라지는 극한에서는, 이 기준은 압력 지지 회전에 대한 KG 기준으로 축소된다.
   • 회전이 사라지는 극한에서는, 이 기준은 대류 불안정에 대한 슈바르츠실트(Schwarzschild) 기준으로 축소된다.
   • 일반적인 경우, 도출된 기준은 마하수 의존성을 명시적으로 드러내지 않는 표준 솔베르그-호일랜드(Solberg-Høiland) 기준과 수학적으로 동일하다.

의의 및 주장
본 논문은 회전이 매우 초음속임에도 불구하고 강착 원반의 유체역학적 안정성이 종종 원래의 레이leigh 기준을 사용하여 설명된다고 주장한다. 마하수를 명시적으로 포함하는 기준을 도출하고 이것이 솔베르그-호일랜드 기준과 동등함을 보여줌으로써, 저자들은 얇은 케플러 원반의 안정성이 단순히 비압축성 레이leigh 조건을 만족하기 때문이 아니라, 솔베르그-호일랜드가 설명하는 원심력과 대류력 사이의 특정한 균형의 결과임을 명확히 한다.

저자들은 KG 기준이 압력 지지 흐름이나 제트(중력이 미미한 경우)에 유용한 정식화(formulation)를 제공하지만, 일반적인 솔베르그-호일랜드 정식화에서는 마하수 의존성이 엔트로피 및 각운동량 기울기에 효과적으로 흡수된다고 결론짓는다. 결과적으로, 등엔트로피 흐oli(isentropic flows, 엔트로피가 일정한 흐름)의 경우, 이 기준은 바로트로픽(barotropic) 유체에 대한 레이leigh(1917)의 분석과 일치하게 원래의 레이leigh 기준으로 축소된다. 이 연구는 높은 마하수가 본래 안정적이었을 압력 지지 흐름을 불안정하게 만들 수 있음을 확인시켜 주며, 이는 천체 제트의 역학 및 초음속 바람과 주변 플라즈마의 상호작용과 관련된 발견이다. 다만, 저자들은 이 불안정성에 의해 촉발되는 비선형 혼합을 완전히 규명하기 위해서는 3D 시뮬레이션이 필요하다고 언급하였다.