차세대 친환경 기저 전력원인 소형 모듈 원자로(SMR, Small Modular Reactor)는 외부 전원이 전면 차단되는 완전 정전(SBO) 비상사태에서도 인위적인 전력 공급이나 냉각 펌프 가동 없이 오직 자연적인 중력과 온도 차이에 의한 부력만으로 노심 잔열을 냉각조(Emergency Coolant Pool)로 영구 배출하는 수동 잔열 제거 계통(PRACS)을 장착합니다. 본 백서에서는 2022 개정 물리·화학·수학·정보를 융합하여, 1차원 폐루프 운동량 보존과 웨이-위그너(Way-Wigner) 방사성 붕괴열 감쇠 곡선을 결합한 1,000Hz 무할당 C# 온디바이스 열수력 시뮬레이터를 공개합니다.

1. 2022 개정 교육과정 교과 연계 분석

  • 물리학 (유체 열팽창, 밀도 차이와 아르키메데스 부력, 마찰 압력 손실): 온도에 따른 냉각재 부피 팽창과 밀도 변화 $\Delta \rho = -\rho_0 \beta \Delta T$, 고온 상승부와 저온 하강부 높이차 $H$에 의한 자연 순환 구동 압두 $\Delta P_{driving} = \rho_0 g \beta (T_{hot} - T_{cold}) H$, 관마찰 손실 $\Delta P_{loss} = K_{loss} \frac{\dot{m}^2}{2 \rho A^2}$.
  • 화학/원자력 (원자핵 반응, 핵분열 생성물 방사성 붕괴열 시간 감쇠): 원자로 비상 정지(SCRAM) 후 핵분열 생성물의 $\beta$선, $\gamma$선 방출 붕괴열 모델 $P_{decay}(t) = 0.066 P_0 (t^{-0.2} - (t + T_{op})^{-0.2})$.
  • 수학 (폐루프 1차원 질량유량 과도응답 상미분방정식 수치해석): 유체 질량유량 $\dot{m}(t)$의 시간 변화율 $\frac{d\dot{m}}{dt} = \frac{A}{L} \left( \Delta P_{driving} - \Delta P_{loss} \right)$에 대한 1차 오일러 및 수렴 적분.
  • 정보 (무할당 실시간 1,000Hz 원자로 열수력 안전 해석기): 1ms 실시간 루프 내에서 가비지 컬렉터 없이 스택 메모리 구조체 `ref struct`로 노심 최고 온도와 피동 냉각 루프 유량을 0-Alloc 모니터링하는 C# 엔진.

2. SMR 수동 자연 순환 열수력 수식 모델

평균 유동 단면적 $A$, 루프 총 길이 $L$인 1차원 폐루프에서 자연 대류 질량유량 $\dot{m}$의 지배방정식은 다음과 같습니다.

$$\frac{L}{A} \frac{d\dot{m}}{dt} = \rho_0 g \beta \Delta T H - \frac{1}{2\rho_0 A^2} \left( \sum K_i + f \frac{L}{D} \right) \dot{m}^2$$

여기서 노심 출구 온도 $T_{hot}$과 열교환기 출구 온도 $T_{cold}$는 노심 발열량 $Q_{decay}(t)$와 유체 열용량 $C_p$의 열평형에 의해 결정됩니다.

$$T_{hot} = T_{cold} + \frac{Q_{decay}(t)}{\dot{m} C_p}$$

펌프가 멈추더라도 노심 온도가 상승하면 부력 $\Delta P_{driving}$이 자동으로 증가하여 순환 유량이 스스로 회복되는 고유 안전성(Inherent Safety)을 수학적으로 증명할 수 있습니다.

