On a Boundary-Initial Value Problem for Fractional Differential Equation with Sequential Caputo derivatives

Resumo Técnico: Problema de Valor de Contorno-Inicial para Equação Diferencial Fracionária com Derivadas Caputo Sequenciais

Autores: Fayziev Yusuf e Jumaeva Shakhnoza.
Tema: Equações Diferenciais Fracionárias (EDFs), Derivadas Caputo Sequenciais, Função Mittag-Leffler Bivariada, Método de Elementos Finitos (FEM).

1. O Problema Investigado

O artigo foca na análise de um problema de valor de contorno-inicial (BVIP) para uma equação diferencial fracionária que envolve derivadas Caputo sequenciais. A equação modelada é:

$$D^\beta_t (D^\alpha_t u(x, t)) + D^\beta_t u(x, t) - u_{xx}(x, t) = f(x, t)$$

Definida no domínio $\Omega = \{(x, t) : 0 < x < 1, 0 < t < T\}$, com as seguintes condições:

• Condições Iniciais: $u(x, 0) = \phi(x)$ e $D^\alpha_t u(x, 0) = \psi(x)$.
• Condições de Contorno: $u(0, t) = u(1, t) = 0$.
• Parâmetros: $0 < \alpha, \beta < 1$. O estudo considera especificamente a região admissível onde $\alpha + \beta > 1$ e $\alpha > \beta$.

O objetivo é encontrar uma solução regular $u(x, t)$ e analisar suas propriedades de existência, unicidade e comportamento numérico. A motivação reside na necessidade de modelar processos com múltiplos efeitos de memória e escalas de tempo, onde a propriedade de semigrupo das derivadas fracionárias clássicas não se aplica, exigindo o uso de operadores sequenciais.

2. Metodologia

A abordagem do artigo divide-se em duas vertentes principais: análise teórica e implementação numérica.

A. Análise Teórica (Solução Analítica)

• Método de Separação de Variáveis (Fourier): O problema é transformado em uma sequência de problemas de valor inicial para equações diferenciais ordinárias (EDOs) fracionárias através da expansão em séries de senos.
• Função Mittag-Leffler Bivariada: A solução exata das EDOs auxiliares é expressa em termos de uma generalização da função Mittag-Leffler, conhecida como Função Mittag-Leffler Bivariada ($E_2$).
• Lemas Auxiliares: Os autores estabelecem propriedades fundamentais desta função, incluindo representações integrais do tipo Hankel, estimativas de convergência e regras de diferenciação/integração fracionária aplicadas a esta função.
• Prova de Existência e Unicidade: Utilizando o método de Fourier e análise de convergência de séries, prova-se que a solução formal converge uniformemente no domínio fechado, garantindo a existência e unicidade da solução regular sob condições de regularidade específicas para os dados iniciais e forçantes (espaços de Hölder).

B. Implementação Numérica

• Reformulação Sequencial: Para evitar a discretização direta de termos de história aninhados (complexos), o operador sequencial é tratado introduzindo uma variável auxiliar $v(x, t) = D^\alpha_t u(x, t)$. Isso converte a equação de ordem superior em um sistema acoplado de primeira ordem.
• Discretização Espacial: Utiliza-se o Método de Elementos Finitos (FEM) com elementos lineares (P1) para a variável espacial.
• Discretização Temporal:
  • Emprega-se uma malha temporal graduada (graded mesh) para lidar com a baixa regularidade (singularidade) da solução próxima a $t=0$.
  • A aproximação das derivadas fracionárias utiliza o Esquema L1 adaptado para malhas não uniformes.
• Algoritmo: Um esquema totalmente discreto é desenvolvido, resolvendo sistemas lineares em cada passo de tempo, atualizando a variável auxiliar $v$ a partir de $u$.

3. Contribuições Principais

1. Solução Analítica Exata: Derivação da solução exata do problema em termos da função Mittag-Leffler bivariada, preenchendo uma lacuna na literatura sobre a resolução analítica de equações com derivadas sequenciais Caputo.
2. Estabelecimento de Propriedades da Função $E_2$: Formulação e prova de lemas cruciais sobre a função Mittag-Leffler bivariada (convergência, limites assintóticos e operações de cálculo fracionário), que são essenciais para a construção da solução.
3. Prova de Bem-Pose: Demonstração rigorosa da existência e unicidade da solução regular para o problema proposto, sob as condições $\alpha + \beta > 1$ e $\alpha > \beta$.
4. Esquema Numérico Robusto: Desenvolvimento de um esquema numérico eficiente que combina a reformulação sequencial com malhas graduadas, superando as dificuldades de precisão associadas às singularidades iniciais típicas de problemas fracionários.

4. Resultados

• Validação Numérica: O esquema foi validado utilizando uma solução manufactured (solução conhecida). Os resultados mostraram alta precisão, com erros máximos na ordem de $10^{-6}$, confirmando a confiabilidade do método.
• Impacto da Malha Graduada: A comparação entre malhas uniformes e graduadas demonstrou que a malha graduada (com parâmetro de graduação $r \approx 1.6$) é superior. Ela restaura a ordem de convergência para valores acima de 1.10 (em contraste com ~0.75 para malhas uniformes) e reduz o erro em um fator de aproximadamente 7.75 para um mesmo número de passos temporais.
• Comportamento Dinâmico:
  • Efeito de $\alpha$ e $\beta$: A análise paramétrica revelou que os parâmetros não apenas alteram a velocidade do processo, mas redefinem fundamentalmente a evolução da solução.
  • Transição de Regime: Observou-se uma transição de um decaimento puramente difusivo para um regime oscilatório amortecido à medida que as ordens fracionárias aumentam.
  • Efeito de "Rebote" (Inércia): Para valores mais altos de $\alpha$ e $\beta$, a energia do sistema decai a zero e apresenta um leve "rebound" (ressalto), indicando uma forma de inércia no modelo, onde o sistema ultrapassa o ponto de equilíbrio antes de estabilizar.
  • Trade-off de Parâmetros: Mapeamentos de calor mostraram que diferentes combinações de $\alpha$ e $\beta$ podem produzir o mesmo resultado em um ponto específico, o que é útil para problemas de identificação inversa.

5. Significância

Este trabalho é significativo por fornecer uma base teórica e numérica sólida para o estudo de equações diferenciais fracionárias com derivadas sequenciais, que são modelos mais realistas para sistemas com múltiplos mecanismos de memória do que as equações fracionárias padrão.

• Avanço Teórico: A conexão explícita entre operadores sequenciais e a função Mittag-Leffler bivariada abre novas portas para a análise de sistemas complexos.
• Aplicabilidade Prática: O esquema numérico proposto é eficiente e preciso, permitindo a simulação de fenômenos físicos complexos (como ondas de difusão-fracionária e equações de telegrafista fracionárias) que antes eram difíceis de tratar numericamente devido às singularidades iniciais.
• Insights Físicos: A descoberta de comportamentos de "inércia" e oscilação em modelos puramente difusivos sugere que a escolha correta das ordens fracionárias sequenciais é crítica para modelar corretamente a dinâmica de materiais viscoelásticos e processos de transporte anômalo.