Harnack inequality for anisotropic fully nonlinear equations with nonstandard growth

1. Problema Investigado

O artigo aborda a regularidade de soluções de viscosidade para uma classe de equações elípticas não lineares completas (fully nonlinear) degeneradas, apresentando crescimento não padrão e anisotropia. A equação principal estudada é dada na forma não divergente:

$$\begin{cases} M^{-}_{\lambda, \Lambda} \left( \left[ |D_i u|^{\frac{p_i-2}{2}} \right] D^2 u \left[ |D_i u|^{\frac{p_i-2}{2}} \right] \right) - \mu \sum_{i} |D_i u|^{p_i-1} \leq c_0 \\ M^{+}_{\lambda, \Lambda} \left( \left[ |D_i u|^{\frac{p_i-2}{2}} \right] D^2 u \left[ |D_i u|^{\frac{p_i-2}{2}} \right] \right) + \mu \sum_{i} |D_i u|^{p_i-1} \geq -c_0 \end{cases}$$

Características principais do problema:

• Anisotropia Degenerada: O operador degenera se qualquer componente individual do gradiente ($\partial_i u$) for zero, não apenas quando o gradiente total é zero (como no Laplaciano $p$-Laplaciano padrão). Isso torna a análise mais complexa do que no caso isotrópico.
• Crescimento Não Padrão: O crescimento do operador varia entre diferentes potências $p_1, \dots, p_n$ (com $2 \leq p_1 \leq \dots \leq p_n$), originando-se de funcionais do tipo $\int F(Dw) dx$ onde $|z|^{p_1} \lesssim F(z) \lesssim |z|^{p_1} + |z|^{p_n}$.
• Soluções de Viscosidade: Diferente de soluções fracas (onde se usa a teoria de De Giorgi-Nash-Moser e desigualdades de Sobolev), as soluções de viscosidade não possuem essas ferramentas diretas, exigindo métodos puramente geométricos e de comparação.

O objetivo central é estabelecer uma Desigualdade de Harnack intrínseca e, consequentemente, a regularidade de Hölder local para essas soluções, sob condições adequadas sobre os expoentes $p_i$.

---

2. Metodologia

Os autores desenvolvem uma abordagem baseada em geometria intrínseca, adaptando técnicas de escalonamento não homogêneo para lidar com a degenerescência anisotrópica. A prova é estruturada em três pilares principais:

A. Escalonamento Intrínseco e Cubos Intrínsecos

Devido à não homogeneidade da equação, o domínio natural de análise não é o cubo euclidiano padrão $Q_r$, mas sim um cubo intrínseco $K_r(M)$, definido como:
$$K_r(M) := \{ |x_i| \leq r^{\alpha_i} M^{\beta_i} \}$$
onde $\alpha_i = 1/p_i$ e $\beta_i = -(p_n - p_i)/p_i$. Este formato ajusta a geometria do domínio à escala natural da solução no nível $M$.

B. Estimativas de Medida (Método do Parabolóide Deslizante)

Para obter estimativas básicas de medida (Lema 3.1), os autores adaptam o método do parabolóide deslizante (introduzido por Cabré e desenvolvido por Savin).

• Adaptação Anisotrópica: Em vez de um parabolóide quadrático padrão, eles utilizam um parabolóide anisotrópico da forma $\phi(x) = \sum K_i |x_i|^{p_i/(p_i-1)}$.
• Desafio de Regularidade: Como essa função não é $C^2$ nos hiperplanos coordenados, eles utilizam um argumento de "restrição a fatias" (restriction to slices) para lidar com a singularidade, assumindo $C^2$ apenas para clareza na demonstração, mas garantindo a validade no sentido de viscosidade.

C. Propriedade de Dobramento (Doubling Property) e Função Barreira

A contribuição mais técnica e inovadora do trabalho é a construção de uma função barreira anisotrópica explícita (Seção 4) para provar a propriedade de dobramento (Lema 5.1).

• A Barreira: Eles constroem uma função da forma $\Phi(x) = (|b_1 x_1|^{a_1} + \dots + |b_n x_n|^{a_n})^{-\gamma}$.
• Condição Estrutural: A existência de parâmetros $a_i, \gamma$ que tornam $\Phi$ uma subsolução válida depende criticamente de uma condição específica sobre os expoentes $p_i$ (Equação 1.6).
• Mecanismo: Ao aplicar o princípio da comparação entre a solução $u$ e a barreira $\Phi$, eles demonstram que se a solução é grande em um ponto, ela deve ser grande em uma vizinhança intrínseca, estabelecendo a propriedade de dobramento.

