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

तकनीकी सारांश: $\mathbb{R}^n$ में डैम्प्ड वेव इक्वेशन्स के विश्लेषणात्मक और गेनबाकर स्पेक्ट्रल समाधान

समस्या का निरूपण
यह अध्ययन एक आयताकार डोमेन $\Omega = [0, \pi]^n$ पर परिभाषित बहु-आयामी डैम्प्ड वेव इक्वेशन (damped wave equation) के समाधानों के अस्तित्व, विशिष्टता और संख्यात्मक सन्निकटन की जांच करता है। नियंत्रित समीकरण एक मिश्रित हाइपरबोलिक-पैराबोलिक ऑपरेटर है जिसमें $|x|^{-\alpha}$ के रूप में एक विलक्षण स्रोत पद (singular source term) शामिल है:
$$\frac{\partial^2 u}{\partial t^2} + \gamma \frac{\partial u}{\partial t} - \Delta u + \frac{u}{|x|^\alpha} = f(x, t)$$
समस्या को कॉची प्रारंभिक स्थितियों (Cauchy initial conditions) ($u(x,0)=\phi(x)$, $u_t(x,0)=\psi(x)$) और गैर-समांगी न्यूमैन सीमा स्थितियों (non-homogeneous Neumann boundary conditions) $\partial u / \partial \nu = g_j$ के तहत निरूपित किया गया है। विश्लेषण उस शासन (regime) पर केंद्रित है जहाँ विलक्षणता पैरामीटर (singularity parameter) $0 < \alpha < n$ संतुष्ट करता है और डैम्पिंग गुणांक $\gamma > 0$ है। विलक्षण विभव (singular potential) की उपस्थिति एकीकरण (integrability) और नियमितता (regularity) के संबंध में महत्वपूर्ण विश्लेषणात्मक चुनौतियाँ पेश करती है, जबकि ऑपरेटर की मिश्रित प्रकृति तरंग प्रसार (wave propagation) और विसारक अपव्यय (diffusive dissipation) दोनों को एक साथ संबोधित करने की आवश्यकता को दर्शाती है।

कार्यप्रणाली
यह शोध एक दोहरी विश्लेषणात्मक और संख्यात्मक पद्धति का उपयोग करता है:

1. विश्लेषणात्मक ढांचा (इगनफंक्शन एक्सपेंशन):

   • समाधान को एक क्षणिक समांगी भाग ($u_h$) और एक स्थिर-अवस्था विशिष्ट भाग ($u_p$) में विभाजित किया जाता है।
   • विशिष्ट समाधान का निर्माण समय में फूरियर श्रृंखला (Fourier series) का उपयोग करके किया जाता है, जो समय-निर्भर PDE को प्रत्येक टेम्पोरल हार्मोनिक $k$ के लिए युग्मित दीर्घवृत्तीय समीकरणों (coupled elliptic equations) की एक प्रणाली में कम कर देता है।
   • जटिल-मान वाले स्थानिक फलनों $U_k(x) = A_k(x) - iB_k(x)$ को पेश करके, प्रणाली को एक जटिल हेल्महोल्ट्ज़-प्रकार के समीकरण में कम किया जाता है: $(L - \lambda_k)U_k = F_k$, जहाँ $L = -\Delta + |x|^{-\alpha}$ और $\lambda_k = k^2 - i\gamma k$ है।
   • समाधान को न्यूमैन लैप्लासियन के इगनफंक्शन्स $\Phi_m(x) = \prod \cos(m_j x_j)$ का उपयोग करके एक दोहरी श्रृंखला के रूप में स्पष्ट रूप से प्रदर्शित किया जाता है। गुणांकों का निर्धारण आंतरिक स्रोतों और सीमा प्रवाह (boundary fluxes) को ध्यान में रखते हुए, ऑर्थोगोनल आधार पर प्रोजेक्शन के माध्यम से किया जाता है।

2. संख्यात्मक ढांचा (गेनबाकर स्पेक्ट्रल कोलोकेशन):

   • पैरामीटर $\lambda=1$ के साथ गेनबाकर बहुपदों (Gegenbauer polynomials) का उपयोग करके एक उच्च-क्रम स्पेक्ट्रल विधि विकसित की गई है (जो इष्टतम नोड वितरण के लिए चेबिशेव बहुपदों में बदल जाती है)।
   • स्थानिक विविक्तीकरण (spatial discretization) डोमेन ज्यामिति को संभालने के लिए टेंसर-प्रोडक्ट ग्रिड पर चेबिशेव-गॉस-लोबैटो (CGL) नोड्स का उपयोग करता है।
   • मूल बिंदु (origin) के पास विलक्षण विभव $|x|^{-\alpha}$ को संख्यात्मक ओवरफ्लो को रोकने और साहसिक व्यवहार (asymptotic behavior) को बनाए रखने के लिए $V_\epsilon(x) = (|x|^2 + \epsilon^2)^{-\alpha/2}$ का उपयोग करके नियमित (regularize) किया जाता है।
   • समय एकीकरण (time integration) के लिए एक द्वितीय-क्रम इम्पलिसिट-एक्सप्लिसिट (IMEX) योजना का उपयोग किया जाता है, जो प्रतिबंधात्मक CFL शर्तों के बिना स्थिरता सुनिश्चित करने के लिए कठोर विसरण पद (stiff diffusion term) को अंतर्निहित (implicitly) और हाइपरबोलिक पदों को स्पष्ट (explicitly) रूप से उपचारित करती है।

