Stabilization of Parabolic Time-Varying PDEs using Certified Reduced-Order Receding Horizon Control

기술 요약: 인증된 축소 차수 예측 제어기를 이용한 포물형 시간 가변 PDE 의 안정화

문제 제기
본 논문은 유한 차수 예측 제어 (RHC, 일명 모델 예측 제어 MPC) 를 사용하여 선형 시간 가변 포물형 편미분방정식 (PDE) 의 안정화 문제를 다룬다. PDE 에 RHC 를 적용할 때의 주요 과제는 겹치는 시간 구간에서 고차원 오픈-루프 최적 제어 문제를 반복적으로 해결하는 데 따른 계산 비용이다. 축소 차수 모델 (ROM) 은 이러한 계산을 가속화할 수 있는 길을 제공하지만, ROM 에서 유도된 제어 입력을 직접 전체 차수 모델 (FOM) 에 적용할 경우 원래 시스템의 폐루프 안정성 및 부분 최적성 보장이 훼손될 수 있다. 또한, 본 논문은 제어 입력의 희소성을 강제하기 위해 제곱된 $\ell_1$-노름을 포함하는 비연속적 비용 함수를 고려하며, 이는 최적화 및 오차 분석의 복잡성을 더한다.

방법론
저자들은 전체 차수 시스템에 대한 엄밀한 안정성 분석과 인증된 모델 축소 전략을 결합한 2 단계 프레임워크를 제안한다.

1. 전체 차수 모델 (FOM) 분석:
   저자들은 힐베르트 공간에서 시간 가변 포물형 시스템에 대한 연속 시간 RHC scheme 의 지수적 안정성과 부분 최적성을 먼저 확립한다. 이전 연구들이 안정성을 보장하기 위해 종단 비용이나 제약 조건에 의존한 것과 달리, 본 프레임워크는 특정 연결 (concatenation) scheme 과 최종 시간 관측성 부등식을 활용한다. 이 접근법은 초기 시간 관측성 추정치를 요구하지 않으면서도 더 약한 단계 비용 (특히 $L^2(\Omega)$-추적 항) 과 비연속적 제어 비용 (희소성 장려 정규화 포함) 의 사용을 가능하게 한다.

2. 축소 차수 모델링 (ROM) 및 오차 추정:
   계산 비용을 해결하기 위해 저자들은 적절한 직교 분해 (POD) 를 사용하는 투영 기반 갈르킨 (Galerkin) 모델 차수 축소를 적용한다. 유한 시간 구간 최적 제어 문제에 대한 엄밀한 사후 (a posteriori) 오차 분석이 개발되었다. 여기에는 상태, 수반 상태, 최적 제어 입력, 그리고 가치 함수에 대한 오차 추정기 도출이 포함된다. 중요하게도, 이러한 추정치는 비연속적 볼록 제어 비용과 부정확한 초기 값을 고려한다.

3. 인증된 ROM-RHC 알고리즘:
   핵심 기여는 적응형 ROM-RHC 알고리즘 (알고리즘 2) 이다. 축소된 제어 입력을 맹목적으로 적용하는 대신, 알고리즘은 유도된 사후 오차 추정기를 사용하여 온라인에서 "완화된 동적 계획법 (Relaxed Dynamic Programming, RDP)" 조건을 검증한다.

   • 각 시간 단계에서 알고리즘은 축소된 문제를 해결하고 ROM 가치 함수 및 오차 추정기에 기반한 성능 상한 ($\bar{\alpha}^r_k$) 을 계산한다.
   • 성능 상한이 사용자가 지정한 임계값 $\tilde{\alpha}$를 만족하면, 축소된 제어 입력이 전체 시스템에 적용된다.
   • 조건이 충족되지 않으면 축소 기저가 정제된다 (예: 전체 차수 해로부터 새로운 POD 스냅샷을 추가). 그 후 단계가 반복된다.
     이 적응형 과정은 축소된 제어 입력이 무한 시간 구간 FOM 문제에 대한 보장된 성능 상한을 가지면서 전체 차수 시스템을 안정화하도록 보장한다.

