Gradient constraints, Born-Infeld, and maximal surfaces via superposition of infinitely many $p$-Laplacians

We study the Dirichlet problem for an infinite series of $p$-Laplacians, $-{\sum }_{p=2}^{\infty }{a}_p{\Delta }_{p}u=fin\Omega ,u=gon\partial \Omega ,$ where $\{a_p\}$ is a sequence of nonnegative numbers whose power series has radius of convergence $\sigma \in (0,\infty ]$. The operator is formally $-\mathrm{div}(A(|\nabla u|)\nabla u)$ with $A$ singular at $|\nabla u|=\sigma$; model cases are the mean curvature operator in Minkowski space and the Born-Infeld operator of nonlinear electrostatics. The radius of convergence forces the gradient constraint $\Vert \nabla u{\Vert }_{\infty }\le \sigma$, so the natural variational problem is constrained. A unique minimizer of the associated energy always exists (for boundary data compatible with the constraint) and always solves a variational inequality, and we identify the saturation flux $\Lambda :={\sum }_{p\ge 2}{a}_p{\sigma }^{p-1}$ as the quantity governing solvability of the equation: if $\Lambda <\infty$ a weak solution exists only if $|{\int }_{E}f|\le \Lambda P(E)$ for every set of finite perimeter $E⋐\Omega$, so for $f\equiv \lambda$ no solution exists once $\lambda >\Lambda h(\Omega )$, $h(\Omega )$ the Cheeger constant of $\Omega$. The threshold is sharp on balls, where the minimizer has a saturation region $\{|\nabla u|=\sigma \}$ of positive measure. When $\Lambda =\infty$, as for Born-Infeld type operators, no such obstruction is present, and the minimizer solves the equation whenever a Lipschitz bound below $\sigma$ is available; for $f\equiv 0$ we obtain such a bound, uniform in the truncations of the series, under a bounded slope condition on the boundary datum. As applications we obtain an isoperimetric bound on the charge densities supported by a nonlinear electrostatics with saturating displacement, and a quantitative convergence rate for the weak field expansion of the Born-Infeld model.