A Variable-Step IMEX-BDF2 Finite Element Method for the Navier–Stokes Equations: Optimal $L^2$ Error Estimates for Velocity and Pressure

तकनीकी सारांश: नेवियर-स्टोक्स समीकरणों के लिए एक परिवर्तनीय-चरण (Variable-Step) IMEX-BDF2 परिमित तत्व विधि (Finite Element Method)

समस्या विवरण
यह शोध पत्र गैर-स्थिर (non-stationary) असंपीड्य नेवियर-स्टोक्स (NS) समीकरणों के संख्यात्मक अनुकरण (numerical simulation) को संबोधित करता है, जो विशेष रूप से उच्च रेनॉल्ड्स संख्या पर अंतर्निहित बहु-पैमाने (multiscale) व्यवहार प्रदर्शित करते हैं। जबकि अनुकूलन योग्य समय-चरण (adaptive time-stepping) को ऐसे प्रवाहों में सटीकता और दक्षता को संतुलित करने के लिए एक मान्यता प्राप्त उपकरण के रूप में उपयोग किया जाता है, द्वितीय-क्रम परिवर्तनीय-चरण विधियों (विशेष रूप से 'बैकवर्ड डिफरेंस फॉर्मूला' या BDF2) के लिए मौजूदा सैद्धांतिक ढांचे में महत्वपूर्ण कमियां मौजूद हैं। परिवर्तनीय-चरण BDF2 योजनाओं पर पिछले अध्ययन या तो केवल प्रथम-क्रम अभिसरण (convergence) स्थापित कर पाए, या गैर-इष्टतम नॉर्म्स (जैसे G-फंक्शनल नॉर्म्स) में त्रुटि अनुमान प्रदान किए, या $L^2$ में इष्टतम अस्थायी दरें प्राप्त कीं लेकिन स्थानिक दरें उप-इष्टतम (suboptimal) रहीं। इसके अतिरिक्त, परिवर्तनीय-चरण सेटिंग में दबाव क्षेत्र (pressure field) के लिए एक पूर्ण, इष्टतम $L^2$ त्रुटि अनुमान साहित्य में अनुपलब्ध था।

कार्यप्रणाली (Methodology)
लेखक एक पूर्णतः विविक्त (fully discrete), परिवर्तनीय-चरण इम्प्लिसिट-एक्सप्लिसिट (IMEX) BDF2 परिमित तत्व विधि (FEM) का प्रस्ताव करते हैं, जिसमें टेलर-हुड (Taylor–Hood) तत्वों ($P_2-P_1$) का उपयोग किया गया है। यह योजना स्थिरता सुनिश्चित करने के लिए संवहन पद (convection term) के लिए एक स्क्यू-सिमेट्रिक (skew-symmetric) रूप का उपयोग करती है। मुख्य पद्धतिगत घटक निम्नलिखित हैं:

1. परिवर्तनीय समय-चरण (Variable Time-Stepping): विधि एक गैर-समान समय विभाजन का उपयोग करती है जिसमें चरण आकार $\tau_k$ और आसन्न चरण अनुपात $r_k = \tau_k/\tau_{k-1}$ होते हैं। एक स्थिरता स्थिति लागू की जाती है जहाँ $0 < r_k \le r_{max} < 4.864$।
2. आरंभीकरण (Initialization): वैश्विक द्वितीय-क्रम सटीकता को बनाए रखने के लिए, पहले समय-चरण की गणना क्रैंक-निकोल्सन (Crank–Nicolson) योजना का उपयोग करके की जाती है, जबकि बाद के चरणों में परिवर्तनीय-चरण BDF2 सूत्र का उपयोग किया जाता है।
3. डिस्क्रीट ऑर्थोगोनल कॉन्वोल्यूशन (DOC) कर्नेल: विश्लेषण गैर-समान समय चरणों को संभालने के लिए DOC कर्नेल का लाभ उठाता है। यह ढांचा प्रारंभिक चरण के प्रभाव को सटीक रूप से मापने और त्रुटि विश्लेषण में समय-चरण के परिवर्तनों के प्रबंधन की अनुमति देता है।
4. त्रुटि विश्लेषण रणनीति:
   • $H^1$ अनुमान: वेग के लिए एक इष्टतम $H^1$ त्रुटि अनुमान गणितीय आगमन (mathematical induction) के माध्यम से स्थापित किया गया है, जिसमें स्पष्ट स्थिरांक (constants) प्राप्त करने के लिए DOC ढांचे का उपयोग किया गया है।
   • $L^2$ अनुमान (वेग): $H^1$ अनुमान को बिना सटीकता खोए इष्टतम $L^2$ अनुमान में अपग्रेड करने के लिए, लेखक एक नेगेटिव-नॉर्म (negative-norm) तर्क का उपयोग करते हैं। यह तकनीक गैर-रेखीय समस्या पर सीधे ऑबिन-निट्शे द्वैतता (Aubin–Nitsche duality) तर्क को लागू करने की जटिलताओं से बचती है।
   • दबाव अनुमान: नेगेटिव-नॉर्म तकनीक और डिस्क्रीट इन्फ-सप (inf-sup) स्थिति का उपयोग करते हुए दबाव के लिए एक इष्टतम $L^2$ त्रुटि अनुमान भी साथ ही प्राप्त किया गया है।

