A semi-Lagrangian method for the polyatomic ESBGK model

तकनीकी सारांश: पॉलीएटमिक ESBGK मॉडल के लिए एक सेमी-लैग्रेंजियन विधि

समस्या विवरण
पॉलीएटमिक गैसों के लिए विरल गैस गतिकी (rarefied gas dynamics) का सटीक अनुकरण करने हेतु ऐसे काइनेटिक मॉडलों की आवश्यकता होती है जो आंतरिक ऊर्जा मोड (घूर्णी और कंपन) को ध्यान में रखते हैं, जो परिवहन गुणांकों (transport coefficients), विश्रांति प्रक्रियाओं (relaxation processes) और गैर-संतुलन प्रभावों को महत्वपूर्ण रूप से प्रभावित करते हैं। जबकि बोल्ट्ज़मैन समीकरण एक मौलिक विवरण प्रदान करता है, इसका कोलिजन ऑपरेटर पांच-गुने समाकलन (five-fold integral) के कारण गणनात्मक रूप से अत्यधिक कठिन है। विश्रांति-प्रकार के मॉडल, जैसे कि BGK और उसका विस्तार, एलिप्सॉइडल सांख्यिकीय ESBGK, एक अधिक सुलभ विकल्प प्रदान करते हैं, जो कोलिजन ऑपरेटर को स्थानीय संतुलन वितरण की ओर एक विश्रांति पद (relaxation term) द्वारा प्रतिस्थापित करते हैं। हालाँकि, मानक BGK मॉडल सही प्रैंड्टल संख्या (Prandtl number) को पुनरुत्पादित करने में विफल रहते हैं, और मौजूदा संख्यात्मक विधियाँ अक्सर निरंतरता सीमा (कम नूसेन संख्या/Knudsen numbers) में विश्रांति पद की कठोरता (stiffness) से जूझती हैं या उनके लिए प्रतिबंधात्मक समय-चरण बाधाएं (time-step constraints) आवश्यक होती हैं। इसके अतिरिक्त, मजबूत सीमा स्थितियों (boundary conditions) और एसिम्प्टोटिक प्रिजर्वेशन गुणों वाले पॉलीएटमिक ESBGK मॉडलों के लिए नियतात्मक (deterministic) विधियाँ सीमित हैं।

कार्यप्रणाली
लेखक पॉलीएटमिक ESBGK मॉडल के लिए एक नई नियतात्मक, ग्रिड-आधारित सेमी-लैग्रेंजियन योजना प्रस्तावित करते हैं। कार्यप्रणाली निम्नलिखित घटकों को एकीकृत करती है:

1. सेमी-लैग्रेंजियन ढांचा: परिवहन पद (transport term) को अभिलक्षणों की विधि (method of characteristics) का उपयोग करके उपचारित किया जाता है। वर्तमान समय चरण से पिछले चरण तक अभिलक्षणों को पीछे ट्रेस करके, यह योजना विशिष्ट कॉउरेंट-फ्रीड्रिक्स-लेवी (CFL) स्थिति से बचती है जो स्पष्ट यूलेरियन परिवहन (explicit Eulerian transport) से जुड़ी होती है, जिससे बड़े समय चरणों की अनुमति मिलती है।
2. विश्रांति का निहित (Implicit) उपचार: कठोर विश्रांति पद को एक निहित ए-स्टेबल लीनियर मल्टीस्टेप विधि (विशेष रूप से दूसरी श्रेणी की सटीकता के लिए BDF2) का उपयोग करके संभाला जाता है। मानक निहित योजनाओं के विपरीत, जिन्हें बड़े गैर-रेखीय सिस्टम को हल करने की आवश्यकता होती है, BGK ऑपरेटर की विशिष्ट संरचना लेखकों को समस्या को पुनर्गठित करने की अनुमति देती है।
3. स्पष्ट पुनर्गठन (Explicit Reformulation): स्थानीय संतुलन वितरण के उच्च-क्रम मोमेंट्स (विशेष रूप से ट्रांसलेशनल और इंटरनल तापमान) को पिछले समय चरण से ज्ञात मात्राओं के फलन के रूप में व्यक्त करके, लेखक एक ऐसी प्रणाली प्राप्त करते हैं जहाँ निहित अपडेट गणनात्मक रूप से सस्ता हो जाता है।
   • द्रव्यमान और संवेग को संरक्षण नियमों के माध्यम से स्पष्ट रूप से अपडेट किया जाता है।
   • ट्रांसलेशनल ($T_{tr}$) और इंटरनल ($T_{int}$) तापमान को प्रत्येक ग्रिड बिंदु पर एक छोटे, दो-आयामी गैर-रेखीय सिस्टम (न्यूटन के मेथड द्वारा हल किया गया) के माध्यम से हल किया जाता है।
   • वेग सहप्रसरण मैट्रिक्स (velocity covariance matrix) को स्पष्ट रूप से प्राप्त किया जाता है।
