तकनीकी सारांश: गतिज समीकरणों के लिए तरंग-कण अपघटन (Wave-Particle Decomposition) I: सिद्धांत और संख्यात्मकता

1. समस्या विवरण

काइनेटिक समीकरणों (जैसे बोल्ट्ज़मैन समीकरण) का संख्यात्मक अनुकरण, मल्टीस्केल व्यवस्थाओं (multiscale regimes) में एक मौलिक चुनौती का सामना करता है। प्रत्यक्ष काइनेटिक विवेकीकरण (जैसे DSMC या नियतात्मक सॉल्वर) निरंतरता सीमा (continuum limit) में गणनात्मक रूप से महंगे होते हैं क्योंकि उन्हें सूक्ष्म पैमानों (औसत मुक्त पथ और टकराव समय) को हल करना पड़ता है, भले ही समाधान साम्यावस्था (equilibrium) के करीब हो। इसके विपरीत, शुद्ध स्थूल (macroscopic) सॉल्वर (यूलर या नेवियर-स्टोक्स) विरल व्यवस्थाओं (rarefied regimes) में विफल हो जाते हैं जहाँ गैर-साम्यावस्था प्रभाव हावी होते हैं।

मौजूदा मल्टीस्केल विधियाँ, जैसे कि यूनिफाइड गैस-काइनेटिक वेव-पार्टिकल (UGKWP) विधि, इस अंतर को पाटने का प्रयास करती हैं, जिसमें समाधान को तरंग (wave) और कण (particle) घटकों में विभाजित किया जाता है। हालाँकि, पारंपरिक UGKWP वैश्विक समय-चरण विभाजन (global time-step splitting) पर निर्भर करता है जहाँ कण अंश को वैश्विक समय चरण $\Delta t$ पर $e^{-\Delta t/\tau}$ टकराव प्रायिकता द्वारा निर्धारित किया जाता है। गैर-समान मेष (non-uniform meshes) पर, एक एकल अत्यधिक परिष्कृत सेल पूरे डोमेन में कण अंश को कृत्रिम रूप से पक्षपाती बना सकता है, जिससे विभिन्न स्थानीय विश्रांति पैमानों (relaxation scales) वाले क्षेत्रों में अक्षमता या अशुद्धता उत्पन्न हो सकती है।

2. कार्यप्रणाली: तरंग-कण अपघटन (WPD)

लेखक एक नया तरंग-कण अपघटन (WPD) प्रस्तावित करते हैं जिसे सीधे नियंत्रित आंशिक अवकल समीकरण (PDE) के स्तर पर तैयार किया गया है, न कि एक संख्यात्मक अपडेट रणनीति के रूप में। मुख्य नवाचार एक स्थानीय काइनेटिक होराइजन (local kinetic horizon), $T$, की परिभाषा है, जो काइनेटिक इतिहास को दो अलग-अलग घटकों में विभाजित करता है।

2.1 सैद्धांतिक सूत्रीकरण

एक काइनेटिक विश्रांति समीकरण (जैसे BGK या शाखोव मॉडल) के अभिलक्षण समाकल समाधान (characteristic integral solution) से शुरू करते हुए, वितरण फलन $f$ को इस प्रकार विभाजित किया जाता है:
$$f(x, \xi, t) = W(x, \xi, t) + P(x, \xi, t)$$
जहाँ:

• तरंग घटक ($W$): स्थानीय काइनेटिक होराइजन $T$ पर विश्लेषणात्मक रूप से संचित विश्राप्ति इतिहास का प्रतिनिधित्व करता है। इसे समय $0$से$T$ तक विश्रांति कर्नेल (relaxation kernel) के समाकल द्वारा परिभाषित किया जाता है।
• कण घटक ($P$): क्षितिज (horizon) के परे जीवित रहने वाले टकराव रहित मुक्त-परिवहन स्मृति (collisionless free-transport memory) का प्रतिनिधित्व करता है। इसे $e^{-\Phi(T)}f(x-\xi T, \xi, t-T)$ द्वारा परिभाषित किया गया है, जहाँ $\Phi(T)$ संचयी टकराव आवृत्ति है।

यह अपघटन एक एकीकृत समीकरण प्रणाली प्रदान करता है:

1. तरंग समीकरण (Wave Equation): एक परिवहन समीकरण जिसमें क्षितिज सीमा पर विनिमय को दर्शाने वाला एक स्रोत पद (source term) होता है।
2. कण समीकरण (Particle Equation): एक परिवहन समीकरण जो क्षितिज सीमा स्रोत द्वारा संचालित होता है।
3. कुल संरक्षण नियम (Total Conservation Law): एक स्रोत-मुक्त स्थूल संरक्षण नियम जहाँ कुल फ्लक्स, तरंग फ्लक्स और कण फ्लक्स का योग है।

