Optimal (partial) transport to non-convex polygonal domains

Resumo Técnico: Transporte Ótimo (Parcial) para Domínios Poligonais Não Convexos

Formulação do Problema
Este artigo investiga a regularidade de mapas de transporte ótimo e de fronteiras livres quando o domínio alvo é um polígono não convexo em $\mathbb{R}^2$. O estudo aborda dois cenários distintos:

1. Transporte Ótimo Completo: Dados dois domínios limitados $\Omega$ e $\Omega^*$ com medida igual, o objetivo é encontrar um plano de transporte que minimize o custo quadrático $\int |x-y|^2 d\gamma$. Pelo teorema de Brenier, isso é equivalente a encontrar um potencial convexo $u$ tal que $(Du)_\# \chi_\Omega = \chi_{\Omega^*}$. O foco está na estrutura do conjunto singular $A = \{x \in \Omega : u \text{ não é diferenciável em } x\}$.
2. Transporte Ótimo Parcial: Dada uma massa $m < \min(|\Omega|, |\Omega^*|)$, o problema busca transportar um subdomínio $U \subset \Omega$ para um subdomínio $V \subset \Omega^*$. A fronteira da região de transporte ativa, $F = \partial U \cap \Omega$, é denominada "fronteira livre". O artigo analisa a regularidade de $F$ quando $\Omega^*$ é um polígono não convexo.

Metodologia
Os autores empregam uma combinação de análise convexa, teoria da medida geométrica e teoria da regularidade para equações de Monge-Ampère. Os componentes metodológicos-chave incluem:

• Transformada de Legendre e Potenciais Duais: A análise depende fortemente da função convexa dual $v = u^*$ (ou $v$ no caso de transporte parcial). A diferenciabilidade de $u$ em $x$ corresponde à diferenciabilidade de $v$ em pontos do subgradiente $\partial^- u(x)$.
• Localização e Conjuntos de Nível Inferior: Os autores utilizam conjuntos de nível inferior centrados $S^c_h[v](y)$ para localizar o comportamento do potencial próximo a singularidades. Ao estabelecer que esses conjuntos de nível inferior permanecem dentro de pequenas vizinhanças da fronteira alvo, eles invocam estimativas quantitativas de convexidade uniforme.
• Decomposição Geométrica de Singularidades: O conjunto singular (ou fronteira livre) é decomposto com base na geometria do subgradiente $\partial^- u(x)$ em relação ao polígono alvo $\Omega^*$. Os autores categorizam pontos $x$ com base se $\partial^- u(x)$ contém um vértice de $\Omega^*$, toca uma aresta ou toca múltiplas arestas.
• Escalonamento de Regularidade: Para regiões onde a geometria alvo é localmente convexa (por exemplo, próximo ao interior de uma aresta), os autores aplicam a teoria da regularidade de Caffarelli e estimativas de Schauder para equações uniformemente elípticas. Eles demonstram que o potencial é $C^\infty$ nessas regiões, implicando que o conjunto singular é uma curva suave.
• Estimativas de Obliquidade: No caso de transporte parcial, a prova de suavidade depende do estabelecimento de uma estimativa de obliquidade (o ângulo entre a normal à fronteira livre e a normal à fronteira alvo é não nulo), o que permite a aplicação de resultados de regularidade de fronteira.

Principais Contribuições e Resultados

• Estrutura do Conjunto Singular (Transporte Completo):
  O artigo prova o Teorema 1.1, que fornece uma descrição precisa do conjunto singular $A$ quando $\Omega^*$ é um polígono não convexo. Os autores mostram que $A$ é ou um conjunto finito ou, afastado de um número finito de pontos $\{x_i\}$, é localmente uma curva unidimensional suave.

  • Os pontos "excepcionais" $\{x_i\}$ correspondem a locais onde o subgradiente $\partial^- u(x)$ toca um vértice do polígono ou toca três ou mais arestas.
  • Perto desses pontos excepcionais, o conjunto singular é uma união de um número finito de curvas Lipschitz.
  • A prova envolve decompor $A$ em conjuntos $\Sigma_1$ (tocando vértices), $\Sigma'_2$ (tocando exatamente duas arestas) e $\Sigma''_2$ (tocando $\ge 3$ arestas). Mostra-se que $\Sigma_1$ e $\Sigma''_2$ são conjuntos finitos, enquanto $\Sigma'_2$ consiste em curvas suaves.

• Regularidade da Fronteira Livre (Transporte Parcial):
  O artigo prova o Teorema 1.2, estabelecendo que a fronteira livre $F$ no problema de transporte ótimo parcial é suave, afastada de um número finito de pontos singulares.

  • Similar ao caso de transporte completo, os pontos singulares na fronteira livre surgem onde o subgradiente toca um vértice de $\Omega^*$ ou múltiplas arestas.
  • Os autores provam que o conjunto de pontos tocando vértices ($F_1$) ou múltiplas arestas ($F_2$) é finito.
  • Para pontos tocando exatamente uma aresta ($F_3$), a fronteira livre é mostrada como sendo suave ($C^\infty$) via estimativas de obliquidade e teoria de regularidade de fronteira.

• Generalizações e Conjecturas:

  • A Observação 1.2 nota que esses resultados se estendem a densidades positivas suaves gerais e, com regularidade ligeiramente mais fraca ($C^{1,\alpha}$), a densidades limitadas superior e inferiormente.
  • A Observação 3.2 indica que o método se aplica a alvos que são uniões de dois domínios convexos limitados em dimensões gerais.
  • A Seção 5 formula duas conjecturas para dimensões superiores ($\mathbb{R}^n$). A Conjectura 5.1 postula que, para transporte completo para um politopo não convexo, o conjunto singular é localmente uma hipersuperfície suave de dimensão $(n-1)$, afastado de um conjunto de dimensão de Hausdorff no máximo $n-2$. A Conjectura 5.2 faz uma afirmação similar para a fronteira livre no transporte parcial.

Significado e Alegações
O artigo alega fornecer uma descrição mais precisa da estrutura de singularidade no transporte ótimo do que anteriormente estabelecida para domínios planares não convexos. Enquanto trabalhos anteriores (por exemplo, Figalli [16]) estabeleceram que conjuntos singulares são variedades $C^1$ afastados de um conjunto enumerável, este trabalho refina esse resultado para alvos poligonais, mostrando que o conjunto singular é suave afastado de um conjunto finito de pontos.

Os autores destacam a relevância prática desses resultados em computações numéricas, especificamente para geração de malhas em domínios poligonais usando algoritmos de transporte ótimo. Ao caracterizar as singularidades, o artigo auxilia na compreensão das limitações e do comportamento de tais algoritmos próximos a cantos não convexos.

O artigo mantém um tom modesto em relação a dimensões superiores, afirmando explicitamente que os resultados para $n=2$ são provados, enquanto a extensão para $n \ge 3$ permanece uma conjectura. Os autores não alegam ter resolvido o caso de dimensões superiores, mas fornecem uma estrutura estrutural e conjecturas específicas para investigação futura.