AIDT 온디바이스 STEM 유도 브릴루앙 산란(SBS) 광음향 포논 결합 서보 엔진

1. 개요: 빛과 극초단 음향파(Hypersonic Phonon)의 코히런트 비선형 상호작용

빛이 매질 내를 진행할 때 열역학적 밀도 요동에 의해 산란되는 현상을 자연 브릴루앙 산란(Spontaneous Brillouin Scattering)이라 합니다. 입사 펌프 광의 강도가 특정 임계값을 초과하면, 전자기장의 그래디언트에 의해 매질을 압축하는 전자왜곡(Electrostriction) 효과가 발생하여 강력한 진행형 굴절률 격자(초음파 음향 포논, \(\\sim 10\\,\\text{GHz}\))가 형성됩니다.

이 동적 음향 격자는 후속 펌프 광을 도플러 천이된 후방 산란파(Stokes 파동)로 강력하게 브래그 반사시키며, 반사된 Stokes 파동과 펌프 파동의 간섭 비트(Beat)가 다시 음향파를 증폭시키는 양성 피드백 루프를 형성합니다. 이것이 유도 브릴루앙 산란(Stimulated Brillouin Scattering, SBS)이며, 수십 킬로미터에 이르는 장거리 광케이블의 온도 및 물리적 변형률(Strain)을 센티미터 단위로 분산 감지하는 핵심 원리입니다.

2. 3파 결합 파동 방정식 및 브릴루앙 주파수 천이 수식 모델

단일 모드 광섬유에서 전방 진행 펌프 파동 \(E_p(z)\), 후방 진행 Stokes 파동 \(E_s(z)\), 음향 밀도 파동 \(\\rho(z)\) 사이의 3파 상호작용(Three-wave interaction) 방정식은 다음과 같습니다:

\\[ \\frac{d I_p}{dz} = -g_B I_p I_s - \\alpha I_p \\]

\\[ -\\frac{d I_s}{dz} = g_B I_p I_s - \\alpha I_s \\]

여기서 \(\\alpha\)는 광섬유 전송 손실 계수, \(g_B(\\nu)\)는 로렌츠 형상의 브릴루앙 이득 스펙트럼(Brillouin Gain Spectrum, BGS)입니다:

\\[ g_B(\\nu) = g_0 \\frac{(\\Delta \\nu_B / 2)^2}{(\\nu - \\nu_B)^2 + (\\Delta \\nu_B / 2)^2} \\]

후방 산란 브릴루앙 주파수 천이(Brillouin Frequency Shift, \(\\nu_B\))는 음속 \(V_a\) 및 코어 굴절률 \(n_{\\text{eff}}\)와 직접 연계됩니다:

\\[ \\nu_B = \\frac{2 n_{\\text{eff}} V_a}{\\lambda_p} \\]

온도 변화 \(\\Delta T\)와 축 변형률 \(\\Delta \\epsilon\)에 따른 주파수 천이는 정밀 선형 계수를 나타냅니다:

\\[ \\nu_B(T, \\epsilon) = \\nu_{B0} + C_T \\Delta T + C_\\epsilon \\Delta \\epsilon \\]

실리카 표준 단일모드 광섬유(SMF-28)에서 \(\\nu_{B0} \\approx 10.85\\,\\text{GHz}\) (1550 nm 기준), \(C_T \\approx 1.07\\,\\text{MHz/}^\\circ\\text{C}\), \(C_\\epsilon \\approx 0.05\\,\\text{MHz/}\\mu\\epsilon\)입니다.

3. 실시간 BOTDA 분산 감지 서보 아키텍처 다이어그램

