Finite element approximation of time-dependent mean field games with nondifferentiable Hamiltonians
Yohance A. P. Osborne, Iain Smears
TL;DR
This work develops a time-dependent mean field game framework with nondifferentiable Hamiltonians by formulating the system as a partial differential inclusion using the Hamiltonian subdifferential. It proposes a monotone finite element method with mass lumping and a novel volume-based stabilization that enforces a discrete maximum principle, enabling stable and convergent space-time discretizations without restrictive time-step constraints. The authors prove existence and uniqueness of weak solutions for the MFG-PDI and establish convergence of the discrete solutions to the continuous weak solution, with strong convergence for the value function and density under standard monotonicity assumptions. A numerical experiment corroborates the theoretical results, displaying expected convergence rates and strong transport-field convergence, highlighting practical applicability to MFGs with nondifferentiable Hamiltonians.
Abstract
The standard formulation of the PDE system of Mean Field Games (MFG) requires the differentiability of the Hamiltonian. However in many cases, the structure of the underlying optimal problem leads to a convex but nondifferentiable Hamiltonian. For time-dependent MFG systems, we introduce a generalization of the problem as a Partial Differential Inclusion (PDI) by interpreting the derivative of the Hamiltonian in terms of the subdifferential set. In particular, we prove the existence and uniqueness of weak solutions to the resulting MFG PDI system under standard assumptions in the literature. We propose a monotone stabilized finite element discretization of the problem, using conforming affine elements in space and an implicit Euler discretization in time with mass-lumping. We prove the strong convergence in $L^2(H^1)$ of the value function approximations, and strong convergence in $L^p(L^2)$ of the density function approximations, together with strong $L^2$-convergence of the value function approximations at the initial time.
