대한민국 해군 도산안창호급(KSS-III) 잠수함 리튬이온(Li-ion) 배터리 팩 오프가스·전압·온도 8채널 센서 융합 및 1,000Hz 0-Alloc C# 셀 열폭주 선제 격리 엔진

1. 리튬이온 셀 전기화학적 등가회로(ECM) 및 Bernardi 열생성 역학

잠수함 추진 모터에 대전류를 방전하는 리튬이온 셀의 전기적 거동은 2차 테브난(Thevenin) 등가회로 모델에 의해 다음과 같이 기술됩니다:

여기서 $OCV$는 개방회로전압, $R_0$는 옴 내부저항, $R_1 C_1$과 $R_2 C_2$는 전기이중층 및 확산 편극 시정수입니다. 이때 셀 내부에서 발생하는 총 열량 $\dot{Q}_{gen}$은 가역적 엔트로피 열과 비가역적 줄열을 결합한 Bernardi 열생성 공식에 의해 결정됩니다:

가역적 반응열 항목 $-IT \frac{\partial OCV}{\partial T}$는 충·방전 상태에 따라 흡열 또는 발열로 작용하지만, 고출력 기동 시에는 내부저항에 의한 줄 발열($I^2 R_0$)과 과전압 편극 손실이 지배적으로 셀 온도를 가속시킵니다.

2. 3단계 아레니우스 열폭주 반응속도론 및 오프가스 생성 메커니즘

비정상 발열 또는 내부 미세 단락(Micro Short-circuit)이 발생할 때 리튬이온 셀은 3단계의 비가역적 연쇄 열분해 반응을 거치며 열폭주로 전이됩니다:

1단계 SEI 분해 시 전해액 유기용매(EC, DMC)가 분해되면서 $H_2, CO, CH_4, C_2H_4$ 등의 가스가 발생하여 셀 내압이 상승하고 안전 벤트(Vent)가 개방됩니다. 이때 방출되는 수소($H_2$)와 일산화탄소($CO$)의 분압 증가 속도는 온도 센서의 열전도 지연(Thermal Lag)보다 수십 초 빠르게 감지되므로, 최선제 조기 경보의 결정적 지표가 됩니다.

3. 오프가스·전압·온도 8채널 다중 지표 열폭주 판별 함수

단일 센서 오작동으로 인한 오차단을 방지하기 위해 8개 채널($H_2, CO, VOC$, 셀 온도 $T$, 셀 전압 $V$, 전압 변화율 $\dot{V}$, 온도 변화율 $\dot{T}$, 모듈 전류 $I$)의 정규화 가중치 판별 지수(Thermal Runaway Risk Index, $\mathcal{R}_{TR}$)를 다음과 같이 정의합니다:

여기서 $\sigma(x) = 1/(1+e^{-x})$는 시그모이드 정규화 함수입니다. $\mathcal{R}_{TR} > 0.85$ 임계치 초과 시 즉각 해당 배터리 스트링의 고속 파이로퓨즈(Pyrofuse) 기폭 신호를 인가하여 10ms 이내에 갈바닉 절연(Galvanic Isolation)을 수행합니다.

4. 1,000Hz 0-Alloc C# 배터리 열폭주 선제 차단 엔진

도산안창호급 잠수함 배터리 관리 시스템(BMS)에 실장되는 아래의 C# 엔진은 8채널 센서 데이터를 입력받아 Bernardi 열생성 모델 및 위험 지수 $\mathcal{R}_{TR}$을 가비지 컬렉션 할당 없이(Zero-Alloc) 1,000Hz(1.0ms) 주기로 산출하는 고신뢰성 국방 아키텍처입니다.

using System;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;
using System.Numerics;

namespace ZeliDesk.Defense.SubmarineESS
{
    /// <summary>
    /// 리튬이온 배터리 모듈 8채널 다중 센서 계측 구조체
    /// </summary>
    [StructLayout(LayoutKind.Sequential)]
    public readonly struct BatterySensorFrame
    {
        public readonly float CellVoltageVolts;    // 셀 전압 [V] (정상: 3.2~4.2)
        public readonly float CellTemperatureC;    // 셀 표면 온도 [C]
        public readonly float ModuleCurrentAmps;   // 모듈 충방전 전류 [A]
        public readonly float HydrogenPpm;         // 수소(H2) 오프가스 농도 [ppm]
        public readonly float CarbonMonoxidePpm;   // 일산화탄소(CO) 농도 [ppm]
        public readonly float VocIndex;            // 휘발성 유기화합물 VOC 지수
        public readonly float AmbientPressureKpa;  // 배터리 격실 내압 [kPa]
        public readonly float CoolantFlowLpm;      // 냉각수 유량 [L/min]

        public BatterySensorFrame(float v, float temp, float current, float h2, float co, float voc, float press, float coolant)
        {
            CellVoltageVolts = v;
            CellTemperatureC = temp;
            ModuleCurrentAmps = current;
            HydrogenPpm = h2;
            CarbonMonoxidePpm = co;
            VocIndex = voc;
            AmbientPressureKpa = press;
            CoolantFlowLpm = coolant;
        }
    }