주요 기여
본 논문은 네 가지 주요 기여를 제공한다:

1. FOM 에 대한 안정성 분석: $L^2(\Omega)$-추적 항과 비연속적 비용을 사용하여 시간 가변 선형 포물형 시스템의 연속 시간 RHC 에 대한 지수적 안정성과 부분 최적성을 확립한다. 이 분석은 초기 시간 관측성이 아닌 최종 시간 관측성 부등식에 의존하며, 부분 최적성 매개변수에 대한 새로운 명시적 상한을 제공한다.
2. 엄밀한 오차 추정: 저자들은 가치 함수 및 비용 함수의 ROM 근사에 대한 사후 오차 추정식을 유도한다. 이러한 추정식은 일반적인 비연속적 볼록 제어 비용과 부정확한 초기 값에 적용 가능하며, 완화된 정규성 가정 하에서 갈르킨 ROM 의 보간 성질과 점근적 수렴을 확립한다.
3. 인증된 ROM-RHC scheme: 새로운 알고리즘이 제안되었으며, 이는 이러한 오차 추정식과 완화된 동적 계획법 원리의 변형을 결합한다. 이 scheme 은 FOM 의 지수적 안정화와 무한 시간 구간 FOM 문제에 대한 부분 최적성을 보장한다. 이 방법은 RDP 부등식이 유한 차수 기저로 만족되도록 ROM 이 구성될 수 있음을 보장한다.
4. 수치적 검증: 지수적으로 불안정한 자유 동역학과 비연속적 $\ell_1$-정규화를 가진 시스템에서 접근법이 검증되었다. 결과는 ROM-RHC 가 전체 차수 방법의 안정성 및 부분 최적성 보장을 유지하면서 상당한 계산 속도 향상을 달성함을 보여준다.

결과
13 개의 작동기와 비연속적 비용 함수를 가진 2 차원 포물형 PDE 에 대한 수치 실험은 이론적 발견을 확인한다:

• 안정성: FOM 및 ROM-RHC scheme 모두 이론적 상한과 일치하는 지수적 감쇠율로 시스템을 성공적으로 안정화한다.
• 정확도: 축소 차수 해는 전체 차수 해에 비해 극히 낮은 상대 오차를 보인다 (예: 비용 함수의 상대 오차는 $10^{-7}$에서 $10^{-8}$ 수준).
• 효율성: ROM-RHC scheme 은 FOM-RHC 에 비해 계산 시간에서 약 10~13 배의 속도 향상을 달성한다. 이는 전체 차수 해가 모델 정제가 필요할 때만 트리거되므로 전체 차수 기울기 평가 횟수가 약 20 배 감소함으로써 달성된다.
• 오차 추정: 사후 오차 추정기는 신뢰할 수 있는 것으로 입증되었으나, 특히 큰 예측 시간 구간과 작은 정규화 매개변수의 경우 실제 오차를 과대평가하는 경향이 있다. 그러나 추정기는 축소 기저 차수가 증가함에 따라 오차의 지수적 감쇠를 정확하게 포착한다.

의의 및 주장
본 논문은 시간 가변 포물형 PDE 의 효율적이고 안정적인 제어를 위한 강력한 프레임워크를 제공한다고 주장한다. 엄밀한 사후 오차 분석과 적응형 모델 축소를 통합함으로써, 이 방법은 이론적 보장을 희생하지 않으면서 고차원 PDE 에 대한 표준 RHC 의 계산적 불가능성을 극복한다. 저자들은 안정성을 보장하기 위해 종단 비용이나 제약 조건이 필요하지 않으므로 비연속적 목적 함수를 포함하는 더 넓은 범위의 문제에 적용 가능하다고 강조한다. 인증 메커니즘은 축소 차수 제어기가 단순한 근사가 아닌 정량화된 성능 상한을 가진 검증된 전체 차수 시스템 안정화기임을 보장한다. 이 작업은 더 복잡한 시스템을 처리하는 한 걸음으로 제시되며, 향후 비선형 포물형 시스템, 매개변수 불확실성, 적응형 시간 구간 선택에 대한 연구가 제안된다.