Higher regularity in congested traffic dynamics.

In this paper, we consider minimizers of integral functionals of the type F (u) : = ∫ Ω [ 1 p (| D u | - 1) + p + f · u ] d x for p > 1 in the vectorial case of mappings u : R n ⊃ Ω → R N with N ≥ 1 . Assuming that f belongs to L n + σ for some σ > 0 , we prove that H (D u) is continuous in Ω for any continuous function H : R Nn → R Nn vanishing on { ξ ∈ R Nn : | ξ | ≤ 1 } . This extends previous results of Santambrogio and Vespri (Nonlinear Anal 73:3832–3841, 2010) when n = 2 , and Colombo and Figalli (J Math Pures Appl (9) 101(1):94–117, 2014) for n ≥ 2 , to the vectorial case N ≥ 1 .

Keywords: 49N60; 35D10; 35J70

Introduction and main result

In this paper, we study $C^1$ -regularity of minimizers of integral functionals of the Calculus of Variations with widely degenerate convex integrands of the form

1.1 $$\begin{array}{c}\hfill \mathcal{F}(u):={\int }_{\mathrm{\Omega }}[F(Du)+f·u]\mathrm{d}x,\end{array}$$

Graph

where $\mathrm{\Omega }\subset {\mathbb{R}}^n$ , $n\ge 2$ , is a bounded domain and $u:\mathrm{\Omega }\to {\mathbb{R}}^N$ , $N\ge 1$ , a possibly vector valued function. We concentrate ourself on the study of the prototype integrand

1.2 $$\begin{array}{c}\hfill F(\xi ):=\frac{1}{p}{(|\xi |-1)}_{+}^p,\end{array}$$

Graph

for some $p>1$ . The datum f is required to belong to $L^{n+\sigma }$ for some $\sigma >0$ . The functional $\mathcal{F}$ with the specific integrand F from (1.2) is the prototype for a class of more general functionals where F is a convex function vanishing inside some convex set, and satisfying specific growth and ellipticity assumptions. For sake of clarity, the results in this paper are stated and proved for the functionals $\mathcal{F}(u)$ as in (1.1)–(1.2). However, we expect our techniques to apply to a general class of integrands with a widely degenerate structure as well. The functional $\mathcal{F}(u)$ and its associated Euler–Lagrange system

1.3 $$\begin{array}{c}\hfill \mathrm{div}((|Du|-1{)}_{+}^{p-1}\frac{\mathit{Du}}{|Du|})=f\end{array}$$

Graph

naturally arise in problems of optimal transport with congestion effects. In fact, minimizing (1.1) with $N=1$ and the integrand from (1.2) is equivalent to the dual minimization problem

1.4 $$\begin{array}{c}\hfill min\{{\int }_{\mathrm{\Omega }}\mathcal{H}(\sigma )\mathrm{d}x:\sigma \in L^q(\mathrm{\Omega },{\mathbb{R}}^n),\mathrm{div}\sigma =f,\sigma ·{\nu }_{\partial \mathrm{\Omega }}=0\},\end{array}$$

Graph

where the integrand

$$\begin{array}{c}\hfill \mathcal{H}(\sigma )=H(|\sigma |),\text{with}H(t)=t+\frac{1}{q}{t}^q\text{and}\frac{1}{p}+\frac{1}{q}=1\end{array}$$

Graph

is the convex conjugate of F, or equivalently $F={\mathcal{H}}^{\ast }$ , and $\sigma$ represents the traffic flow. The function $g(t)=H^{\prime }(t)$ models the congestion effect. Note that $\sigma \mapsto g(\sigma )$ is increasing and $g(0)=1>0$ , so that moving in an empty street has nonzero cost. As shown in [[4]] the unique minimizer $\sigma (x)$ of (1.4) is given by $D_{\xi }F(Du(x))$ . We refer to [[2]–[4], [6], [36]] and the references therein for detailed motivations and for the physical meaning of the regularity of minimizers. It would be interesting to investigate if there are applications of the vectorial problem. However, our main motivation to consider the very degenerate system (1.3) was from a mathematical point of view.

In connection with congested traffic dynamic problems the regularity of minimizers, as well as the regularity of weak solutions of the associated autonomous Euler–Lagrange system has been an active field of research in recent years. For instance, in [[2], [4], [9], [24]] Lipschitz regularity of minimizers has been established under suitable assumptions on the datum f.

At this point it is worthwhile to observe that, in general, no more than Lipschitz regularity can be expected for solutions of equations or systems as in (1.3). Indeed when $f=0$ , every 1-Lipschitz continuous function solves (1.3). On the other hand, in the scalar case $N=1$ assuming $f\in L^{n+\sigma }$ for some $\sigma >n$ , it was shown by Santambrogio and Vespri [[36]] for $n=2$ and Figalli and Colombo [[10]] for $n\ge 2$ , that the composition of an arbitrary continuous function vanishing on the set $\{|\mathrm{\nabla }u|\le 1\}$ with $\mathrm{\nabla }u$ is continuous.

Our aim in this paper is to investigate $C^1$ -regularity of minimizers in the vectorial case $N\ge 1$ . In general, regularity in the vectorial case is much more delicate and minimizers may be irregular although the integrand is smooth, cf. [[15], [38]]. In this respect, regularity can be expected only for integrands with special structure. The first result in this direction has been obtained by Uhlenbeck [[40]] for the p-Laplace system when $p\ge 2$ . She proved that weak solutions are of class $C^{1,\alpha }$ . The scalar case had previously been established by Ural'ceva [[41]], while the case $p\in (1,2)$ was obtained by Tolksdorf [[39]]. As already mentioned, we cannot hope for a $C^{1,\alpha }$ -regularity result for the elliptic system (1.3), since Lipschitz continuity is optimal. However, we are able to establish in the vectorial setting that the composition $\mathcal{K}(Du)$ is continuous for any continuous function $\mathcal{K}:{\mathbb{R}}^{\mathit{Nn}}\to \mathbb{R}$ vanishing on $\{\xi \in {\mathbb{R}}^{\mathit{Nn}}:|\xi |\le 1\}$ . This phenomenon is somewhat reminiscent of comparable results for the Stefan problem, in which the continuity of the energy cannot been shown, but the temperature shows a logarithmic type continuity, cf. [[16], [31]].

Statement of the main result

Before formulating the main regularity result, we need to introduce a few notations. The natural energy space to deal with (local) minimizers of the integral functional $\mathcal{F}$ is the Sobolev space $W^{1,p}(\mathrm{\Omega },{\mathbb{R}}^N)$ . Then, (local) minimizers in $W_{\mathrm{loc}}^{1,p}(\mathrm{\Omega },{\mathbb{R}}^N)$ of the functional $\mathcal{F}$ are weak solutions of the Euler–Lagrange system

1.5 $$\begin{array}{c}\hfill \mathrm{div}\mathbf{A}(Du)=f\end{array}$$

Graph

and vice versa, where

$$\begin{array}{c}\hfill \mathbf{A}(\xi ):=h(|\xi |)\xi ,\text{with}h(t):=\frac{{(t-1)}_{+}^{p-1}}{t}\text{for}t\in {\mathbb{R}}_{+},\end{array}$$

Graph

for some $p>1$ . A function $u\in W_{\mathrm{loc}}^{1,p}(\mathrm{\Omega },{\mathbb{R}}^N)$ is a weak solution of the Euler–Lagrange system (1.5) if and only if

$$\begin{array}{c}\hfill {\int }_{\mathrm{\Omega }}\mathbf{A}(Du)·D\phi \mathrm{d}x=-{\int }_{\mathrm{\Omega }}f·\phi \mathrm{d}x\end{array}$$

Graph

holds true for any testing function $\phi \in C_0^{\infty }(\mathrm{\Omega },{\mathbb{R}}^N)$ . Our main result proves the continuity of the composition $\mathcal{K}(Du)$ explained above.

Theorem 1.1

Let $p>1$ , $f\in L^{n+\sigma }(\mathrm{\Omega },{\mathbb{R}}^N)$ for some $\sigma >0$ and $u\in W_{\mathrm{loc}}^{1,p}(\mathrm{\Omega },{\mathbb{R}}^N)$ be a weak solution of (1.5) in $\mathrm{\Omega }$ . Then,

$$\begin{array}{c}\hfill \mathcal{K}(Du)\in C^0(\mathrm{\Omega })\end{array}$$

Graph

for any continuous function $\mathcal{K}:{\mathbb{R}}^{\mathit{Nn}}\to \mathbb{R}$ vanishing on $\{\xi \in {\mathbb{R}}^{\mathit{Nn}}:|\xi |\le 1\}$ .

By carefully tracing the dependence of constants on the parameter $\delta$ in the proof of Theorem 3.6, one could determine an explicit modulus of continuity of $\mathcal{G}(Du)$ , where $\mathcal{G}$ is defined in (2.2). However, it is not clear if $\mathcal{G}(Du)$ is Hölder continuous in general. For a different very degenerate elliptic equation a counterexample to Hölder continuity is provided in [[11]].

Theorem 1.1 can be regarded as the vectorial analog of the regularity results of Santambrogio and Vespri [[36], Theorem 11] and of Figalli and Colombo [[11], Theorem 1.1] as far as the model type integral functional is considered. The vectorial case cannot be treated with the methods from [[11], [36]], since these are tailored to the scalar case. Nevertheless, some steps in our proof are similar, for example, the approximation procedure by some sequence of uniformly elliptic problems. The main difference in our proof is that we are establishing a variant of DiBenedetto's and Manfredi's proofs of $C^{1,\alpha }$ -regularity of minimizers to p-energy type functionals [[17], [32]]. It is also inspired by the arguments from DiBenedetto and Friedman's pioneering results on $C^{1,\alpha }$ -regularity for parabolic p-Laplacian systems [[18], [20]]. Roughly speaking, our strategy is the adaptation of De Giorgi's approach to the level of gradients in combination with Campanato type comparison arguments. The past has shown that De Giorgi's approach is extremely flexible. Therefore, we expect that our approach can be transferred to larger classes of widely degenerate functionals in the vectorial case. However, and in order to keep the individual steps as simple as possible, we limit ourselves to treating the model case.

Strategy of the proof

Concerning the overall strategy of proof a few words are in order. First, we observe that weak solutions of (1.5) are Lipschitz continuous. This has been proved in [[2], [4], [9]]. Moreover, functionals as in (1.1) fit into the broader context of asymptotically convex functionals, i.e. functionals having a p-Laplacian type structure only at infinity. This class of functionals has been widely studied, since the local Lipschitz regularity result by Chipot and Evans [[8]]. In particular we mention generalizations allowing super- and sub-quadratic growth [[28], [30], [35]], lower order terms [[34]]. Extensions to various other settings can be found in the non-complete list [[12]–[14], [23]–[27], [37]].

The proof of Theorem 1.1 is divided into several steps and starts by an approximation procedure. Indeed, by replacing h(t) by $h_{\epsilon }(t):=h(t)+\epsilon$ for $\epsilon >0$ and considering instead of (1.5) the Dirichlet-problem on a ball compactly contained in $\mathrm{\Omega }$ associated to the regularized coefficients $h_{\epsilon }$ and with Dirchlet boundary datum u we obtain a sequence of approximating more regular mappings $u_{\epsilon }$ . In particular, $u_{\epsilon }$ has second weak derivatives in $L_{\text{loc}}^2$ . In Sect. 3.1 we summarize the most important properties, i.e. uniform energy bounds, uniform quantitative interior $L^{\infty }$ -gradient bounds, uniform quantitative higher differentiability $W^{2,2}$ -estimates, and finally strong $L^p$ -convergence of ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })\to {\mathcal{G}}_{\delta }(Du)$ in the limit $\delta \downarrow 0$ . The nonlinear mapping ${\mathcal{G}}_{\delta }:{\mathbb{R}}^{\mathit{Nn}}\to {\mathbb{R}}^{\mathit{Nn}}$ with $\delta \in (0,1]$ is defined by

$$\begin{array}{c}\hfill {\mathcal{G}}_{\delta }(\xi ):=\frac{{(|\xi |-1-\delta )}_{+}}{|\xi |}\xi ,\text{for}\xi \in {\mathbb{R}}^{\mathit{Nn}}.\end{array}$$

Graph

Observe that ${\mathcal{G}}_{\delta }$ vanishes on the larger set $\{|\xi |\le 1+\delta \}$ . The reason for considering ${\mathcal{G}}_{\delta }$ is that on the complement of $\{|\xi |>1+\delta \}$ the system (1.5) behaves non-degenerate in the sense that the vector field $\mathbf{A}$ admits a uniform ellipticity bound from below, of course, with constants depending on $\delta$ . This point of view has already been exploited in [[11], [36]]. As a first main result we prove that ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ is Hölder continuous uniformly with respect to $\epsilon$ . However, the constants in the quantitative estimate, i.e. the Hölder exponent and the Hölder norm, may blow up when $\delta \downarrow 0$ . We distinguish between two different regimes: the degenerate and non-degenerate regime. The degenerate regime is characterized by the fact that the measure of those points in a ball in which $|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })|$ is far from its supremum is large, while the non-degenerate regime is characterized by the opposite. In the non-degenerate regime we compare $u_{\epsilon }$ with a solution of a linearized system. This allows us to derive a quantitative $L^2$ -excess-improvement for ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })$ on some smaller ball (see Proposition 3.4). This step utilizes a suitable comparison estimate and the higher integrability of $u_{\epsilon }$ . On the smaller ball we are again in the non-degenerate regime, so that the argument can be iterated yielding a Campanato-type estimate for the $L^2$ -excess of ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })$ . In the degenerate regime we establish that $U_{\epsilon }:=(|D{u}_{\epsilon }{|-1-\delta )}_{+}^2$ is a subsolution to a linear uniformly elliptic equation with measurable coefficients; of course the ellipticity constants depend on $\delta$ and blow up as $\delta \downarrow 0$ . At this stage a De Giorgi type argument allows a reduction of the modulus of ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ on some smaller ball (see Proposition 3.5). However, on this smaller scale it is not clear whether or not we are in the degenerate or non-degenerate regime. Therefore one needs to distinguish between these two regimes again. In the non-degenerate regime we can conclude as above, while in the degenerate regime the reduction of the modulus of ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ applies again. This argument can be iterated as long as we stay in the degenerate regime. However, if at a certain scale the switching from degenerate to non-degenerate occurs, the above Campanato type decay applies. If no switching occurs, we have at any scale of the iteration process a reduction of the modulus of ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ . This, however, shows that the supremum of $|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })|$ —and hence also the one of $|{\mathcal{G}}_{2\delta }(D{u}_{\epsilon })|$ —on shrinking concentric balls converges to 0. Altogether this leads to a quantitative Hölder estimate for ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })$ which remains stable under the already established convergence $u_{\epsilon }\to u$ as $\epsilon \downarrow 0$ . The final step consits in passing to the limit $\delta \downarrow 0$ and conclude that $\mathcal{G}(Du):=\frac{{(|Du|-1)}_{+}}{|Du|}Du$ is continuous. This can be achieved by an application of Ascoli–Arzela's theorem. It is here where we loose control on the quantitative Hölder exponent. At this point the continuity of $\mathcal{K}(Du)$ for any continuous function $\mathcal{K}$ vanishing on $\{|\xi |\le 1\}$ is an immediate consequence.

Notation and preliminary results

Notation

For the open ball of radius $\varrho >0$ and center $x_o\in {\mathbb{R}}^n$ we write $B_{\varrho }(x_o)\subset {\mathbb{R}}^n$ . The mean value of a function $v\in L^1(B_{\varrho }(x_o),{\mathbb{R}}^k)$ is defined by

Graph

If the center is clear from the context we omit the reference to the center and write $B_{\varrho }$ respectively ${(v)}_{\varrho }$ for short. For the standard scalar product on Euclidean spaces ${\mathbb{R}}^k$ as well as the space ${\mathbb{R}}^{\mathit{kn}}$ of $k\times n$ matrices, we use the notation $\xi ·\eta$ . Finally, we use the notion $\mathrm{\nabla }u$ for the gradient of a scalar function u, while we use Du for a vector field u.

Throughout this paper we abbreviate

2.1 $$\begin{array}{c}\hfill g(t):=\frac{{(t-1)}_{+}^p}{t}\text{for}t\in {\mathbb{R}}_{+}\end{array}$$

Graph

and

2.2 $$\begin{array}{c}\hfill \mathcal{G}(\xi ):=\frac{{(|\xi |-1)}_{+}}{|\xi |}\xi \text{for}\xi \in {\mathbb{R}}^k,k\in \mathbb{N}.\end{array}$$

Graph

Observe that

$$\begin{array}{c}\hfill g(|\xi |)\xi ={|\mathcal{G}(\xi )|}^{p-1}\mathcal{G}(\xi )\text{for}\text{any}\xi \in {\mathbb{R}}^k.\end{array}$$

Graph

Moreover, for $\delta \in (0,1]$ we define

2.3 $$\begin{array}{c}\hfill {\mathcal{G}}_{\delta }(\xi ):=\frac{{(|\xi |-1-\delta )}_{+}}{|\xi |}\xi \text{for}\xi \in {\mathbb{R}}^k,k\in \mathbb{N}\end{array}$$

Graph

and note that ${\mathcal{G}}_0\equiv \mathcal{G}$ .

Generic constants are denoted by c. They may vary from line to line. Relevant dependencies on parameters and special constants will be suitably emphasized using parentheses or subscripts.

Algebraic inequalities

In this section, we summarize the relevant algebraic inequalities that will be needed later on. The first lemma follows from an elementary computation.

Lemma 2.1

For $\eta ,\zeta \in {\mathbb{R}}_{\ne 0}^k$ , $k\in \mathbb{N}$ we have

$$\begin{array}{c}\hfill |\frac{\eta }{|\eta |}-\frac{\zeta }{|\zeta |}|\le \frac{2}{|\eta |}|\eta -\zeta |.\end{array}$$

Graph

The next lemma can be deduced as in [[29], Lemma 8.3].

Lemma 2.2

For any $\alpha >0$ , there exists a constant $c=c(\alpha )$ such that, for all $\eta ,\zeta \in {\mathbb{R}}_{\ne 0}^k$ , $k\in \mathbb{N}$ , we have

$$\begin{array}{c}\hfill \frac{1}{c}{||\eta |}^{\alpha -1}\eta -{|\zeta |}^{\alpha -1}\zeta |\le (|\eta |+|\zeta |{)}^{\alpha -1}|\eta -\zeta |\le {c||\eta |}^{\alpha -1}\eta -{|\zeta |}^{\alpha -1}\zeta |.\end{array}$$

Graph

Lemma 2.3

Let $\delta \ge 0$ and $\eta ,\zeta \in {\mathbb{R}}^k$ , $k\in \mathbb{N}$ . Then, for ${\mathcal{G}}_{\delta }$ as defined in (2.3) we have

$$\begin{array}{c}\hfill |{\mathcal{G}}_{\delta }(\eta )-{\mathcal{G}}_{\delta }(\zeta )|\le 3|\eta -\zeta |.\end{array}$$

Graph

Moreover, if $\delta >0$ and $|\eta |\ge 1+\delta$ there holds

$$\begin{array}{c}\hfill |\eta -\zeta |\le (1+\frac{2}{\delta })|\mathcal{G}(\eta )-\mathcal{G}(\zeta )|.\end{array}$$

Graph

Proof

We distinguish between different cases. If $|\eta |,|\zeta |\le 1+\delta$ the inequality holds trivially since ${\mathcal{G}}_{\delta }(\eta )=0={\mathcal{G}}_{\delta }(\zeta )$ . If $|\eta |,|\zeta |>1+\delta$ , we apply Lemma 2.1 and obtain

$$\begin{array}{cc}\hfill |{\mathcal{G}}_{\delta }(\eta )-{\mathcal{G}}_{\delta }(\zeta )|& =|\frac{|\eta |-1-\delta }{|\eta |}\eta -\frac{|\zeta |-1-\delta }{|\zeta |}\zeta |\hfill \\ \hfill & \le |\eta -\zeta |+(1+\delta )|\frac{\eta }{|\eta |}-\frac{\zeta }{|\zeta |}|\hfill \\ \hfill & \le 3|\eta -\zeta |.\hfill \end{array}$$

Graph

If $|\eta |>1+\delta$ and $|\zeta |\le 1+\delta$ , we have

$$\begin{array}{c}\hfill |{\mathcal{G}}_{\delta }(\eta )-{\mathcal{G}}_{\delta }(\zeta )|=|{\mathcal{G}}_{\delta }(\eta )|=|\eta |-1-\delta \le |\eta |-|\zeta |\le |\eta -\zeta |.\end{array}$$

Graph

The case when $|\eta |\le 1+\delta$ and $|\zeta |>1+\delta$ is similar, we just have to interchange the role of $\eta$ and $\zeta$ . Joining the three cases gives the first assertion of the Lemma.

Now, we come to the proof of the second assertion. First, we consider the case $|\zeta |\le 1$ in which $\mathcal{G}(\zeta )=0$ . In this case we have

$$\begin{array}{cc}\hfill \frac{|\eta -\zeta |}{|\mathcal{G}(\eta )-\mathcal{G}(\zeta )|}& =\frac{|\eta -\zeta |}{|\mathcal{G}(\eta )|}=\frac{|\eta -\zeta |}{|\eta |-1}\le \frac{|\eta |+1}{|\eta |-1}=1+\frac{2}{|\eta |-1}\le 1+\frac{2}{\delta }.\hfill \end{array}$$

Graph

Next, we consider the case $|\zeta |>1$ . Recall that by assumption $|\eta |\ge 1+\delta$ . We start with the observation that $\mathcal{G}:{\mathbb{R}}^k\backslash \{|\eta |\le 1\}\to {\mathbb{R}}^k\backslash \{0\}$ is a one to one mapping whose inverse mapping is given by ${\mathcal{G}}^{-1}(\tilde{\xi })=\frac{|\tilde{\xi }|+1}{|\tilde{\xi }|}\tilde{\xi }$ . We let $\tilde{\eta }:=\mathcal{G}(\eta )$ and $\tilde{\zeta }:=\mathcal{G}(\zeta )$ and estimate with Lemma 2.1

$$\begin{array}{c}\hfill \frac{|{\mathcal{G}}^{-1}(\tilde{\eta })-{\mathcal{G}}^{-1}(\tilde{\zeta })|}{|\tilde{\eta }-\tilde{\zeta }|}=\frac{\left|\frac{|\tilde{\eta }|+1}{|\tilde{\eta }|},\tilde{\eta },-,\frac{|\tilde{\zeta }|+1}{|\tilde{\zeta }|},\tilde{\zeta }\right|}{|\tilde{\eta }-\tilde{\zeta }|}=\frac{\left|\tilde{\eta },+,\frac{\tilde{\eta }}{|\tilde{\eta }|},-,\tilde{\zeta },-,\frac{\tilde{\zeta }}{|\tilde{\zeta }|}\right|}{|\tilde{\eta }-\tilde{\zeta }|}\le 1+\frac{2}{|\tilde{\eta }|}\le 1+\frac{2}{\delta }.\end{array}$$

Graph

In the second to last estimate we used $|\tilde{\eta }|>\delta$ . Using in the previous inequality the definition of $\tilde{\eta }$ and $\tilde{\zeta }$ the claim immediately follows. $\square$

Lemma 2.4

There exists a constant $c=c(p)$ such that for any $a>1$ and $b\ge 0$ we have

$$\begin{array}{c}\hfill |h(b)-h(a)|b\le c(p)\frac{{[a-1+{(b-1)}_{+}]}^{p-1}}{a-1}|b-a|.\end{array}$$

Graph

Proof

We apply Lemma 2.2 with $\alpha =p-1>0$ to obtain

$$\begin{array}{cc}\hfill |h(a)-h(b)|b& =|h(a)(b-a)+{(a-1)}^{p-1}-{(b-1)}_{+}^{p-1}|\hfill \\ \hfill & \le {h(a)|b-a|+|(a-1)}^{p-1}-{(b-1)}_{+}^{p-1}|\hfill \\ \hfill & \le {h(a)|b-a|+c[(a-1)+(b-1)}_{+}{]}^{p-2}|b-a|\hfill \\ \hfill & \le c(p)\frac{{[a-1+{(b-1)}_{+}]}^{p-1}}{a-1}|b-a|.\hfill \end{array}$$

Graph

This proves the claim. $\square$

Lemma 2.5

For $a>1$ we have

$$\begin{array}{c}\hfill |h^{\prime }(a)|\le \frac{p{(a-1)}^{p-2}}{a}.\end{array}$$

Graph

Moreover, for $a,b>1$ there holds

$$\begin{array}{c}\hfill |h^{\prime }(b)b-h^{\prime }(a)a|\le c(p)[{(a-1)}^{p-3}+{(b-1)}^{p-3}]|b-a|.\end{array}$$

Graph

Proof

By direct computation we have for $a>1$ that

2.4 $$\begin{array}{cc}\hfill h^{\prime }(a)& =\frac{(p-1){(a-1)}^{p-2}a-{(a-1)}^{p-1}}{a^2}=\frac{{(a-1)}^{p-2}[p-2+\frac{1}{a}]}{a},\hfill \end{array}$$

Graph

from which the first claim immediately follows. We now turn our attention to the second claim. We may assume that $1<a<b$ ; otherwise we interchange the role of a and b. In view of (2.4) we find

