Trefftz methods with evanescent plane waves

यहाँ WAVES 2026 में प्रस्तुत शोध पत्र "Trefftz methods with evanescent plane waves" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

यह शोध पत्र हेल्महोल्ट्ज़ समीकरण ($\Delta u + \kappa^2 u = 0$) को हल करने के लिए उपयोग किए जाने वाले शास्त्रीय ट्रेफ़्त्ज़ विधियों (Trefftz methods) में निहित संख्यात्मक अस्थिरता (numerical instability) को संबोधित करता है।

• संदर्भ: ट्रेफ़्त्ज़ विधियाँ समाधानों को उन आधार कार्यों (basis functions) का उपयोग करके अनुमानित करती हैं जो PDE के सटीक समाधान होते हैं। मानक दृष्टिकोण प्रोपैगेटिव प्लेन वेव्स (PPWs) ($e^{i\kappa d \cdot x}$ जहाँ $d \in \mathbb{R}^n$) का उपयोग करता है।
• समस्या: जबकि PPWs तीव्र अभिसरण दर (convergence rates) प्रदान करते हैं, वे गंभीर संख्यात्मक अस्थिरता से ग्रस्त होते हैं। सैद्धांतिक परिणाम (Theorem 1) प्रदर्शित करते हैं कि कुछ हेल्महोल्ट्ज़ समाधानों (विशेष रूप से जिनमें "इवेनेसेंट मोड्स" शामिल होते हैं) के लिए, PPWs का उपयोग करके सटीक अनुमान लगाने के लिए अत्यधिक बड़े यूक्लिडियन नॉर्म वाले गुणांकों (coefficients) की आवश्यकता होती है।
• परिणाम: परिमित-परिशुद्धता अंकगणित (finite-precision arithmetic) में, ये बड़े गुणांक अपरिहार्य संख्यात्मक निरसन (numerical cancellation) और इल-कंडीशनिंग (ill-conditioning) की ओर ले जाते हैं, जिससे रैखिक बीजगणित तकनीकों (जैसे प्रीकंडीशनिंग) के बावजूद समाधान का स्थिर प्रतिनिधित्व करना असंभव हो जाता है।

2. कार्यप्रणाली (Methodology)

लेखक मानक PPW आधार को इवेनेसेंट प्लेन वेव्स (EPWs) से बदलने और अल्ट्रावीक वेरिएशनल फॉर्मूलेशन (UWVF) के भीतर एक विशिष्ट सैंपलिंग रणनीति लागू करने का प्रस्ताव करते हैं।

A. इवेनेसेंट प्लेन वेव्स (EPWs)

प्रोपैगेशन वेक्टर $d$ को वास्तविक इकाई गोले (real unit sphere) तक सीमित रखने के बजाय, यह विधि $d \in \mathbb{C}^n$ की अनुमति देती है जहाँ $d \cdot d = 1$ है।

• परिभाषा: एक EPW को $e^{i\kappa d \cdot x}$ के रूप में परिभाषित किया जाता है, जहाँ $d = \zeta p + i \eta e$ है।
  • $p$: प्रोपैगेशन दिशा ($p \in S^{n-1}$)।
  • $e$: इवेनेसेंस दिशा ($e \in S^{n-1}$, $p$ के लंबवत)।
  • $\eta$: इवेनेसेंस स्ट्रेंथ ($\eta \ge 0$)।
  • $\zeta = \sqrt{1+\eta^2}$: अपीयरेंट वेवनंबर स्केलिंग।
• व्यवहार: एक तरंग दिशा $p$ में वेवनंबर $\kappa\zeta$ के साथ दोलन करती है और दिशा $e$ में दर $\kappa\eta$ के साथ तेजी से घटती (decay) है।
• सैद्धांतिक आधार: PPWs के विपरीत, किसी भी सीमित डोमेन में हेल्महोल्ट्ज़ समाधान को EPWs के निरंतर सुपरपोजिशन (Herglotz transform के माध्यम से) के रूप में दर्शाया जा सकता है, जिसके गुणांक सीमित (bounded) होते हैं। यह सुनिश्चित करता है कि सटीक अनुमान के लिए अस्थिर, बड़े गुणांकों की आवश्यकता नहीं है।

