A fully diagonalized spectral method on the unit ball

Resumo Técnico: Um Método Espectral Totalmente Diagonalizado no Balão Unitário

Enunciado do Problema
O artigo aborda a solução numérica de um problema de valor de contorno de Dirichlet (BVP) para uma equação do tipo Schrödinger estacionária no balão unitário $B_d \subset \mathbb{R}^d$. A equação específica é:
$$-\Delta u + \lambda(1 - \|x\|^2)^\kappa u = f(x), \quad x \in B_d$$
sujeita à condição de contorno $u(x) = 0$ na esfera $\partial B_d = S^{d-1}$, onde $\lambda > 0$ e $\kappa$ é um inteiro não negativo. A formulação variacional padrão envolve encontrar $u \in H^1_0(B_d)$ tal que uma forma bilinear específica seja igual a um funcional linear.

Metodologia
O autor propõe um método espectral que evita a solução de grandes sistemas lineares densos tipicamente associados à discretização de equações diferenciais parciais (PDEs) no balão. A metodologia procede da seguinte forma:

1. Estrutura Variacional e Produto Interno de Sobolev:
   A forma bilinear associada ao problema variacional é reconhecida como definindo um produto interno em $H^1_0(B_d)$. Para impor a condição de contorno $u=0$ na esfera, as funções de teste são escolhidas do espaço $R^d_n = \{ (1-\|x\|^2)p(x) : p(x) \in \Pi^d_n \}$. Um produto interno de Sobolev específico $\langle \cdot, \cdot \rangle_S$ é definido no espaço polinomial $\Pi^d$ por:
   $$\langle u, v \rangle_S = \lambda \int_{B_d} (1-\|x\|^2)^{\kappa+2} u v \, dx + \int_{B_d} \nabla((1-\|x\|^2)u) \cdot \nabla((1-\|x\|^2)v) \, dx$$
   Este produto interno surge naturalmente da formulação variacional.

2. Construção de Polinômios Ortogonais de Sobolev:
   O cerne do método é a construção de uma base para $R^d_n$ consistindo de polinômios mutuamente ortogonais em relação a $\langle \cdot, \cdot \rangle_S$.

   • Estrutura: A base é construída usando harmônicos esféricos $Y^{n-2j}_\nu$ e polinômios univariados $q_j$ dependentes da variável radial $r = \|x\|$. Os polinômios assumem a forma $R^n_{j,\nu}(x) = q_j(2\|x\|^2 - 1)Y^{n-2j}_\nu(x)$.
   • Redução ao Caso Univariado: É provado que a condição de ortogonalidade de Sobolev multivariada reduz-se a um produto interno de Sobolev univariado para os polinômios radiais $q_j$. Este produto interno univariado envolve um termo $L^2$ ponderado e um termo envolvendo a derivada de $(1-t)q(t)$.

3. Fórmulas de Conexão e Coeficientes Recursivos:
   O artigo estabelece uma conexão entre estes novos polinômios ortogonais de Sobolev e polinômios de Jacobi clássicos (especificamente $P^{(1, \beta)}_k$).

   • Os polinômios de Sobolev $q^{(\beta)}_k$ são expressos como uma combinação linear de polinômios de Jacobi clássicos.
   • Crucialmente, a conexão é "bandada" (banded): $P^{(1, \beta)}_k$ está relacionado a $q^{(\beta)}_k$ e apenas aos $\kappa+1$ polinômios precedentes.
   • Um algoritmo é fornecido para calcular os coeficientes de conexão recursivamente. Para o caso específico $\kappa=0$, uma representação de forma fechada e uma relação de recorrência não linear para os coeficientes são derivadas, permitindo que sejam expressos via frações contínuas.

4. Solução Totalmente Diagonalizada:
   A solução $u(x)$ é expandida na base de polinômios ortogonais de Sobolev multiplicada pelo fator de contorno $(1-\|x\|^2)$. Como a base é ortogonal em relação ao produto interno natural do operador, a projeção de Galerkin leva a um sistema diagonal. Os coeficientes de Fourier-Sobolev da solução são derivados recursivamente dos coeficientes de Fourier do termo de origem $f(x)$ na base ortogonal clássica, eliminando a necessidade de inverter uma matriz.

Principais Contribuições

• Diagonalização: A principal contribuição é um método espectral que resulta em um sistema totalmente diagonalizado para o problema de Dirichlet no balão. Isso contrasta com métodos espectrais padrão que frequentemente resultam em matrizes densas e mal condicionadas.
• Polinômios Ortogonais de Sobolev: O artigo fornece uma teoria construtiva para polinômios ortogonais de Sobolev multivariados no balão unitário em relação a um produto interno envolvendo gradientes e um peso específico $(1-\|x\|^2)^\kappa$.
• Algoritmos Recursivos: Estabelece fórmulas de conexão explícitas entre estes polinômios de Sobolev e polinômios de Jacobi clássicos, permitindo o cálculo recursivo dos coeficientes de expansão sem resolver sistemas lineares.
• Verificação Numérica: Um experimento numérico para o caso $d=2, \kappa=0$ demonstra a precisão do método. A análise de erro mostra convergência exponencial (super-exponencial, dada a analiticidade da solução) conforme o grau da aproximação aumenta.

Significância e Alegações
O autor afirma que esta abordagem aproveita a estrutura variacional do problema para simplificar significativamente a implementação numérica. Ao utilizar polinômios ortogonais de Sobolev, o método evita o gargalo computacional de resolver grandes sistemas lineares, que é uma limitação comum em métodos espectrais para PDEs no balão. O artigo afirma que as fórmulas de conexão derivadas permitem a derivação eficiente e recursiva dos coeficientes da solução, tornando o método "totalmente diagonalizado". Os resultados numéricos confirmam que o método é confiável e preciso, alcançando alta precisão com aproximações de grau relativamente baixo para soluções suaves.