Error estimates for finite element discretizations of the instationary Navier-Stokes equations

기술 요약: 비정상 나비에-스토크스 방정식의 유한 요소 이산화에 대한 오차 추정

문제 제기
본 연구는 동질 디리클레 (미끄럼 없음) 경계 조건을 갖는 2 차원 비정상 (시간 의존적) 나비에-스토크스 방정식의 수치 해석을 다룬다. 해당 방정식은 시간 방향으로는 불연속 갤러킨 (dG) 방법과 공간 방향으로는 inf-sup 안정 유한 요소를 결합한 완전 이산화 기법을 사용하여 이산화된다. 주요 목적은 완전 이산 해, 특히 $L^\infty(I; L^2(\Omega))$ 노름에 대한 엄밀한 오차 추정을 유도하는 것이다.

최근 문헌에서 비정상 스토크스 문제에 대한 최적 근사형 오차 추정 (특히 $L^\infty(I; L^2(\Omega))$, $L^2(I; H^1(\Omega))$, 및 $L^2(I; L^2(\Omega))$ 노름에서) 이 확립되었으나, 이러한 결과를 비선형 나비에-스토크스 방정식으로 확장하는 것은 상당한 도전을 제기한다. 기존 문헌에서는 종종 $L^\infty(I; L^2(\Omega))$ 및 $L^2(I; H^1(\Omega))$ 오차를 결합된 방식으로 추정하여, $H^1$ 공간 오차의 영향으로 인해 $L^2$ 공간 오차 추정의 차수가 감소하는 결과를 초래했다. 본 논문은 이러한 열화 없이 최적 수렴 차수를 달성하기 위해 $L^\infty(I; L^2(\Omega))$ 오차를 분리하는 것을 목표로 한다.

방법론
저자들은 dG 시간 이산화에 맞춰 조정된 오차 분할, 쌍대성 논증, 그리고 안정성 분석의 조합을 활용한다.

1. 정규성 분석: 본 논문은 먼저 연속 나비에-스토크스 해에 필요한 정규성 속성을 확립한다. 특히 비선형 항 $(u \cdot \nabla)u$ 를 분석한다. 볼록 영역과 충분한 데이터 정규성 ($f \in L^2(I; L^2(\Omega)^2)$, $u_0 \in V$) 에 대한 가정 하에, 저자들은 해가 $H^2$ 공간 정규성을 가지며 비선형 항이 임의의 $s < \infty$에 대해 $L^s(I; L^2(\Omega)^2)$에 속함을 보인다. 이 정규성은 오차 분석에서 비선형성을 처리하는 데 필수적이다.

2. 이산화 기법:

   • 시간: $q$ 차수 (특히 $q=0$ 및 $q=1$) 의 불연속 갤러킨 방법이 사용된다. 이는 시간 노드에서의 연속성 요구 사항 없이 조각별 다항식 근사를 허용한다.
   • 공간: 속도와 압력에 대해 inf-sup 안정 유한 요소 쌍 (예: 테일러-후드 또는 미니 요소) 이 사용된다. 이산 발산 자유 공간 $V_h$가 활용되며, 표준 형태가 $V_h$에서 반대칭성을 결여하고 있으므로 안정성을 보장하기 위해 삼선형 형태의 반대칭 변형 $\hat{c}(\cdot, \cdot, \cdot)$이 도입된다.

3. 안정성 도구: 핵심적인 기술적 기여는 특별히 조정된 이산 그론월 보조정리 (Lemma 4.11) 의 개발이다. 우변에 제곱 노름을 요구하는 일반적인 이산 그론월 보조정리와 달리, 이 버전은 선형 항 (제곱되지 않은 기여) 을 수용하며 시간 방향에서 $L^1$인 우변을 처리하도록 설계되었다. 이를 통해 최종 시간 $T$에 대해 지수적으로 증가하지 않는 안정성 경계를 얻을 수 있으며, 이는 연속 에너지 추정과 부합한다.

4. 오차 분석 전략:

   • 오차 분할: 오차 $u - u_{kh}$는 두 부분으로 분할된다: $\xi = u - \tilde{u}_{kh}$ (스토크스 사영의 오차) 와 $\eta_{kh} = \tilde{u}_{kh} - u_{kh}$ (나머지).
   • 스토크스 성분: 오차 $\xi$는 스토크스 문제에 대한 기존 최적 근사 결과를 사용하여 상한이 설정된다.
   • 나머지 성분: 오차 $\eta_{kh}$는 쌍대성 논증을 사용하여 추정된다. 연속 해와 이산 해의 평균 ($u_{ukh} = \frac{1}{2}(u + u_{kh})$) 을 중심으로 선형화된 이산 쌍대 문제가 구성된다.
   • 쌍대 안정성: 이산 쌍대 문제에 대한 안정성 결과가 $L^\infty(I; L^2(\Omega)) \cap L^2(I; H^1(\Omega))$ 및 더 강한 노름 ($L^\infty(I; H^1(\Omega))$) 에서 유도된다. 이러한 안정성 결과는 이전에 확립된 주 해의 정규성과 이산 그론월 보조정리에 의존한다.

