위성 온도 진동의 비선형 해석

초록

위성의 열평형을 단일 온도 변수로 모델링하고, 태양 복사·지구 복사·내부 발열을 입력, 흑체 복사를 출력으로 하는 비선형 미분방정식을 도출하였다. 정성·섭동·수치 해석을 통해 온도가 주기적인 한계궤도(리밋 사이클)로 수렴함을 확인하고, 파라미터에 따라 (i) 주기보다 긴 시간에 느리게 수렴하거나 (ii) 주기보다 짧은 시간에 급격히 수렴하는 두 경우가 존재함을 제시한다. 평균화된 해가 정확히 풀리는 경우도 논의한다.

상세 분석

본 논문은 소형 인공위성을 등온체로 가정하고, 외부에서 들어오는 열원(태양 복사, 행성 복사, 알베도 반사)과 내부에서 발생하는 전자기·기계적 손실을 Q_in(t)라 두었다. 열손실은 표면적 A와 방사율 ε를 갖는 흑체 복사 Q_out = εσA T^4 로 표현한다. 이때 온도 T(t)만을 독립 변수로 하는 비선형 1차 상미분방정식

  C dT/dt = Q_in(t) – εσA T^4

이 도출된다. 여기서 C는 위성의 열용량이며, Q_in(t)는 궤도에 따라 주기적으로 변하는 함수이다. 방정식은 T^4 항 때문에 비선형이며, Q_in(t)의 주기성으로 인해 강제 진동이 발생한다.

비정성 분석에서는 T(t)와 Q_in(t) 사이의 위상공간을 살펴, 모든 궤도 초기조건이 유한한 구역 안으로 수렴하고, 그 안에서 유일한 폐곡선(리밋 사이클)이 존재함을 보인다. 리밋 사이클의 존재와 안정성은 Poincaré–Bendixson 정리를 이용해 증명되며, 시스템이 강제 진동에 의해 주기적인 궤도를 끌어당긴다.

섭동 해법에서는 Q_in(t)를 평균값 Q̄와 작은 진동 성분 q(t)로 분해하고, T = T_0 + ε θ(t) 형태의 전개를 시도한다. 평균 방정식 C dT_0/dt = Q̄ – εσA T_0^4 은 정확히 해석적으로 풀 수 있으며, 정적 평형 온도 T_eq를 제공한다. 작은 진동 θ(t)는 선형화된 방정식 dθ/dt = –(4εσA T_eq^3 / C) θ + q(t)/C 로 기술되며, 감쇠율 λ = 4εσA T_eq^3 / C 가 양수이면 θ는 지수적으로 감쇠한다. λ가 궤도 주기 τ와 비교될 때 두 가지 동역학적 경우가 나타난다.

(i) λ τ ≪ 1 인 경우, 감쇠가 매우 느리므로 온도는 평균값 주변에서 천천히 진동하고, 전체 수렴 시간은 한 주기보다 길다. 이때 평균 방정식이 전체 동역학을 거의 대변한다.

(ii) λ τ ≫ 1 인 경우, 감쇠가 급격히 일어나 온도가 몇 초 안에 리밋 사이클에 도달한다. 실제 궤도 주기보다 짧은 시간에 수렴하므로 초기 조건에 대한 민감도가 크게 감소한다.

수치 시뮬레이션은 Runge–Kutta 4차 방법을 이용해 다양한 파라미터 집합에 대해 T(t)를 계산하고, 위 두 경우의 전이 구간을 확인한다. 시뮬레이션 결과는 이론적 예측과 일치하며, 특히 내부 발열이 큰 경우 λ가 크게 증가해 빠른 수렴이 관찰된다.

결론적으로, 위성 온도는 비선형 방사 손실과 주기적 열입력의 상호작용으로 인해 안정적인 주기적 궤도에 수렴한다. 파라미터에 따라 수렴 속도가 달라지므로, 설계 단계에서 열용량 C와 방사율 ε를 적절히 조정하면 온도 변동을 최소화하거나 원하는 온도 프로파일을 구현할 수 있다.