Optimal (partial) transport to non-convex polygonal domains

Technische Samenvatting: Optimale (Partiële) Transport naar Niet-Convexe Polygonaal Domeinen

Probleemformulering
Dit artikel onderzoekt de regulariteit van optimale transportafbeeldingen en vrije randen wanneer het doeldomein een niet-convex polygoon is in $\mathbb{R}^2$. De studie behandelt twee verschillende scenario's:

1. Volledige Optimaal Transport: Gegeven twee begrensde domeinen $\Omega$ en $\Omega^*$ met gelijke maat, is het doel een transportplan te vinden dat de kwadratische kosten $\int |x-y|^2 d\gamma$ minimaliseert. Volgens de stelling van Brenier is dit equivalent aan het vinden van een convex potentiaal $u$ zodat $(Du)_\# \chi_\Omega = \chi_{\Omega^*}$. De focus ligt op de structuur van de singuliere verzameling $A = \{x \in \Omega : u \text{ is niet differentieerbaar in } x\}$.
2. Optimaal Partiële Transport: Gegeven een massa $m < \min(|\Omega|, |\Omega^*|)$, zoekt het probleem naar het transporteren van een subdomein $U \subset \Omega$ naar een subdomein $V \subset \Omega^*$. De rand van het actieve transportgebied, $F = \partial U \cap \Omega$, wordt de "vrije rand" genoemd. Het artikel analyseert de regulariteit van $F$ wanneer $\Omega^*$ een niet-convex polygoon is.

Methodologie
De auteurs maken gebruik van een combinatie van convexe analyse, meettheorie van maat en regulariteitstheorie voor Monge-Ampère-vergelijkingen. Belangrijke methodologische componenten zijn:

• Legendre-transformatie en Duale Potentialen: De analyse steunt zwaar op de duale convexe functie $v = u^*$ (of $v$ in het geval van partiële transport). De differentieerbaarheid van $u$ in $x$ komt overeen met de differentieerbaarheid van $v$ in punten van de subgradiënt $\partial^- u(x)$.
• Lokalisatie en Subniveausets: De auteurs maken gebruik van gecentreerde subniveausets $S^c_h[v](y)$ om het gedrag van de potentiaal in de buurt van singulariteiten te lokaliseren. Door vast te stellen dat deze subniveausets binnen kleine omgevingen van de doelrand blijven, maken ze gebruik van kwantitatieve uniforme convexiteitsschattingen.
• Geometrische Decompositie van Singulariteiten: De singuliere verzameling (of vrije rand) wordt gedecomponeerd op basis van de geometrie van de subgradiënt $\partial^- u(x)$ ten opzichte van het doelpolygoon $\Omega^*$. De auteurs categoriseren punten $x$ op basis van of $\partial^- u(x)$ een hoekpunt van $\Omega^*$ raakt, één rand raakt, of meerdere randen raakt.
• Regulariteit-Bootstrapping: Voor gebieden waar de doelgeometrie lokaal convex is (bijvoorbeeld in de buurt van het interieur van een rand), passen de auteurs de regulariteitstheorie van Caffarelli en Schauder-schattingen voor uniform elliptische vergelijkingen toe. Ze tonen aan dat de potentiaal $C^\infty$ is in deze gebieden, wat impliceert dat de singuliere verzameling een gladde kromme is.
• Schattingen voor Obliquiteit: In het geval van partiële transport rust het bewijs van gladheid op het vaststellen van een obliquiteitsschatting (de hoek tussen de normaal op de vrije rand en de normaal op de doelrand is niet nul), wat de toepassing van randregulariteitsresultaten mogelijk maakt.

Belangrijkste Bijdragen en Resultaten

