가이드필드가 있는 상대론적 자기 재결합 연구

초록

본 논문은 상대론적 전자·양전자 두 유체 모델을 이용해, 가이드 필드가 존재하는 2차원 자기 재결합의 비선형 전개를 시뮬레이션한다. 가이드 필드가 강해질수록 방출 에너지의 구성이 엔탈피 중심에서 전자기 파동(Poynting) 중심으로 전환되며, 전하 분리와 가이드 필드 압축 현상이 두드러진다. 재결합을 유지하는 주요 메커니즘은 두 유체의 관성 효과이며, 전통적인 저항성보다 관성 항이 전기장 생성에 더 크게 기여한다.

상세 분석

이 연구는 상대론적 두 유체 방정식에 맥스웰 방정식을 결합한 완전한 형태를 사용하여, 전자와 양전자가 각각 독립적으로 움직이는 상황을 구현하였다. 초기 조건은 Harris‑like 전류 시트에 가이드 필드 B_G 를 추가한 형태이며, 전류 시트 두께 L 에 비해 전자 관성 길이와 라모어 반경을 충분히 작게 설정해 물리적 스케일을 확보하였다. 시뮬레이션은 Lax‑Wendroff 스킴에 인공 점성을 도입해 수치적 안정성을 유지하면서도 실제 물리적 점성 효과를 모사하도록 설계되었다.

가이드 필드가 없는 경우(Run 1)와 비교했을 때, B_G/B_0 = 0.5 이상인 경우에는 재결합 영역 주변에 강한 전하 분리(양전자는 상부, 전자는 하부)와 B_y 압축이 관찰된다. 전하 분리 정도는 ρ_c/(γ_p n_p+γ_e n_e)≈0.5에 달해, 비중성 층이 넓은 영역에 걸쳐 존재한다는 점이 특징이다. 이는 기존의 비상대론적 혹은 단일 유체 MHD 모델에서는 거의 무시되는 현상으로, 두 유체 간의 상대 속도 차이가 전기장을 유도하고, 결과적으로 E×B 흐름을 재결합 구역으로 끌어들인다.

재결합 전기장 E_y 의 구성 요소를 Ohm’s law 형태로 분해한 결과, 대류 항(−hv×B)와 관성 항(특히 advective inertia)이 주된 기여를 한다는 것이 확인되었다. 인위적인 마찰 저항(τ_fr)도 존재하지만, 로렌츠 인자 γ≈1.6 정도가 이를 억제하여 전체 E_y 의 약 1/4~1/3 정도만을 담당한다. 중앙 확산층에서는 인공 점성 항이 advective inertia 를 대체하며, 이는 실제 플라즈마에서의 미세한 입자 충돌이나 파동‑입자 상호작용을 근사적으로 표현한 것으로 해석될 수 있다.

에너지 흐름을 Poynting 플럭스와 플라즈마 엔탈피·운동 에너지로 분해한 결과, 가이드 필드가 강해질수록 B_y 에 의한 Poynting 플럭스 비중이 급격히 증가한다. σ_y,in (가이드 필드의 마그네틱 에너지 비율)이 0.1 수준에서 0.07 정도의 재결합 속도 R을 보이며, σ_y,in≈0.05 일 때는 Poynting 플럭스가 전체 에너지 전달의 절반 이상을 차지한다. 반면, 가이드 필드가 없는 경우에는 엔탈피 플럭스가 주된 에너지 운반 메커니즘이며, 전자·양전자 쌍의 가열이 주요 에너지 변환 경로가 된다.

이러한 결과는 상대론적 RMHD 모델이 전하 분리와 관성 효과를 충분히 포착하지 못함을 시사한다. 두 유체 모델은 전하 비중성 가정이 깨지는 상황에서도 안정적인 재결합을 유지할 수 있음을 보여주며, 특히 고마그네틱화(σ≫1) 환경에서 가이드 필드가 에너지 변환 스펙트럼을 어떻게 재편성하는지를 명확히 밝힌다.