심우주 우주선 차폐 B-UHMWPE/텅스텐 다층재 감쇠 & 1,000Hz 0-Alloc C# 선량당량 모니터링 서보 엔진

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

namespace ZeliDesk.Stem.RadiationProtection
{
    /// <summary>
    /// 우주선 방사선 검출 센서 데이터 (스택 할당 0-Alloc)
    /// </summary>
    [StructLayout(LayoutKind.Sequential, Pack = 8)]
    public readonly struct SpaceRadiationSensorData
    {
        public readonly double FastNeutronFluxCm2S;  // 고속 중성자 플럭스 (/cm^2*s)
        public readonly double ThermalNeutronFluxCm2S;// 열중성자 플럭스 (/cm^2*s)
        public readonly double GammaDoseRateMicroGyH; // 감마 흡수 선량률 (uGy/h)
        public readonly double ProtonFluxCm2S;        // 고에너지 양성자 플럭스

        public SpaceRadiationSensorData(double fastN, double thN, double gamma, double proton)
        {
            FastNeutronFluxCm2S = fastN;
            ThermalNeutronFluxCm2S = thN;
            GammaDoseRateMicroGyH = gamma;
            ProtonFluxCm2S = proton;
        }
    }

    /// <summary>
    /// 선량률 상태 및 능동 차폐 보조 제어 명령
    /// </summary>
    public ref struct RadiationShieldCommand
    {
        public double ActiveShieldMotorPositionMm; // 텅스텐/B-PE 셔터 전진 변위 (mm)
        public double AuxiliaryMagneticCoilCurrentA;
        public double TotalEquivalentDoseRateMicroSvH;
        public bool RadiationStormAlert;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public RadiationShieldCommand(double posMm, double coilA, double doseSvH, bool alert)
        {
            ActiveShieldMotorPositionMm = posMm;
            AuxiliaryMagneticCoilCurrentA = coilA;
            TotalEquivalentDoseRateMicroSvH = doseSvH;
            RadiationStormAlert = alert;
        }
    }

    /// <summary>
    /// 1,000Hz 무할당 우주 방사선 복합 차폐 서보 엔진
    /// </summary>
    public ref struct SpaceRadiationShieldEngine
    {
        private const double FluxToDoseFastNeutron = 3.5e-4;   // 플럭스 -> uSv/h 변환 인자
        private const double FluxToDoseThermalNeutron = 4.2e-5;
        private const double FluxToDoseProton = 2.1e-4;
        private const double StormThresholdMicroSvH = 500.0;   // SPE 폭풍 경보 임계치 (500 uSv/h)

        private double _accumulatedDoseMicroSv;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void ExecuteMonitoringCycle(
            in SpaceRadiationSensorData sensor,
            double dtSeconds,
            out RadiationShieldCommand command)
        {
            // 1. 방사선 성분별 유효 선량당량률 계산 (uSv/h)
            double neutronDoseRate = sensor.FastNeutronFluxCm2S * FluxToDoseFastNeutron +
                                     sensor.ThermalNeutronFluxCm2S * FluxToDoseThermalNeutron;
            double gammaDoseRate = sensor.GammaDoseRateMicroGyH * 1.0; // 감마선 wR = 1
            double protonDoseRate = sensor.ProtonFluxCm2S * FluxToDoseProton;

            double totalDoseRate = neutronDoseRate + gammaDoseRate + protonDoseRate;

            // 2. 누적 피폭 선량 계산
            _accumulatedDoseMicroSv += totalDoseRate * (dtSeconds / 3600.0);

            // 3. 방사선 폭풍(SPE) 경보 및 능동 차폐기 구동 명령 산출
            bool alert = totalDoseRate > StormThresholdMicroSvH;

            double shutterPosMm = 0.0;
            double coilCurrentA = 0.0;

            if (alert)
            {
                // 위험도에 비례하여 보조 차폐벽 전진 및 초전도 편향 자기장 코일 여자
                double severity = (totalDoseRate - StormThresholdMicroSvH) / StormThresholdMicroSvH;
                shutterPosMm = Math.Clamp(50.0 + severity * 50.0, 50.0, 150.0);
                coilCurrentA = Math.Clamp(20.0 + severity * 40.0, 20.0, 120.0);
            }
            else
            {
                shutterPosMm = 0.0;
                coilCurrentA = 0.0;
            }

            command = new RadiationShieldCommand(shutterPosMm, coilCurrentA, totalDoseRate, alert);
        }
    }
}