$$\begin{array}{c}\hfill |h^{\prime }(b)b-h^{\prime }{(a)a|\le (b-1)}^{p-2}|\frac{1}{a}-\frac{1}{b}|+|(p-2)+\frac{1}{a}||{(b-1)}^{p-2}-{(a-1)}^{p-2}|.\end{array}$$

Graph

For the first term we have

$$\begin{array}{c}\hfill {(b-1)}^{p-2}|\frac{1}{a}-\frac{1}{b}|=\frac{{(b-1)}^{p-2}}{\mathit{ab}}|a-b|\le {(b-1)}^{p-3}|a-b|.\end{array}$$

Graph

For the second term we estimate

$$\begin{array}{cc}& |(p-2)+\frac{1}{a}||{(b-1)}^{p-2}-{(a-1)}^{p-2}|\hfill \\ \hfill & \le p|{(b-1)}^{p-2}-{(a-1)}^{p-2}|\hfill \\ \hfill & \le p|p-2|\underset{t\in [a,b]}{max}{(t-1)}^{p-3}|b-a|\hfill \\ \hfill & \le {p|p-2|[(a-1)}^{p-3}+{(b-1)}^{p-3}]|b-a|.\hfill \end{array}$$

Graph

This completes the proof of the lemma. $\square$

Lemma 2.6

For any $t\in {\mathbb{R}}_{+}$ we have

$$\begin{array}{c}\hfill g{(t)}^2\le h(t){(t-1)}_{+}^p\end{array}$$

Graph

and

$$\begin{array}{c}\hfill g{(t)}^2+g^{\prime }{(t)}^2{t}^2\le \frac{p^2}{p-1}[h(t)+h^{\prime }(t)t]{(t-1)}_{+}^p.\end{array}$$

Graph

Proof

The first assertion can be achieved, since for $t\in {\mathbb{R}}_{+}$ we have

$$\begin{array}{cc}\hfill g{(t)}^2& =\frac{{(t-1)}_{+}^{2p}}{t^2}\le \frac{{(t-1)}_{+}^{2p-1}}{t}=h(t){(t-1)}_{+}^p.\hfill \end{array}$$

Graph

For the second claim, we first observe that both sides are zero for $t\le 1$ . Therefore, it remains to consider $t>1$ . Recalling (2.4) we compute

2.5 $$\begin{array}{c}\hfill h(t)+h^{\prime }(t)t=(p-1){(t-1)}^{p-2}.\end{array}$$

Graph

Moreover, we have

$$\begin{array}{c}\hfill g^{\prime }(t)=\frac{{(t-1)}^{p-1}[(p-1)t+1]}{t^2},\end{array}$$

Graph

so that

$$\begin{array}{cc}\hfill g{(t)}^2+g^{\prime }{(t)}^2{t}^2& =\frac{{(t-1)}^{2p}}{t^2}+\frac{{(t-1)}^{2p-2}{[(p-1)t+1]}^2}{t^2}\hfill \\ \hfill & =\frac{{(t-1)}^{2p-2}}{t^2}[{(t-1)}^2+[(p-1)t+1{]}^2]\hfill \\ \hfill & \le \frac{{(t-1)}^{2p-2}}{t^2}[(t-1)+[(p-1)t+1]{]}^2\hfill \\ \hfill & =p^2{(t-1)}^{2p-2}.\hfill \end{array}$$

Graph

Dividing the previous inequality by $h(t)+h^{\prime }(t)t$ we infer that

$$\begin{array}{cc}\hfill \frac{g{(t)}^2+g^{\prime }{(t)}^2{t}^2}{h(t)+h^{\prime }(t)t}& \le \frac{p^2{(t-1)}^{2p-2}}{(p-1){(t-1)}^{p-2}}=\frac{p^2}{p-1}{(t-1)}^p\hfill \end{array}$$

Graph

holds, proving the second claimed inequality. $\square$

Bilinear forms

For $\epsilon \in [0,1]$ we define

$$\begin{array}{c}\hfill h_{\epsilon }(t):=h(t)+\epsilon ,\text{for}t\in {\mathbb{R}}_{\ge 0}.\end{array}$$

Graph

We note that $h_{\epsilon }\in C^1({\mathbb{R}}_{\ge 0})$ for $p>2$ , while for $p\le 2$ , we have $h_{\epsilon }\in W^{1,1}({\mathbb{R}}_{>0})\cap C^1([0,1)\cap (1,\infty ))$ . For $p=2$ we additionally have $h_{\epsilon }\in W^{1,\infty }({\mathbb{R}}_{\ge 0})$ . In any case, $h_{\epsilon }\equiv \epsilon$ on the interval [0, 1]. Moreover, we let

2.6 $$\begin{array}{c}\hfill {\mathbf{A}}_{\epsilon }(\xi ):=h_{\epsilon }(|\xi |)\xi ,\text{for}\xi \in {\mathbb{R}}^{\mathit{Nn}}.\end{array}$$

Graph

For $\xi \in {\mathbb{R}}^{\mathit{Nn}}\backslash \{0\}$ with $|\xi |\ne 1$ if $1<p<2$ we define the bilinear forms

$$\begin{array}{c}\hfill {\mathcal{A}}_{\epsilon }(\xi )(\eta ,\zeta ):=h_{\epsilon }(|\xi |)\eta ·\zeta +h_{\epsilon }^{\prime }(|\xi |)|\xi |\sum _{i,j=1}^N\sum _{\alpha ,\beta ,\gamma =1}^n\frac{{\xi }_{\alpha }^i{\eta }_{\alpha \gamma }^i{\xi }_{\beta }^j{\zeta }_{\beta \gamma }^j}{{|\xi |}^2}\text{for}\eta ,\zeta \in {\mathbb{R}}^{N{n}^2}\end{array}$$

Graph

and

2.7 $$\begin{array}{c}\hfill {\mathcal{B}}_{\epsilon }(\xi )(\eta ,\zeta ):=h_{\epsilon }(|\xi |)\eta ·\zeta +h_{\epsilon }^{\prime }(|\xi |)|\xi |\sum _{i,j=1}^N\sum _{\alpha ,\beta =1}^n\frac{{\xi }_{\alpha }^i{\eta }_{\alpha }^i{\xi }_{\beta }^j{\zeta }_{\beta }^j}{{|\xi |}^2}\text{for}\eta ,\zeta \in {\mathbb{R}}^{\mathit{Nn}}\end{array}$$

Graph

and

2.8 $$\begin{array}{c}\hfill {\mathcal{C}}_{\epsilon }(\xi )(\eta ,\zeta ):=h_{\epsilon }(|\xi |)\eta ·\zeta +h_{\epsilon }^{\prime }(|\xi |)|\xi |\sum _{i=1}^{N_1}\sum _{\alpha ,\beta =1}^n\frac{{\xi }_{\alpha }^i{\eta }_{\alpha }{\xi }_{\beta }^i{\zeta }_{\beta }}{{|\xi |}^2}\text{for}\eta ,\zeta \in {\mathbb{R}}^n.\end{array}$$

Graph

Observe that all forms are symmetric in the arguments $\eta$ and $\zeta$ . Due to the special structure of $h_{\epsilon }$ and $h_{\epsilon }^{\prime }$ the compositions ${\mathcal{A}}_{\epsilon }(Dv)$ , ${\mathcal{B}}_{\epsilon }(Dv)$ and ${\mathcal{C}}_{\epsilon }(Dv)$ are well definined for $v\in W^{1,p}$ as integrable functions. Therefore integral calculations involving these quantities make sense.

The next Lemma provides the relevant ellipticity and boundedness properties of the bilinear forms ${\mathcal{A}}_{\epsilon }(\xi )$ , ${\mathcal{B}}_{\epsilon }(\xi )$ and ${\mathcal{C}}_{\epsilon }(\xi )$ . The following abbreviations

$$\begin{array}{cc}\hfill \mathit{\lambda }(t)& :=min\{h(t),(p-1){(t-1)}^{p-2}\}\text{for}t>1,\hfill \end{array}$$

Graph

and

$$\begin{array}{cc}\hfill \mathbf{\Lambda }(t)& :=max\{h(t),(p-1){(t-1)}^{p-2}\}\text{for}t>1,\hfill \end{array}$$

Graph

and $\mathit{\lambda }(t)=0=\mathbf{\Lambda }(t)$ for $0\le t\le 1$ prove to be useful in the formulation of the Lemma.

Lemma 2.7

Let $\epsilon \in [0,1]$ and $\xi \in {\mathbb{R}}^{\mathit{Nn}}\backslash \{0\}$ . The bilinear form ${\mathcal{A}}_{\epsilon }(\xi )$ defined above satisfies

$$\begin{array}{c}\hfill [\epsilon +{\mathit{\lambda }(|\xi |)]|\zeta |}^2\le {\mathcal{A}}_{\epsilon }(\xi )(\zeta ,\zeta )\le [\epsilon +{\mathbf{\Lambda }(|\xi |)]|\zeta |}^2\text{for}\text{any}\zeta \in {\mathbb{R}}^{N{n}^2}.\end{array}$$

Graph

The analogous estimates hold for the bilinear form ${\mathcal{B}}_{\epsilon }$ and any $\eta ,\zeta \in {\mathbb{R}}^{\mathit{Nn}}$ , as well as for ${\mathcal{C}}_{\epsilon }$ and any $\eta ,\zeta \in {\mathbb{R}}^n$ .

Proof

If $|\xi |\le 1$ the inequality holds trivially. Therefore, it remains to consider the case $|\xi |>1$ . We first establish the lower bound. If $h_{\epsilon }^{\prime }(|\xi |)\ge 0$ we omit the second term in the definition of ${\mathcal{A}}_{\epsilon }$ and obtain

$$\begin{array}{cc}\hfill {\mathcal{A}}_{\epsilon }(\xi )(\zeta ,\zeta )& \ge h_{\epsilon }{(|\xi |)|\zeta |}^2\ge [\epsilon +{\mathit{\lambda }(|\xi |)]|\zeta |}^2,\hfill \end{array}$$

Graph

while for $h_{\epsilon }^{\prime }(|\xi |)<0$ we use Cauchy-Schwarz inequality and (2.5) to conclude

$$\begin{array}{cc}\hfill {\mathcal{A}}_{\epsilon }(\xi )(\zeta ,\zeta )& \ge [\epsilon +h(|\xi |)+h^{\prime }{(|\xi |)|\xi |]|\zeta |}^2\hfill \\ \hfill & =[\epsilon +(p-1){(|\xi |-1)}^{p-2}{]|\zeta |}^2\ge [\epsilon +{\mathit{\lambda }(|\xi |)]|\zeta |}^2.\hfill \end{array}$$

Graph

Now, we turn our attention to the upper bound. If $h_{\epsilon }^{\prime }(|\xi |)\ge 0$ we have

$$\begin{array}{cc}\hfill {\mathcal{A}}_{\epsilon }(\xi )(\zeta ,\zeta )& \le [h_{\epsilon }(|\xi |)+h_{\epsilon }^{\prime }{(|\xi |)|\xi |]|\zeta |}^2\hfill \\ \hfill & =[\epsilon +(p-1){(|\xi |-1)}^{p-2}{]|\zeta |}^2\le [\epsilon +{\mathbf{\Lambda }(|\xi |)]|\zeta |}^2,\hfill \end{array}$$

Graph

while for $h_{\epsilon }^{\prime }(|\xi |)<0$ we obtain

$$\begin{array}{c}\hfill {\mathcal{A}}_{\epsilon }(\xi )(\zeta ,\zeta )\le h_{\epsilon }{(|\xi |)|\zeta |}^2\le [\epsilon +{\mathbf{\Lambda }(|\xi |)]|\zeta |}^2.\end{array}$$

Graph

This proves the claim for the bilinear form ${\mathcal{A}}_{\epsilon }$ . The corresponding estimates for ${\mathcal{B}}_{\epsilon }$ and ${\mathcal{C}}_{\epsilon }$ follow in the same way. $\square$

It should also be mentioned that the coercive symmetric bilinear forms fulfill Cauchy-Schwarz inequality. In particular, we have

$$\begin{array}{c}\hfill |{\mathcal{C}}_{\epsilon }(\xi )(\eta ,\zeta )|\le \sqrt{{\mathcal{C}}_{\epsilon }(\xi )(\eta ,\eta )}\sqrt{{\mathcal{C}}_{\epsilon }(\xi )(\zeta ,\zeta )}\text{for}\text{any}\eta ,\zeta \in {\mathbb{R}}^n.\end{array}$$

Graph

In the next Lemma we put together the monotonicity and growth properties of the vector field ${\mathbf{A}}_{\epsilon }$ .

Lemma 2.8

Let $\epsilon \in [0,1]$ and $\xi ,\tilde{\xi }\in {\mathbb{R}}^k$ with $|\xi |>1$ . Then, we have

$$\begin{array}{c}\hfill |{\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi )|\le c(p)[\epsilon +\frac{[(|\xi |-1)+(|\tilde{\xi }{|-1)}_{+}{]}^{p-1}}{|\xi |-1}]|\tilde{\xi }-\xi |\end{array}$$

Graph

and