B. सैंपलिंग रणनीति (2D रेसिपी)

एक निश्चित "डिग्री ऑफ फ्रीडम" (DOF) बजट $P$ के साथ एक असतत (discrete) EPW स्पेस बनाने के लिए, लेखक एक विशिष्ट सैंपलिंग रेसिपी प्रस्तावित करते हैं:

1. पैरामीटर्स: एक फूरियर ट्रंकेशन स्तर $L = P/4$ को परिभाषित करें।
2. सैंपलिंग डोमेन: $Y = [0, 2\pi) \times \{\pm 1\} \times [0, 1]$ पर एक समान माप (uniform measure) से बिंदु $(\theta, \phi, \xi)$ लें।
   • $\theta$: प्रोपैगेशन दिशा के लिए ध्रुवीय कोण (Polar angle)।
   • $\phi$: क्षय (decay) दिशा के लिए चिह्न (Sign)।
   • $\xi$: इवेनेसेंस स्ट्रेंथ के लिए पैरामीटर।
3. मैपिंग: नमूनों को भौतिक मापदंडों में बदलें:
   • $\zeta_p = \max\left(1, \frac{2L}{\kappa \text{diam}(K)}\xi_p\right)$
   • $\eta_p = \sqrt{\zeta_p^2 - 1}$
4. अनुकूलनशीलता (Adaptivity): यदि DOF बजट कम है ($2L \le \kappa \text{diam}(K)$), तो विधि डिफ़ॉल्ट रूप से शुद्ध PPWs ($\eta=0$) का उपयोग करती है। जैसे-जैसे $P$ बढ़ता है, आधार स्वचालित रूप से जटिल समाधान विशेषताओं को पकड़ने के लिए अधिक इवेनेसेंट मोड्स को शामिल करता है।

C. संख्यात्मक कार्यान्वयन (UWVF के साथ रेगुलाइजेशन)

इस विधि को अल्ट्रावीक वेरिएशनल फॉर्मूलेशन (UWVF) का उपयोग करके लागू किया गया है:

• ओवरसैंपलिंग (Oversampling): आधार कार्यों के बीच निकट-रैखिक निर्भरता (near-linear dependence) को संभालने के लिए, सिस्टम को ओवरसैंपल किया जाता है (उदाहरण के लिए, ट्रायल फंक्शन की तुलना में 10% अधिक टेस्ट फंक्शन)।
• रेगुलाइजेशन (Regularization): रैखिक प्रणाली के डायगोनल ब्लॉक मैट्रिक्स $D$ को ट्रंकेटेड सिंगुलर वैल्यू डीकंपोजिशन (SVD) का उपयोग करके रेगुलाइज किया जाता है। छोटे सिंगुलर मानों (थ्रेशोल्ड $\epsilon$ से नीचे) को हटा दिया जाता है, और एक स्यूडो-इनवर्स (pseudo-inverse) की गणना की जाती है। यह स्थानीय रूप से (एलिमेंट-वाइज) और समानांतर (parallel) में किया जाता है।
• सिस्टम: समाधान प्राप्त करने के लिए अंतिम रैखिक प्रणाली $(I - D^\dagger_\epsilon C)u_\epsilon = D^\dagger_\epsilon b$ को हल किया जाता है।

3. मुख्य योगदान

1. सैद्धांतिक औचित्य: यह शोध पत्र PPWs की सैद्धांतिक सीमा (Theorem 1) को पुख्ता करता है और Herglotz transform गुणांकों की बाध्यता (boundedness) का लाभ उठाकर EPWs को एक स्थिर विकल्प के रूप में स्थापित करता है।
2. व्यावहारिक एल्गोरिदम: यह 2D में (और 3D के लिए संदर्भ के साथ) EPW मापदंडों को चुनने के लिए एक ठोस, सरल "रेसिपी" प्रदान करता है जो मेष आकार (mesh size) और वेवनंबर के आधार पर प्रोपैगेटिव और इवेनेसेंट मोड्स के बीच संतुलन बनाता है।
3. स्थिरीकरण तकनीक: यह इल-कंडीशन्ड रैखिक प्रणालियों को सटीकता खोए बिना मजबूती से हल करने के लिए ओवरसैंपलिंग के साथ एलिमेंट-वाइज SVD रेगुलाइजेशन को जोड़ता है।