4. आयामी कमी (Chu-Reduction): 3D वेग और आंतरिक ऊर्जा प्रभावों को बनाए रखते हुए 2D स्थानिक सिमुलेशन को सुगम बनाने के लिए, लेखक एक चू-रिडक्शन (Chu-reduction) का उपयोग करते हैं। यह वितरण फलन की लंबवत वेग घटक और आंतरिक ऊर्जा चर पर निर्भरता को कम करता है, जिसके परिणामस्वरूप तीन रिड्यूस्ड वितरण फलन ($f_1, f_2, f_3$) प्राप्त होते हैं जो युग्मित BGK-प्रकार के समीकरणों का पालन करते हैं।
5. सीमा स्थितियाँ (Boundary Conditions): शोध पत्र नियतात्मक सॉल्वर के लिए विशिष्ट सीमा उपचार प्रस्तुत करता है:
   • डिफ्यूसिव-रिफ्लेक्टिव वॉल्स: पुनर्गठन चरण (कोलिजन-पूर्व) के ठीक बाद सीधे लागू किया जाता है, जो द्रव्यमान संरक्षण सुनिश्चित करता है।
   • इनफ्लो/आउटफ्लो: यूलर समीकरण एक्सट्रपोलेशन से प्राप्त सबसोनिक दबाव-चालित स्थितियाँ, जो कथित तौर पर नियतात्मक BGK सॉल्वरों के लिए पहले रिपोर्ट नहीं की गई थीं।

प्रमुख योगदान
शोध पत्र पाँच प्राथमिक योगदानों को रेखांकित करता है:

1. नया विवेकीकरण (Discretization): पॉलीएटमिक ESBGK मॉडल का एक सेमी-लैग्रेंजियन विवेकीकरण जिसमें कोलिजन ऑपरेटर का एक स्पष्ट उपचार शामिल है। इसमें दो तापमानों के लिए एक गैर-रेखीय सिस्टम को हल करना और एक स्पष्ट वेग सहप्रसरण मैट्रिक्स व्युत्पन्न करना शामिल है, जो पिछले मोनोएटमिक दृष्टिकोणों का सामान्यीकरण करता है।
2. अभिसरण प्रमाण (Convergence Proof): एक कठोर प्रमाण जो यह दर्शाता है कि प्रथम-क्रम की योजना कंप्रेसिबल नेवियर-स्टोक्स-फूरियर समीकरणों के एक सुसंगत विवेकीकरण की एसिम्प्टोटिक रूप से अभिसरण करती है, जो सही परिवहन गुणांकों (श्यानता और ऊष्मा चालन) को प्राप्त करती है और मोनोएटमिक ESBGK के पिछले प्रमाणों का सामान्यीकरण करती है।
3. सीमा स्थिति कार्यान्वयन: डिफ्यूज़-रिफ्लेक्टिव और प्रेशर-ड्रिवन इनफ्लो/आउटफ्लो सीमा स्थितियों का एक अभिनव, कुशल कार्यान्वयन जो पुनरावृत्ति (iterative) सॉल्वर या सीमाओं तक एक्सट्रपोलेशन से बचता है।
4. व्यापक बेंचमार्किंग: फूरियर और कुएट प्रवाह परीक्षण मामलों के लिए BGK, ESBGK, पॉलीएटमिक ESBGK, DSMC और नेवियर-स्टोक्स मॉडलों का एक व्यवस्थित तुलनात्मक अध्ययन, जो विभिन्न प्रवाह व्यवस्थाओं में मॉडल की निष्ठा का आकलन करने के लिए एक सुसंगत संदर्भ ढांचा स्थापित करता है।
5. जटिल ज्यामिति सिमुलेशन: गैर-तुच्छ ज्यामिति और समय-निर्भर सीमाओं में मजबूती प्रदर्शित करते हुए, स्थिर और गतिशील ऑरिफिस फ्लो कॉन्फ़िगरेशन के सफल अनुप्रयोग।