$$\begin{array}{cc}\hfill ({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·(\tilde{\xi }-\xi )& \ge [\epsilon +\frac{min\{1,p-1\}}{2^{p+1}}\frac{{(|\xi |-1)}^p}{|\xi |(|\xi |+|\tilde{\xi }|)}]{|\tilde{\xi }-\xi |}^2.\hfill \end{array}$$

Graph

Proof

The first inequality results from the following chain of inequalities

$$\begin{array}{cc}\hfill |{\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi )|& \le h_{\epsilon }(|\xi |)|\tilde{\xi }-\xi |+|h(|\tilde{\xi }|)-h(|\xi |)||\tilde{\xi }|\hfill \\ \hfill & \le h_{\epsilon }(|\xi |)|\tilde{\xi }-\xi |+c\frac{[(|\xi |-1)+(|\tilde{\xi }{|-1)}_{+}{]}^{p-1}}{|\xi |-1}||\tilde{\xi }|-|\xi ||\hfill \\ \hfill & \le c(p)[\epsilon +\frac{[(|\xi |-1)+(|\tilde{\xi }{|-1)}_{+}{]}^{p-1}}{|\xi |-1}]|\tilde{\xi }-\xi |.\hfill \end{array}$$

Graph

From the third last line to the second last we used Lemma 2.4. For the proof of the second inequality we abbreviate

$$\begin{array}{c}\hfill {\xi }_s:=\xi +s(\tilde{\xi }-\xi ),\text{for}s\in [0,1].\end{array}$$

Graph

Keeping this in mind we compute

$$\begin{array}{cc}\hfill ({\mathbf{A}}_{\epsilon }(\tilde{\xi })& -{\mathbf{A}}_{\epsilon }(\xi ))·(\tilde{\xi }-\xi )\hfill \\ \hfill & ={\int }_0^1\frac{\mathrm{d}}{\mathrm{d}s}{\mathbf{A}}_{\epsilon }({\xi }_s)\mathrm{d}s·(\tilde{\xi }-\xi )\hfill \\ \hfill & ={\int }_0^1[h_{\epsilon }(|{\xi }_s|)(\tilde{\xi }-\xi )+\frac{h_{\epsilon }^{\prime }(|{\xi }_s|)}{|{\xi }_s|}{\xi }_s·(\tilde{\xi }-\xi ){\xi }_s]\mathrm{d}s·(\tilde{\xi }-\xi )\hfill \\ \hfill & ={\int }_0^1{\mathcal{B}}_{\epsilon }({\xi }_s)(\tilde{\xi }-\xi ,\tilde{\xi }-\xi )\mathrm{d}s\hfill \\ \hfill & \ge {\int }_0^1[\epsilon +\mathit{\lambda }(|{\xi }_s|)]\mathrm{d}s|\tilde{\xi }{-\xi |}^2\hfill \\ \hfill & \ge [\epsilon +\frac{min\{1,p-1\}}{|\xi |+|\tilde{\xi }|}{\int }_0^1(|{\xi }_s{|-1)}_{+}^{p-1}\mathrm{d}s]{|\tilde{\xi }-\xi |}^2.\hfill \end{array}$$

Graph

In turn we used Lemma 2.7 and the elementary inequality $(|{\xi }_s{|-1)}_{+}\le |{\xi }_s|\le |\xi |+|\tilde{\xi }|$ . Now, we distinguish whether or not $|\tilde{\xi }|\le |\xi |$ . If $|\tilde{\xi }|\le |\xi |$ , then

$$\begin{array}{c}\hfill |{\xi }_s|\ge (1-s)|\xi |-s|\tilde{\xi }|\ge (1-2s)|\xi |>1\forall s\in [0,\frac{|\xi |-1}{2|\xi |}).\end{array}$$

Graph

For $s\in [0,\frac{|\xi |-1}{4|\xi |}]$ this implies a bound from below in the form

$$\begin{array}{cc}\hfill (|{\xi }_s{|-1)}_{+}& =|{\xi }_s|-1\ge (1-2s)|\xi |-1\ge [1-\frac{|\xi |-1}{2|\xi |}]|\xi |-1=\frac{1}{2}(|\xi |-1).\hfill \end{array}$$

Graph

So it follows that

$$\begin{array}{c}\hfill {\int }_0^1(|{\xi }_s{|-1)}_{+}^{p-1}\mathrm{d}s\ge {\int }_0^{\frac{|\xi |-1}{4|\xi |}}(|{\xi }_s{|-1)}_{+}^{p-1}\mathrm{d}s=\frac{1}{2^{p+1}}\frac{{(|\xi |-1)}^p}{|\xi |}\end{array}$$

Graph

holds true. In the case that $|\tilde{\xi }|>|\xi |$ , we estimate $|{\xi }_s|$ from below by

$$\begin{array}{c}\hfill |{\xi }_s|\ge s|\tilde{\xi }|-(1-s)|\xi |>(2s-1)|\xi |>1\forall s\in (\frac{|\xi |+1}{2|\xi |},1].\end{array}$$

Graph

Therefore, for $s\in [\frac{3|\xi |+1}{4|\xi |},1]$ we obtain

$$\begin{array}{c}\hfill (|{\xi }_s{|-1)}_{+}=|{\xi }_s|-1\ge (2s-1)|\xi |-1\ge [\frac{3|\xi |+1}{2|\xi |}-1]|\xi |-1=\frac{1}{2}(|\xi |-1).\end{array}$$

Graph

This yields

$$\begin{array}{c}\hfill {\int }_0^1(|{\xi }_s{|-1)}_{+}^{p-1}\mathrm{d}s\ge {\int }_{\frac{3|\xi |+1}{4|\xi |}}^1(|{\xi }_s{|-1)}_{+}^{p-1}\mathrm{d}s\ge \frac{1}{2^{p+1}}\frac{{(|\xi |-1)}^p}{|\xi |}.\end{array}$$

Graph

Inserting this above, we obtain the second claim of the Lemma. $\square$

Lemma 2.9

Let $\epsilon ,\delta \in (0,1]$ and $\xi ,\tilde{\xi }\in {\mathbb{R}}^{\mathit{Nn}}$ . Then, we have

$$\begin{array}{c}\hfill {\epsilon }^{\frac{1}{2}}|\xi -\tilde{\xi }{|}^2+|{\mathcal{G}}_{\delta }(\xi )-{\mathcal{G}}_{\delta }(\tilde{\xi }){|}^p\le {\epsilon }^{\frac{1}{2}}{|\xi |}^2+c{\epsilon }^{-\frac{1}{2}}({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·(\tilde{\xi }-\xi )\end{array}$$

Graph

for a constant $c=c(p,\delta )$ .

Proof

If $|\xi |,|\tilde{\xi }|\le 1+\delta$ , we have ${\mathcal{G}}_{\delta }(\xi )=0={\mathcal{G}}_{\delta }(\tilde{\xi })$ . Therefore, the desired inequality follows from the second inequality in Lemma 2.8 after omitting the positive second term in the bracket on the right-hand side.

Therefore, it remains to consider the case where either $|\xi |>1+\delta$ or $|\tilde{\xi }|>1+\delta$ . We again distinguish two cases and start with $|\xi |\ge |\tilde{\xi }|$ . Note that this implies $|\xi |>1+\delta$ . If $p\ge 2$ we use Lemma 2.3 to conclude

$$\begin{array}{c}\hfill |{\mathcal{G}}_{\delta }(\tilde{\xi })-{\mathcal{G}}_{\delta }(\xi ){|}^p\le 3^p|\tilde{\xi }{-\xi |}^p\le 2^{p-2}{3}^p{|\xi |}^{p-2}{|\tilde{\xi }-\xi |}^2,\end{array}$$

Graph

while in the case $p<2$ we use Young's inequality and Lemma 2.3 to obtain

$$\begin{array}{cc}\hfill |{\mathcal{G}}_{\delta }(\tilde{\xi })-{\mathcal{G}}_{\delta }(\xi ){|}^p& ={|\xi |}^{\frac{p(2-p)}{2}}{|\xi |}^{\frac{p(p-2)}{2}}|{\mathcal{G}}_{\delta }(\tilde{\xi })-{\mathcal{G}}_{\delta }(\xi ){|}^p\hfill \\ \hfill & \le \frac{1}{2}{\epsilon }^{\frac{1}{2}}{|\xi |}^p+c{\epsilon }^{-\frac{2-p}{2p}}{|\xi |}^{p-2}|{\mathcal{G}}_{\delta }(\tilde{\xi })-{\mathcal{G}}_{\delta }(\xi ){|}^2\hfill \\ \hfill & \le \frac{1}{2}{\epsilon }^{\frac{1}{2}}{|\xi |}^2+c{\epsilon }^{-\frac{1}{2}}{|\xi |}^{p-2}{|\tilde{\xi }-\xi |}^2,\hfill \end{array}$$

Graph

for a constant $c=c(p)$ . Combining both cases, taking into account the elementary inequalities $\frac{1}{|\xi |}\le \frac{2}{|\xi |+|\tilde{\xi }|}$ and $|\xi |\le (1+\frac{1}{\delta })(|\xi |-1)$ , and finally applying Lemma 2.8, we obtain

$$\begin{array}{cc}\hfill {\epsilon }^{\frac{1}{2}}{|\tilde{\xi }-\xi |}^2+|{\mathcal{G}}_{\delta }(\xi )-{\mathcal{G}}_{\delta }(\tilde{\xi }){|}^p& \le {\epsilon }^{\frac{1}{2}}|\tilde{\xi }{-\xi |}^2+\frac{1}{2}{\epsilon }^{\frac{1}{2}}{|\xi |}^2+\frac{c{(|\xi |-1)}_{+}^p}{{\epsilon }^{\frac{1}{2}}|\xi |(|\tilde{\xi }|+|\xi |)}|\tilde{\xi }-\xi {|}^2\hfill \\ \hfill & \le \frac{1}{2}{\epsilon }^{\frac{1}{2}}{|\xi |}^2+c{\epsilon }^{-\frac{1}{2}}({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·(\tilde{\xi }-\xi ),\hfill \end{array}$$

Graph

where $c=c(p,\delta )$ . This proves the claimed inequality in this case. In the remaning case, i.e. $|\tilde{\xi }|>|\xi |$ , we obtain a similar estimate. We only have to replace on the right-hand side ${|\xi |}^2$ by $|\tilde{\xi }{|}^2$ . Then, we use the estimate $|\tilde{\xi }{|}^2\le 2(|\tilde{\xi }{-\xi |}^2+{|\xi |}^2)$ and absorb ${\epsilon }^{\frac{1}{2}}{|\tilde{\xi }-\xi |}^2$ by Lemma 2.8 into the second term on the right-hand side. In this way, we obtain

$$\begin{array}{cc}\hfill {\epsilon }^{\frac{1}{2}}{|\tilde{\xi }-\xi |}^2+|{\mathcal{G}}_{\delta }(\xi )-{\mathcal{G}}_{\delta }(\tilde{\xi }){|}^p& \le \frac{1}{2}{\epsilon }^{\frac{1}{2}}{|\tilde{\xi }|}^2+c{\epsilon }^{-\frac{1}{2}}({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·(\tilde{\xi }-\xi )\hfill \\ \hfill & \le {\epsilon }^{\frac{1}{2}}{|\xi |}^2+c{\epsilon }^{-\frac{1}{2}}({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·(\tilde{\xi }-\xi ),\hfill \end{array}$$

Graph

proving the claim also in this case. $\square$

In the following Lemma we quantify the remainder term in the linearization of ${\mathbf{A}}_{\epsilon }$ . In the application, it can be assumed that the linearization only takes place in points $\xi$ with $|\xi |>1$ in a quantifiable way. The precise statement is

Lemma 2.10

Let $\epsilon \in [0,1]$ , and $\xi ,\tilde{\xi }\in {\mathbb{R}}^{\mathit{Nn}}$ with $|\xi |\ge 1+\frac{1}{4}\mu$ and $|\xi |,|\tilde{\xi }|\le 1+2\mu$ for some $\mu >0$ . Then, we have

$$\begin{array}{c}\hfill |{\mathcal{B}}_{\epsilon }(\xi )(\tilde{\xi }-\xi ,\zeta )-({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·\zeta |\le c(p){\mu }^{p-3}|\tilde{\xi }{-\xi |}^2|\zeta |,\forall \zeta \in {\mathbb{R}}^{\mathit{Nn}}.\end{array}$$

Graph

Proof

We distinguish two cases. We start with the case $|\xi -\tilde{\xi }|\le \frac{1}{8}\mu$ . For $s\in [0,1]$ we write ${\xi }_s:=\xi +s(\tilde{\xi }-\xi )$ . Note that

2.9 $$\begin{array}{c}\hfill |{\xi }_s|\ge |\xi |-s|\tilde{\xi }-\xi |\ge 1+\frac{1}{4}\mu -\frac{1}{8}\mu =1+\frac{1}{8}\mu \forall s\in [0,1].\end{array}$$

Graph

Similarly to the computations in the proof of Lemma 2.8 we have

$$\begin{array}{cc}\hfill ({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·\zeta & ={\int }_0^1\frac{\mathrm{d}}{\mathrm{d}s}{\mathbf{A}}_{\epsilon }(\xi +s(\tilde{\xi }-\xi ))·\zeta \mathrm{d}s\hfill \\ \hfill & ={\int }_0^1[h_{\epsilon }(|{\xi }_s|)(\tilde{\xi }-\xi )+\frac{h_{\epsilon }^{\prime }(|{\xi }_s|)}{|{\xi }_s|}{\xi }_s·(\tilde{\xi }-\xi ){\xi }_s]·\zeta \mathrm{d}s\hfill \\ \hfill & ={\int }_0^1{\mathcal{B}}_{\epsilon }({\xi }_s)(\tilde{\xi }-\xi ,\zeta )\mathrm{d}s.\hfill \end{array}$$

Graph

This allows us to re-write

$$\begin{array}{cc}\hfill |{\mathcal{B}}_{\epsilon }(\xi )(\tilde{\xi }-\xi ,\zeta )-({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·\zeta |& =|{\int }_0^1[{\mathcal{B}}_{\epsilon }(\xi )-{\mathcal{B}}_{\epsilon }({\xi }_s)](\tilde{\xi }-\xi ,\zeta )\mathrm{d}s|.\hfill \end{array}$$

Graph

We decompose and estimate the integrand appearing on the right-hand side and obtain in this way

$$\begin{array}{cc}\hfill |[{\mathcal{B}}_{\epsilon }(\xi )-{\mathcal{B}}_{\epsilon }({\xi }_s)](\tilde{\xi }-\xi ,\zeta )|& \le \mathbf{I}+\mathbf{II},\hfill \end{array}$$

Graph

where

$$\begin{array}{cc}\hfill \mathbf{I}& :=|h(|\xi |)-h(|{\xi }_s|)||\tilde{\xi }-\xi ||\zeta |\hfill \\ \hfill \mathbf{II}& :=|(h^{\prime }(|\xi |)|\xi |\frac{{\xi }_{\alpha }^i{\xi }_{\beta }^j}{{|\xi |}^2}-h^{\prime }(|{\xi }_s|)|{\xi }_s|\frac{{({\xi }_s)}_{\alpha }^i{({\xi }_s)}_{\beta }^j}{|{\xi }_s{|}^2}){(\tilde{\xi }-\xi )}_{\alpha }^i{\zeta }_{\beta }^j|\hfill \end{array}$$

Graph

In view of Lemma 2.4 we have

$$\begin{array}{cc}\hfill |h_{\epsilon }(|\xi |)-h_{\epsilon }(|{\xi }_s|)|& \le c(p)\frac{[|\xi |-1+(|{\xi }_s{|-1)}_{+}{]}^{p-1}}{|{\xi }_s|(|\xi |-1)}|\xi -{\xi }_s|\hfill \\ \hfill & \le c(p){\mu }^{p-3}|\xi -\tilde{\xi }|,\hfill \end{array}$$

Graph

which immediately implies

$$\begin{array}{c}\hfill \mathbf{I}\le c(p){\mu }^{p-3}|\xi -\tilde{\xi }{|}^2|\zeta |.\end{array}$$

Graph

For the second term we use Lemma 2.1, Lemma 2.5, the assumptions on $\xi ,\tilde{\xi }$ and (2.9) to obtain

$$\begin{array}{cc}\hfill \mathbf{II}& \le [\sum _{\alpha ,\beta =1}^n\sum _{i,j=1}^N(h^{\prime }(|\xi |)|\xi |\frac{{\xi }_{\alpha }^i{\xi }_{\beta }^j}{{|\xi |}^2}-h^{\prime }(|{\xi }_s|)|{\xi }_s|\frac{{({\xi }_s)}_{\alpha }^i{({\xi }_s)}_{\beta }^j}{|{\xi }_s{|}^2}{)}^2{]}^{\frac{1}{2}}|\tilde{\xi }-\xi ||\zeta |\hfill \\ \hfill & =|h^{\prime }(|\xi |)|\xi |\frac{\xi \otimes \xi }{{|\xi |}^2}-h^{\prime }(|{\xi }_s|)|{\xi }_s|\frac{{\xi }_s\otimes {\xi }_s}{|{\xi }_s{|}^2}||\tilde{\xi }-\xi ||\zeta |\hfill \\ \hfill & \le [|h^{\prime }(|\xi |)||\xi ||\frac{\xi \otimes \xi }{{|\xi |}^2}-\frac{{\xi }_s\otimes {\xi }_s}{|{\xi }_s{|}^2}|+|h^{\prime }(|\xi |)|\xi |-h^{\prime }(|{\xi }_s|)|{\xi }_s||]|\tilde{\xi }-\xi ||\zeta |\hfill \\ \hfill & \le [2|h^{\prime }(|\xi |)|\frac{|\xi |+|{\xi }_s|}{|\xi |}|\xi -{\xi }_s|+|h^{\prime }(|\xi |)|\xi |-h^{\prime }(|{\xi }_s|)|{\xi }_s||]|\tilde{\xi }-\xi ||\zeta |\hfill \\ \hfill & \le c(p)[{(|\xi |-1)}^{p-3}\frac{|\xi |+|{\xi }_s|}{|\xi |}+{(|\xi |-1)}^{p-3}+(|{\xi }_s{|-1)}^{p-3}]|\tilde{\xi }{-\xi |}^2|\zeta |\hfill \\ \hfill & \le c(p){\mu }^{p-3}|\tilde{\xi }{-\xi |}^2|\zeta |.\hfill \end{array}$$

Graph

Inserting the preceding estimates above, we find that

$$\begin{array}{cc}\hfill |{\mathcal{B}}_{\epsilon }(\xi )(\tilde{\xi }-\xi ,\zeta )-({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·\zeta |& \le c(p){\mu }^{p-3}|\tilde{\xi }{-\xi |}^2|\zeta |.\hfill \end{array}$$

Graph

At this stage it remains to consider the case $|\xi -\tilde{\xi }|>\frac{1}{8}\mu$ . Note that

$$\begin{array}{c}\hfill {\mathcal{B}}_{\epsilon }(\xi )(\tilde{\xi }-\xi ,\zeta )=\epsilon (\tilde{\xi }-\xi )·\zeta +{\mathcal{B}}_0(\xi )(\tilde{\xi }-\xi ,\zeta ).\end{array}$$

Graph

Therefore, by Lemma 2.7 and the assumption $\mu <8|\tilde{\xi }-\xi |$ we obtain

$$\begin{array}{cc}& |{\mathcal{B}}_{\epsilon }(\xi )(\tilde{\xi }-\xi ,\zeta )-({\mathbf{A}}_{\epsilon }(\tilde{\xi })-{\mathbf{A}}_{\epsilon }(\xi ))·\zeta |\hfill \\ \hfill & =|{\mathcal{B}}_0(\xi )(\tilde{\xi }-\xi ,\zeta )-(h(|\tilde{\xi }|)\tilde{\xi }-h(|\xi |)\xi )·\zeta |\hfill \\ \hfill & \le [\mathbf{\Lambda }(|\xi |)|\tilde{\xi }-\xi |+(|\tilde{\xi }{|-1)}_{+}^{p-1}+{(|\xi |-1)}^{p-1}]|\zeta |\hfill \\ \hfill & \le c(p)[{\mu }^{p-2}|\tilde{\xi }-\xi |+{\mu }^{p-1}]|\zeta |\hfill \\ \hfill & \le c(p){\mu }^{p-3}|\tilde{\xi }{-\xi |}^2|\zeta |.\hfill \end{array}$$

Graph

Joining both cases yields the claim. $\square$

Lemma 2.11

For any $\epsilon \in [0,1]$ , any ball $B_R\subset {\mathbb{R}}^n$ and any $v\in W_{\mathrm{loc}}^{2,2}(B_R,{\mathbb{R}}^N)$ we have

$$\begin{array}{c}\hfill |D[g(|Dv|)Dv]{|}^2\le \frac{2p^2}{p-1}{\mathcal{A}}_{\epsilon }(Dv)(D^{2}v,D^{2}v){(|Dv|-1)}_{+}^p\text{a.e.}\text{in}{B}_R,\end{array}$$

Graph

where g is defined in (2.1).

Proof

For $\alpha \in \{1,\cdots ,n\}$ we compute

$$\begin{array}{cc}\hfill |D_{\alpha }[g(|Dv|)Dv]{|}^2& =|g(|Dv|)D_{\alpha }Dv+g^{\prime }(|Dv|)D_{\alpha }|Dv|Dv{|}^2\hfill \\ \hfill & \le 2[g{(|Dv|)}^2|D_{\alpha }Dv{|}^2+g^{\prime }{(|Dv|)}^2{|Dv|}^2|D_{\alpha }|Dv|{|}^2].\hfill \end{array}$$

Graph

Summing with respect to $\alpha \in \{1,\cdots ,n\}$ , we obtain

$$\begin{array}{cc}\hfill |D[g(|Dv|)Dv]{|}^2& =\sum _{\alpha =1}^n|D_{\alpha }[g(|Dv|)Dv]{|}^2\hfill \\ \hfill & \le 2[g{(|Dv|)}^2|D^2{v|}^2+g^{\prime }{(|Dv|)}^2{|Dv|}^2|\mathrm{\nabla }|Dv|{|}^2].\hfill \end{array}$$

Graph

In the preceding inequality we replace the term $g{(|Dv|)}^2$ by $h(|Dv|){(|Dv|-1)}_{+}^p$ , which is possible by an application of the first inequality in Lemma 2.6. Moreover, we would also like to replace $g^{\prime }{(|Dv|)}^2{|Dv|}^2$ by $h^{\prime }{(|Dv|)|Dv|(|Dv|-1)}_{+}^p$ . To this aim we have to distinguish two cases. If $h^{\prime }(|Dv|)\ge 0$ , we use the second inequality in Lemma 2.6 to replace $g^{\prime }{(|Dv|)}^2{|Dv|}^2$ by $[h(|Dv|)+h^{\prime }{(|Dv|)|Dv|](|Dv|-1)}_{+}^p$ . Thereby, we may omit the positive term $g{(|Dv|)}^2$ on the left-hand side. Inserting this above and using also Kato's inequality in the form $|\mathrm{\nabla }|Dv||\le |D^{2}v|$ , we obtain

$$\begin{array}{cc}& |D[g(|Dv|)Dv]{|}^2\hfill \\ \hfill & \le \frac{p^2}{p-1}{(|Dv|-1)}_{+}^p[h(|Dv|)|D^2{v|}^2+[h(|Dv|)+h^{\prime }(|Dv|)|Dv|]|\mathrm{\nabla }|Dv|{|}^2]\hfill \\ \hfill & \le \frac{2p^2}{p-1}{(|Dv|-1)}_{+}^p[h(|Dv|)|D^2{v|}^2+h^{\prime }(|Dv|)|Dv||\mathrm{\nabla }|Dv|{|}^2]\hfill \\ \hfill & \le \frac{2p^2}{p-1}{\mathcal{A}}_{\epsilon }(Dv)(D^{2}v,D^{2}v){(|Dv|-1)}_{+}^p.\hfill \end{array}$$

Graph

Otherwise, if $h^{\prime }(|Dv|)<0$ we use Kato's inequality twice and again the second inequality in Lemma 2.6 to obtain

$$\begin{array}{cc}\hfill |D[g(|Dv|)Dv]{|}^2& \le 2[g{(|Dv|)}^2+g^{\prime }{(|Dv|)}^2{|Dv|}^2]{|D^{2}v|}^2\hfill \\ \hfill & \le \frac{2p^2}{p-1}{(|Dv|-1)}_{+}^p[h(|Dv|)+h^{\prime }(|Dv|)|Dv|]|D^2{v|}^2\hfill \\ \hfill & \le \frac{2p^2}{p-1}{(|Dv|-1)}_{+}^p[h(|Dv|)|D^2{v|}^2+h^{\prime }(|Dv|)|Dv||\mathrm{\nabla }|Dv|{|}^2]\hfill \\ \hfill & \le \frac{2p^2}{p-1}{\mathcal{A}}_{\epsilon }(Dv)(D^{2}v,D^{2}v){(|Dv|-1)}_{+}^p.\hfill \end{array}$$

Graph

This proves the asserted inequality also in the second case. $\square$

Proof of Theorem 1.1

In this section we will prove Theorem 1.1 under the hypothesis that Propositions 3.4 and 3.5 below are true. The remainder of the paper is then devoted to the proof of those two propositions.

Here and in the following we denote by $u\in W^{1,p}(\mathrm{\Omega },{\mathbb{R}}^N)$ a weak solution of (1.5). We first observe that $u\in W_{\mathrm{loc}}^{1,\infty }(\mathrm{\Omega },{\mathbb{R}}^N)$ ; cf [[2], [4]] for the scalar case and [[9]] for the vectorial case. Therefore, we may always assume that Du is locally bounded in $\mathrm{\Omega }$ .

Regularization

The first step in the proof consists in the construction of more regular approximating solutions. To this aim we consider a fixed ball $B_R\equiv B_R(y_o)⋐\mathrm{\Omega }$ . We let $\epsilon \in (0,1]$ and $\mathfrak{p}:=max\{p,2\}$ and recall the definition of the regularized coefficients ${\mathbf{A}}_{\epsilon }$ from (2.6). By $u_{\epsilon }\in u+W_0^{1,\mathfrak{p}}(B_R,{\mathbb{R}}^N)$ we denote the unique weak solution of the regularized elliptic system

3.1 $$\begin{array}{c}\hfill \left\{\begin{array}{cc}\mathrm{div}{\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })=f,& \text{in}{B}_R,\hfill \\ [5pt]u_{\epsilon }=u,& \text{on}\partial B_R.\hfill \end{array}\right)\end{array}$$

Graph

The weak formulation of (3.1) is

3.2 $$\begin{array}{c}\hfill {\int }_{B_R}{\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })·D\phi \mathrm{d}x=-{\int }_{B_R}f·\phi \mathrm{d}x\forall \phi \in W_0^{1,\mathfrak{p}}(B_R,{\mathbb{R}}^N).\end{array}$$

Graph

Note that $u_{\epsilon }\in W_{\mathrm{loc}}^{1,\infty }(B_R,{\mathbb{R}}^N)\cap W_{\mathrm{loc}}^{2,2}(B_R,{\mathbb{R}}^N)$ . This can be retrieved similarly as in [[33]].

Lemma 3.1

For any $\epsilon \in (0,1]$ we have $u_{\epsilon }\in W_{\mathrm{loc}}^{1,\infty }(B_R,{\mathbb{R}}^N)\cap W_{\mathrm{loc}}^{2,2}(B_R,{\mathbb{R}}^N)$ . Moreover, for any ball $B_{2\varrho }(x_o)⋐B_R$ the (uniform with respect to $\epsilon$ ) quantitative $L^{\infty }$ -gradient bound

Graph

with $c=c(n,N,p,\sigma ,\Vert f{\Vert }_{L^{n+\sigma }(B_R)})$ and the quantitative $W^{2,2}$ -estimate

$$\begin{array}{cc}\hfill {\int }_{B_{\varrho }(x_o)}{|D^2{u}_{\epsilon }|}^2\mathrm{d}x& \le \frac{c}{\epsilon {\varrho }^2}[{\int }_{B_{2\varrho }(x_o)}(|D{u}_{\epsilon }{|}^p+\epsilon |D{u}_{\epsilon }{|}^2)\mathrm{d}x+\frac{{\varrho }^2}{\epsilon }{\int }_{B_{2\varrho }(x_o)}{|f|}^2\mathrm{d}x]\hfill \end{array}$$

Graph

with $c=c(n,N,p)$ hold true.

Our first observation is a uniform energy bound for $D{u}_{\epsilon }$ .

Lemma 3.2

There exists a constant $c=c(n,p)$ such that for any $\epsilon \in (0,1]$ we have

$$\begin{array}{cc}\hfill {\int }_{B_R}(|D{u}_{\epsilon }{|}^p+\epsilon {|D{u}_{\epsilon }|}^2)\mathrm{d}x& \le c{\int }_{B_R}{(|Du|}^p+{\epsilon |Du|}^2+1)\mathrm{d}x+c{R}^n{\Vert f\Vert }_{L^n(B_R)}^{\frac{p}{p-1}}.\hfill \end{array}$$

Graph

Proof

The desired estimate can be deduced with a standard argument by testing the weak formulation (3.2) with the test-function $\phi :=u_{\epsilon }-u$ . Indeed, we have

$$\begin{array}{cc}& {\int }_{B_R}[(|D{u}_{\epsilon }{|-1)}_{+}^p+\epsilon {|D{u}_{\epsilon }|}^2]\mathrm{d}x\hfill \\ \hfill & \le {\int }_{B_R}{h}_{\epsilon }(|D{u}_{\epsilon }|)D{u}_{\epsilon }·D{u}_{\epsilon }\mathrm{d}x\hfill \\ \hfill & ={\int }_{B_R}{h}_{\epsilon }(|D{u}_{\epsilon }|)D{u}_{\epsilon }·Du\mathrm{d}x-{\int }_{B_R}f·(u_{\epsilon }-u)\mathrm{d}x.\hfill \end{array}$$

Graph

By Young's inequality we obtain for the first integral on the right-hand side

$$\begin{array}{cc}\hfill {\int }_{B_R}{h}_{\epsilon }(|D{u}_{\epsilon }|)D{u}_{\epsilon }·Du\mathrm{d}x& \le {\int }_{B_R}[(|D{u}_{\epsilon }{|-1)}_{+}^{p-1}|Du|+\epsilon |D{u}_{\epsilon }||Du|]\mathrm{d}x\hfill \\ \hfill & \le \frac{1}{4}{\int }_{B_R}(|D{u}_{\epsilon }{|-1)}_{+}^p\mathrm{d}x+c{\int }_{B_R}{|Du|}^p\mathrm{d}x\hfill \\ \hfill & +\frac{\epsilon }{2}{\int }_{B_R}(|D{u}_{\epsilon }{|}^2+{|Du|}^2)\mathrm{d}x.\hfill \end{array}$$

Graph

The second integral on the right-hand side is estimated with Hölder's and Sobolev's inequality, so that

$$\begin{array}{cc}& |{\int }_{B_R}f·(u_{\epsilon }-u)\mathrm{d}x|\hfill \\ \hfill & \le [{\int }_{B_R}{|f|}^n\mathrm{d}x{]}^{\frac{1}{n}}[{\int }_{B_R}{|u_{\epsilon }-u|}^{\frac{n}{n-1}}\mathrm{d}x{]}^{\frac{n-1}{n}}\hfill \\ \hfill & \le c[{\int }_{B_R}{|f|}^n\mathrm{d}x{]}^{\frac{1}{n}}{\int }_{B_R}|D{u}_{\epsilon }-Du|\mathrm{d}x\hfill \\ \hfill & \le c|B_R{|}^{\frac{p-1}{p}}[{\int }_{B_R}{|f|}^n\mathrm{d}x{]}^{\frac{1}{n}}[{\int }_{B_R}[(|D{u}_{\epsilon }{|-1)}_{+}^p+{|Du|}^p+1]\mathrm{d}x{]}^{\frac{1}{p}}\hfill \\ \hfill & \le c{R}^n{\Vert f\Vert }_{L^n(B_R)}^{\frac{p}{p-1}}+\frac{1}{4}{\int }_{B_R}[(|D{u}_{\epsilon }{|-1)}_{+}^p+{|Du|}^p+1]\mathrm{d}x,\hfill \end{array}$$

Graph

where $c=c(n,p)$ . We insert these inequalities above and reabsorb the terms containing $D{u}_{\epsilon }$ from the right-hand side into the left. In this way, we get

$$\begin{array}{cc}& {\int }_{B_R}[(|D{u}_{\epsilon }{|-1)}_{+}^p+\epsilon {|D{u}_{\epsilon }|}^2]\mathrm{d}x\hfill \\ \hfill & \le c{\int }_{B_R}{[|Du|}^p+{\epsilon |Du|}^2+1]\mathrm{d}x+c{R}^n{\Vert f\Vert }_{L^n(B_R)}^{\frac{p}{p-1}},\hfill \end{array}$$

Graph

with a constant $c=c(n,p)$ . The desired uniform energy bound can easily be deduced from the preceding inequality. $\square$

The next lemma ensures strong convergence of the approximating solutions, in the sense that ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ strongly converges to ${\mathcal{G}}_{\delta }(Du)$ in $L^p$ .

Lemma 3.3

Let $\delta \in (0,1]$ and $u_{\epsilon }$ with $\epsilon \in (0,1]$ be the unique solution of the Dirichlet problem (3.1). Then, we have

$$\begin{array}{c}\hfill {\mathcal{G}}_{\delta }(D{u}_{\epsilon })\to {\mathcal{G}}_{\delta }(Du)\text{in}{L}^p(B_R,{\mathbb{R}}^{\mathit{Nn}})\text{as}\epsilon \downarrow 0.\end{array}$$

Graph

Proof

Testing (3.2) and the weak formulation of (1.5) with $\phi :=u_{\epsilon }-u$ we have

$$\begin{array}{c}\hfill {\int }_{B_R}({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(Du))·(D{u}_{\epsilon }-Du)\mathrm{d}x=\epsilon {\int }_{B_R}Du·(Du-D{u}_{\epsilon })\mathrm{d}x.\end{array}$$

Graph

Using Lemma 2.9 and Young's inequality we find

$$\begin{array}{cc}& {\int }_{B_R}[{\epsilon }^{\frac{1}{2}}{|D{u}_{\epsilon }-Du|}^2+|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })-{\mathcal{G}}_{\delta }(Du){|}^p]\mathrm{d}x\hfill \\ \hfill & \le {\epsilon }^{\frac{1}{2}}{\int }_{B_R}{|Du|}^2\mathrm{d}x+c{\epsilon }^{-\frac{1}{2}}{\int }_{B_R}({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(Du))·(D{u}_{\epsilon }-Du)\mathrm{d}x\hfill \\ \hfill & \le {\epsilon }^{\frac{1}{2}}{\int }_{B_R}{|Du|}^2\mathrm{d}x+c{\epsilon }^{\frac{1}{2}}{\int }_{B_R}|Du||D{u}_{\epsilon }-Du|\mathrm{d}x\hfill \\ \hfill & \le {\epsilon }^{\frac{1}{2}}{\int }_{B_R}|D{u}_{\epsilon }{-Du|}^2\mathrm{d}x+c{\epsilon }^{\frac{1}{2}}{\int }_{B_R}{|Du|}^2\mathrm{d}x.\hfill \end{array}$$

Graph

Re-absorbing the first integral from the right into the left-hand side, we obtain

$$\begin{array}{cc}\hfill {\int }_{B_R}|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })-{\mathcal{G}}_{\delta }(Du){|}^p\mathrm{d}x& \le c{\epsilon }^{\frac{1}{2}}{\int }_{B_R}{|Du|}^2\mathrm{d}x.\hfill \end{array}$$

Graph

The integral on the right-hand side is finite, since $Du\in L_{\mathrm{loc}}^{\infty }(\mathrm{\Omega },{\mathbb{R}}^{\mathit{Nn}})$ . Therefore, the preceding inequality implies strong convergence ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })\to {\mathcal{G}}_{\delta }(Du)$ in $L^p(B_R,{\mathbb{R}}^{\mathit{Nn}})$ as $\epsilon \downarrow 0$ . $\square$

Hölder-continuity of Gδ(Duε)

We recall the definition of ${\mathcal{G}}_{\delta }$ from (2.3). In this subsection we will prove that ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ is locally Hölder continuous in $B_R$ for any $\delta \in (0,1]$ . This will be achieved in Theorem 3.6. Thereby, it is essential that all constants are independent of $\epsilon$ .

The proof of Theorem 3.6 relies on the distinction between two different regimes—the degenerate and non-degenerate regime. In the non-degenerate regime we will prove an excess-decay estimate for ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ (see Proposition 3.4), while in the degenerate regime we establish a reduction of the modulus of ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ (see Proposition 3.5). The precise setup is as follows. We consider a ball $B_R⋐\mathrm{\Omega }$ and denote by $u_{\epsilon }$ the unique weak solution of the Cauchy–Dirichlet problem (3.1). For $0<r_o<R$ we let $r_1:=\frac{1}{2}(R+r_o)$ . Then, by Lemma 3.1 and Lemma 3.2 we have (uniform with respect to $\epsilon$ ) boundedness of $D{u}_{\epsilon }$ on $B_{r_1}$ . More precisely, there exists a constant

3.3 $$\begin{array}{c}\hfill M=M(n,N,p,R-r_o,{\Vert Du\Vert }_{L^p(B_R)},...)\end{array}$$

Graph

independent of $\epsilon >0$ , such that $\Vert D{u}_{\epsilon }{\Vert }_{L^{\infty }(B_{r_1})}\le M$ . We may assume that $M\ge 3$ . Now, we consider a center $x_o\in B_{r_o}$ and a radius $\varrho \le r_1$ such that $B_{2\varrho }(x_o)\subset B_{r_1}$ . On this ball we have

3.4 $$\begin{array}{c}\hfill \underset{B_{2\varrho }(x_o)}{sup}|D{u}_{\epsilon }|\le 1+\delta +\mu \end{array}$$

Graph

for some $\mu >0$ such that

3.5 $$\begin{array}{c}\hfill 1+\delta +\mu \le M.\end{array}$$

Graph

Next, for $\nu \in (0,1)$ we define the super-level set of $|D{u}_{\epsilon }|$ by

$$\begin{array}{c}\hfill E_{\varrho }^{\nu }(x_o):=\{x\in B_{\varrho }(x_o):|D{u}_{\epsilon }(x)|-1-\delta >(1-\nu )\mu \}.\end{array}$$

Graph

The definition of the super-level set allows us to distinguish between the degenerate regime which is characterized by the measure condition $|B_{\varrho }(x_o)\backslash E_{\varrho }^{\nu }(x_o)|\ge \nu |B_{\varrho }(x_o)|$ and the non-degenerate regime which is characterized by the reversed inequality. Roughly speaking, in the degenerate regime the set of points with $|D{u}_{\epsilon }|$ small has large measure, while in the non-degenerate regime the set of points with $|D{u}_{\epsilon }|$ small is small in measure. We start with the latter one. In the following we abbreviate $\beta :=\frac{\sigma }{n+\sigma }$ .

Proposition 3.4

Let $\epsilon ,\delta \in (0,1]$ and

3.6 $$\begin{array}{c}\hfill 0<\delta <\mu .\end{array}$$

Graph

Then, there exist $\nu ={\nu (n,N,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta )\in (0,\frac{1}{4}]$ and $\widehat{\varrho }=\widehat{\varrho }(n,N,p,\sigma ,$ ${\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta )\in (0,1]$ such that there holds: Whenever $B_{2\varrho }(x_o)\subset B_{r_1}$ is a ball with radius $\varrho \le \widehat{\varrho }$ and center $x_o\in B_{r_o}$ , and $u_{\epsilon }$ is the unique weak solution of the Dirichlet problem (3.1) and hypothesis (3.4) and (3.5) and the measure condition

3.7 $$\begin{array}{c}\hfill |B_{\varrho }(x_o)\backslash E_{\varrho }^{\nu }(x_o)|<\nu |B_{\varrho }(x_o)|\end{array}$$

Graph

are satisfied, then the limit

3.8 $$\begin{array}{c}\hfill {\mathrm{\Gamma }}_{x_o}:=\underset{r\downarrow 0}{lim}({\mathcal{G}}_{2\delta }(D{u}_{\epsilon }){)}_{x_o,r}\end{array}$$

Graph

exists, and the excess decay estimate

3.9

Graph

holds true. Moreover, we have

$$\begin{array}{c}\hfill |{\mathrm{\Gamma }}_{x_o}|\le \mu .\end{array}$$

Graph

The statement for the degenerate regime is as follows.

Proposition 3.5

Let $\epsilon ,\delta \in (0,1]$ , $\mu >0$ and $\nu \in (0,\frac{1}{4}]$ . Then, there exist constants $\kappa ={\kappa (n,p,\sigma ,\Vert f\Vert }_{n+\sigma },M,\delta ,\nu )\in [2^{-\beta /2},1)$ and $c_o=c_o{(n,p,\Vert f\Vert }_{n+\sigma },M,\delta ,\nu )\ge 1$ such that the there holds: Whenever $B_{2\varrho }(x_o)\subset B_{r_1}$ is a ball with center $x_o\in B_{r_o}$ , and $u_{\epsilon }$ is the unique weak solution of the Dirichlet problem (3.1) and hypothesis (3.4) and (3.5) and the measure condition

3.10 $$\begin{array}{c}\hfill |B_{\varrho }(x_o)\backslash E_{\varrho }^{\nu }(x_o)|\ge \nu |B_{\varrho }(x_o)|\end{array}$$

Graph

are satisfied, then, either

$$\begin{array}{c}\hfill {\mu }^2<c_o{\varrho }^{\beta },\end{array}$$

Graph

or

$$\begin{array}{c}\hfill \underset{B_{\varrho /2}(x_o)}{sup}|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })|\le \kappa \mu \end{array}$$

Graph

hold true.

We postpone the proofs of Proposition 3.4 and 3.5 to Chapters 4–6 and continue with formulating the main result of this subsection.

Theorem 3.6

Let $\epsilon ,\delta \in (0,1]$ and $u_{\epsilon }$ be the unique weak solution of the Dirichlet problem (3.1) in $B_R$ . Then, ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ is Hölder continuous in $B_{r_o}$ for any $0<r_o<R$ with Hölder-exponent ${\alpha }_{\delta }\in (0,1)$ and a Hölder constant $c_{\delta }$ both depending on ${n,N,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M$ and $\delta$ .

Proof

By

$$\begin{array}{c}\hfill \nu ={\nu (n,N,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta )\in (0,\frac{1}{4}]\end{array}$$

Graph

and

$$\begin{array}{c}\hfill \widehat{\varrho }=\widehat{\varrho }{(n,N,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta )\in (0,1]\end{array}$$

Graph

we denote the constants from Proposition 3.4 and by

$$\begin{array}{cc}& \kappa ={\kappa (n,p,\sigma ,\Vert f\Vert }_{n+\sigma },M,\delta ,\nu )\in [2^{-\beta /2},1),\hfill \\ \hfill & c_o=c_o{(n,p,\Vert f\Vert }_{n+\sigma },M,\delta ,\nu )\ge 1\hfill \end{array}$$

Graph

we denote the ones from Proposition 3.5. The dependence of $\nu$ on the structural parameters implies that $\kappa$ depends on ${n,N,p,\sigma ,\Vert f\Vert }_{n+\sigma },M$ and $\delta$ . Finally, we let $\mu =M-1-\delta$ and

$$\begin{array}{c}\hfill {\varrho }_{\ast }=min\{\widehat{\varrho },[\frac{{(\kappa \mu )}^2}{c_o}{]}^{\frac{1}{\beta }}\}.\end{array}$$

Graph

We consider a ball $B_{2\varrho }(x_o)\subset B_{r_1}$ with center $x_o\in B_{r_o}$ and $\varrho \le {\varrho }_{\ast }$ as described above. On this ball we have (3.3)–(3.5) satisfied. Our aim now is to prove that ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })$ is Hölder continuous in $B_{r_o}$ with Hölder exponent

$$\begin{array}{c}\hfill \alpha :=-\frac{log\kappa }{log2}\in (0,\frac{\beta }{2}].\end{array}$$

Graph

In turn, substituting $2\delta$ by $\delta$ , this proves the claim of the proposition. We proceed in two steps.

Step 1: We prove that the limit

$$\begin{array}{c}\hfill {\mathrm{\Gamma }}_{x_o}:=\underset{r\downarrow 0}{lim}({\mathcal{G}}_{2\delta }(D{u}_{\epsilon }){)}_{x_o,r}\end{array}$$

Graph

exists and that

3.11

Graph

holds true for a constant $c=c(n,p,\sigma ,\Vert f{\Vert }_{n+\sigma },M,\delta )$ . To this aim, we define for $i\in {\mathbb{N}}_0$ radii

$$\begin{array}{c}\hfill {\varrho }_i:=2^{-i}\varrho \text{and}{\mu }_i:={\kappa }^i\mu \end{array}$$

Graph

and observe that

3.12 $$\begin{array}{c}\hfill {\mu }_i={\kappa }^i\mu \le 2^{-\alpha i}\mu =(\frac{{\varrho }_i}{\varrho }{)}^{\alpha }\mu \text{for}\text{any}i\in \mathbb{N}.\end{array}$$

Graph

Now, suppose that assumption (3.10) holds on $B_{\varrho }(x_o)$ . Then, Proposition 3.5 yields that either ${\mu }^2<c_o{\varrho }^{\beta }$ or

$$\begin{array}{c}\hfill \underset{B_{{\varrho }_1}(x_o)}{sup}|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })|\le \kappa \mu ={\mu }_1.\end{array}$$

Graph

Note that the first alternative cannot happen, since it would imply

$$\begin{array}{c}\hfill {\mu }^2<c_o{\varrho }^{\beta }\le c_o{\varrho }_{\ast }^{\beta }\le {\kappa }^2{\mu }^2.\end{array}$$

Graph

Hence, we conclude that (3.4) holds on $B_{{\varrho }_1}(x_o)$ with $\mu ={\mu }_1$ . If the measure condition (3.10) is satisfied with $\varrho ={\varrho }_1$ and $\mu ={\mu }_1$ , then a second application of Proposition 3.5 yields that either ${\mu }_1^2<c_o{\varrho }_1^{\beta }$ or

$$\begin{array}{c}\hfill \underset{B_{{\varrho }_2}(x_o)}{sup}|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })|\le \kappa {\mu }_1={\mu }_2.\end{array}$$

Graph

As before, the first alternative cannot happen, since it would imply

$$\begin{array}{c}\hfill {\mu }_1^2<c_o{\varrho }_1^{\beta }\le 2^{-\beta }{(\kappa \mu )}^2\le {\kappa }^2{\mu }_1^2.\end{array}$$

Graph

Assume now that (3.10) is satisfied for $i=1,\cdots ,i_o-1$ up to some $i_o\in \mathbb{N}$ , i.e. that (3.10) holds true on the balls $B_{{\varrho }_i}(x_o)$ with $\mu ={\mu }_i$ . Then, we iteratively conclude that

3.13 $$\begin{array}{c}\hfill \underset{B_{{\varrho }_i}(x_o)}{sup}|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })|\le {\mu }_i,\text{for}i=0,\cdots ,i_o.\end{array}$$

Graph

Now assume that (3.10) fails to hold for some $i_o\in {\mathbb{N}}_0$ . If ${\mu }_{i_o}>\delta$ , the hypothesis of Proposition 3.4 are satisfied on $B_{{\varrho }_{i_o}}(x_o)$ and we conclude that the limit

$$\begin{array}{c}\hfill {\mathrm{\Gamma }}_{x_o}:=\underset{r\downarrow 0}{lim}({\mathcal{G}}_{2\delta }(D{u}_{\epsilon }){)}_{x_o,r}\end{array}$$

Graph

exists and that

Graph

Moreover, we have

3.14 $$\begin{array}{c}\hfill |{\mathrm{\Gamma }}_{x_o}|\le {\mu }_{i_o}.\end{array}$$

Graph

Therefore, we obtain from the preceding inequality and (3.12) that

Graph

holds true for any $0<r\le {\varrho }_{i_o}$ . For a radius $r\in ({\varrho }_{i_o},\varrho ]$ there exists $i\in \{0,\cdots ,i_o\}$ such that ${\varrho }_{i+1}<r\le {\varrho }_i$ . Using (3.13), (3.14) and (3.12) we obtain

Graph

Combining the preceding two inequalities, we have shown (3.11) provided ${\mu }_{i_o}>\delta$ .

In the case ${\mu }_{i_o}\le \delta$ , we have ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })=0$ on $B_{{\varrho }_{i_o}}(x_o)$ . Combining this with (3.13) and keeping in mind that ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })\le {\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ , we obtain

3.15 $$\begin{array}{c}\hfill \underset{B_{{\varrho }_i}(x_o)}{sup}|{\mathcal{G}}_{2\delta }(D{u}_{\epsilon })|\le {\mu }_i,\text{for}\text{any}i\in \{0,1,\cdots \}.\end{array}$$

Graph

In the final case when (3.10) holds for any $i\in \mathbb{N}$ , then (3.13) is satisfied for any $i\in \mathbb{N}$ and hence we obtain (3.15) also in this case. (3.15), however, implies

$$\begin{array}{c}\hfill {\mathrm{\Gamma }}_{x_o}:=\underset{r\downarrow 0}{lim}({\mathcal{G}}_{2\delta }(D{u}_{\epsilon }){)}_{x_o,r}=0.\end{array}$$

Graph

For $r\in (0,\varrho ]$ we find $i\in \{0,1,\cdots \}$ such that ${\varrho }_{i+1}<r\le {\varrho }_i$ . Then, (3.15) and (3.12) imply

Graph

This establishes (3.11) in the remaining cases.

Step 2: Here, we prove that the Lebesgue representative $x\mapsto {\mathrm{\Gamma }}_x$ of ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })$ is Hölder continuous in $B_{r_o}$ . The proof is standard once the excess decay (3.11) is established. For convenience of the reader we give the details. We consider $x_1,x_2\in B_{r_o}$ . If $r:=|x_1-x_2|\le {\varrho }_{\ast }$ we define $\tilde{x}:=\frac{1}{2}(x_1+x_2)$ and obtain from Step 1 that

Graph

This can be re-written in the form

$$\begin{array}{c}\hfill |{\mathrm{\Gamma }}_{x_1}-{\mathrm{\Gamma }}_{x_2}|\le c(\frac{|x_1-x_2|}{{\varrho }_{\ast }}{)}^{\alpha }\mu .\end{array}$$

Graph

Otherwise, if $r=|x_1-x_2|>{\varrho }_{\ast }$ , we trivially have

$$\begin{array}{c}\hfill |{\mathrm{\Gamma }}_{x_1}-{\mathrm{\Gamma }}_{x_2}|\le 2\mu \le 2(\frac{|x_1-x_2|}{{\varrho }_{\ast }}{)}^{\alpha }\mu .\end{array}$$

Graph

Together, this establishes that the Lebesgue representative $x\mapsto {\mathrm{\Gamma }}_x$ of ${\mathcal{G}}_{2\delta }(D{u}_{\epsilon })$ is Hölder continuous in $B_{r_o}$ with Hölder exponent $\alpha$ . Note that $\alpha$ admits the same dependencies as $\kappa$ , i.e. $\alpha =\alpha (n,N,p,\sigma ,\Vert f{\Vert }_{n+\sigma },M,\delta )$ . This finishes the proof of the theorem. $\square$

Continuity of G(Du)

In this subsection it is important that all estimates are independent of $\epsilon$ . More precisely, constants might depend on $\delta$ , but are independent of $\epsilon$ .

Proof of Theorem 1.1

We let $\epsilon \in (0,1]$ and consider a fixed ball $B_R\equiv B_R(y_o)⋐\mathrm{\Omega }$ . By $u_{\epsilon }$ we denote the weak solution to (3.1) constructed in Sect. 3.1. Next, we fix $\delta \in (0,1]$ and $r\in (0,R)$ . From Theorem 3.6 we know that ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ is Hölder continuous in ${\overline{B}}_r$ with Hölder-exponent ${\alpha }_{\delta }\in (0,1)$ and constant $c_{\delta }$ , both depending at most on ${n,N,p,\sigma ,\Vert f\Vert }_{n+\sigma },M$ , and $\delta$ . From Lemma 3.3 we know that ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })\to {\mathcal{G}}_{\delta }(Du)$ in $L^p(B_R,{\mathbb{R}}^{\mathit{Nn}})$ as $\epsilon \downarrow 0$ . This implies that there exists a subsequence ${\epsilon }_i\downarrow 0$ as $i\to \infty$ such that ${\mathcal{G}}_{\delta }(D{u}_{{\epsilon }_i})\to {\mathcal{G}}_{\delta }(Du)$ a.e. in $B_R$ . On the other hand, by Ascoli-Arzelà's Theorem we conclude that ${\mathcal{G}}_{\delta }(D{u}_{{\epsilon }_i})$ converges uniformly on compact subsets of $B_R$ . Therefore the limit ${\mathcal{G}}_{\delta }(Du)$ is Hölder continuous in ${\overline{B}}_r$ with Hölder-exponent ${\alpha }_{\delta }\in (0,1)$ and constant $c_{\delta }$ . In particular, ${\mathcal{G}}_{\delta }(Du)$ is continuous in ${\overline{B}}_r$ for any $\delta \in (0,1]$ . Moreover, we have ${\mathcal{G}}_{\delta }(Du)\to \mathcal{G}(Du)$ uniformly in ${\overline{B}}_r$ as $\delta \downarrow 0$ . Indeed

