A dichotomy of finite element spaces and its application to an energy-conservative scheme for the regularized long wave equation

तकनीकी सारांश: परिमित तत्व स्थानों (Finite Element Spaces) का एक द्वैत और नियमितीकृत दीर्घ तरंग समीकरण (Regularized Long Wave Equation) के लिए एक ऊर्जा-संरक्षण योजना में इसका अनुप्रयोग

समस्या विवरण
यह शोधपत्र गैररेखीय विसरणी तरंग समीकरणों (nonlinear dispersive wave equations) के ऊर्जा-संरक्षण गैलरकिन विविक्तीकरणों (Galerkin discretizations) में देखे गए एक विशिष्ट घटना को संबोधित करता है: परिमित तत्व स्थान के बहुपद क्रम (polynomial degree) के आधार पर अभिसरण व्यवहार (convergence behavior) में एक असामान्य द्वैत। अनुकूल अभिसरण दरें तब प्राप्त की जाती हैं जब विषम बहुपद डिग्री वाले निरंतर लैग्रेंज परिमित तत्व स्थानों का उपयोग किया जाता है, जबकि सम-डिग्री वाले बहुपद कम अभिसरण दरें प्रदर्शित करते हैं। लेखक यह प्रदर्शित करने का लक्ष्य रखते हैं कि यह व्यवहार परिमित तत्व स्थानों की अपनी संरचना के प्रति अंतर्निहित है, न कि किसी विशिष्ट संख्यात्मक योजनाओं का कोई कृत्रिम प्रभाव। इसके अतिरिक्त, यह शोधपत्र नियमितीकृत दीर्घ तरंग (RLW) समीकरण (जिसे बेंजामिन-बोना-माहॉ या BBM समीकरण भी कहा जाता है) के लिए इस गुण के निहितार्थों की जांच करता है, जो द्रव्यमान (mass), आवेग (impulse), और ऊर्जा संरक्षण नियमों को धारण करने वाला एक मौलिक मॉडल है।

कार्यप्रणाली
विश्लेषण एक समान आवधिक ग्रिड (uniform periodic grid) पर परिभाषित एक-आयामी निरंतर लैग्रेंज परिमित तत्व स्थानों ($V^k_h$) के ढांचे के भीतर किया गया है। कार्यप्रणाली का मुख्य हिस्सा परिमित तत्व स्थान $V^k_h = S^1_h \oplus S^2_h$ का प्रत्यक्ष-योग अपघटन (direct-sum decomposition) है, जहाँ $S^1_h$ शीर्ष-आधारित फलनों (vertex-based functions) से बना है और $S^2_h$ बबल फलनों (bubble functions) से बना है।

1. $L^2$-प्रक्षेप (Projections) का विश्लेषण: लेखक अवकलजों के $L^2$-प्रक्षेप $P$ का गहन विश्लेषण करते हैं। वे सिद्ध करते हैं कि पर्याप्त सुचारू आवधिक फलनों के लिए, त्रुटि पद $\|P[(Pu - u)_x]\|$ विशेष रूप से विषम बहुपद डिग्री ($k = 2\kappa + 1$) के लिए एक "सुपर-एप्रोक्सिमेशन" (super-approximation) गुण प्रदर्शित करता है, जो $O(h^{k+1})$ की कोटि प्राप्त करता है। इसके विपरीत, सम डिग्री ($k = 2\kappa$) के लिए, त्रुटि मानक अनुकूल कोटि $O(h^k)$ पर बनी रहती है। यह अंतर ऑर्थोगोनल बहुपदों और विशिष्ट आधार फलनों (जैसे $\psi$ और $\phi_j$) के गुणों का उपयोग करते हुए शीर्ष और बबल उप-स्थानों पर प्रक्षेपों के विस्तृत अनुमानों के माध्यम से प्राप्त किया गया है।
2. अर्ध-विविक्त सूत्रीकरण (Semi-Discrete Formulation): RLW समीकरण के लिए एक संशोधित गैलरकिन परिमित तत्व विधि प्रस्तावित की गई है। इस योजना को द्रव्यमान और एक त्रिघात ऊर्जा फलन (cubic energy functional) दोनों को संरक्षित करने के लिए डिज़ाइन किया गया है। अर्ध-विविक्त सूत्रीकरण को गैररेखीय और विसरणी पदों के $L^2$-प्रक्षेप को शामिल करते हुए साधारण अवकल समीकरणों (ODEs) की एक प्रणाली के रूप में व्यक्त किया गया है।
3. त्रुटि विश्लेषण: प्रस्तावित अर्ध-विविक्त सूत्रीकरण के लिए सैद्धांतिक 'अ प्रायोरी' (a priori) त्रुटि सीमाएं स्थापित की गई हैं। प्रमाण खंड 2 में स्थापित अवकलजों के $L^2$-प्रक्षेप की द्वैत संपत्ति पर निर्भर करता है। विश्लेषण $u$ और उसके व्युत्पन्न $u_x$ के लिए विशिष्ट अभिसरण दरों को प्राप्त करने हेतु विषम और सम $k$ के बीच अंतर करता है।
4. आवेग सन्निकटन (Impulse Approximation): जबकि प्रस्तावित योजना द्रव्यमान और ऊर्जा को संरक्षित करती है, यह आवेग फलन को सख्ती से संरक्षित नहीं करती है। लेखक आवेग सन्निकटन की सटीकता का विश्लेषण करते हैं, यह दिखाते हुए कि यह बहुपद डिग्री पर निर्भर उच्च सटीकता के साथ संरक्षित रहता है।
5. पूर्णतः विविक्त योजना (Fully Discrete Scheme): समय में ऊर्जा फलन को संरक्षित करते हुए इस अर्ध-विविक्त प्रणाली को एकीकृत करने के लिए, लेखक रिलैक्सेशन रनगे-कुट्टा (Relaxation Runge–Kutta) विधियों का उपयोग करते हैं। ये विधियाँ सटीक विविक्त ऊर्जा के संरक्षण को सुनिश्चित करने के लिए एक रिलैक्सेशन पैरामीटर $\gamma_n$ के माध्यम से समय चरण (time step) को समायोजित करती हैं।

