기술 요약: 고충실도 태양 복사압 모델링을 통한 GPU 가속 궤도 전파

문제 정의

정확한 궤도 전파는 우주 임무 설계, 기동 계획 및 추진제 최적화, 특히 우주 임무의 복잡성과 수가 증가함에 따라 매우 중요합니다. 우주선 궤도에 영향을 미치는 지배적인 비중력 섭동 중 하나는 태양 복사압(Solar Radiation Pressure, SRP)입니다. 단순화된 해석적 모델(예: Cannonball, N-Plate)은 계산 비용은 낮지만, 자기 차폐(self-occlusion), 다중 반사, 관절형 태양전지판 운동과 같은 복잡한 기하학적 효과를 포착하지 못합니다. 이러한 단순화는 장기 전파 시 누적 오차를 발생시킵니다.

반대로, 고충실도 광선 추적(ray-tracing) 모델은 상세한 우주선 메쉬(mesh)와 광자 상호작용을 시뮬레이션하여 정확한 SRP 추정치를 제공합니다. 그러나 이러한 모델은 수천에서 수백만 번의 SRP 평가가 필요한 수치 궤도 전파기에 통합될 때 계산 비용이 지나치게 높아집니다. 또한, 기존 도구들은 전파 과정 중에 우주선 기하 구조(예: 움직이는 태양전지판)를 동적으로 업데이트하거나, 동일한 조건 하에서 보간 기반 근사치와 온라인 고충실도 계산을 체계적으로 비교하는 기능이 부족한 경우가 많습니다

방법론

저자들은 기존의 HiFi-SoRaP 프레임워크를 확장하여, 고충실도 SRP 모델을 통합하기 위한 세 가지 상호 보완적인 전략을 다루는 완전한 궤도 전파 환경을 구축했습니다:

1. Vulkan 기반 GPU 가속: 기존의 OpenGL 기반 광선 추적 구현을 Vulkan 컴퓨트 파이프라인으로 포팅했습니다. 이러한 전환은 메모리 관리, 동기화 및 컴퓨트 스케줄링에 대한 명시적인 제어를 제공합니다. 주요 최적화 사항은 다음과 같습니다:

   • GPU 기반 리덕션(Reduction): GPU와 CPU 사이의 데이터 전송 오버명을 최소화하기 위해 서브그룹 연산을 사용하여 GPU에서 직접 SRP 기여도를 누적합니다.
   • 수치적 일관성: 수치적 안정성을 유지하기 위해 Kahan 보상 합산(Kahan compensated summation)을 포함합니다.
   • 병렬 실행: 독립적인 컴퓨트 셰이더를 사용하여 이산화된 이미지 평면상의 픽셀당 광선을 처리합니다.

2. 가동형 태양전지판 모델링: 우주선 모델을 루트가 있는 유향 트리 표현($G = (N, E, L)$)을 사용하여 관절 구조를 지원하도록 확장했습니다. 이를 통해 각 SRP 평가 전에 기하 구조를 동적으로 업데이트할 수 있습니다. 알고리즘은 트리를 순회하여 가동 조인트를 식출하고, 태양 방향에 맞게 패널 법선(normal)을 정렬하기 위한 회전(각도 제한 적용)을 계산하며, 관련 메쉬 서브트리에 강체 변환을 적용합니다.

3. 보간 기반 근사(SPAD): 정확도와 비용 간의 트레이드오프를 평가하기 위해, 프레임워크는 고충실도 GPU 모델을 통해 생성된 오프라인 데이터셋으로부터 SRP 힘을 검색하는 사전 계산된 SRP 데이터셋(SPAD)을 활용합니다.

실험 설정:
프레임워크는 두 가지 우주선 모델, 즉 기하학적으로 단순한 "BoxWing" 위성(212개 삼각형)과 고충실도 Parker Solar Probe(PSP) 모델(39,257개 삼각형)을 사용하여 검증되었습니다. 시뮬레이션은 명시적 Runge-Kutta 7/8(RK78) 적분기를 사용하여 헬리오센트릭 타원 궤도(장반경 0.557 AU, 이심률 0.682)에서 수행되었습니다. 연구에서는 다음을 비교했습니다:

• 직접 RayTrace GPU (OpenGL vs. Vulkan).
• SPAD 보간 (다양한 샘플링 그리드 사용).
• 고정형 vs. 가동형 태양전지판 구성.

주요 기여

본 논문은 궤도 전파 분야에 세 가지 주요 기여를 합니다:

1. Vulkan 기반 고충실도 SRP 엔진: 원래의 OpenGL 기반 방식의 수치적 정확성을 유지하면서 직접 광선 추적 SRP 평가를 가속화하는 새로운 구현을 제공합니다.
2. 동적 기하 구조 통합: 정적 기하 구조 가정에서 벗어나, 전파 중 관절형 태양전지판의 연속적인 방향 변화를 고려하도록 SRP 모델을 확장했습니다.
3. 체계적 비교 프레임워크: 해석적 모델, 보간 기반 접근 방식(SPAD), 그리고 온라인 고충실도 광선 추적을 동일한 전파 조건 하에서 직접 비교할 수 있는 통합된 환경을 제공합니다.

