대한민국 해군 정조대왕급 이지스 구축함 차세대 한국형 수직발사체계(KVLS-II) 대형 로켓 모터 점화 시 초음속 플룸 가스 팽창파 충격 및 1,000Hz 0-Alloc C# 배기 덕트 벤트 밸브 과압 해소 제어 엔진

1. KVLS-II 초음속 플룸 가스 유동 및 1차원 압축성 충격파 역학

고체 로켓 모터 노즐 출구(단면적 $A_e$)에서 분출되는 초음속 연소가스의 배출 압력 $P_e$, 마하수 $M_e \approx 2.5$, 가스 비열비 $\gamma \approx 1.25$ 하에서 플레넘 챔버로 진입하는 1차원 충격파 관통 거동은 질량, 운동량, 에너지 보존 법칙에 의해 결정됩니다:

초음속 제트가 U-Duct 바닥면에 충돌하여 형성되는 정체 압력(Stagnation Pressure $P_0$)은 레일리 초음속 피토관 공식에 의해 정압 $P_1$ 대비 10배 이상 증폭됩니다:

2. 플레넘 챔버 압력 상승 및 초크 벤트(Choked Vent) 밸브 모델

체적 $V_{pl}$의 플레넘 챔버 내부에서 로켓 가스 유입량 $\dot{m}_{in}$, U-Duct 배기 유량 $\dot{m}_{duct}$, 그리고 상부 과압 방출 밸브 배출 유량 $\dot{m}_{vent}$ 간의 에너지 수지식에 따른 압력 변화율은 다음과 같습니다:

플레넘 압력 $P_{pl}$이 대기압 $P_{atm}$ 대비 임계 압력비($P_{pl}/P_{atm} \ge 1.85$)를 초과하므로, 벤트 밸브를 통과하는 유동은 항상 음속 초크 상태(Choked Flow)를 유지합니다:

과압 허용 한계($P_{\text{limit}} = 1.2\,\text{MPa}$)를 사수하기 위해, 제어기는 실측 압력 상승률 $\frac{dP_{pl}}{dt}$에 선제 비례하여 벤트 도어 단면적 $A_{vent}(t)$를 $15\,\text{ms}$ 이내로 완전 전개합니다.

3. 1,000Hz 0-Alloc C# Kvls2GasReliefEngine.cs 구현

아래의 C# 구현체는 플레넘 챔버 4개 압력 센서, 연소가스 온도계, 발사 시퀀스 이벤트를 1ms 주기로 수신하여 가비지 컬렉션(GC) 0% 무할당으로 유압 벤트 댐퍼 개도율을 결정론적으로 지령하는 온디바이스 엔진입니다.

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

namespace ZeliDesk.Defense.NavyWeapons
{
    /// <summary>
    /// KVLS-II 플레넘 챔버 텔레메트리
    /// </summary>
    [StructLayout(LayoutKind.Sequential)]
    public struct KvlsPlenumTelemetry
    {
        public Vector4 ChamberPressuresMpa; // 4개 센서 실측 압력 [MPa] (정격 안전한계 1.2 MPa)
        public float GasTemperatureK;       // 연소가스 온도 [K]
        public float RocketThrustKn;        // 로켓 모터 실측 추력 [kN]
        public float IgnitionTimeSec;       // 점화 후 경과 시간 [s]
        public bool RocketIgnited;          // 점화 개시 플래그
    }

    /// <summary>
    /// 1,000Hz 과압 벤트 밸브 출력 지령
    /// </summary>
    [StructLayout(LayoutKind.Sequential)]
    public struct KvlsVentCommand
    {
        public float VentDamperOpeningRatio;// 벤트 밸브 개도율 [0.0 ~ 1.0]
        public float EstimatedChamberPressure;// 챔버 평균 압력 [MPa]
        public float PressureRiseRateMpaPerS; // 압력 상승률 dP/dt [MPa/s]
        public float ChokedFlowMassKgPerSec;  // 벤트 배출 초크 유량 [kg/s]
        public bool OverpressureAlert;        // 과압 위험 경보 (1.0 MPa 초과 시)
    }

    /// <summary>
    /// 1,000Hz 무할당(Zero-Alloc) KVLS-II 핫가스 과압 릴리프 엔진
    /// </summary>
    public ref struct Kvls2GasReliefEngine
    {
        private const float MaxAllowedPressureMpa = 1.20f;  // 최대 허용 압력
        private const float TargetPressureMpa = 0.85f;      // 목표 유지 압력
        private const float MaxVentAreaM2 = 0.45f;          // 벤트 최대 개구 면적
        private const float Gamma = 1.25f;                  // 고체 추진제 비열비
        private const float GasConstantR = 320.0f;          // 가스 정수 J/(kg.K)

        private float _lastMeanPressure;
        private float _filteredOpening;

        public Kvls2GasReliefEngine(float initialOpening = 0.0f)
        {
            _lastMeanPressure = 0.1013f; // 초기 대기압
            _filteredOpening = initialOpening;
        }

        /// <summary>
        /// 1,000Hz (dt = 0.001s) 초고속 핫가스 과압 방출 밸브 제어 루프
        /// </summary>
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void Update(in KvlsPlenumTelemetry telemetry, ref KvlsVentCommand cmd)
        {
            const float dt = 0.001f;

            // 1. 4개 센서 평균 압력 및 압력 상승률(dP/dt) 산출
            float meanP = (telemetry.ChamberPressuresMpa.X + telemetry.ChamberPressuresMpa.Y +
                           telemetry.ChamberPressuresMpa.Z + telemetry.ChamberPressuresMpa.W) * 0.25f;
            cmd.EstimatedChamberPressure = meanP;

            float dpDt = (meanP - _lastMeanPressure) / dt;
            cmd.PressureRiseRateMpaPerS = dpDt;
            _lastMeanPressure = meanP;

            cmd.OverpressureAlert = (meanP >= 1.0f);

            // 2. 점화 시 선제적 피드포워드 + 압력 피드백 밸브 제어
            float targetOpening = 0.0f;
            if (telemetry.RocketIgnited)
            {
                // 점화 초기 0.2초간 급격한 압력파에 대비한 선제 개방
                float feedforward = (telemetry.IgnitionTimeSec <= 0.35f) ? 0.65f : 0.25f;
                
                // 압력 오차 및 상승률 비례 게인
                float pErr = MathF.Max(0.0f, meanP - TargetPressureMpa);
                float feedback = (pErr * 1.5f) + (MathF.Max(0.0f, dpDt) * 0.04f);

                targetOpening = Math.Clamp(feedforward + feedback, 0.0f, 1.0f);
            }

            // 고속 유압 액추에이터 동특성 필터링 (슬루율 제한)
            _filteredOpening += (targetOpening - _filteredOpening) * 0.12f;
            cmd.VentDamperOpeningRatio = _filteredOpening;

            // 3. 초크 유동 배출 질량 유량 계산
            float actualArea = MaxVentAreaM2 * cmd.VentDamperOpeningRatio;
            float tGas = MathF.Max(300.0f, telemetry.GasTemperatureK);
            float sonicFactor = MathF.Sqrt(Gamma / (GasConstantR * tGas)) * 
                                MathF.Pow(2.0f / (Gamma + 1.0f), (Gamma + 1.0f) / (2.0f * (Gamma - 1.0f)));

            // P [Pa] = meanP * 1e6
            cmd.ChokedFlowMassKgPerSec = 0.92f * actualArea * (meanP * 1e6f) * sonicFactor;
        }
    }
}

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