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

Technische Samenvatting: Analytische en Gegenbauer Spectrale Oplossingen voor Gedempte Golfvergelijkingen met Singuliere Potentialen in $\mathbb{R}^n$

Probleemformulering
Deze studie onderzoekt de existentie, uniciteit en numerieke benadering van oplossingen voor een meerdimensionale gedempte golfvergelijking gedefinieerd op een rechthoekig domein $\Omega = [0, \pi]^n$. De leidende vergelijking is een gemengde hyperbolisch-parabolische operator die een singuliere bronterm van de vorm $|x|^{-\alpha}$ bevat:
$$\frac{\partial^2 u}{\partial t^2} + \gamma \frac{\partial u}{\partial t} - \Delta u + \frac{u}{|x|^\alpha} = f(x, t)$$
Het probleem wordt geformuleerd onder Cauchy-beginvoorwaarden ($u(x,0)=\phi(x)$, $u_t(x,0)=\psi(x)$) en niet-homogene Neumann-randvoorwaarden ($\partial u / \partial \nu = g_j$). De analyse richt zich op het regime waar de singulariteitsparameter voldoet aan $0 < \alpha < n$ en de dempingscoëfficiënt $\gamma > 0$. De aanwezigheid van de singuliere potentiaal introduceert aanzienlijke analytische uitdagingen met betrekking tot integreerbaarheid en regulariteit, terwijl het gemengde karakter van de operator vereist dat zowel golfvoortplanting als diffusieve dissipatie simultaan worden geadresseerd.

Methodologie
Het artikel hanteert een duale analytische en numerieke benadering:

1. Analytisch Kader (Eigenfunctie-expansie):

   • De oplossing wordt ontleed in een transiënte homogene component ($u_h$) en een stationaire particuliere component ($u_p$).
   • De particuliere oplossing wordt geconstrueerd met behulp van een Fourierreeks in de tijd, waardoor de tijdsafhankelijke PDE wordt gereduceerd tot een systeem van gekoppelde elliptische vergelijkingen voor elke temporele harmonische $k$.
   • Door complexe-waardige ruimtelijke functies $U_k(x) = A_k(x) - iB_k(x)$ te introduceren, wordt het systeem gereduceerd tot een complexe Helmholtz-type vergelijking: $(L - \lambda_k)U_k = F_k$, waarbij $L = -\Delta + |x|^{-\alpha}$ en $\lambda_k = k^2 - i\gamma k$.
   • De oplossing wordt expliciet weergegeven als een dubbele reeks met behulp van de eigenfuncties $\Phi_m(x) = \prod \cos(m_j x_j)$ van de Neumann-Laplace-operator. De coëfficiënten worden bepaald via projectie op de orthogonale basis, waarbij rekening wordt gehouden met zowel interne bronnen als randfluxen.

2. Numeriek Kader (Gegenbauer Spectrale Collocatie):

   • Er wordt een hogere-orde spectraal methode ontwikkeld met behulp van Gegenbauer-polynomen met parameter $\lambda=1$ (wat reduceert tot Chebyshev-polynomen voor een optimale node-distributie).
   • Ruimtelijke discretisatie maakt gebruik van Chebyshev-Gauss-Lobatto (CGL) nodes op een tensorproduct-grid om het domein-geometrie te behanden.
   • De singuliere potentiaal $|x|^{-\alpha}$ wordt nabij de oorsprong geregulariseerd met behulp van $V_\epsilon(x) = (|x|^2 + \epsilon^2)^{-\alpha/2}$ om numerieke overflow te voorkomen terwijl het asymptotische gedrag behouden blijft.
   • Tijdintegratie wordt uitgevoerd met behulp van een tweede-orde Implicit-Explicit (IMEX) schema, waarbij de stijve diffusieterm impliciet wordt behandeld en de hyperbolische termen expliciet, om stabiliteit te garanderen zonder restrictieve CFL-condities.

Belangrijkste Bijdragen en Theoretische Resultaten
Het artikel stelt rigoureuze wiskundige eigenschappen van de oplossing vast:

• Convergentie van Coëfficiënten: Er wordt bewezen dat de expansiecoëfficiënten asymptotisch vervallen als $O(|m|^{-n-\alpha})$, wat ervoor zorgt dat de reeks goed gedefinieerd is ondanks de singulariteit, mits $0 < \alpha < n$.
• Stabiliteit en Niet-Resonantie: De dempingscoëfficiënt $\gamma > 0$ blijkt cruciaal. Het zorgt ervoor dat de complexe spectrale noemer $|\mu_m - \lambda_k| \geq \gamma k > 0$ is, wat deling door nul voorkomt en pure resonantie-fenomenen elimineert.
• Regulariteit en Sobolev-ruimten:
  • De oplossing $u$ is begrensd en wel-gesteld in $L^2(\Omega)$.
  • De ruimtelijke gradiënt $\nabla u$ blijft kwadratisch integreerbaar (d.w.z. $u \in H^1(\Omega)$) alleen onder de striktere voorwaarde $0 < \alpha < n-2$. Voor $n-2 \leq \alpha < n$ vertoont de gradiënt een singulariteit bij de oorsprong.
• Globale Stabiliteit: Een expliciete norm-schatting is afgeleid die aantoont dat de $L^2$-norm van de oplossing exponentieel vervalt als $e^{-\gamma t/2}$, wat globale stabiliteit en energie-dissipatie bevestigt.
• Foutanalyse: Theoretische foutgrenzen voor de Gegenbauer-collocatiemethode zijn vastgesteld. De ruimtelijke discretiefout convergeert algebraïsch met een snelheid van $O(N^{-(n-\alpha-1/2)})$, wat de regulariteitsverlies reflecteert die door de singulariteit wordt geïnduceerd. Voor zwakke singulariteiten wordt bijna-spectrale convergentie waargenomen.

Numerieke Resultaten
Uitgebreide numerieke experimenten zijn uitgevoerd in één, twee en drie dimensies.

• Validatie: De simulaties valideren de theoretische convergentiesnelheden. Voor zwakke singulariteiten ($\alpha = 0.3$) bereikt de methode een hoge-orde convergentie (snelheden $> 3.5$). Voor sterke singulariteiten nabij de kritische drempel ($\alpha \approx n - 0.5$) degradeert de convergentiesnelheid gracieus naar de theoretische algebraïsche limiet.
• Energiegedrag: Numerieke plots bevestigen de exponentiële afname van energie over de tijd, consistent met de dempingsparameter $\gamma$.
• Dimensionale Schaling: De resultaten tonen aan dat de methode robuust en stabiel blijft over dimensies $n=1, 2, 3$, waarbij de computationele kosten schalen als $O(N^n)$.

Betekenis
Het artikel beweert een gat in de literatuur te vullen door een expliciete analytische representatie te bieden voor $n$-dimensionale gemengde hyperbolisch-parabolische operatoren onder singuliere potentialen, een probleem waarbij standaard existentie-theorema's vaak tekortschieten in het bieden van constructieve oplossingen. Het werk demonstreert dat de dempingscoëfficiënt fungeert als een fundamentele stabilisator tegen resonantie in de aanwezigheid van singulariteiten. Bovendien biedt de geïntegreerde analytisch-numerieke benadering een robuust model voor het karakteriseren van golfvoortplanting en diffusie in niet-homogene media met gelokaliseerde singulariteiten. De auteurs positioneren deze bevindingen als direct toepasbaar op fysieke velden zoals thermoelasticiteit, vloeistof-structuur interactie en elektromagnetische golfvoortplanting in dissipatieve media.