Analysis of a finite element method for second order uniformly elliptic PDEs in non-divergence form

weifeng Qiu 의 논문 "Analysis of a Finite Element Method for Second Order Uniformly Elliptic PDEs in Non-Divergence Form"에 대한 상세한 기술적 요약입니다.

1. 문제 제기

이 논문은 $\Omega \subset \mathbb{R}^d$ ($d=2, 3$) 인 Lipschitz 다면체 영역에서 정의된 두 가지 유형의 2 차 편미분방정식 (PDE) 에 대한 수치해석을 다룹니다:

1. 비발산 형태의 선형 타원형 PDE:
   $$A : D^2u + b \cdot \nabla u + cu = f \quad \text{in } \Omega, \quad u = 0 \text{ on } \partial\Omega$$
   여기서 $A$는 대칭적이고 균일하게 양정치인 행렬이며, $b, c$는 유계 계수입니다.

2. 균일 타원형 Hamilton-Jacobi-Bellman (HJB) 방정식:
   $$\sup_{\alpha \in \Lambda} [A_\alpha : D^2u + b_\alpha \cdot \nabla u + c_\alpha u - f_\alpha] = 0 \quad \in \Omega, \quad u = 0 \text{ on } \partial\Omega$$
   여기서 $\Lambda$는 인덱스 집합이며, 계수들은 균일 타원성과 유계 조건을 만족합니다.

주요 과제:

• 비발산 형태: 2 차 항 $A:D^2u$를 $A$나 영역의 높은 정칙성을 가정하지 않고는 테스트 함수로 쉽게 이동시킬 수 없기 때문에, 표준적인 약형 공식화 (적분변환에 기반) 는 직접 적용할 수 없습니다.
• 해 공간: 분석은 점성해 (viscosity solutions) 가 아닌 Sobolev 공간 $W^{2,p}(\Omega)$의 **강해 (strong solutions)**에 초점을 맞춥니다.
• 영역 기하학: 이 논문은 기존 문헌에서 종종 요구되는 $C^{1,1}$ 또는 볼록 영역보다 더 일반적인 Lipschitz 다면체 영역, 특히 2 차원의 비볼록 다각형을 처리하는 것을 목표로 합니다.
• 계수 정칙성: 이 방법은 최적의 수렴 속도를 유지하면서 (특히 HJB 방정식의 경우) 불연속 계수를 수용해야 합니다.

2. 방법론

저자는 선형 PDE 와 HJB 방정식 모두에 대한 통합된 유한요소법 (FEM) 을 제안합니다.

A. 이산화 및 연산자

• 메쉬: 준균일 (quasi-uniform), 심플렉스, 정합 (conforming) 삼각분할 $\mathcal{T}_h$.
• 유한요소 공간: $V_h = H^1_0(\Omega) \cap P_r(\mathcal{T}_h)$, 여기서 다항식 차수 $r \geq 2$.
• 이산 기울기 ($\nabla_h$): 이 방법은 메쉬 요소에서의 적분변환과 내부 면을 가로지르는 점프 항을 통해 정의된 특정 이산 기울기 연산자를 사용합니다. 이 연산자는 기저 함수의 $C^1$ 연속성을 요구하지 않고 헤시안 $D^2u$를 근사할 수 있게 합니다.
  $$(\nabla_h v)_i = \partial_{h, x_i} v$$
• 가중 함수: 비발산 형태를 처리하기 위해 방정식에 가중 함수 $\gamma$ (HJB 의 경우 $\gamma_\alpha$) 를 곱합니다.
  • 선형 경우: $\gamma(x) = \frac{\text{Tr} A(x)}{|A(x)|^2}$.
  • 이는 방정식을 이산 연산자에 적합하도록 변환합니다: $\tilde{A} : D^2 u + \dots = \tilde{f}$.

B. 수치적 기법

• 선형 경우 (식 14): 모든 $v_h \in V_h$에 대해 다음을 만족하는 $u_h \in V_h$를 찾습니다:
  $$\int_\Omega (\tilde{A} : \nabla_h(\nabla u_h) + \tilde{b} \cdot \nabla u_h + \tilde{c} u_h) v_h \, dx = \int_\Omega \tilde{f} v_h \, dx$$
• HJB 경우 (식 16): 모든 $v_h \in V_h$에 대해 다음을 만족하는 $u_h \in V_h$를 찾습니다:
  $$\left( \sup_{\alpha \in \Lambda} \gamma_\alpha [A_\alpha : \nabla_h(\nabla u_h) + b_\alpha \cdot \nabla u_h + c_\alpha u_h - f_\alpha], v_h \right)_\Omega = 0$$
  참고: $\Lambda$가 단일 원소를 포함하면 HJB 기법은 정확히 선형 기법으로 축소됩니다.

C. 분석적 프레임워크

분석은 다음과 같이 정의된 이산 $W^{2,p}$-노름에 의존합니다:
$$\|v\|_{W^{2,p}_h(\Omega)} = \|D^2_h v\|_{L^p(\Omega)} + \left( \sum_{F \in \mathcal{F}_h^I} h_F^{1-p} \|\llbracket \nabla v \rrbracket\|_{L^p(F)}^p \right)^{1/p}$$
주요 분석 도구는 다음과 같습니다:

1. 전역 $W^{2,p}$ 추정: Lipschitz 다면체에서 연속 문제의 강해에 대한 잘-제시됨 (well-posedness) 을 증명합니다.
2. 이산 컴팩트성: 이산 $W^{2,p}_h$ 노름에서 유계인 수열이 $W^{1,p}$에서 수렴하는 부분수열을 가진다는 것을 확립합니다.
3. 유일성 논증: 이전 연구에서 사용된 특수한 쌍대성 논증 (duality arguments) 을 피하고, 강해의 유일성을 이용하여 안정성을 증명합니다 (모순을 통해).
4. 상确의 가측성: HJB 방정식의 경우, 인덱스 집합 $\Lambda$에 대한 상确이 가측이며 가산 부분집합 $\tilde{\Lambda}$로 근사될 수 있음을 증명합니다 (가정 6).

3. 주요 기여

A. 통합 프레임워크

이 논문은 단일 FEM 공식화가 선형 비발산 PDE 와 완전 비선형 HJB 방정식 모두를 효과적으로 풀 수 있음을 보여줍니다. 인덱스 집합이 단일원소일 때 HJB 기법은 자연스럽게 선형 기법으로 퇴화합니다.

B. 영역 및 계수 가정의 완화

• 영역: 분석은 $C^{1,1}$ 또는 엄격하게 볼록한 영역에서의 최적 수렴 결과를 Lipschitz 다면체 영역으로 확장합니다.
  • $d=3$의 경우, $\Omega$는 볼록해야 합니다.
  • $d=2$의 경우, $\Omega$는 비볼록일 수 있습니다 ($p$에 대한 제한이 있음).
• 계수: 이 방법은 가산 조밀 부분집합을 포함하는 특정 근사 조건 (가정 6) 을 만족하는 한, HJB 방정식에서 불연속 계수를 허용합니다. 이는 많은 기존 DG 방법에서 요구되는 연속성 가정 ($A \in C^0$) 을 완화합니다.

C. $W^{2,p}$에서의 최적 수렴

이 논문은 광범위한 $p$ 범위에 대해 이산 $W^{2,p}_h$-노름에서 최적 수렴 속도를 증명합니다:

• 볼록 다면체 ($d=2,3$): $p \in (1, 2]$.
• 비볼록 다각형 ($d=2$): $p \in (4/3 - \epsilon_1, 4/3 + \epsilon_2)$.
  이는 $L^2$ 또는 $H^2$ ($p=2$) 에서만 수렴을 보장하는 기존 방법들에 비해 중요한 개선입니다.

D. 새로운 안정성 증명

연속 문제의 잘-제시됨에 의존하는 특수한 쌍대성 논증에 의존한 이전 연구들 (예: Feng et al.) 과 달리, 이 논문은 이산 컴팩트성과 강해의 유일성을 사용하여 안정성을 확립합니다. 이 접근법은 비볼록 영역과 불연속 계수에 대해 더 견고합니다.

4. 주요 결과

1. 잘-제시됨 (Theorems 3.3, 4.1):

   • 지정된 $A, b, c$ 및 영역 조건 하에서 선형 PDE 와 HJB 방정식 모두에 대해 $W^{2,p}(\Omega) \cap W^{1,p}_0(\Omega)$ 내의 강해의 존재성과 유일성이 증명됩니다.
   • HJB 방정식의 경우, Cordes 조건이 가산 부분집합 가정을 통해 불연속 계수를 허용하도록 일반화됩니다.

2. 안정성 (Theorems 3.2, 4.3):

   • 이산 해 $u_h$는 안정성 추정을 만족합니다: $\|u_h\|_{W^{2,p}_h} \leq C \|f\|_{L^p}$.
   • 이 안정성은 영역 기하학에 의존하는 특정 $p$ 범위에서 성립합니다.

3. 최적 수렴 (Theorems 3.3, 4.3):

   • 오차 추정이 유도됩니다:
     $$\|u - u_h\|_{W^{1,p}} + \|u - u_h\|_{W^{2,p}_h} \leq C \left( \inf_{\chi_h \in V_h} (\|u - \chi_h\|_{W^{1,p}} + \|u - \chi_h\|_{W^{2,p}_h}) + \|P_h(D^2u) - D^2u\|_{L^p} \right)$$
   • 이는 정확한 해가 충분히 매끄럽다면, 이 방법이 이산 $W^{2,p}$ 노름에서 최적 수렴 속도를 달성함을 의미합니다.

5. 의의

• 이론과 실무의 연결: 이 논문은 비발산 형태 PDE 를 푸는 데 표준 유한요소 공간 ($C^0$ 연속성) 을 사용하는 것에 대한 엄밀한 이론적 기반을 제공합니다. 이는 이산화하기 notoriously 어려운 문제입니다.
• 비볼록성 처리: 분석을 비볼록 2 차원 다각형으로 확장함으로써, 이 방법은 영역이 거의 완벽하게 볼록하지 않은 실제 공학 및 물리 문제의 더 넓은 범위에 적용 가능해집니다.
• HJB 응용: 계수 연속성 가정의 완화는 계수가 불연속일 수 있는 확률적 제어 문제 (예: 전환 체제) 에 이 방법을 적합하게 만들며, 이전의 $C^0$-정합 방법들이 실패했던 시나리오를 다룹니다.
• 통합 접근법: 하나의 수치 기법이 선형 및 완전 비선형 경우 모두를 처리함을 보여줌으로써 복잡한 시스템에 대한 구현 및 이론적 분석을 단순화합니다.

요약하자면, Weifeng Qiu 의 연구는 완화된 정칙성 요구사항을 가진 일반 다면체 영역에서 강노름으로 최적 수렴을 달성하는 견고하고 통합된 유한요소 프레임워크를 수립함으로써 비발산 형태 방정식에 대한 수치 PDE 의 최첨단 기술을 발전시켰습니다.