[논문 리뷰] Temporal evolution of magnetic molecular shocks I. Moving grid simulations

P. Lesaffre, J.-P. Chièze외 2명|ArXiv.org|2004. 07. 28.

Astrophysics and Star Formation Studies참고 문헌 45인용 수 33

한 줄 요약

이 논문은 평균 자유경로까지 해상도를 확보하는 적응형 이동 격자와 은밀한 스킴을 사용하여 시간에 따라 변하는 1차원 다유체 자기유체역학적 시뮬레이션을 통해 자기장이 있는 분자 충격파를 연구한다. 이는 장시간에 걸친 안정된 C형 충격파의 존재를 검증하고, 음속점 형성 이전에 C형 충격파의 최대 속도가 발생함을 밝히며, H₂ 분해/재결합과 충돌성 이온화에 의해 유도되는 진동 현상이 자기장에 의해 억제됨을 규명한다.

ABSTRACT

We present time-dependent 1D simulations of multifluid magnetic shocks with chemistry resolved down to the mean free path. They are obtained with an adaptive moving grid implemented with an implicit scheme. We examine a broad range of parameters relevant to conditions in dense molecular clouds, with preshock densities between 10^3 and 10^5 cm-3, velocities between 10 and 40 km/s, and three different scalings for the transverse magnetic field: B=0,0.1,1 μG \sqrt{n.cm3}. We first use this study to validate the results of Chièze, Pineau des Forêts & Flower (1998), in particular the long delays necessary to obtain steady C-type shocks, and we provide evolutionary time-scales for a much greater range of parameters. We also present the first time-dependent models of dissociative shocks with a magnetic precursor, including the first models of stationary CJ shocks in molecular conditions. We find that the maximum speed for steady C-type shocks is reached before the occurrence of a sonic point in the neutral fluid, unlike previously thought. As a result, the maximum speed for C-shocks is lower than previously believed. Finally, we find a large amplitude bouncing instability in J-type fronts near the H2 dissociation limit (u ~ 25-30 km/s), driven by H2 dissociation/reformation. At higher speeds, we find an oscillatory behaviour of short period and small amplitude linked to collisional ionisation of H. Both instabilities are suppressed after some time when a magnetic field is present. In a companion paper, we use the present simulations to validate a new semi-analytical construction method for young low-velocity magnetic shocks based on truncated steady-state models.

연구 동기 및 목표

• 높은 공간적·시간적 해상도로 조밀한 분자운에서 다유체 자기장 충격파의 시간적 진화를 모델링하기 위해.
• Chièze 등(1998)이 제시한 시간 스트레칭 방법에 의해 예측된 진화 시간스케일과 안정상태 행동을 검증하기 위해.
• 자기장 영향 하에서 분해성 및 부분적으로 이온화되는 충격파, 즉 CJ형 및 정적 충격파의 형성 과정을 조사하기 위해.
• H₂ 화학 반응, 이온화 및 냉각이 진동 불안정성의 원동력이 되는 메커니즘과 자기장이 이를 억제하는 방식을 조사하기 위해.
• 시간에 따라 변화하는 데이터를 기반으로 젊고 저속의 자기장이 있는 충격파에 대한 반분석 모델을 구축하기 위한 기초를 마련하기 위해.

제안 방법

• 충격파를 평균 자유경로까지 해상도를 확보하기 위해 적응형 이동 격자 알고리즘과 은밀한 시간 적분 스킴을 구현하기 위해.
• 중성, 이온, 전자 유체를 포함하는 다유체 MHD 프레임워크를 사용하며, 화학 반응과 에너지 방정식을 상호 연결함으로써.
• 세부적인 H₂ 준위 분포와 냉각을 포함하며, H₂, CO, H₂O, OH 분자 냉각에 대해 LVG 처리를 고려한 복합 냉각 모델을 적용함으로써.
• 밀도 10³–10⁵ cm⁻³, 속도 10–40 km s⁻¹, 횡방향 자기장은 B ∝ √n 비례로 스케일링된 광범위한 파라미터 공간에서 충격파를 시뮬레이션함으로써.
• 화학 평형을 가정하지 않고 비평형 화학 반응과 냉각을 직접 처리함으로써.
• 충격파 전면과 급격한 변화가 일어나는 영역에 집중된 격자 세분화를 통해 불연속성과 전구조를 정밀하게 해상도 확보함으로써.

실험 결과

연구 질문

• RQ1 실제 분자운 조건에서 C형 충격파가 안정상태에 도달하기 위해 필요한 진화 시간스케일은 무엇인가?
• RQ2 자기장 존재 시 C형 및 J형 충격파의 형성과 안정성, 특히 분해 영역에서의 영향은 어떠한가?
• RQ3 충격파 전면의 진동 불안정성은 어떤 물리적 메커니즘이 유도하며, 유입 속도와 밀도에 따라 어떻게 달라지는가?
• RQ4 중성 유체에서 음속점 형성 이전에 안정된 C형 충격파의 최대 속도가 발생하는 지점은 어디인가?
• RQ5 H₂ 분해/재결합과 수소의 충돌성 이온화가 일시적인 진동을 유도하는 정도는 어떠하며, 이러한 현상은 자기장에 의해 억제될 수 있는가?

주요 결과

• 안정된 C형 충격파는 매우 긴 시간스케일이 필요하며, 최대 10⁵ 년 이상이 소요되며, 이는 이전의 시간 스트레칭 결과를 검증하고, 안정상태 모델이 초기 충격 진화를 반영하지 못할 수 있음을 확인한다.
• 안정된 C형 충격파의 최대 속도는 중성 유체에서 음속점 형성 이전에 발생하며, 이는 이전의 가정과 정반대되며, 이는 이전에 추정된 것보다 낮은 C형 충격파 속도 상한선이 존재함을 시사한다.
• 두 가지의 별개의 진동 불안정성이 확인됨: H₂ 분해/재결합에 의해 유도되는 큰 진폭의 반사 불안정성(약 u ≈ 25–30 km s⁻¹), 그리고 높은 속도에서 수소의 충돌성 이온화와 관련된 짧은 주기·작은 진폭의 진동.
• 자기장이 존재할 경우, 특히 자기장 전구가 형성되면서 이 진동이 시간이 지남에 따라 억제됨을 확인하여, 자기장이 충격파의 구조를 안정화시키는 것으로 나타남.
• 이 시뮬레이션은 분자 조건에서 정적 CJ형 충격파의 첫 번째 시간에 따라 변화하는 모델과 자기장 전구가 존재하는 분해성 충격파의 모델을 제공함.
• 격자 해상도와 수치적 파rameter는 진동 진폭과 주기를 영향을 미치며, 이는 수치적 오류가 불안정성 특성 분석에 영향을 줄 수 있음을 시사하며, 정확한 예측을 위해서는 선형 안정성 분석이 필요함.