[논문 리뷰] Guaranteed optimal reachability control of reaction-diffusion equations using one-sided Lipschitz constants and model reduction

Adrien Le Coënt, Laurent Fribourg|arXiv (Cornell University)|2019. 07. 25.

Model Reduction and Neural Networks참고 문헌 39인용 수 4

한 줄 요약

이 논문은 명시적 오일러 시간 이산화와 한쪽 경계 리프시츠(One-Sided Lipschitz, OSL) 상수를 사용하여 1D 반응-확산 편미분방정식에 대한 보장된 최적 제어 방법을 제시한다. OSL 상수가 음수일 경우, 시간 영역 길이에 대해 수렴성과 선형 오차 경계를 증명한다. 또한 더 큰 시스템에서의 계산 효율성을 향상시키기 위해 모델 축소 기법을 도입한다.

ABSTRACT

We show that, for any spatially discretized system of reaction-diffusion, the approximate solution given by the explicit Euler time-discretization scheme converges to the exact time-continuous solution, provided that diffusion coefficient be sufficiently large. By "sufficiently large", we mean that the diffusion coefficient value makes the one-sided Lipschitz constant of the reaction-diffusion system negative. We apply this result to solve a finite horizon control problem for a 1D reaction-diffusion example. We also explain how to perform model reduction in order to improve the efficiency of the method.

연구 동기 및 목표

• 보장된 도달 가능성과 함께 유한 시간 영역 최적 제어 문제를 해결하기 위한 비선형 반응-확산 PDE에 대한 도전 과제를 해결한다.
• 시스템의 OSL 상수가 음수일 경우, 시간 이산화된 해의 수치적 안정성과 수렴성을 보장한다.
• 명시적 오일러 적분과 모델 순서 축소를 조합하여 계산 효율성이 높은 방법을 개발한다.
• 이전 연구에서의 지수적 경계와는 대조적으로, 시간 영역 길이에 따라 선형적으로 증가하는 엄밀한 오차 경계를 제공한다.
• 투영 기반 모델 축소를 통해 더 큰 공간 이산화에 대한 최적 제어 적용을 가능하게 한다.

제안 방법

• 반응-확산 PDE를 공간적으로 반연속화하여 상미분방정식의 시스템으로 변환한다.
• 수치적 해를 구하기 위해 명시적 오일러 시간 적분을 사용하며, 시스템의 OSL 상수가 음수일 경우 수렴이 보장된다.
• 한쪽 경계 리프시츠(One-Sided Lipschitz, OSL) 상수를 활용하여 시간 영역 T에 대해 선형인 사전 오차 경계를 유도한다.
• 적절한 직교 분해(POD)를 사용하여 모델 축소를 수행하고, 고차원 시스템을 저차원 부분공간으로 투영한다.
• 최적 제어 전략을 저차원 시스템에서 합성한 후, 오차가 유한하게 유지되는 방식으로 전체 시스템에 적용한다.
• 투영 기반 오차 추정을 사용하여 전체 시스템과 축소 시스템의 궤적 간의 이격도를 유한하게 제한하며, 오차는 $ \frac{K_2 \sigma}{|\lambda_{h_1}|} $ 비례한다.

실험 결과

연구 질문

• RQ1 특정 조건 하에서 명시적 오일러 시간 이산화가 반응-확산 PDE의 정확한 해로 수렴할 수 있는가?
• RQ2 한쪽 경계 리프시츠 상수가 PDE의 유한 시간 최적 제어에서 오차 경계에 어떤 영향을 미치는가?
• RQ3 반응-확산 시스템에 대해 모델 축소를 어떻게 적용하여 계산 비용을 줄이면서도 제어 성능을 유지할 수 있는가?
• RQ4 축소된 시스템과 전체 시스템 궤적 간의 오차 경계를 정량적으로 유한하게 제한할 수 있는가? 이는 보장된 제어 합성에 기여하는가?
• RQ5 제안된 방법이 이전 접근 방식의 지수적 스케일링과는 대조적으로 시간 영역 길이에 대해 선형 오차 스케일링을 보이는가?

주요 결과

• 확산 계수가 충분히 크면 OSL 상수가 음수가 되어 명시적 오일러 스킴이 반응-확산 PDE의 정확한 해로 수렴한다.
• 최적 제어 해의 오차 경계는 이전 연구에서 발견된 지수적 의존성 $ O(e^{L_f T}) $ 와는 대조적으로 시간 영역 T에 대해 선형적으로 증가한다.
• 투영 기반 모델 축소를 통해 시스템 크기를 크게 줄일 수 있다. 예를 들어 $ M_2 = 10 $ 에서 $ M_1 = 5 $ 로 축소되며, 격자 크기는 약 $ 15^5 \approx 7.6 \times 10^5 $ 배 감소하여 대규모 제어를 가능하게 한다.
• 축소 오차는 $ \frac{K_2 \sigma}{|\lambda_{h_1}|} $ 로 유한하게 제한되며, 여기서 $ \sigma $ 는 POD 절단 오차를 나타내며, 이 경계는 투영된 궤적이 목표 근처의 알려진 이웃 영역 내에 머무르도록 보장한다.
• 1D 이중 안정 반응-확산 시스템에 대한 실험 결과는, 축소 모델에서 유도된 제어 전략이 전체 시스템을 목표 상태로 성공적으로 이끌며, 시뮬레이션 결과가 전체 시스템과 밀도 있게 일치함을 보여준다.
• 후행 시뮬레이션 결과는 이론적 오차 경계가 보수적으로 평가되었음을 보여주며, 실제 궤적 간의 거리는 상한 경계보다 훨씬 작다.