[논문 리뷰] Sunspot simulations with MURaM -- I. Parameter study using potential field initial conditions

Markus Schmassmann, N. Bello González외 2명|arXiv (Cornell University)|2026. 01. 30.

Solar and Space Plasma Dynamics인용 수 0

한 줄 요약

이 논문은 MURaM 코드의 잠재필드 초기 조건을 사용하여 태양흑점 시뮬레이션의 매개변수 연구를 수행하고, 초기 자기장 강도, 플럭스 내용, 박스 크기 및 해상도를 다변량으로 변화시켜 흑점과 펜럼브라 형성 및 관측과의 비교를 평가한다.

ABSTRACT

Context. Existing sunspot simulations fail to reproduce the observed magnetic field distribution due to an artificially increased $B_{hor}$ at the upper boundary. Aims. We explore alternative ways to better reproduce the magnetic and dynamic properties of observed sunspots. Methods. We use the radiative MHD code MURaM. As initial conditions, we placed a potential magnetic field into small-scale dynamo simulations and used potential-field extrapolation at the top. Results. We find that: (1) Simulations with increasing initial magnetic field strengths (20, 40, 80, and 160 kG) show larger spots, umbrae, and penumbrae. (2) The penumbral-to-spot sizes are smaller than those measured in observed sunspots. (3) In none of the runs are pure Evershed (radially outward) flows. Instead, bi-directional flows with inflows in the inner penumbra and outflows in the outer penumbra were measured, similar to early observations of penumbra formation for runs with $\ge80$ kG at 96/32 km resolution, whereas runs with 40 kG or less showed pure inflows. (4) Simulations with 160 kG at 32/16 km resolution contain filaments with bi-directional and Evershed flows. (5) Simulations with fluxes $>10^{22}$ Mx show unrealistically strong fields in the umbra. (6) The best runs with 160 kG and $10^{22}$ Mx give realistic radial profiles of $B_z$ and $B_r$, although stronger fields than observed. (7) Increasing the width of the box and reducing the overall flux by subtracting a uniform opposing vertical field have little influence on internal spot dynamics and fields, but change the mean vertical field outside the spot. Conclusions. Simulations of small ($10^{22}$) sunspots with an initial potential field and intensified bottom magnetic field strength best reproduce observations of the initial stages of sunspot formation. Numerical resolution may be critical for achieving fully developed penumbrae.

연구 동기 및 목표

• 이전 설정과 비교하여 잠재필드 초기 조건이 흑점의 자기·유동 특성을 더 현실적으로 만들어 주는지 평가한다.
• 초기 바닥 자기장 세기(B0), 총 플럭스, 박스 크기, 해상도가 흑점/펜럼브라의 크기 및 역학에 어떻게 영향을 미치는지 결정한다.
• 펜럼브라 형성 중 관측 벤치마크와 비교하여 시뮬레이션된 자기장 프로필과 흐름 패턴을 비교한다.
• 관측된 에버셰드(Evershed)-유사 유출과 양방향 펜럼브라 흐름을 재현하는 조건을 식별한다.

제안 방법

• 열방사 MHD 코드 MURaM을 사용하고, 하단 경계는 개방, 상단 경계는 잠재 필드 조건으로 설정한다.
• 하단에서의 초기 Bz(r)를 잠재장으로 보정하기 위해 가우시안 혹은 균일하게 반대 방향의 플럭스 Bopp로 초기 축 자기장을 B0로 가하고 설정한다.
• 실험에 따라 매개변수를 변화: B0 = 20, 40, 80, 160 kG; FGauss = 1e22 Mx (또는 Bopp를 사용하는 Ftot); 박스 너비는 49.152 또는 98.304 Mm; 해상도는 96/32 km 또는 32/16 km; 그리고 여러 Bopp 값(0–300 G).
• τ=1에서 방위평균된 양(I, Bz, Br, γ)을 계산하고, 강도 등에서의 암흑점/검은점 경계를 강도 등으로부터 추출하며, 반경 방향 흐름 vz, vr 및 부호 혼합 거동을 분석한다.
• 활동의 시간평균 적합을 수행하여 대표 상태 값을 얻고, 관측 흑점과 비교한다(예: NOAA AR 11591).
• 고해상도 실행에서 비회색(non-grey) 복사전달을 포함하여 전달 모델링 및 측정 장비 효과를 가능하게 한다.

Figure 1: Quadrants of bolometric intensity maps for four selected simulation runs with different initial magnetic field strengths. From top-left clock-wise: $B_{0}=20$ kG, 40 kG, 80 kG and 160 kG. The magnetic flux is the same for all four runs: $F=10^{22}$ Mx.

실험 결과

연구 질문

• RQ1 잠재필드 초기 조건에서 바닥 자기장을 강화하는 시작이 흑점과 펜럼브라 특성을 전통적 설정보다 더 관측에 근접하게 만드는가?
• RQ2 초기 자기장 세기, 총 플럭스, 도메인 크기 및 수치 해상도가 펜럼브라 형성 중 흑점의 크기, 펜럼브라 범위, 흐름 패턴(Evershed-유사 유출 대 유입)에 어떻게 영향을 미치는가?
• RQ3 시뮬레이션된 흑점의 자기장 프로필(Bz, Br)과 그 반경 의존성은 관측 벤치마크와 어떻게 비교되는가?
• RQ4 시뮬레이션이 양방향 펜럼브라 흐름과 현실적인 펜럼브라 필라멘트를 언제 어떤 조건에서 재현하는가?
• RQ5 고해상도나 더 큰 도메인이 펜럼브라 발달 및 플럭스 농축의 차이를 완화하거나 악화시키는가?

주요 결과

• 펜럼브라 범위는 초기 바닥 자기장 강도 B0가 증가함에 따라 커지며, 20 kG에서 펜럼브라가 없고 160 kG에서 길고 가느다란 필라멘트가 형성된다.
• FGauss > 1e22 Mx는 나중 단계에서 암흑점의 수직 자기장이 비현실적으로 강해지게 한다.
• BI-directional 펜럼브라 흐름은 B0 ≥ 80 kG인 시뮬레이션에서 나타나며, 펜럼브라 내부로는 유입이, 외부로는 유출이 관찰되며, 해상도가 높아질수록 이 현상이 두드러진다.
• 고해상도 실행(32/16 km)은 양방향 및 에버셰드 유사 흐름을 가진 필라멘트를 만들어내는 반면, 저해상도 실행은 더 강한 유입과 덜 현실적인 펜럼브라 패턴을 보인다.
• 1e22 Mx를 넘는 플럭스는 암흑점의 Bz를 과도하게 강하게 만들 수 있어, 이러한 설정 하에서 현실적인 시뮬레이션에 대한 상한을 시사한다.
• 160 kG, 1e22 Mx(및 고해상도 관련 실행)에서 radius에 따른 Bz와 Br의 최적 매칭 프로파일이 얻어지지만, 전체 흐름은 여전히 태양 관측보다 강할 수 있다.
• 더 넓은 계산 박스나 균일 반대 수직장 제거는 내부 점의 역학에는 큰 영향을 주지 않지만 외부의 평균 수직장에는 영향을 주어 수정한다.

Figure 2: Azimuthally averaged intensity $I/I_{\textrm{qs}}$ (solid lines) as a function of radius in Mm and fraction of the spot radius. The colours indicate the simulations listed in the inset legend. For the older runs (with names ending with 'o', red), the bolometric intensity is shown, and for