A Distributed SOS Program For Local Stability Analysis of Polynomial PDEs in the PIE Representation

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 양자역학, 유체역학, 파동 전파 등 다양한 물리 현상을 모델링하는 데 편미분방정식 (PDE) 이 널리 사용됩니다. 특히 비선형 PDE 의 안정성 분석은 중요한 과제입니다.
• 기존 방법의 한계:
  • 전통적인 PDE 표현은 경계 조건 (BCs) 과 연속성 제약을 상태 변수에 직접 부과합니다.
  • 리아푸노프 함수 (Lyapunov Function, LF) 를 구성할 때, 경계 조건을 처리하기 위해 부분적분 (integration by parts) 이나 포인카레 부등식 (Poincaré inequality) 등을 사용해야 하며, 이는 특정 형태의 PDE 에만 적용 가능합니다.
  • 경계 조건이 변경되거나 (예: 디리클레에서 뉴먼으로), 비선형 항이 복잡해지면 안정성 분석이 매우 어려워지거나 기존 방법이 적용되지 않습니다.
• 핵심 문제: 일반적인 다항식 PDE 에 대해 경계 조건을 명시적으로 처리하지 않고, 체계적이고 자동화된 국소 안정성 (Local Stability) 분석을 수행할 수 있는 프레임워크가 필요합니다.

2. 제안된 방법론 (Methodology)

이 논문은 PDE 의 상태 진화를 **기본 상태 (Fundamental State)**의 진화로 재해석하는 부분적분방정식 (Partial Integral Equation, PIE) 프레임워크를 기반으로 합니다.

A. 기본 상태 (Fundamental State) 와 PIE 표현

• 기본 상태 정의: PDE 상태 $u(t, s)$의 최고 차수 공간 미분 ($v(t) = \partial^n u / \partial s^n$) 을 '기본 상태'로 정의합니다.
• 장점: 기본 상태 $v$는 $L_2$ 공간에 속하며, 원래 PDE 상태와 달리 경계 조건이나 연속성 제약이 없습니다.
• PIE 변환: 선형 PDE 의 경우, 기본 상태에서 PDE 상태로의 매핑은 부분적분 (PI) 연산자로 표현됩니다. 본 논문은 이를 다항식 PDE로 확장합니다.

B. 분산 다항식 (Distributed Polynomials) 과 텐서 -PI 연산자

• 분산 다항식: 기본 상태 $v$에 대한 다항식 동역학을 분산 다항식으로 표현합니다. 이는 $v$의 텐서 곱 ($v^{\otimes k}$) 에 작용하는 새로운 텐서 -PI (Tensor-PI) 연산자를 사용하여 파라미터화됩니다.
• 기능성 -PI 연산자 (Functional-PI Operators): 리아푸노프 함수나 반대수적 집합 (semi-algebraic sets, 예: $L_2$-볼) 을 표현하기 위해 스칼라 값으로 매핑하는 연산자를 정의합니다.
• 벡터화 (Vectorization): 스칼라 T-PI 연산자의 적분을 분산 다항식의 선형 형태로 변환하는 '벡터화' 기법을 도입하여, 등식 제약을 쉽게 부과할 수 있게 합니다.

C. 분산 SOS (Sum-of-Squares) 프로그래밍

• 분산 SOS 다항식: 양수성 (positivity) 을 보장하기 위해 분산 다항식을 SOS 형태로 파라미터화합니다. 이는 유한 차원 다항식의 SOS 프레임워크를 무한 차원 공간 ($L_2$) 으로 확장한 것입니다.
• 국소 안정성 검증:
  1. 리아푸노프 함수 (V): 분산 SOS 다항식으로 정의합니다.
  2. 리 미분 (Lie Derivative): PIE 동역학을 따라가는 $V$의 시간 미분 ($\dot{V}$) 을 계산합니다.
  3. SOS 프로그램 구성: $V > 0$이고 $\dot{V} < 0$ (또는 $-\dot{V}$가 SOS 다항식) 인 조건을 만족하는지 확인하는 최적화 문제를 풉니다. 이는 $L_2$-볼 내부에서의 국소 안정성을 보장합니다.