$$\begin{array}{cc}\hfill |{\mathcal{G}}_{\delta }(Du)-\mathcal{G}(Du)|& =|{(|Du|-1-\delta )}_{+}\frac{\mathit{Du}}{|Du|}-{(|Du|-1)}_{+}\frac{\mathit{Du}}{|Du|}|\hfill \\ \hfill & =|{(|Du|-1-\delta )}_{+}-{(|Du|-1)}_{+}|\le \delta \hfill \end{array}$$

Graph

in ${\overline{B}}_r$ . As the uniform limit of a sequence of continuous functions, $\mathcal{G}(Du)$ itself is continuous on ${\overline{B}}_r$ . Observe that $\mathcal{G}(Du)$ is also uniformly continuous on ${\overline{B}}_r$ .

Now, let $\mathcal{K}:{\mathbb{R}}^{\mathit{Nn}}\to \mathbb{R}$ be any continuous function vanishing on $\{\xi \in {\mathbb{R}}^{\mathit{Nn}}:|\xi |\le 1\}$ . Since $u\in W_{\mathrm{loc}}^{1,\infty }(\mathrm{\Omega },{\mathbb{R}}^N)$ , we find $M>0$ such that $|Du|\le M$ on ${\overline{B}}_r$ . By $\omega :{\mathbb{R}}_{+}\to {\mathbb{R}}_{+}$ we denote the modulus of continuity of $\mathcal{K}$ on $\{\xi \in {\mathbb{R}}^{\mathit{Nn}}:|\xi |\le M\}$ , i.e. for any $\xi ,\zeta \in {\mathbb{R}}^{\mathit{Nn}}$ with $|\xi |,|\zeta |\le M$ we have $|\mathcal{K}(\xi )-\mathcal{K}(\zeta )|\le \omega (|\xi -\zeta |)$ . Next, given $\epsilon \in (0,1)$ we choose $\delta >0$ such that

3.16 $$\begin{array}{c}\hfill |\mathcal{G}(Du(x))-\mathcal{G}(Du(y))|<\epsilon \text{for}\text{any}x,y\in {\overline{B}}_r\text{with}|x-y|<\delta .\end{array}$$

Graph

We now distinguish two cases. First, we assume $|Du(x)|\le 1+\sqrt{\epsilon }$ . If $|Du(x)|\ge 1$ we use $\mathcal{K}=0$ on $\{|\xi |\le 1\}$ to conclude

$$\begin{array}{c}\hfill |\mathcal{K}(Du(x))|=|\mathcal{K}(Du(x))-\mathcal{K}(\frac{Du(x)}{|Du(x)|})|\le \omega (|Du(x)-\frac{Du(x)}{|Du(x)|}|)\le \omega (\sqrt{\epsilon }).\end{array}$$

Graph

The preceding estimate trivially holds if $|Du(x)|\le 1$ . Moreover, by (3.16) we have

$$\begin{array}{cc}\hfill (|Du(y)|-1{)}_{+}& =|\mathcal{G}(Du(y))|\le |\mathcal{G}(Du(y))-\mathcal{G}(Du(x))|+|\mathcal{G}(Du(x))|\hfill \\ \hfill & \le \epsilon +(|Du(x)|-1{)}_{+}\le \epsilon +\sqrt{\epsilon }\le 2\sqrt{\epsilon }.\hfill \end{array}$$

Graph

This implies $|Du(y)|\le 1+2\sqrt{\epsilon }$ . Similarly as above we conclude

$$\begin{array}{c}\hfill |\mathcal{K}(Du(y))|\le \omega (2\sqrt{\epsilon }).\end{array}$$

Graph

Combining the estimates from above, we end up with

$$\begin{array}{c}\hfill |\mathcal{K}(Du(x))-\mathcal{K}(Du(y))|\le \omega (\sqrt{\epsilon })+\omega (2\sqrt{\epsilon })\le 2\omega (2\sqrt{\epsilon }).\end{array}$$

Graph

Now, we consider the case $|Du(x)|>1+\sqrt{\epsilon }$ . Here, Lemma 2.3 and (3.16) imply

$$\begin{array}{c}\hfill |Du(x)-Du(y)|\le (1+\frac{2}{\sqrt{\epsilon }})|\mathcal{G}(Du(x))-\mathcal{G}(Du(y))|\le \epsilon +2\sqrt{\epsilon }\le 3\sqrt{\epsilon },\end{array}$$

Graph

proving

$$\begin{array}{c}\hfill |\mathcal{K}(Du(x))-\mathcal{K}(Du(y))|\le \omega (3\sqrt{\epsilon }).\end{array}$$

Graph

Hence, $\mathcal{K}(Du)$ is continuous on ${\overline{B}}_r$ . Since $B_r⋐B_R⋐\mathrm{\Omega }$ were arbitrary, we have shown that $\mathcal{K}(Du)$ is continuous in $\mathrm{\Omega }$ . This completes the proof of the theorem. $\square$

As mentioned above, we have now finished the proof of the main theorem Theorem 1.1 under the condition that Propositions 3.4 and 3.5 are true. The rest of the paper is now devoted to the proof of those two propositions.

Conclusions from the differentiated system

The main integral inequality for second derivatives

Throughout this subsection we assume as a general requirement that $u_{\epsilon }:B_R\to {\mathbb{R}}^N$ is a weak solution to the regularized system (3.1). Instead of $u_{\epsilon }$ , we write u for the sake of simplicity. In contrast, we will continue to use the subscript $\epsilon$ in the notation for the coefficients $h_{\epsilon }$ and its associated bilinear forms, such as ${\mathcal{C}}_{\epsilon }$ . We recall that the bilinear forms have been defined in Sect. 2.3.

For some index $\beta =1,\cdots ,n$ we differentiate the regularized system (3.1) with respect to $x_{\beta }$ and obtain

4.1 $$\begin{array}{c}\hfill {\int }_{B_R}{D}_{\beta }[h_{\epsilon }(|Du|)Du]·D\phi \mathrm{d}x={\int }_{B_R}f·D_{\beta }\phi \mathrm{d}x,\end{array}$$

Graph

for any $\phi \in W_0^{1,p}(B_R,{\mathbb{R}}^N)$ . We have

$$\begin{array}{cc}\hfill D_{\beta }[h_{\epsilon }(|Du|)D_{\alpha }{u}^i]& =h_{\epsilon }(|Du|)D_{\alpha }{D}_{\beta }{u}^i+h_{\epsilon }^{\prime }(|Du|)|Du|\frac{D_{\alpha }{u}^i{D}_{\gamma }{u}^j{u}_{x_{\gamma }{x}_{\beta }}^j}{{|Du|}^2}\hfill \\ \hfill & =[h_{\epsilon }(|Du|){\delta }^{\mathit{ij}}{\delta }_{\gamma \alpha }+h_{\epsilon }^{\prime }(|Du|)|Du|\frac{D_{\alpha }{u}^i{D}_{\gamma }{u}^j}{{|Du|}^2}]u_{x_{\gamma }{x}_{\beta }}^j.\hfill \end{array}$$

Graph

In (4.1) we choose the testing function $\phi =\zeta \varphi (|Du|)D_{\beta }u$ , where $\zeta \in C_0^1(B_R)$ is non-negative and $\varphi \in W_{\text{loc}}^{1,\infty }({\mathbb{R}}_{\ge 0},{\mathbb{R}}_{\ge 0})$ is non-decreasing. Note that

$$\begin{array}{cc}\hfill D_{\alpha }\phi & =\zeta \varphi (|Du|)D_{\alpha }{D}_{\beta }u+\zeta {\varphi }^{\prime }(|Du|)D_{\alpha }|Du|D_{\beta }u+D_{\alpha }\zeta \varphi (|Du|)D_{\beta }u.\hfill \end{array}$$

Graph

The resulting equations are then summed with respect to $\beta$ from 1 to n. This leads to

$$\begin{array}{cc}& {\int }_{B_R}[h_{\epsilon }(|Du|)|D^2{u|}^2+h_{\epsilon }^{\prime }(|Du|)|Du||\mathrm{\nabla }|Du|{|}^2]\varphi (|Du|)\zeta \mathrm{d}x\hfill \\ \hfill & +{\int }_{B_R}[h_{\epsilon }(|Du|)|Du||\mathrm{\nabla }|Du|{|}^2+h_{\epsilon }^{\prime }(|Du|)|Du\mathrm{\nabla }|Du|{|}^2]{\varphi }^{\prime }(|Du|)\zeta \mathrm{d}x\hfill \\ \hfill & +{\int }_{B_R}[\underset{={\mathcal{C}}_{\epsilon }(Du)(\mathrm{\nabla }|Du|,\mathrm{\nabla }\zeta )|Du|}{\underbrace{h_{\epsilon }(|Du|)|Du|\mathrm{\nabla }|Du|·\mathrm{\nabla }\zeta +h_{\epsilon }^{\prime }(|Du|)Du\mathrm{\nabla }|Du|·Du\mathrm{\nabla }\zeta }}]\varphi (|Du|)\mathrm{d}x\hfill \\ \hfill & ={\int }_{B_R}f·D_{\beta }[\zeta \varphi (|Du|)D_{\beta }u]\mathrm{d}x.\hfill \end{array}$$

