On Runge-Kutta convolution quadrature based fractional variational integrators

Resumo Técnico: Integradores Variacionais Fracionários Baseados em Quadratura de Convolução de Runge-Kutta

Enunciado do Problema
Equações diferenciais fracionárias são essenciais para modelar sistemas com memória de longo prazo, como materiais viscoelásticos e difusão anômala. Embora o cálculo fracionário estenda os princípios variacionais para incluir sistemas dissipativos (amortecimento fracionário), a construção de integradores numéricos de alta ordem para esses sistemas permanece um desafio. Trabalhos anteriores por [23] combinaram a quadratura de convolução de diferenciação retroativa (BDFCQ) com métodos variacionais de Galerkin para a parte conservativa do Lagrangiano. No entanto, essa abordagem sofre de um efeito de "saturação": apesar do potencial de ordem de precisão superior da BDFCQ, os integradores variacionais fracionários (FVIs) resultantes são limitados à precisão de segunda ordem. Essa limitação surge porque a estrutura de passos múltiplos da BDFCQ opera apenas nos valores nodais principais, falhando em utilizar a informação dos estágios internos inerente às discretizações de Galerkin de ordem superior da ação conservativa. Consequentemente, há uma necessidade de integradores variacionais fracionários onde a discretização da parte fracionária seja estruturalmente compatível com a estrutura de estágios internos da discretização conservativa.

Metodologia
Os autores propõem uma nova classe de Integradores Variacionais Fracionários (FVIs) que combinam a Quadratura de Convolução de Runge-Kutta (RKCQ) com métodos variacionais de Galerkin de ordem superior. A metodologia baseia-se em três pilares:

1. Estrutura Variacional Aumentada: Seguindo o princípio de Hamilton restrito [22], o sistema é modelado pelo dobro das variáveis de estado ($x$ e $y$) para incorporar o amortecimento fracionário. O funcional da ação é dividido em uma parte conservativa ($L_{con}$) e uma parte fracionária ($L_{fra}$).
2. Estratégia de Discretização:
   • Parte Conservativa: Aproximada usando métodos variacionais de Galerkin de ordem superior com interpolação polinomial de grau $s$ e quadratura de Lobatto.
   • Parte Fracionária: Aproximada usando RKCQ baseada em métodos Lobatto IIIC. Diferente da BDFCQ, a RKCQ avalia a integral de convolução nos estágios internos de Runge-Kutta ($t_k + c_i h$). Isso permite que a aproximação da derivada fracionária se alinhe perfeitamente com os nós internos usados na discretização de Galerkin da ação conservativa.
3. Derivação das Equações: Os autores derivam as equações de Euler-Lagrange Discretas Fracionárias (FDEL) aplicando o princípio de Hamilton restrito discreto. Chave para esta derivação são novas identidades discretas estabelecidas para RKCQ:
   • Uma propriedade de semigrupo discreto para os pesos de convolução (Lema 4.1).
   • Fórmulas de integração por partes discretas assimétricas (Lemmas 4.2 e 4.3), que são necessárias para lidar com a natureza não local das derivadas fracionárias dentro da estrutura variacional.

Principais Contribuições
O artigo apresenta quatro contribuições primárias:

1. Construção de Integradores Compatíveis: Os autores constroem FVIs onde a estrutura de estágios da RKCQ (especificamente Lobatto IIIC) coincide com os nós internos da discretização de Galerkin. Essa compatibilidade estrutural evita a saturação de ordem observada em abordagens anteriores baseadas em BDFCQ.
2. Fundamentos Teóricos: O artigo estabelece as identidades discretas necessárias (propriedade de semigrupo e integração por partes) para derivar as equações FDEL dentro do contexto de RKCQ.
3. Análise de Convergência:
   • Teórica: Uma prova rigorosa de convergência de segunda ordem é fornecida para o esquema Lobatto IIIC de 2 estágios aplicado a um Lagrangiano quadrático com $\beta = 2\alpha \in (0, 1]$.
   • Numérica: Experimentos numéricos demonstram que o método atinge precisão de ordem superior (4ª e 6ª ordem) para esquemas Lobatto IIIC de 3 e 4 estágios, respectivamente, quando aplicado a problemas suaves.
4. Esquema de Ponto Médio: Uma breve discussão e implementação de um integrador fracionário baseado em ponto médio (MIDCQ) são incluídas, que também atinge precisão de segunda ordem.

Resultados
Os esquemas propostos foram testados em dois modelos: um oscilador harmônico forçado bidimensional e a equação de Bagley-Torvik fracionária.

• Precisão: Para o esquema Lobatto IIIC de 2 estágios, os resultados numéricos confirmaram a convergência de segunda ordem teórica. Para os esquemas de 3 e 4 estágios, as ordens de convergência observadas foram 4 e 6, respectivamente, correspondendo às expectativas teóricas para os métodos de Galerkin subjacentes.
• Decaimento de Energia: Os integradores reproduziram com sucesso o comportamento esperado de decaimento de energia de sistemas com amortecimento fracionário. Em simulações de longo prazo, o FVI proposto manteve melhores características de decaimento de energia comparado aos métodos clássicos RK4.
• Dependência de Regularidade: O estudo destacou que a convergência de alta ordem é contingente na regularidade da solução. Para problemas com força descontínua (ex: funções degrau), a ordem de convergência saturou em primeira ordem, consistente com os requisitos de regularidade da teoria de RKCQ subjacente. No entanto, para forças suaves, a precisão de ordem superior total foi realizada.
• Comparação com BDFCQ: Os resultados confirmam que a abordagem RKCQ supera a limitação de ordem da abordagem BDFCQ ao utilizar plenamente a informação dos estágios internos.

Significância e Alegações
O artigo alega que a principal significância deste trabalho reside na compatibilidade estrutural entre a aproximação da derivada fracionária e a discretização conservativa. Ao utilizar RKCQ baseada em Lobatto IIIC, o método incorpora naturalmente os mesmos nós de estágio usados pelo Lagrangiano discreto de Galerkin. Isso permite que tanto os termos conservativos quanto os fracionários sejam avaliados em uma grade comum, evitando assim a saturação de ordem que aflige abordagens baseadas em passos múltiplos como a BDFCQ.

Os autores observam modestamente que, embora evidências numéricas apoiem a convergência de 4ª e 6ª ordem para esquemas de mais estágios, uma prova teórica rigorosa para essas ordens superiores permanece um problema em aberto. Adicionalmente, o custo computacional da implementação atual é de $O(N^2)$ devido à dependência histórica dos pesos de convolução, embora os autores reconheçam que técnicas de convolução rápida poderiam potencialmente reduzir isso para $O(N \log^2 N)$ em trabalhos futuros. O artigo conclui que esta abordagem fornece um arcabouço robusto para a construção de integradores de alta ordem e preservação de estrutura para sistemas lagrangianos com amortecimento fracionário.