에너지 안전을 보장하는 조기 종료 가능한 포트‑해밀토니안 병렬 결합

본 논문은 대규모 로봇 시스템의 실시간 병렬 시뮬레이션 및 제어에 필수적인 파티셔닝(time‑stepping) 기법에서 발생하는 에너지 비보존 문제를 해결하기 위해, 포트‑해밀토니안(Port‑Hamiltonian) 서브시스템들을 스캐터링(파동) 변수와 Douglas‑Rachford(DR) 분할 스키마에 결합한 새로운 인터페이스 알고리즘을 제안한다. 논문은 먼저 포트‑해밀토니안 시스템을 전통적인 (effort, flow) 쌍 (e_i, f_i) 대신 파동 변수 (a_i, b_i) 로 변환한다. 여기서 a_i = e_i + γ f_i, b_i = e_i – γ f_i (γ>0) 로 정의하고, 역변환을 통해 전력은 P_i = ½(‖a_i‖² – ‖b_i‖²) 로 표현된다. 이 변환은 파동 변수의 노름 차이가 바로 전력임을 명시함으로써, 손실 없는 인터커넥션을 단순한 직교 변환 a = P b (PᵀP=I) 로 기술할 수 있게 만든다. 따라서 파동 변수 공간에서의 인터페이스 연산은 에너지 중립(energy‑neutral)이며, 중간 단계에서 에너지 주입이 발생하지 않는다. 다음으로, 각 서브시스템에 대해 1‑스텝 통합기(예: symplectic Euler, implicit midpoint 등)를 적용했을 때, 파동 변수 입출 관계 b_i = S_n^i(a_i) 라는 고정된 포트 매핑이 정의된다. 매크로 타임스텝 n에서 전체 시스템은 두 제약식 b = S_n(a) 와 a = P b 를 동시에 만족해야 하는 고정점 문제로 귀결된다. 이를 직접 풀면 전역적인 비선형 연립방정식이 되지만, 논문은 Douglas‑Rachford 스플리팅을 이용해 내부 반복(inner‑loop)으로 해결한다. 구체적으로, ζ = col(a, b) 를 정의하고, 직교 투영 Π_C(·) 를 통해 a = P b 제약을 강제한다. DR 업데이트는 ζ^{k+1}=ζ^{k}+Π_C(2y^{k}−ζ^{k})−y^{k} (y^{k}=col(a^{k}, S_n(a^{k}))) 형태이며, 각 반복 단계에서 서브시스템 포트 매핑 S_n^i는 병렬로 평가된다. 초기값은 이전 매크로 스텝의 출력 파동 b^{0} 로 설정하고, K_n 번의 반복을 수행하거나 수렴 기준 ε 이하가 되면 조기 종료한다. 내부 반복이 끝난 뒤에는 최종 입파동 a_n = P b̂_n,⋆ 를 사용해 각 서브시스템의 실제 상태와 출력 파동을 업데이트한다. 핵심 이론적 기여는 두 가지 정리이다. 첫 번째 정리(Theorem 1)는 “증강 저장 부등식”을 제시한다. 여기서는 각 서브시스템 통합기가 자체적으로 패시브(에너지 감소가 포트 전력보다 작음)라는 조건(Condition 1)과 포트 매핑 S_n이 강하게 비팽창적이라는 조건(Condition 2)을 가정한다. 이 두 조건 하에, 내부 반복 횟수 K_n이 어떠한 값이든, 전체 매크로 스텝에 대해 H̃_{n+1} ≤ H̃_n 이 성립한다. 여기서 H̃_n = Σ_i H_i(x_i^n) + ½‖b^{k}−Π_C b^{k}‖² 로 정의된 증강 저장 함수는 파동 변수의 불일치(즉, 인터페이스 잔차)를 에너지 형태로 포함한다. 따라서 제한된 내부 반복만 수행해도 시스템은 이산적으로 패시브하며, 에너지 안전성을 보장한다. 두 번째 정리(Theorem 2)는 K_n → ∞ 일 때 DR 반복이 고정점을 정확히 찾아, b = S_n(P b) 를 만족하게 함으로써 모놀리식 암시적 통합과 동일한 업데이트를 얻는다고 증명한다. 이는 제안된 방법이 “anytime” 특성을 가지면서도, 반복 횟수를 늘리면 정확도가 단조 증가한다는 것을 의미한다. 실험에서는 2‑차원 커플드 진동자 시스템을 사용해 제안된 알고리즘을 검증하였다. 각 서브시스템에 대해 symplectic Euler 통합기를 적용하고, γ와 P를 적절히 튜닝하였다. 결과는 다음과 같다. (1) 에너지 보존 오류는 double precision 환경에서 10⁻¹⁴ 수준으로, 정밀도 한계에 도달했다. (2) RMS 상태 오차는 K_n=0(명시적 스캐터링 인터페이스)에서 시작해 K_n=2, 4, 8 로 증가할수록 단조 감소했으며, K_n=∞(즉 충분히 큰 반복)에서는 모놀리식 해와 거의 일치했다. (3) 내부 반복은 각 서브시스템의 포트 매핑을 병렬로 평가하므로, 전체 연산 비용은 K_n에 비례하지만, 병렬 하드웨어에서 효율적으로 스케일한다. 논문의 마지막 섹션에서는 제한 사항과 향후 연구 방향을 논의한다. 현재는 선형 직교 인터커넥션(PᵀP=I)과 고정된 임피던스 파라미터 γ에 기반하고 있으며, 비선형 또는 비직교 인터페이스에 대한 확장은 추가적인 수학적 분석이 필요하다. 또한, 조건 2(강한 비팽창성)가 모든 고차 통합기에 자동으로 만족되는지는 아직 완전히 밝혀지지 않았다. 향후 연구에서는 비선형 포트 매핑에 대한 일반화, 적응형 γ 튜닝, 그리고 실시간 MPC와 같은 고속 제어 응용에 대한 하드웨어 구현을 목표로 한다. 전반적으로, 이 논문은 파동 변수와 DR 스플리팅을 결합해 “조기 종료 가능한 에너지 안전 인터페이스”를 구현함으로써, 고속 로봇 시뮬레이션 및 제어 분야에서 실시간 병렬 계산을 가능하게 하는 중요한 이론적·실용적 기여를 제공한다.