AIDT 온디바이스 STEM 페로브스카이트 태양전지 이온 이동 및 J-V 히스테리시스 벡터 서보 엔진

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.Photovoltaics
{
    public readonly ref struct PerovskiteSolarState
    {
        public readonly float CurrentDensityJMaCm2;
        public readonly float InterfaceIonChargeQ;
        public readonly float CompensatedVoltageV;
        public readonly float PowerDensityMwCm2;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public PerovskiteSolarState(float j, float qIon, float vComp, float power)
        {
            CurrentDensityJMaCm2 = j;
            InterfaceIonChargeQ = qIon;
            CompensatedVoltageV = vComp;
            PowerDensityMwCm2 = power;
        }
    }

    public ref struct PerovskiteHysteresisServo
    {
        private const float JscNominal = 24.5f; // mA/cm^2
        private const float VocNominal = 1.15f; // V
        private const float IdealityFactor = 1.6f;
        private const float Vt = 0.02585f;      // thermal voltage at 300K

        private float _accumulatedIonPolarizationV;
        private float _ionRelaxationTimeSec;

        public PerovskiteHysteresisServo(float tauIonSec = 0.5f)
        {
            _accumulatedIonPolarizationV = 0.0f;
            _ionRelaxationTimeSec = MathF.Max(0.01f, tauIonSec);
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public PerovskiteSolarState UpdateServo(float appliedV, float scanRateVPerS, float dt = 0.0001f)
        {
            // Ion migration follows dynamic charging: dV_ion/dt = (V_applied - V_ion)/tau
            float targetPolarization = (appliedV - (VocNominal * 0.5f)) * 0.25f;
            _accumulatedIonPolarizationV += (targetPolarization - _accumulatedIonPolarizationV) * (dt / _ionRelaxationTimeSec);

            // Effective internal voltage considering ion shielding
            float effectiveV = appliedV - _accumulatedIonPolarizationV;

            // Diode characteristic equation with ion modified barrier
            float jDark = 1e-7f * (MathF.Exp(effectiveV / (IdealityFactor * Vt)) - 1.0f);
            float currentJ = MathF.Max(0.0f, JscNominal - jDark);

            // True unperturbed voltage compensated for MPPT tracking
            float compensatedV = appliedV + _accumulatedIonPolarizationV;
            float powerMw = currentJ * appliedV;

            return new PerovskiteSolarState(currentJ, _accumulatedIonPolarizationV, compensatedV, powerMw);
        }
    }
}