광자 도플러 속도계로 입자 속도와 크기 분포 복원

초록

본 논문은 가스 매질 내에서 발생하는 금속 충격 유출물(ejecta)의 입자 속도와 크기 분포를 Photon Doppler Velocimetry(PDV) 신호를 이용해 동시에 추정하는 방법을 제시한다. Radiative Transfer Equation(RTE)을 수치적으로 풀어 다중 산란 효과를 포함한 PDV 스펙트로그램을 시뮬레이션하고, 이를 실험 데이터와 비교해 초기 입자 크기 분포를 역추정한다.

상세 분석

이 연구는 기존에 PDV가 단일 산란 가정 하에 입자 속도 분포만을 제공한다는 한계를 넘어, 다중 산란을 정량적으로 기술하는 Radiative Transfer Equation(RTE)을 도입함으로써 입자 크기 정보까지 추출할 수 있음을 증명한다. 저자들은 먼저 충격 파동에 의해 금속 표면에서 발생하는 미세 제트(micro‑jet)와 그 파편화 과정을 설명하고, 초기 입자 크기‑속도 분포 g(a,v)를 크기와 속도가 독립적인 h(a)·j(v) 형태로 가정한다. 이후 가스 매질과의 상호작용을 KIVA‑II 기반의 drag force 모델과 Weber 수에 의한 파편화 모델을 결합한 Phénix 코드로 시뮬레이션한다. 중요한 점은 입자 집단을 ‘수치 입자(numerical particle)’로 대표하여 질량 보존을 유지하면서도 계산 효율성을 확보한 것이다. PDV 신호는 입자들의 도플러 이동에 의해 발생하는 주파수 변화를 검출하고, 이를 Short‑Term Fourier Transform(STFT)으로 시간‑주파수 스펙트로그램 S(t,ω)으로 변환한다. 다중 산란 상황에서는 scattered field가 개별 입자들의 합이 아니라 RTE의 해인 specific intensity I(r,Ω,ν)로 표현되며, 이때 입자 크기 분포는 Mie scattering 효율을 통해 extinction coefficient와 phase function에 직접 영향을 미친다. 저자들은 Phénix 시뮬레이션 결과를 입력으로 RTE를 수치해석(Finite‑Volume 혹은 Discrete‑Ordinate 방법)하여 가스 종류(진공, 헬륨, 공기)별로 예상되는 스펙트로그램을 생성한다. 실험 데이터와의 반복적 비교(iterative fitting)를 통해 초기 h(a) 파라미터를 조정하고, 최종적으로 입자 평균 반경이 수백 나노미터 수준임을 확인한다. 이 과정에서 입자 속도 분포는 Asay foil 측정값을 그대로 사용했으며, 크기 분포만을 가변 파라미터로 두어 역문제의 차원을 크게 줄였다. 결과적으로 PDV는 기존의 ‘속도 전용’ 진단을 넘어, 가스 매질 내에서 입자 크기와 동역학을 동시에 파악할 수 있는 강력한 도구로 재정의된다.