Optimal (partial) transport to non-convex polygonal domains

Resumen Técnico: Transporte Óptimo (Parcial) a Dominios Poligonales No Convexos

Formulación del Problema
Este artículo investiga la regularidad de los mapas de transporte óptimo y las fronteras libres cuando el dominio objetivo es un polígono no convexo en $\mathbb{R}^2$. El estudio aborda dos escenarios distintos:

1. Transporte Óptimo Completo: Dados dos dominios acotados $\Omega$ y $\Omega^*$ con medida igual, el objetivo es encontrar un plan de transporte que minimice el costo cuadrático $\int |x-y|^2 d\gamma$. Por el teorema de Brenier, esto es equivalente a encontrar un potencial convexo $u$ tal que $(Du)_\# \chi_\Omega = \chi_{\Omega^*}$. El enfoque se centra en la estructura del conjunto singular $A = \{x \in \Omega : u \text{ no es diferenciable en } x\}$.
2. Transporte Óptimo Parcial: Dada una masa $m < \min(|\Omega|, |\Omega^*|)$, el problema busca transportar un subdominio $U \subset \Omega$ a un subdominio $V \subset \Omega^*$. La frontera de la región de transporte activa, $F = \partial U \cap \Omega$, se denomina "frontera libre". El artículo analiza la regularidad de $F$ cuando $\Omega^*$ es un polígono no convexo.

Metodología
Los autores emplean una combinación de análisis convexo, teoría de la medida geométrica y teoría de regularidad para ecuaciones de Monge-Ampère. Los componentes metodológicos clave incluyen:

• Transformada de Legendre y Potenciales Duales: El análisis depende en gran medida de la función convexa dual $v = u^*$ (o $v$ en el caso de transporte parcial). La diferenciabilidad de $u$ en $x$ corresponde a la diferenciabilidad de $v$ en puntos del subgradiente $\partial^- u(x)$.
• Localización y Conjuntos de Nivel Inferior: Los autores utilizan conjuntos de nivel inferior centrados $S^c_h[v](y)$ para localizar el comportamiento del potencial cerca de las singularidades. Al establecer que estos conjuntos de nivel inferior permanecen dentro de vecindades pequeñas de la frontera objetivo, invocan estimaciones cuantitativas de convexidad uniforme.
• Descomposición Geométrica de Singularidades: El conjunto singular (o la frontera libre) se descompone basándose en la geometría del subgradiente $\partial^- u(x)$ en relación con el polígono objetivo $\Omega^*$. Los autores categorizan los puntos $x$ según si $\partial^- u(x)$ contiene un vértice de $\Omega^*$, toca un borde o toca múltiples bordes.
• Mejora de la Regularidad (Bootstrapping): Para regiones donde la geometría objetivo es localmente convexa (por ejemplo, cerca del interior de un borde), los autores aplican la teoría de regularidad de Caffarelli y las estimaciones de Schauder para ecuaciones uniformemente elípticas. Demuestran que el potencial es $C^\infty$ en estas regiones, lo que implica que el conjunto singular es una curva suave.
• Estimaciones de Oblicuidad: En el caso de transporte parcial, la prueba de suavidad se basa en establecer una estimación de oblicuidad (el ángulo entre la normal a la frontera libre y la normal a la frontera objetivo es distinto de cero), lo que permite la aplicación de resultados de regularidad de frontera.

Contribuciones y Resultados Clave

• Estructura del Conjunto Singular (Transporte Completo):
  El artículo demuestra el Teorema 1.1, que proporciona una descripción precisa del conjunto singular $A$ cuando $\Omega^*$ es un polígono no convexo. Los autores muestran que $A$ es o bien un conjunto finito o, alejado de un número finito de puntos $\{x_i\}$, es localmente una curva unidimensional suave.

  • Los puntos "excepcionales" $\{x_i\}$ corresponden a ubicaciones donde el subgradiente $\partial^- u(x)$ toca un vértice del polígono o toca tres o más bordes.
  • Cerca de estos puntos excepcionales, el conjunto singular es una unión de un número finito de curvas Lipschitz.
  • La prueba implica descomponer $A$ en conjuntos $\Sigma_1$ (tocando vértices), $\Sigma'_2$ (tocando exactamente dos bordes) y $\Sigma''_2$ (tocando $\ge 3$ bordes). Se demuestra que $\Sigma_1$ y $\Sigma''_2$ son conjuntos finitos, mientras que $\Sigma'_2$ consiste en curvas suaves.

• Regularidad de la Frontera Libre (Transporte Parcial):
  El artículo demuestra el Teorema 1.2, estableciendo que la frontera libre $F$ en el problema de transporte óptimo parcial es suave alejada de un número finito de puntos singulares.

  • Similar al caso de transporte completo, los puntos singulares en la frontera libre surgen donde el subgradiente toca un vértice de $\Omega^*$ o múltiples bordes.
  • Los autores prueban que el conjunto de puntos tocando vértices ($F_1$) o múltiples bordes ($F_2$) es finito.
  • Para los puntos tocando exactamente un borde ($F_3$), se demuestra que la frontera libre es suave ($C^\infty$) mediante estimaciones de oblicuidad y teoría de regularidad de frontera.

• Generalizaciones y Conjeturas:

  • La Observación 1.2 señala que estos resultados se extienden a densidades positivas suaves generales y, con una regularidad ligeramente más débil ($C^{1,\alpha}$), a densidades acotadas superior e inferiormente.
  • La Observación 3.2 indica que el método se aplica a objetivos que son uniones de dos dominios convexos acotados en dimensiones generales.
  • La Sección 5 formula dos conjeturas para dimensiones superiores ($\mathbb{R}^n$). La Conjetura 5.1 postula que, para el transporte completo a un politopo no convexo, el conjunto singular es localmente una hipersuperficie suave de dimensión $(n-1)$ alejada de un conjunto de dimensión de Hausdorff a lo sumo $n-2$. La Conjetura 5.2 hace una afirmación similar para la frontera libre en el transporte parcial.

Significado y Afirmaciones
El artículo afirma proporcionar una descripción más precisa de la estructura de singularidades en el transporte óptimo que la establecida previamente para dominios planares no convexos. Mientras que trabajos anteriores (por ejemplo, Figalli [16]) establecieron que los conjuntos singulares son variedades $C^1$ alejados de un conjunto numerable, este trabajo refina ese resultado para objetivos poligonales, mostrando que el conjunto singular es suave alejado de un conjunto finito de puntos.

Los autores destacan la relevancia práctica de estos resultados en cálculos numéricos, específicamente para la generación de mallas en dominios poligonales utilizando algoritmos de transporte óptimo. Al caracterizar las singularidades, el artículo ayuda a comprender las limitaciones y el comportamiento de dichos algoritmos cerca de esquinas no convexas.

El artículo mantiene un tono modesto respecto a las dimensiones superiores, declarando explícitamente que los resultados para $n=2$ están demostrados, mientras que la extensión a $n \ge 3$ permanece como una conjetura. Los autores no afirman haber resuelto el caso de dimensiones superiores, sino que proporcionan un marco estructural y conjeturas específicas para futuras investigaciones.