주요 기여 및 결과
본 논문은 충분히 작은 이산화 매개변수 $k$(시간 간격) 와 $h$(메쉬 크기) 에 대해 다음과 같은 주요 결과를 확립한다:

1. $L^2(I; H^1(\Omega))$ 오차 추정:
   $f \in L^2(I; L^2(\Omega)^2)$ 및 $u_0 \in V$에 대해, 오차는 다음을 만족한다:
   $$\|u - u_{kh}\|_{L^2(I; H^1(\Omega))} \leq C(\sqrt{k} + h)$$
   이 결과는 오차 분할 접근법과 표준 안정성 논증을 사용하여 유도된다.

2. $L^\infty(I; L^2(\Omega))$ 오차 추정 (주요 결과):
   $f \in L^\infty(I; L^2(\Omega)^2)$ 및 $u_0 \in V \cap H^2(\Omega)^2$에 대해, 본 논문은 최적 근사형 추정을 증명한다:
   $$\|u - u_{kh}\|_{L^\infty(I; L^2(\Omega))} \leq C \ln\left(\frac{T}{k}\right) \left( \inf_{\chi_{kh} \in V_{kh}} \|u - \chi_{kh}\|_{L^\infty(I; L^2(\Omega))} + \|u - R^S_h(u, p)\|_{L^\infty(I; L^2(\Omega))} \right)$$
   여기서 $R^S_h$는 정상 스토크스 리츠 사영이다.
   특정 정규성 가정 하에, 이는 명시적 수렴 속도 $O(\ln(T/k)^2 (k + h^2))$를 산출한다. 이는 $H^1$ 노름으로부터의 차수 감소 없이 시간과 공간에서 최적 차수를 달성함으로써 이전 추정 (예: $O(\sqrt{k} + h)$) 보다 개선된 것이다.

3. $L^2(I; L^2(\Omega))$ 오차 추정:
   $f \in L^2(I; L^2(\Omega)^2)$ 및 $u_0 \in V$에 대해, 오차는 다음을 만족한다:
   $$\|u - u_{kh}\|_{L^2(I; L^2(\Omega))} \leq C(k + h^2)$$
   이는 $L^\infty$ 추정과 유사한 쌍대성 기법을 사용하지만 쌍대 문제에 대한 더 강력한 안정성 결과를 활용하여 유도된다.

4. 유일성과 안정성:
   본 논문은 충분히 작은 $k$와 $h$에 대해 완전 이산 해의 유일성을 증명하고, 이산 해에 대한 $L^\infty(I; H^1(\Omega))$ 및 이산 스토크스 연산자 $A_h$와 관련된 강력한 안정성 경계를 확립한다.

의의
저자들은 본 연구의 주요 의의가 $L^\infty(I; L^2(\Omega))$ 노름에서 비선형 나비에-스토크스 방정식으로 최적 근사형 오차 추정을 확장하는 데 있다고 주장한다. 이 특정 노름은 시간 방향의 점별 상태 제약을 받는 최적 제어 문제의 분석에 중요하며, 여기서 해당 라그랑주 승수는 시간 방향의 측도이다.

본 논문은 이전 접근법들이 종종 $L^\infty(I; L^2(\Omega))$ 및 $L^2(I; H^1(\Omega))$ 오차를 결합하여 $L^2$ 공간 오차에 대한 비최적 수렴 속도를 초래했다고 강조한다. 정교한 쌍대성 논증과 조정된 이산 그론월 보조정리를 통해 $L^\infty(I; L^2(\Omega))$ 오차를 분리함으로써, 저자들은 유한 요소 공간의 근사 능력과 일치하는 (로그 인자까지) 최적 수렴 차수를 달성한다. 특히 정규성 분석을 통한 비선형 항 처리와 특정 쌍대성 설정을 포함하여 개발된 기법들은 낮은 정규성 데이터를 갖는 완전 이산 나비에-스토크스 문제를 분석하기 위한 견고한 프레임워크를 제공한다.