[논문 리뷰] Advanced modeling for the HIT-SI Experiment

Alan A. Kaptanoglu, Thomas Benedett외 3명|arXiv (Cornell University)|2020. 03. 01.

Magnetic confinement fusion research참고 문헌 51인용 수 10

한 줄 요약

이 논문은 HIT-SI 스페로메이크 실험의 시뮬레이션 정확도를 향상시키기 위해 PSI-Tet 및 NIMROD 코드에 이중 온도 자기유체역학(MHD) 모델을 도입한다. 전자와 이온 온도를 별도로 진화시키는 모델은 실험 데이터와의 일치도를 향상시켜, 특히 토로이드 전류, 밀도 프로파일 및 전류 중심 대칭성에서 개선된 결과를 보이며, 높은 주입기 주파수에서의 큰 점성 및 압축성 가열을 드러낸다. 주요 기여는 저밀도, 고성능 영역에서 예측 능력 향상으로, 성능 향상은 ≈5이다.

ABSTRACT

A two-temperature magnetohydrodynamic (MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet code and used to model plasma dynamics in the HIT-SI experiment. When compared with single-temperature Hall-MHD, the two-temperature Hall-MHD model demonstrates improved qualitative agreement with experimental measurements, including: far-infrared interferometry, ion Doppler spectroscopy, Thomson scattering, and magnetic probe measurements. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5-68.5 kHz range. At all frequencies the two-temperature models result in increased toroidal current, lower chord-averaged density, and symmetrization of the current centroid, relative to single-temperature simulations. Both codes produce higher average temperatures and toroidal currents as the injector frequency is increased. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and indicate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for the artificial diffusivity, and Dirichlet wall temperature and density. Artificial diffusivity and the density boundary condition both significantly modify the plasma density profiles, leading to larger average temperatures, higher toroidal current, and increased relative density fluctuations at low diffusivity and low wall density. High power, low density simulations at 14.5 kHz achieve sufficiently high gain (G = 5) to generate significant volumes of closed flux lasting 1-2 injector periods.

연구 동기 및 목표

• 전자와 이온 온도의 별도 진화를 통합함으로써 HIT-SI 스페로메이크 실험의 시뮬레이션 정밀도를 향상시키기.
• 에너지 분할과 플라즈마 역학을 포괄하지 못하는 단일 온도 홀-MHD 모델의 한계를 해결하기.
• 간섭계측, 도플러 스펙트로스코피, 토머슨 산란 및 자기 측정계를 포함한 다수의 실험 진단 결과와의 비교를 통해 이중 온도 모델을 검증하기.
• 주입기 주파수, 인공 확산계수, 벽 경계 조건이 플라즈마 성능과 안정성에 미치는 영향을 조사하기.
• 주입기의 회로 모델 통합 및 열 운반 폐쇄 모델 개선을 통해 향후 검증 가능성을 높이기.

제안 방법

• PSI-Tet 코드에 전자 및 이온에 대해 별도 에너지 방정식을 푸는 이중 온도 홀-MHD 모델을 구현한다.
• NIMROD 코드를 사용해 상호 검증을 수행하며, 동일한 이중 온도 수식을 적용하여 코드 간 일관성을 평가한다.
• 주입기 주파수(14.5–68.5 kHz), 인공 확산계수(50–1000 m²/s²), 벽 온도 및 밀도에 대한 디리클레 경계 조건을 포함한 파라미터 스캔을 수행한다.
• 점성 및 응축성 가열을 정량화하기 위해 에너지 균형과 열역학적 유량을 계산하며, 특히 고주파수에서의 영향을 집중 분석한다.
• NIM로드에서 시뮬레이션 안정화를 위해 초확산성(Dh∇⁴n)을 적용하고, PSI-Tet의 명시적 확산성과 비교한다.
• PSI-Tet에서 실험적 파형과 주입기 기하구조를 사용하여 전류 중심 비대칭성과 파형 반응에 미치는 영향을 평가한다.

실험 결과

연구 질문

• RQ1 전자와 이온 온도의 별도 진화가 HIT-SI 시뮬레이션에서 실험 측정 결과와의 일치도에 어떤 영향을 미치는가?
• RQ2 주입기 주파수가 이중 온도 MHD 모델에서 토로이드 전류, 온도 및 밀도 프로파일에 미치는 영향은 무엇인가?
• RQ3 인공 확산계수와 벽 경계 조건(온도 및 밀도)이 플라즈마 밀도 프로파일, 온도 및 전류 중심 위치에 미치는 영향은 무엇인가?
• RQ4 14.5 kHz 주입기 주파수에서 저밀도, 고출력 시뮬레이션은 성능 향상 G ≈ 5를 달성하고 1–2 주입기 주기 동안 폐쇄된 플럭스를 유지할 수 있는가?
• RQ5 PSI-Tet는 주입기 유도 전류 중심 비대칭성을 재현하지만 NIMROD는 그렇지 않은 이유는 무엇이며, 주입기 기하구조는 이 현상에 어떤 영향을 미치는가?

주요 결과

• 이중 온도 모델은 단일 온도 모델에 비해 더 높은 체적 평균 이온 온도와 더 낮은 연직 평균 밀도를 생성하며, 높은 주입기 주파수에서 ⟨Ti⟩ > ⟨Te⟩를 보인다.
• PSI-Tet와 NIMROD 시뮬레이션 모두 이중 온도 모델링으로 토로이드 전류가 증가하고 전류 중심이 대칭화되며, 특히 주입기 기하구조를 해상도로 표현한 PSI-Tet에서 두드러진다.
• 14.5 kHz에서 저밀도, 고출력 시뮬레이션은 성능 향상 G ≈ 5를 달성하고 50–100 µs(1–2 주입기 주기) 동안 폐쇄된 플럭스를 유지하며 고성능 영역으로의 전이를 시사한다.
• PSI-Tet에서 인공 확산계수를 1000에서 50 m²/s²로 감소시키면 연직 평균 밀도가 약 20% 감소하고, 평균 온도와 토로이드 전류가 증가하며 상대적 밀도 변동이 증폭된다.
• 벽 온도 스캔 결과, 더 강한 열압력이 벽 근처에 존재함에 따라 전류 중심이 내측으로 이동하고 ⟨β⟩가 감소하며, 이는 압력 균형 효과와 일치한다.
• 이중 온도 모델은 주입기 임피던스와 체적 평균 온도가 주입기 주파수에 대해 초선형적으로 증가하는 경향을 드러내며, 고주파수에서 에너지 결합 효율이 향상됨을 시사한다.