[논문 리뷰] Density profile of a self-gravitating polytropic turbulent fluid in the context of ensembles of molecular clouds

Sava Donkov, Ivan Zh. Stefanov외 2명|arXiv (Cornell University)|2020. 12. 11.

Astrophysics and Star Formation Studies참고 문헌 67인용 수 7

한 줄 요약

이 논문은 구형 대칭 분자의 자기중력적, 다이내믹한 유체에 대해, 핵부근에서 등온(Γ = 1)에서 딱딱한 다이내믹 상태(Γ > 1)로의 상태 방정정식 전이를 가정하여 밀도 프로파일 방정식을 유도한다. 총 단위 질량당 에너지가 유동을 따라 보존되며, 베르누이 유사 방정식이 도출되며, 그 중 하나의 해는 Γ = 4/3일 때 밀도 프로파일 ρ ∝ r⁻³를 얻는다. 이는 고밀도 영역에서 관측되거나 시뮬레이션된 거듭제곱 꼬리와 일치한다.

ABSTRACT

We obtain an equation for the density profile in a self-gravitating polytropic spherically symmetric turbulent fluid with an equation of state $p_{ m gas}\propto ho^\Gamma$. This is done in the framework of ensembles of molecular clouds represented by single abstract objects as introduced by Donkov et al. (2017). The adopted physical picture is appropriate to describe the conditions near to the cloud core where the equation of state changes from isothermal (in the outer cloud layers) with $\Gamma=1$ to one of `hard polytrope' with exponent $\Gamma>1$. On the assumption of steady state, as the accreting matter passes through all spatial scales, we show that the total energy per unit mass is an invariant with respect to the fluid flow. The obtained equation reproduces the Bernoulli equation for the proposed model and describes the balance of the kinetic, thermal and gravitational energy of a fluid element. We propose as well a method to obtain approximate solutions in a power-law form which results in four solutions corresponding to different density profiles, polytropic exponents and energy balance equations for a fluid element. One of them, a density profile with slope $-3$ and polytropic exponent $\Gamma=4/3$, matches with observations and numerical works and, in particular, leads to a second power-law tail of the density distribution function in dense, self-gravitating cloud regions.

연구 동기 및 목표

• 핵부근에서 중력과 비이상 기체 거동이 지배하는 분자운의 밀도 구조를 모델링하기 위해.
• 밀도 확률 밀도 함수(PDF)의 두 번째 거듭제곱 꼬리가 고밀도 영역에서 관측되는 이유를 설명하기 위해.
• 유동하는 자기중력적 유체 층에서 운동에너지, 열에너지, 중력에너지가 균형을 이루는 정적인 해석적 프레임워크를 개발하기 위해.
• 이론적 밀도 프로파일을 관측 및 수치 결과에서 관찰된 프리스타르 및 프로스타르 핵의 거듭제곱 꼬리 기울기와 연결하기 위해.

제안 방법

• 구형 대칭과 정적 유동 조건 하에서 밀도 프로파일에 대한 비선형 적분 방정식을 수립한다.
• 외곽부에서는 등온(Γ = 1)에서 핵부근에서는 Γ > 1로 전이되는 다이내믹 상태 방정식 pgas ∝ ρ^Γ을 적용한다.
• 유동을 따라 단위 질량당 총 에너지가 보존되며, 운동에너지, 열에너지, 중력에너지가 균형을 이루는 베르누이 방정식과 유사한 방정식을 도출한다.
• 다양한 에너지 지배 조건을 가정한 근사 해를 도출하기 위해 ρ(l) ∝ l⁻ᵖ 형태의 거듭제곱 꼬리 가정을 제안한다.
• 유도된 시스템을 해결하여 운동에너지, 열에너지, 중력에너지 항의 상대적 중요도에 따라 네 가지의 별개의 해를 식별한다.
• 유동의 정적 상태 가정을 검증하기 위해, 유착 시간스케일 τ₀와 천체역학적 시간스케일을 비교하여 τ₀ ≪ τₐ임을 보여준다.

실험 결과

연구 질문

• RQ1 자기중력적, 다이내믹한 난류 유체가 분자운 핵부근에서 정적 유동 조건을 만족할 때 어떤 밀도 프로파일이 도출되는가?
• RQ2 등온에서 딱딱한 다이내믹 상태 방정식으로의 전이가 에너지 균형과 최종 밀도 구조에 어떤 영향을 미치는가?
• RQ3 관측된 밀도 PDF의 두 번째 거듭제곱 꼬리(기울기 q ≈ -1)가 일관된 에너지 균형 모델로 설명될 수 있는가?
• RQ4 단위 질량당 총 에너지 보존이 반경 방향 밀도 프로파일을 결정하는 데 어떤 역할을 하는가?
• RQ5 운동에너지, 열에너지, 중력에너지의 상대 기여가 최종 밀도 프로파일에 어떻게 영향을 미치는가?

주요 결과

• 단위 질량당 총 에너지가 유동을 따라 보존되며, 운동에너지, 열에너지, 중력에너지가 균형을 이루는 베르누이 유사 방정식이 도출된다.
• 거듭제곱 꼬리 가정에 기반한 네 가지 해석적 해가 유도되며, 각각 운동에너지, 열에너지, 중력에너지 항의 지배 정도에 해당한다.
• 한 해는 Γ = 4/3일 때 밀도 프로파일 ρ ∝ r⁻³를 얻으며, 이는 고밀도 영역에서 관측되거나 시뮬레이션된 거듭제곱 꼬리와 일치한다.
• 이 r⁻³ 프로파일은 PDF 기울기 q ≈ -1에 해당하며, 고밀도 영역에서 관측된 두 번째 거듭제곱 꼬리와 일치한다(예: Schneider et al. 2015c; Kritsuk et al. 2011).
• 시간스케일 분석을 통해 모델은 (준)정적 상태 가정을 정당화하며, 핵 질량 증가 시간스케일 τ₀가 천체역학적 시간스케일보다 훨씬 짧음을 보여준다.
• p = 3 및 Γ = 4/3인 해는 밀도 PDF의 두 번째 거듭제곱 꼬리를 기술하는 것으로 해석되며, 밀도 높은 하위기구에 느린 유착으로 인해 발생한다.