A priori error estimates for optimal control problems governed by the transient Stokes equations and subject to state constraints pointwise in time

기술적 요약: 점별 시간 상태 제약이 있는 과도 스토크스 방정식의 최적 제어에 대한 사전 오차 추정

문제 정식화
본 논문은 과도 스토크스 방정식에 의해 지배되는 상태 제약 최적 제어 문제를 다룹니다. 목적은 원하는 프로파일로부터의 상태 편차와 제어에 대한 정규화 항을 포함하는 2 차 비용 함수를 최소화하는 것입니다. 시스템은 다음 조건을 따릅니다:

1. 상태 방정식: 동질적인 초기 조건 및 경계 조건을 갖는 과도 스토크스 방정식.
2. 제약 조건: 공간과 시간에서 제어 변수 $q(t,x)$ 에 대한 점별 상한.
3. 상태 제약 조건: 시간의 점별 상한을 만족해야 하는 상태의 $L^2(\Omega)$ 함수에 대한 제약. 구체적으로, 모든 $t \in [0, T]$ 에 대해 $\int_\Omega u(t,x) \cdot w(x) \, dx \leq \beta$를 만족해야 합니다.

영역 $\Omega \subset \mathbb{R}^d$ ($d=2,3$) 은 볼록한 다각형 또는 다면체 영역으로 가정합니다.

방법론 및 이산화
저자들은 공간 및 시간 근사를 결합한 완전 이산화 전략을 사용합니다:

• 공간 이산화: 속도와 압력 공간에 대해 Inf-sup 안정 유한 요소 (예: Taylor-Hood 또는 MINI 요소) 를 사용합니다. 분석은 표기법을 단순화하기 위해 이산 발산 자유 형식을 활용합니다.
• 시간 이산화: 0 차 불연속 갈러킨 (dG) 방법 (시간에 대해 조각별 상수) 을 적용합니다.
• 제어 이산화: 두 가지 접근 방식을 분석합니다:
  1. 변분 이산화: 제어 공간은 사전에 이산화되지 않으며, 상태 방정식만 이산화됩니다.
  2. 완전 이산화: 제어는 상태 및 쌍대 변수의 격자와 일치하는 공간과 시간 모두에서 조각별 상수 함수를 사용하여 이산화됩니다.

분석은 $L^\infty(I; L^2(\Omega))$ 노름에서 과도 스토크스 방정식에 대해 최근 확립된 최접근 유형 오차 추정에 의존합니다. 중요한 기술적 요구 사항은 이산 Slater 점의 존재를 보장하기 위해 공간 메쉬 크기 $h$와 시간 단계 $k$ 사이의 약한 결합 조건, 구체적으로 $h \to 0$일 때 $\ln(T/k)h \to 0$이어야 한다는 조건입니다.

주요 기여 및 이론적 결과

1. 최적 제어의 정칙성: 본 논문은 최적 제어 $\bar{q}$에 대한 개선된 정칙성을 확립합니다. 데이터에 대한 추가 가정 없이, $\bar{q}$가 $L^2(I; H^1(\Omega)^d) \cap C(\bar{I}; L^2(\Omega)^d)$에 속함이 보입니다. 특정 조건 하에 (예: 제어 제약의 부재 또는 가중 함수 $w$의 특정 성질), 상태 제약과 관련된 라그랑주 승수가 디랙 기여를 포함하지 않음이 보여지며, 이는 쌍대 상태와 제어에 대한 더 높은 정칙성을 의미합니다.
2. 사전 오차 추정: 주요 이론적 결과는 $L^2(I \times \Omega)$ 노름에서 최적 제어에 대한 사전 오차 추정입니다. 완전 이산 문제에 대해, 연속 최적 제어 $\bar{q}$와 이산 해 $\bar{q}_\sigma$ 사이의 오차는 다음을 만족합니다:
   $$\|\bar{q} - \bar{q}_\sigma\|_{L^2(I \times \Omega)} \leq C \frac{\ln(T/k)}{\sqrt{\alpha}} (k^{1/2} + h)$$
   여기서 $C$는 $k$와 $h$에 무관하며, $\alpha$는 정규화 매개변수입니다. 이 추정은 제어 제약이 있는 경우에 성립하며, 해의 고차 정칙성을 사전에 가정하는 대신 최적성 체계에서 유도된 정칙성에 의존합니다.
3. 이산 잘-제정성: 저자들은 연속 조건에서 유도된 이산 Slater 조건을 활용하여 변분 이산화 및 완전 이산 최적 제어 문제에 대한 해의 존재성과 유일성을 증명합니다.

수치 결과
이론적 발견은 2 차원 및 3 차원 수치 실험으로 보완됩니다:

• 매끄러운 데이터: 매끄러운 해석적 데이터를 가진 예시는 공간에서 $O(h)$의 이론적 수렴 차수와 시간에서 약 $O(k^{1/2})$의 수렴 차수를 보여줍니다. 결과는 시간 단계가 메쉬 크기에 비해 충분히 세분화되지 않으면 공간 이산화 오차가 지배적일 수 있음을 강조합니다.
• 낮은 정칙성 데이터: 거친 데이터 (부드럽지 않은 원하는 상태 및 가중 함수) 를 가진 예시는 시간에서 감소된 수렴 속도 (약 $O(k^{0.6})$) 를 보여주며, 이는 초수렴 효과가 실제 속도를 가릴 수 있는 매끄러운 경우보다 이론적 기대치인 $O(k^{1/2})$와 더 잘 부합합니다.
• 제어 제약: 세 번째 예시는 활성 제어 제약을 포함합니다. 수치 결과는 제어 제약의 존재가 비제약 경우에 대해 유도된 수렴 속도를 저하시키지 않음을 확인합니다.

의의 및 범위
본 논문은 점별 시간 상태 제약이 있는 과도 스토크스 최적 제어 문제에 대한 오차 추정을 유도하는 것이 주요 기여라고 주장하며, 이 설정은 이전에 그러한 엄격한 분석이 부족했습니다. 이 작업은 포물형 방정식 (열 방정식) 에 대한 기존 결과를 스토크스 체계로 확장합니다. 저자들은 그들의 분석이 최적 제어에 대한 강한 정칙성 가정을 요구하지 않으며, 대신 최적성 조건으로부터 필요한 정칙성을 유도한다고 강조합니다. 유도된 오차 추정은 특히 정규화 매개변수 $\alpha$, 시간 단계 $k$, 메쉬 크기 $h$ 간의 상호작용과 관련하여 예시에서 관찰된 수치적 성능에 대한 이론적 기초를 제공합니다. 본 논문은 연속 및 이산 문제를 분석하지만, 정적 나비에 - 스토크스 방정식에 대한 오차 추정을 유도하거나 상태 제약 하의 유체 제어의 수학적 프레임워크를 넘어선 새로운 응용을 제안하지는 않는다고 명시적으로 언급합니다.