Discrete variational calculus for double-bracket dissipation

1. Formulação do Problema

O artigo aborda o desafio de simular numericamente sistemas mecânicos sujeitos a dissipação de duplo parêntese. Trata-se de sistemas dissipativos onde a energia diminui ao longo da trajetória, mas estruturas geométricas específicas — especificamente as órbitas coadjuntas da álgebra de Lie subjacente — devem permanecer invariantes.

• Contexto: Integradores geométricos padrão (como integradores simpléticos) excel em preservar propriedades geométricas (simpleticidade, mapas de momento) para sistemas conservativos. No entanto, adicionar dissipação geralmente destrói esses invariantes.
• Desafio Específico: Em sistemas como satélites com amortecedores, fluidos geofísicos e física de plasmas, a dinâmica evolui no dual de uma álgebra de Lie ($\mathfrak{g}^*$). Embora a energia seja dissipada, o sistema deve permanecer estritamente em uma órbita coadjunta específica (um conjunto de nível de invariantes de Casimir).
• Limitação dos Métodos Existentes: Integradores de propósito geral (por exemplo, Runge-Kutta) frequentemente falham em preservar a órbita coadjunta exatamente, levando a um desvio nos invariantes de Casimir e a um comportamento de longo prazo qualitativamente incorreto, mesmo que sejam precisos de alta ordem.

2. Metodologia

Os autores desenvolvem um Integrador Variacional Discreto especificamente adaptado para sistemas Euler-Poincaré e Lie-Poisson forçados com dissipação de duplo parêntese.

A. Fundamentação Teórica

1. Dissipação de Duplo Parêntese: O artigo utiliza o quadro estabelecido por Bloch et al. (1996), onde a força dissipativa é definida por meio de um mapa $\phi: \mathfrak{g} \to \mathfrak{g}$. As equações contínuas de movimento são:
   $$\dot{\mu} = \text{ad}^*_{\frac{\delta H}{\delta \mu}} \mu + \text{ad}^*_{\phi(\frac{\delta H}{\delta \mu})} \mu$$
   Esta formulação garante que, enquanto o Hamiltoniano $H$ diminui ($\dot{H} \leq 0$), as funções de Casimir $C$ são conservadas ($\dot{C} = 0$), mantendo a trajetória na órbita coadjunta.

2. Princípio Variacional Discreto:

   • Os autores definem um Lagrangiano Discreto $L_d: G \to \mathbb{R}$ no grupo de Lie $G$.
   • Eles introduzem um princípio variacional generalizado discreto com variações específicas $\delta w_k$ que incorporam o mapa dissipativo $\phi_d$.
   • A ação discreta é minimizada sujeita a variações da forma:
     $$\delta w_k = -\eta_k w_k + w_k (\text{Ad}_{\phi_d(w_k)} \eta_{k+1})$$
     onde $\eta_k$ são variações na álgebra de Lie que se anulam nas fronteiras.

3. Mapas de Retração:
   Para construir algoritmos práticos, os autores utilizam mapas de retração $\tau: \mathfrak{g} \to G$ (difeomorfismos locais analíticos que mapeiam a álgebra de Lie para o grupo de Lie, por exemplo, o mapa Exponencial ou o mapa de Cayley).

   • O momento discreto $M_k$ é atualizado via:
     $$M_{k+1} = \text{Ad}^*_{w_k \phi_d(w_k)} M_k$$
   • Esta regra de atualização é derivada discretizando o princípio variacional contínuo usando o mapa de retração para relacionar incrementos do grupo ($w_k$) com velocidades na álgebra de Lie ($\xi_k$).

B. Construção do Algoritmo

O método Duplo Parêntese Discreto (DDB) proposto procede da seguinte forma:

1. Inicializar o momento discreto $M_0 = \mu(t_0)$.
2. Para cada passo de tempo $k$:
   • Resolver para o elemento da álgebra de Lie $\xi_k$ tal que $M_k = (d\tau^{-1}_{h\xi_k})^* \frac{\partial l}{\partial \xi}(\xi_k)$.
   • Calcular o incremento do grupo $w_k = \tau(h\xi_k)$.
   • Atualizar o momento: $M_{k+1} = \text{Ad}^*_{w_k \phi_d(w_k)} M_k$.
