열복사 점원 모델로 푸에리어 이상 현상 해명하기

초록

**
본 논문은 소수의 점‑형 열원으로 구성된 간단한 모델을 이용해 파이오니어 10·11호의 열복사에 의한 추진력을 계산하고, 파라미터 민감도 분석을 통해 열효과가 관측된 비정상 가속도의 상당 부분을 설명할 수 있음을 보인다.

**

상세 분석

**
이 연구는 파이오니어 탐사선의 복잡한 구조를 완전한 유한요소(FEM) 모델로 재현하기 어려운 현실을 인정하고, 대신 “점‑형 열원” 접근법을 채택한다. 각 열원은 등방성 혹은 Lambertian(코사인‑법칙) 방사 특성을 갖으며, 방출 전력 W와 위치 (x₀,y₀,z₀)만 지정하면 시간 평균 포인팅 벡터 S를 통해 복사압과 반작용력을 직접 계산한다. 방사압은 p_rad = α S·n /c 로 표현되며, α = 1은 완전 흡수, α = 2는 완전 확산 반사를 의미한다.

핵심은 복사된 광선이 다른 표면에 닿아 그림자를 만들거나 재반사되는 과정을 적분함으로써 전체 추진력을 구한다는 점이다. 저자들은 이 과정을 수치적으로 구현하기 위해 표면을 여러 개의 작은 면소자로 분할하고, 각 면에 일정 수(1, 4, 16, 64, 144)의 점‑형 열원을 배치한다. 테스트 케이스에서는 1 m² 평면에서 방사된 열이 다른 평면에 미치는 힘을 계산했으며, 점‑형 열원의 수가 늘어날수록 정확도가 급격히 향상돼 16 개 이상이면 오차가 1 % 이하로 수렴함을 확인했다.

모델 적용 범위는 주로 RTG(방사성 동력원), 메인 버스, 안테나 디시 세 부분이다. 전체 전력 소비량은 잘 알려져 있으므로, 각 부품에 할당된 전력을 기반으로 열원 강도를 추정한다. 주요 가정으로는 (1) 탐사선이 스핀 안정화돼 축 외 방향의 힘이 평균적으로 소멸한다, (2) 외부 표면의 온도 구배는 미미하고, (3) 재료의 광학적 특성(반사율·흡수율)은 시간에 따라 크게 변하지 않는다. 이러한 가정은 민감도 분석을 통해 검증되었으며, 전력 배분 비율이나 α 값의 변동이 최종 가속도 추정에 미치는 영향이 정량적으로 제시된다.

결과적으로, 저자들은 열복사에 의한 가속도가 a_th ≈ (5 ~ 8) × 10⁻¹⁰ m/s² 로, 관측된 a_P ≈ 8.5 × 10⁻¹⁰ m/s² 와 비교적 근접함을 보여준다. 특히 RTG에서 방출된 열이 안테나 디시와 버스에 반사·흡수되는 복합 효과가 전체 가속도의 약 30 % ~ 40 % 를 차지한다는 점을 강조한다. 이는 이전 연구에서 제시된 “열효과가 전체 가속도의 1/3 정도를 설명한다”는 주장과 일치한다.

하지만 모델은 (가) 표면 손상·오염에 따른 방사율 변동, (나) 내부 구조물의 미세한 열전도 경로, (다) 장기적인 방사성 붕괴에 따른 전력 감소 등을 정밀히 반영하지 못한다는 한계가 있다. 또한, 점‑형 열원 배치가 실제 구조와 완전히 일치하지 않으므로, 절대적인 오차 범위는 ± 0.5 × 10⁻¹⁰ m/s² 정도로 추정된다. 이러한 불확실성을 감안하면, 열복사만으로 전체 이상 가속도를 완전히 설명하기는 아직 어려우며, 추가적인 시스템 오류나 새로운 물리 현상의 가능성을 완전히 배제할 수는 없다.

**