Analytical and Gegenbauer spectral solutions of damped wave equations with singular potentials in ℝn

Sintesi Tecnica: Soluzioni Spettrali Analitiche e di Gegenbauer per Equazioni d'Onda Smorzate con Potenziali Singolari in $\mathbb{R}^n$

Formulazione del Problema
Questo studio investiga l'esistenza, l'unicità e l'approssimazione numerica di soluzioni di un'equazione d'onda smorzata multidimensionale definita su un dominio rettangolare $\Omega = [0, \pi]^n$. L'equazione governante è un operatore misto iperbolico-parabolico che incorpora un termine sorgente singolare della forma $|x|^{-\alpha}$:
$$\frac{\partial^2 u}{\partial t^2} + \gamma \frac{\partial u}{\partial t} - \Delta u + \frac{u}{|x|^\alpha} = f(x, t)$$
Il problema è formulato sotto condizioni iniziali di Cauchy ($u(x,0)=\phi(x)$, $u_t(x,0)=\psi(x)$) e condizioni al contorno di Neumann non omogenee ($\partial u / \partial \nu = g_j$). L'analisi si concentra nel regime in cui il parametro di singolarità soddisfa $0 < \alpha < n$ e il coefficiente di smorzamento $\gamma > 0$. La presenza del potenziale singolare introduce significative sfide analitiche riguardanti l'integrabilità e la regolarità, mentre la natura mista dell'operatore richiede di affrontare simultaneamente la propagazione delle onde e la dissipazione diffusiva.

Metodologia
Il documento impiega un duplice approccio analitico e numerico:

1. Framework Analitico (Espansione in Autofunzioni):

   • La soluzione è decomposta in una parte omogenea transitoria ($u_h$) e una parte particolare a regime stazionario ($u_p$).
   • La soluzione particolare è costruita mediante una serie di Fourier nel tempo, riducendo la PDE dipendente dal tempo a un sistema di equazioni ellittiche accoppiate per ogni armonica temporale $k$.
   • Introducendo funzioni spaziali a valori complessi $U_k(x) = A_k(x) - iB_k(x)$, il sistema viene ridotto a un'equazione di tipo Helmholtz complessa: $(L - \lambda_k)U_k = F_k$, dove $L = -\Delta + |x|^{-\alpha}$ e $\lambda_k = k^2 - i\gamma k$.
   • La soluzione è rappresentata esplicitamente come una doppia serie utilizzando le autofunzioni $\Phi_m(x) = \prod \cos(m_j x_j)$ del Laplaciano di Neumann. I coefficienti sono determinati tramite proiezione su una base ortogonale, tenendo conto sia delle sorgenti interne che dei flussi di confine.

2. Framework Numerico (Collocazione Spettrale di Gegenbauer):

   • Viene sviluppato un metodo spettrale ad alto ordine utilizzando i polinomi di Gegenbauer con parametro $\lambda=1$ (che si riduce ai polinomi di Chebyshev per una distribuzione ottimale dei nodi).
   • La discretizzazione spaziale utilizza nodi di Chebyshev-Gauss-Lobatto (CGL) su una griglia prodotto tensoriale per gestire la geometria del dominio.
   • Il potenziale singolare $|x|^{-\alpha}$ viene regolarizzato vicino all'origine usando $V_\epsilon(x) = (|x|^2 + \epsilon^2)^{-\alpha/2}$ per prevenire l'overflow numerico pur preservando il comportamento asintotico.
   • L'integrazione temporale è eseguita utilizzando uno schema Implicit-Explicit (IMEX) del secondo ordine, trattando il termine di diffusione rigido (stiff) in modo implicito e i termini iperbolici in modo esplicito per garantire la stabilità senza condizioni CFL restrittive.

Contributi Chiave e Risultati Teorici
Il documento stabilisce rigorose proprietà matematiche della soluzione:

• Convergenza dei Coefficienti: Si dimostra che i coefficienti di espansione decadono asintoticamente come $O(|m|^{-n-\alpha})$, garantendo che la serie sia ben definita nonostante la singolarità, a condizione che $0 < \alpha < n$.
• Stabilità e Non-Risonanza: Il coefficiente di smorzamento $\gamma > 0$ si rivela critico. Esso assicura che il denominatore spettrale complesso $|\mu_m - \lambda_k| \geq \gamma k > 0$, prevenendo la divisione per zero ed eliminando i fenomeni di pura risonanza.
• Regolarità e Spazi di Sobolev:
  • La soluzione $u$ è dimostrata essere limitata e ben posta in $L^2(\Omega)$.
  • Il gradiente spaziale $\nabla u$ rimane quadrato-integrabile (ovvero $u \in H^1(\Omega)$) solo sotto la condizione più stringente $0 < \alpha < n-2$. Per $n-2 \leq \alpha < n$, il gradiente presenta una singolarità nell'origine.
• Stabilità Globale: Viene derivata una stima esplicita della norma che mostra come la norma $L^2$ della soluzione decada esponenzialmente come $e^{-\gamma t/2}$, confermando la stabilità globale e la dissipazione di energia.
• Analisi dell'Errore: Sono stabiliti limiti di errore teorici per il metodo di collocazione di Gegenbauer. L'errore di discretizzazione spaziale converge algebricamente con un tasso di $O(N^{-(n-\alpha-1/2)})$, riflettendo la perdita di regolarità indotta dalla singolarità. Per singolarità deboli, si osserva una convergenza quasi-spettrale.

Risultati Numerici
Sono stati condotti esperimenti numerici esaustivi in una, due e tre dimensioni.

• Validazione: Le simulazioni validano i tassi di convergenza teorici. Per singolarità deboli ($\alpha = 0.3$), il metodo raggiunge una convergenza di alto ordine (tassi $> 3.5$). Per singolarità forti vicino alla soglia critica ($\alpha \approx n - 0.5$), il tasso di convergenza degrada con grazia verso il limite algebrico teorico.
• Comportamento dell'Energia: I grafici numerici confermano il decadimento esponenziale dell'energia nel tempo, coerentemente con il parametro di smorzamento $\gamma$.
• Scalabilità Dimensionale: I risultati dimostrano che il metodo rimane robusto e stabile attraverso le dimensioni $n=1, 2, 3$, con costi computazionali che scalano come $O(N^n)$.

Significatività
Il saggio sostiene di colmare una lacuna nella letteratura fornendo una rappresentazione analitica esplicita per operatori misti iperbolici-parabolici $n$-dimensionali sotto potenziali singolari, un problema in cui i teoremi di esistenza standard spesso non riescono a fornire soluzioni costruttive. Il lavoro dimostra che il coefficiente di smorzamento agisce come un fondamentale stabilizzatore contro la risonanza in presenza di singolarità. Inoltre, l'approccio integrato analitico-numerico offre un modello robusto per caratterizzare la propagazione delle onde e la diffusione in mezzi non omogenei con singolarità localizzate. Gli autori posizionano questi risultati come direttamente applicabili a campi fisici quali la termoelasticità, l'interazione fluido-struttura e la propagazione di onde elettromagnetiche in mezzi dissipativi.