극저온 터보펌프 캐비테이션 인듀서 역류 와류 방지 서보 엔진

1. 유체역학 및 열역학적 배경: 유효 흡입 수두(NPSHa)와 열역학적 효과(B-factor)

터보펌프 입구에서의 유효 흡입 수두(Net Positive Suction Head, NPSHa)는 정압, 증기압, 동압 수두로 정의됩니다:

극저온 유체(액체산소, 액체수소)는 상온 물과 달리 기화 시 막대한 잠열을 흡수하여 기포 주위 액체 온도를 급격히 낮추는 열역학적 효과(Thermodynamic Effect, $B$-factor)를 나타냅니다:

이 열적 지연 효과로 인해 캐비테이션 기포의 붕괴 주기가 불규칙해지며, 특정 유량비($Q/Q_d < 0.85$)에서 날개 전방으로 고속 선회 역류 와류가 뿜어져 나와 터보펌프 축계를 파괴하는 저주파 캐비테이션 불안정성을 유발합니다.

2. 공학적 모델링: 레일리-플레셋 기포 진동과 고주파 맥동 억제

단일 기포의 반경 변화율은 레일리-플레셋(Rayleigh-Plesset) 비선형 미분방정식을 따릅니다:

고주파 동압 센서 어레이로부터 $100 \sim 500\text{ Hz}$ 대역의 캐비테이션 비대칭 회전 섭동 에너지를 0.1ms마다 FFT 없이 직접 시간 도메인 적분하여, 산화제 탱크의 헬륨 가압 압력을 실시간 증압하고 인듀서 하우징 경계층 흡입 밸브를 개방합니다.

3. 10,000Hz 0-Alloc C# 극저온 인듀서 캐비테이션 서보

펌프 입구 압력 맥동을 매 0.1ms마다 수집하여 탱크 가압 밸브 및 바이패스 인젝터를 제어하는 C# 무할당 알고리즘입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.RocketTurbopump
{
    public readonly ref struct TurbopumpTelemetryFrame
    {
        public readonly double TimestampSec;
        public readonly double InletStaticPressureBar;   // 인듀서 입구 정압 (bar)
        public readonly double LoxSaturationVaporPressureBar; // 액체산소 포화 증기압 (bar)
        public readonly double HighFreqPressurePulsationRms; // 고주파 압력 맥동 실효값 (bar)
        public readonly double InducerRotationalSpeedRpm;    // 터보펌프 회전수 (RPM)
        public readonly double VolumetricFlowRateM3Sec;      // 산화제 체적 유량 (m^3/s)

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public TurbopumpTelemetryFrame(
            double timestampSec,
            double inletStaticPressureBar,
            double loxSaturationVaporPressureBar,
            double highFreqPressurePulsationRms,
            double inducerRotationalSpeedRpm,
            double volumetricFlowRateM3Sec)
        {
            TimestampSec = timestampSec;
            InletStaticPressureBar = inletStaticPressureBar;
            LoxSaturationVaporPressureBar = loxSaturationVaporPressureBar;
            HighFreqPressurePulsationRms = highFreqPressurePulsationRms;
            InducerRotationalSpeedRpm = inducerRotationalSpeedRpm;
            VolumetricFlowRateM3Sec = volumetricFlowRateM3Sec;
        }
    }

    public ref struct CavitationMitigationCommand
    {
        public double HeliumPressurizationValveRatio;  // 산화제 탱크 헬륨 가압 밸브 개도율 (0.0 ~ 1.0)
        public double InducerBleedBypassRatio;         // 역류 와류 억제 바이패스 밸브 개도율 (0.0 ~ 1.0)
        public double EstimatedCavitationNumberSigma;  // 현재 캐비테이션 수 sigma
        public bool CavitationSurgeAlert;              // 캐비테이션 서지 위험 경보 플래그
        public double TurbopumpHealthMargin;           // 펌프 수명 안전 마진

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void Reset()
        {
            HeliumPressurizationValveRatio = 0.0;
            InducerBleedBypassRatio = 0.0;
            EstimatedCavitationNumberSigma = 1.0;
            CavitationSurgeAlert = false;
            TurbopumpHealthMargin = 1.0;
        }
    }

    public sealed class CryogenicInducerServo
    {
        private const double CriticalSigmaThreshold = 0.042; // 캐비테이션 서지 발생 임계 시그마
        private const double PulsationRmsThresholdBar = 0.35; // 허용 압력 맥동 한계 (bar)
        private const double LoxDensityKgM3 = 1141.0;         // 액체산소 밀도 (kg/m^3)

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void ComputeMitigationCycle(
            in TurbopumpTelemetryFrame frame,
            ref CavitationMitigationCommand command)
        {
            command.Reset();

            // 1. 블레이드 팁 속도 및 캐비테이션 수 sigma 계산
            double tipSpeedMs = (frame.InducerRotationalSpeedRpm * 2.0 * Math.PI / 60.0) * 0.075; // r=75mm
            double dynamicHead = 0.5 * LoxDensityKgM3 * tipSpeedMs * tipSpeedMs * 1.0e-5; // bar 환산
            double netPressureBar = Math.Max(0.001, frame.InletStaticPressureBar - frame.LoxSaturationVaporPressureBar);
            double sigma = netPressureBar / Math.Max(1.0, dynamicHead);
            command.EstimatedCavitationNumberSigma = sigma;

            // 2. 캐비테이션 서지 및 맥동 진동 판정
            bool surgeRisk = (sigma < CriticalSigmaThreshold) || 
                             (frame.HighFreqPressurePulsationRms > PulsationRmsThresholdBar);
            command.CavitationSurgeAlert = surgeRisk;

            // 3. 10,000Hz 헬륨 증압 및 역류 와류 능동 바이패스 제어
            if (surgeRisk)
            {
                double sigmaDeficit = Math.Max(0.0, CriticalSigmaThreshold - sigma);
                command.HeliumPressurizationValveRatio = Math.Clamp(0.5 + (sigmaDeficit * 15.0), 0.5, 1.0);
                command.InducerBleedBypassRatio = Math.Clamp(frame.HighFreqPressurePulsationRms / 0.8, 0.2, 0.85);
                command.TurbopumpHealthMargin = Math.Clamp(sigma / CriticalSigmaThreshold * 0.6, 0.1, 0.6);
            }
            else
            {
                command.HeliumPressurizationValveRatio = 0.15; // 평상시 기저 가압
                command.InducerBleedBypassRatio = 0.0;
                command.TurbopumpHealthMargin = Math.Clamp(sigma / 0.1, 0.7, 1.0);
            }
        }
    }
}

4. 물리 교육 및 로켓 추진 실무 연계 방안

물리 교사 및 항공우주 연구원은 본 시뮬레이터를 기반으로 다음과 같은 융합 수업을 진행할 수 있습니다:

• 베르누이 정리와 증기압 상변화: 고속 유동에서 정압 감소가 유발하는 비등(Boiling) 현상과 캐비테이션 기포 파괴 충격파($> 1000\text{ bar}$) 모델링
• 극저온 추진제의 열역학적 지연 효과: 액체산소 기화 잠열이 국소 온도를 냉각시켜 증기압을 낮춤으로써 기포 성장을 지연시키는 열적 안정화 탐구
• 10,000Hz 0-Alloc 발사체 온디바이스 추진 제어: 우주 발사체의 치명적 POGO 진동을 원천 봉쇄하는 초고속 C# 임베디드 펌프 제어 아키텍처 실습