Graph

Now, we compute the right-hand side.

$$\begin{array}{cc}& {\int }_{B_R}f·D_{\beta }[\zeta \varphi (|Du|)D_{\beta }u]\mathrm{d}x\hfill \\ \hfill & ={\int }_{\mathrm{\Omega }}f·D_{\beta }{D}_{\beta }u\varphi (|Du|)\zeta \mathrm{d}x+{\int }_{\mathrm{\Omega }}f·{\varphi }^{\prime }(|Du|)Du\mathrm{\nabla }|Du|\zeta \mathrm{d}x\hfill \\ \hfill & +{\int }_{\mathrm{\Omega }}f·\varphi (|Du|)Du\mathrm{\nabla }\zeta \mathrm{d}x\hfill \\ \hfill & ={\mathbf{R}}_1+{\mathbf{R}}_2+{\mathbf{R}}_3,\hfill \end{array}$$

Graph

with the obvious meaning of ${\mathbf{R}}_i$ . For the first term, we have by Young's inequality

$$\begin{array}{c}\hfill {\mathbf{R}}_1\le \tau {\int }_{B_R}{h}_{\epsilon }(|Du|)|D^2{u|}^2\varphi (|Du|)\zeta \mathrm{d}x+\frac{1}{\tau }{\int }_{B_R}\frac{{|f|}^2\varphi (|Du|)}{h_{\epsilon }(|Du|)}\zeta \mathrm{d}x\end{array}$$

Graph

for any $\tau \in (0,1)$ . Similarly, we get

$$\begin{array}{c}\hfill {\mathbf{R}}_2\le \tau {\int }_{B_R}{h}_{\epsilon }(|Du|)|Du||\mathrm{\nabla }|Du|{|}^2{\varphi }^{\prime }(|Du|)\zeta \mathrm{d}x+\frac{1}{\tau }{\int }_{B_R}\frac{{|f|}^2{\varphi }^{\prime }(|Du|)|Du|}{h_{\epsilon }(|Du|)}\zeta \mathrm{d}x\end{array}$$

Graph

and

$$\begin{array}{c}\hfill {\mathbf{R}}_3\le {\int }_{\mathrm{\Omega }}|f|\varphi (|Du|)|Du||\mathrm{\nabla }\zeta |\mathrm{d}x.\end{array}$$

Graph

Inserting this above and re-absorbing the terms containing second derivatives from the right-hand side into the left, we obtain

$$\begin{array}{cc}& {\int }_{B_R}[(1-\tau )h_{\epsilon }(|Du|)|D^2{u|}^2+h_{\epsilon }^{\prime }(|Du|)|Du||\mathrm{\nabla }|Du|{|}^2]\varphi (|Du|)\zeta \mathrm{d}x\hfill \\ \hfill & +{\int }_{B_R}[(1-\tau )h_{\epsilon }(|Du|)|\mathrm{\nabla }|Du|{|}^2+h_{\epsilon }^{\prime }(|Du|)|Du|\frac{|Du\mathrm{\nabla }|Du|{|}^2}{{|Du|}^2}]{\varphi }^{\prime }(|Du|)|Du|\zeta \mathrm{d}x\hfill \\ \hfill & +{\int }_{B_R}{\mathcal{C}}_{\epsilon }(Du)(\mathrm{\nabla }|Du|,\mathrm{\nabla }\zeta )\varphi (|Du|)|Du|\mathrm{d}x\hfill \\ \hfill & \le \frac{1}{\tau }{\int }_{B_R}{|f|}^2[\frac{\varphi (|Du|)}{h_{\epsilon }(|Du|)}+\frac{{\varphi }^{\prime }(|Du|)|Du|}{h_{\epsilon }(|Du|)}]\zeta \mathrm{d}x+{\int }_{B_R}|f|\varphi (|Du|)|Du||D\zeta |\mathrm{d}x\hfill \end{array}$$

Graph

for any non-negative function $\zeta \in C_0^1(B_R)$ . In the preceding inequality the parameter $\tau$ is at our disposal. We choose $\tau =\frac{1}{2}min\{1,p-1\}$ . For the first factor, i.e. the term $[\cdots ]$ , in the integrand of the first integral on the left-hand side we have

$$\begin{array}{cc}& (1-\tau )h_{\epsilon }(|Du|)|D^2{u|}^2+h_{\epsilon }^{\prime }(|Du|)|Du||\mathrm{\nabla }|Du|{|}^2-\frac{1}{2}{\mathcal{A}}_{\epsilon }(Du)(D^{2}u,D^{2}u)\hfill \\ \hfill & =\frac{1}{2}max\{0,2-p\}{h}_{\epsilon }(|Du|)|D^2{u|}^2+\frac{1}{2}{h}_{\epsilon }^{\prime }(|Du|)|Du||D|Du|{|}^2\ge 0.\hfill \end{array}$$

Graph

In fact, if $h_{\epsilon }^{\prime }(|Du|)\ge 0$ the inequality is obvious. If otherwise $h_{\epsilon }^{\prime }(|Du|)<0$ , which can only happen if $p<2$ and $|Du|>1$ , the result follows by an application of Kato's inequality and (2.4). Indeed

$$\begin{array}{cc}& \frac{1}{2}max\{0,2-p\}{h}_{\epsilon }(|Du|)|D^2{u|}^2+\frac{1}{2}{h}_{\epsilon }^{\prime }(|Du|)|Du||\mathrm{\nabla }|Du|{|}^2\hfill \\ \hfill & \ge \frac{1}{2}[(2-p)h(|Du|)+h^{\prime }(|Du|)|Du|]|D^2{u|}^2\hfill \\ \hfill & =\frac{1}{2}[(2-p)\frac{{(|Du|-1)}_{+}^{p-1}}{|Du|}+\frac{{(|Du|-1)}_{+}^{p-2}[(p-2)|Du|+1]}{|Du|}]{|D^{2}u|}^2\hfill \\ \hfill & =\frac{{(|Du|-1)}_{+}^{p-2}}{2|Du|}[(2-p)(|Du|-1)+(p-2)|Du|+1]|D^2{u|}^2\hfill \\ \hfill & =(p-1)\frac{{(|Du|-1)}_{+}^{p-2}}{2|Du|}{|D^{2}u|}^2\ge 0.\hfill \end{array}$$

Graph

For the term in brackets of the second integral on the left-hand side a similar computation applies. The result of the calculation is

$$\begin{array}{cc}& (1-\tau )h_{\epsilon }(|Du|)|D|Du|{|}^2+h_{\epsilon }^{\prime }(|Du|)|Du|\frac{|Du\mathrm{\nabla }|Du|{|}^2}{{|Du|}^2}\hfill \\ \hfill & -\frac{1}{2}{\mathcal{C}}_{\epsilon }(Du)(\mathrm{\nabla }|Du|,\mathrm{\nabla }|Du|)\ge 0.\hfill \end{array}$$

Graph

Using the last and second last inequality above, we obtain an inequality which can be interpreted in two ways. On the one hand it can be seen as an energy inequality for the second derivatives of u. On the other hand—by discarding on the left-hand side the two non-negative terms containing second derivatives—the inequality implies that |Du| is a subsolution of an elliptic equation with measurable coefficients.

Lemma 4.1

Let $\epsilon \in (0,1]$ and $u=u_{\epsilon }$ a weak solution to the regularized system (3.1) on $B_R$ . Then, for any non-decreasing function $\varphi \in W_{\text{loc}}^{1,\infty }({\mathbb{R}}_{\ge 0},{\mathbb{R}}_{\ge 0})$ and any non-negative testing function $\zeta \in C_0^1(B_R)$ we have

$$\begin{array}{cc}& {\int }_{B_R}[{\mathcal{A}}_{\epsilon }(Du)(D^{2}u,D^{2}u)\varphi (|Du|)+{\mathcal{C}}_{\epsilon }(Du)(\mathrm{\nabla }|Du|,\mathrm{\nabla }|Du|){\varphi }^{\prime }(|Du|)|Du|]\zeta \mathrm{d}x\hfill \\ \hfill & +2{\int }_{B_R}{\mathcal{C}}_{\epsilon }(Du)(\mathrm{\nabla }|Du|,\mathrm{\nabla }\zeta )\varphi (|Du|)|Du|\mathrm{d}x\hfill \\ \hfill & \le c{\int }_{B_R}{|f|}^2[\frac{\varphi (|Du|)}{h_{\epsilon }(|Du|)}+\frac{{\varphi }^{\prime }(|Du|)|Du|}{h_{\epsilon }(|Du|)}]\zeta \mathrm{d}x+2{\int }_{B_R}|f|\varphi (|Du|)|Du||\mathrm{\nabla }\zeta |\mathrm{d}x,\hfill \end{array}$$

Graph

where $c:=\frac{4}{min\{1,p-1\}}$ . $\square$

Subsolution to an elliptic equation

We start by showing that $(|D{u}_{\epsilon }{|-1-\delta )}_{+}^2$ is a sub-solution of a certain elliptic equation. More precisely

Lemma 4.2

Let $\epsilon \in (0,1]$ and $u_{\epsilon }\in W^{1,p}(B_R,{\mathbb{R}}^N)$ be a weak solution of the regularized system (3.1) satisfying (3.4) and (3.5) on $B_{\varrho }(x_o)\subset B_{r_o}⋐B_R$ . Then, the function

4.2 $$\begin{array}{c}\hfill U_{\epsilon }:=(|D{u}_{\epsilon }|-1-\delta {)}_{+}^2\end{array}$$

Graph

is a sub-solution of a linear elliptic equation on $B_{\varrho }(x_o)$ , in the sense that

4.3 $$\begin{array}{c}\hfill {\int }_{B_{\varrho }(x_o)}{A}_{\alpha \beta }{D}_{\alpha }{U}_{\epsilon }{D}_{\beta }\zeta \mathrm{d}x\le c{\int }_{B_{\varrho }(x_o)}{|f|}^2\zeta \mathrm{d}x+c\mu {\int }_{B_{\varrho }(x_o)}|f||D\zeta |\mathrm{d}x,\end{array}$$

Graph

holds true for any non-negative test function $\zeta \in C_0^1(B_{\varrho }(x_o))$ and with a universal constant $c=c(p,M,\delta )$ . The coefficients $A_{\alpha \beta }$ are given by

$$\begin{array}{c}\hfill A_{\alpha \beta }{\eta }_{\alpha }{\zeta }_{\beta }=|D{u}_{\epsilon }|{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\eta ,\zeta )\text{for}\eta ,\zeta \in {\mathbb{R}}^n,\end{array}$$

Graph

where ${\mathcal{C}}_{\epsilon }$ is the bilinear form defined in (2.8).

Proof

We apply Lemma 4.1 with $\varphi (t)={(t-1-\delta )}_{+}$ . Due to Lemma 2.7 the first two integrals on the left-hand side are non-negative and therefore can be discarded. In this way, we obtain with $c=\frac{4}{min\{1,p-1\}}$ that

4.4 $$\begin{array}{cc}\hfill \mathbf{L}& :={\int }_{B_{\varrho }(x_o)}|D{u}_{\epsilon }|{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }|D{u}_{\epsilon }|,\mathrm{\nabla }\zeta )\varphi (|D{u}_{\epsilon }|)\mathrm{d}x\hfill \\ \hfill & \le c{\int }_{B_{\varrho }(x_o)}{|f|}^2[\frac{\varphi (|D{u}_{\epsilon }|)}{h_{\epsilon }(|D{u}_{\epsilon }|)}+\frac{{\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{h_{\epsilon }(|D{u}_{\epsilon }|)}]\zeta \mathrm{d}x\hfill \\ \hfill & +{\int }_{B_{\varrho }(x_o)}|f|\varphi (|D{u}_{\epsilon }|)|D{u}_{\epsilon }||D\zeta |\mathrm{d}x=:\mathbf{R}\hfill \end{array}$$

Graph

holds true for any non-negative $\zeta \in C_0^1(B_{\varrho }(x_o))$ . To proceed further we compute

$$\begin{array}{c}\hfill \frac{\varphi (t)}{h_{\epsilon }(t)}\le \frac{\varphi (t)}{h(t)}=\frac{t{(t-1-\delta )}_{+}}{{(t-1)}^{p-1}}\le \frac{t^2}{{\delta }^{p-1}}\end{array}$$

Graph

and

$$\begin{array}{c}\hfill \frac{{\varphi }^{\prime }(t)t}{h_{\epsilon }(t)}\le \frac{{\varphi }^{\prime }(t)t}{h(t)}=\frac{{\mathit{\chi }}_{\{t>1+\delta \}}{t}^2}{{(t-1)}^{p-1}}\le \frac{t^2}{{\delta }^{p-1}}.\end{array}$$

Graph

The above inequalities together with the general assumptions (3.4), (3.5) allow to estimate

$$\begin{array}{c}\hfill \frac{\varphi (|D{u}_{\epsilon }|)}{h_{\epsilon }(|D{u}_{\epsilon }|)}+\frac{{\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{h_{\epsilon }(|D{u}_{\epsilon }|)}\le \frac{2M^2}{{\delta }^{p-1}}.\end{array}$$

Graph

Moreover, we have $\varphi (|D{u}_{\epsilon }|)|D{u}_{\epsilon }|\le M\mu$ . In this way, we obtain for the right-hand side in (4.4) the estimate

$$\begin{array}{cc}\hfill \mathbf{R}& \le c{\int }_{B_R}{|f|}^2\zeta \mathrm{d}x+c\mu {\int }_{B_R}|f||D\zeta |\mathrm{d}x,\hfill \end{array}$$

Graph

where $c=c(p,M,\delta )$ . Now, we consider the left-hand side in (4.4). Observe that $\mathrm{\nabla }{U}_{\epsilon }=2\varphi (|D{u}_{\epsilon }|)\mathrm{\nabla }|D{u}_{\epsilon }|$ , so that by the linearity of ${\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })$ with respect to the first variable we have

$$\begin{array}{cc}\hfill {\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }|D{u}_{\epsilon }|,\mathrm{\nabla }\zeta )\varphi (|D{u}_{\epsilon }|)& =\frac{1}{2}{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }{U}_{\epsilon },\mathrm{\nabla }\zeta ).\hfill \end{array}$$

Graph

For the left-hand side this has the consequence that

$$\begin{array}{cc}\hfill \mathbf{L}& =\frac{1}{2}{\int }_{B_{\varrho }(x_o)}|D{u}_{\epsilon }|{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }{U}_{\epsilon },\mathrm{\nabla }\zeta )\mathrm{d}x=\frac{1}{2}{\int }_{B_{\varrho }(x_o)}{A}_{\alpha \beta }(x)D_{\alpha }{U}_{\epsilon }{D}_{\beta }\zeta \mathrm{d}x\hfill \end{array}$$

Graph

holds true. Here, we have taken into account the definition of the coefficients $A_{\alpha \beta }$ . Altogether we have shown the claim (4.3). $\square$

The coefficients $A_{\alpha \beta }$ in (4.3) can be explicitly written as

$$\begin{array}{c}\hfill A_{\alpha \beta }:=h_{\epsilon }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|[{\delta }_{\alpha \beta }+\frac{h_{\epsilon }^{\prime }(|D{u}_{\epsilon }|)D_{\alpha }{u}_{\epsilon }·D_{\beta }{u}_{\epsilon }}{h_{\epsilon }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}].\end{array}$$

Graph

They are only degenerate elliptic due to the factor $h_{\epsilon }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|$ which, for $\epsilon =0$ , vanishes on the set $\{|Du|\le 1\}$ . On the other hand $U_{\epsilon }$ has its support in the set $B_R\cap \{|D{u}_{\epsilon }|\ge 1+\delta \}$ . This allows us to modify the coefficients on $B_R\cap \{|D{u}_{\epsilon }|\le 1+\delta \}$ . This idea will lead us to an energy estimate for $U_{\epsilon }$ in the next lemma.

Lemma 4.3

Let $\epsilon \in (0,1]$ and $u_{\epsilon }\in W^{1,p}(B_R,{\mathbb{R}}^N)$ be a weak solution of the regularized system (3.1) satisfying (3.4) and (3.5) on $B_{\varrho }(x_o)\subset B_{r_o}⋐B_R$ and denote by $U_{\epsilon }$ the function defined in (4.2). Then, for any $k>0$ and any $\tau \in (0,1)$ we have

$$\begin{array}{cc}& {\int }_{B_{\tau \varrho }(x_o)}{|D{(U_{\epsilon }-k)}_{+}|}^2\mathrm{d}x\hfill \\ \hfill & \le \frac{c}{{(1-\tau )}^2{\varrho }^2}{\int }_{B_{\varrho }(x_o)}{(U_{\epsilon }-k)}_{+}^2\mathrm{d}x+c{\Vert f\Vert }_{n+\sigma }^2|B_{\varrho }(x_o)\cap \{U_{\epsilon }>k\}{|}^{1-\frac{2}{n}+\frac{2\beta }{n}},\hfill \end{array}$$

Graph

where $c=c(n,p,M,\delta )$ and $\beta =\frac{\sigma }{n+\sigma }$ .

Proof

We let $A_{\alpha \beta }$ be the coefficients defined in Lemma 4.2 and extend them from $B_R\cap \{|D{u}_{\epsilon }|>1+\delta \}$ to the complement $B_R\cap \{|D{u}_{\epsilon }|\le 1+\delta \}$ by letting $A_{\alpha \beta }\equiv {\delta }_{\alpha \beta }$ . The new coefficients ${\tilde{A}}_{\alpha \beta }$ are thus defined by

$$\begin{array}{c}\hfill {\tilde{A}}_{\alpha \beta }(x):=\left\{\begin{array}{cc}{\delta }_{\alpha \beta },& \text{on}\{x\in B_R:|D{u}_{\epsilon }(x)|\le 1+\delta \},\hfill \\ [7pt]A_{\alpha \beta }(x),& \text{on}\{x\in B_R:|D{u}_{\epsilon }(x)|>1+\delta \}.\hfill \end{array}\right)\end{array}$$

Graph

From this definition and Lemma 4.2 we observe that $U_{\epsilon }$ is a weak sub-solution also with the modified coefficients. More precisely, we have that

4.5 $$\begin{array}{c}\hfill {\int }_{B_{\varrho }(x_o)}{\tilde{A}}_{\alpha \beta }{D}_{\alpha }{U}_{\epsilon }{D}_{\beta }\zeta \mathrm{d}x\le c{\int }_{B_{\varrho }(x_o)}{|f|}^2\zeta \mathrm{d}x+c\mu {\int }_{B_{\varrho }(x_o)}|f||D\zeta |\mathrm{d}x,\end{array}$$

Graph

for any non-negative $\zeta \in C_0^1(B_R)$ .

We now investigate the upper bound and ellipticity of the coefficients ${\tilde{A}}_{\alpha \beta }$ . We will show that there exist $0<\lambda \le \mathrm{\Lambda }<\infty$ both depending at most on p, M and $\delta$ such that

$$\begin{array}{c}\hfill {\lambda |\zeta |}^2\le {\tilde{A}}_{\alpha \beta }(x){\zeta }_{\alpha }{\zeta }_{\beta }\le \mathrm{\Lambda }{|\zeta |}^2\end{array}$$

Graph

for any $x\in B_{r_o}$ and $\zeta \in {\mathbb{R}}^n$ . We start with the former one. On the set where $|D{u}_{\epsilon }|\le 1+\delta$ the upper bound holds with $\mathrm{\Lambda }=1$ , while on the set where $|D{u}_{\epsilon }|>1+\delta$ we have from Lemma 2.7 that

$$\begin{array}{cc}\hfill {\tilde{A}}_{\alpha \beta }{\zeta }_{\alpha }{\zeta }_{\beta }& =|D{u}_{\epsilon }|{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\zeta ,\zeta )\hfill \\ \hfill & \le |D{u}_{\epsilon }|[\epsilon +\mathbf{\Lambda }(|D{u}_{\epsilon }{|)]|\zeta |}^2\hfill \\ \hfill & =[\epsilon |D{u}_{\epsilon }|+max\{(|D{u}_{\epsilon }{|-1)}^{p-1},(p-1)|D{u}_{\epsilon }|(|D{u}_{\epsilon }{|-1)}^{p-2}\}]{|\zeta |}^2\hfill \\ \hfill & \le [M+\frac{p{M}^p}{\delta }]{|\zeta |}^2,\hfill \end{array}$$

Graph

which proves the claim with $\mathrm{\Lambda }=\mathrm{\Lambda }(p,M,\delta )=M+\frac{p{M}^p}{\delta }$ . Similarly, on the set where $|D{u}_{\epsilon }|\le 1+\delta$ the ellipticity holds with $\lambda =1$ , while on the set where $|D{u}_{\epsilon }|>1+\delta$ we have from Lemma 2.7 that

$$\begin{array}{cc}\hfill {\tilde{A}}_{\alpha \beta }{\zeta }_{\alpha }{\zeta }_{\beta }& =|D{u}_{\epsilon }|{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\zeta ,\zeta )\hfill \\ \hfill & \ge |D{u}_{\epsilon }|[\epsilon +\mathit{\lambda }(|D{u}_{\epsilon }{|)]|\zeta |}^2\hfill \\ \hfill & \ge min\{(|D{u}_{\epsilon }{|-1)}^{p-1},(p-1)|D{u}_{\epsilon }|(|D{u}_{\epsilon }{|-1)}^{p-2}\}{|\zeta |}^2\hfill \\ \hfill & \ge min\{1,p-1\}{\delta }^{p-1}{|\zeta |}^2,\hfill \end{array}$$

Graph

which proves the claim with $\lambda =\lambda (p,\delta )=min\{1,p-1\}{\delta }^{p-1}$ .

Now, the claimed energy estimate follows in a standard way by choosing in (4.5) a test-function of the form $\zeta ={\eta }^2{(U_{\epsilon }-k)}_{+}$ with a cut-off function $\eta \in C_0^1(B_{\varrho }(x_o))$ with $\eta \equiv 1$ on $B_{\tau \varrho }(x_o)$ and $|\mathrm{\nabla }\eta |\le \frac{2}{\tau \varrho }$ , cf. [[19], Chapter 10.1]. $\square$

Energy estimates

Here, we assume that the hypothesis of Proposition 3.4 are in force. Our starting point is again Lemma 4.1. This time we keep the two non-negative terms containing the quadratic forms ${\mathcal{A}}_{\epsilon }$ and ${\mathcal{C}}_{\epsilon }$ on the left-hand side. For any non-decreasing function $\varphi \in W_{\text{loc}}^{1,\infty }({\mathbb{R}}_{\ge 0},{\mathbb{R}}_{\ge 0})$ and any non-negative function $\zeta ={\eta }^2\in C_0^1(B_{\varrho }(x_o))$ we have

