Spectral Koopman Approach for Reconstructing State-space Geometry of Cislunar Restricted 3-Body Problem

기술 요약: 지구-달 제한 3체 문제의 상태 공간 기하학 재구성을 위한 스펙트럼 코업만(Koopman) 접근법

문제 정의
지구, 달, 태양의 중력적 상호작용에 의해 지배되는 시즐루나(cislunar) 환경은 리브레이션 포인트(libration points), 불변 다양체(invariant manifolds), 그리고 혼돈적 수송(chaotic transport)을 특징으로 하는 고도의 비선형 동역학계입니다. 이 영역에 대한 전통적인 미션 설계 및 궤적 분석은 원형 제한 3체 문제(CR3BP)와 수치 최적화 기법(예: 미분 보정, 멀티 슈팅)에 의존합니다. 그러나 이러한 기존 방식들은 세 가지 지속적인 한계에 직면해 있습니다:

1. 충실도 불일치(Fidelity Mismatches): 단순화된 모델과 고충실도 실제 물리 현상 사이의 차이.
2. 계산 비용(Computational Cost): 항법 업데이트가 발생할 때마다 막대한 재계산이 필요함.
3. 확장성 문제(Scalability Issues): 장기 호라이즌의 불확실성 전파를 위한 몬테카를로 방법의 낮은 확장성.

또한, 표준적인 국소 선형화 기법은 공명(resonance)이나 수송 채널와 같은 강력한 비선형 구조에 의해 조직된 전역적 위상 공간 기하학을 포착하는 데 실패합니다. 연산자 이론적 접근 방식(특히 코업만 연산자)은 비선형 시스템에 대한 전역적인 선형 표현을 제공하지만, 기존의 데이터 기반 근사법(예: 확장 동적 모드 분해, EDMD)은 정교한 기저(basis) 선택이 필요하며 고유한 전역 스펙트럼을 정확하게 회복하는 데 어려움을 겪습니다.

방법론
저자들은 임의의 사전 정의된 딕셔너리(dictionary) 선택 없이 평면 CR3BP의 전역 위상 공간 기하학을 재구성하기 위해, **경로 적분 공식화 기반의 코업만 스펙트럼(path integral formulation of the Koopman spectrum)**을 사용하는 새로운 접근 방식을 제안합니다.

1. 코업만 연산자 프레임워크: 시스템 동역학 $\dot{x} = f(x)$는 관측량(observables)의 진화를 통해 분석됩니다. 코업만 연산자 $U_t$는 비선형 흐름에 대한 선형 기술을 제공합니다. 주요 고유함수 $\phi_\lambda$는 $\frac{\partial \phi_\lambda}{\partial x}f = \lambda \phi_\lambda$ 방정식을 만족합니다.
2. 경로 적분 공식화: 유한 차원 부공간에서 연산자를 근사하는 대신, 저자들은 주요 고유함수를 직접 계산합니다. 벡터장을 선형($Ax$) 및 비선형($F_n(x)$) 성분으로 분해함으로써, 비선형 부분인 고유함수$h_\lambda(x)$는 다음과 같이 경로 적분으로 표현됩니다:
   $$h_\lambda(x) = \int_0^\infty e^{-\lambda \tau} w_\lambda^\top F_n(s_\tau(x)) d\tau$$
   여기서 $s_\tau(x)$는 시스템 궤적을 나타냅니다.
3. 수치적 구현: 고유함수는 유한 시간 호라이즌 $T$ 동안 적분을 평가함으로써 근사됩니다. 불안정 고유함수( $\lambda_u > 0$ 와 관련된 경우)는 순방향 궤적을 사용하여 계산되며, 안정 고유함수( $\lambda_s < 0$ 와 관련된 경우)는 역방향 궤적을 사용하여 계산됩니다.
4. 기하학 재구성: 저자들은 이러한 주요 고유함수들의 결합 제로-레벨 곡선(joint zero-level curves)이 시스템의 평형점(라그랑주 점)과 관련된 불변 다양체(안정, 불안정, 중심 다양체)에 대응한다는 성질을 활용합니다.

핵심 기여

