토카막 대붕괴 완화 및 도주 전자(Runaway Electron) 파쇄 펠릿(SPI) 100㎲ 긴급 주입 서보 엔진

1. 물리적 배경: 토카막 전류 붕괴(CQ)와 도주 전자 눈사태

토카막 플라스마 안정성 한계를 넘어서면 수 밀리초 내에 플라스마 온도가 $10\text{ keV}$에서 $10\text{ eV}$ 이하로 곤두박질치는 열적 붕괴(TQ)가 발생합니다. 플라스마 저항률 $\eta \propto T_e^{-3/2}$이 급증함에 따라 렌츠의 법칙($\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$)에 의해 수십 $\text{V/m}$에 달하는 강력한 토로이달 유도 전기장 $E_\parallel$이 형성됩니다:

전기장 $E_\parallel$이 콘월-드라이서 임계 전기장 $E_c$를 초과하면 쿨롱 충돌 저항이 속도 증가에 따라 감소하는 상대론적 영역에 진입하여 도주 전자가 발생하며, 이 전자가 배경 전자를 2차 충돌 이온화시키는 코너-헤이스티(Connor-Hastie) 눈사태(Avalanche) 증폭이 발생합니다:

수 메가암페어($\sim \text{MA}$)의 고에너지(수십 $\text{MeV}$) 도주 전자 빔이 토카막 진공용기 및 베릴륨/텅스텐 대면재(PFC)에 국소 충돌하면 용기 관통 및 냉각수 폭발 사고를 유발합니다.

2. 공학적 모델링: 파쇄 펠릿(SPI) 기화 및 불순물 혼합 역학

도주 전자 빔의 형성을 차단하려면 플라스마 밀도를 로젠블루스(Rosenbluth) 밀도 $n_R$ 이상으로 급격히 증가시켜야 합니다. 극저온($\sim 5\text{ K}$) 고체 네온/중수소 펠릿을 초속 $200 \sim 800\text{ m/s}$로 가속한 뒤 완충 튜브 끝단의 쐐기(Shredder)에 충돌시켜 마이크로 파편 조각군으로 파쇄 분사합니다:

파쇄된 펠릿 조각들이 플라스마 중심부까지 침투하여 복사 냉각을 유도하고 전자 밀도를 $10^{21}\text{ m}^{-3}$ 이상으로 높여 임계 전기장 $E_c$를 $E_\parallel$ 이상으로 끌어올림으로써 눈사태를 소멸시킵니다.

3. 10,000Hz 0-Alloc C# 토카막 대붕괴 완화 SPI 긴급 서보

