Lipschitz Regularity of Viscosity Solutions to the Infinity Laplace Equation

Xiao Han, Fang Liu^*

School of Mathematics and Statistics, Nanjing University of Science and Technology, Nanjing, China.

DOI: 10.4236/jamp.2023.1110197                

Abstract

In this paper, we study the viscosity solutions of the Neumann problem

in a bounded C² domain Ω, where Δ^N_∞ is called the normalized infinity Laplacian. The normalized infinity Laplacian was first studied by Peres, Shramm, Sheffield and Wilson from the point of randomized theory named tug-of-war, which has wide applications in optimal mass transportation, financial option price problems, digital image processing, physical engineering, etc. We give the Lipschitz regularity of the viscosity solutions of the Neumann problem. The method we adopt is to choose suitable auxiliary functions as barrier functions and combine the perturbation method and viscosity solutions theory.

Keywords

Normalized Infinity Laplacian, Viscosity Solution, Lipschitz Regularity

Keywords:

1. Introduction

In this paper, we study the Lipschitz regularity of the viscosity solutions of the Neumann problem

$\{\begin{array}{ll}{\Delta }_{\infty }^{N}u+\beta \left|Du\right|+\xi \left(x\right)\cdot Du+\eta \left(x\right)u=g\left(x\right)\hfill & \text{in}\Omega \mathrm{,}\hfill \\ \frac{\partial u}{\partial \stackrel{\to }{n}}=0\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$ (1)

where Ω is a bounded $C^2$ domain in ${ℝ}^{\mathbb{N}}$ , ${ℝ}^{\mathbb{N}}$ denotes a set of $\mathbb{N}$ -dimensional Euclidean space, $\beta \in ℝ$ , $\partial \Omega$ denotes the boundary of Ω, $\stackrel{\to }{n}\left(x\right)$ is the unit exterior normal to the domain Ω at $x\in \partial \Omega$ , $\xi \left(x\right)\mathrm{<mo>\; :\; </mo>}\Omega \to {ℝ}^{\mathbb{N}}$ , $\eta \left(x\right)\mathrm{<mo>\; :\; </mo>}\Omega \to ℝ$ are continuous in $\bar{\Omega }$ (the closure of Ω), $g\left(x\right)$ is a bounded function in $\bar{\Omega }$ and

${\Delta }_{\infty }^{N}u=\langle D^{2}u\frac{Du}{\left|Du\right|}\mathrm{,}\frac{Du}{\left|Du\right|}\rangle \mathrm{,}$ (2)

is called the normalized infinity Laplacian.

The infinity Laplace equation

${\Delta }_{\infty }u={\displaystyle \underset{i=1}{\overset{n}{\sum }}}{\displaystyle \underset{j=1}{\overset{n}{\sum }}}\frac{\partial u}{\partial x_i}\frac{{\partial }^{2}u}{\partial x_i\partial x_j}\frac{\partial u}{\partial x_j}=0,$

is the Euler-Lagrange equation associated with $L^{\infty }$ -variational problem related to the absolutely minimizing Lipschitz extensions, which was first studied by Aronsson [1] [2] [3] [4] . The infinity Laplacian has attracted more and more attention because it is highly degenerate and has no variational structure. It has been widely used in the Monge-Kantorovich mass transfer problem in [5] , digital image processing in [6] [7] and financial mathematics in [8] .

${\Delta }_{\infty }^{N}u+\xi \left(x\right)\cdot Du$ is called the infinity Laplacian with a transport term related to tug-of-war. López, Navarro and Rossi [9] gave an explanation from the point of tug-of-war game. Let us briefly recall the game: let F be the final payoff function defined in a narrow strip around the boundary $\partial \Omega$ . The tug-of-war game with a transport term is played with two stages. First the players toss an unfair coin, which has head probability $0<C\left(\epsilon \right)<1$ and tail probability $1-C\left(\epsilon \right)$ . If the players have obtained a head, then they toss a new (fair) coin and the winner moves the token to any new position $x^1\in {\bar{B}}_{\epsilon }\left(x^0\right)$ . But if in the first (unfair) coin toss, the players obtain a tail, the token is moved to $x^0+\xi \left(x^0\right)\epsilon$ . Note that there is no strategies of the players involved if they get a tail in the first coin toss. The game continues until the first time the token arrives to $x^{\varsigma }\in {ℝ}^{\mathbb{N}}\backslash \Omega$ and then Player I earns $F\left(x^{\varsigma }\right)$ , Player II earns $-F\left(x^{\varsigma }\right)$ , where F is the extension of $f$ from $\partial \Omega$ to a small strip ${\Gamma }_{\epsilon }=\left\{x\in {ℝ}^{\mathbb{N}}\backslash \Omega \mathrm{<mo>\; :\; </mo>}\text{dist}\left(x\mathrm{,}\partial \Omega \right)<\epsilon {\Vert \xi \Vert }_{\infty }\right\}$ and gives the final payoff of the game. López, Navarro and Rossi found a viscosity solution to

$\{\begin{array}{ll}-{\Delta }_{\infty }^{N}u-\langle Du\mathrm{,}\xi \rangle =0\hfill & \text{in}\Omega \mathrm{,}\hfill \\ u=f\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$ (3)

where $f$ is a Lipschitz continuous function. They obtained the existence and uniqueness of a viscosity solution by a L^p-approximation procedure when ξ is a continous gradient vector field. They also proved the stability of the unique solution with respect to ξ. In addition, when ξ is Lipschitz continuous but not necessarily a gradient, they proved that the problem (3) has a viscosity solution. Some kinds of modified tug-of-war have received a lot of attention, such as [10] - [19] .

${\Delta }_{\infty }^{N}u+\beta \left|Du\right|$ is called the β-biased infinity Laplacian, which was first introduced by Peres, Pete and Somersille when modelling the stochastic game named biased tug-of-war in [17] . They investigated the random game with a final payoff function and a running payoff function. It’s a zero sum game with two players in which the earnings of one of them are the losses of the other. Armstrong, Smart and Somersille [20] studied the mixed Dirichlet-Neumann boundary value problem

$\{\begin{array}{ll}-{\Delta }_{\infty }^{N}u-\beta \left|Du\right|=g\hfill & \text{in}\Omega \mathrm{,}\hfill \\ u=f\hfill & \text{on}{\Gamma }_D\mathrm{,}\hfill \\ \frac{\partial u}{\partial \stackrel{\to }{n}}=0\hfill & \text{on}{\Gamma }_N\mathrm{,}\hfill \end{array}$

where ${\Gamma }_D\cup {\Gamma }_N=\partial \Omega$ is a partition of $\partial \Omega$ with ${\Gamma }_D$ nonempty and closed. They obtained existence, uniqueness and stability results for the boundary-value problem. Liu and Yang [21] established the existence of the principal Dirichlet eigenvalue based on the comparison principle. They also established the Harnack inequality and the Lipschitz regularity of a nonnegative viscosity supersolution to the β-biased equation

${\Delta }_{\infty }^{N}u+\beta \left|Du\right|+\lambda a\left(x\right)u=g\left(x\right)\mathrm{,}$

where $\lambda \in ℝ$ and $\lambda \ge 0$ , the weight function $a\left(x\right)$ is positive in $\bar{\Omega }$ and $a\left(x\right)\in C\left(\Omega \right)\cap L^{\infty }\left(\Omega \right)$ . The key of their method is to choose suitable exponential cones as barrier functions.

For the case $\beta =0$ , Lu and Wang [22] [23] studied the inhomogeneous infinity Laplace equation

${\Delta }_{\infty }^{N}u=g\text{in}\Omega \mathrm{.}$

They showed existence and uniqueness of the viscosity solutions of the Dirichlet problem under the intrinsic condition that g does not change its sign from the PDE's methods. Patrizi [24] studied the following Neumann problem

$\{\begin{array}{ll}{\Delta }_{\infty }^{N}u+\xi \left(x\right)\cdot Du+\eta \left(x\right)u=0\hfill & \text{in}\Omega \mathrm{,}\hfill \\ \frac{\partial u}{\partial \stackrel{\to }{n}}=0\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$

and showed the Lipschitz regularity in the whole $\bar{\Omega }$ of the viscosity solutions and obtained the existence of the principal eigenvalue.

Aronsson [25] obtained the specific form of a viscosity solution to the infinity Laplace equation ( ${\Delta }_{\infty }u=0$ ) in two-dimensional space: $u\left(x_1\mathrm{,}{x}_2\right)=x_1^{\frac{4}{3}}-x_2^{\frac{4}{3}}$ . Thus, the regularity of infinity harmonic functions (viscosity solutions to ${\Delta }_{\infty }u=0$ ) is at most $C^{\mathrm{1,}\frac{1}{3}}$ . In [26] , the $C^1$ regularity of infinity harmonic functions was proved by Savin in dimension two. Later, Evans and Savin [27] established the $C^{\mathrm{1,}\alpha }$ regularity of infinity harmonic functions for some $\alpha >0$ in dimension two. For $\mathbb{N}\ge 3$ , Evans and Smart [28] [29] proved that infinity harmonic functions in ${ℝ}^{\mathbb{N}}$ are differentiable everywhere.

In this paper, we study the Lipschitz regularity of viscosity solutions of the Neumann problem (1). The main result can be summarized as the following theorem.

Theorem 1 Assume that Ω is a bounded domain of class $C^2$ , $\beta \in ℝ$ , $\xi \left(x\right)\mathrm{<mo>\; :\; </mo>}\Omega \to {ℝ}^{\mathbb{N}}$ , $\eta \left(x\right)\mathrm{<mo>\; :\; </mo>}\Omega \to ℝ$ are continuous in $\bar{\Omega }$ , $g$ is a bounded function in $\bar{\Omega }$ . If $u\in C\left(\bar{\Omega }\right)$ is a viscosity solution of

$\{\begin{array}{ll}{\Delta }_{\infty }^{N}u+\beta \left|Du\right|+\xi \left(x\right)\cdot Du+\eta \left(x\right)u=g\left(x\right)\hfill & \text{in}\Omega \mathrm{,}\hfill \\ \frac{\partial u}{\partial \stackrel{\to }{n}}=0\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$ (4)

then there exists a constant $C_0$ depending on $\Omega \mathrm{,}\beta \mathrm{,}{\Vert g\Vert }_{\infty }$ and ${\Vert u\Vert }_{\infty }$ such that

$\left|u\left(x\right)-u\left(y\right)\right|\le C_0\left|x-y\right|\mathrm{,}\forall x\mathrm{,}y\in \bar{\Omega }\mathrm{.}$ (5)

2. Definitions of the Viscosity Solutions

Since the normalized infinity Laplacian ${\Delta }_{\infty }^N$ is singular at the points where the gradient vanishes, we give a proper explanation to the operator by the viscosity solutions theory according to Crandall, Ishii and Lions [30] .

We denote $S\left(\mathbb{N}\right)$ as the set of symmetric matrices on ${ℝ}^{\mathbb{N}\times \mathbb{N}}$ and define $\Vert X\Vert$ in ${ℝ}^{\mathbb{N}\times \mathbb{N}}$ by letting $\Vert X\Vert =\sup \left\{\left|X\theta \right|\mathrm{<mo>\; |\; </mo>}\theta \in {ℝ}^{\mathbb{N}}\mathrm{,}\left|\theta \right|\le 1\right\}$ .

Denote $\sigma \mathrm{<mo>\; :\; </mo>}{ℝ}^{\mathbb{N}}\to S\left(\mathbb{N}\right)$ :

$\sigma \left(p\right)\mathrm{<mo>\; :\; </mo>}=\frac{p\otimes p}{{\left|p\right|}^2}\mathrm{,}$

where $\otimes$ denotes the tensor product.

Then we get

${\Delta }_{\infty }^{N}u=\text{tr}\left(\sigma \left(Du\right)D^{2}u\right)\mathrm{,}$

for any $u\in C^2\left(\Omega \right)$ .

It is easy to check that the following properties are valid.

(1) $\sigma \left(p\right)$ is homogeneous of order 0, i.e., for any $a\in ℝ$ and $p\in {ℝ}^{\mathbb{N}}$ , one has

$\sigma \left(ap\right)=\sigma \left(p\right)\mathrm{.}$

(2) For all $p\in {ℝ}^{\mathbb{N}}$ , one has

$0\le \sigma \left(p\right)\le I_{\mathbb{N}}\mathrm{,}$

where $I_{\mathbb{N}}$ denotes the identity matrix in ${ℝ}^{\mathbb{N}\times \mathbb{N}}$ .

(3) $\sigma \left(p\right)$ is idempotent, i.e.,

${\left(\sigma \left(p\right)\right)}^2=\sigma \left(p\right)\mathrm{.}$

Suppose that Ω is a bounded $C^2$ domain. Obviously, one has the interior sphere condition and the uniform exterior sphere condition, i.e.,

(Ω1) For $\forall x\in \partial \Omega$ , there exist $R>0$ and $y\in \Omega$ for which $\left|x-y\right|=R$ and $B_R\left(y\right)\subset \Omega$ .

(Ω2) For $\forall x\in \partial \Omega$ , there exists $r>0$ such that $B_r\left(x+r\stackrel{\to }{n}\left(x\right)\right)\cap \Omega =\varnothing$ .

From (Ω2), one has

$\langle y-x\mathrm{,}\stackrel{\to }{n}\left(x\right)\rangle \le \frac{1}{2r}{\left|y-x\right|}^2\mathrm{,}x\in \partial \Omega \mathrm{,}y\in \bar{\Omega }\mathrm{.}$ (6)

Due to the $C^2$ -regularity of Ω, we obtain the existence of a neighborhood of $\partial \Omega$ in $\bar{\Omega }$ on which the distance to the boundary

$d\left(x\right)\mathrm{<mo>\; :\; </mo>}=\inf \left\{\left|x-y\right|\mathrm{,}y\in \partial \Omega \right\}\mathrm{,}x\in \bar{\Omega }$

is of class $C^2$ . Without loss of generality, we assume that $\left|Dd\left(x\right)\right|\le 1$ on $\bar{\Omega }$ .

The $\text{USC}\left(\bar{\Omega }\right)$ denotes the set of upper semicontinuous functions on $\bar{\Omega }$ and the $\text{LSC}\left(\bar{\Omega }\right)$ denotes the set of lower semicontinuous functions on $\bar{\Omega }$ . We define $B\mathrm{<mo>\; :\; </mo>}\partial \Omega \times ℝ\times {ℝ}^{\mathbb{N}}\to ℝ$ . Now we give the definitions of the viscosity solutions of the Neumann problem according to [30] [31] .

Definition 1 Any function $u\in \text{USC}\left(\bar{\Omega }\right)$ (resp., $u\in \text{LSC}\left(\bar{\Omega }\right)$ ) is called a viscosity subsolution (resp., viscosity supersolution) of

$\{\begin{array}{ll}{\Delta }_{\infty }^{N}u+\beta \left|Du\right|+\xi \left(x\right)\cdot Du+\eta \left(x\right)u=g\left(x\right)\hfill & \text{in}\Omega \mathrm{,}\hfill \\ B\left(x\mathrm{,}u\mathrm{,}Du\right)=0\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$ (7)

if the following conditions hold:

(1) For every $x_0\in \Omega$ , for all $\phi \in C^2\left(\bar{\Omega }\right)$ , such that $u-\phi$ has a local strict maximum (resp., strict minimum) at $x_0$ with $u\left(x_0\right)=\phi \left(x_0\right)$ and $D\phi \left(x_0\right)\ne 0$ , one has

${\Delta }_{\infty }^N\phi \left(x_0\right)+\beta \left|D\phi \left(x_0\right)\right|+\xi \left(x_0\right)\cdot D\phi \left(x_0\right)+\eta \left(x_0\right)u\left(x_0\right)\ge \left(\text{resp}\mathrm{.,}\le \right)g\left(x_0\right)\mathrm{.}$

If $u\equiv k$ ( $k$ is a constant) in a neighborhood of $x_0$ , then

$\eta \left(x_0\right)k\ge \left(\text{resp}\mathrm{.,}\le \right)g\left(x_0\right)\mathrm{.}$

(2) For every $x_0\in \partial \Omega$ , for all $\phi \in C^2\left(\bar{\Omega }\right)$ , such that $u-\phi$ has a local maximum (resp., minimum) at $x_0$ with $u\left(x_0\right)=\phi \left(x_0\right)$ and $D\phi \left(x_0\right)\ne 0$ , one has

$\begin{array}{l}\min \{-{\Delta }_{\infty }^N\phi \left(x_0\right)-\beta \left|D\phi \left(x_0\right)\right|-\xi \left(x_0\right)\cdot D\phi \left(x_0\right)-\eta \left(x_0\right)u\left(x_0\right)\\ +g\left(x_0\right),B\left(x_0,u\left(x_0\right),D\phi \left(x_0\right)\right)\}\le 0.\end{array}$

(resp.,

$\begin{array}{l}\max \{-{\Delta }_{\infty }^N\phi \left(x_0\right)-\beta \left|D\phi \left(x_0\right)\right|-\xi \left(x_0\right)\cdot D\phi \left(x_0\right)-\eta \left(x_0\right)u\left(x_0\right)\\ +g\left(x_0\right),B\left(x_0,u\left(x_0\right),D\phi \left(x_0\right)\right)\}\ge 0.\end{array}$ )

If $u\equiv k$ ( $k$ is a constant) in a neighborhood of $x_0$ in $\bar{\Omega }$ , then

$\min \left\{-\eta \left(x_0\right)k+g\left(x_0\right)\mathrm{,}B\left(x_0\mathrm{,}k\mathrm{,0}\right)\right\}\le 0.$

(resp.,

$\max \left\{-\eta \left(x_0\right)k+g\left(x_0\right)\mathrm{,}B\left(x_0\mathrm{,}k\mathrm{,0}\right)\right\}\ge 0.$ )

We call that $u$ is a viscosity solution if $u$ is both a viscosity supersolution and a viscosity subsolution.

The definition of the viscosity solutions can be also given by semijets ${\bar{J}}^{\mathrm{2,}+}u\left(x_0\right)$ and ${\bar{J}}^{\mathrm{2,}-}u\left(x_0\right)$ according to [32] .

Definition 2 The second-order superjet of $u$ at $x_0$ is defined to be the set

$J^{2,+}u\left(x_0\right)=\left\{\left(D\phi \left(x\right),D^2\phi \left(x\right)\right):\phi \text{is}{C}^2\text{and}u-\phi \text{has}\text{a}\text{local}\text{maximum}\text{at}{x}_0\right\},$

whose closure is defined as

$\begin{array}{l}{\bar{J}}^{2,+}u\left(x_0\right)=\{\left(p,X\right)\in {ℝ}^{\mathbb{N}}\times S\left(\mathbb{N}\right):\exists \left(x_n,p_n,X_n\right)\in \Omega \times {ℝ}^{\mathbb{N}}\times S\left(\mathbb{N}\right)\\ \text{such}\text{that}\left(p_n,X_n\right)\in J^{2,+}u\left(x_n\right)\text{and}\\ \left(x_n,u\left(x_n\right),p_n,X_n\right)\to \left(x_0,u\left(x_0\right),p,X\right)\},\end{array}$

and the second-order subjet of $u$ at $x_0$ is defined to be the set

$J^{2,-}u\left(x_0\right)=\left\{\left(D\phi \left(x\right),D^2\phi \left(x\right)\right):\phi \text{is}{C}^2\text{and}u-\phi \text{has}\text{a}\text{local}\text{minimum}\text{at}{x}_0\right\},$

whose closure is defined as

$\begin{array}{l}{\bar{J}}^{2,-}u\left(x_0\right)=\{\left(p,X\right)\in {ℝ}^{\mathbb{N}}\times S\left(\mathbb{N}\right):\exists \left(x_n,p_n,X_n\right)\in \Omega \times {ℝ}^{\mathbb{N}}\times S\left(\mathbb{N}\right)\\ \text{such}\text{that}\left(p_n,X_n\right)\in J^{2,-}u\left(x_n\right)\text{and}\\ \left(x_n,u\left(x_n\right),p_n,X_n\right)\to \left(x_0,u\left(x_0\right),p,X\right)\}.\end{array}$

Next we give the definitions of viscosity solutions by semijets.

Definition 3 Any function $u\in \text{USC}\left(\bar{\Omega }\right)$ (resp., $u\in \text{LSC}\left(\bar{\Omega }\right)$ ) is called a viscosity subsolution (resp., viscosity supersolution) of

$\{\begin{array}{ll}{\Delta }_{\infty }^{N}u+\beta \left|Du\right|+\xi \left(x\right)\cdot Du+\eta \left(x\right)u=g\left(x\right)\hfill & \text{in}\Omega \mathrm{,}\hfill \\ B\left(x\mathrm{,}u\mathrm{,}Du\right)=0\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$

if the following conditions hold:

(1) For every $x_0\in \Omega$ , $\forall \left(p\mathrm{,}X\right)\in {\bar{J}}^{\mathrm{2,}+}u\left(x_0\right)$ (resp., $\left(p\mathrm{,}X\right)\in {\bar{J}}^{\mathrm{2,}-}u\left(x_0\right)$ ) and $p\ne 0$ , one has

$\frac{1}{{\left|p\right|}^2}\langle Xp,p\rangle +\beta \left|p\right|+\xi \left(x_0\right)\cdot p+\eta \left(x_0\right)u\left(x_0\right)\ge \left(\text{resp}.,\le \right)g\left(x_0\right).$

If $u\equiv k$ ( $k$ is a constant) in a neighborhood of $x_0$ , then

$\eta \left(x_0\right)k\ge \left(\text{resp}\text{.,}\le \right)g\left(x_0\right)\mathrm{.}$

(2) For every $x_0\in \partial \Omega$ , $\forall \left(p\mathrm{,}X\right)\in {\bar{J}}^{\mathrm{2,}+}u\left(x_0\right)$ (resp., $\left(p\mathrm{,}X\right)\in {\bar{J}}^{\mathrm{2,}-}u\left(x_0\right)$ ) and $p\ne 0$ , one has

$\min \left\{-\frac{1}{{\left|p\right|}^2}\langle Xp,p\rangle -\beta \left|p\right|-\xi \left(x_0\right)\cdot p-\eta \left(x_0\right)u\left(x_0\right)+g\left(x_0\right),B\left(x_0,u\left(x_0\right),p\right)\right\}\le 0.$

(resp.,

$\max \left\{-\frac{1}{{\left|p\right|}^2}\langle Xp,p\rangle -\beta \left|p\right|-\xi \left(x_0\right)\cdot p-\eta \left(x_0\right)u\left(x_0\right)+g\left(x_0\right),B\left(x_0,u\left(x_0\right),p\right)\right\}\ge 0.$ )

If $u\equiv k$ ( $k$ is a constant) in a neighborhood of $x_0$ in $\bar{\Omega }$ , then

$\min \left\{-\eta \left(x_0\right)k+g\left(x_0\right)\mathrm{,}B\left(x_0\mathrm{,}k\mathrm{,0}\right)\right\}\le 0.$

(resp.,

$\max \left\{-\eta \left(x_0\right)k+g\left(x_0\right)\mathrm{,}B\left(x_0\mathrm{,}k\mathrm{,0}\right)\right\}\ge 0.$ )

We call that $u$ is a viscosity solution if $u$ is both a viscosity supersolution and a viscosity subsolution.

3. Lipschitz Regularity of Viscosity Solutions

In this section, we show the Lipschitz regularity of the viscosity solutions of the Neumann problem (1).

Theorem 2 Assume that Ω is a bounded domain of class $C^2$ , $\beta \in ℝ$ , $\xi \left(x\right)\mathrm{<mo>\; :\; </mo>}\Omega \to {ℝ}^{\mathbb{N}}$ , $\eta \left(x\right)\mathrm{<mo>\; :\; </mo>}\Omega \to ℝ$ are continuous in $\bar{\Omega }$ , $g$ and $h$ are bounded functions in $\bar{\Omega }$ . Let $u\in \text{USC}\left(\bar{\Omega }\right)$ be a viscosity subsolution of

$\{\begin{array}{ll}{\Delta }_{\infty }^{N}u+\beta \left|Du\right|+\xi \left(x\right)\cdot Du+\eta \left(x\right)u=g\left(x\right)\hfill & \text{in}\Omega \mathrm{,}\hfill \\ \frac{\partial u}{\partial \stackrel{\to }{n}}=0\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$

and $v\in \text{LSC}\left(\bar{\Omega }\right)$ be a viscosity supersolution of

$\{\begin{array}{ll}{\Delta }_{\infty }^{N}v+\beta \left|Dv\right|+\xi \left(x\right)\cdot Dv+\eta \left(x\right)v=h\left(x\right)\hfill & \text{in}\Omega \mathrm{,}\hfill \\ \frac{\partial v}{\partial \stackrel{\to }{n}}=0\hfill & \text{on}\partial \Omega \mathrm{,}\hfill \end{array}$

with $u$ and $v$ bounded, or $v\ge 0$ and $v$ bounded. If $m={\max }_{\bar{\Omega }}\left(u-v\right)\ge 0$ , then there exists $C_0>0$ such that

$u\left(x\right)-v\left(y\right)\le m+C_0\left|x-y\right|\mathrm{,}\forall x\mathrm{,}y\in \bar{\Omega }\mathrm{,}$ (8)

where $C_0$ depends on $\Omega ,\beta ,\mathbb{N},{\Vert \xi \Vert }_{\infty }\mathrm{,}{\Vert \eta \Vert }_{\infty }\mathrm{,}{\Vert g\Vert }_{\infty }\mathrm{,}{\Vert h\Vert }_{\infty }\mathrm{,}{\Vert v\Vert }_{\infty }\mathrm{,}m$ and ${\Vert u\Vert }_{\infty }$ or ${\sup }_{\bar{\Omega }}u$ .

Proof. We set

$\Psi \left(x\right)=PQ\left|x\right|-P{\left(Q\left|x\right|\right)}^2\mathrm{,}$

and

$\psi \left(x\mathrm{,}y\right)=m+{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(x\right)+d\left(y\right)\right)}\Psi \left(x-y\right)\mathrm{,}$

where $M$ is a fixed constant, $P$ and $Q$ are two positive constants to be chosen later.

If $Q\left|x\right|\le \frac{1}{4}$ , then

$\Psi \left(x\right)\ge \frac{3}{4}PQ\left|x\right|\mathrm{.}$ (9)

Define

${\Delta }_Q\mathrm{<mo>\; :\; </mo>}=\left\{\left(x\mathrm{,}y\right)\in {ℝ}^{\mathbb{N}}\times {ℝ}^{\mathbb{N}}\mathrm{<mo>\; :\; </mo>}\left|x-y\right|\le \frac{1}{4Q}\right\}\mathrm{.}$

Fix $P$ such that

$\underset{{\bar{\Omega }}^2}{\max }\left(u\left(x\right)-v\left(y\right)\right)\le m+\frac{P}{8}{\text{e}}^{-2\left(\left|\beta \right|+1\right)M{d}_0}\mathrm{,}$ (10)

where $d_0={\max }_{x\in \bar{\Omega }}d\left(x\right)$ . If we take $Q$ large enough, there holds

$u\left(x\right)-v\left(y\right)-\psi \left(x\mathrm{,}y\right)\le \mathrm{0,}\left(x\mathrm{,}y\right)\in {\Delta }_Q\cap {\bar{\Omega }}^2\mathrm{.}$

Step 1. Suppose by contradiction that for each $Q$ there exists a point $\left(\bar{x}\mathrm{,}\bar{y}\right)\in {\Delta }_Q\cap {\bar{\Omega }}^2$ such that

$u\left(\bar{x}\right)-v\left(\bar{y}\right)-\psi \left(\bar{x},\bar{y}\right)=\underset{{\Delta }_Q\cap {\bar{\Omega }}^2}{\max }\left(u\left(x\right)-v\left(y\right)-\psi \left(x,y\right)\right)>0.$

Here we have dropped the dependence of $\bar{x}\mathrm{,}\bar{y}$ on $Q$ for simplicity of notations.

If $v\ge 0$ , we obtain that $\Psi \left(x-y\right)$ is non-negative in ${\Delta }_Q$ and $m\ge 0$ by the inequality (9). Then $u\left(\bar{x}\right)>0$ .

Clearly $\bar{x}\ne \bar{y}$ . For any $x\mathrm{,}y\in \bar{\Omega }$ with $\left|x-y\right|=\frac{1}{4Q}$ , we get

$\begin{array}{c}u\left(x\right)-v\left(y\right)\le m+\frac{P}{8}{\text{e}}^{-2\left(\left|\beta \right|+1\right)M{d}_0}\\ \le m+\frac{P}{2}{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(x\right)+d\left(y\right)\right)}Q\left|x-y\right|\\ \le \psi \left(x,y\right).\end{array}$

Thus, $\left(\bar{x},\bar{y}\right)\in \text{int}\left({\Delta }_Q\right)\cap {\bar{\Omega }}^2$ .

Next we compute the derivatives of $\psi$ at $\left(\bar{x},\bar{y}\right)$ ,

$\begin{array}{l}{D}_x\psi \left(\bar{x},\bar{y}\right)={\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}PQ\{-\left(\left|\beta \right|+1\right)M\left|\bar{x}-\bar{y}\right|\left(1-Q\left|\bar{x}-\bar{y}\right|\right)Dd\left(\bar{x}\right)\begin{array}{c}\\ \end{array}\\ +\left(1-2Q\left|\bar{x}-\bar{y}\right|\right)\frac{\bar{x}-\bar{y}}{\left|\bar{x}-\bar{y}\right|}\},\end{array}$

and

$\begin{array}{l}{D}_y\psi \left(\bar{x},\bar{y}\right)={\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}PQ\{-\left(\left|\beta \right|+1\right)M\left|\bar{x}-\bar{y}\right|\left(1-Q\left|\bar{x}-\bar{y}\right|\right)Dd\left(\bar{y}\right)\begin{array}{c}\\ \end{array}\\ -\left(1-2Q\left|\bar{x}-\bar{y}\right|\right)\frac{\bar{x}-\bar{y}}{\left|\bar{x}-\bar{y}\right|}\}.\end{array}$

For large $Q$ , one has

$0<{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}PQ\left(\frac{1}{2}-\left(\left|\beta \right|+1\right)M\left|\bar{x}-\bar{y}\right|\right)\le \left|D_x\psi \left(\bar{x},\bar{y}\right)\right|\le 2PQ,$ (11)

and

$0<{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}PQ\left(\frac{1}{2}-\left(\left|\beta \right|+1\right)M\left|\bar{x}-\bar{y}\right|\right)\le \left|D_y\psi \left(\bar{x},\bar{y}\right)\right|\le 2PQ.$ (12)

By the inequality (6), if $\bar{x}\in \partial \Omega$ , one has

$\begin{array}{l}\langle D_x\psi \left(\bar{x},\bar{y}\right),\stackrel{\to }{n}\left(x\right)\rangle \\ ={\text{e}}^{-\left(\left|\beta \right|+1\right)Md\left(\bar{y}\right)}PQ\{\left(\left|\beta \right|+1\right)M\left|\bar{x}-\bar{y}\right|\left(1-Q\left|\bar{x}-\bar{y}\right|\right)\begin{array}{c}\\ \end{array}\\ +\left(1-2Q\left|\bar{x}-\bar{y}\right|\right)\langle \frac{\bar{x}-\bar{y}}{\left|\bar{x}-\bar{y}\right|},\stackrel{\to }{n}\left(\bar{x}\right)\rangle \}\\ \ge {\text{e}}^{-\left(\left|\beta \right|+1\right)Md\left(\bar{y}\right)}PQ\left\{\frac{3}{4}\left(\left|\beta \right|+1\right)M\left|\bar{x}-\bar{y}\right|-\left(1-2Q\left|\bar{x}-\bar{y}\right|\right)\frac{\left|\bar{x}-\bar{y}\right|}{2r}\right\}\\ \ge \frac{1}{2}{\text{e}}^{-\left(\left|\beta \right|+1\right)Md\left(\bar{y}\right)}PQ\left|\bar{x}-\bar{y}\right|\left(\frac{3}{2}M\left(\left|\beta \right|+1\right)-\frac{1}{r}\right)>0,\end{array}$

where $r$ is the radius in the uniform exterior sphere condition (Ω2) and we have chosen $M>\frac{2}{3\left(\left|\beta \right|+1\right)r}$ .

Similarly, if $\bar{y}\in \partial \Omega$ , one has

$\langle -D_y\psi \left(\bar{x},\bar{y}\right)\mathrm{,}\stackrel{\to }{n}\left(\bar{y}\right)\rangle \le \frac{1}{2}{\text{e}}^{-\left(\left|\beta \right|+1\right)Md\left(\bar{x}\right)}PQ\left|\bar{x}-\bar{y}\right|\left(-\frac{3}{2}M\left(\left|\beta \right|+1\right)+\frac{1}{r}\right)<0.$

Since $u$ is a viscosity subsolution and $v$ is a viscosity supersolution, we obtain

$\begin{array}{l}\text{tr}\left(\sigma \left(D_x\psi \left(\bar{x},\bar{y}\right)\right)X\right)+\beta \left|D_x\psi \left(\bar{x},\bar{y}\right)\right|+\xi \left(\bar{x}\right)\cdot D_x\psi \left(\bar{x},\bar{y}\right)+\eta \left(\bar{x}\right)u\left(\bar{x}\right)\ge g\left(\bar{x}\right),\\ \text{if}\left(D_x\psi \left(\bar{x},\bar{y}\right),X\right)\in {\bar{J}}^{2,+}u\left(\bar{x}\right),\end{array}$

and

$\begin{array}{l}\text{tr}\left(\sigma \left(D_y\psi \left(\bar{x},\bar{y}\right)\right)Y\right)+\beta \left|D_y\psi \left(\bar{x},\bar{y}\right)\right|-\xi \left(\bar{y}\right)\cdot D_y\psi \left(\bar{x},\bar{y}\right)+\eta \left(\bar{y}\right)v\left(\bar{y}\right)\le h\left(\bar{y}\right),\\ \text{if}\left(-D_y\psi \left(\bar{x},\bar{y}\right),Y\right)\in {\bar{J}}^{2,-}v\left(\bar{y}\right).\end{array}$

Then the inequalities (11) and (12) hold for any maximum point $\left(\bar{x},\bar{y}\right)\in {\Delta }_Q\cap {\bar{\Omega }}^2$ , provided $Q$ is large enough.

Step 2. For every $\epsilon >0$ , there exist $X\mathrm{,}Y\in S\left(\mathbb{N}\right)$ such that $\left(D_x\psi \left(\bar{x},\bar{y}\right)\mathrm{,}X\right)\in {\bar{J}}^{\mathrm{2,}+}u\left(\bar{x}\right)$ , $\left(-D_y\psi \left(\bar{x},\bar{y}\right)\mathrm{,}Y\right)\in {\bar{J}}^{\mathrm{2,}-}v\left(\bar{y}\right)$ and

$\left(\begin{array}{cc}X& 0\\ 0& -Y\end{array}\right)\le D^2\psi \left(\bar{x},\bar{y}\right)+\epsilon {\left(D^2\psi \left(\bar{x},\bar{y}\right)\right)}^2\mathrm{.}$ (13)

Next we estimate the right-hand side of the inequality (13):

$\begin{array}{c}{D}^2\psi \left(\bar{x},\bar{y}\right)=\Psi \left(\bar{x}-\bar{y}\right)D^2\left({\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\right)+D\left({\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\right)\\ \otimes D\left(\Psi \left(\bar{x}-\bar{y}\right)\right)+D\left(\Psi \left(\bar{x}-\bar{y}\right)\right)\otimes D\left({\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\right)\\ +{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}{D}^2\left(\Psi \left(\bar{x}-\bar{y}\right)\right).\end{array}$

We denote

$A_1:=\Psi \left(\bar{x}-\bar{y}\right)D^2\left({\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\right),$

$\begin{array}{l}{A}_2:=D\left({\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\right)\otimes D\left(\Psi \left(\bar{x}-\bar{y}\right)\right)\\ +D\left(\Psi \left(\bar{x}-\bar{y}\right)\right)\otimes D\left({\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\right),\end{array}$

$A_3:={\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}{D}^2\left(\Psi \left(\bar{x}-\bar{y}\right)\right).$

One has

$A_1\le CQ\left|\bar{x}-\bar{y}\right|\left(\begin{array}{cc}{I}_{\mathbb{N}}& 0\\ 0& I_{\mathbb{N}}\end{array}\right)\mathrm{,}$ (14)

and

$A_2\le CQ\left(\begin{array}{cc}{I}_{\mathbb{N}}& 0\\ 0& I_{\mathbb{N}}\end{array}\right)+CQ\left(\begin{array}{cc}{I}_{\mathbb{N}}& -I_{\mathbb{N}}\\ -I_{\mathbb{N}}& I_{\mathbb{N}}\end{array}\right)\mathrm{.}$ (15)

Indeed, for $\rho \mathrm{,}\tau \in {ℝ}^{\mathbb{N}}$ , we have

$\begin{array}{l}\langle A_2\left(\rho ,\tau \right),\left(\rho ,\tau \right)\rangle \\ =2\left(\left|\beta \right|+1\right)M{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\{\langle Dd\left(\bar{x}\right)\otimes D\Psi \left(\bar{x}-\bar{y}\right)\left(\tau -\rho \right),\rho \rangle \\ +\langle Dd\left(\bar{y}\right)\otimes D\Psi \left(\bar{x}-\bar{y}\right)\left(\tau -\rho \right),\tau \rangle \}\\ \le CQ\left(\left|\rho \right||+\left|\tau \right|\right)\left|\tau -\rho \right|\\ \le CQ\left({\left|\rho \right|}^2+{\left|\tau \right|}^2\right)+CQ{\left|\tau -\rho \right|}^2,\end{array}$

where $C$ denotes various positive constants independent of $Q$ .

Now we are ready to estimate A₃. For $D^2\left(\Psi \left(\bar{x}-\bar{y}\right)\right)$ , one has

$D^2\left(\Psi \left(\bar{x}-\bar{y}\right)\right)=\left(\begin{array}{cc}{D}^2\Psi \left(\bar{x}-\bar{y}\right)& -D^2\Psi \left(\bar{x}-\bar{y}\right)\\ -D^2\Psi \left(\bar{x}-\bar{y}\right)& D^2\Psi \left(\bar{x}-\bar{y}\right)\end{array}\right)\mathrm{,}$

and the Hessian matrix of $\Psi \left(x\right)$ is

$D^2\Psi \left(x\right)=\frac{PQ}{\left|x\right|}\left(I_{\mathbb{N}}-\frac{x\otimes x}{{\left|x\right|}^2}\right)-2P{Q}^2{I}_{\mathbb{N}}\mathrm{.}$ (16)

Denoting

$\epsilon =\frac{\left|\bar{x}-\bar{y}\right|}{2PQ{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}}\mathrm{,}$ (17)

we obtain

$\begin{array}{l}\epsilon A_1^2\le CQ{\left|\bar{x}-\bar{y}\right|}^3{I}_{2\mathbb{N}}\mathrm{,}\epsilon A_2^2\le CQ\left|\bar{x}-\bar{y}\right|I_{2\mathbb{N}}\mathrm{,}\\ \epsilon \left(A_1{A}_2+A_2{A}_1\right)\le CQ{\left|\bar{x}-\bar{y}\right|}^2{I}_{2\mathbb{N}}\mathrm{,}\\ \epsilon \left(A_1{A}_3+A_3{A}_1\right)\le CQ\left|\bar{x}-\bar{y}\right|I_{2\mathbb{N}}\mathrm{,}\epsilon \left(A_2{A}_3+A_3{A}_2\right)\le CQ{I}_{2\mathbb{N}}\mathrm{,}\end{array}$ (18)

where

$I_{2\mathbb{N}}=\left(\begin{array}{cc}{I}_{\mathbb{N}}& 0\\ 0& I_{\mathbb{N}}\end{array}\right)\mathrm{.}$

By the inequalities (14), (15), (18) and

${\left(D^2\left(\Psi \left(\bar{x}-\bar{y}\right)\right)\right)}^2=\left(\begin{array}{cc}2{\left(D^2\Psi \left(\bar{x}-\bar{y}\right)\right)}^2& -2{\left(D^2\Psi \left(\bar{x}-\bar{y}\right)\right)}^2\\ -2{\left(D^2\Psi \left(\bar{x}-\bar{y}\right)\right)}^2& 2{\left(D^2\Psi \left(\bar{x}-\bar{y}\right)\right)}^2\end{array}\right),$

one obtains

$\left(\begin{array}{cc}X& 0\\ 0& -Y\end{array}\right)\le O\left(Q\right)\left(\begin{array}{cc}{I}_{\mathbb{N}}& 0\\ 0& I_{\mathbb{N}}\end{array}\right)+\left(\begin{array}{cc}B& -B\\ -B& B\end{array}\right)\mathrm{,}$ (19)

where

$B=CQ{I}_{\mathbb{N}}+{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\left[D^2\Psi \left(\bar{x}-\bar{y}\right)+\frac{\left|\bar{x}-\bar{y}\right|}{PQ}{\left(D^2\Psi \left(\bar{x}-\bar{y}\right)\right)}^2\right]\mathrm{.}$

Thus, we can rewrite the inequality (19) as

$\left(\begin{array}{cc}\tilde{X}& 0\\ 0& -\tilde{Y}\end{array}\right)\le \left(\begin{array}{cc}B& -B\\ -B& B\end{array}\right)\mathrm{,}$ (20)

where $\tilde{X}=X-O\left(Q\right)I_{\mathbb{N}}$ , $\tilde{Y}=Y+O\left(Q\right)I_{\mathbb{N}}$ .

Multiplying on the left of the inequality (20) by the non-negative symmetric matrix

$\left(\begin{array}{cc}\sigma \left(D_x\psi \left(\bar{x},\bar{y}\right)\right)& 0\\ 0& \sigma \left(D_y\psi \left(\bar{x},\bar{y}\right)\right)\end{array}\right),$

one has

$\begin{array}{l}\text{tr}\left(\sigma \left(D_x\psi \left(\bar{x},\bar{y}\right)\right)\tilde{X}\right)-\text{tr}\left(\sigma \left(D_y\psi \left(\bar{x},\bar{y}\right)\right)\tilde{Y}\right)\\ \le \text{tr}\left(\sigma \left(D_x\psi \left(\bar{x},\bar{y}\right)\right)B\right)+\text{tr}\left(\sigma \left(D_y\psi \left(\bar{x},\bar{y}\right)\right)B\right).\end{array}$ (21)

We aim to get the estimate on the right side of the inequality (21). Next we define

$0\le H:=\frac{\left(\bar{x}-\bar{y}\right)\otimes \left(\bar{x}-\bar{y}\right)}{{\left|\bar{x}-\bar{y}\right|}^2}\le I_{\mathbb{N}},$

and compute $\text{tr}\left(HB\right)$ . Since $\frac{x\otimes x}{{\left|x\right|}^2}$ is idempotent, one has

${\left(D^2\Psi \left(x\right)\right)}^2=\frac{P^2{Q}^2}{{\left|x\right|}^2}\left(1-4Q\left|x\right|\right)\left(I_{\mathbb{N}}-\frac{x\otimes x}{{\left|x\right|}^2}\right)+4P^2{Q}^4{I}_{\mathbb{N}}\mathrm{.}$

For large $Q$ , since $\text{tr}H=1$ and $4Q\left|\bar{x}-\bar{y}\right|\le 1$ , we get

$\begin{array}{c}\text{tr}\left(HB\right)=CQ+{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}\left(-2P{Q}^2+4P{Q}^3\left|\bar{x}-\bar{y}\right|\right)\\ \le CQ-{\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}P{Q}^2\le -C{Q}^2.\end{array}$

Thus, we can write $D_x\psi \left(\bar{x},\bar{y}\right)$ as

$D_x\psi \left(\bar{x},\bar{y}\right)={\text{e}}^{-\left(\left|\beta \right|+1\right)M\left(d\left(\bar{x}\right)+d\left(\bar{y}\right)\right)}PQ\left(v_1+v_2\right)\mathrm{,}$

where

$v_1=-\left(\left|\beta \right|+1\right)M\left|\bar{x}-\bar{y}\right|\left(1-Q\left|\bar{x}-\bar{y}\right|\right)Dd\left(\bar{x}\right)\mathrm{,}$

and

$v_2=\left(1-2Q\left|\bar{x}-\bar{y}\right|\right)\frac{\bar{x}-\bar{y}}{\left|\bar{x}-\bar{y}\right|}\mathrm{.}$

Therefore,

$\sigma \left(D_x\psi \left(\bar{x},\bar{y}\right)\right)=\frac{v_1\otimes v_1}{{\left|v_1+v_2\right|}^2}+\frac{v_1\otimes v_2+v_2\otimes v_1}{{\left|v_1+v_2\right|}^2}+\frac{v_2\otimes v_2}{{\left|v_1+v_2\right|}^2}\mathrm{.}$

Since $Q\left|\bar{x}-\bar{y}\right|\le \frac{1}{4}$ , for large $Q$ , one has

$\frac{1}{4}=\frac{1}{2}-\frac{1}{4}\le \left|v_2\right|-\left|v_1\right|\le \left|v_1+v_2\right|\le \left|v_1\right|+\left|v_2\right|\le 2,$

and

$\Vert B\Vert \le \frac{CQ}{\left|\bar{x}-\bar{y}\right|}\mathrm{.}$

Then

$\left|\text{tr}\left(\frac{v_1\otimes v_1}{{\left|v_1+v_2\right|}^2}B\right)\right|\le C{\left|\bar{x}-\bar{y}\right|}^2\Vert B\Vert \le CQ\left|\bar{x}-\bar{y}\right|\mathrm{,}$

$\left|\text{tr}\left(\frac{v_1\otimes v_2+v_2\otimes v_1}{{\left|v_1+v_2\right|}^2}B\right)\right|\le C\left|\bar{x}-\bar{y}\right|\Vert B\Vert \le CQ\mathrm{,}$

and

$\text{tr}\left(\frac{v_2\otimes v_2}{{\left|v_1+v_2\right|}^2}B\right)=\frac{\text{tr}\left(HB\right)}{{\left|v_1+v_2\right|}^2}\le -C{Q}^2\mathrm{.}$

We conclude that

$\text{tr}\left(\sigma \left(D_x\psi \left(\bar{x},\bar{y}\right)\right)B\right)\le O\left(Q\right)-C{Q}^2\mathrm{.}$

Similarly, we can get the following estimate

$\text{tr}\left(\sigma \left(D_y\psi \left(\bar{x},\bar{y}\right)\right)B\right)\le O\left(Q\right)-C{Q}^2\mathrm{.}$

Therefore, by the inequality (21) one has

$\text{tr}\left(\sigma \left(D_x\psi \left(\bar{x},\bar{y}\right)\right)\tilde{X}\right)-\text{tr}\left(\sigma \left(D_y\psi \left(\bar{x},\bar{y}\right)\right)\tilde{Y}\right)\le O\left(Q\right)-C{Q}^2\mathrm{.}$

Step 3. By the definition of $\tilde{X}$ and $\tilde{Y}$ and the fact that $u$ , $v$ are respectively viscosity subsolution and viscosity supersolution, one has

$\begin{array}{c}g\left(\bar{x}\right)-\eta \left(\bar{x}\right)u\left(\bar{x}\right)\le \text{tr}\left(\sigma \left(D_x\psi \right)X\right)+\beta \left|D_x\psi \right|+\xi \left(\bar{x}\right)\cdot D_x\psi \\ \le \text{tr}\left(\sigma \left(D_x\psi \right)\tilde{X}\right)+O\left(Q\right)+\beta \left|D_x\psi \right|+\xi \left(\bar{x}\right)\cdot D_x\psi \\ \le \text{tr}\left(\sigma \left(D_y\psi \right)Y\right)+O\left(Q\right)-C{Q}^2+\beta \left|D_x\psi \right|+\xi \left(\bar{x}\right)\cdot D_x\psi \\ \le -\beta \left|D_y\psi \right|+\xi \left(\bar{y}\right)\cdot D_y\psi -\eta \left(\bar{y}\right)v\left(\bar{y}\right)+h\left(\bar{y}\right)+O\left(Q\right)\\ -C{Q}^2+\beta \left|D_x\psi \right|+\xi \left(\bar{x}\right)\cdot D_x\psi .\end{array}$

According to the inequalities (11) and (12), one gets

$g\left(\bar{x}\right)-h\left(\bar{y}\right)-\eta \left(\bar{x}\right)u\left(\bar{x}\right)+\eta \left(\bar{y}\right)v\left(\bar{y}\right)\le O\left(Q\right)-C{Q}^2\mathrm{.}$ (22)

If $u$ and $v$ are both bounded, the left-hand side of the inequality (22) is bounded from below by $-{\Vert g\Vert }_{\infty }-{\Vert h\Vert }_{\infty }-{\Vert \eta \Vert }_{\infty }\left({\Vert u\Vert }_{\infty }+{\Vert v\Vert }_{\infty }\right)$ . Otherwise, if $v$ is non-negative and bounded, then $u\left(\bar{x}\right)\ge 0$ and that quantity is greater than $-{\Vert g\Vert }_{\infty }-{\Vert h\Vert }_{\infty }-{\Vert \eta \Vert }_{\infty }\left(\sup u+{\Vert v\Vert }_{\infty }\right)$ . On the other hand, the right-hand side of the inequality (22) goes to $-\infty$ as $Q\to +\infty$ . Hence, taking $Q$ large enough, we can obtain a contradiction and this concludes the proof.

Theorem 1 is an immediate consequence of Theorem 2.

Proof of Theorem 1. Since $u\left(x\right)\in C\left(\bar{\Omega }\right)$ is a viscosity solution of the problem (4), $u\left(x\right)$ is both a viscosity subsolution and a viscosity supersolution of the problem (4). Thus, we have $m=0$ . Since $u$ is bounded, by Theorem 2, we can get immediately the Lipschizt estimate (5).

4. Conclusion

In this paper, we establish the Lipschitz regularity of the problem (1) arsing from the generalized random tug-of-war game. The Lipschitz regularity is an indispensable part and an important issue in the study of PDEs. The Lipschitz regularity also plays a vital role in applications, such as image processing, financial problems and physical engineering.

Acknowledgements

We thank the anonymous referees for the careful reading of the manuscript and useful suggestions and comments.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Aronsson, G. (1965) Minimization Problems for the Functional supxF (x, f(x), f’(x)). Arkiv for Matematik, 6, 33-53. https://doi.org/10.1007/BF02591326
[2]	Aronsson, G. (1966) Minimization Problems for the Functional supxF (x, f(x), f’(x)) II. Arkiv for Matematik, 6, 409-431. https://doi.org/10.1007/BF02590964
[3]	Aronsson, G. (1967) Extension of Functions Satisfying Lipschitz Conditions. Arkiv for Matematik, 6, 551-561. https://doi.org/10.1007/BF02591928
[4]	Aronsson, G. (1969) Minimization Problems for the Functional supxF (x, f(x), f’(x)) III. Arkiv for Matematik, 7, 509-512. https://doi.org/10.1007/BF02590888
[5]	Evans, L.C. and Gangbo, W. (1999) Differential Equations Methods for the Monge-Kantorovich Mass Transfer Problem. Memoirs of the American Mathematical Society, 137, 1-75. https://doi.org/10.1090/memo/0653
[6]	Elion, C. and Vese, L.A. (2007) An Image Decomposition Model Using the Total Variation and the Infinity Laplacian. Proceedings Computational Imaging V, Vol. 6498, 261-270. https://doi.org/10.1117/12.716079
[7]	Elmoataz, A., Toutain, M. and Tenbrinck, D. (2015) On the p-Laplacian and ∞-Laplacian on Graphs with Applications in Image and Data Processing. SIAM Journal on Imaging Sciences, 8, 2412-2451. https://doi.org/10.1137/15M1022793
[8]	Atar, R. and Budhiraja, A. (2010) A Stochastic Differential Game for the Inhomogeneous ∞-Laplace Equation. The Annals of Probability, 38, 498-531. https://doi.org/10.1214/09-AOP494
[9]	López-Soriano, R., Navarro-Climent, J.C. and Rossi, J.D. (2013) The Infinity Laplacian with a Transport Term. Journal of Mathematical Analysis and Applications, 398, 752-765. https://doi.org/10.1016/j.jmaa.2012.09.030
[10]	Charro, F., Garca-Azorero, J. and Rossi, J.D. (2009) A Mixed Problem for the Infinity Laplacian via Tug-of-War Games. Calculus of Variations and Partial Differential Equations, 34, 307-320. https://doi.org/10.1007/s00526-008-0185-2
[11]	Garca-Azorero, J., Manfredi, J., Peral, I. and Rossi, J.D. (2007) The Neumann Problem for the ∞-Laplacian and the Monge-Kantorovich Mass Transfer Problem. Nonlinear Analysis: Theory, Methods and Applications, 66, 349-366. https://doi.org/10.1016/j.na.2005.11.030
[12]	Gómze, I. and Rossi, J.D. (2013) Tug-of-War Games and the Infinity Laplacian with Spatial Dependence. Communications on Pure and Applied Analysis, 12, 1959-1983. https://doi.org/10.3934/cpaa.2013.12.1959
[13]	Liu, F. and Yang, X.P. (2012) Solutions to an Inhomogeneous Equation Involving Infinity-Laplacian. Nonlinear Analysis: Theory, Methods and Applications, 75, 5693-5701. https://doi.org/10.1016/j.na.2012.05.017
[14]	Liu, F. (2022) The Eigenvalue Problem for a Class of Degenerate Operators Related to the Normalized p-Laplacian. Discrete and Continuous Dynamical Systems-B, 27, 2701-2720. https://doi.org/10.3934/dcdsb.2021155
[15]	Liu, F., Tian, L. and Zhao, P.B. (2020) A Weighted Eigenvalue Problem of the Degenerate Operator Associated with Infinity Laplacian. Nonlinear Analysis: Theory, Methods and Applications, 200, Article ID: 112001. https://doi.org/10.1016/j.na.2020.112001
[16]	Nyström, K. and Parviainen, M. (2014) Tug-of-War, Market Manipulation, and Option Pricing. Mathematical Finance, 27, 279-312. https://doi.org/10.1111/mafi.12090
[17]	Peres, Y., Pete, G. and Somersille, S. (2010) Biased Tug-of-War, the Biased Infinity Laplacian, and Comparison with Exponential Cones. Calculus of Variations and Partial Differential Equations, 38, 541-564. https://doi.org/10.1007/s00526-009-0298-2
[18]	Peres, Y. and Sheffield, S. (2008) Tug-of-War with Noise: A Game-Theoretic View of the p-Laplacian. Duke Mathematical Journal, 145, 91-120. https://doi.org/10.1215/00127094-2008-048
[19]	Peres, Y., Schramm, O., Sheffield, S. and Wilson, D. (2009) Tug-of-War and the Infinity Laplacian. Journal of the American Mathematical Society, 22, 167-210. https://doi.org/10.1090/S0894-0347-08-00606-1
[20]	Armstrong, S.N., Smart, C.K. and Somersille, S.J. (2011) An Infinity Laplace Equation with Gradient Term and Mixed Boundary Conditions. Proceedings of the American Mathematical Society, 139, 1763-1776. https://doi.org/10.1090/S0002-9939-2010-10666-4
[21]	Liu, F. and Yang, X.P. (2021) A Weighted Eigenvalue Problem of the Biased Infinity Laplacian. Nonlinearity, 34, 1197-1237. https://doi.org/10.1088/1361-6544/abd85d
[22]	Lu, G.Z. and Wang, P.Y. (2008) A PDE Perspective of the Normalized Infinity Laplacian. Communications in Partial Differential Equations, 33, 1788-1817. https://doi.org/10.1080/03605300802289253
[23]	Lu, G.Z. and Wang, P.Y. (2008) Inhomogeneous Infinity Laplace Equation. Advances in Mathematics, 217, 1838-1868. https://doi.org/10.1016/j.aim.2007.11.020
[24]	Patrizi, S. (2011) The Principal Eigenvalue of the ∞-Laplacian with the Neumann Boundary Condition. ESAIM: Control, Optimisation and Calculus of Variations, 17, 575-601. https://doi.org/10.1051/cocv/2010019
[25]	Aronsson, G. (1984) On Certain Singular Solutions of the Partial Differential Equation u2xuxx + 2uxuyuxy + u2yuyy = 0. Manuscripta Mathematica, 47, 133-151. https://doi.org/10.1007/BF01174590
[26]	Savin, O. (2005) C1 Regularity for Infinity Harmonic Functions in Two Dimensions. Archive for Rational Mechanics and Analysis, 176, 351-361. https://doi.org/10.1007/s00205-005-0355-8
[27]	Evans, L.C. and Savin, O. (2008) C1,α Regularity for Infinity Harmonic Functions in Two Dimensions. Calculus of Variations and Partial Differential Equations, 32, 325-347. https://doi.org/10.1007/s00526-007-0143-4
[28]	Evans, L.C. and Smart, C.K. (2011) Adjoint Methods for the Infinity Laplacian Partial Differential Equation. Archive for Rational Mechanics and Analysis, 201, 87-113. https://doi.org/10.1007/s00205-011-0399-x
[29]	Evans, L.C. and Smart, C.K. (2011) Everywhere Differentiability of Infinity Harmonic Functions. Calculus of Variations and Partial Differential Equations, 42, 289-299. https://doi.org/10.1007/s00526-010-0388-1
[30]	Crandall, M.G., Ishii, H. and Lions, P.L. (1992) User’s Guide to Viscosity Solutions of Second Order Partial Differential Equations. Bulletin of the American Mathematical Society, 27, 1-67. https://doi.org/10.1090/S0273-0979-1992-00266-5
[31]	Birindelli, I. and Demengel, F. (2007) Eigenvalue, Maximum Principle and Regularity for Fully Nonlinear Homogeneous Operators. Communications on Pure and Applied Analysis, 6, 335-366. https://doi.org/10.3934/cpaa.2007.6.335
[32]	Crandall, M.G., Evans, L.C. and Gariepy, R.F. (2001) Optimal Lipschitz Extensions and the Infinity Laplacian. Calculus of Variations and Partial Differential Equations, 13, 123-139. https://doi.org/10.1007/s005260000065