방사광 X선 극저온 나노프로브 금속유기골격체(MOF) 초임계 이산화탄소 분자 포집 흡착 서보 C# SIMD 엔진

3. Zero-Allocation C# SIMD MOF 흡착 제어 서보 엔진

using System;
using System.Runtime.CompilerServices;
using System.Runtime.Intrinsics;
using System.Runtime.Intrinsics.X86;

namespace ZeliDesk.Stem.Synchrotron.Mof
{
    public readonly ref struct MofAdsorptionState
    {
        public readonly float AdsorbedAmountMolKg;
        public readonly float LatticeStrainPct;
        public readonly float PiezoValveVoltage;
        public readonly bool StableAdsorption;

        public MofAdsorptionState(float amount, float strain, float voltage, bool stable)
        {
            AdsorbedAmountMolKg = amount;
            LatticeStrainPct = strain;
            PiezoValveVoltage = voltage;
            StableAdsorption = stable;
        }
    }

    public sealed class MofCarbonCaptureServoEngine
    {
        private const float SampleRateHz = 100000.0f;
        private const float Dt = 1.0f / SampleRateHz;

        // MOF 재료 물리 상수
        private const float Qmax = 8.5f;          // mol/kg
        private const float K0 = 1.2e-6f;         // 1/kPa
        private const float DeltaH = 38000.0f;    // J/mol (흡착열)
        private const float R = 8.314f;           // J/(mol*K)
        private const float N = 0.85f;            // Freundlich exponent
        private const float KappaAd = 0.0018f;    // 흡착 유도 변형 계수 (% per mol/kg)

        // PID 제어 파라미터
        private float _integralError = 0.0f;
        private float _prevError = 0.0f;
        private const float Kp = 45.0f;
        private const float Ki = 250.0f;
        private const float Kd = 0.015f;

        [MethodImpl(MethodImplOptions.AggressiveInlining | MethodImplOptions.AggressiveOptimization)]
        public unsafe MofAdsorptionState ControlAdsorptionCycle(
            float measuredBraggAngleDeg,
            float referenceBraggAngleDeg,
            float chamberPressureKpa,
            float cryogenicTempK)
        {
            // 1. 회절 각 변위 및 격자 변형률(Strain) 계산
            float deltaThetaRad = (measuredBraggAngleDeg - referenceBraggAngleDeg) * (MathF.PI / 180.0f);
            float refThetaRad = referenceBraggAngleDeg * (MathF.PI / 180.0f);

            // delta_d / d = -cot(theta) * delta_theta
            float cotTheta = 1.0f / MathF.Tan(refThetaRad);
            float latticeStrainPct = -cotTheta * deltaThetaRad * 100.0f;

            // 2. 현재 열역학 조건에서 이론 흡착량 계산 (SIMD 가속 평준화)
            float arrheniusFactor = MathF.Exp(DeltaH / (R * cryogenicTempK));
            float kP = K0 * arrheniusFactor * chamberPressureKpa;
            float kPPowN = MathF.Pow(MathF.Max(kP, 1.0e-8f), N);
            float theoreticalAdsorption = Qmax * (kPPowN / (1.0f + kPPowN));

            // 3. 목표 격자 변형률 대비 오차 산출 (목표 변형률: 0.015% 미만 유지)
            float targetStrainPct = KappaAd * theoreticalAdsorption;
            float strainError = targetStrainPct - latticeStrainPct;

            // 4. PID 연산
            _integralError += strainError * Dt;
            _integralError = Math.Clamp(_integralError, -0.5f, 0.5f);
            float derivative = (strainError - _prevError) * SampleRateHz;
            _prevError = strainError;

            float controlSignal = (Kp * strainError) + (Ki * _integralError) + (Kd * derivative);

            // 피에조 주입 밸브 제어 전압 (0 ~ 100V 범위)
            float piezoVoltage = Math.Clamp(50.0f + controlSignal, 0.0f, 100.0f);
            bool isStable = MathF.Abs(strainError) < 0.002f;

            return new MofAdsorptionState(
                amount: theoreticalAdsorption,
                strain: latticeStrainPct,
                voltage: piezoVoltage,
                stable: isStable
            );
        }
    }
}