Rates of convergence of finite element approximations of second-order mean field games with nondifferentiable Hamiltonians

1. 연구 배경 및 문제 정의 (Problem Definition)

• 주제: 평균장 게임 (Mean Field Games, MFG) 시스템의 수치 해법, 특히 **미분 불가능한 해밀토니안 (Nondifferentiable Hamiltonians)**을 포함하는 정적 (Stationary) 2 차 MFG 시스템에 대한 유한 요소법 (FEM) 의 오차 분석입니다.
• 수학적 모델:
  • 상태 공간 $\Omega$에서 플레이어의 분포 밀도 $m$과 가치 함수 (Value function) $u$를 구하는 연립 편미분 방정식 (PDE) 시스템입니다.
  • HJB 방정식: $-\nu\Delta u + H(x, \nabla u) = F[m](x)$
  • KFP 방정식 (포함 관계): $-\nu\Delta m - G(x) \in \text{div}(m \partial_p H(x, \nabla u))$
  • 여기서 $H$는 해밀토니안이며, $\partial_p H$는 $p$에 대한 $H$의 **부분 미분 (Subdifferential)**입니다.
• 핵심 난제:
  • 기존 MFG 연구들은 대부분 해밀토니안 $H$가 미분 가능하다고 가정하여 KFP 방정식의 대류 항 (advective term) 을 명확하게 정의했습니다.
  • 그러나 실제 응용 (최소 탈출 시간 문제, 뱅크 - 뱅크 제어 등) 에서는 최적 피드백 제어가 유일하지 않아 $H$가 미분 불가능할 수 있으며, 이때 $\partial_p H$는 집합 값 함수가 됩니다.
  • 이로 인해 KFP 방정식은 편미분 포함 (Partial Differential Inclusion, PDI) 형태로 변하며, 수치 해법에서 선택 가능한 측정 가능 함수 (measurable selections) 간의 정량적 제어 (quantitative control) 가 어렵고, 밀도 함수 $m$의 정칙성 (regularity) 이 매우 낮아질 수 있습니다.

2. 방법론 (Methodology)

저자들은 **안정화된 유한 요소법 (Stabilized Finite Element Methods)**을 사용하여 위 문제를 해결하고 수렴 속도를 증명하기 위해 다음과 같은 접근법을 취했습니다.

• 규칙화 (Regularization) 전략:
  • 미분 불가능한 해밀토니안 $H$를 Moreau-Yosida 규칙화를 통해 $C^{1,1}$ 클래스의 미분 가능한 함수 $H_\lambda$로 근사합니다.
  • 이를 통해 원래의 PDI 시스템 (1.1) 을 규칙화된 PDE 시스템 (1.3) 으로 변환하여 수치 해석을 용이하게 합니다.
• 이중 오차 분석 (Dual Error Analysis):
  • 전체 오차를 두 단계로 나누어 분석합니다:
    1. 규칙화 오차: 원래 해 $(u, m)$과 규칙화된 해 $(u_\lambda, m_\lambda)$ 사이의 차이.
    2. 이산화 오차: 규칙화된 연속 해 $(u_\lambda, m_\lambda)$와 규칙화된 FEM 근사해 $(u_{k,\lambda}, m_{k,\lambda})$ 사이의 차이.
• 안정화 기법:
  • **Discrete Maximum Principle (DMP)**을 만족하도록 설계된 안정화 항 (Stabilization term, $D_k$) 을 도입합니다. 이는 밀도 함수 $m$의 비음성 (non-negativity) 을 보장하고, 해밀토니안의 미분 불가능성으로 인한 수치적 불안정성을 제어합니다.
  • 반-매끄러움 (Semismoothness) 성질: 미분 불가능한 해밀토니안 $H$에 대한 Lemma 2.1 (반-매끄러움 결과) 을 활용하여, 해밀토니안의 비선형 결합을 선형화하고 오차 항을 제어합니다.
• 이중성 논증 (Duality Argument):
  • 밀도 함수 $m$의 $L^2$ 오차 한계를 얻기 위해 Aubin-Nitsche 이중성 논증을 변형하여 적용합니다. 이는 밀도 함수의 낮은 정칙성 ($H^1$만 보장됨) 을 극복하기 위해 필수적입니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