    /// <summary>
    /// 배터리 안전 제어 및 격리 지령 구조체
    /// </summary>
    [StructLayout(LayoutKind.Sequential)]
    public struct BatteryProtectionCommand
    {
        public float ThermalRunawayRisk; // 열폭주 위험 지수 R_TR (0.0 ~ 1.0)
        public float HeatGenerationWatts; // 실시간 발열량 Q_gen [W]
        public float VoltageDropRateVps; // 전압 강하율 dV/dt [V/s]
        public float TempRiseRateCps;    // 온도 상승율 dT/dt [C/s]
        public bool TriggerPyrofuseTrip; // 파이로퓨즈 고속 강제 격리 지령
        public bool EmergencyCoolingMax; // 비상 냉각 펌프 최대 가동
        public int WarningLevel;         // 위험 단계: 0=정상, 1=주의, 2=경고, 3=비상차단
    }

    /// <summary>
    /// 1,000Hz 무할당(Zero-Alloc) KSS-III 리튬이온 배터리 열폭주 선제 차단 엔진
    /// </summary>
    public ref struct Kss3BatteryProtectionEngine
    {
        private const float H2ThresholdPpm = 25.0f;     // H2 경보 임계치
        private const float COThresholdPpm = 15.0f;     // CO 경보 임계치
        private const float TempRiseRateLimitCps = 1.2f;// 초당 1.2도 상승 한계
        private const float VoltageDropLimitVps = -0.15f;// 초당 0.15V 급락 한계
        private const float InternalResistanceOhm = 0.0018f; // 셀 내부저항 1.8 mOhm

        private float _lastVoltage;
        private float _lastTemp;
        private bool _isInitialized;

        public Kss3BatteryProtectionEngine(float initialVoltage = 3.85f, float initialTemp = 25.0f)
        {
            _lastVoltage = initialVoltage;
            _lastTemp = initialTemp;
            _isInitialized = true;
        }

        /// <summary>
        /// 1,000Hz (dt = 0.001s) 주기 열폭주 징후 판별 및 보호 지령 산출 루프
        /// </summary>
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void Update(in BatterySensorFrame sensor, float dt, ref BatteryProtectionCommand cmd)
        {
            if (!_isInitialized)
            {
                _lastVoltage = sensor.CellVoltageVolts;
                _lastTemp = sensor.CellTemperatureC;
                _isInitialized = true;
            }

            // 1. 전압 및 온도 미분율 계산 (초당 변화율)
            float invDt = 1.0f / dt;
            float dV_dt = (sensor.CellVoltageVolts - _lastVoltage) * invDt;
            float dT_dt = (sensor.CellTemperatureC - _lastTemp) * invDt;
            _lastVoltage = sensor.CellVoltageVolts;
            _lastTemp = sensor.CellTemperatureC;

            cmd.VoltageDropRateVps = dV_dt;
            cmd.TempRiseRateCps = dT_dt;

            // 2. Bernardi 열생성량 Q_gen = I^2 * R0 + 엔트로피항
            float i = sensor.ModuleCurrentAmps;
            float jouleHeat = i * i * InternalResistanceOhm;
            float entropyHeat = -i * (sensor.CellTemperatureC + 273.15f) * (-0.00022f); // dOCV/dT 근사
            cmd.HeatGenerationWatts = jouleHeat + entropyHeat;

            // 3. 오프가스 및 상태 지표 시그모이드 정규화
            float sH2 = Sigmoid((sensor.HydrogenPpm - H2ThresholdPpm) * 0.15f);
            float sCO = Sigmoid((sensor.CarbonMonoxidePpm - COThresholdPpm) * 0.2f);
            float sTdot = Sigmoid((dT_dt - TempRiseRateLimitCps) * 2.5f);
            float sVdot = Sigmoid((-dV_dt - (-VoltageDropLimitVps)) * 10.0f);
            float sTempAbs = Sigmoid((sensor.CellTemperatureC - 65.0f) * 0.3f);

            // 4. 가중치 융합 위험 지수 R_TR 산출 (가스 40%, 미분율 40%, 절대온도 20%)
            float risk = (0.25f * sH2) + (0.15f * sCO) + (0.25f * sTdot) + (0.15f * sVdot) + (0.20f * sTempAbs);
            cmd.ThermalRunawayRisk = Math.Clamp(risk, 0.0f, 1.0f);

            // 5. 다단계 보호 상태 기계 (FSM)
            if (cmd.ThermalRunawayRisk > 0.80f || sensor.HydrogenPpm > 100.0f || sensor.CellTemperatureC > 85.0f)
            {
                cmd.WarningLevel = 3; // 비상 차단
                cmd.TriggerPyrofuseTrip = true;
                cmd.EmergencyCoolingMax = true;
            }
            else if (cmd.ThermalRunawayRisk > 0.55f || sensor.HydrogenPpm > 40.0f)
            {
                cmd.WarningLevel = 2; // 경고
                cmd.TriggerPyrofuseTrip = false;
                cmd.EmergencyCoolingMax = true;
            }
            else if (cmd.ThermalRunawayRisk > 0.30f || sensor.CellTemperatureC > 45.0f)
            {
                cmd.WarningLevel = 1; // 주의
                cmd.TriggerPyrofuseTrip = false;
                cmd.EmergencyCoolingMax = false;
            }
            else
            {
                cmd.WarningLevel = 0; // 정상
                cmd.TriggerPyrofuseTrip = false;
                cmd.EmergencyCoolingMax = false;
            }
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        private static float Sigmoid(float x)
        {
            if (x >= 10.0f) return 1.0f;
            if (x <= -10.0f) return 0.0f;
            return 1.0f / (1.0f + MathF.Exp(-x));
        }
    }
}

5. 2022 개정 화학·물리·정보 교육과정 연계표