Improved Arrow-Hurwicz method for Stationary Inductionless Magnetohydrodynamics System

तकनीकी सारांश: स्थिर इंडक्शनलेस मैग्नेटोहाइड्रोडायनामिक्स के लिए उन्नत एरो-हर्वीज़ विधि

समस्या विवरण
यह शोध पत्र स्थिर इंडक्शनलेस मैग्नेटोहाइड्रोडायनामिक्स (IMHD) प्रणाली के संख्यात्मक सन्निकटन (numerical approximation) को संबोधित करता है। यह प्रणाली उन परिदृश्यों का मॉडल करती है जहाँ चालक तरल के भीतर धाराओं द्वारा उत्पन्न चुंबकीय क्षेत्र, बाहरी चुंबकीय क्षेत्र $\mathbf{B}$ की तुलना में नगण्य होता है। नियंत्रित समीकरण वेग $\mathbf{u}$, किनेमैटिक दबाव $p$, धारा घनत्व $\mathbf{J}$, और विद्युत विभव $\phi$ को युग्मित (couple) करते हैं। यह प्रणाली लोरेन्त्ज़ बल पद (Lorentz force term) द्वारा संशोधित नेवियर-स्टोक्स समीकरणों द्वारा अभिलक्षित है, जो ओम और फैराडे के नियमों से प्राप्त धारा घनत्व समीकरण के साथ युग्मित है, साथ ही असंपीड्यता (incompressibility) और आवेश संरक्षण बाधाओं को भी शामिल करती है।

यद्यपि विभिन्न संख्यात्मक दृष्टिकोण मौजूद हैं, जिनमें मोनोलिथिक परिमित तत्व विधियाँ (monolithic finite element methods) और पूर्णतः युग्मित पुनरावृत्ति योजनाएं शामिल हैं, वे अक्सर प्रत्येक पुनरावृत्ति पर बड़े, युग्मित सैडल-पॉइंट रैखिक प्रणालियों को हल करने की आवश्यकता रखते हैं। लेखक उल्लेख करते हैं कि मौजूदा विकेंद्रीकृत (decoupled) विधियाँ, जैसे कि शिया और यांग [28] द्वारा IMHD पर लागू किया गया मानक एरो-हर्वीज़ (AH) क्रम, स्थिरता या अभिसरण दक्षता (convergence efficiency) के संबंध में सीमाओं का सामना कर सकते हैं, विशेष रूप से लैग्ड वेरिएबल्स (lagged variables) से उत्पन्न होने वाले मिश्रित पदों (mixed terms) के प्रबंधन के संबंध में।

कार्यप्रणाली
लेखक एरो-हर्वीज़ विधि पर आधारित एक नवीन, विकेंद्रीकृत पुनरावृत्ति परिमित तत्व योजना प्रस्तावित करते हैं। यह दृष्टिकोण स्थिरता और अभिसरण को बढ़ाने के लिए मानक AH ढांचे में विशिष्ट संशोधनों को पेश करता है:

1. विकेंद्रीकरण रणनीति (Decoupling Strategy): विधि चार क्षेत्र चरों ($\mathbf{u}, p, \mathbf{J}, \phi$) के समाधान को प्रत्येक पुनरावृत्ति चरण में अलग-अलग रैखिक प्रणालियों में विभाजित करती है।
2. संवेग और निरंतरता (Momentum and Continuity): वेग-दबाव उप-प्रणाली [25] से प्रेरित एक AH योजना का उपयोग करती है। इसमें वर्तमान और पिछले पुनरावृत्तियों के बीच वेग लैप्लासियनों के अंतर ($\frac{1}{\rho_1}(\nabla(\mathbf{u}^{n+1}-\mathbf{u}^n), \nabla\mathbf{v})$) से जुड़े एक दंड पद (penalty term) को जोड़ना शामिल है और एक दंड पैरामीटर $\gamma_1$ के माध्यम से असंपीड्यता बाधा को शिथिल (relax) करना शामिल है।
3. धारा घनत्व नवाचार (Current Density Innovation): प्राथमिक पद्धतिगत योगदान धारा घनत्व समीकरण है। क्लासिकल स्ट्रीमलाइन डिफ्यूजन विधि [5] से प्रेरित होकर, लेखक नए दंड पदों को पेश करते हैं जो $\frac{1}{\rho_2}(\mathbf{J}^{n+1}-\mathbf{J}^n, \mathbf{K})_{\mathbf{B}}$ के रूप में होते हैं, जहाँ $(\cdot, \cdot)_{\mathbf{B}}$ चुंबकीय क्षेत्र के क्रॉस उत्पाद के साथ एक क्वाज़ी-इनर प्रोडक्ट (quasi-inner product) है। ये पद गैर-ऋणात्मक ऊर्जा और विपाश्विक (dissipative) पद उत्पन्न करने के लिए डिज़ाइन किए गए हैं, जो प्रभावी रूप से उन प्रतिकूल मिश्रित पदों को नियंत्रित करते हैं जो वेरिएबल्स के लैग होने पर उत्पन्न होते हैं।
4. परिमित तत्व विविक्तीकरण (Finite Element Discretization): यह योजना अनुरूप परिमित तत्व स्थानों का उपयोग करती है: वेग और दबाव के लिए टेलर-हुड तत्व ($P_2/P_1$), और धारा घनत्व एवं विद्युत विभव के लिए क्रमशः रविआट-थॉमस तत्व ($RT_k$) और विच्छिन्न (discontinuous) पीसवाइज लीनियर तत्व ($P_{k,dc}$)। यह विकल्प विच्छिन्न आवेश संरक्षण गुण ($\nabla \cdot \mathbf{J}_h = 0$) सुनिश्चित करता है।

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

