A Residual-Certified Product-Integration Method for Nonlinear Caputo–Volterra Equations

Resumo Técnico: Um Método de Integração de Produto com Certificação de Resíduo para Equações de Caputo–Volterra Não Lineares

Enunciado do Problema
O artigo aborda a solução numérica de problemas de valor inicial de Caputo não lineares, que podem ser reformulados como equações integrais de Volterra com núcleos fracamente singulares do tipo Abel. Especificamente, o estudo foca na equação:
$$C D^\alpha_{0+} u(t) = F(t, u(t)), \quad u(0) = u_0, \quad 0 < \alpha < 1$$
onde $F$ satisfaz uma condição de Lipschitz local. O principal desafio abordado é resolver o sistema algébrico não linear resultante da discretização, particularmente quando a constante de Lipschitz global é grande (excedendo o limiar exigido para os teoremas padrão de ponto fixo de Banach) e quando a solução apresenta singularidades iniciais. O artigo busca determinar se uma regra de integração de produto de ordem superior pode ser implementada de forma causal, fornecendo simultaneamente um limite rigoroso e computável para soluções não lineares inexatas.

Metodologia
Os autores propõem um esquema de integração de produto implícito usando interpolação linear por partes do termo não linear $F(t, u(t))$ em uma malha não uniforme.

1. Discretização: O núcleo de Abel é integrado exatamente contra o interpolante linear. Isso resulta em um sistema não linear triangular inferior onde o valor no nó atual $U_i$ aparece implicitamente.
2. Condição de Solvabilidade: Diferente das regras de extremidade esquerda estritamente causais, o termo do nó atual introduz um coeficiente diagonal $\kappa_i L_i$. O método baseia-se em um teste de solvabilidade local da malha: se $\mu_h = \max_i (\kappa_i L_i) < 1$, o sistema discreto possui uma solução única. Esta condição é mais fraca do que a condição global padrão $S = L^* T^\alpha / \Gamma(\alpha+1) < 1$.
3. Solvers Iterativos: Tanto as iterações de Picard globais quanto as iterações de Picard nodais (sequenciais) são provadas como convergentes sob a condição local $\mu_h < 1$, mesmo quando a constante global $S > 1$.
4. Certificação de Resíduo: Um componente fundamental é a derivação de um limite de erro algébrico componente a componente. Para qualquer vetor discreto inexato $W$, o erro $|W - U_h|$ é limitado por $(I - C_h)^{-1} |W - T_h W|$, onde $C_h$ é uma matriz triangular inferior não negativa derivada das constantes de Lipschitz e dos pesos de quadratura. Este limite é computado via uma única substituição direta (forward substitution), separando o erro devido à iteração não linear incompleta do erro de discretização.
5. Análise de Convergência:
   • Para soluções suaves, o método alcança convergência de segunda ordem em malhas suficientemente finas.
   • Para soluções com singularidades iniciais do tipo potência ($|g''(t)| \le C t^{\sigma-2}$), o erro em uma malha graduada ($t_j = T(j/M)^r$) é de ordem $O(M^{-\min\{r\sigma, 2\}})$.

Principais Contribuições
O artigo afirma que sua novidade não reside na regra de quadratura em si (que é clássica) ou no uso de malhas graduadas, mas na combinação dos seguintes elementos:

• Teste de Diagonal Exata: Uma condição explícita e local da malha para a existência e unicidade da solução discreta que não requer que a constante de Lipschitz global seja pequena.
• Argumento de Contração Ponderada: Uma prova de que as iterações de Picard globais e nodais convergem sob a condição local $\mu_h < 1$.
• Certificado de Resíduo de Grade Única: Um limite componente a componente computável para o erro algébrico de qualquer iterado inexato. Diferente de estimadores baseados em extrapolação, este não requer uma malha refinada ou expansões de erro assintóticas.
• Separação de Erro: Uma estimativa de dois níveis que separa explicitamente o erro introduzido pelo término da iteração não linear do erro de discretização.
• Limite de Resolvente Independente da Malha: Uma estimativa de Volterra em blocos provando que o limite do resolvente permanece independente do tamanho da malha, fechando a prova de convergência.

Resultados
Experimentos numéricos validam os achados teóricos através de quatro modelos:

1. Confiabilidade do Certificado: O limite certificado excede consistentemente o erro algébrico observado para iterações de Picard globais, mesmo em casos onde a constante de Lipschitz global $S$ é significamente maior que 1 (por exemplo, $S \approx 7.7$). A "efetividade" do estimador (razão entre o limite e o erro) melhora conforme as iterações progridem.
2. Taxas de Convergência:
   • Em benchmarks suaves, o método proposto alcança convergência de segunda ordem, superando o método Adams PECE fracionário (ordem $\approx 1+\alpha$) e a regra de extremidade esquerda explícita (ordem 1).
   • Em problemas com singularidades iniciais, malhas uniformes resultam em ordens de convergência reduzidas (ex: $\approx 1.2$), enquanto malhas graduadas recuperam com sucesso a taxa de segunda ordem.
3. Problemas Não Fabricados: O método performa robustamente em um problema não linear sem uma solução exata conhecida, com contagens de iteração diminuindo conforme a malha é refinada.
4. Custo Computacional: O método implícito é computacionalmente mais caro que as bases explícitas (Adams PECE, extremidade esquerda) devido às resoluções não lineares e construção de pesos. No entanto, ele fornece um critério de parada rigoroso para o solver não linear e uma precisão superior em comparação com a regra de extremidade esquerda.

Significância e Alegações
O artigo posiciona modestamente sua contribuição como um refinamento teórico e algorítmico de técnicas existentes de integração de produto, em vez de uma descoberta de novas regras de quadratura. Sua significância reside em:

• Estender a aplicabilidade da integração de produto implícita para casos onde a contração global falha ($S > 1$) ao utilizar propriedades locais da malha.
• Fornecer um critério de parada rigoroso e computável para solvers não lineares que não depende de suposições assintóticas ou computações adicionais em malhas refinadas.
• Oferecer um framework unificado que separa erros algébricos e de discretização, permitindo uma estimativa de erro confiável para equações de Volterra fracionárias não lineares.

Os autores declaram explicitamente que não reivindicam o primeiro uso de malhas graduadas, transformações de suavização ou regras de integração de produto para estas equações, mas sim a integração específica de um teste de solvabilidade diagonal, provas de contração ponderada e um certificado de resíduo de grade única.