A fully diagonalized spectral method on the unit ball

Resumen Técnico: Un Método Espectral Completamente Diagonalizado en la Bola Unitaria

Planteamiento del Problema
El artículo aborda la solución numérica de un problema de valor en la frontera de Dirichlet (BVP) para una ecuación de tipo Schrödinger estacionaria en la bola unitaria $B_d \subset \mathbb{R}^d$. La ecuación específica es:
$$-\Delta u + \lambda(1 - \|x\|^2)^\kappa u = f(x), \quad x \in B_d$$
sujeta a la condición de contorno $u(x) = 0$ en la esfera $\partial B_d = S^{d-1}$, donde $\lambda > 0$ y $\kappa$ es un entero no negativo. La formulación variacional estándar consiste en encontrar $u \in H^1_0(B_d)$ tal que una forma bilineal específica sea igual a un funcional lineal.

Metodología
El autor propone un método espectral que evita la resolución de grandes sistemas lineales densos típicamente asociados con la discretización de ecuaciones diferenciales parciales (PDE) en la bola. La metodología procede de la siguiente manera:

1. Marco Variacional y Producto Interno de Sobolev:
   Se reconoce que la forma bilineal asociada al problema variacional define un producto interno en $H^1_0(B_d)$. Para imponer la condición de contorno $u=0$ en la esfera, se eligen las funciones de prueba del espacio $R^d_n = \{ (1-\|x\|^2)p(x) : p(x) \in \Pi^d_n \}$. Se define un producto interno de Sobolev específico $\langle \cdot, \cdot \rangle_S$ en el espacio polinomial $\Pi^d$ como:
   $$\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 producto interno surge naturalmente de la formulación variacional.

2. Construcción de Polinomios Ortogonales de Sobolev:
   El núcleo del método es la construcción de una base para $R^d_n$ compuesta por polinomios mutuamente ortogonales con respecto a $\langle \cdot, \cdot \rangle_S$.

   • Estructura: La base se construye utilizando armónicos esféricos $Y^{n-2j}_\nu$ y polinomios univariados $q_j$ que dependen de la variable radial $r = \|x\|$. Los polinomios toman la forma $R^n_{j,\nu}(x) = q_j(2\|x\|^2 - 1)Y^{n-2j}_\nu(x)$.
   • Reducción al Caso Univariado: Se demuestra que la condición de ortogonalidad de Sobolev multivariada se reduce a un producto interno de Sobolev univariado para los polinomios radiales $q_j$. Este producto interno de Sobolev univariado involucra un término $L^2$ ponderado y un término que implica la derivada de $(1-t)q(t)$.

3. Fórmulas de Conexión y Coeficientes Recursivos:
   El artículo establece una conexión entre estos nuevos polinomios ortogonales de Sobolev y los polinomios de Jacobi clásicos (específicamente $P^{(1, \beta)}_k$).

   • Los polinomios de Sobolev $q^{(\beta)}_k$ se expresan como una combinación lineal de polinomios de Jacobi clásicos.
   • Crucialmente, la conexión es "bandeada": $P^{(1, \beta)}_k$ está relacionado con $q^{(\beta)}_k$ y solo con los $\kappa+1$ polinomios precedentes.
   • Se proporciona un algoritmo para calcular los coeficientes de conexión recursivamente. Para el caso específico de $\kappa=0$, se deriva una representación de forma cerrada y una relación de recurrencia no lineal para los coeficientes, lo que permite expresarlos mediante fracciones continuas.

4. Solución Completamente Diagonalizada:
   La solución $u(x)$ se expande en la base de los polinomios ortogonales de Sobolev multiplicados por el factor de frontera $(1-\|x\|^2)$. Debido a que la base es ortogonal con respecto al producto interno natural del operador, la proyección de Galerkin conduce a un sistema diagonal. Los coeficientes de Fourier-Sobolev de la solución se derivan recursivamente de los coeficientes de Fourier del término fuente $f(x)$ en la base ortogonal clásica, eliminando la necesidad de invertir una matriz.

Contribuciones Clave

• Diagonalización: La principal contribución es un método espectral que resulta en un sistema completamente diagonalizado para el problema de Dirichlet en la bola. Esto contrasta con los métodos espectrales estándar que a menudo producen matrices completas y mal condicionadas.
• Polinomios Ortogonales de Sobolev: El artículo proporciona una teoría constructiva para polinomios ortogonales de Sobolev multivariados en la bola unitaria con respecto a un producto interno que involucra gradientes y un peso específico $(1-\|x\|^2)^\kappa$.
• Algoritmos Recursivos: Establece fórmulas de conexión explícitas entre estos polinomios de Sobolev y los polinomios de Jacobi clásicos, permitiendo el cálculo recursivo de los coeficientes de expansión sin resolver sistemas lineales.
• Verificación Numérica: Un experimento numérico para el caso $d=2, \kappa=0$ demuestra la precisión del método. El análisis de error muestra una convergencia exponencial (super-exponencial, dada la analiticidad de la solución) a medida que aumenta el grado de la aproximación.

Significancia y Reivindicaciones
El autor afirma que este enfoque aprovecha la estructura variacional del problema para simplificar significativamente la implementación numérica. Al utilizar polinomios ortogonales de Sobolev, el método evita el cuello de botella computacional de resolver grandes sistemas lineales, que es una limitación común en los métodos espectrales para las PDE en la bola. El artículo sostiene que las fórmulas de conexión derivadas permiten la derivación eficiente y recursiva de los coeficientes de la solución, haciendo que el método sea "completamente diagonalizado". Los resultados numéricos confirman que el método es fiable y preciso, logrando una alta precisión con grados de aproximación relativamente bajos para soluciones suaves.