• स्थिरता और परिबद्धता (Stability and Boundedness): एक "लघु डेटा" स्थिति के तहत (विशेष रूप से, रेनॉल्ड्स और कपलिंग नंबरों के सापेक्ष स्रोत पदों पर एक सीमा), लेखक सिद्ध करते हैं कि पुनरावृत्ति अनुक्रम समान रूप से परिबद्ध रहता है।
• ज्यामितीय अभिसरण (Geometric Convergence): विश्लेषण प्रदर्शित करता है कि विधि ज्यामितीय (रैखिक) अभिसरण प्राप्त करती है। प्रमाण एक विशिष्ट ऊर्जा नॉर्म (energy norm) के निर्माण और यह दिखाने पर निर्भर करता है कि त्रुटि अनुक्रम एक गुणांक $\lambda \in (0,1)$ के साथ संकुचित होता है।
• त्रुटि अनुमान (Error Estimates): शोध पत्र त्रुटि सीमाओं को व्युत्पन्न करता है जो विविक्तीकरण त्रुटि (मेश आकार $h$ पर निर्भर) और पुनरावृत्ति त्रुटि (पुनरावृत्ति गणना $n$ पर निर्भर) को जोड़ती है।
• विकेंद्रीकरण लाभ (Decoupling Advantage): पूर्णतः युग्मित दृष्टिकोणों के विपरीत, प्रस्तावित विधि सहायक समीकरणों के लिए सममित और धनात्मक-निश्चित (symmetric and positive-definite) गुणांक मैट्रिसेस वाली रैखिक प्रणालियों का परिणाम देती है, जो कुशल समाधान रणनीतियों को सुगम बनाती है।

संख्यात्मक परिणाम
लेखक फ्रीफेम++ (FreeFem++) लाइब्रेरी का उपयोग करके अपनी सैद्धांतिक खोजों को संख्यात्मक प्रयोगों के माध्यम से मान्य करते हैं:

• सुचारू समाधान परीक्षण (Smooth Solution Test): एक यूनिट स्क्वायर पर एक निर्मित सुचारू समाधान के साथ, योजना सभी चरों ($\mathbf{u}, p, \mathbf{J}, \phi$) के लिए अपेक्षित इष्टतम अभिसरण दर प्राप्त करती है और मेश रिफाइनमेंट के बावजूद लगभग समान पुनरावृत्ति संख्या बनाए रखती है।
• विलक्षण समाधान परीक्षण (Singular Solution Test): एक एल-आकार (L-shaped) डोमेन में, जहाँ री-एंट्रेंट कॉर्नर (re-entrant corners) हैं, विधि वेग और धारा घनत्व क्षेत्रों में विलक्षणताओं (singularities) को सफलतापूर्वक कैप्चर करती है। परिणाम दिखाते हैं कि प्रस्तावित योजना [28] की AH विधि की तुलना में तेजी से अभिसरित होती है, विशेष रूप से दंड पैरामीटर $\rho$ के बड़े मानों के लिए, जहाँ [28] की विधि अस्थिर हो जाती है।
• लिड-ड्रिवन कैविटी फ्लो (Lid-Driven Cavity Flows): विधि का परीक्षण स्थिर और गैर-स्थिर चुंबकीय क्षेत्रों के साथ 2D और 3D लिड-ड्रिवन कैविटी समस्याओं के लिए किया गया है। परिणाम उच्च रेनॉल्ड्स और हार्टमैन संख्याओं को संभालने की योजना की क्षमता को प्रदर्शित करते हैं, जो जटिल प्रवाह संरचनाओं को कैप्चर करते हैं। पुनरावृत्ति गणनाएँ उचित रहती हैं, और विधि 3D समस्याओं तक प्रभावी ढंग से स्केल होती है।

महत्व और दावे
शोध पत्र दावा करता है कि प्रस्तावित "उन्नत एरो-हर्वीज़" विधि स्थिर IMHD के लिए मौजूदा पुनरावृत्ति योजनाओं की तुलना में एक महत्वपूर्ण प्रगति प्रदान करती है। प्राथमिक महत्व धारा घनत्व समीकरण में नए दंड पदों की शुरूआत में है, जिसे लेखक तर्क देते हैं कि सटीकता से समझौता किए बिना मिश्रित पदों को नियंत्रित करने और ज्यामितीय अभिसरण सुनिश्चित करने के लिए केंद्रीय है।

लेखक इस बात पर जोर देते हैं कि उनका दृष्टिकोण सभी चार क्षेत्र चरों के समाधान को प्रभावी ढंग से अलग करता है, जिससे एक कम्प्यूटेशनल संरचना प्राप्त होती है जो पूर्णतः युग्मित या पूर्व में अध्ययन की गई विकेंद्रीकृत विधियों से भिन्न और अधिक कुशल है। उनके संख्यात्मक परीक्षण पुष्टि करते हैं कि योजना कुशल और मजबूत है, जो सुचारू और विलक्षण दोनों प्रकार के समाधानों तथा जटिल 3D ज्यामिति को संभालने में सक्षम है, जबकि आवेश संरक्षण के भौतिक गुण को बनाए रखती है। यह कार्य सुझाव देता है कि यह विधि औद्योगिक MHD प्रक्रियाओं, जैसे पंप और जनरेटर के अनुकरण के लिए एक व्यवहार्य विकल्प प्रदान करती है, जहाँ उच्च हार्टमैन संख्या और जटिल सीमा स्थितियाँ प्रचलित हैं।