Radial transport of electric current by electromagnetic microturbulence in tokamaks

기술 요약: 전자기 미세 난류에 의한 토카막 내 전기 전류의 반경 방향 수송

문제 제기
토카막 플라즈마에서 토로이달 각운동량의 난류 수송은 회전 프로파일을 결정하며, 이는 다시 가둠 성능과 안정성에 영향을 미친다. 이온의 수송 기여는 정전기적 한계 내에서 잘 이해되어 있는 반면, 전자의 기여는 작은 전자 질량으로 인해 종종 무시되곤 한다. 그러나 토로이달 각운동량의 반경 방향 플럭스는 토로이달 전류의 반경 방향 플럭스와 본질적으로 연결되어 있다. 전류 플럭스는 전하 대 질량 비에 비례하여 스케일링되므로, 겉보기에 미미해 보이는 전자 운동량 플럭스가 상당한 난류 전류 플럭스에 해당할 수 있다. 이 전류 플럭스의 발산(divergence)은 토로이달 전류 밀도를 재분배하여 안전 인자($q$) 프로파일을 수정할 가능성이 있다. 이러한 효과는 전자기적 요동(특히 마이크로티어링 모드와 키네틱 벌룬 모드)이 지배적이 되는 고베타($\beta$) 플라즈마에서 특히 중요하지만, 운동량 수송에서의 전자기적 항과 전류 재분배 사이의 연결 고리에 대한 탐구는 아직 충분히 이루어지지 않았다.

방법론
이를 조사하기 위해 저자들은 두 가지 국소 자이로키네틱 코드인 GENE와 CGYRO에 완전한 전자기 운동량 플럭스 진단법을 구현하였다. 기존의 GENE 및 CGYRO 구현에는 섭동된 전자기장($\delta A_\parallel$ 및 $\delta B_\parallel$)의 기여가 누락되어 있었다. 저자들은 입자 운동과 전자기장(구체적으로 맥스웰 응력 및 평행 전류 섭동 관련 항)의 기여를 분리하여 토로이달 각운동량 플럭스의 전체 식을 유도하였다.

이 새로운 진단법의 타당성은 엄격한 선형 및 비선형 교차 코드 벤치마크를 통해 검증되었다:

• 선형 벤치마크: 상하 대칭성을 깨뜨리기 위해 비제로 벌룬 각도를 가진 단일 선형 모드의 시뮬레이션을 수행하였다. 결과 분석 결과, 입자 플럭스와 다양한 운동량 플럭스 성분이 GENE와 CGYRO 사이에서 10% 이내(기울도가 높을수록 5% 미만으로 개선)의 일치도를 보였다.
• 비선형 벤치마크: 회전 전단(flow shear)에 의존하지 않고 운동량 플럭스를 생성하기 위해 상하 비대칭 자기 평형 상태에서 시뮬레이션을 수행하였다. 입자, 열, 운동량의 시간 평균 정상 상태 플럭스는 두 코드 간에 5% 이내로 일치하였다.

검증된 이 진단법들을 다양한 플라즈마 베타($\beta_e$) 값 범위에서의 마이크로티어링 모드(MTM) 및 키네틱 벌룬 모드(KBM) 구동 난류에 대한 비선형 시뮬레이션에 적용하였다.

주요 결과

1. MTM 구동 난류:

   • MTM이 지배하는 난류에서 총 운동량 수송은 열 수송에 비해 비효율적이다(낮은 프란틀 수, $Pr \approx 0.04$).
   • 정전기적 기대와는 반대로, 운동량 플럭스는 전자의 난류 맥스웰 응력($\Pi^{\delta A_\parallel}_{a,EM}$) 기여에 의해 지배된다.
   • 전자 운동량 플럭스는 $\beta_e \gtrsim 0.025$에서 부스트랩 전류로부터 유도된 기준 척도($\Pi_{BS}$)를 초과한다. 이는 이와 관련된 난류 전류 플럭스가 부스트랩 전류를 유지하는 충돌 과정과 경쟁할 만큼 크다는 것을 시사한다.
   • 이 결과는 운동량 플럭스가 회전 전단에 의해 구동되든 토로이달 마하 수에 의해 구동되든 상관없이 유효하다.

2. KBM 구동 난류:

   • 이온 온도 구배(ITG)에서 KBM 지배 난류로의 전이 과정에서 열 및 총 운동량 플럭스가 모두 증가하며, 프란틀 수가 현저히 커진다.
   • 총 수송은 여전히 정전기적 이온 기여에 의해 지배된다.
   • 그러나, 전이 과정 동안 전자 운동량 플럭스가 급격히 상승하며, 이는 거의 전적으로 맥스웰 응력 항에 의한 것이다.
   • KBM 지배 사례에서도 전자 운동량 플럭스는 부스트랩 기반 기준 척도($\Pi_{BS}$)를 초과한다.
   • 결정적으로, 맥스웰 응력 항을 유지하는 것이 필수적이다. 이 항이 없다면 전자 운동량 플럭스가 과소평가되어, 상당한 전류 재분배의 가능성을 놓치게 될 것이다.

의의 및 주장
본 논문은 전자기적 요동가 토카막의 난류 운동량 수송을 질적으로 변화시킨다고 주장한다. 주요 발견은 맥스웰 응력에 의해 매개되는 전자 운동량 플럭스가 부스트랩 전류와 경쟁할 만큼 상당한 수준의 난류 전류 플럭스를 구동할 수 있다는 점이다.

저자들은 다음과 같이 단언한다:

• 고베타 플라즈마(구형 토카막 또는 STEP와 같은 미래 장치)에서 전자기적 운동량 수송을 무시하는 것은 회전 프로파일뿐만 아니라 전류 및 안전 인자 프로파일의 주요 차수 진화에 대한 잘못된 예측을 초래할 수 있다.
• 이 효과는 총 토로이달 각운동량 수송은 비효율적인 상태(회전 구배를 보존함)를 유지하면서도, 동시에 상당한 토로이달 전류 재분배를 일으키는 영역을 허용한다.
• 본 연구는 이 효과의 수치적 존재를 확립하고 이를 측정하기 위한 진단법을 검증하였으며, 이것이 중요해지는 특정 파라미터 영역(고 $\beta_e$, MTM/KBM 지배 영역)을 식별하였다.

저자들은 국소적 플럭스 크기가 중요성을 시사하더라도, 실제 전류 프로파일의 수정 여부는 플럭스의 부호, 반경 방향 변화 및 전역적 발산에 달려 있으며 이는 국소 시뮬레이션만으로는 결정될 수 없음을 언급하며, 전역적 함의에 대해서는 신중한 태도를 유지하였다. 그들은 이러한 역학에 대한 물리적 이론을 개발하고 그 전역적 결과를 결정하는 것이 향후 연구 과제라고 결론지었다.