3. 1,000Hz 0-Alloc C# SMR 자연 순환 시뮬레이션 코어

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Education.StemNuclear.SmrThermalHydraulics
{
    public readonly record struct SmrReactorState(
        float ThermalPowerMw,          // 원자로 초기 정격 열출력 (MW)
        float TimeAfterScramSeconds,   // 정지 후 경과 시간 (초)
        float EmergencyPoolTempCelsius // 비상 냉각수조 온도 (°C)
    );

    public ref struct SmrCoolingTelemetry
    {
        public float DecayHeatKw;
        public float NaturalCirculationMassFlowKgPerS;
        public float CoreInletTempCelsius;
        public float CoreOutletTempCelsius;
        public float PeakFuelCladdingTempCelsius;
        public bool PassiveSafetyMarginMaintained;
    }

    public static class SmrPassiveCoolingCore
    {
        private const float GravityG = 9.80665f;
        private const float WaterDensityKgM3 = 750f; // 300°C 고압 경수 밀도
        private const float ThermalExpansionBeta = 1.8e-3f; // 1/K
        private const float SpecificHeatCp = 5200f; // J/(kg*K)
        private const float ThermalCenterHeightDiffH = 4.5f; // 노심-열교환기 수직 고도차 (m)
        private const float LoopFlowAreaA = 0.25f; // 유로 단면적 (m^2)
        private const float LoopLengthL = 18.0f; // 루프 총 길이 (m)
        private const float HydraulicResistanceK = 35.0f; // 마찰 및 형상 손실 계수
        private const float CladdingSafeLimitCelsius = 650.0f;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void SimulateStep(
            in SmrReactorState state,
            ref float currentFlowRateKgS, // 스택 관리 유량 상태
            float dtSec,
            ref SmrCoolingTelemetry telemetry)
        {
            // 1. Way-Wigner 핵분열 생성물 붕괴열 모델 산출
            float t = MathF.Max(1.0f, state.TimeAfterScramSeconds);
            float decayFraction = 0.066f * (MathF.Pow(t, -0.2f) - MathF.Pow(t + 1e7f, -0.2f));
            float qDecayWatts = (state.ThermalPowerMw * 1e6f) * decayFraction;
            telemetry.DecayHeatKw = qDecayWatts / 1000f;

            // 2. 코어 입/출구 온도 예측
            float flow = MathF.Max(currentFlowRateKgS, 0.5f);
            float deltaT = qDecayWatts / (flow * SpecificHeatCp);
            float tInlet = state.EmergencyPoolTempCelsius + 15f;
            float tOutlet = tInlet + deltaT;

            telemetry.CoreInletTempCelsius = tInlet;
            telemetry.CoreOutletTempCelsius = tOutlet;
            telemetry.PeakFuelCladdingTempCelsius = tOutlet + deltaT * 0.4f;

            // 3. 열 부력 구동 압두 및 루프 저항 계산
            float drivingHead = WaterDensityKgM3 * GravityG * ThermalExpansionBeta * deltaT * ThermalCenterHeightDiffH;
            float flowLoss = 0.5f / (WaterDensityKgM3 * LoopFlowAreaA * LoopFlowAreaA) * HydraulicResistanceK * flow * flow;

            // 4. 질량유량 미분 적분 (0-Alloc)
            float dFlowDt = (LoopFlowAreaA / LoopLengthL) * (drivingHead - flowLoss);
            currentFlowRateKgS = MathF.Max(0.1f, currentFlowRateKgS + dFlowDt * dtSec);

            telemetry.NaturalCirculationMassFlowKgPerS = currentFlowRateKgS;
            telemetry.PassiveSafetyMarginMaintained = telemetry.PeakFuelCladdingTempCelsius < CladdingSafeLimitCelsius;
        }
    }
}

4. 교육적 가치 및 미래 원자력 안전 연계

본 교구는 후쿠시마 원전 사고의 교훈인 전원 상실(SBO) 위험을 어떻게 혁신적인 수동 안전 시스템(PRACS)과 자연 부력 열수력학으로 원천 극복하는지를 물리학의 열대류 및 베르누이 방정식과 직접 연계하여 가르칩니다. 학생들은 비상 정지 후 시간의 흐름에 따라 붕괴열이 감쇠하는 동안 노심과 냉각조 사이의 자연 대류가 스스로 평형 유량을 찾아가는 과정을 온디바이스 그래프로 확인하며 첨단 에너지 공학의 안전 철학을 체득합니다.