[논문 리뷰] Identification and characterization of current sheets in collisionless plasma turbulence

Amirhassan Chatraee Azizabadi, Neeraj Jain외 1명|arXiv (Cornell University)|2020. 09. 08.

Ionosphere and magnetosphere dynamics참고 문헌 68인용 수 13

한 줄 요약

이 연구는 Zhdankin 등(2013)의 MHD 기반 전류판 식별 알고리즘을 복사되지 않는 플라즈마 난류의 2D 하이브리드-키네틱 시뮬레이션에 확장하여, 노이즈가 있는 키네틱 환경에서의 전류판 자동 탐지 및 특성 분석을 가능하게 한다. 주요 발견은 전류판이 이온 관성 길이 이하에서 격자 해상도 수준으로 두껍게 얇아지며, 갈라짐 모드 및 전자 스케일 물리에 대한 잠재적 불안정성의 징후를 보이며, 그 역학적 해석을 위한 3D 시뮬레이션에서 전자 관성 효과가 필요하다는 점이다.

ABSTRACT

The properties of current sheets forming in a ion-kinetically turbulent collisionless plasma are investigated by utilizing the results of two-dimensional hybrid-kinetic numerical simulations. For this sake the algorithm proposed by Zhdankin et al. (2013) for the analysis of current sheets forming in MHD-turbulent plasmas, was extended to analyse the role and propertes of current sheets formating in a much noisier kinetically turbulent plasma. The applicability of this approach to the analysis of kinetically-turbulent plasmas is verified. Invesigated are, e.g., the effects of the choice of parameters on the current sheet recognition, viz. the threshold current density, the minimum current density and of the local regions around current density peaks. The main current sheet properties are derived, their peak current density, the peak current carrier velocity (mainly electrons), the thickness and length of the current sheets, i.e. also their aspect ratio (length/thickness). By varying the grid resolution of the simulations it is shown that, as long as the electron inertia is not taken into account, the current sheets thin down well below ion inertial length scale until numerical (grid-resolution based) dissipation stops any the further thinning.

연구 동기 및 목표

• 노이즈가 있는 키네틱 스케일 플라즈마 난류에서 전류판을 자동으로 식별하기 위한 알고리즘 개발 및 검증.
• 하이브리드-키네틱 시뮬레이션에서 전류판의 피크 전류 밀도, 두께, 길이, 길이 대비 두께 비율 등의 특성 특성화.
• 알고리즘 파라미터(기준 전류 밀도, 국소 영역 크기, 최소 전류 밀도)가 식별 신뢰성에 미치는 영향 평가.
• 수치 해상도가 전류판 얇아짐에 미치는 영향과 전자 스케일 물리가 불안정성 발생에 미치는 역할 평가.

제안 방법

• 복사되지 않는 플라즈마 난류의 2D 하이브리드-키네틱 시뮬레이션에 Zhdankin 등(2013)의 MHD 전류판 식별 알고리즘을 적용.
• 정의된 국소 영역 크기 n 내에서 기준 전류 밀도 Jthr 이상의 국소 전류 밀도 최대값을 탐지하기 위해 파이썬 기반 코드 구현.
• 피크에 연결된 연속성을 확보하기 위해 최소 전류 밀도 기준 Jmin,i = Jmax,i / 2를 사용하여 전류판 경계 정의.
• 전류 밀도의 반최대값 너비를 측정하기 위해 헤시안 행렬의 주된 고유벡터 방향을 이용해 전류판 두께 계산.
• 전류판의 2차원 횡단면 상에서 임의의 두 점 간 최대 거리로 전류판 길이 정의.
• 수치 수렴성과 전류판 얇아짐의 한계를 테스트하기 위해 시뮬레이션 격자 해상도와 입자 수 다양화.

실험 결과

연구 질문

• RQ1 기준 전류 밀도 Jthr, 국소 영역 크기 n, 최소 전류 밀도 Jmin,i 등의 알고리즘 파라미터 선택이 키네틱 난류에서 전류판 식별에 미치는 영향은 어떻게 되는가?
• RQ2 2D 하이브리드-키네틱 플라즈마 난류에서 형성된 전류판의 통계적 특성(두께, 길이, 길이 대비 두께 비율)은 무엇인가?
• RQ3 전자 관성 없이도 전류판이 이온 관성 길이 스케일 이하로 얇아지는 정도는 어느 정도인가?
• RQ4 전류판 얇아짐이 불안정성으로 이어질 수 있으며, 전자 스케일 물리와 전자-이온 유동의 역할은 무엇인가?
• RQ5 수치 해상도와 격자 기반 소산이 관측된 전류판 형태에 어떤 영향을 미치는가?

주요 결과

• 전자 관성이 포함되지 않은 2D 하이브리드-키네틱 시뮬레이션에서 전류판은 이온 관성 길이 이하로 격자 해상도 한계까지 얇아진다.
• 전자 스케일 물리가 없음에도 불구하고, 전류판은 평균 길이 대비 두께 비율 약 20를 보이며, 이온 스케일 두께에서도 유지된다.
• 알고리즘 성능은 파라미터 선택에 민감하며, 신뢰성 있는 탐지를 위해 Jthr, n, Jmin,i의 사례별 캘리브레이션 필요.
• 피크 전류 밀도와 전자 지배 전류 운반 입자의 속도는 전류판 역학의 주요 지표이다.
• 길이 대비 두께 비율 20 이상인 얇은 전류판은 갈라짐 유형 불안정성, 예를 들어 전자-유동 불안정성과 같은 유형에 취약하다.
• 이론적 추정에 따르면 전자 스케일 전류판에서 전자-이온 유동 비율 |uez/uiz| ≫ 1이 되며, 이는 불안정성의 주요 자유 에너지원을 나타낸다.