3. O mapa $\phi_d$ é construído a partir do $\phi$ contínuo usando a retração: $\phi_d(w_k) = \tau(h\phi(\xi_k))$.

3. Contribuições Principais

1. Preservação Exata das Órbitas Coadjuntas: A contribuição primária é um integrador geométrico que preserva exatamente as órbitas coadjuntas na presença de dissipação de duplo parêntese. Diferentemente dos métodos padrão onde os invariantes de Casimir desviam, este método garante que a solução numérica permaneça estritamente na variedade definida pelas condições iniciais.
2. Derivação Variacional: O método é derivado de um princípio variacional discreto, garantindo que herde a estabilidade e as propriedades de preservação de estrutura dos integradores variacionais.
3. Flexibilidade via Mapas de Retração: O quadro permite o uso de diferentes mapas de retração (Exponencial, Cayley, etc.), possibilitando a construção de integradores com diferentes custos computacionais e ordens de precisão (por exemplo, esquemas simétricos de segunda ordem).
4. Aplicação ao Corpo Rígido: A teoria é aplicada explicitamente ao corpo rígido dissipativo, fornecendo uma formulação matricial concreta para $SO(3)$.

4. Resultados e Simulações Numéricas

Os autores validam o método usando o corpo rígido dissipativo como caso de teste, comparando o método DDB (usando tanto retrações Exponenciais quanto de Cayley) contra:

• Runge-Kutta de 4ª Ordem (RK4)
• Lobatto IIIC de 3 estágios
• DOP853 (8ª ordem, usado como referência)

Principais Descobertas:

• Preservação de Casimir: Em todas as simulações, os métodos DDB preservaram o invariante de Casimir (a magnitude do momento angular $\|M\|$) até a precisão da máquina. Em contraste, RK4 e Lobatto IIIC mostraram desvio significativo no Casimir, fazendo com que a trajetória deixasse a esfera coadjunta.
• Comportamento Assintótico:
  • Cenário A (Equilíbrios Isolados): Para um corpo rígido com momentos de inércia distintos, o sistema converge para a rotação em torno do eixo de máxima inércia (energia mínima). Os métodos DDB convergiram para o ponto limite correto, enquanto métodos não geométricos desviaram devido a erros de Casimir.
  • Cenário B (Equilíbrios de Círculo Máximo): Para um corpo simétrico (dois momentos de inércia iguais), o sistema converge para um círculo máximo de equilíbrios. Os métodos DDB convergiram corretamente para um ponto neste círculo, enquanto outros métodos desviaram da variedade.
• Eficiência: Os métodos DDB (especificamente com retração de Cayley) foram computacionalmente eficientes, frequentemente mais rápidos que o método Lobatto IIIC de alta ordem, enquanto forneciam precisão qualitativa superior.
• Convergência: Os métodos demonstraram convergência para a solução de referência, com o erro no momento final sendo significativamente menor para os métodos DDB em simulações de longo tempo devido à preservação da geometria subjacente.

5. Significado

Este trabalho preenche uma lacuna crítica na integração numérica geométrica:

• Fidelidade Física: Para sistemas como satélites com amortecedores internos ou confinamento de plasma, a conservação do momento angular (ou o Casimir equivalente) é fisicamente obrigatória. Integradores padrão violam isso, levando a resultados não físicos ao longo de tempos de simulação longos.
• Robustez: Os integradores propostos são robustos para integração de longo prazo, capturando corretamente a estabilidade assintótica do sistema sem desvio artificial.
• Generalizabilidade: O quadro não se limita a corpos rígidos. Os autores notam que esta abordagem pode ser estendida para álgebroides de Lie e fluidos ideais, oferecendo um caminho para simular dinâmica de fluidos complexa com dissipação enquanto preserva os invariantes de Casimir de dimensão infinita das equações de Euler.

Em conclusão, o artigo apresenta um integrador geométrico rigoroso, baseado em variacionais, que lida com sucesso com a dissipação de duplo parêntese, garantindo que simulações numéricas de sistemas mecânicos dissipativos permaneçam geometricamente fiéis à física contínua.