The (Limited) Effect of Viscosity in Multiphase Turbulent Mixing

문제 제기
다상 기체는 성간매질 (ISM), 은하권계면 (CGM), 은하단내부매질 (ICM) 과 같은 천체물리학적 환경에서 보편적으로 존재합니다. 이러한 시스템의 역학은 난류, 복사 냉각, 그리고 수송 과정 간의 상호작용에 의해 지배됩니다. 자기장이 난류 복사 혼합층 (TRMLs) 에서 켈빈 - 헬름홀츠 불안정성 (KHI) 을 억제하고 혼합을 변화시키는 역할은 연구되어 왔으나, 점성의 구체적인 영향은 여전히 논쟁의 대상입니다. 완전히 이온화된 플라즈마는 매우 높은 점성을 가질 것으로 예상되지만 (Spitzer 1962), 천체물리학적 시스템의 순 유효 점성은 불확실합니다. 중요한 질문은 점성이 난류와 혼합을 억제하여 난류 속도와 냉각률 사이의 "보편적" 스케일링 관계 (Fielding et al. 2020; Tan et al. 2021) 가 성립하는 데 필요한 조건을 무효화하는지 여부이며, 특히 카를로비츠 수 ($Ka < 1$) 가 층류 흐름을 시사하는 영역에서 그렇습니다.

방법론
저자들은 점성의 효과를 분리하기 위해 Athena 코드를 사용한 이상화된 수치 시뮬레이션을 수행했습니다. 연구는 두 가지 주요 단계로 진행되었습니다:

1. 단열 2 차원 시뮬레이션: KHI 를 억제하는 데 필요한 임계 점도 ($\eta_{Crit}$) 를 결정하기 위해, 저자들은 다양한 과밀도 ($\chi = \rho_{cold}/\rho_{hot}$) 와 마하 수 ($M$) 를 가진 전단 흐름을 시뮬레이션했습니다. 경계면의 인공적 효과를 피하기 위해 과밀도 의존성을 테스트하기 위해 평면 슬랩 설정을, 마하 수 의존성을 테스트하기 위해 평면 시트 설정을 사용했습니다. 저자들은 점성 시간 척도와 KHI 성장 시간 척도를 비교하여 $\eta_{Crit}$에 대한 분석적 추정을 도출했으며, 다중 밀도 유체에서의 운동량 차이를 고려하기 위해 운동 점도에 조화 평균을 적용했습니다.
2. 3 차원 난류 복사 혼합층 (TRMLs): $\chi=100$인 과밀도에 대해 유도된 임계 점도를 사용하여, 저자들은 명시적 점도를 포함한 TRML 의 3 차원 시뮬레이션을 수행했습니다. 그들은 $\eta_{Crit}$에 대한 점도 비율 (0% 에서 50% 까지) 과 섭동 파장을 조정하여 다름코ehler 수 ($Da = t_{turb}/t_{cool}$) 를 변화시켰습니다. 이를 통해 약한 냉각 ($Da < 1$) 과 강한 냉각 ($Da > 1$) 영역 모두를 탐색할 수 있었습니다. 시뮬레이션에는 복사 냉각 (Townsend 알고리즘) 이 포함되었으나 명시적 열전도는 포함되지 않았으며, 이는 $Ka < 1$을 보장했습니다.

주요 기여 및 결과