संख्यात्मक परिणाम
लेखक अपने तरीके को कई परीक्षण मामलों के माध्यम से मान्य करते हैं:

• अभिसरण: योजना $L_1$ त्रुटि नॉर्म्स में घोषित सटीकता क्रम (प्रथम और द्वितीय क्रम) प्राप्त करती है।
• फूरियर और कुएट प्रवाह: फ्री-मॉलिकुलर, ट्रांजिशन, स्लिप और कंटिनम व्यवस्थाओं में सिमुलेशन दिखाते हैं कि:
  • सभी मॉडल (BGK, ESBGK, पॉलीएटमिक ESKG) वेग प्रोफाइल के लिए कंटिनम और फ्री-मॉलिकुलर व्यवस्थाओं में सैद्धांतिक सीमाओं से मेल खाते हैं।
  • पॉलीएटमिक ESKG मॉडल तापमान प्रोफाइल और हीट फ्लक्स को कैप्चर करने में मानक BGK मॉडल की तुलना में काफी बेहतर प्रदर्शन करता है, जो सही प्रैंड्टल संख्या और आंतरिक ऊर्जा विश्रांति को ध्यान में रखते हुए DSMC परिणामों के साथ सही ढंग से मेल खाता है।
  • मानक BGK मॉडल कंटिनम सीमा में निश्चित प्रैंड्टल संख्या 1 के कारण शियर स्ट्रेस और ऊष्मा चालन को एक साथ कैप्चर करने में विफल रहता है।
• ऑरिफिस फ्लो:
  • कंटिनम व्यवस्था ($Kn \approx 4 \times 10^{-3}$) में, पॉलीएटमिक ESKG परिणाम दबाव और वेग के लिए नेवियर-स्टोक्स (Star-CCM+) डेटा के करीब आते हैं, जो तापमान समझौते के मामले में BGK वेरिएंट्स में सर्वश्रेष्ठ है।
  • ट्रांजिशन व्यवस्था ($Kn \approx 9 \times 10^{-2}$) में, विधि DSMC सिमुलेशन के साथ अच्छी तरह से सहमत होती है, जहाँ नेवियर-स्टोक्स अमान्य है।
  • विभिन्न दबाव अनुपातों के लिए मास फ्लो रेट DSMC और नेवियर-स्टोक्स संदर्भ डेटा दोनों के अनुरूप हैं।
• मूविंग ऑरिफिस: विधि एक समय-निर्भर ऑरिफिस ऊंचाई का सफलतापूर्वक अनुकरण करती है, जो प्रवाह विस्तार (flow expansion), पुनरावृत्ति (recirculation) और सुपरसोनिक विस्तार क्षेत्र में शॉक डायमंड्स जैसे जटिल प्रवाह फीचर्स को कैप्चर करती है।

महत्व और दावे
शोध पत्र का दावा है कि प्रस्तावित विधि पॉलीएटमिक काइनेटिक मॉडलों के अनुकरण के लिए एक अत्यधिक कुशल और संख्यात्मक रूप से स्थिर नियतात्मक ढांचा प्रदान करती है। सेमी-लैग्रेजियन दृष्टिकोण को कोलिजन टर्म के एक निहित, फिर भी गणनात्मक रूप से सस्ते उपचार के साथ जोड़कर, यह विधि काइनेटिक ट्रांसपोर्ट की प्रतिबंधात्मक समय-चरण बाधाओं को समाप्त करती है और कंटिनम सीमा में कठोरता (stiffness) को संभालती है। यह योजना एसिम्प्टोटिक प्रिजर्विंग (जैसे ही $Kn \to 0$ होता है, यूलर समीकरणों में अभिसरण करती है) और स्टिफली एक्यूरेट (सही परिवहन गुणांकों के साथ नेवियर-स्टोक्स समीकरणों को पुनः प्राप्त करती है) है।

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