प्रमुख योगदान और परिणाम

• सैद्धांतिक द्वैत: प्राथमिक योगदान यह प्रमाण है कि $L^2$-प्रक्षेप के अवकलजों का सुपर-एप्रोक्सिमेशन गुण विशेष रूप से विषम-डिग्री लैग्रेंज परिमित तत्व स्थानों के लिए ही विशिष्ट है। विषम $k$ के लिए, $\|P[(Pu - u)_x]\| = O(h^{k+1})$; सम $k$ के लिए, यह $O(h^k)$ है।
• RLW के लिए अभिसरण दरें: इस गुण को RLW समीकरण के लिए ऊर्जा-संरक्षण योजना पर लागू करते हुए, लेखक स्थापित करते हैं कि:
  • विषम $k$ के लिए: संख्यात्मक समाधान अनुकूल कोटि $O(h^{k+1})$ के साथ अभिसरित होता है, और व्युत्पन्न सन्निकटन $u_x$ भी $O(h^{k+1})$ की कोटि के साथ अभिसरित होता है (सुपर-कन्वर्जेंस)।
  • सम $k$ के लिए: अभिसरण उप-अनुकूल (suboptimal) है, जिसमें समाधान $O(h^k)$ पर और व्युत्पन्न $O(h^{k-1})$ पर अभिसरित होता है।
• आवेग संरक्षण: शोधपत्र यह प्रदर्शित करता है कि हालांकि संशोधित योजना आवेग फलन को सख्ती से संरक्षित नहीं करती है, फिर भी आवेग सन्निकटन में त्रुटि विषम $k$ के लिए $O(h^{2k})$ और सम $k$ के लिए $O(h^{2k-2})$ के क्रम की है।
• संख्यात्मक सत्यापन: व्यापक संख्यात्मक प्रयोग सैद्धांतिक भविष्यवाणियों की पुष्टि करते हैं। त्रुटि $\|P[(Pu - u)_x]\|$ और RLW समाधान के लिए अभिसरण दरें प्राप्त सीमाओं से मेल खाती हैं। इसके अतिरिक्त, सिमुलेशन दिखाते हैं कि द्रव्यमान और ऊर्जा मशीन परिशुद्धता (machine precision) तक संरक्षित हैं, जबकि आवेग त्रुटियां छोटी रहती हैं और अनुमानित अभिसरण पैटर्न का पालन करती हैं।

महत्व
शोधपत्र का दावा है कि ऊर्जा-संरक्षण योजनाओं के लिए अभिसरण व्यवहार में देखा गया यह द्वैत एक संख्यात्मक विसंगति नहीं है, बल्कि उपयोग किए जाने वाले परिमित तत्व स्थानों का एक मौलिक गुण है। $L^2$-प्रक्षेप के अवकलजों के माध्यम से बहुपद डिग्री की समता (parity) को सीधे अभिसरण दर से जोड़कर, यह कार्य पूर्ववर्ती अनुभवजन्य अवलोकनों (जैसे [28]) के लिए एक सैद्धांतिक स्पष्टीकरण प्रदान करता है। परिणाम सुझाव देते हैं कि त्रिघात ऊर्जा फलनों और व्युत्पन्न प्रक्षेपों से जुड़ी समस्याओं के लिए, विषम-डिग्री बहुपदों का चयन केवल एक अनुमानित विधि (heuristic) नहीं बल्कि इष्टतम सटीकता प्राप्त करने के लिए एक आवश्यकता है। RLW समीकरण के लिए प्रस्तावित ऊर्जा-संरक्षण योजना एक ठोस अनुप्रयोग के रूप में कार्य करती है, जो यह प्रदर्शित करती है कि यदि उपयुक्त परिमित तत्व स्थान का चयन किया जाए, तो उच्च-क्रम सटीकता और मौलिक भौतिक मात्राओं (द्रव्यमान और ऊर्जा) के संरक्षण को एक साथ प्राप्त किया जा सकता है।