One dimensional high-order moment models with realistic collisions for nonequilibrium ion transport in weakly ionized plasmas

कल्पना कीजिए एक ऐसी दुनिया की जहाँ अदृश्य कण एक चमकते कोहरे में नृत्य कर रहे हैं, जो आपके शहर के नियॉन साइन से लेकर भविष्य के अंतरिक्ष यानों के इंजनों तक, सब कुछ संचालित कर रहे हैं। यह प्लाज्मा का क्षेत्र है, पदार्थ की चौथी अवस्था, जहाँ परमाणुओं से उनके इलेक्ट्रॉन छिन लिए जाते हैं और वे आवेशित आयनों और मुक्त इलेक्ट्रॉनों का एक अराजक सूप बन जाते हैं। इनमें से कई चमकते बादलों में, आयन (भारी, धनावेशित हिस्से) केवल आलस में नहीं तैरते; उन्हें विद्युत क्षेत्रों द्वारा धकेला जाता है, जिससे वे भीड़ भरे आर्केड में बम्पर कारों की तरह तटस्थ गैस परमाणुओं से टकराते हैं। आमतौर पर, वैज्ञानिक इन कणों की गति की भविष्यवाणी करने के लिए सरल नियमों का उपयोग करते हैं, यह मानते हुए कि वे सभी लगभग एक ही गति और तापमान पर चल रहे हैं, जैसे कि एक शांत भीड़ सड़क पर चल रही हो। लेकिन कभी-कभी, विद्युत धक्का इतना शक्तिशाली होता है, या हवा इतनी पतली होती है, कि आयन एक उन्माद में आ जाते हैं। वे जंगली अलग-अलग गति से चलने लगते हैं, जिनमें से कुछ आगे निकल जाते हैं जबकि अन्य पीछे रह जाते हैं, जिससे एक अराजक, एकतरफा भीड़ बन जाती है जो सरल नियमों को तोड़ देती है। इस "गैर-संतुलन" (non-equilibrium) अराजकता को समझना महत्वपूर्ण है क्योंकि यदि हम गणित में गलती करते हैं, तो हम उपग्रहों के लिए बेहतर इलेक्ट्रिक थ्रस्टर्स या अधिक कुशल औद्योगिक प्लाज्मा उपकरण डिजाइन नहीं कर सकते।

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

उन्मादी आयनों की कहानी

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

लेखकों ने इन आयनों का वर्णन करने के लिए नए नियमों, या "मोमेंट मॉडल्स" (moment models) को बनाकर इसे ठीक करने का निर्णय लिया। एक "मोमेंट" को भीड़ के व्यवहार के सारांश के रूप में सोचें। सबसे सरल मॉडल (3M) बस यह पूछता है: "वहाँ कितने बच्चे हैं, उनकी औसत गति क्या है, और उनका औसत तापमान क्या है?" यह एक भीड़ का वर्णन करने जैसा है कि, "वे 5 मील प्रति घंटे की गति से चल रहे हैं।" लेकिन एक अराजक प्लाज्मा में, यह पर्याप्त नहीं है। आयन औसतन 5 मील प्रति घंटे की गति से चल सकते हैं, लेकिन आधा हिस्सा स्थिर हो सकता है जबकि दूसरा आधा हिस्सा 10 मील प्रति घंटे की गति से दौड़ रहा हो। इसे पकड़ने के लिए, लेखकों ने अपने मॉडलों में अधिक "मोमेंट्स" जोड़े। उन्होंने एक 4M मॉडल जोड़ा जो यह जाँचता है कि क्या भीड़ एक दिशा में दूसरी दिशा की तुलना में अधिक गर्म है (anisotropy), एक 5M मॉडल जो देखता है कि गति वितरण कितना "तिरछा" (skewed) है (heat flux), और अंत में, एक शानदार 6M मॉडल जो सब कुछ जाँचता है, जिसमें यह भी शामिल है कि वितरण कितना "नुकीला" (kurtosis) है।

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

अपने नए "स्मार्ट लेंसों" का परीक्षण करने के लिए, लेखकों ने दो बहुत अलग प्लाज्मा परिदृश्यों के सिमुलेशन चलाए। पहला एक चमकते बॉक्स की तरह दो तैरती हुई दीवारों के बीच फंसा हुआ प्लाज्मा था। उन्होंने इसका परीक्षण 0.05 mTorr (लगभग निर्वात) से 50 mTorr (थोड़ा घना) के दबाव सीमा में किया। दूसरा परिदृश्य एक डायरेक्ट-करंट (DC) डिस्चार्ज था, जो एक नियॉन ट्यूब की तरह है जिसमें एक मजबूत विद्युत क्षेत्र आयनों को एक दीवार की ओर धकेलता है। उन्होंने इसे 200 से 500 mTorr की विस्तृत दबाव सीमा में चलाया।

परिणाम चौंकाने वाले थे। जब उन्होंने अपने नए मॉडलों की तुलना "गोल्ड स्टैंडर्ड" सिमुलेशन के विरुद्ध की जो प्रत्येक कण को ट्रैक करते हैं (जिसे PIC-MCC कहा जाता है), तो पुराने 3M और 4M मॉडल संघर्ष करते दिखे, विशेष रूप से दीवारों के पास शीथ (sheath) क्षेत्रों में जहाँ विद्युत क्षेत्र सबसे मजबूत होता है। वे अक्सर तापमान और गति को गलत बताते थे क्योंकि वे आयनों के एकतरफा आंदोलन को नहीं समझ पाते थे। BGK मॉडल, जो उस "तत्काल शांत होने" के अनुमान पर निर्भर करते हैं, और भी खराब थे, जो अक्सर भौतिकी को पूरी तरह से पकड़ने में विफल रहे।

हालाँकि, नया 6M मॉडल, जिसमें लंबवत तापमान (perpendicular temperature) और हीट फ्लक्स शामिल है, एक विजेता साबित हुआ। इसने कण-ट्रैकिंग सिमुलेशन के साथ लगभग पूरी तरह से मेल खाया। इसने सही ढंग से भविष्यवाणी की कि कम-दबाव वाले क्षेत्रों में, आयन अपनी गति वितरण में लंबी, भारी पूंछ (tails) विकसित करते हैं (जिसका अर्थ है कि कुछ आयन बहुत तेज़ गति से चल रहे हैं), और उच्च-दबाव वाले क्षेत्रों में, वे जिस दिशा में चल रहे हैं और जिस दिशा में घूम रहे हैं, उसके बीच तापमान के अंतर को विकसित करते हैं। सबसे महत्वपूर्ण बात यह है कि 6M मॉडल बिना किसी "शोर" या सांख्यिकीय धुंधलेपन के, आयन गति वितरण के वास्तविक आकार को पुनर्गठित कर सका, और यह साधारण फ्लूइड मॉडलों की गति के बराबर था।

लेखकों ने पाया कि 5M मॉडल भी बहुत अच्छा था, लेकिन 6M मॉडल का थोड़ा बढ़त थी, विशेष रूप से उच्च-दबाव वाले शीथ में जहाँ लंबवत तापमान (आयनों की "घूमने वाली" ऊर्जा) महत्वपूर्ण होने लगती है। उन्होंने दिखाया कि इन उच्च-क्रम के विवरणों को अनदेखा करने से यह बताने में त्रुटि होती है कि आयन दीवारों से कितनी ऊर्जा के साथ टकराते हैं, जो कंप्यूटर चिप्स को उकेरने (etching) या अंतरिक्ष यान के थ्रस्टर्स को डिजाइन करने जैसी चीजों के लिए महत्वपूर्ण है।

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