AIDT 온디바이스 STEM 방사광 X선 오페란도(Operando) 나노촉매 미세유체 반응기 차압 및 가스 크로마토그래피 동기화 서보 제어기

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.Synchrotron.OperandoCatalysis
{
    /// <summary>
    /// 오페란도 미세유체 셀 차압 및 X선 분광 계측 상태
    /// </summary>
    public readonly ref struct OperandoCellTelemetryState
    {
        public readonly double DifferentialPressureBar; // 입출구 차압 Delta P [bar]
        public readonly double InletPressureBar;          // 인입 가스 절대압 [bar]
        public readonly double WhiteLinePeakEv;          // XANES 백색선 피크 에너지 [eV]
        public readonly double CellTemperatureCelsius;   // 반응기 온도 [degC]

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public OperandoCellTelemetryState(double diffPress, double inletPress, double peakEv, double temp)
        {
            DifferentialPressureBar = diffPress;
            InletPressureBar = inletPress;
            WhiteLinePeakEv = peakEv;
            CellTemperatureCelsius = temp;
        }
    }

    /// <summary>
    /// 압전 배압 밸브 구동 및 촉매 상태 출력
    /// </summary>
    public readonly ref struct OperandoPressureControlOutput
    {
        public readonly double PiezoValveCommandVolt; // 배압 밸브 구동 전압 [V]
        public readonly double EstimatedOxidationState;// 촉매 금속 유효 산화수 (+1.0 ~ +4.0)
        public readonly double MembraneStressRatio;   // 멤브레인 허용 응력 대비 사용률 (0.0~1.0)
        public readonly bool IsPressureSafe;          // 차압 정상(<0.10 bar) 유지 여부

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public OperandoPressureControlOutput(double volt, double oxState, double stressRatio, bool safe)
        {
            PiezoValveCommandVolt = volt;
            EstimatedOxidationState = oxState;
            MembraneStressRatio = stressRatio;
            IsPressureSafe = safe;
        }
    }

    /// <summary>
    /// 10,000Hz 오페란도 반응기 차압 서보 및 분광 동기화 엔진
    /// </summary>
    public ref struct OperandoFlowCellServoEngine
    {
        private const double MaxSafeDiffPressBar = 0.10; // 0.10 bar 초과 시 멤브레인 위험
        private readonly double _referenceEdgeEv;
        private readonly double _oxStateSensitivity;
        private readonly double _kp;
        private readonly double _ki;
        private readonly double _kd;
        private readonly double _valveGainVoltPerBar;
        private readonly double _maxValveVolt;

        private double _integral;
        private double _prevErrorBar;

        public OperandoFlowCellServoEngine(
            double referenceEdgeEv,
            double oxStateSensitivity,
            double kp,
            double ki,
            double kd,
            double valveGainVoltPerBar,
            double maxValveVolt)
        {
            _referenceEdgeEv = referenceEdgeEv;
            _oxStateSensitivity = oxStateSensitivity;
            _kp = kp;
            _ki = ki;
            _kd = kd;
            _valveGainVoltPerBar = valveGainVoltPerBar;
            _maxValveVolt = maxValveVolt;
            _integral = 0.0;
            _prevErrorBar = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public OperandoPressureControlOutput Update(
            in OperandoCellTelemetryState state,
            double targetDiffPressBar,
            double dtSeconds)
        {
            // 1. 차압 오차 연산 (목표 차압 통상 0.02 bar)
            double errorBar = state.DifferentialPressureBar - targetDiffPressBar;

            // 2. PID 제어
            _integral += errorBar * dtSeconds;
            double derivative = (errorBar - _prevErrorBar) / dtSeconds;
            _prevErrorBar = errorBar;

            double compBar = _kp * errorBar + _ki * _integral + _kd * derivative;
            double valveVolt = Math.Clamp(compBar * _valveGainVoltPerBar, 0.0, _maxValveVolt);

            // 3. SiN 멤브레인 응력 사용률 추정
            double stressRatio = Math.Clamp(Math.Abs(state.DifferentialPressureBar) / MaxSafeDiffPressBar, 0.0, 1.5);
            bool isSafe = stressRatio <= 1.0;

            // 4. XANES 피크 편위 기반 촉매 산화수 실시간 추정
            double dEv = state.WhiteLinePeakEv - _referenceEdgeEv;
            double oxState = Math.Clamp(1.0 + dEv * _oxStateSensitivity, 0.0, 4.0);

            return new OperandoPressureControlOutput(valveVolt, oxState, stressRatio, isSafe);
        }
    }
}