• 임계 점도 스케일링: 이 연구는 고과밀도 한계에서 KHI 를 억제하는 데 필요한 임계 점도가 과밀도에 비례하여 선형적으로 스케일링됨 ($\eta_{Crit} \propto \chi$) 을 확인했습니다. 이는 차갑고 밀도가 높은 매질이 운동량의 대부분을 운반하므로 불안정성을 억제하기 위해 더 강한 마찰이 필요하기 때문입니다. 저자들은 또한 높은 마하 수에서 압축성 효과가 점성 없이도 KHI 를 안정화시켜, 불안정성이 억제되는 임계 마하 수 ($M_{Crit}$) 를 정의함을 검증했습니다.
• **약한 냉각 영역 ($Da < 1$):** 이 영역에서는 냉각 시간 척도가 난류 혼합 시간보다 깁니다. 저자들은 점성이 난류를 강력하게 억제하여 흐름을 난류에서 층류로 전환시킨다고 발견했습니다. 결과적으로, 비점성 유체에 대해 유도된 난류 속도 ($u'$) 와 표면 밝기 ($Q$) 사이의 보편적 관계는 무효화됩니다. 그러나 역설적으로 총 볼로메트릭 광도는 거의 영향을 받지 않습니다. 저자들은 이를 "병목 효과"로 귀인하는데, 긴 냉각 시간이 혼합 메커니즘 (확산 vs 대류) 에 관계없이 전체 에너지 손실을 제한하기 때문입니다.
• 강한 냉각 영역 ($Da > 1$): 이 영역에서는 냉각 시간 척도가 점성 시간 척도보다 짧습니다. 저자들은 점성이 시스템 역학에 제한적인 영향을 미친다는 것을 발견했습니다. 점성이 존재함에도 불구하고, 냉각 과정은 계면에서의 온도와 속도 구배를 날카롭게 만들어 점성의 평탄화 효과를 효과적으로 상쇄합니다. 이로 인해 KHI 가 발달하여 비점성 사례와 유사한 난류 수준을 유지할 수 있습니다. 결과적으로 $u'$와 $Q$ 사이의 스케일링 관계는 유지되며, 시스템은 효과적으로 비점성처럼 행동합니다.
• 점성 대 냉각 경쟁: 저자들은 냉각이 점성을 지배하는 조건에 대한 분석적 조건 ($Da > \zeta$, 여기서 $\zeta$는 임계 점도의 비율) 을 도출했습니다. 그들의 수치 결과는 이 예측과 일치합니다. 강한 냉각 영역에서, 점성 실행 (약간 더 넓은 전면으로 인해) 에서의 중간 온도 기체의 양이 증가하여 비점성 실행에 비해 표면 밝기가 약간 증가합니다 (최대 약 2 배). 그러나 근본적인 스케일링 법칙은 유효합니다.

의의 및 함의
이 논문은 다상 난류 혼합에서 점성의 역할이 단순하게 예상된 것보다 더 제한적이라고 주장하며, 특히 은하 유출, CGM, ICM 등 많은 천체물리학적 환경을 특징짓는 강한 냉각 영역에서 그렇습니다.

• 관측적 진단: 이 발견은 선 폭 (난류) 과 방출도 (냉각) 를 연관 짓는 확립된 이론적 예측이 점성이 중요할 수 있는 시스템에서도 관측 데이터에 적용될 수 있음을 시사합니다. 단, 시스템이 강한 냉각 영역에 있어야 합니다.
• 서브그리드 모델: 이 결과는 대규모 우주론적 시뮬레이션에 사용되는 현재 다상 서브그리드 모델에 대한 물리적 근거를 제공합니다. 이러한 모델들은 종종 질량 전달률과 냉각 처방을 알리기 위해 비점성, 소규모 시뮬레이션에 의존합니다. 저자들은 냉각 지배 영역이 난류와 혼합의 점성 억제를 효과적으로 무효화하기 때문에 점성이 존재할 때에도 이러한 처방이 유효함을 입증했습니다.
• 천체물리학적 맥락: 이 연구는 점성이 대규모 형태 (예, 램 압력에 의해 벗겨진 꼬리) 에는 영향을 미칠 수 있지만, 빠른 냉각 한계에서는 에너지 소산이나 질량 전달률을 근본적으로 바꾸지 않는다고 강조합니다. 이는 은하 헤일로와 은하단에서 다상 기체의 진화를 해석하는 데 확신을 줍니다.

저자들은 그들의 이상화된 설정과 관련하여 주의점을 지적했는데, 구체적으로는 온도 의존성 스피처 점도 대신 일정한 점도를 사용한 점, 이방성 점성을 유도하는 자기장의 부재, 그리고 명시적 열전도의 부재가 그것입니다. 그들은 이러한 요소들을 천체물리학적 플라즈마에서의 수송 과정에 대한 이해를 더욱 정교하게 하기 위한 향후 연구의 방향으로 제안합니다.