प्रमुख योगदान
इस कार्य के प्राथमिक योगदान सैद्धांतिक हैं, जो परिवर्तनीय-चरण BDF2 विधियों के विश्लेषण में विशिष्ट कमियों को भरते हैं:

• इष्टतम $L^2$ त्रुटि अनुमान: यह शोध पत्र वेग के लिए $L^2$ नॉर्म में स्थान और समय दोनों में पूर्णतः इष्टतम ए प्रियोरी (a priori) त्रुटि अनुमान स्थापित करता है। यह पिछले परिवर्तनीय-चरण विश्लेषणों की तुलना में एक महत्वपूर्ण प्रगति है जो उप-इष्टतम स्थानिक दरों या गैर-$L^2$ नॉर्म्स तक सीमित थे।
• दबाव अभिसरण: यह कार्य परिवर्तनीय-चरण BDF2 सेटिंग में दबाव क्षेत्र के लिए पहला इष्टतम $L^2$ त्रुटि अनुमान प्रदान करता है, जो कि पहले साहित्य में उपलब्ध नहीं था।
• आरंभीकरण विश्लेषण: लेखक कठोरता से प्रदर्शित करते हैं कि प्रारंभिक चरण के लिए क्रैंक-निकोल्सन योजना का उपयोग करने से परिवर्तनीय-चरण BDF2 विधि की समग्र द्वितीय-क्रम अस्थायी सटीकता बनी रहती है।
• स्पष्ट स्थिरांक: विश्लेषण मेष आकार (mesh size) और समय-चरण बाधाओं के लिए स्पष्ट ऊपरी सीमाएं प्रदान करता है, जो विशिष्ट नियमितता मान्यताओं के तहत त्रुटि अनुमानों की वैधता सुनिश्चित करता है।

परिणाम
सैद्धांतिक निष्कर्षों का समर्थन FreeFEM++ का उपयोग करके किए गए संख्यात्मक प्रयोगों द्वारा किया गया है:

• सटीकता सत्यापन: सटीक समाधानों के साथ 2D लैटिस-वोर्टेक्स समस्या पर परीक्षण अनुमानित अभिसरण दरों की पुष्टि करते हैं। यह योजना ग्रेडेड और रैंडम टाइम ग्रिड दोनों पर वेग और दबाव के लिए द्वितीय-क्रम अभिसरण प्राप्त करती है।
• स्थानिक अभिसरण (Spatial Convergence): स्थिर समय चरणों के साथ, विधि वेग के लिए $L^2$ नॉर्म में तृतीय-क्रम अभिसरण और $H^1$ नॉर्म में द्वितीय-क्रम अभिसरण प्रदर्शित करती है, तथा दबाव के लिए $L^2$ नॉर्म में द्वितीय-क्रम अभिसरण दिखाती है।
• मजबूती (Robustness): योजना सैद्धांतिक सीमा 4.864 से अधिक के अनुपात होने पर भी अपेक्षित सटीकता बनाए रखती है (7 के अनुपात तक परीक्षण किया गया), जो यह सुझाव देता है कि सैद्धांतिक स्थिति व्यावहारिक रूप से स्थिरता के लिए पर्याप्त है लेकिन अनिवार्य नहीं है।
• अनुकूलन योग्य अनुप्रयोग (Adaptive Application): $Re=1000$ पर लिड-ड्रिवन कैविटी फ्लो सिमुलेशन एक अनुकूलन योग्य समय-चरण रणनीति की प्रभावशीलता को प्रदर्शित करता है। विधि सफलतापूर्वक छोटे चरणों के साथ तीव्र प्रवाह संक्रमणों को हल करती है और जैसे ही प्रवाह स्थिर अवस्था की ओर बढ़ता है, चरण आकार को बढ़ा देती है, जिससे सुचारू चरण-आकार विकास बना रहता है।

महत्व
यह शोध पत्र नेवियर-स्टोक्स समीकरणों के लिए परिवर्तनीय-चरण BDF2 विधियों के सिद्धांत में एक उल्लेखनीय अंतर को भरने का दावा करता है। योजना के लिए—जिसमें दबाव क्षेत्र भी शामिल है—एक पूर्ण और इष्टतम त्रुटि चित्रण (error portrait) प्रदान करके, यह इस व्यापक रूप से उपयोग किए जाने वाले परिवर्तनीय-चरण दृष्टिकोण के सैद्धांतिक आधार को सुदृढ़ करता है। यह कार्य प्रमाणित करता है कि परिवर्तनीय-चरण IMEX-BDF2-FEM का विश्लेषण इष्टतम सटीकता प्राप्त करने के लिए कठोरता से किया जा सकता है, जिससे यह बहु-पैमाने वाले प्रवाह घटनाओं के अनुकरण के लिए एक विश्वसनीय उपकरण बन जाता है जहाँ अनुकूलन योग्य रिज़ॉल्यूशन महत्वपूर्ण है। लेखक उल्लेख करते हैं कि ये परिणाम दीर्घकालिक स्थिरता अनुमानों और युग्मित मल्टीफिजिक्स सिस्टम (coupled multiphysics systems) के लिए भविष्य के विस्तार का आधार रखते हैं।