토로이달 자속 코일, 미르노프 코일(Mirnov Coil), 경X선(Hard X-ray) 검출 신호를 0.1ms 주기로 수집하여 도주 전자 급발작 전조를 판정하고 고압 가스 솔레노이드 밸브를 초고속 트리거하는 C# 무할당 엔진입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.FusionTokamak
{
    public readonly ref struct TokamakDisruptionFrame
    {
        public readonly double TimestampSec;
        public readonly double PlasmaCurrentMA;         // 플라스마 전류 (MA)
        public readonly double LoopVoltageVolts;        // 토로이달 루프 전압 (V)
        public readonly double HardXRayCountRate;       // 경X선 계수율 (counts/ms)
        public readonly double MirnovMagneticFluctuation; // dB/dt 자장 섭동 진폭 (T/s)
        public readonly double CoreElectronTempKeV;     // 중심 전자 온도 (keV)

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public TokamakDisruptionFrame(
            double timestampSec,
            double plasmaCurrentMA,
            double loopVoltageVolts,
            double hardXRayCountRate,
            double mirnovMagneticFluctuation,
            double coreElectronTempKeV)
        {
            TimestampSec = timestampSec;
            PlasmaCurrentMA = plasmaCurrentMA;
            LoopVoltageVolts = loopVoltageVolts;
            HardXRayCountRate = hardXRayCountRate;
            MirnovMagneticFluctuation = mirnovMagneticFluctuation;
            CoreElectronTempKeV = coreElectronTempKeV;
        }
    }

    public ref struct SpiServoCommand
    {
        public bool TriggerGasPuff;                     // 고압 헬륨/수소 발사 밸브 개방
        public double SolenoidFireDelayMicroSec;       // 발사 지연 시간 (us)
        public double PelletBarrelPressureBar;         // 가속 총열 압력 (bar)
        public double EstimatedAvalancheGrowthRate;    // 도주 전자 눈사태 성장률 (s^-1)
        public double MitigationMargin;                 // 대붕괴 억제 안전 마진

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void Reset()
        {
            TriggerGasPuff = false;
            SolenoidFireDelayMicroSec = 0.0;
            PelletBarrelPressureBar = 0.0;
            EstimatedAvalancheGrowthRate = 0.0;
            MitigationMargin = 1.0;
        }
    }

    public sealed class TokamakRunawayElectronSpiServo
    {
        private const double CriticalLoopVoltageThreshold = 18.5; // TQ 직후 임계 유도 전압 (V)
        private const double MirnovMhdInstabilityLimit = 42.0;    // MHD 섭동 한계 (T/s)
        private const double HardXRayAvalancheTrigger = 1250.0;   // HXR 도주 전자 계수율 (cps/ms)
        private const double BaseBarrelPressureBar = 60.0;        // SPI 가속 기본 가스압 (bar)

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void EvaluateDisruptionCycle(
            in TokamakDisruptionFrame frame,
            ref SpiServoCommand command)
        {
            command.Reset();

            // 1. 코어 열적 붕괴(TQ) 조기 감지: 전자 온도 급락 및 자장 섭동 급증
            bool isThermalQuenchImminent = (frame.CoreElectronTempKeV < 0.15) || 
                                           (frame.MirnovMagneticFluctuation > MirnovMhdInstabilityLimit);

            // 2. 상대론적 도주 전자 눈사태 성장률 추정
            double electricFieldRatio = Math.Max(0.0, (frame.LoopVoltageVolts / CriticalLoopVoltageThreshold) - 1.0);
            double growthRate = electricFieldRatio * 850.0 + (frame.HardXRayCountRate / 100.0);
            command.EstimatedAvalancheGrowthRate = growthRate;

            // 3. 파쇄 펠릿(SPI) 100us 긴급 발사 트리거 판정
            if (isThermalQuenchImminent || frame.HardXRayCountRate > HardXRayAvalancheTrigger || frame.LoopVoltageVolts > CriticalLoopVoltageThreshold)
            {
                command.TriggerGasPuff = true;
                // 급격한 눈사태 성장 시 솔레노이드 디스크 트리거를 50us 내 즉시 개방
                command.SolenoidFireDelayMicroSec = Math.Max(20.0, 100.0 - (growthRate * 0.05));
                // 펠릿 파쇄 도달 속도를 높이기 위해 총열 압력을 100 bar까지 승압
                command.PelletBarrelPressureBar = Math.Clamp(BaseBarrelPressureBar + (growthRate * 0.02), 60.0, 120.0);
                command.MitigationMargin = 0.0; // 비상 억제 루틴 활성화
            }
            else
            {
                command.TriggerGasPuff = false;
                command.SolenoidFireDelayMicroSec = 0.0;
                command.PelletBarrelPressureBar = BaseBarrelPressureBar;
                command.MitigationMargin = Math.Clamp(1.0 - (frame.LoopVoltageVolts / CriticalLoopVoltageThreshold), 0.0, 1.0);
            }
        }
    }
}

4. 교육 및 실무 시뮬레이션 적용 방안

초·중·고 교사 및 핵융합 연구자들은 본 시뮬레이터를 통해 다음과 같은 물리 탐구 활동을 진행할 수 있습니다:

• 패러데이 법칙과 플라스마 붕괴 전압 측정: 플라스마 전류가 수 밀리초 내에 급감할 때 유도 기전력이 진공용기와 도주 전자에 미치는 영향 분석
• 상대론적 전자 눈사태의 동역학: 고에너지 전자의 쿨롱 충돌 단면적 감소와 상대론적 질량 증가 효과 연계 탐구
• KSTAR/ITER 분산 완화 기술(DMS) 제어: 네온/중수소 펠릿 침투 깊이 및 대면재 보호를 위한 초정밀 0-Alloc 임베디드 제어 루프 검증