3. 주요 기여 (Key Contributions)

1. 다항식 PDE 를 위한 새로운 PIE 표현: 선형 및 2 차 PDE 를 넘어, 일반적인 1 차원 다항식 PDE 를 기본 상태의 분산 다항식으로 표현하는 체계를 정립했습니다.
2. 텐서 -PI 및 기능성 -PI 연산자 도입: PDE 의 비선형 항과 리아푸노프 함수의 적분 항을 체계적으로 파라미터화할 수 있는 새로운 연산자 클래스를 정의했습니다.
3. 분산 SOS 프로그래밍 프레임워크: 분산 SOS 다항식을 정의하고, 이를 통해 PDE 의 국소 안정성을 검증하는 SOS 프로그램을 구성했습니다. 이 방법은 경계 조건을 동역학 표현 내부에 내재화하여 처리하므로, 다양한 경계 조건에 대해 일관된 분석이 가능합니다.
4. 이론적 증명: PIE 해와 원래 PDE 해 사이의 일대일 대응 (bijection) 을 증명하여, PIE 상에서의 안정성 분석이 원래 PDE 의 안정성을 보장함을 보였습니다.

4. 수치 실험 및 결과 (Results)

• 대상 시스템: Fisher 방정식 (반응 - 확산 방정식) 을 예제로 사용했습니다.
  • $u_t = u_{ss} + \alpha u - \beta u^2$, 경계 조건: $u(t, 0)=u(t, 1)=0$.
• 실험 설정:
  • 파라미터 $(\alpha, \beta) = (5, -1)$로 설정.
  • 다양한 초기 조건 크기 ($r$, $L_2$-볼 반지름) 에 대해 안정성 영역을 검증.
• 결과:
  • 안정성 반지름: 제안된 SOS 프로그램을 통해 $r \approx 4.048$까지의 $L_2$-볼에서 국소 안정성이 검증되었습니다. 시뮬레이션 결과 $r > 4.05$에서는 불안정성이 관측되어, 얻어진 하한이 매우 정밀 (tight) 함을 확인했습니다.
  • 감쇠율 (Decay Rate): 초기 상태 크기에 따른 지수 감쇠율 $\lambda$의 하한을 계산했습니다. $r$이 작아질수록 (선형 영역에 가까워질수록) $\lambda$는 선형화된 시스템의 고유값 ($\approx 4.869$) 에 수렴하는 것을 확인했습니다.
  • 시뮬레이션: 계산된 감쇠율과 상한선 ($M e^{-\lambda t} \|u_0\|$) 이 실제 PDE 시뮬레이션 결과와 잘 일치함을 Fig. 1 을 통해 입증했습니다.

5. 의의 및 결론 (Significance)

• 일반성: 이 방법은 특정 PDE 유형이나 경계 조건에 국한되지 않고, 1 차원 다항식 PDE 전반에 적용 가능한 범용적인 안정성 분석 도구입니다.
• 자동화 가능성: 경계 조건 처리를 위한 수동적인 부분적분이나 특수한 부등식 적용이 불필요하며, SDP (Semidefinite Programming) 솔버를 통해 자동화된 안정성 검증이 가능합니다.
• 확장성: 현재는 1 차원 스칼라 PDE 에 국한되어 있지만, 이 프레임워크는 결합된 PDE (coupled PDEs) 및 고차원 문제로 확장될 수 있는 잠재력을 가지고 있습니다.
• 실용성: PIETOOLS 와 같은 오픈소스 소프트웨어 패키지에 통합되어, 제어 이론 및 공학 분야에서 비선형 PDE 시스템의 설계 및 분석에 직접 활용될 수 있습니다.

요약하자면, 이 논문은 PDE 의 복잡하고 까다로운 경계 조건 문제를 '기본 상태'와 'PIE'라는 새로운 수학적 표현으로 우회하여, 다항식 SOS 기법을 통해 체계적이고 자동화된 국소 안정성 분석을 가능하게 한 획기적인 연구입니다.