• Structuur van de Singuliere Verzameling (Volledige Transport):
  Het artikel bewijst Stelling 1.1, die een nauwkeurige beschrijving geeft van de singuliere verzameling $A$ wanneer $\Omega^*$ een niet-convex polygoon is. De auteurs tonen aan dat $A$ óf een eindige verzameling is, óf, op afstand van een eindig aantal punten $\{x_i\}$, lokaal een gladde eendimensionale kromme is.

  • De "uitzonderlijke" punten $\{x_i\}$ komen overeen met locaties waar de subgradiënt $\partial^- u(x)$ een hoekpunt van het polygoon raakt of drie of meer randen raakt.
  • In de buurt van deze uitzonderlijke punten is de singuliere verzameling een vereniging van een eindig aantal Lipschitz-krommen.
  • Het bewijs houdt in dat $A$ wordt gedecomponeerd in verzamelingen $\Sigma_1$ (hoekpunten rakend), $\Sigma'_2$ (precies twee randen rakend) en $\Sigma''_2$ ($\ge 3$ randen rakend). Er wordt aangetoond dat $\Sigma_1$ en $\Sigma''_2$ eindige verzamelingen zijn, terwijl $\Sigma'_2$ bestaat uit gladde krommen.

• Regulariteit van de Vrije Rand (Partiële Transport):
  Het artikel bewijst Stelling 1.2, waarin wordt vastgesteld dat de vrije rand $F$ in het optimale partiële transportprobleem glad is, op afstand van een eindig aantal singuliere punten.

  • Net als in het geval van volledige transport ontstaan de singuliere punten op de vrije rand waar de subgradiënt een hoekpunt van $\Omega^*$ raakt of meerdere randen raakt.
  • De auteurs bewijzen dat de verzameling punten die hoekpunten ($F_1$) of meerdere randen ($F_2$) raken, eindig is.
  • Voor punten die precies één rand raken ($F_3$), wordt aangetoond dat de vrije rand glad is ($C^\infty$) via obliquiteitsschattingen en randregulariteitstheorie.

• Generalisaties en Conjecturen:

  • Opmerking 1.2 merkt op dat deze resultaten uitbreiden naar algemene positieve gladde dichtheden en, met iets zwakkere regulariteit ($C^{1,\alpha}$), naar dichtheden die boven en onder begrensd zijn.
  • Opmerking 3.2 geeft aan dat de methode van toepassing is op doelen die verenigingen zijn van twee begrensde convexe domeinen in algemene dimensies.
  • Sectie 5 formuleert twee conjecturen voor hogere dimensies ($\mathbb{R}^n$). Conjectuur 5.1 stelt dat voor volledige transport naar een niet-convex polytoop, de singuliere verzameling lokaal een gladde $(n-1)$-dimensionale hypervlak is, op afstand van een verzameling met Hausdorff-dimensie ten hoogste $n-2$. Conjectuur 5.2 doet een vergelijkbare bewering voor de vrije rand in partiële transport.

Betekenis en Claims
Het artikel claimt een nauwkeurigere beschrijving te geven van de singulariteitsstructuur in optimaal transport dan eerder vastgesteld voor niet-convexe planaire domeinen. Hoewel eerder werk (bijvoorbeeld Figalli [16]) vaststelde dat singuliere verzamelingen $C^1$-variëteiten zijn op afstand van een aftelbare verzameling, verfijnt dit werk dat resultaat voor polygonaal doelen, en toont het aan dat de singuliere verzameling glad is op afstand van een eindige verzameling punten.

De auteurs benadrukken de praktische relevantie van deze resultaten in numerieke berekeningen, specifiek voor mesh-generatie op polygonaal domeinen met behulp van optimal transport-algoritmen. Door de singulariteiten te karakteriseren, helpt het artikel bij het begrijpen van de beperkingen en het gedrag van dergelijke algoritmen in de buurt van niet-convexe hoeken.

Het artikel behoudt een bescheiden toon met betrekking tot hogere dimensies, en stelt expliciet dat de resultaten voor $n=2$ bewezen zijn, terwijl de uitbreiding naar $n \ge 3$ een conjectuur blijft. De auteurs claimen niet het hogere-dimensionale geval opgelost te hebben, maar bieden wel een structureel raamwerk en specifieke conjecturen voor toekomstig onderzoek.