배경지향 슬리렌 측정에서 밀도 추정의 불확실성 정량화

초록

본 논문은 배경지향 슬리렌(BOS) 실험에서 점 패턴의 변위를 교차상관 기반 PIV 불확실성 추정 기법으로 정량화하고, 이를 포아송 방정식 풀이와 결합해 밀도 필드의 지역·즉시 불확실성 한계를 제공한다. 합성 가우시안 밀도장과 초음속 웨지 흐름 실험을 통해 변위 불확실성 전파와 경계조건의 영향을 검증하였다.

상세 분석

본 연구는 BOS 측정 체인의 모든 단계에 걸쳐 불확실성을 체계적으로 추적한다는 점에서 기존 밀도 측정 방법과 차별화된다. 먼저, BOS에서 얻어지는 점 패턴 변위는 PIV와 동일한 교차상관 절차를 사용하므로, PIV에서 개발된 직접형 불확실성 추정법(이미지 매칭, 상관통계, 상관모멘트)을 그대로 적용한다. 저자는 간접형 방법의 보정 의존성을 피하기 위해 직접형 방법만을 선택했으며, 특히 상관통계(CS)와 상관모멘트(MC) 기법을 중심으로 변위 표준편차를 계산한다. 변위 불확실성 σ_u는 광학 파라미터(배율 M, 점패턴‑밀도구배 거리 Δz)와 결합되어 투영 밀도 구배 σ_∇ρ = (σ_u)/(M·Δz·K·n₀) 형태로 전파된다(여기서 K는 Gladstone‑Dale 상수, n₀는 주변 굴절률).

밀도 구배를 얻은 뒤에는 포아송 방정식 ∇²ρ = –∇·(∇ρ) 를 이산화하여 선형 시스템 Aρ = b 로 변환한다. 저자는 2차 중심차분 스키마를 사용해 라플라시안 연산자 A와 소스항 b를 희소 행렬 형태로 구현함으로써 대규모 이미지에서도 계산 효율을 확보한다. 불확실성 전파는 공분산 행렬 Σ_b 를 A⁻¹·Σ_b·(A⁻¹)ᵀ 로 전파하는 방식으로 수행되며, 경계조건(디리클레)에 대한 라벨 행렬 L을 포함해 경계점의 확정값을 반영한다. 이 과정에서 공분산의 대각 원소 √(Σ_ρ,ii) 가 각 격자점의 밀도 불확실성 σ_ρ,i 로 해석된다.

합성 실험에서는 가우시안 형태의 밀도장을 사용해 실제 변위와 추정 변위, 그리고 전파된 밀도 불확실성을 비교하였다. 결과는 전파된 σ_ρ가 실제 오차와 높은 상관관계를 보이며, 특히 밀도 구배가 급격히 변하는 영역 주변에서 σ_ρ가 크게 증가함을 확인한다. 실험적 검증으로는 초음속 유동이 웨지 위를 흐르는 경우를 선택했으며, 충격파와 경계층에서 발생하는 급격한 밀도 변화가 전체 시야에 걸쳐 불확실성을 확대시키는 현상을 관찰했다. 특히 포아송 솔버를 이용한 경우, 디리클레 경계가 지정된 영역으로부터 멀어질수록 σ_ρ가 단조 증가하는 패턴을 보였는데, 이는 라플라시안 연산자의 전역적 특성에 기인한다.

이와 같이 변위 단계에서부터 밀도 재구성 단계까지 일관된 선형·확률 모델을 적용함으로써, BOS 데이터에 대한 신뢰 구간을 실시간으로 제공할 수 있다. 또한, 불확실성 전파가 사용된 밀도 적분 방식(포아송 vs. 다른 적분 스키마)에 따라 결과가 달라짐을 보여, 향후 연구에서는 비선형 적분 기법이나 베이지안 역학 모델과의 결합 가능성을 제시한다.