graph TD A[협선폭 1550nm DFB 레이저] --> B[광 커플러 펌프/프로브 분기] B --> C[EOM 광 강도 펄스 생성기 펌프 광 입사 z=0] B --> D[광 주파수 변조기 RF 스위퍼 프루브 광 입사 z=L] C & D --> E[광섬유 내 역방향 충돌 및 SBS 음향 포논 격자 증폭] E --> F[광 수신 광검출기 BGS 스펙트럼 수집] F --> G[10,000Hz 0-Alloc C# 로렌츠 피팅 서보 엔진] G --> H[거리별 온도 Delta T & 변형률 Delta eps 분산 매핑] H --> D

4. 고성능 0-Alloc C# SBS 로렌츠 피크 추적 서보 구현

BGS 이득 스펙트럼으로부터 중심 천이 주파수 \(\\nu_B\)를 피팅하고 온도와 변형률을 실시간 추출하는 0-Alloc C# 엔진입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.BrillouinScattering
{
    public readonly record struct SbsParameters(
        double CenterFrequencyHz, // 기본 nu_B0 ~ 10.85 GHz
        double GainLinewidthHz,   // Delta nu_B ~ 30 MHz
        double TempCoefficient,   // Hz / K (예: 1.07e6)
        double StrainCoefficient  // Hz / microstrain (예: 50e3)
    );

    public ref struct SbsServoState
    {
        public double EstimatedNuB;
        public double PeakGain;
        public double TemperatureCelsius;
        public double StrainMicrostrain;
        public bool IsConverged;
    }

    public static class StimulatedBrillouinServo
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void FitLorentzianPeak(
            in SbsParameters p,
            ReadOnlySpan freqScanHz,
            ReadOnlySpan measuredGain,
            double baselineTempC,
            ref SbsServoState state)
        {
            if (freqScanHz.Length < 3 || freqScanHz.Length != measuredGain.Length)
            {
                state.IsConverged = false;
                return;
            }

            // 1. 최대 이득 위치 탐색 (초기 추정)
            double maxG = -1.0;
            int maxIdx = 0;
            for (int i = 0; i < measuredGain.Length; i++)
            {
                if (measuredGain[i] > maxG)
                {
                    maxG = measuredGain[i];
                    maxIdx = i;
                }
            }

            // 2. 포물선 삼각 보간으로 서브-빈(Sub-bin) 중심 주파수 추정
            if (maxIdx > 0 && maxIdx < measuredGain.Length - 1)
            {
                double y1 = measuredGain[maxIdx - 1];
                double y2 = measuredGain[maxIdx];
                double y3 = measuredGain[maxIdx + 1];

                double denom = 2.0 * (2.0 * y2 - y1 - y3);
                double delta = (denom != 0.0) ? (y1 - y3) / denom : 0.0;

                double stepHz = freqScanHz[1] - freqScanHz[0];
                state.EstimatedNuB = freqScanHz[maxIdx] + delta * stepHz;
                state.PeakGain = y2 + 0.125 * (y1 - y3) * delta;
                state.IsConverged = true;
            }
            else
            {
                state.EstimatedNuB = freqScanHz[maxIdx];
                state.PeakGain = maxG;
                state.IsConverged = false;
            }

            // 3. 온도 천이 분리 계산 (변형률 구속 상태 가정 시)
            double deltaNu = state.EstimatedNuB - p.CenterFrequencyHz;
            state.TemperatureCelsius = baselineTempC + (deltaNu / p.TempCoefficient);
            state.StrainMicrostrain = deltaNu / p.StrainCoefficient;
        }
    }
}

5. 교육 연수 및 국방 인프라 방호 적용 지침

방통대 직무연수 및 국가 중요시설 광학 경계망

물리/정보 교과 연계: 음향 포논-광자 비선형 결합, 브래그 산란 조건, 로렌츠 선형 피팅 및 수십 킬로미터 광섬유 OTDR 신호 처리 원리 교육.
국방 해저 케이블 및 군사기지 울타리 감시: 수중 해저 통신망의 도청 침입 시도에 따른 미세 굽힘 변형률 감지 및 잠수함 해안 침투 음향 진동 모니터링.