• 전역 선형 표현: 본 논문은 코업만 프레임워크가 국소 선형화의 한계를 극복하고 CR3BP의 전역 선형 표현을 생성할 수 있음을 입증합니다.
• 고유함수의 직접 계산: 경로 적분 방식을 활용함으로써, 저자들은 EDMD와 같은 방식에서 요구되는 수동적인 딕셔너리 설계 과정을 우회하여 시계열 데이터로부터 주요 고유함수를 직접 계산합니다.
• 불변 다양체 식별: 본 연구는 계산된 고유함수의 제로-레벨 집합이 위상 공간의 특성을 인코딩하며, 특히 다섯 개의 라그랑주 점($L_1$ ~ $L_5$)과 관련된 안정, 불안정, 중심 다양체를 식별함을 확립합니다.
• 수송 채널 시각화: 이 접근법은 한 평형점의 불안정 다양체가 다른 평형점의 안정 다양체와 어떻게 정렬되어 수송 채널을 형성하는지를 보여줌으로써, 시즐루나 환경의 "보이지 않는 고속도로"를 성공적으로 시각화합니다.

결과
제안된 방법은 평면 4D CR3BP 모델(지구-달 시스템)에 적용되었습니다. 시뮬레이션 결과는 다음과 같습니다:

• 제로 속도 곡선(ZVCs)의 재구성: 이 방법은 야코비 상수(Jacobi constant)에 의해 정의된 에너지 허용 영역과 금지 구역을 복구합니다.
• 평형점에서의 국소 기하학:
  • $L_1$ 및 $L_2$: 고유함수는 수송 관문을 명확하게 구분합니다. 불안정 다양체는 외부로 뻗어 나가며(지구 탈출 또는 외부 진입 허용), 안정 다양체는 내부로 휘어지며(궤적을 포획) 작용합니다.
  • $L_3$: 외부 영역으로 누출되는 취약한 다양체를 나타내며, 강력한 수송 통로를 형성하지 못합니다.
  • $L_4$ 및 $L_5$: 중심 고유함수가 지배하며, 이는 선형 안정성을 가진 삼각형 점(triangular points)과 일치하는 폐곡선 형태의 타폴로드(tadpole) 및 트로잔(Trojan) 궤적을 생성합니다.
• 전역 위상 공간 분할: 고유함수의 제로-레벨 곡선은 위상 공간을 스펙트럼적으로 분할합니다. 에너지 기반의 야코비 표면이 운동이 가능한 곳을 정의한다면, 코업만 곡선은 궤적이 해당 영역 내에서 어떻게 점근적으로 조직되는지를 드러내어 분리선(separatrix)과 불변 경계를 식별합니다.
• 수송 메커니즘: 결과는 $L_1$과 $L_2$ 사이의 수송 채널을 보여주며, 여기서 $L_1$의 불안정 다양체와 $L_2$의 안정 다양체가 불변 튜브를 형성하여 지구와 달 영역 간의 이동을 용이하게 합니다.

의의 및 주장
본 논문은 스펙트럼 코업만 접근법이 개별 궤적에서 관측량의 진화로 초점을 전환함으로써 시즐루나 역학을 위한 강력한 도구를 제공한다고 주장합니다. 주요 의의는 다음과 같습니다:

1. 전역 기하학 복구: CR3BP의 전역 위상 공간 기하학을 시계열 데이터로부터 재구성하여, 국소적 방법이 놓치는 비선형 구조를 포착합니다.
2. 다운스트림 작업 지원: 불확실성 전파, 도달 가능성 분석(reachability analysis), 제어 설계의 토대를 제공합니다. 코업만 연산자는 선형이기 때문에, 저자들은 비선형 CR3BP에 대한 불확실성 분석 및 제어 전략 설계가 전통적인 비선형 방식에 비해 더 단순해지고 확장 가능해질 것이라고 단언합니다.
3. 데이터 기반 발견: 우주 역학 분야에서 데이터 기반 발견의 새로운 방향을 제시하며, 연산자의 선형성을 통해 본질적으로 비선형적인 시스템에 대한 제어 전략의 분석 및 설계를 가능하게 합니다.

저자들은 현재의 연구가 CR3BP의 기하학적 구조를 복구하는 데 집중하고 있으나, 이 프레임워크가 복잡한 다체 환경에서의 실시간 궤적 최적화 및 강건 제어(robust control)를 위한 미래 응용 분야에 대한 높은 가능성을 지니고 있다고 결론짓습니다.