• 수렴 속도 증명 (Theorem 5.1):

  • 일반 유계 다면체 Lipschitz 영역에서, FEM 근사해 $(u_k, m_k)$와 정확한 해 $(u, m)$ 사이의 오차에 대한 **사전 오차 한계 (A priori error bound)**를 증명했습니다.
  • 결과:
    $$\|u - u_k\|_{H^1(\Omega)} + \|m - m_k\|_{L^2(\Omega)} \lesssim h_k^{\gamma/3}$$
    • 여기서 $h_k$는 메쉬 크기, $\gamma \in (0, 1]$는 타원형 정칙성 (elliptic regularity) 지수입니다.
    • 볼록 영역 (Convex domain) 인 경우 $\gamma=1$이 되어 $O(h^{1/3})$의 수렴 속도를 얻습니다.
  • 의미: 이는 미분 불가능한 해밀토니안을 가진 MFG 시스템에 대해 유한 요소법의 수렴 속도를 정량적으로 규명한 최초의 결과 중 하나입니다.

• 규칙화 및 이산화의 결합 (Corollary 5.2):

  • 실제 계산에서는 미분 불가능한 포함 관계를 직접 푸는 것보다 규칙화된 PDE 시스템을 푸는 것이 더 쉽습니다.
  • 저자들은 규칙화된 FEM 해 $(u_{k, \lambda_k}, m_{k, \lambda_k})$가 원래 해에 대해 다음과 같은 오차 한계를 가진다는 것을 보였습니다:
    $$\|u - u_{k, \lambda_k}\|_{H^1} + \|m - m_{k, \lambda_k}\|_{L^2} \lesssim \lambda_k^{1/2} + \lambda_k^{-1} h_k^\gamma$$
  • 최적의 $\lambda_k \sim h_k^{2\gamma/3}$를 선택하면 위와 동일한 $O(h^{\gamma/3})$의 수렴 속도를 얻으며, 이는 실제 계산에서 규칙화 기법의 유효성을 입증합니다.

• 밀도 함수의 $L^2$ 수렴:

  • 밀도 함수의 기울기 ($\nabla m$) 에 대한 $H^1$ 수렴은 부분 미분 집합의 불연속성으로 인해 여전히 열린 문제이나, 본 논문은 밀도 함수 자체에 대한 $L^2$ 수렴 속도를 확립했습니다.

4. 의의 및 중요성 (Significance)

1. 이론적 확장: 기존에 미분 가능한 해밀토니안에만 국한되었던 MFG 수치 해석 이론을 **비미분 가능 (Nondifferentiable)**한 경우로 확장했습니다. 이는 최적 제어 이론에서 발생하는 실제적인 문제 (예: 뱅크 - 뱅크 제어) 를 MFG 프레임워크 내에서 엄밀하게 다룰 수 있는 기반을 마련했습니다.
2. 정량적 오차 분석: 기존 연구들이 단순히 수렴성 (Convergence) 만을 보였던 것과 달리, 구체적인 **오차 상한 (Error bounds)**과 **수렴 속도 (Rates of convergence)**를 제시했습니다.
3. 실용적 적용 가능성: 규칙화 (Regularization) 기법이 수치 해법에서 계산의 용이성을 제공할 뿐만 아니라, 이론적으로도 통제 가능한 오차 범위 내에서 정확한 해를 근사할 수 있음을 보였습니다. 이는 실제 공학 및 경제학 모델링에서 MFG 시스템을 구현할 때 중요한 지침이 됩니다.
4. 수치 기법의 견고성: 안정화 기법 (Stabilization) 과 이산 최대 원리 (DMP) 를 결합하여 밀도 함수의 물리적 성질 (비음성) 을 보존하면서도 수렴성을 증명하는 강력한 프레임워크를 제시했습니다.

결론

이 논문은 미분 불가능한 해밀토니안을 갖는 평균장 게임 시스템에 대해, Moreau-Yosida 규칙화와 안정화된 유한 요소법을 결합하여 $O(h^{1/3})$ (볼록 영역 기준) 의 수렴 속도를 가진 수치 해법을 제안하고 엄밀하게 증명했습니다. 이는 비선형 편미분 포함 문제를 다루는 수치 해석 분야에서 중요한 이론적 진전이며, 실제 응용 문제에서의 수치 시뮬레이션 신뢰도를 높이는 데 기여합니다.