D. Estimativa $L^\epsilon$ e Iteração

Combinando as estimativas de medida e a propriedade de dobramento, os autores derivam uma estimativa de decaimento algébrico para os conjuntos de nível (Teorema 6.1), utilizando o Lema de Cobertura de Vitali adaptado para cubos intrínsecos. Isso permite provar a desigualdade de Harnack final.

---

3. Contribuições Chave

1. Primeira Desigualdade de Harnack para Soluções de Viscosidade: O artigo estabelece a primeira desigualdade de Harnack para soluções de viscosidade de operadores anisotrópicos degenerados com crescimento não padrão. Resultados anteriores existiam apenas para soluções fracas ou para casos isotrópicos.
2. Construção de Barreira Anisotrópica Explícita: A construção da função barreira $\Phi$ é uma novidade técnica significativa. Ela supera a dificuldade de não ter uma estrutura variacional (como no caso de soluções fracas) e a ausência de desigualdades de Sobolev adequadas para soluções de viscosidade neste contexto.
3. Condição Estrutural Otimizada (Relativa): Os autores identificam uma condição específica sobre a diferença entre os expoentes $p_1$ e $p_n$ (Equação 1.6) necessária para a validade da desigualdade. Eles discutem como essa condição se relaciona com condições conhecidas para limitabilidade e regularidade de soluções fracas.
4. Generalização do Método de Geometria Intrínseca: O trabalho estende a técnica de geometria intrínseca (originalmente para o Laplaciano $p$-parabólico) para o caso elíptico completamente não linear com anisotropia extrema, onde cada direção tem uma taxa de crescimento diferente.

---

4. Resultados Principais

O resultado central é o Teorema 1.1, que afirma:

Seja $u \in C(\Omega)$ uma solução de viscosidade não negativa satisfazendo (1.1). Sob a condição de que os expoentes $p_i$ satisfaçam:
$$\frac{p_n - 1}{p_1 - 1} \leq \frac{p_n}{p} \frac{n}{n - \lambda/\Lambda}$$
(onde $p$ é a média harmônica dos $p_i$), existem constantes universais $C_0, R_0$ tais que:

1. Limite Superior: Se $u(0) > 0$, então $\sup_{K_r(u(0))} u \leq C_0 u(0)$.
2. Limite Inferior: Se $u(0) > 0$, então $\inf_{K_r(C_0 u(0))} u \geq \frac{1}{C_0} u(0)$.

Corolário de Regularidade (Corolário 1.1):
Como consequência direta da desigualdade de Harnack, prova-se que a solução $u$ é localmente Hölder contínua em $\Omega$ com um expoente universal $\alpha \in (0, 1)$.

---

5. Significado e Impacto

• Preenchimento de Lacunas Teóricas: Antes deste trabalho, a regularidade de Hölder e a desigualdade de Harnack para soluções de viscosidade de operadores anisotrópicos degenerados com crescimento não padrão eram problemas abertos. As técnicas clássicas de De Giorgi falhavam devido à falta de estrutura variacional e desigualdades de Sobolev adequadas para este contexto.
• Robustez Metodológica: A demonstração de que é possível construir barreiras explícitas e adaptar o método de parabolóide deslizante para geometrias altamente anisotrópicas abre caminho para o estudo de outras classes de equações não lineares degeneradas.
• Aplicabilidade: O resultado é relevante para modelos físicos e matemáticos onde a difusão ou condutividade varia drasticamente dependendo da direção e da magnitude do gradiente, como em materiais compósitos anisotrópicos ou fluidos não newtonianos complexos.
• Limitações e Futuro: Os autores notam que a condição sobre o "gap" entre $p_1$ e $p_n$ é mais restritiva do que a condição ótima conhecida para a limitabilidade de soluções fracas, sugerindo que a construção da barreira pode ser refinada no futuro para relaxar essa condição.

Em resumo, o artigo representa um avanço significativo na teoria de regularidade de EDPs não lineares, unindo geometria intrínseca, análise de viscosidade e construção de barreiras para resolver um problema de alta complexidade técnica.