$$\begin{array}{cc}& {\int }_{B_{\varrho }(x_o)}[{\mathcal{A}}_{\epsilon }(D{u}_{\epsilon })(D^2{u}_{\epsilon },D^2{u}_{\epsilon })\varphi (|D{u}_{\epsilon }|)\hfill \\ \hfill & +{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }|D{u}_{\epsilon }|,\mathrm{\nabla }|D{u}_{\epsilon }|){\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|]{\eta }^2\mathrm{d}x\hfill \\ \hfill & \le 4|{\int }_{B_{\varrho }(x_o)}{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }|D{u}_{\epsilon }|,\mathrm{\nabla }\eta )\varphi (|D{u}_{\epsilon }|)|D{u}_{\epsilon }|\eta \mathrm{d}x|\hfill \\ \hfill & +c(p){\int }_{B_{\varrho }(x_o)}{|f|}^2[\frac{\varphi (|D{u}_{\epsilon }|)}{h_{\epsilon }(|D{u}_{\epsilon }|)}+\frac{{\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{h_{\epsilon }(|D{u}_{\epsilon }|)}]{\eta }^2\mathrm{d}x\hfill \\ \hfill & +4{\int }_{B_{\varrho }(x_o)}|f|\varphi (|D{u}_{\epsilon }|)|D{u}_{\epsilon }||\mathrm{\nabla }\eta |\eta \mathrm{d}x\hfill \\ \hfill & =:4\mathbf{I}+c(p)\mathbf{II}+4\mathbf{III}\hfill \end{array}$$

Graph

with the obvious meaning of $\mathbf{I}$ – $\mathbf{III}$ . For $\mathbf{III}$ we have

$$\begin{array}{cc}\hfill 4\mathbf{III}& \le 2{\int }_{B_{\varrho }(x_o)}{[|f|}^2\frac{{\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{h_{\epsilon }(|D{u}_{\epsilon }|)}{\eta }^2+h_{\epsilon }(|D{u}_{\epsilon }|)\frac{{\varphi }^2(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{{\varphi }^{\prime }(|D{u}_{\epsilon }|)}{|\mathrm{\nabla }\eta |}^2]\mathrm{d}x.\hfill \end{array}$$

Graph

Moreover, we estimate the integral $\mathbf{I}$ by Cauchy-Schwarz inequality and obtain

$$\begin{array}{cc}\hfill 4\mathbf{I}& \le {\int }_{B_{\varrho }(x_o)}{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }|D{u}_{\epsilon }|,\mathrm{\nabla }|D{u}_{\epsilon }|){\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|{\eta }^2\mathrm{d}x\hfill \\ \hfill & +4{\int }_{B_{\varrho }(x_o)}{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }\eta ,\mathrm{\nabla }\eta )\frac{{\varphi }^2(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{{\varphi }^{\prime }(|D{u}_{\epsilon }|)}\mathrm{d}x\hfill \\ \hfill & \le {\int }_{B_{\varrho }(x_o)}{\mathcal{C}}_{\epsilon }(D{u}_{\epsilon })(\mathrm{\nabla }|D{u}_{\epsilon }|,\mathrm{\nabla }|D{u}_{\epsilon }|){\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|{\eta }^2\mathrm{d}x\hfill \\ \hfill & +4{\int }_{B_{\varrho }(x_o)}[\epsilon +\mathbf{\Lambda }(|D{u}_{\epsilon }|)]\frac{{\varphi }^2(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{{\varphi }^{\prime }(|D{u}_{\epsilon }|)}{|\mathrm{\nabla }\eta |}^2\mathrm{d}x.\hfill \end{array}$$

Graph

From the second to last inequality we used the upper bound from Lemma 2.7 to estimate the second integral. Inserting the results above and re-absorbing the first term from the right into the left, we find that

4.6 $$\begin{array}{cc}& {\int }_{B_{\varrho }(x_o)}{\mathcal{A}}_{\epsilon }(D{u}_{\epsilon })(D^2{u}_{\epsilon },D^2{u}_{\epsilon })\varphi (|D{u}_{\epsilon }|){\eta }^2\mathrm{d}x\hfill \\ \hfill & \le 4{\int }_{B_{\varrho }(x_o)}[\epsilon +\mathbf{\Lambda }(|D{u}_{\epsilon }|)+h_{\epsilon }(|D{u}_{\epsilon }|)]\frac{{\varphi }^2(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{{\varphi }^{\prime }(|D{u}_{\epsilon }|)}{|\mathrm{\nabla }\eta |}^2\mathrm{d}x\hfill \\ \hfill & +c{\int }_{B_{\varrho }(x_o)}{|f|}^2\frac{\varphi (|D{u}_{\epsilon }|)+{\varphi }^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|}{h_{\epsilon }(|D{u}_{\epsilon }|)}{\eta }^2\mathrm{d}x\hfill \end{array}$$

Graph

holds true with a constant $c=c(p)$ . We now choose $\varphi (t):={(t-1)}_{+}^p\tilde{\varphi }(t)$ , where $\tilde{\varphi }\in W_{\text{loc}}^{1,\infty }({\mathbb{R}}_{\ge 0},{\mathbb{R}}_{\ge 0})$ is non-decreasing. Note that

$$\begin{array}{c}\hfill {\varphi }^{\prime }(t)={(t-1)}_{+}^{p-1}[p\tilde{\varphi }(t)+{(t-1)}_{+}{\tilde{\varphi }}^{\prime }(t)].\end{array}$$

Graph

For $t\in [0,1+2\mu ]$ we compute

$$\begin{array}{cc}& [\epsilon +\mathbf{\Lambda }(t)+h_{\epsilon }(t)]\frac{{\varphi }^2(t)t}{{\varphi }^{\prime }(t)}\hfill \\ \hfill & \le 2[\epsilon +\mathbf{\Lambda }(t)]\frac{{(t-1)}_{+}^{p+1}t{\tilde{\varphi }}^2(t)}{p\tilde{\varphi }(t)+{(t-1)}_{+}{\tilde{\varphi }}^{\prime }(t)}\hfill \\ \hfill & \le 2[\epsilon {(t-1)}_{+}^{p+1}t+max\{{(t-1)}_{+}^{2p},(p-1){(t-1)}_{+}^{2p-1}t\}]\frac{{\tilde{\varphi }}^2(t)}{p\tilde{\varphi }(t)+{(t-1)}_{+}{\tilde{\varphi }}^{\prime }(t)}\hfill \\ \hfill & \le c[\epsilon {\mu }^{p+2}+{\mu }^{2p}]\frac{{\tilde{\varphi }}^2(t)}{p\tilde{\varphi }(t)+{(t-1)}_{+}{\tilde{\varphi }}^{\prime }(t)}\hfill \\ \hfill & \le c(p,M,\delta )\frac{{\mu }^{2p}{\tilde{\varphi }}^2(t)}{p\tilde{\varphi }(t)+{(t-1)}_{+}{\tilde{\varphi }}^{\prime }(t)}.\hfill \end{array}$$

Graph

In turn we used $\delta \le \mu \le M$ from (3.5) and (3.6), which implies on the one hand $t\le 1+2\mu \le (2+\frac{1}{\delta })\mu$ , and on the other hand ${\mu }^2\le max\{{\delta }^{2-p},M^{2-p}\}{\mu }^p$ . Next, we compute

$$\begin{array}{cc}\hfill \frac{\varphi (t)+{\varphi }^{\prime }(t)t}{h_{\epsilon }(t)}& \le \frac{\varphi (t)+{\varphi }^{\prime }(t)t}{h(t)}\hfill \\ \hfill & =[{(t-1)}_{+}\tilde{\varphi }(t)+pt\tilde{\varphi }(t)+{(t-1)}_{+}{\tilde{\varphi }}^{\prime }(t)t]t\hfill \\ \hfill & \le c(p,M,\delta )[\tilde{\varphi }(t)+{\tilde{\varphi }}^{\prime }(t)t].\hfill \end{array}$$

Graph

Due to assumptions (3.4) and (3.6) we know that $|D{u}_{\epsilon }|\le 1+2\mu$ on $B_{\varrho }(x_o)$ . This allows us to use the preceding estimates in (4.6) to bound the right-hand side from above. Moreover, by Lemma 2.11 the left-hand side in (4.6) can be estimated from below. Proceeding in this way we obain

4.7 $$\begin{array}{cc}& {\int }_{B_{\varrho }(x_o)}|D[g(|D{u}_{\epsilon }|)D{u}_{\epsilon }]{|}^2\tilde{\varphi }(|D{u}_{\epsilon }|){\eta }^2\mathrm{d}x\hfill \\ \hfill & \le c{\int }_{B_{\varrho }(x_o)}\frac{{\mu }^{2p}{\tilde{\varphi }}^2(|D{u}_{\epsilon }|)}{p\tilde{\varphi }(|D{u}_{\epsilon }|)+(|D{u}_{\epsilon }{|-1)}_{+}{\tilde{\varphi }}^{\prime }(|D{u}_{\epsilon }|)}{|\mathrm{\nabla }\eta |}^2\mathrm{d}x\hfill \\ \hfill & +c{\int }_{B_{\varrho }(x_o)}{|f|}^2[\tilde{\varphi }(|D{u}_{\epsilon }|)+{\tilde{\varphi }}^{\prime }(|D{u}_{\epsilon }|)|D{u}_{\epsilon }|]{\eta }^2\mathrm{d}x,\hfill \end{array}$$

Graph

for any $\eta \in C_0^1(B_{\varrho }(x_o))$ . The constant c depends only on p, M, and $\delta$ .

Different concrete choices of $\tilde{\varphi }$ in (4.7) result in two important energy inequalities. The first one is

Lemma 4.4

Let $\epsilon \in (0,1]$ and $u_{\epsilon }\in W^{1,p}(B_R,{\mathbb{R}}^N)$ be a weak solution of the regularized system (3.1) such that hypotheses (3.4), (3.5) and (3.6) are in force on $B_{\varrho }(x_o)\subset B_{r_o}⋐B_R$ . Then, for any $\tau \in (0,1)$ there holds

$$\begin{array}{c}\hfill {\int }_{B_{\tau \varrho }(x_o)}|D[g(|D{u}_{\epsilon }|)D{u}_{\epsilon }]{|}^2\mathrm{d}x\le c[\frac{{\mu }^{2p}}{{\varrho }^2{(1-\tau )}^2}+{\varrho }^{-\frac{2n}{n+\sigma }}{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2]|B_{\varrho }|\end{array}$$

Graph

for some universal constant $c=c(n,p,M,\delta )$ .

Proof

We apply inequality (4.7) with the choice $\tilde{\varphi }\equiv 1$ . The cut-off function $\eta \in C_0^1(B_{\varrho }(x_o))$ is chosen such that $\eta \equiv 1$ in $B_{\tau \varrho }(x_o)$ , $0\le \eta \le 1$ , and $|\mathrm{\nabla }\eta |\le \frac{2}{(1-\tau )\varrho }$ . This leads us to

$$\begin{array}{cc}\hfill {\int }_{B_{\tau \varrho }(x_o)}|D[g(|D{u}_{\epsilon }|)D{u}_{\epsilon }]{|}^2\mathrm{d}x& \le c{\int }_{B_{\varrho }(x_o)}[{\mu }^{2p}{|\mathrm{\nabla }\eta |}^2+{|f|}^2]\mathrm{d}x\hfill \\ \hfill & \le \frac{c{\mu }^{2p}}{{(1-\tau )}^2}{\varrho }^{n-2}+c{\varrho }^{n-\frac{2n}{n+\sigma }}{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2\hfill \end{array}$$

Graph

with a constant $c=c(n,p,M,\delta )$ , which is the claimed energy estimate. $\square$

The second energy estimate is

Lemma 4.5

Let $\nu \in (0,\frac{1}{4}]$ , $\epsilon \in (0,1]$ and $u_{\epsilon }\in W^{1,p}(B_R,{\mathbb{R}}^N)$ be a weak solution of the regularized system (3.1) such that hypotheses (3.4), (3.5), (3.6) and (3.7) are in force on $B_{\varrho }(x_o)\subset B_{r_o}⋐B_R$ . Then, for any $\tau \in (0,1)$ we have

$$\begin{array}{c}\hfill {\int }_{E_{\tau \varrho }^{\nu }(x_o)}|D[g(|D{u}_{\epsilon }|)D{u}_{\epsilon }]{|}^2\mathrm{d}x\le c[\frac{{\mu }^{2p}\nu }{{\varrho }^2{(1-\tau )}^2}+\frac{{\varrho }^{-\frac{2n}{n+\sigma }}}{\nu }{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2]|B_{\varrho }|,\end{array}$$

Graph

for a constant $c=c(n,p,M,\delta )$ .

Proof

This time we choose

$$\begin{array}{c}\hfill \tilde{\varphi }(t)={(t-1-\delta -k)}_{+}^2\text{with}k:=(1-2\nu )\mu \end{array}$$

Graph

in inequality (4.7), and obtain

$$\begin{array}{cc}& {\int }_{B_{\varrho }(x_o)}{\eta }^2|D[g(|D{u}_{\epsilon }|)D{u}_{\epsilon }]{|}^2(|D{u}_{\epsilon }|-1-\delta -k{)}_{+}^2\mathrm{d}x\hfill \\ \hfill & \le c{\int }_{B_{\varrho }(x_o)}\frac{{\mu }^{2p}(|D{u}_{\epsilon }|-1-\delta -k{)}_{+}^3}{p(|D{u}_{\epsilon }|-1-\delta -k)+2(|D{u}_{\epsilon }|-1)}{|\mathrm{\nabla }\eta |}^2\mathrm{d}x\hfill \\ \hfill & +c{\int }_{B_{\varrho }(x_o)}{|f|}^2|D{u}_{\epsilon }|(|D{u}_{\epsilon }|-1-\delta -k{)}_{+}{\eta }^2\mathrm{d}x.\hfill \end{array}$$

Graph

On $B_{\varrho }(x_o)\cap \{|D{u}_{\epsilon }|>1+\delta +k\}$ we have

$$\begin{array}{c}\hfill (|D{u}_{\epsilon }|-1-\delta -k{)}_{+}\le \mu -k=\mu -(1-2\nu )\mu =2\nu \mu ,\end{array}$$

Graph

and

$$\begin{array}{c}\hfill p(|D{u}_{\epsilon }|-1-\delta -k)+2(|D{u}_{\epsilon }|-1)\ge 2(\delta +k)\ge 2k\ge \mu ,\end{array}$$

Graph

since $\nu \le \frac{1}{4}$ . Again, we choose $\eta \in C_0^1(B_{\varrho }(x_o))$ to be a non-negative cut-off function with $\eta \equiv 1$ in $B_{\tau \varrho }(x_o)$ , $0\le \eta \le 1$ , and $|\mathrm{\nabla }\eta |\le \frac{2}{(1-\tau )\varrho }$ . This, together with the fact that $|D{u}_{\epsilon }|\le 1+2\mu$ on $B_{\varrho }(x_o)$ , allows us to estimate the right-hand side in the above inequality. Indeed, we have

$$\begin{array}{cc}& {\int }_{B_{\tau \varrho }(x_o)}|D[g(|D{u}_{\epsilon }|)D{u}_{\epsilon }]{|}^2(|D{u}_{\epsilon }|-1-\delta -k{)}_{+}^2\mathrm{d}x\hfill \\ \hfill & \le c[\frac{{\nu }^3{\mu }^{2p+2}}{{(1-\tau )}^2}{\varrho }^{n-2}+\nu {\mu }^2{\int }_{B_{\varrho }(x_o)}{|f|}^2\mathrm{d}x].\hfill \end{array}$$

Graph

Therefore it remains to estimate the left-hand side from below. The integral has to be taken only on the set of points $x\in B_{\varrho }(x_o)$ with $|D{u}_{\epsilon }(x)|-1-\delta >k=(1-2\nu )\mu$ . We shrink this set to those points satisfying the stronger condition $|D{u}_{\epsilon }(x)|-1-\delta >(1-\nu )\mu >k$ , i.e. to $E_{\tau \varrho }^{\nu }(x_o)$ . On this set we have

$$\begin{array}{c}\hfill |D{u}_{\epsilon }|-1-\delta -k\ge (1-\nu )\mu -(1-2\nu )\mu =\nu \mu .\end{array}$$

Graph

Inserting this above we conclude that

$$\begin{array}{c}\hfill {\nu }^2{\mu }^2{\int }_{E_{\tau \varrho }^{\nu }(x_o)}|D[g(|D{u}_{\epsilon }|)D{u}_{\epsilon }]{|}^2\mathrm{d}x\le c[\frac{{\nu }^3{\mu }^{2p+2}}{{\varrho }^2{(1-\tau )}^2}+\nu {\mu }^2{\varrho }^{-\frac{2n}{n+\sigma }}{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2]|B_{\varrho }|\end{array}$$

Graph

holds true. This proves the claim. $\square$

The non-degenerate regime

The aim of this section is to prove Proposition 3.4. Throughout this section we presume the following general assumptions. For given $\epsilon \in (0,1]$ we denote by $u_{\epsilon }\in W^{1,p}(B_R,{\mathbb{R}}^N)$ the unique weak solution of the regularized system (3.1). Moreover, we assume that for some $\delta \in (0,1]$ and $\mu >\delta$ and a ball $B_{2\varrho }(x_o)\subset B_{r_1}$ with $\varrho \le 1$ assumptions (3.4)–(3.6) are in force. We denote by

5.1

Graph

the $L^2$ -excess of $D{u}_{\epsilon }$ on $B_{\varrho }(x_o)$ , i.e. the $L^2$ -mean square deviation of $D{u}_{\epsilon }$ from its mean value ${(D{u}_{\epsilon })}_{x_o,\varrho }$ .

Higher integrability

An ingredient in the proof of Proposition 3.4 is the following higher integrability result.

Lemma 5.1

Under the general assumptions of Sect. 5 there exist $\vartheta =\vartheta (n,p,\sigma ,M,\delta )\in (0,min\{\frac{1}{2},\frac{n+\sigma }{2}-1\}]$ and $c=c(n,p,M,\delta )$ such for any $\xi \in {\mathbb{R}}^{\mathit{Nn}}$ satisfying

$$\begin{array}{c}\hfill 1+\delta +\frac{1}{4}\mu \le |\xi |\le 1+\delta +\mu ,\end{array}$$

Graph

we have

Graph

Proof

We consider a ball $B_s(z_o)\subset B_{\varrho }(x_o)$ . We test the weak form (3.2) of the elliptic system by the testing function

$$\begin{array}{c}\hfill \phi :={\eta }^{2}w,\text{where}w:=u_{\epsilon }-{(u_{\epsilon })}_{z_o,s}-\xi (x-z_o)\end{array}$$

Graph

and $\eta \in C_0^1(B_s(z_o))$ is a standard cut-off function with $\eta \equiv 1$ in $B_{s/2}(z_o)$ , $0\le \eta \le 1$ , and $|\mathrm{\nabla }\eta |\le \frac{4}{s}$ . We obtain

$$\begin{array}{cc}\hfill 0& ={\int }_{B_s(z_o)}[{\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })·D\phi +f·\phi ]dx\hfill \\ \hfill & ={\int }_{B_s(z_o)}[({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(\xi ))·D\phi +f·\phi ]\mathrm{d}x\hfill \\ \hfill & ={\int }_{B_s(z_o)}[({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(\xi ))·[{\eta }^{2}Dw+2\eta \mathrm{\nabla }\eta \otimes w]+f·\phi ]\mathrm{d}x.\hfill \end{array}$$

Graph

We use the monotonicity of ${\mathbf{A}}_{\epsilon }$ from Lemma 2.8 in order to estimate the first term from below. Due to our assumption on $\xi$ and (3.4) we have $|\xi |+|D{u}_{\epsilon }|\le 5|\xi |$ and therefore obtain

$$\begin{array}{cc}\hfill {\int }_{B_s(z_o)}& {\eta }^2({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(\xi ))·Dw\mathrm{d}x\ge [\epsilon +\lambda \frac{{(|\xi |-1)}^p}{{|\xi |}^2}]{\int }_{B_s(z_o)}{\eta }^2{|Dw|}^2\mathrm{d}x,\hfill \end{array}$$

Graph

where $\lambda =\frac{1}{5·2^{p+1}}min\{1,p-1\}$ . To bound the second integral from above we use the structural upper bound from Lemma 2.8. Moreover, we observe that $(|D{u}_{\epsilon }{|-1)}_{+}\le 4(|\xi |-1)$ due to our assumption on $\xi$ . This allows us to estimate

$$\begin{array}{cc}& 2|{\int }_{B_s(z_o)}\eta ({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(\xi ))·\mathrm{\nabla }\eta \otimes w\mathrm{d}x|\hfill \\ \hfill & \le c[\epsilon +{(|\xi |-1)}^{p-2}]{\int }_{B_s(z_o)}\eta |Dw||w||\mathrm{\nabla }\eta |dx\hfill \\ \hfill & \le [\epsilon +\frac{1}{2}\lambda \frac{{(|\xi |-1)}^p}{{|\xi |}^2}]{\int }_{B_s(z_o)}{\eta }^2{|Dw|}^2\mathrm{d}x\hfill \\ \hfill & +c[\epsilon +{|\xi |}^2{(|\xi |-1)}^{p-4}]{\int }_{B_s(z_o)}{|w|}^2{|\mathrm{\nabla }\eta |}^2\mathrm{d}x.\hfill \end{array}$$

Graph

for a constant $c=c(p)$ . Re-absorbing terms on the left-hand side, we find that

$$\begin{array}{cc}& \frac{1}{2}\lambda \frac{{(|\xi |-1)}^p}{{|\xi |}^2}{\int }_{B_s(z_o)}{\eta }^2{|Dw|}^2\mathrm{d}x\hfill \\ \hfill & \le c[1+{|\xi |}^2{(|\xi |-1)}^{p-4}]{\int }_{B_s(z_o)}{|w|}^2{|D\eta |}^2\mathrm{d}x+{\int }_{B_s(z_o)}|f||w|\mathrm{d}x,\hfill \end{array}$$

Graph

where $c=c(p)$ . Due to the assumption on $\xi$ , the assumption (3.5), and the particular choice of $\eta$ , we conclude with an application of Hölder's and Sobolev–Poincaré's inequality a reverse Hölder inequality of the form

Graph

with a constant $c=c(n,p,M,\delta )$ . The dependence of c upon M only occurs in the sub-quadratic case of $p<2$ . The claim, i.e. the higher integrability, now follows with Gehring's lemma, since $Dw=D{u}_{\epsilon }-\xi$ , cf. [[1], Theorem 3.22], [[22], Theorem 2.4] and [[42], Theorem 3.3]. Note that $\vartheta$ can always be diminished if necessary. $\square$

Comparison with a linear system

In this section we will consider the weak solution $v\in u_{\epsilon }+W_0^{1,2}(B_{\varrho /2}(x_o),{\mathbb{R}}^N)$ of the linear elliptic system

5.2 $$\begin{array}{c}\hfill {\int }_{B_{\varrho /2}(x_o)}{\mathcal{B}}_{\epsilon }({(D{u}_{\epsilon })}_{x_o,\varrho /2})(Dv,D\phi )\mathrm{d}x=0,\end{array}$$

Graph

for any $\phi \in W_0^{1,2}(B_{\varrho /2}(x_o),{\mathbb{R}}^N)$ as comparison function to our solution $u_{\epsilon }$ of the regularized elliptic system (3.1). Recall that ${\mathcal{B}}_{\epsilon }$ has been defined in (2.7).

Lemma 5.2

Let the general assumptions of Sect. 5 be in force and assume that

5.3 $$\begin{array}{c}\hfill 1+\delta +\frac{1}{4}\mu \le |{(D{u}_{\epsilon })}_{x_o,\varrho }|\le 1+\delta +\mu .\end{array}$$

Graph

Then, there exists $\vartheta =\vartheta (n,p,M,\sigma ,\delta )\in (0,min\{\frac{1}{2},\frac{n+\sigma }{2}-1\}]$ and $c=c(n,p,M,\delta$ , ${\Vert f\Vert }_{L^{n+\sigma }(B_R)})$ such that

Graph

Here, v is the unique weak solution of the Dirichlet problem (5.2) and $\beta =\frac{\sigma }{n+\sigma }$ .

Proof

Throughout the proof we omit the reference to the center $x_o$ and write $B_{\varrho }$ instead of $B_{\varrho }(x_o)$ . Moreover, we abbreviate $\xi :={(D{u}_{\epsilon })}_{\varrho }$ . Using the weak form (3.2) of the elliptic system we obtain

$$\begin{array}{cc}\hfill 0& ={\int }_{B_{\varrho /2}}{\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })·D\phi \mathrm{d}x+{\int }_{B_{\varrho /2}}f·\phi \mathrm{d}x\hfill \\ \hfill & ={\int }_{B_{\varrho /2}}({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(\xi ))·D\phi \mathrm{d}x+{\int }_{B_{\varrho /2}}f·\phi \mathrm{d}x,\hfill \end{array}$$

Graph

for any $\phi \in W_0^{1,2}(B_{\varrho /2},{\mathbb{R}}^N)$ . Using also the fact that v is a weak solution of the linear elliptic system (5.2), we find that

$$\begin{array}{cc}& {\int }_{B_{\varrho /2}}{\mathcal{B}}_{\epsilon }(\xi )(D{u}_{\epsilon }-Dv,D\phi )\mathrm{d}x\hfill \\ \hfill & ={\int }_{B_{\varrho /2}}{\mathcal{B}}_{\epsilon }(\xi )(D{u}_{\epsilon },D\phi )\mathrm{d}x\hfill \\ \hfill & ={\int }_{B_{\varrho /2}}{\mathcal{B}}_{\epsilon }(\xi )(D{u}_{\epsilon }-\xi ,D\phi )\mathrm{d}x\hfill \\ \hfill & ={\int }_{B_{\varrho /2}}[{\mathcal{B}}_{\epsilon }(\xi )(D{u}_{\epsilon }-\xi ,D\phi )-({\mathbf{A}}_{\epsilon }(D{u}_{\epsilon })-{\mathbf{A}}_{\epsilon }(\xi ))·D\phi ]\mathrm{d}x-{\int }_{B_{\varrho /2}}f·\phi \mathrm{d}x\hfill \\ \hfill & \le c(p){\mu }^{p-3}{\int }_{B_{\varrho /2}}|D{u}_{\epsilon }{-\xi |}^2|D\phi |\mathrm{d}x+({\int }_{B_{\varrho /2}}{|f|}^2\mathrm{d}x{)}^{\frac{1}{2}}({\int }_{B_{\varrho /2}}{|\phi |}^2\mathrm{d}x{)}^{\frac{1}{2}}.\hfill \end{array}$$

Graph

Here we used from the second to last line Lemma 2.10. This is possible since (3.4) and (5.3) are in force. Since $u_{\epsilon }-v\in W_0^{1,2}(B_{\varrho /2},{\mathbb{R}}^N)$ , we may choose the testing function $\phi =u_{\epsilon }-v$ . Together with the bound from below from Lemma 2.7 and Hölder's and Poincaré's inequality this leads us to

$$\begin{array}{cc}& \gamma {\mu }^{p-2}{\int }_{B_{\varrho /2}}{|D{u}_{\epsilon }-Dv|}^2\mathrm{d}x\hfill \\ \hfill & \le {\int }_{B_{\varrho /2}}{\mathcal{B}}_{\epsilon }(\xi )(D{u}_{\epsilon }-Dv,D{u}_{\epsilon }-Dv)\mathrm{d}x\hfill \\ \hfill & \le c(p){\mu }^{p-3}{\int }_{B_{\varrho /2}}|D{u}_{\epsilon }{-\xi |}^2|D{u}_{\epsilon }-Dv|\mathrm{d}x\hfill \\ \hfill & +({\int }_{B_{\varrho /2}}{|f|}^2\mathrm{d}x{)}^{\frac{1}{2}}({\int }_{B_{\varrho /2}}{|u_{\epsilon }-v|}^2\mathrm{d}x{)}^{\frac{1}{2}}\hfill \\ \hfill & \le c(p){\mu }^{p-3}({\int }_{B_{\varrho /2}}|D{u}_{\epsilon }{-\xi |}^4\mathrm{d}x{)}^{\frac{1}{2}}({\int }_{B_{\varrho /2}}{|D{u}_{\epsilon }-Dv|}^2\mathrm{d}x{)}^{\frac{1}{2}}\hfill \\ \hfill & +c(n)\varrho ({\int }_{B_{\varrho /2}}{|f|}^2\mathrm{d}x{)}^{\frac{1}{2}}({\int }_{B_{\varrho /2}}{|D{u}_{\epsilon }-Dv|}^2\mathrm{d}x{)}^{\frac{1}{2}}\hfill \end{array}$$

Graph

for a constant $\gamma =\gamma (p,\delta )>0$ . We divide both sides by

$$\begin{array}{c}\hfill \gamma {\mu }^{p-2}[{\int }_{B_{\varrho /2}}{|D{u}_{\epsilon }-Dv|}^2\mathrm{d}x{]}^{\frac{1}{2}},\end{array}$$

Graph

square the result and finally take means. This implies

Graph

with a constant $c=c(n,p,M,\delta )$ . Here we have also used $\delta <\mu \le M$ . At this stage arrived, we want to reduce the integrability exponent on the right-hand side from 4 to $2(1+\vartheta )$ , where $\vartheta =\vartheta (n,p,M,\delta )\in (0,min\{\frac{1}{2},\frac{n+\sigma }{2}-1\}]$ is the integrability exponent from the higher integrability Lemma 5.1. This is possible since $|D{u}_{\epsilon }|$ and $|\xi |$ are bounded by $(2+\frac{1}{\delta })\mu$ on account of (3.4) and (3.6). Then, the application of the higher integrability lemma yields

Graph

where $c=c(n,p,M,\delta )$ . Inserting this inequality above and noting that $\varrho \le 1$ finishes the proof of the lemma. $\square$

The following a priori estimate for solutions to linear elliptic systems can be inferred from [[5]] once the ellipticity conditions for the quadratic form ${\mathcal{B}}_{\epsilon }({(D{u}_{\epsilon })}_{x_o,\varrho })$ are established; see also [[21], Theorem 2.3].

Lemma 5.3

Let the general assumptions of Sect. 5 be in force and assume that (5.3) holds true. Then, the weak solution $v\in W^{1,2}(B_{\varrho /2}(x_o),{\mathbb{R}}^N)$ of the linear elliptic system (5.2) satisfies $v\in W_{\text{loc}}^{2,2}(B_{\varrho /2}(x_o),{\mathbb{R}}^N)$ and there exists a constant $c_o=c_o(n,N,p,\delta )$ such that for any $\tau \in (0,\frac{1}{2}]$ we have

Graph

Proof

As mentioned before, the a priori estimate is standard. The constant $c_o$ depends on the dimensions n, N and the ellipticity constant and the upper bound of the quadratic form ${\mathcal{B}}_{\epsilon }({(D{u}_{\epsilon })}_{x_o,\varrho })$ . Due to assumption (5.3) and Lemma 2.7 these quantities only depend on p and $\delta$ . $\square$

Exploiting the measure theoretic information

The aim of this subsection is to convert the measure theoretic information (3.7) into a lower bound for the mean value of $D{u}_{\epsilon }$ and smallness of the excess.

Lemma 5.4

Let the general assumptions of Sect. 5 be in force. Furthermore, assume that (3.7) holds for some $\nu \in (0,\frac{1}{4}]$ . Then there exists a constant $c=c(n,p,M,\delta )$ such that for any $\tau \in [\frac{1}{2},1)$ there holds

$$\begin{array}{c}\hfill \mathrm{\Phi }(x_o,\tau \varrho )\le c{\mu }^2[\frac{{\nu }^{\frac{2}{n}}}{{(1-\tau )}^2}+\frac{{\varrho }^{2\beta }}{\nu }{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2].\end{array}$$

Graph

Proof

Throughout the proof we omit the reference to the center $x_o$ and write $B_{\varrho }$ instead of $B_{\varrho }(x_o)$ . We define $\zeta \in {\mathbb{R}}^{\mathit{Nn}}$ by

$$\begin{array}{c}\hfill {|\zeta |}^{p-1}\zeta =({|\mathcal{G}(D{u}_{\epsilon })|}^{p-1}\mathcal{G}(D{u}_{\epsilon }){)}_{\tau \varrho }\end{array}$$

Graph

and let

$$\begin{array}{c}\hfill \tilde{\zeta }:={\mathcal{G}}^{-1}(\zeta ).\end{array}$$

Graph

Note that $|\zeta |\le \delta +\mu$ by (3.4) and $|\tilde{\zeta }|\le 1+\delta +\mu$ . Due to the minimality of the integral average ${(D{u}_{\epsilon })}_{\tau \varrho }$ with respect to the mapping $\xi \mapsto {\int }_{B_{\tau \varrho }}{|D{u}_{\epsilon }-\xi |}^2\mathrm{d}x$ , we have

Graph

We recall that $|D{u}_{\epsilon }|,|\tilde{\zeta }|\le 1+\delta +\mu$ and hence by (3.6) we have $|D{u}_{\epsilon }|,|\tilde{\zeta }|\le (2+\frac{1}{\delta })\mu$ . Due to assumption (3.7) we therefore obtain for the second integral

$$\begin{array}{c}\hfill \mathbf{II}\le \frac{c(\delta ){\mu }^2}{|B_{\tau \varrho }|}|B_{\varrho }\backslash E_{\varrho }^{\nu }|\le \frac{c(\delta )\nu {\mu }^2}{{\tau }^n}.\end{array}$$

Graph

For the estimate of $\mathbf{I}$ we first note that $|D{u}_{\epsilon }|\ge 1+\frac{3}{2}\delta$ on $E_{\tau \varrho }^{\nu }$ since $\nu \in (0,\frac{1}{4}]$ and $\mu \ge \delta$ . Therefore, the application of Lemma 2.3 yields

$$\begin{array}{c}\hfill \mathbf{I}\le \frac{c(\delta )}{|B_{\tau \varrho }|}{\int }_{E_{\tau \varrho }^{\nu }}|\mathcal{G}(D{u}_{\epsilon })-\zeta {|}^2\mathrm{d}x.\end{array}$$

Graph

Next, we note that

$$\begin{array}{c}\hfill |\mathcal{G}(D{u}_{\epsilon })|+|\zeta |\ge |\mathcal{G}(D{u}_{\epsilon })|=(|D{u}_{\epsilon }{|-1)}_{+}>\delta +(1-\nu )\mu >\frac{1}{2}\mu \text{on}{E}_{\tau \varrho }^{\nu }.\end{array}$$

Graph

Using this information, Lemma 2.2, the choice of $\zeta$ , and Poincaré's inequality we obtain

Graph

We once again decompose the domain of integration into $E_{\tau \varrho }^{\nu }$ and $B_{\tau \varrho }\backslash E_{\tau \varrho }^{\nu }$ . Subsequently applying Hölder's inequality and taking into account assumption (3.7) leads us to

$$\begin{array}{cc}\hfill \mathbf{I}& \le \frac{c(p,\delta ){\varrho }^2}{{\mu }^{2p-2}|B_{\tau \varrho }|}[{\int }_{E_{\tau \varrho }^{\nu }}|D[g(D{u}_{\epsilon })D{u}_{\epsilon }]{|}^2\mathrm{d}x+{\nu }^{\frac{2}{n}}{\int }_{B_{\tau \varrho }\backslash E_{\tau \varrho }^{\nu }}|D[g(D{u}_{\epsilon })D{u}_{\epsilon }]{|}^2\mathrm{d}x].\hfill \end{array}$$

Graph

We note that $\tau \ge \frac{1}{2}$ and hence $|B_{\tau \varrho }|\ge c(n){\varrho }^n$ . For the first integral we use Lemma 4.5, while for the second one we use Lemma 4.4 and the assumption $\mu \ge \delta$ . In this way we obtain

$$\begin{array}{cc}\hfill \mathbf{I}& \le \frac{c{\mu }^2}{{(1-\tau )}^2}[\nu +{\nu }^{\frac{2}{n}}]+\frac{c{\varrho }^{2-\frac{2n}{n+\sigma }}}{\nu {\mu }^{2p-2}}{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2\hfill \\ \hfill & \le c{\mu }^2[\frac{{\nu }^{\frac{2}{n}}}{{(1-\tau )}^2}+\frac{{\varrho }^{\frac{2\sigma }{n+\sigma }}}{\nu }{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2],\hfill \end{array}$$

Graph

for a constant $c=c(n,p,M,\delta )$ . Inserting this above yields the desired estimate. $\square$

Lemma 5.5

Let the general assumptions of Sect. 5 be in force. Then, for any $\theta \in (0,\frac{1}{64}]$ there exist $\nu =\nu (n,p,M,\delta ,\theta )\in (0,\frac{1}{4}]$ and ${\varrho }_o={\varrho }_o{(n,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},$ $M,\delta ,\theta )\in (0,1]$ such that the smallness assumption $\varrho \le {\varrho }_o$ and the measure theoretic hypothesis (3.7) imply

5.4 $$\begin{array}{c}\hfill |{(D{u}_{\epsilon })}_{x_o,\varrho }|\ge 1+\delta +\frac{1}{2}\mu \text{and}\mathrm{\Phi }(x_o,\varrho )\le \theta {\mu }^2.\end{array}$$

Graph

Proof

We let $\tau \in [\frac{1}{2},1)$ , $\nu \in (0,\frac{1}{4}]$ and ${\varrho }_o\in (0,1]$ . Consider $B_{\varrho }(x_o)\subset B_R$ with $\varrho \le {\varrho }_o$ . For convenience in notation we omit the reference to the center $x_o$ . Using the minimality of ${(Du)}_{\varrho }$ with respect to the mapping $\xi \mapsto {\int }_{B_{\varrho }}{|Du-\xi |}^2\mathrm{d}x$ and decomposing the domain of integration into $B_{\tau \varrho }$ and $B_{\varrho }\backslash B_{\tau \varrho }$ , we obtain

Graph

For the first integral we use Lemma 5.4 and obtain

$$\begin{array}{c}\hfill \mathbf{I}={\tau }^n\mathrm{\Phi }(\tau \varrho )\le c{\mu }^2[\frac{{\nu }^{\frac{2}{n}}}{{(1-\tau )}^2}+\frac{{\varrho }^{2\beta }}{\nu }{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2],\end{array}$$

Graph

where $c=c(n,p,M,\delta )$ . For the second integral we use $|D{u}_{\epsilon }|\le 1+\delta +\mu \le c(\delta )\mu$ and get

$$\begin{array}{cc}\hfill \mathbf{II}& \le \frac{4{(1+\delta +\mu )}^2|B_{\varrho }\backslash B_{\tau \varrho }|}{|B_{\varrho }|}\le c(\delta ){\mu }^2(1-{\tau }^n)\le c(n,\delta ){\mu }^2(1-\tau ),\hfill \end{array}$$

Graph

so that

$$\begin{array}{c}\hfill \mathrm{\Phi }(\varrho )\le c{\mu }^2[\frac{{\nu }^{\frac{2}{n}}}{{(1-\tau )}^2}+(1-\tau )+\frac{{\varrho }^{2\beta }}{\nu }{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2],\end{array}$$

Graph

for a constant $c=c(n,p,M,\delta )$ . Now, we first choose $\tau \in [\frac{1}{2},1)$ in dependence on $n,p,M,\delta$ and $\theta$ in such a way that $c(1-\tau )\le \frac{1}{3}\theta$ . Subsequently, we choose $\nu \in (0,\frac{1}{4}]$ in dependence on $n,p,M,\delta$ and $\theta$ such that

$$\begin{array}{c}\hfill \nu \le min\{(\frac{\theta {(1-\tau )}^2}{3c}{)}^{\frac{n}{2}},\frac{\delta }{4(1+\delta )}\}.\end{array}$$

Graph

Finally, we choose ${\varrho }_o\in (0,1]$ such that

$$\begin{array}{c}\hfill {\varrho }_o^{2\beta }\le \frac{\nu \theta }{3c(1+{\Vert f\Vert }_{L^{n+\sigma }(B_R)}^2)}.\end{array}$$

Graph

In this way we obtain (5.4) ${}_2$ .

To prove (5.4) ${}_1$ , we first observe that the measure theoretic assumption (3.7) implies

$$\begin{array}{c}\hfill |E_{\varrho }^{\nu }|>(1-\nu )|B_{\varrho }|.\end{array}$$

Graph

Hence, due to the definition of the set $E_{\varrho }^{\nu }$ , we obtain

$$\begin{array}{cc}\hfill {\int }_{B_{\varrho }}|D{u}_{\epsilon }|\mathrm{d}x& \ge {\int }_{E_{\varrho }^{\nu }}|D{u}_{\epsilon }|\mathrm{d}x\ge (1+\delta +(1-\nu )\mu )|E_{\varrho }^{\nu }|\hfill \\ \hfill & \ge (1-\nu )(1+\delta +(1-\nu )\mu )|B_{\varrho }|.\hfill \end{array}$$

Graph

On the other hand, due to (5.4) ${}_2$ , we have

Graph

so that

$$\begin{array}{c}\hfill |{(D{u}_{\epsilon })}_{\varrho }|\ge (1-\nu )(1+\delta +(1-\nu )\mu )-\sqrt{\theta }\mu .\end{array}$$

Graph

Due to the choice of $\nu$ and the fact that $n\ge 2$ we have $\nu \le {(\frac{1}{3}\theta )}^{\frac{n}{2}}\le \frac{1}{3}\theta \le \frac{1}{2}\sqrt{\theta }$ and $\nu \le \frac{\delta }{4(1+\delta )}$ . Together with the assumptions $\delta \le \mu$ and $\theta \le \frac{1}{64}$ we obtain

$$\begin{array}{cc}& -\nu (1+\delta )+{(1-\nu )}^2\mu -\sqrt{\theta }\mu -\frac{1}{2}\mu \ge -\frac{1}{4}\delta +(\frac{1}{2}-2\nu -\sqrt{\theta })\mu \hfill \\ \hfill & \ge [\frac{1}{4}-2\sqrt{\theta }]\mu \ge 0.\hfill \end{array}$$

Graph

Inserting this above yields the claim (5.4) ${}_1$ and finishes the proof of the lemma. $\square$

Proof of Proposition 3.4

Our aim in this subsection is to prove Proposition 3.4. We start with an excess-decay estimate for the excess $\mathrm{\Phi }(x_o,\varrho )$ of $D{u}_{\epsilon }$ .

Lemma 5.6

Assume that the general hypotheses of Sect. 5 are in force. Let $\tau \in (0,\frac{1}{2}]$ and $\vartheta =\vartheta (n,p,M,\delta )\in (0,\frac{1}{2}]$ be the exponent from Lemma 5.2. If

5.5 $$\begin{array}{c}\hfill |{(D{u}_{\epsilon })}_{x_o,\varrho }|\ge 1+\delta +\frac{1}{4}\mu \text{and}\mathrm{\Phi }(x_o,\varrho )\le {\tau }^{\frac{n+2}{\vartheta }}{\mu }^2,\end{array}$$

Graph

hold true, then we have the quantitative excess decay estimate

$$\begin{array}{c}\hfill \mathrm{\Phi }(x_o,\tau \varrho )\le c_{\ast }[{\tau }^2\mathrm{\Phi }(x_o,\varrho )+{\tau }^{-n}{\varrho }^{2\beta }{\mu }^2]\end{array}$$

Graph

with a constant $c_{\ast }=c_{\ast }{(n,N,p,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta )$ .

Proof

Throughout the proof we omit the reference to the center $x_o$ and write $B_{\varrho }$ instead of $B_{\varrho }(x_o)$ . By $v\in u_{\epsilon }+W_0^{1,2}(B_{\varrho },{\mathbb{R}}^N)$ we denote the unique weak solution of the linear elliptic system (5.2). For $\tau \in (0,\frac{1}{2}]$ , we have

Graph

In view of Lemma 5.3 we deduce

Graph

where $c=c(n,N,p,\delta )$ . Inserting this above and applying Lemma 5.2 and assumption (5.5) ${}_2$ , we end up with

Graph

Note that the constant $c_{\ast }$ depends on ${n,N,p,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M$ and $\delta$ . $\square$

Proof of Proposition 3.4

By $\vartheta =\vartheta (n,p,\sigma ,M,\delta )\in (0,min\{\frac{1}{2},\frac{n+\sigma }{2}-1\}]$ we denote the constant from Lemma 5.2 and by $c_{\ast }=c_{\ast }{(n,N,p,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta )$ the one from Lemma 5.6. For $\beta =\frac{\sigma }{n+\sigma }\in (0,1)$ we define $\tau \in (0,\frac{1}{8}]$ by

$$\begin{array}{c}\hfill \tau :=min\{\frac{1}{8},2^{-\frac{1}{\beta }},{(10c_{\ast })}^{-\frac{1}{2(1-\beta )}}\}.\end{array}$$

Graph

For the particular choice $\theta ={\tau }^{\frac{n+2}{\vartheta }}$ we let ${\varrho }_o={\varrho }_o{(n,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta )\in (0,1]$ be the radius from Lemma 5.5. Finally, we define

$$\begin{array}{c}\hfill \widehat{\varrho }:=min\{{\varrho }_o,{(2c_{\ast })}^{-\frac{1}{2\beta }}{\tau }^{1+\frac{2n+2}{\beta \vartheta }}\},\end{array}$$

Graph

so that $\widehat{\varrho }$ depends on ${n,N,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M,\delta$ . In the following we consider a ball $B_{2\varrho }(x_o)\subset B_{r_1}$ with $\varrho \le \widehat{\varrho }$ . As before, we omit the reference to the center $x_o$ and write $B_{\varrho }$ instead of $B_{\varrho }(x_o)$ . By $\nu =\nu (n,p,M,\delta ,\theta =\frac{1}{2}{\tau }^{\frac{n+2}{\vartheta }})\in (0,\frac{1}{4}]$ we denote the constant from Lemma 5.5 and assume that (3.7) is satisfied for this particular choice of $\nu$ . Note that by our choice of $\tau$ the parameter $\nu$ depends on ${n,N,p,\sigma ,\Vert f\Vert }_{L^{n+\sigma }(B_R)},M$ and $\delta$ . From Lemma 5.5 applied with $\theta ={\tau }^{\frac{n+2}{\vartheta }}$ we infer that

5.6 $$\begin{array}{c}\hfill |{(Du)}_{\varrho }|\ge 1+\delta +\frac{1}{2}\mu \text{and}\mathrm{\Phi }(\varrho )\le {\tau }^{\frac{n+2}{\vartheta }}{\mu }^2.\end{array}$$

Graph

By induction we shall prove that for any $i\in \mathbb{N}$ we have

$$\begin{array}{cc}& \mathrm{\Phi }({\tau }^i\varrho )\le {\tau }^{\frac{n+2}{\vartheta }}{\tau }^{2\beta i}{\mu }^2{(\text{I})}_i\hfill \end{array}$$

Graph

and

$$\begin{array}{cc}& |{(Du)}_{{\tau }^i\varrho }|\ge 1+\delta +[\frac{1}{2}-\frac{1}{8}\sum _{j=0}^{i-1}{2}^{-j}]\mu .{(\text{II})}_i\hfill \end{array}$$

Graph

For $i=1$ we can apply Lemma 5.6, since (5.6) ensures that the assumptions of the lemma are satisfied. Then, (I) ${}_1$ follows from Lemma 5.6, (5.6) ${}_2$ and our choices of $\tau$ and $\widehat{\varrho }$ , since

$$\begin{array}{cc}\hfill \mathrm{\Phi }(\tau \varrho )& \le c_{\ast }[{\tau }^2\mathrm{\Phi }(\varrho )+{\tau }^{-n}{\varrho }^{2\beta }{\mu }^2]\hfill \\ \hfill & \le \frac{1}{2}{\tau }^{2\beta }\mathrm{\Phi }(\varrho )+\frac{c_{\ast }}{{\tau }^n}{\varrho }^{2\beta }{\mu }^2\hfill \\ \hfill & \le {\tau }^{\frac{n+2}{\vartheta }}{\tau }^{2\beta }[\frac{1}{2}+\frac{c_{\ast }{\widehat{\varrho }}^{2\beta }}{{\tau }^{\frac{n+2}{\vartheta }+n+2\beta }}]{\mu }^2\hfill \\ \hfill & \le {\tau }^{\frac{n+2}{\vartheta }}{\tau }^{2\beta i}{\mu }^2.\hfill \end{array}$$

Graph

For the proof of (II) ${}_1$ we use (5.6) ${}_2$ and ${\tau }^{\frac{n+2}{\vartheta }}\le {\tau }^{n+2}$ to obtain

Graph

so that

$$\begin{array}{c}\hfill |{(Du)}_{\tau \varrho }-{(Du)}_{\varrho }|\le \tau \mu \le \frac{1}{8}\mu .\end{array}$$

Graph

Together with (5.6) ${}_1$ this implies (II) ${}_1$ .

Now, we consider $i>1$ and prove (I) ${}_i$ and (II) ${}_i$ assuming that (I) ${}_{i-1}$ and (II) ${}_{i-1}$ hold. From (I) ${}_{i-1}$ and (II) ${}_{i-1}$ we observe that the assumptions of Lemma 5.6 as formulated in (5.5) are satisfied on $B_{{\tau }^{i-1}\varrho }$ . Therefore, applying the lemma with ${\tau }^{i-1}\varrho$ instead of $\varrho$ , recalling the choices of $\tau$ and $\widehat{\varrho }$ and joining the result with (I) ${}_{i-1}$ yields

$$\begin{array}{cc}\hfill \mathrm{\Phi }({\tau }^i\varrho )& \le c_{\ast }[{\tau }^2\mathrm{\Phi }({\tau }^{i-1}\varrho )+{\tau }^{-n}{({\tau }^{i-1}\varrho )}^{2\beta }{\mu }^2]\hfill \\ \hfill & \le \frac{1}{2}{\tau }^{2\beta }\mathrm{\Phi }({\tau }^{i-1}\varrho )+\frac{c_{\ast }}{{\tau }^n}{({\tau }^{i-1}\varrho )}^{2\beta }{\mu }^2\hfill \\ \hfill & \le {\tau }^{\frac{n+2}{\vartheta }}{\tau }^{2\beta i}[\frac{1}{2}+\frac{c_{\ast }{\widehat{\varrho }}^{2\beta }}{{\tau }^{\frac{n+2}{\vartheta }+n+2\beta }}]{\mu }^2\hfill \\ \hfill & \le {\tau }^{\frac{n+2}{\vartheta }}{\tau }^{2\beta i}{\mu }^2.\hfill \end{array}$$

Graph

This proves (I) ${}_i$ . Moreover, from (I) ${}_{i-1}$ and ${\tau }^{\frac{n+2}{\vartheta }}\le {\tau }^{n+2}$ we obtain

Graph

so that

$$\begin{array}{c}\hfill {|(Du)}_{{\tau }^i\varrho }-{(Du)}_{{\tau }^{i-1}\varrho }|\le {\tau }^{\beta (i-1)}\tau \mu \le \frac{1}{8}{2}^{-(i-1)}\mu ,\end{array}$$

Graph

by our choice of $\tau$ . Together with (II) ${}_{i-1}$ , this proves (II) ${}_i$ .

We now come to the proof of (3.8) and (3.9). For $i\in \mathbb{N}$ we obtain from the minimizing property of the mean value, Lemma 2.3, (I) ${}_i$ , (5.6) and our choice of $\tau$ that

5.7

Graph

This allows us to compute

Graph

Given $j<k$ , we use the preceding inequality to conclude that

5.8 $$\begin{array}{cc}\hfill |({\mathcal{G}}_{2\delta }(Du){)}_{{\tau }^j\varrho }-({\mathcal{G}}_{2\delta }(Du){)}_{{\tau }^k\varrho }|& \le \sum _{i=j+1}^k|({\mathcal{G}}_{2\delta }(Du){)}_{{\tau }^i\varrho }-({\mathcal{G}}_{2\delta }(Du){)}_{{\tau }^{i-1}\varrho }|\hfill \\ \hfill & \le {\tau }^{\frac{n}{2}+1}\sum _{i=j+1}^k{\tau }^{\beta (i-1)}\mu \le {\tau }^{\frac{n}{2}+1}\frac{{\tau }^{\beta j}}{1-{\tau }^{\beta }}\mu \hfill \\ \hfill & \le 2{\tau }^{\frac{n}{2}+1}{\tau }^{\beta j}\mu .\hfill \end{array}$$

Graph

This shows that ${({({\mathcal{G}}_{2\delta }(Du))}_{{\tau }^i\varrho })}_{i=1}^{\infty }$ is a Cauchy sequence and therefore the limit

$$\begin{array}{c}\hfill {\mathrm{\Gamma }}_{x_o}:=\underset{i\to \infty }{lim}({\mathcal{G}}_{2\delta }(Du){)}_{{\tau }^i\varrho }\end{array}$$

Graph

exists. Passing to the limit $k\to \infty$ in (5.8) yields

$$\begin{array}{c}\hfill |({\mathcal{G}}_{2\delta }(Du){)}_{{\tau }^j\varrho }-{\mathrm{\Gamma }}_{x_o}|\le 2{\tau }^{\frac{n}{2}+1}{\tau }^{\beta j}\mu \text{for}\text{any}j\in \mathbb{N}.\end{array}$$

Graph

Joining this with (5.7), we find

Graph

For $r\in (0,\varrho ]$ there exists $j\in {\mathbb{N}}_0$ such that ${\tau }^{j+1}\varrho <r\le {\tau }^j\varrho$ . Then, we obtain from the last inequality

Graph

This implies

Graph

so that also

$$\begin{array}{c}\hfill {\mathrm{\Gamma }}_{x_o}=\underset{r\downarrow 0}{lim}({\mathcal{G}}_{2\delta }(Du){)}_r.\end{array}$$

Graph

Finally, due to assumption (3.4) we have $|{({\mathcal{G}}_{2\delta }(Du))}_r|\le \mu$ for any $0<r\le \varrho$ , which implies $|{\mathrm{\Gamma }}_{x_o}|\le \mu$ . This finishes the proof of Proposition 3.4. $\square$

The degenerate regime

Our aim in this section is to prove Proposition 3.5, which treats the degenerate regime. The proof relies on a De Giorgi type reduction argument reducing the supremum of $U_{\epsilon }=(|D{u}_{\epsilon }{|-1-\delta )}_{+}^2$ under the measure theoretic assumption (3.10). The starting point is the energy estimate for $U_{\epsilon }$ from Lemma 4.3.

As in Section 5, we first formulate the general assumptions. For $\epsilon \in (0,1]$ we denote by $u_{\epsilon }\in W^{1,p}(B_R,{\mathbb{R}}^N)$ the unique weak solution to the Dirichlet problem (3.1) associated to the regularized system. We assume that (3.4) is in force for some $\mu ,\delta >0$ on some ball $B_{2\varrho }(x_o)\subset B_{r_1}⋐B_R$ . Let $U_{\epsilon }:=(|D{u}_{\epsilon }{|-1-\delta )}_{+}^2$ denote the function defined in (4.2). Note that (3.4) implies

$$\begin{array}{c}\hfill \underset{B_{2\varrho }(x_o)}{sup}{U}_{\epsilon }\le {\mu }^2.\end{array}$$

Graph

Moreover, we set $\beta :=\frac{\sigma }{n+\sigma }\in (0,1)$ .

We start by a De Giorgi type lemma for $U_{\epsilon }$ , which can for instance be deduced as in [[19], Chap. 10, Proposition 4.1] by the use of the energy estimate from Lemma 4.3. For the readers convenience we provide the proof in the appendix Sect. 7.

Lemma 6.1

(Reducing the supremum) Assume that the general assumptions of Sect. 6 are in force and let $\theta \in (0,1)$ . Then, there exists $\tilde{\nu }=\tilde{\nu }{(n,p,\Vert f\Vert }_{n+\sigma },M,\delta )\in (0,1)$ such that the measure theoretic assumption

$$\begin{array}{c}\hfill |\{x\in B_{\varrho }(x_o):U_{\epsilon }(x)>(1-\theta ){\mu }^2\}|<\tilde{\nu }|B_{\varrho }(x_o)|,\end{array}$$

Graph

implies that either

$$\begin{array}{c}\hfill {\mu }^2<\frac{{\varrho }^{\beta }}{\theta },\end{array}$$

Graph

or

$$\begin{array}{c}\hfill U_{\epsilon }\le (1-\frac{1}{2}\theta ){\mu }^2\text{in}{B}_{\varrho /2}(x_o)\end{array}$$

Graph

hold true. $\square$

The proof of the next Lemma can be deduced as in [[19], Chap. 10, Proposition 5.1] utilizing the energy estimate from Lemma 4.3; see also Sect. 7.

Lemma 6.2

Assume that the general assumptions of Sect. 6 are in force and assume that (3.10) is satisfied for some $\nu \in (0,1)$ . Then, for any $i_{\ast }\in \mathbb{N}$ we either have

$$\begin{array}{c}\hfill {\mu }^2<2^{i_{\ast }}{\varrho }^{\beta }/\nu \end{array}$$

Graph

or

$$\begin{array}{c}\hfill |\{x\in B_{\varrho }(x_o):U_{\epsilon }(x)>(1-2^{-i_{\ast }}\nu ){\mu }^2\}|<\frac{c_{\ast }}{\nu \sqrt{i_{\ast }}}|B_{\varrho }(x_o)|\end{array}$$

Graph

for a constant $c_{\ast }=c_{\ast }{(n,p,\Vert f\Vert }_{n+\sigma },M,\delta )$ .

Now, we have all the prerequisites at hand to provide the

Proof of Proposition 3.5

Let $\tilde{\nu }\in (0,1)$ and $c_{\ast }$ be the constants from Lemmas 6.1 and 6.2. Note that both depend on ${n,p,\Vert f\Vert }_{n+\sigma },M$ and $\delta$ . Choose $i_{\ast }\in \mathbb{N}$ such that

$$\begin{array}{c}\hfill i_{\ast }\ge (\frac{c_{\ast }}{\nu \tilde{\nu }}{)}^2.\end{array}$$

Graph

Then $i_{\ast }$ depends ${n,p,\Vert f\Vert }_{n+\sigma },M,\delta$ and $\nu$ . Lemma 6.2 implies that either ${\mu }^2<2^{i_{\ast }}{\varrho }^{\beta }/\nu$ , or

$$\begin{array}{c}\hfill |\{x\in B_{\varrho }(x_o):U_{\epsilon }(x)>(1-2^{-i_{\ast }}\nu ){\mu }^2\}|<\frac{c_{\ast }}{\nu \sqrt{i_{\ast }}}|B_{\varrho }(x_o)|\le \tilde{\nu }|B_{\varrho }(x_o)|.\end{array}$$

Graph

In the first case the proposition is proved with $c_o=2^{i_{\ast }}/\nu$ , while in the second case we may apply Lemma 6.1 with $\theta =2^{-i_{\ast }}\nu$ . Therefore either ${\mu }^2<2^{i_{\ast }}{\varrho }^{\beta }/\nu$ or

$$\begin{array}{c}\hfill U_{\epsilon }\le (1-2^{-(i_{\ast }+1)}\nu ){\mu }^2\text{in}{B}_{\varrho /2}(x_o).\end{array}$$

Graph

The first alternative coincides with the first alternative above, while the second one implies the sup-bound for ${\mathcal{G}}_{\delta }(D{u}_{\epsilon })$ for any $\kappa \ge \sqrt{1-2^{-(i_{\ast }+1)}\nu }$ since $U_{\epsilon }={|{\mathcal{G}}_{\delta }(D{u}_{\epsilon })|}^2$ . Therefore we my choose $\kappa \in [2^{-\beta /2},1)$ as required. $\square$

Acknowledgements

R. Giova and A. Passarelli di Napoli have been partially supported by the Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM). A. Passarelli di Napoli has been partially supported by Università degli Studi di Napoli Federico II through the Project FRA (000022-75-2021-FRA-PASSARELLI)

We would like to thank N. Liao for discussions concerning the modulus of continuity.

Funding Information

Open Access funding enabled and organized by Projekt DEAL.

Appendix

Here, we provide the proofs of Lemmas 6.1 and 6.2. We first state a well known iteration lemma which can be found for instance in [[19], Chap. 9.15.1, Lemma 15.1].

Lemma 7.1

Let ${(Y_i)}_{i\in {\mathbb{N}}_0}$ be a sequence of non-negative numbers satisfying

$$\begin{array}{c}\hfill Y_{i+1}\le \kappa b^i{Y}_i^{1+\gamma }\text{for}\text{all}i\in {\mathbb{N}}_0\end{array}$$

Graph

with some positive constants $\kappa ,\gamma$ and $b>1$ . If

$$\begin{array}{c}\hfill Y_0\le {\kappa }^{-\frac{1}{\gamma }}{b}^{-\frac{1}{{\gamma }^2}},\end{array}$$

Graph

then $Y_i\to 0$ as $i\to \infty$ .

The next lemma can be seen as a discrete version of the isoperimetric inequality, cf. [[19], Chap. 10.5.1, inequality (5.4)].

Lemma 7.2

Let $v\in W^{1,1}(B_{\varrho }(x_o))$ and $k,\ell \in \mathbb{R}$ with $k<\ell$ . Then, there exists a constant c depending on n such that

$$\begin{array}{c}\hfill (\ell -k)|B_{\varrho }(x_o)\cap \{v<k\}|\le \frac{c{\varrho }^{n+1}}{|B_{\varrho }(x_o)\cap \{v>\ell \}|}{\int }_{B_{\varrho }(x_o)\cap \{k<v<\ell \}}|Dv|\mathrm{d}x.\end{array}$$

Graph

Proof of Lemma 6.1

For $i\in \mathbb{N}$ we let

$$\begin{array}{c}\hfill B_i:=B_{{\varrho }_i}(x_o),\text{where}{\varrho }_i:=\frac{1}{2}\varrho (1+2^{-i}),\end{array}$$

Graph

and consider for the levels

$$\begin{array}{c}\hfill k_i:=(1-\frac{1}{2}\theta -2^{-(i+1)}\theta ){\mu }^2\end{array}$$

Graph

the normalized measure of the associated super-level sets

$$\begin{array}{c}\hfill {\mathit{Y}}_i:=\frac{|A_i|}{|B_i|},\text{where}{A}_i:=\{x\in B_i:U_{\epsilon }(x)>k_i\}.\end{array}$$

Graph

In view of Hölder's, Sobolev's inequality and Lemma 4.3 we now estimate[1]

$$\begin{array}{cc}\hfill {\int }_{B_{i+1}}{(U_{\epsilon }-k_i)}_{+}^2\mathrm{d}x& \le |A_i{|}^{\frac{2}{n}}[{\int }_{B_{i+1}}{(U_{\epsilon }-k_i)}_{+}^{\frac{2n}{n-2}}\mathrm{d}x{]}^{\frac{n-2}{n}}\hfill \\ \hfill & \le c|A_i{|}^{\frac{2}{n}}{\int }_{B_{i+1}}[|D{(U_{\epsilon }-k_i)}_{+}{|}^2+\frac{{(U_{\epsilon }-k_i)}_{+}^2}{{\varrho }_{i+1}^2}]\mathrm{d}x\hfill \\ \hfill & \le \frac{c{4}^i}{{\varrho }^2}|A_i{|}^{\frac{2}{n}}{\int }_{B_i}{(U_{\epsilon }-k_i)}_{+}^2\mathrm{d}x+{c\Vert f\Vert }_{n+\sigma }^2{|A_i|}^{1+\frac{2\beta }{n}}\hfill \\ \hfill & \le \frac{c{4}^i{\theta }^2{\mu }^4}{{\varrho }^2}|A_i{|}^{1+\frac{2}{n}}+c{|A_i|}^{1+\frac{2\beta }{n}},\hfill \end{array}$$

Graph

where $c=c(n,p,\Vert f{\Vert }_{n+\sigma },M,\delta )$ . On the other hand, by the definition of $k_i$ we have

$$\begin{array}{c}\hfill |A_{i+1}|\le \frac{4^{i+2}}{{\theta }^2{\mu }^4}{\int }_{B_{i+1}}{(U_{\epsilon }-k_i)}_{+}^2\mathrm{d}x.\end{array}$$

Graph

Here we used that $U_{\epsilon }-k_i\ge k_{i+1}-k_i$ on $A_{i+1}$ . Joining the preceding two inequalities and recalling the definition of ${\mathit{Y}}_i$ shows

$$\begin{array}{cc}\hfill {\mathit{Y}}_{i+1}& =\frac{|A_{i+1}|}{|B_{i+1}|}\hfill \\ \hfill & \le \frac{c{4}^{2i}{|A_i|}^{1+\frac{2\beta }{n}}}{|B_i{|}^{1+\frac{2\beta }{n}}}[\frac{|A_i{|}^{\frac{2(1-\beta )}{n}}{|B_i|}^{\frac{2\beta }{n}}}{{\varrho }^2}+\frac{|B_i{|}^{\frac{2\beta }{n}}}{{\theta }^2{\mu }^4}]\hfill \\ \hfill & \le c{4}^{2i}{\mathit{Y}}_i^{1+\frac{2\beta }{n}}[1+\frac{{\varrho }^{2\beta }}{{\theta }^2{\mu }^4}].\hfill \end{array}$$

Graph

If ${\mu }^2<{\varrho }^{\beta }/\theta$ , the lemma is proved. Otherwise, the term in brackets on the right-hand side is bounded by 2, so that

$$\begin{array}{c}\hfill {\mathit{Y}}_{i+1}\le c{4}^{2i}{\mathit{Y}}_i^{1+\frac{2}{n}},\end{array}$$

Graph

where $c=c(n,p,\Vert f{\Vert }_{n+\sigma },M,\delta )$ . Therefore, Lemma 7.1 ensures that $Y_i\to 0$ in the limit $i\to \infty$ , provided that ${\mathit{Y}}_0\le 2^{-n^2}{c}^{-\frac{n}{2}}=:\tilde{\nu }$ and hence $U_{\epsilon }\le (1-\frac{1}{2}\theta ){\mu }^2$ on $B_{\varrho /2}(x_o)$ . $\square$

Proof of Lemma 6.2

For convenience in notation we omit the reference to the center $x_o$ . For $i\in \{0,\cdots ,i_{\ast }\}$ we consider the super-level sets

$$\begin{array}{c}\hfill A_i:=\{x\in B_{\varrho }:U_{\epsilon }(x)>k_i\},\text{where}{k}_i:=(1-2^{-i}\nu ){\mu }^2.\end{array}$$

Graph

Then, due to assumption (3.10) we have $|B_{\varrho }\backslash A_i|\ge \nu |B_{\varrho }|$ for any $i\in \{0,\cdots ,i_{\ast }\}$ . Applying Lemma 7.2 to $U_{\epsilon }$ on $B_{\varrho }$ with $k=k_i$ and $\ell =k_{i+1}$ , we obtain with a constant $c=c(n)$ that

$$\begin{array}{cc}\hfill \frac{\nu {\mu }^2}{2^{i+1}}|A_{i+1}|& \le \frac{c{\varrho }^{n+1}}{|B_{\varrho }\backslash A_i|}{\int }_{A_i\backslash A_{i+1}}|D{U}_{\epsilon }|\mathrm{d}x\hfill \\ \hfill & \le \frac{c\varrho }{\nu }[{\int }_{B_{\varrho }}|D{(U_{\epsilon }-k_i)}_{+}{|}^2\mathrm{d}x{]}^{\frac{1}{2}}(|A_i|-|A_{i+1}|{)}^{\frac{1}{2}}.\hfill \end{array}$$

Graph

In view of the energy estimate from Lemma 4.3 and the definition of $k_i$ we have

$$\begin{array}{cc}\hfill {\int }_{B_{\varrho }}{|D{(U_{\epsilon }-k_i)}_{+}|}^2\mathrm{d}x& \le \frac{c}{{\varrho }^2}{\int }_{B_{2\varrho }}|{(U_{\epsilon }-k_i)}_{+}{|}^2\mathrm{d}x+c{\Vert f\Vert }_{n+\sigma }^2{\varrho }^{n-2+2\beta }\hfill \\ \hfill & \le \frac{c{\nu }^2{\mu }^4}{2^{2i}{\varrho }^2}[1+(\frac{2^i{\varrho }^{\beta }}{\nu {\mu }^2}{)}^2]|B_{\varrho }|\hfill \end{array}$$

Graph

with a constant $c=c(n,p,\Vert f{\Vert }_{n+\sigma },M,\delta )$ . If ${\mu }^2<2^i{\varrho }^{\beta }/\nu$ , the lemma is proved. Otherwise, the term in brackets on the right-hand side is bounded by 2. Inserting this above yields

$$\begin{array}{cc}\hfill |A_{i+1}{|}^2& \le \frac{c}{{\nu }^2}(|A_i|-|A_{i+1}|)|B_{\varrho }|,\hfill \end{array}$$

Graph

for $c=c(n,p,\Vert f{\Vert }_{n+\sigma },M,\delta )$ . Now, we add these inequalities for $i=0,1,\cdots ,i_{\ast }-1$ and obtain

$$\begin{array}{c}\hfill i_{\ast }|A_{i_{\ast }}{|}^2\le \frac{c}{{\nu }^2}|B_{\varrho }|\sum _{i=0}^{i_{\ast }-1}(|A_i|-|A_{i+1}|)\le \frac{c}{{\nu }^2}|A_0||B_{\varrho }|\le \frac{c}{{\nu }^2}{|B_{\varrho }|}^2.\end{array}$$

Graph

Therefore, we have

$$\begin{array}{c}\hfill |A_{i_{\ast }}|\le \frac{c}{\nu \sqrt{i_{\ast }}}|B_{\varrho }|.\end{array}$$

Graph

This proves the assertion of the lemma. $\square$

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

1 Björn A, Björn J. Nonlinear Potential Theory on Metric Spaces. EMS Tracts in Mathematics. 2011: Zürich; European Mathematical Society (EMS). 1231.31001 2 Brasco L. Global gradient estimates for solutions to a certain degenerate elliptic equation. Nonlinear Anal. 2011; 74; 2: 516-531. 2733227. 1202.35103 3 Brasco L, Carlier G. On certain anisotropic elliptic equations arising in congested optimal transport: local gradient bounds. Adv. Calc. Var. 2014; 7: 379-407. 3296496. 1305.49067 4 Brasco L, Carlier G, Santambrogio F. Congested traffic dynamics, weak flows and very degenerate elliptic equations. J. Math. Pures Appl. (9). 2010; 93; 6: 652-671. 2651987. 1193.35075 5 Campanato S. Equazioni ellittiche del II ordine e spazi. Ann. Mat. Pura Appl. 1965; 69: 321-381. 192168. 0145.36603 6 Carlier G, Jimenez C, Santambrogio F. Optimal transportation with traffic congestion and Wardrop equilibria. SIAM J. Control Optim. 2008; 47; 3: 1330-1350. 2407018. 1206.90209 7 Carlier G, Santambrogio F. A continuous theory of traffic congestion and Wardrop equilibria. Jo. Math. Sci. 2012; 181; 6: 792-804. 2870230. 1262.90016 8 Chipot M, Evans LC. Linearization at infinity and Lipschitz estimates for certain problems in calculus of variations. Proc. R. Soc. Edinburgh Sect. A. 1986; 102: 291-303. 0602.49029 9 Clop A, Giova R, Hatami F, Passarelli di Napoli A. Very degenerate elliptic equations under almost critical Sobolev regularity. Forum Math. 2020; 32; 6: 1515-1537. 4168702. 1460.35120 Colombo M, Figalli A. An excess-decay result for a class of degenerate elliptic equations. Discret. Contin. Dyn. Syst. Ser. S. 2014; 7; 4: 631-652. 3177756. 1286.35114 Colombo M, Figalli A. Regularity results for very degenerate elliptic equations. J. Math. Pures Appl. (9). 2014; 101; 1: 94-117. 3133426. 1282.35175 Cupini G, Giannetti F, Giova R, Passarelli di Napoli A. Higher integrability for minimizers of asymptotically convex integrals with discontinuous coefficients. Nonlinear Anal. 2017; 154: 7-24. 3614641. 1360.49027 Cupini G, Giannetti F, Giova R, Passarelli di Napoli A. Regularity results for vectorial minimizers of a class of degenerate convex integrals. J. Differ. Equ. 2018; 265; 9: 4375-4416. 3843304. 1411.49021 Cupini G, Guidorzi M, Mascolo E. Regularity of minimizers of vectorial integrals with growth. Nonlinear Anal. 2003; 54; 4: 591-616. 1983438. 1027.49032 De Giorgi E. Un esempio di estremali discontinue per un problema variazionale di tipo ellitico. Boll. UMI. 1968; 4: 135-137. 0155.17603 DiBenedetto D. Continuity of weak solutions to certain singular parabolic equations. Ann. Mat. Pura Appl. 1982; 130; 4: 131-176. 663969. 0503.35018 DiBenedetto D. local regularity of weak solutions of degenerate elliptic equations. Nonlinear Anal. 1983; 7; 8: 827-850. 709038. 0539.35027 DiBenedetto E. Degenerate Parabolic Equation. Universitext. 1993: New York; Springer DiBenedetto E. Partial Differential Equations. Cornerstones. 20102: Boston; Birkhäuser Boston Ltd. 1188.35001 DiBenedetto E, Friedman A. Hölder estimates for nonlinear degenerate parabolic systems. J. Reine Angew. Math. 1985; 357: 1-22. 783531. 0549.35061 Duzaar F, Grotowski JF. Optimal interior partial regularity for nonlinear elliptic systems: the method of A-harmonic approximation. Manuscripta Math. 2000; 103; 3: 267-298. 1802484. 0971.35025 Duzaar F, Grotowski JF, Kronz M. Partial and full boundary regularity for minimizers of functionals with nonquadratic growth. J. Convex Anal. 2004; 11; 2: 437-476. 2158914. 1066.49022 Eleuteri M, Marcellini P, Mascolo E. Lipschitz estimates for systems with ellipticity conditions at infinity. Ann. Mat. Pura Appl. 2016; 195: 1575-1603. 3537963. 1354.35035 Fonseca I, Fusco N, Marcellini P. An existence result for a nonconvex variational problem via regularity. ESAIM Control Optim. Calc. Var. 2002; 7: 69-95. 1925022. 1044.49011 Foss M. Global regularity for almost minimizers of nonconvex variational problems. Ann. Mat. Pura Appl. (4). 2008; 187; 2: 263-321. 2372803. 1223.49041 Foss M, Passarelli di Napoli A, Verde A. Global Morrey regularity results for asymptotically convex variational problems. Forum Math. 2008; 20: 921-953. 2445123. 1147.49029 Foss M, Passarelli di Napoli A, Verde A. Global Lipschitz regularity for almost minimizers of asymptotically convex variational problems. Ann. Mat. Pura Appl. (4). 2010; 189: 127-162. 2556762. 1180.49039 Giaquinta M, Modica G. Remarks on the regularity of the minimizers of certain degenerate functionals. Manuscripta Math. 1986; 57: 55-99. 866406. 0607.49003 Giusti E. Direct Methods in the Calculus of Variations. 2003: Tuck Link; World Scientific Publishing Company. 1028.49001 Leone C, Passarelli di Napoli A, Verde A. Lipschitz regularity for some asymptotically subquadratic problems. Nonlinear Anal. 2007; 67: 1532-1539. 2323300. 1113.49042 Liao, N: On the logarithmic type boundary modulus of continuity for the Stefan problem. (preprint) (2021) Manfredi JJ. Regularity for minima of functionals with p-growth. J. Differ. Equ. 1988; 76; 2: 203-212. 969420. 0674.35008 Marcellini P. Everywhere regularity for a class of elliptic systems without growth conditions. Ann. Scuola Norm. Sup. Pisa Cl. Sci.(4). 1996; 23; 1: 1-25. 1401415. 0922.35031 Passarelli di Napoli A, Verde A. A regularity result for asymptotically convex problems with lower order terms. J. Convex Anal. 2008; 15: 131-148. 2389007. 1245.35009 Raymond JP. Lipschitz regularity for some asymptotically convex problems. Proc. R. Soc. Edinburgh Sect. A. 1991; 117: 59-73. 1096219. 0725.49012 Santambrogio F, Vespri V. Continuity in two dimensions for a very degenerate elliptic equation. Nonlinear Anal. 2010; 73: 3832-3841. 2728558. 1202.35107 Scheven C, Schmidt T. Asymptotically regular problems. I. Higher integrability. J. Differ. Equ. 2010; 248: 745-791. 2578447. 1190.49044 Sverák V, Yan X. Non-Lipschitz minimizers of smooth uniformly convex variational integrals. Proc. Natl. Acad. Sci. USA. 2002; 99; 24: 15269-15276. 1946762. 1106.49046 Tolksdorf P. Everywhere-regularity for some quasilinear systems with a lack of ellipticity. Ann. Mat. Pura Appl. 1983; 134; 4: 241-266. 736742. 0538.35034 Uhlenbeck K. Regularity for a class of nonlinear elliptic systems. Acta Math. 1977; 138: 219-240. 474389. 0372.35030 Ural'ceva NN. Degenerate quasilinear elliptic systems. Zap. Naučn. Sem. Leningrad. Otdel. Mat. Inst. Steklov. (LOMI). 1968; 7: 184-222. 244628 Zatorska-Goldstein A. Very weak solutions of nonlinear subelliptic equations. Ann. Acad. Sci. Fenn. Math. 2005; 30; 2: 407-436. 2173373. 1082.35063 If $n=2$ we may choose any exponent $>2$ instead of $\frac{2n}{n-2}$ .

By Verena Bögelein; Frank Duzaar; Raffaella Giova and Antonia Passarelli di Napoli

Reported by Author; Author; Author; Author