The weak Galerkin method for a class of Gross-Pitaevskii type eigenvalue problems

तकनीकी सारांश: ग्रॉस-पिटावस्की प्रकार की समस्याओं के एक वर्ग के लिए वीक गैलेरकिन विधि (The Weak Galerkin Method)

समस्या विवरण
यह शोध पत्र ग्रॉस-पिटावस्की प्रकार के मॉडलों से उत्पन्न होने वाली गैररेखीय आइगेनवैल्यू (eigenvalue) समस्याओं के संख्यात्मक समाधान को संबोधित करता है। इस समस्या को $H^1_0(\Omega)$ सोबोलेव स्पेस (Sobolev space) पर एक $L^2$-नॉर्म बाधा के अधीन एक ऊर्जा फलन $E(v)$ के वेरियेशनल मिनिमाइजेशन (variational minimization) के रूप में तैयार किया गया है। ऊर्जा फलन में एक सममित मैट्रिक्स $A(x)$ और एक विभव (potential) $V(x)$ से संबंधित एक रैखिक भाग और एक फलन $F(v^2)$ द्वारा परिभाषित एक गैररेखीय भाग शामिल है। विशिष्ट लक्ष्य ग्राउंड स्टेट आइगेनपेयर $(\lambda, u)$ का सन्निकटन (approximation) करना है, जहाँ $u$ मिनिमाइज़र है और $\lambda$ संगत आइगेनवैल्यू है।

कार्यप्रणाली (Methodology)
लेखक गैररेखीय आइगेनवैल्यू समस्या को विविक्त (discretize) करने के लिए वीक गैलेरकिन (WG) परिमित तत्व विधि (finite element method) का उपयोग करते हैं। कार्यप्रणाली निम्नलिखित चरणों के माध्यम से आगे बढ़ती है:

1. विविक्तीकरण (Discretization): एक WG परिमित तत्व स्थान $V_h$ का निर्माण खंडों (elements) पर पीसवाइज पॉलिनोमियल (piecewise polynomials) का उपयोग करके किया जाता है। एक पूर्णांक $k \ge 1$ के लिए, विविक्त फलन $v_h = \{v_0, v_b\}$ प्रत्येक तत्व पर डिग्री $k$ के एक पॉलिनोमियल $v_0$ और किनारों पर डिग्री $k-1$ के एक पॉलिनोमियल $v_b$ से बना होता है।
2. वीक ग्रेडिएंट (Weak Gradient): $L^2$ प्रोजेक्शन का उपयोग करते हुए एक विविक्त वीक ग्रेडिएंट ऑपरेटर $\nabla_w$ को परिभाषित किया जाता है, जो तत्वों की सीमाओं के पार विच्छिन्न (discontinuous) सन्निकटन की अनुमति देता है।
3. स्थिरीकरण (Stabilization): विविक्त फलनों के आंतरिक और सीमा घटकों के बीच अंतर को नियंत्रित करने के लिए एक स्थिरीकरण द्विपद रूप (bilinear form) $s(v_h, w_h)$ पेश किया जाता है।
4. गैररेखीय योजना (Nonlinear Scheme): निरंतर गैररेखीय ऑपरेटर को एक विविक्त ऑपरेटर $A_{u_h}^h$ द्वारा प्रतिस्थापित किया जाता है। WG योजना एक विविक्त आइगेनपेयर $(u_h, \lambda_h)$ की खोज करती है जो विविक्त ऊर्जा फलन $E_h$ से संबंधित विविक्त वेरियेशनल समीकरण को संतुष्ट करती है।
5. पोस्ट-प्रोसेसिंग (Post-Processing): आइगेनवैल्यू के निचले स्तर (lower bounds) प्राप्त करने के लिए, लेखक एक पोस्ट-प्रोसेसिंग तकनीक का उपयोग करते हैं जो WG सन्निकटन को अनुरूप लैग्रेंज तत्व सन्निकटन (conforming Lagrange element approximations - रैखिक और द्विघात) के साथ जोड़ती है।

मुख्य योगदान और सैद्धांतिक परिणाम
यह शोध पत्र एक कठोर त्रुटि विश्लेषण (error analysis) प्रदान करता है और WG सन्निकटन के लिए विशिष्ट सीमाएँ स्थापित करता है:

• ऊर्जा के लिए निचला स्तर (Lower Bound for Energy): लेखक यह सिद्ध करते हैं कि WG योजना द्वारा उत्पादित विविक्त ऊर्जा $E_h(u_h)$, सटीक ग्राउंड स्टेट ऊर्जा $E(u)$ के लिए एक निचला स्तर (lower bound) है, अर्थात $E(u) - E_h(u_h) \ge 0$।
• त्रुटि अनुमान (Error Estimates): गैररेखीयता $F$ (मान्यताएं A1–A5) और ग्राउंड स्टेट $u \in H^{k+1}(\Omega)$ की नियमितता पर विशिष्ट मान्यताओं के तहत, यह शोध पत्र पर्याप्त छोटे मेश आकार $h$ के लिए निम्नलिखित अभिसरण दरें (convergence rates) प्राप्त करता है:
  • $L^2$-नॉर्म त्रुटि: $\|u - u_h\| \lesssim h^{k+1-2\varepsilon}$.
  • विविक्त $H^1$-नॉर्म त्रुटि: $\|u - u_h\|_{1,h} \lesssim h^{k-\varepsilon}$.
  • आइगेनवैल्यू त्रुटि: $|\lambda - \lambda_h| \lesssim h^{k-\varepsilon}$.
    यहाँ, $\varepsilon \in (0, 1)$ एक छोटा स्थिरांक है जो स्थिरीकरण पद (stabilization term) से उत्पन्न होता है।
• आइगेनवैल्यू के लिए निचला स्तर: अनुरूप तत्वों के सन्निकटन के साथ WG सन्निकटन को जोड़कर, लेखक एक सुधारा गया आइगेनवैल्यू $\tilde{\lambda}_h$ निर्मित करते हैं जो एक निचला स्तर गुण संतुष्ट करता है: $0 \le \lambda - \tilde{\lambda}_h \lesssim h^{2-2\varepsilon}$ (के लिए $k=1$)।
• आइगेनवैल्यू के लिए ऊपरी स्तर (Upper Bound for Eigenvalue): इसी प्रकार, रैखिक अनुरूप तत्वों का उपयोग करके आइगेनवैल्यू के लिए एक ऊपरी स्तर व्युत्पन्न किया गया है, जो दिखाता है कि $0 \le \tilde{\lambda}_h - \lambda \lesssim h^2$।

संख्यात्मक परिणाम (Numerical Results)
सैद्धांतिक निष्कर्षों को विभिन्न विभव $V(x,y)$ और गैररेखीयताओं $f(u^2)$ के साथ यूनिट स्क्वायर डोमेन पर संख्यात्मक प्रयोगों के माध्यम से सत्यापित किया गया है।

• अभिसरण दरें (Convergence Rates): प्रयोग सैद्धांतिक अभिसरण क्रमों की पुष्टि करते हैं। $k=1$ के लिए, देखा गया क्रम दोनों आइगेनवैल्यू और ऊर्जा के लिए लगभग 2.0 है। $k=2$ के लिए, देखा गया क्रम लगभग 4.0 है।
• निचले स्तर का सत्यापन: संख्यात्मक डेटा लगातार दिखाता है कि गणना की गई WG ऊर्जा और पोस्ट-प्रोसेस्ड आइगेनवैल्यू, संदर्भ मानों (जो अत्यधिक परिष्कृत मेश द्वारा अनुमानित हैं) से कम हैं, जो सैद्धांतिक निचले स्तर के दावों का समर्थन करते हैं।
• मजबूती (Robustness): विधि गैररेखीयता पैरामीटर $\beta$ और विभिन्न विभव विन्यासों के विभिन्न मानों के बीच स्थिरता प्रदर्शित करती है।

महत्व
यह शोध पत्र यह दावा करता है कि वीक गैलेरिक विधि न केवल गैररेखीय आइगेनवैल्यू समस्याओं के लिए सटीक सन्निकटन प्रदान करती है, बल्कि स्वाभाविक रूप से ग्राउंड स्टेट ऊर्जा के लिए निचले स्तर (lower bounds) भी उत्पन्न करती है। यह गुण भौतिकी और इंजीनियरिंग अनुप्रयोगों में विशेष रूप से मूल्यवान है जहाँ ऊर्जा स्तरों को ओवरएस्टिमेट (overestimate) करने की तुलना में अंडरएस्टिमेट (underestimate) करना अक्सर अधिक सुरक्षित या वांछनीय होता है। इसके अलावा, आइगेनवैल्यू के निचले और ऊपरी स्तर प्राप्त करने के लिए पोस्ट-प्रोसेसिंग तकनीक विकसित करना, सटीक समाधान की आवश्यकता के बिना त्रुटि अनुमान और समाधान सत्यापन के लिए एक व्यावहारिक उपकरण प्रदान करता है। यह कार्य कठोर त्रुटि विश्लेषण के साथ रैखिक से गैररेखीय स्पेक्ट्रल समस्याओं तक WG विधियों की प्रयोज्यता का विस्तार करता है।