[논문 리뷰] Non-axisymmetric wind-accretion simulations. II. Density gradients

M. Ruffert|arXiv (Cornell University)|1999. 03. 19.

Astro and Planetary Science인용 수 14

한 줄 요약

이 연구는 3차원 비축대칭 풍속 가속에 대해 전단 밀도 기울기가 있는 매질에서 압축성 유체역학 시뮬레이션을 수행하였으며, 높은 해상도를 확보하기 위해 중첩 격자 구조를 사용한 조각별 파라볼릭 방법(PPM)을 적용하였다. 주요 발견은 질량 가속도율은 균일한 매질과 비교해 거의 변화가 없지만, 각운동량 가속도는 불안정성에 의해 강하게 조절되며, 안정된 단계에서는 가용 각운동량의 약 70%까지 가속도가 이루어지고, 불안정한 단계에서는 변동성이 크며 평균가를 감소시키는 것으로 나타났다. 특히 기울기가 큰 경우 후행 흐름이 유지되는 경향이 있다.

ABSTRACT

The hydrodynamics of a variant of classical Bondi-Hoyle-Lyttleton accretion is investigated: a totally absorbing sphere moves at various Mach numbers (3 and 10) relative to a medium, which is taken to be an ideal gas having a density gradient (of 3%, 20% or 100% over one accretion radius) perpendicular to the relative motion. Similarly to the 3D models published previously, both with velocity gradients and without, the models with a density gradient presented here exhibit non-stationary flow patterns, although the Mach cone remains fairly stable. The accretion rates of mass, linear and angular momenta do not fluctuate as strongly as published previously for 2D models. No obvious trend of the dependency of mass accretion rate fluctuations on the density gradient can be discerned. The average specific angular momentum accreted is roughly between zero and 70% of the total angular momentum available in the accretion cylinder in the cases where the average is prograde. Due to the large fluctuations during accretion, the average angular momentum of some models is retrograde by up to 25%. Small gradients hardly influence the average accretion rates as compared to accretion from a homogeneous medium, while very large ones succeed to dominate and form an accretion flow in which the sense of rotation is not inverted.

연구 동기 및 목표

• 횡방향 밀도 기울기가 3차원 유체역학적 풍속 가속도에 미치는 영향을 조사하기 위해.
• 마하 수, 등온지수 및 기울기 강도가 질량 및 각운동량 가속도율에 미치는 영향을 규명하기 위해.
• 밀도 기울기가 있는 경우와 없는 경우의 가속도 역학을 비교하여 안정성 및 시간 평균 각운동량 이행을 평가하기 위해.
• 밀도 기울기가 흐름의 구조 및 각운동량 수송에 지배적인지 평가하며, 특히 속도 기울기 모델과 비교하여 분석하기 위해.

제안 방법

• 고해상도를 확보하기 위해 조각별 파라볼릭 방법(PPM)을 사용한 수치적 유체역학 해법을 적용하였다.
• 9단계의 격자 중첩을 통해 최소 격자 크기 δ = 1/256 R_A로 고해상도를 확보하여 가속체 주변에서 정밀한 해석을 수행하였다.
• 가속체는 반사율이 100%인 구형으로 모델링되었으며, R⋆ = 0.02 R_A의 반지름을 가지며, 수치적 특이점을 피하기 위해 부드러운 잠재력이 적용되었다.
• 마하 수(M∞ = 1.4, 3, 10), 밀도 기울기 강도(ερ = 0.03, 0.20, 1.00), 등온지수(γ = 5/3, 4/3, 1.01)를 변화시켜 시뮬레이션을 수행하였다.
• 시간에 따른 질량 및 각운동량 성분을 계산하였으며, 변동성을 분석하기 위해 통계적 측도(평균, 표준편차)를 사용하였다.
• 가속도 실린더의 이론적 각운동량을 기준으로 하여 가용 각운동량의 비율을 평가하였다.

실험 결과

연구 질문

• RQ1 횡방향 밀도 기울기가 3차원 풍속 가속도에서 안정성 및 시간 평균 질량 가속도율에 어떤 영향을 미치는가?
• RQ2 가속도 실린더의 가용 각운동량 중 얼마나 많은 비율이 가속도되는가? 이는 기울기 강도와 마하 수에 따라 어떻게 달라지는가?
• RQ3 밀도 기울기 모델과 이전에 연구된 속도 기울기 모델 간의 질량 및 각운동량 가속도율 변동성은 어떻게 비교되는가?
• RQ4 강한 밀도 기울기(ερ = 100%)가 흐름을 안정화시키고 후행 각운동량 가속도를 방지하는가?
• RQ5 최대 각운동량 가속도의 분석적 추정치가 일시적인 안정 단계에서의 가속도 행동을 예측하는 데 사용될 수 있는가?

주요 결과

• 작은 밀도 기울기(ερ = 0.03)의 경우, 일시적인 준안정 단계가 발생하여 특정 각운동량의 가속도가 이론적 최대값의 약 10% 이내로 유지되며, 이는 분석적 추정치와 양호한 일치를 보였다.
• 불안정성이 발생한 후 평균 특정 각운동량 가속도는 크게 감소하였으며, 후행 경우에서는 변동성이 평균을 최대 25% 감소시켰다.
• 평균 특정 각운동량 가속도는 가용 총량의 0%에서 70% 사이에서 변동하였으며, 대부분의 모델에서 전행 값이 지배적이었다.
• 질량 가속도율의 변동성은 이전 2차원 시뮬레이션보다 뚜렷하지 않았으며, 기울기 강도 증가와의 명백한 추세도 관찰되지 않았다.
• 극단적인 경우 ερ = 1.00에서는 가속된 각운동량이 항상 전행이었으며, 배경 각운동량이 크기 때문에 불안정성의 상대적 영향이 작았다.
• ερ = 1.00 모델에서는 입구 흐름의 특정 각운동량이 가속체 표면에서 케플러 한계를 초과하지만, 가속된 각운동량은 여전히 전행이면서 안정된 상태를 유지하였다.