Dynamic Optimal Transport with Optimal Preferential Paths

Resumo Técnico: Transporte Ótimo Dinâmico com Caminhos Preferenciais Ótimos

Formulação do Problema
Este trabalho investiga um problema de transporte ótimo dinâmico definido em um domínio $\Omega \subset \mathbb{R}^d$ composto por uma região de volume (bulk) e uma curva regular, não interceptante e de dimensão inferior $\Gamma \subset \Omega$. O modelo descreve o transporte de massa onde a população pode se mover dentro do volume $\Omega$ e ao longo da curva $\Gamma$, com a possibilidade de troca de massa entre as duas regiões a um custo. A massa total do sistema é conservada.

A dinâmica é governada por um par de equações de continuidade acopladas. Sejam $\rho_t$ e $\mu_t$ as distribuições de massa dependentes do tempo no volume e na curva, respectivamente, com seus respectivos momentos $J_t$ e $V_t$. A troca de massa é modelada por um termo de fluxo $f_t$ (troca normal) e termos de contorno $G^0_t, G^1_t$ (troca tangencial nas extremidades da curva). O sistema está sujeito a funções de mobilidade não lineares $m_\Omega$ e $m_\Gamma$, que regulam a dinâmica no volume e na curva, respectivamente.

O objetivo é minimizar um funcional de ação $A_\Gamma$ que agrega a energia cinética no volume, a energia cinética na curva e o custo de transferência de massa entre eles. O funcional é definido como:
$$A_\Gamma = \int_\Omega \frac{|J_t|^2}{m_\Omega(\rho_t)} d\lambda_\Omega + \alpha_1 \int_\Gamma \frac{|V_t|^2}{m_\Gamma(\mu_t)} d\lambda_\Gamma + \alpha_2 \int_\Gamma \frac{|f_t|^2}{m_\Gamma(\mu_t)} d\lambda_\Gamma + \dots$$
onde $\alpha_1, \alpha_2 > 0$ são parâmetros que controlam os custos relativos de transporte na curva versus no volume, e o custo de entrar ou sair da curva.

O artigo aborda dois cenários distintos:

1. Curva Fixa: A geometria de $\Gamma$ é fixa, e o objetivo é encontrar a evolução ótima de $(\rho_t, \mu_t, J_t, V_t, f_t)$.
2. Curva Variável: A forma de $\Gamma$ também é uma variável no problema de otimização. O objetivo é encontrar tanto a curva ótima quanto a dinâmica de transporte ótima.

Metodologia
Os autores empregam o método direto do cálculo de variações para estabelecer a existência de minimizadores.

• Existência para Curvas Fixas: Os autores definem um espaço de tuplas admissíveis que satisfazem as equações de continuidade acopladas em um sentido fraco (distribucional). Eles provam que o funcional de ação é convexo e semicontínuo inferiormente em relação à convergência tênue das medidas. Ao estabelecer estimativas de Hölder para as medidas dependentes do tempo e utilizar um teorema de Arzelà-Ascoli generalizado, eles demonstram a compacidade das sequências minimizantes, provando assim a existência de um minimizador.
• Análise Assintótica: O artigo analisa o comportamento do sistema conforme os parâmetros de custo $\alpha_2$ e $\alpha_3$ (que controlam a troca de massa) tendem a zero ou ao infinito.
  • Quando $\alpha \to \infty$, o custo de troca torna-se proibitivo, e o problema se desacopla em dois problemas de transporte ótimo independentes em $\Omega$ e $\Gamma$ (desde que as condições de compatibilidade de massa sejam atendidas).
  • Quando $\alpha \to 0$, a troca torna-se gratuita, e o problema converge para um problema de transporte acoplado com custo de troca zero.
