반사열이 파일론 탐사선 가속에 미치는 영향 모델링

초록

**
본 논문은 파일론 10·11 탐사선의 열복사에 의한 미세 가속도를 보다 정확히 계산하기 위해 확산·정반사 두 가지 반사 모델을 포함한 새로운 방식을 제시한다. 점광원과 Phong‑쉐이딩 기법을 이용해 복잡한 구조의 열반사를 정량화하고, 파라미터 스터디를 통해 열가속도의 상한·하한을 추정한다. 결과는 기존 연구와 일치하며, 열효과가 관측된 비정상 가속도의 상당 부분을 설명함을 확인한다.

**

상세 분석

**
논문은 파일론 10·11 탐사선의 비정상 가속도(≈ 8.74 × 10⁻¹⁰ m s⁻²)를 열복사 현상으로 설명하려는 일련의 연구 흐름에 기반한다. 기존 연구에서는 라마르틴 점광원을 이용해 방사능을 추정했으나, 반사 효과를 단순히 흡수·반사 계수 α(1~2)로 대체하는 수준에 머물렀다. 본 연구는 이를 넘어, 광학에서 널리 쓰이는 Phong‑쉐이딩 모델을 차용해 확산 반사(k_d)와 정반사(k_s, α) 두 성분을 명시적으로 계산한다.

핵심 수식은 방사 Poynting 벡터 S(x) = W/(π‖x−x₀‖²) · n·(x−x₀) · (x−x₀)/‖x−x₀‖ 로, 여기서 W는 방출 전력, n은 표면 법선이다. 이 벡터를 기반으로 표면에 입사하는 에너지 흐름 E_ilum와 방사압 p_rad을 정의하고, 반사된 복사장을 S_rd, S_rs 로 각각 확산·정반사 성분에 대해 전개한다. 특히 정반사 항은 (r·(x−x′))^α 형태의 ‘shininess’ 지수를 도입해 거울성 정도를 조절한다.

모델 구현 단계는 다음과 같다. (1) 파일론 구조를 고전적인 주요 부품(고이득 안테나, 장비함, 두 개의 라디오‑열 발생기(RTG))으로 단순화하고, 각 부품을 다수의 라마르틴 점광원으로 대체한다. (2) 각 점광원에서 방출된 복사장을 계산하고, 다른 표면에 입사하는 부분을 적분해 E_ilum을 구한다. (3) Phong‑반사식을 적용해 반사된 복사장을 구하고, 다시 입사·반사 과정을 반복한다. (4) 전체 복사 압력을 표면에 대해 적분해 총 힘을 얻고, 이를 질량으로 나누어 가속도를 산출한다.

파라미터 스터디에서는 k_d, k_s, α 값을 현실적인 범위(예: k_d = 0.60.9, k_s = 0.10.4, α = 10~100)로 변동시켜 열가속도의 상·하한을 탐색한다. 결과는 열가속도가 초기 전력(~2500 W)에서 약 30 %~70 % 수준으로 감소하며, 전력 감소율(≈ 1 %/년)과 일치하는 시간적 감쇠를 보인다. 또한, 반사 효과를 포함했을 때 가속도 상한이 약 5 × 10⁻¹⁰ m s⁻², 하한이 2 × 10⁻¹⁰ m s⁻²로, 관측값의 대부분을 포괄한다는 점을 확인한다.

열전도 효과도 간단히 평가했는데, RTG와 장비함 사이, 장비함과 안테나 사이의 전도율을 고려한 결과 전도에 의한 전력 전달은 1~4 W 수준으로 전체 열복사에 비해 무시할 수 있는 수준임을 보여준다.

이러한 접근법은 (1) 복잡한 구조의 다중 반사 경로를 정량화, (2) 계산 효율성을 유지하면서도 물리적 파라미터를 명시적으로 다룰 수 있다는 장점을 갖는다. 따라서 파일론 탐사선뿐 아니라, 유사한 열복사‑반사 문제가 있는 장거리 우주선 설계·분석에 적용 가능성이 크다.

**