2.2 हाइड्रोडायनामिक संरचना

वेव ऑपरेटर एक चैपन-एनस्कोग (Chapman–Enskog) विस्तार को स्वीकार करता है। परिणामी वेव फ्लक्स में होराइजन-भारित गुणांक ($A_0$ और $A_1$) होते हैं जो यूलर और नेवियर-स्टोक्स फ्लक्स को भारित करते हैं।

• निरंतरता सीमा में (स्थानीय होराइजन $T \gg$ विश्रांति समय $\tau$), गुणांक इकाई के करीब पहुँच जाते हैं, और वेव फ्लक्स पूर्ण नेवियर-स्टोक्स फ्लक्स को पुनः प्राप्त करता है, जबकि कण घटक लुप्त हो जाता है।
• विरल सीमा में ($T \ll \tau$), वेव फ्लक्स शून्य हो जाता है, और परिवहन पूरी तरह से कण घटक द्वारा ले जाया जाता है।

2.3 संख्यात्मक एल्गोरिदम

इस विधि को एक संरक्षित मैक्रो-माइक्रो फाइनाइट-वॉल्यूम स्कीम के रूप में लागू किया गया है:

• कुल अपडेट: संरक्षित चरों को वेव फ्लक्स और कण फ्लक्स के योग का उपयोग करके अपडेट किया जाता है।
• वेव फ्लक्स: नेवियर-स्टोक्स स्तर पर एक गैस-काइनेटिक स्कीम (GKS) का उपयोग करके, होराइजन-भारित गुणांकों को शामिल करते हुए, मूल्यांकन किया जाता है।
• कण फ्लक्स: या तो एक नियतात्मक डिस्क्रीट-ऑर्डिनेट ($S_N$) विधि या एक मोंटे कार्लो (MC) प्रतिनिधित्व का उपयोग करके मूल्यांकन किया जाता है।
• युग्मन (Coupling): कण घटक को समाकल समाधान के टकराव रहित कारक के माध्यम से गतिशील रूप से उत्पन्न किया जाता है। स्थानीय होराइजन $T$ को सेल-स्केल विकास समय ($T = \text{CFL}_l \Delta t_{i,loc}$) द्वारा परिभाषित किया जाता है, जो यह सुनिश्चित करता है कि विभाजन वैश्विक समय चरण के बजाय मेष-स्केल विश्रांति प्रक्रिया के लिए स्थानीय है।

3. प्रमुख योगदान

1. PDE-स्तर का अपघटन: पिछले तरीकों के विपरीत, जो वैश्विक समय-चरण प्रायिकताओं के आधार पर समाधान को विभाजित करते हैं, यह कार्य एक स्थानीय काइनेटिक होराइजन का उपयोग करके PDE स्तर पर वेव-पार्टिकल स्प्लिट को परिभाषित करता है। यह सुनिश्चित करता है कि विभाजन विकास काल और विश्रांति समय ($T/\tau$) के स्थानीय अनुपात द्वारा नियंत्रित है।
2. भौतिक रूप से सुसंगत कण जनसंख्या: कण घटक को एक टकराव रहित कारक द्वारा उत्पन्न एक भिन्नात्मक काइनेटिक जनसंख्या के रूप में प्राप्त किया गया है, न कि एक हस्ताक्षरित गणितीय अवशेष (signed mathematical residual) के रूप में जिसे साम्यावस्था अवस्था को घटाकर प्राप्त किया जाता है।
3. परिस्थिति-अनुकूल स्केलिंग (Regime-Adaptive Scaling): यह विधि सक्रिय काइनेटिक डिग्री ऑफ फ्रीडम (DoF) को स्थानीय भौतिक व्यवस्था के साथ गतिशील रूप से स्केल करने के लिए एक कठोर तंत्र प्रदान करती है। निरंतरता क्षेत्रों में, कण जनसंख्या तेजी से कम हो जाती है; विरल क्षेत्रों में, यह प्रभावी हो जाती है।
4. एकीकृत ढांचा: यह फॉर्मूलेशन एक एकल प्रणाली प्रदान करता है जो पूरे क्नुडसन स्पेक्ट्रम (Knudsen spectrum) में मान्य है, जिसमें एक स्रोत-मुक्त कुल संरक्षण नियम, एक वेव समीकरण और एक कण समीकरण शामिल है।

4. संख्यात्मक परिणाम

लेखक एक, दो और तीन-आयामी संख्यात्मक बेंचमार्क के माध्यम से ढांचे को मान्य करते हैं:

• 1D बेंचमार्क:

  • स्मूथ सटीकता: विरल, संक्रमण और निकट-निरंतरता व्यवस्थाओं में सुचारू आवधिक प्रवाह (smooth periodic flows) में द्वितीय-क्रम स्थानिक सटीकता प्रदर्शित की गई।
  • सोड शॉक ट्यूब (Sod Shock Tube): विभिन्न क्नुडसन संख्याओं में घनत्व, वेग और तापमान के लिए UGKS संदर्भ समाधानों के साथ सहमति दिखाई। इसमें कोई कृत्रिम संख्यात्मक प्रसार (numerical diffusion) नहीं आया।
  • कुएट प्रवाह (Couette Flow): दीवार-बद्ध प्रवाहों में $Kn=1$से$10^{-4}$ तक वेग स्लिप और घनत्व विविधताओं को सफलतापूर्वक कैप्चर किया।
  • नॉर्मल शॉक संरचना: मजबूत शॉक ($M=3, 8, 10$) के लिए स्थिर गैर-साम्यावस्था परतों और उच्च-क्रम मोमेंट्स को सटीक रूप से हल किया।

• 2D बेंचमार्क:

  • लिड-ड्रिवन कैविटी (Lid-Driven Cavity): विरल ($Kn=1, 0.075$) और निरंतर ($Re=1000, 3200$) दोनों व्यवस्थाओं में प्रदर्शन को मान्य किया। विधि ने मशीन परिशुद्धता तक द्रव्यमान संरक्षण बनाए रखा और निरंतरता सीमा में मानक श्यान भंवर संरचनाओं (viscous vortex structures) को पुनः प्राप्त किया।
  • हाइपरसोनिक सिलेंडर फ्लो: $M_\infty=8$ पर $Kn=1, 10^{-2}, 10^{-3}$ के लिए परीक्षण किया गया। नियतात्मक ($S_N$) और स्टोकेस्टिक (MC) दोनों कार्यान्वयन संदर्भ समाधानों के साथ सुसंगत रहे। परिणामों ने निरंतरता सीमा के करीब पहुँचने पर काइनेटिक डिग्री ऑफ फ्रीडम के सक्रिय न्यूनीकरण को रेखांकित किया।

• 3D बेंचमार्क:

  • X38 कॉन्फ़िगरेशन: $M_\infty=5$ और $AoA=15^\circ$ पर जटिल 3D ज्यामिति के लिए $Kn=10^{-1}, 10^{-2}, 10^{-4}$ पर लागू किया गया। विधि ने जटिल ज्यामिति को सफलतापूर्वक संभाला, जो वेव-प्रधान और कण-प्रधान क्षेत्रों के बीच सुचारू संक्रमण दिखाया।

• दक्षता: सिलेंडर फ्लो परीक्षणों में, वैश्विक-विभाजन विधियों (UGKWP) की तुलना में WPD-MC (स्थानीय) विधि ने प्रभावी काइनेटिक डिग्री ऑफ फ्रीडम और वॉल टाइम में महत्वपूर्ण कमी दिखाई, विशेष रूप से $Kn=10^{-3}$ पर जहाँ प्रभावी DoF को ~88 के कारक से कम कर दिया गया था।

5. महत्व और दावे

लेखक दावा करते हैं कि प्रस्तावित वेव-पार्टिकल अपघटन मल्टीस्केल प्रवाह सिमुलेशन के लिए एक मजबूत और परिस्थिति-अनुकूल ढांचा प्रदान करता है। इसका प्राथमिक महत्व निम्नलिखित है:

• वैश्विक समय चरणों से अलगाव: काइनेटिक होराइजन के माध्यम से स्थानीय रूप से स्प्लिट को परिभाषित करके, यह विधि गैर-समान मेष पर वैश्विक समय-चरण बाधाओं के कारण होने वाले कृत्रिम पक्षपात से बचती है।
• एसिम्प्टोटिक प्रिजर्वेशन (Asymptotic Preservation): औपचारिक विश्लेषण और संख्यात्मक परीक्षण पुष्टि करते हैं कि यह योजना निरंतरता क्षेत्र में सही नेवियर-स्टोक्स सीमा को पुनः प्राप्त करती है और सूक्ष्म पैमानों को हल किए बिना विरल क्षेत्रों में काइनेटिक निरंतरता बनाए रखती है।
• भौतिक सुसंगतता: कण घटक को एक भौतिक रूप से सुसंगत भिन्नात्मक जनसंख्या के रूप में व्याख्या किया गया है, जो हस्ताक्षरित अवशेषों की तुलना में गैर-साम्यावस्था परिवहन का अधिक स्वाभाविक प्रतिनिधित्व प्रदान करता है।

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