결과

1. 수치적 정확도 및 Vulkan 성능

• 수치적 일관성: Vulkan 구현은 계산된 SRP 가속도의 상대적 차이를 $5 \times 10^{-4}$ 미만으로 유지함으로써 기존 OpenGL 방식의 정확도를 보존합니다.
• 속도 향상:
  • 개별 SRP 계산의 경우, GPU 기반 리덕션을 적용한 Vulkan 구현은 복잡한 PSP 모델에서 기존 OpenGL 구현 대비 9.4배, 단순한 BoxWing 모델에서 4.4배의 속도 향상을 달축했습니다.
  • 전체 궤도 전파(20 회전)의 경우, PSP 모델에 대해 총 실행 시간이 최대 15.2배 단축되었습니다 (약 5시간 18분에서 약 20분 53초로 감소).
  • 성능 이득은 광선 추적 계산이 실행 시간을 지배하는 복잡한 기하 구조에서 가장 컸으며, 이는 고정된 통신 오버헤드의 영향을 최소화했기 때문입니다.

2. 보간(SPAD) vs. 직접 평가

• 정확도: 20x10 샘플링 그리드를 사용한 SPAD 기반 보간은 직접 RayTrace GPU 평가의 궤도를 높은 충실도로 재현했습니다. 원 궤도에서 위치 차이는 약 0.12 km였으며, 이는 N-Plate 모델의 12.06 km와 대조됩니다.
• 효율성: 보간은 궤도 전파의 계산 비용을 몇 시간에서 밀리초 단위로 줄였습니다 (수 차수의 속도 향상).
• 한계: 정확도는 샘플링 그리드에 따라 달라집니다. 이중 선형 보간이 필요한 이동된 그리드(SPAD 20x10S)는 샘플링 지점과 정렬된 그리드보다 큰 오차를 보였으나, 여전히 해석적 모델보다는 정확했습니다. 연속적인 자세 변화가 보간에 가장 까다로운 시나리오로 확인되었습니다.

3. 가동형 태양전지판의 영향

• 정확도: 패널 운동을 무시하면 특히 면적 질량비가 큰 우주선의 경우 상당한 오차가 발생합니다.
  • PSP의 경우, SRP 힘의 평균 상대 오차는 약 7.5%였으며, 이는 20 회전 후 약 27,590 km의 궤도 편차를 초래했습니다.
  • BoxWing 모델의 경우, 평균 상대 오차는 2.75 (275%)였으며, 이는 약 317만 km의 궤도 편차를 발생시켰습니다.
• 계산 오버헤드: 가동형 패널을 포함하면 광선-삼각형 교차 횟수(ray-triangle intersections)가 증가하여 광선 추적 작업량이 늘어납니다.
  • PSP의 경우, 오버헤드는 완만했습니다 (전파 시간 약 8% 증가).
  • BoxWing의 경우, 움직이는 대형 배열로 인해 교차 횟수가 크게 변하여 오버헤드가 더 높았습니다 (전파 시간 약 55% 증가).

의의 및 가이드라인

본 논문은 SRP 모델링 전략의 선택이 계산 효율성, 수치적 정확도 및 물리적 충실도 사이의 균형을 맞춰야 한다고 결론짓습니다:

• SPAD 보간은 우주선의 자세가 사전 계산된 데이터셋의 범위 내에 있는 경우, 일상적인 장기 궤도 전파 및 임무 분석에 권장됩니다. 이는 매우 유리한 정확도 대비 성능 비율을 제공합니다.
• 직접 고충실도 광선 추적(Vulkan에 의해 가속됨)은 최대의 물리적 충실도가 필요하거나, 참조 데이터셋을 생성할 때, 또는 조명 조건이 보간의 유효 범위를 벗어날 때 정당화됩니다. Vulkan 구현은 이전에는 너무 비용이 많이 들었던 복잡한 우주선에 대해서도 이를 가능하게 합니다.
• 가동형 패널 모델링은 관절형 태양전지판 배열을 가진 우주선에 필수적입니다. 패널 운동을 무시하면 장기적인 궤도 오차가 크게 발생할 수 있으며, 시뮬레이션의 현실성 개선을 고려할 때 관련 계산 비용은 수용 가능한 수준입니다.

저자들은 실용적인 워크플로를 제안합니다: 직접 GPU 광선 추적(필요한 경우 관절 구조 포함)을 통해 고충실도 참조 데이터를 생성하고, 일상적인 전파에는 SPAD 기반 보간을 사용하며, 보간의 유효성이 훼손될 때만 직접 평가로 전환하는 방식입니다.