• Existência para Curvas Variáveis: Para lidar com a otimização sobre a forma da curva, os autores introduzem um termo de regularização ao funcional de ação para evitar autointerseções e garantir regularidade. Eles utilizam a Energia de Ponto-Tangente ($E_{tp}^p$), originalmente definida para curvas fechadas, adaptada aqui para curvas abertas. Esta energia penaliza autointerseções e garante que as curvas com energia finita sejam subvariedades de dimensão um imersas. O funcional de regularização total $R(\gamma)$ inclui a Energia de Ponto-Tangente, um termo logarítmico para evitar que a curva se feche (garantindo condições de contorno bem definidas) e uma penalidade de comprimento.
  • A compacidade é estabelecida mostrando que a Energia de Ponto-Tangente limita a norma Sobolev-Slobodeckij $W^{2-1/p, p}$ da curva, que se imerge compactamente em $C^{1,\alpha}$.
  • A existência de um minimizador para o problema conjunto (curva e medidas) é provada combinando a compacidade das medidas (Proposição 4.6) com a semicontinuidade inferior da ação e do funcional de regularização.

Esquema Numérico
Os autores desenvolvem um algoritmo numérico baseado em uma abordagem primal-dual e uma formulação de Lagrange Aumentada (ALG2).

• Discretização: O domínio espaço-tempo é discretizado usando um método de elementos finitos com elementos lineares por partes para a curva e elementos de prisma para o espaço-tempo.
• Algoritmo: Para uma curva fixa, o algoritmo resolve um problema de ponto de sela alternando as atualizações das variáveis duais (potenciais) e das variáveis primais (densidades e momentos).
• Curva Variável: Para otimizar a curva, o método de Lagrange Aumentada é incorporado dentro de um loop de descida de gradiente projetado. O gradiente do funcional objetivo em relação aos pontos de controle da curva é aproximado usando diferenças centrais. A Energia de Ponto-Tangente é discretizada usando uma regra de quadratura específica para lidar com singularidades.

Principais Contribuições e Resultados

1. Formulação Rigorosa: O artigo fornece uma definição rigorosa de equações de continuidade acopladas para o transporte de massa entre um domínio de volume e uma curva de dimensão inferior, incluindo mecanismos de troca de massa e mobilidades não lineares.
2. Existência de Minimizadores: Estabelece a existência de minimizadores para o problema de transporte dinâmico tanto para caminhos preferenciais fixos quanto variáveis.
3. Limites Assintóticos: O estudo caracteriza o comportamento limite do sistema acoplado conforme os custos de troca desaparecem ou tornam-se infinitos, ligando o modelo ao transporte ótimo clásso não acoplado com mobilidades.
4. Regularização via Energia de Ponto-Tangente: O trabalho estende a aplicação da Energia de Ponto-Tangente para curvas abertas dentro de um arcabouço de transporte ótimo, provando que esta energia impõe a regularidade necessária (não autointerseção e regularidade $C^1$) para a existência de soluções no cenário de curva variável.
5. Implementação Numérica: Os autores apresentam um esquema numérico funcional capaz de computar geodésicas tanto para curvas fixas quanto otimizadas, demonstrando a influência dos parâmetros de custo e da geometria na dinâmica de transporte.

Significância
O artigo afirma contribuir para a teoria do transporte ótimo dinâmico ao generalizar o framework de Benamou-Brenier para incluir estruturas de dimensão inferior (curvas) dentro do domínio. Isso modela cenários onde o transporte é mais rápido ou mais barato ao longo de caminhos específicos (ex: rodovias), mas requer um custo para acessá-los. Ao permitir que o próprio caminho seja otimizado, o trabalho aborda problemas de planejamento ótimo onde a infraestrutura (o caminho) é determinada pelas distribuições iniciais e finais da população. O uso de mobilidades não lineares e o tratamento rigoroso da geometria da curva via Energia de Ponto-Tangente distinguem este trabalho de estudos anteriores que frequentemente dependiam de dualidade convexa ou focavam apenas em fronteiras fixas. Os resultados numéricos ilustram como a interação entre custos de transporte e a geometria do caminho dita a evolução da massa, oferecendo uma ferramenta para modelar redes de transporte complexas.