4. संख्यात्मक परिणाम

लेखक UWVF का उपयोग करके तीन परीक्षण मामलों के विरुद्ध इस विधि को मान्य करते हैं:

• पॉइंट सोर्स (फंडामेंटल सॉल्यूशन):
  • सेटअप: एक सिंगुलैरिटी (singularity) जो सीमा के पास स्थित है ($\kappa = 128$)।
  • परिणाम: जैसे-जैसे DOFs बढ़ते हैं, अस्थिरता के कारण PPW अभिसरण (convergence) रुक जाता है। EPWs समान DOFs के लिए PPWs की तुलना में $10^6$ गुना कम त्रुटि प्राप्त करते हैं।
• वेवनंबर का प्रभाव:
  • सेटअप: एप्रोक्सिमेशन स्पेस के आकार को रैखिक रूप से स्केल करते हुए ($\text{P} = 4\kappa$) $\kappa$ को बढ़ाना।
  • परिणाम: $\kappa$ बढ़ने के साथ EPW त्रुटि कम होती है, जबकि PPW त्रुटि उच्च रहती है या विचलन (diverge) करती है।
• स्कैटरिंग समस्या (Scattering Problem):
  • सेटअप: कोनों और एक कैविटी वाले नॉन-कॉन्वेक्स, साउंड-सॉफ्ट बाधा द्वारा स्कैटरिंग ($\kappa = 16$)।
  • परिणाम: PPW अभिसरण जल्दी ही स्थिर (stagnate) हो जाता है। EPWs समाधान को सफलतापूर्वक हल करते हैं, जो काफी उच्च सटीकता प्राप्त करते हैं, और कोनों की सिंगुलैरिटी तथा कैविटी प्रभावों को प्रभावी ढंग से संभालते हैं।

5. महत्व

यह कार्य उच्च-आवृत्ति तरंग प्रसार सिमुलेशन (high-frequency wave propagation simulation) में एक महत्वपूर्ण प्रगति का प्रतिनिधित्व करता है। विशुद्ध रूप से प्रोपैगेटिव बेस से इवेनेसेंट प्लेन वेव्स की ओर शिफ्ट होकर, लेखक उस मौलिक संख्यात्मक बाधा को दूर करते हैं जिसने शास्त्रीय ट्रेफ़्त्ज़ विधियों की सटीकता को सीमित किया था।

• स्थिरता: यह विधि उन क्षेत्रों में उच्च-सटीक समाधानों की अनुमति देती है जहाँ शास्त्रीय ट्रेफ़्त्ज़ विधियाँ संख्यात्मक शोर (numerical noise) के कारण विफल हो जाती हैं।
• दक्षता: यह अत्यधिक बड़े मेष या बहुपद डिग्री (polynomial degrees) की आवश्यकता के बिना उच्च सटीकता प्राप्त करता है, जिससे ट्रेफ़्त्ज़ स्कीम्स की विशेषता वाली घातांकीय अभिसरण दर (exponential convergence rates) बनी रहती है।
• अनुप्रयोग: प्रस्तावित सैंपलिंग रणनीति और रेगुलाइजेशन इस विधि को जटिल ज्यामिति, उच्च वेवनंबर और सिंगुलर समाधानों (जैसे पॉइंट सोर्स या स्कैटरिंग कॉर्नर) के लिए मजबूत बनाता है।

संक्षेप में, यह शोध पत्र सिद्ध करता है कि इवेनेसेंट प्लेन वेव्स केवल एक सैद्धांतिक जिज्ञासा नहीं हैं, बल्कि हेल्महोल्ट्ज़ समीकरण के स्थिर, उच्च-क्रम ट्रेफ़्त्ज़ विविक्तकरण (discretizations) के लिए एक व्यावहारिक आवश्यकता हैं।