비선형 베이지안 도플러 토모그래피: 흐름과 온도 동시 재구성

초록

본 논문은 비선형 가우시안 프로세스 토모그래피(GPT)와 라플라스 근사를 결합한 베이지안 프레임워크를 제시한다. 전체 도플러 전방 모델을 유지하면서 로그‑가우시안 프로세스 사전분포를 이용해 저방사 영역에서의 속도 발산을 억제하고, 방사율·이온 온도·플로우 속도를 동시에 복원한다. 합성 팬텀 데이터와 RT‑1 장치의 코히어런스 이미징 스펙트로스코피(CIS) 실험 데이터를 통해 방법의 정확성과 실용성을 검증하였다.

상세 분석

이 연구는 도플러 스펙트럼 이미징의 비선형 역문제를 베이지안 관점에서 해결한다. 먼저 물리량(방사율 e, 이온 온도 Ti, 플로우 속도 vi)을 공간적으로 연관된 가우시안 프로세스(GP)로 모델링하고, 커널은 위치에 따라 변하는 길이 스케일 ℓ(r)을 갖는 Gibbs 커널을 채택해 플라즈마 내부의 비균질성을 반영한다. 전방 모델은 방사율과 온도에 의한 도플러 브로드닝, 속도에 의한 도플러 시프트를 완전하게 포함하는 비선형 함수 g(f)이며, 이를 직접 사용함으로써 기존의 선형화 가정이나 작은 시프트 가정을 탈피한다. 비선형성으로 인해 사후분포는 비가우시안이 되므로 라플라스 근사를 적용해 최빈값(모드)와 헤시안의 역을 이용해 가우시안 근사(posterior≈N(μ̂,Σ̂))를 얻는다. 특히 속도 재구성에 로그‑가우시안 프로세스 사전(p(v)∝exp(−½(v−μ)ᵀK⁻¹(v−μ)))을 도입해 방사율이 낮아 도플러 정보가 약해지는 영역에서도 물리적으로 타당한 속도 값을 유지한다. 최적화는 뉴턴‑라프슨 방식으로 수행되며, 단계별로 방사율(로그‑방사율 ê)과 온도·속도 결합 변수(ã, v̂)를 순차적으로 추정·마진화한다. 하이퍼파라미터(σ_f,ℓ(r) 등)는 증거 최적화 혹은 교차 검증을 통해 결정한다. 합성 팬텀 실험에서는 실제와 비교해 평균 제곱 오차가 5 % 이하로 감소했으며, RT‑1 실험에서는 토로이달 플로우와 온도 프로파일이 기존 방법보다 더 선명하게 재현되었다. 이 프레임워크는 강한 플로우와 큰 온도 변동이 존재하는 플라즈마 진단뿐 아니라, 천체 물리, 대기 레이더, 의료 초음파 등 다양한 도플러 기반 영상에 적용 가능성을 제시한다.