प्रमुख योगदान और सैद्धांतिक परिणाम
यह शोध समाधान के कठोर गणितीय गुणों को स्थापित करता है:

• गुणांकों का अभिसरण (Convergence of Coefficients): यह सिद्ध किया गया है कि विस्तार गुणांक $O(|m|^{-n-\alpha})$ के रूप में स्पर्शोन्मुख रूप से (asymptotically) घटते हैं, जिससे यह सुनिश्चित होता है कि यदि $0 < \alpha < n$ हो, तो श्रृंखला सुपरिभाषित है, भले ही विलक्षणता मौजूद हो।
• स्थिरता और गैर-अनुनाद (Stability and Non-Resonance): यह दिखाया गया है कि डैम्पिंग गुणांक $\gamma > 0$ महत्वपूर्ण है। यह सुनिश्चित करता है कि जटिल स्पेक्ट्रल डिनोमिनेटर $|\mu_m - \lambda_k| \geq \gamma k > 0$, जिससे शून्य से विभाजन रुकता है और शुद्ध अनुनाद (pure resonance) की घटनाओं को समाप्त करता है।
• नियमितता और सोबोलेव स्पेस (Regularity and Sobolev Spaces):
  • समाधान $u$ को $L^2(\Omega)$ में परिबद्ध (bounded) और सुव्यवस्थित (well-posed) दिखाया गया है।
  • स्थानिक ग्रेडिएंट $\nabla u$ केवल सख्त स्थिति $0 < \alpha < n-2$ के तहत ही वर्ग-एकीकृत (square-integrable) रहता है (अर्थात $u \in H^1(\Omega)$)। यदि $n-2 \leq \alpha < n$ है, तो ग्रेडिएंट मूल बिंदु पर विलक्षणता प्रदर्शित करता है।
• वैश्विक स्थिरता (Global Stability): एक स्पष्ट नॉर्म अनुमान प्राप्त किया गया है जो दिखाता है कि समाधान का $L^2$-नॉर्म $e^{-\gamma t/2}$ के रूप में तेजी से (exponentially) घटता है, जो वैश्विक स्थिरता और ऊर्जा अपव्यय की पुष्टि करता है।
• त्रुटि विश्लेषण (Error Analysis): गेनबाकर कोलोकेशन विधि के लिए सैद्धांतिक त्रुटि सीमाएं स्थापित की गई हैं। स्थानिक विविक्तीकरण त्रुटि $O(N^{-(n-\alpha-1/2)})$ की दर से बीजगणितीय रूप से अभिसरित होती है, जो विलक्षणता द्वारा प्रेरित नियमितता हानि को दर्शाती है। कमजोर विलक्षणताओं के लिए, निकट-स्पेक्ट्रल अभिसरण (near-spectral convergence) देखा जाता है।

संख्यात्मक परिणाम
एक, दो और तीन आयामों में व्यापक संख्यात्मक प्रयोग किए गए।

• सत्यापन: सिमुलेशन सैद्धांतिक अभिसरण दरों को मान्य करते हैं। कमजोर विलक्षणताओं ($\alpha = 0.3$) के लिए, विधि उच्च-क्रम अभिसरण (दर $> 3.5$) प्राप्त करती है। महत्वपूर्ण दहलीज के पास मजबूत विलक्षणताओं ($\alpha \approx n - 0.5$) के लिए, अभिसरण दर सैद्धांतिक बीजगणितीय सीमा तक सुचारू रूप से गिर जाती है।
• ऊर्जा व्यवहार: ऊर्जा के समय के साथ होने वाले घातांतीय क्षय (exponential decay) की पुष्टि संख्यात्मक प्लॉट्स द्वारा की जाती है, जो डैम्पिंग पैरामीटर $\gamma$ के अनुरूप है।
• आयामी स्केलिंग (Dimensional Scaling): परिणाम प्रदर्शित करते हैं कि विधि $n=1, 2, 3$ आयामों में मजबूत और स्थिर रहती है, जिसमें कम्प्यूटेशनल लागत $O(N^n)$ के रूप में स्केल होती है।

महत्व
यह शोध दावा करता है कि यह साहित्य में एक अंतराल को भरता है, जो विलक्षण विभवों के तहत $n$-आयामी मिश्रित हाइपरबोलिक-पैराबोलिक ऑपरेटरों के लिए एक स्पष्ट विश्लेषणात्मक प्रतिनिधित्व प्रदान करता है—एक ऐसी समस्या जहाँ मानक अस्तित्व प्रमेय अक्सर रचनात्मक समाधान प्रदान करने में विफल रहते हैं। यह कार्य दर्शाता है कि डैम्पिंग गुणांक विलक्षणता की उपस्थिति में अनुनाद के विरुद्ध एक मौलिक स्टेबलाइज़र के रूप में कार्य करता है। इसके अतिरिक्त, एकीकृत विश्लेषणात्मक-संख्यात्मक दृष्टिकोण गैर-समांगी माध्यमों में स्थानीयकृत विलक्षणताओं के साथ तरंग प्रसार और विसरण को चरित्रित करने के लिए एक मजबूत मॉडल प्रदान करता है। लेखक इन निष्कर्षों को थर्मोइलास्टिसिटी, फ्लुइड-स्ट्रक्चर इंटरैक्शन और विसरक माध्यमों में इलेक्ट्रोमैग्नेटिक वेव प्रोपेगेशन जैसे भौतिक क्षेत्रों में सीधे लागू करने योग्य स्थिति में रखते हैं।