Geometric numerical discretization of electromagnetic quasineutral models

समस्या विवरण
चुंबकित प्लाज्मा का संख्यात्मक सिमुलेशन संलयन ऊर्जा (fusion energy), अंतरिक्ष मौसम और खगोल भौतिकी के लिए आवश्यक है, फिर भी यह लंबाई और समय के विशाल पृथक्करण (separation of scales) के कारण चुनौतियों का सामना करता है। जबकि गैर-मैक्सवेलियन वितरण (non-Maxwellian distributions) और तरंग-कण अंतःक्रियाओं (wave-particle interactions) को पकड़ने के लिए विलोस-मैक्सवेल (Vlasov-Maxwell) प्रणाली पर आधारित गतिज मॉडल (kinetic models) आवश्यक हैं, डेबाई लंबाई (Debye length) और उच्च-आवृत्ति वाली विद्युत चुम्बकीय तरंगों (जैसे लैंगमुइर और प्रकाश तरंगें) को हल करना कई अनुप्रयोगों के लिए अत्यधिक गणनात्मक रूप से महंगा है। अर्ध-तटस्थ (quasineutral) मॉडल, जो आवेश तटस्थता ($\rho=0$) मानते हैं और इन छोटे पैमानों को फ़िल्टर कर देते हैं, एक अधिक कुशल विकल्प प्रदान करते हैं। हालाँकि, अर्ध-तटस्थ प्रणालियों के लिए मौजूदा संख्यात्मक दृष्टिकोण अक्सर उन 'एसिम्प्टोटिक-प्रिजर्विंग' (asymptotic-preserving) स्कीम्स पर निर्भर करते हैं जिन्हें डेबाई लंबाई के शून्य होने पर उच्च-आवृत्ति मोड को दबाने के लिए संख्यात्मक डैम्पिंग (numerical damping) की आवश्यकता होती है, या उनमें आधुनिक संरचना-संरक्षण (structure-preserving) कण-इन-सेल (PIC) विधियों में पाई जाने वाली ज्यामितीय संरचना (जैसे ऊर्जा संरक्षण और सटीक विचलन प्रतिबंध) का अभाव होता है।

कार्यप्रणाली
लेखक पूर्णतः गतिज, अर्ध-तटस्थ विलोस-मैक्सवेल समीकरणों को हल करने के लिए ज्यामितीय विद्युत चुम्बकीय PIC ढांचे, GEMPICX का विस्तार करते हैं। यह कार्यप्रणाली निम्नलिखित स्तंभों पर आधारित है:

• ज्यामितीय ढांचा (Geometric Framework): यह योजना द्वैत ग्रिड (प्राथमिक और द्वैत मेश) पर मिमेटिक फाइनाइट डिफरेंस (mimetic finite differences) का उपयोग करती है। क्षेत्र चर (field variables) को असतत स्थानों (नोड्स, किनारों, फलकों, आयतनों) पर परिभाषित किया गया है जो एक असतत डी रॅम कॉम्प्लेक्स (discrete de Rham complex) बनाते हैं, जो चुंबकीय क्षेत्र की विचलन-मुक्त (divergence-free) स्थिति जैसी टोपोलॉजिकल गुणों के सटीक संरक्षण को सुनिश्चित करता है।
• वैरियेशनल सिद्धांत (Variational Principle): इस विविक्तीकरण (discretization) को ट्रोंसी और कैंपोरियल (2015) द्वारा दिए गए लैग्रेंजियन पर आधारित एक असतत क्रिया सिद्धांत (discrete action principle) से प्राप्त किया गया है। यह सुनिश्चित करता है कि कणों और क्षेत्रों के लिए अर्ध-असतत गति के समीकरण (semi-discrete equations of motion) अंतर्निहित हैमिल्टनियन संरचना के अनुरूप हैं।
• अर्ध-तटस्थ सूत्रीकरण (Quasineutral Formulation): अन्य एसिम्प्टोटिक-प्रिजर्विंग दृष्टिकोणों के विपरीत, जो पूर्ण विलोस-मैक्सवेल प्रणाली को हल करते हैं और $\epsilon_0 \to 0$ की सीमा लेते हैं, यह विधि सीधे अर्ध-तटस्थ सीमा में व्युत्पन्न की गई है। फलस्वरूप, विस्थापन धारा (displacement current) पद अनुपस्थित है, और विद्युत क्षेत्र को प्रत्येक समय चरण में समय के व्युत्पन्न (time derivative) के माध्यम से विकसित नहीं किया जा सकता है। इसके बजाय, विद्युत क्षेत्र को प्रत्येक समय चरण में एक विविक्त कर्ल-कर्ल समीकरण (Equation 2.25) को हल करके अप्रत्यक्ष रूप से (implicitly) प्राप्त किया जाता है, जो गैर-विलक्षण (non-singular) और धनात्मक निश्चित (positive definite) है।
• धारा विचलन सुधार (Current Divergence Correction): मशीन परिशुद्धता (machine precision) पर $\nabla \cdot \mathbf{J} = 0$ की अर्ध-तटस्थ स्थिति को संतुष्ट करने के लिए (जो सांख्यिकीय शोर के कारण कण निक्षेपण (particle deposition) द्वारा स्वतः संतुष्ट नहीं होती है), एक लैग्रेंज मल्टीप्लायर पेश किया गया है। यह एक "विचलन-सुधारक" (divergence-corrector) विद्युत क्षेत्र ($\mathbf{E}^\star = -\nabla \psi^\star$) के रूप में प्रकट होता है, जिसे मुख्य समय चरण से पहले कण वेगों को समायोजित करने के लिए एक पॉइसन-समान समीकरण (Equation 6.8) के माध्यम से हल किया जाता है।
• टाइम-स्टेपिंग (Time-Stepping): एल्गोरिदम एक द्वितीय-क्रम सटीक स्प्लिट स्कीम का उपयोग करता है। इसमें धारा विचलन को ठीक करना, वेक्टर क्षमता $\mathbf{A}$ (कूलम्ब गेज $\nabla \cdot \mathbf{A} = 0$ का उपयोग करके) के लिए हल करना, विद्युत क्षेत्र $\mathbf{E}$ के लिए हल करना, और अंत में बोरिस पुशर (Boris pusher) का उपयोग करके कणों की स्थिति और वेग को अपडेट करना शामिल है।

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

1. GEMPICX का विस्तार: यह शोध पत्र पूर्णतः गतिज इलेक्ट्रॉनों और आयनों के लिए अर्ध-तटस्थ सीमा हेतु मिमेटिक फाइनाइट डिफरेंस GEMPIC ढांचे के पहले विस्तार को प्रस्तुत करता है।
2. प्रत्यक्ष अर्ध-तटस्थ व्युत्पत्ति: यह योजना किसी पूर्ण सॉल्वर की एसिम्प्टोटिक निरंतरता पर निर्भर रहने के बजाय सीधे अर्ध-तटस्थ सीमा में व्युत्पन्न की गई है। यह सुनिश्चित करता है कि संख्यात्मक समाधान में स्पूरियस (spurious) विद्युत चुम्बकीय और लैंगमुइर तरंगें बिना किसी संख्यात्मक डैम्पिंग के अनुपस्थित रहें।
3. संरचना संरक्षण: अर्ध-असतत सूत्रीकरण गणितीय रूप से चुंबकीय और गतिज ऊर्जाओं के योग के संरक्षण को सिद्ध करता है। यह मशीन परिशुद्धता के माध्यम से धारा घनत्व के शून्य-विचलन की आवश्यकता को भी बनाए रखता है।
4. अप्रत्यक्ष विद्युत क्षेत्र सॉल्वर: लेखक विद्युत क्षेत्र प्राप्त करने के लिए एक सुदृढ़ रैखिक प्रणाली (Equation 5.12) व्युत्पन्न करते हैं, जो क्षेत्र अपडेट में विद्युत विभव (electric potentials) की आवश्यकता को समाप्त करती है और विस्थापन धारा की अनुपस्थिति में मानक कर्ल-कर्ल ऑपरेटर से जुड़ी विलक्षणता (singularity) की समस्याओं से बचती है।
5. स्केलेबिलिटी (Scalability): कार्यान्वयन को MPI समानांतरकरण के साथ AMReX फ्रेमवर्क में एकीकृत किया गया है, जो HPC संसाधनों पर मजबूत और कमजोर स्केलिंग क्षमताओं का प्रदर्शन करता है।

परिणाम
ढांचे को आवधिक सीमा स्थितियों (periodic boundary conditions) के साथ श्रृंखला परीक्षणों के माध्यम से मान्य किया गया:

• विक्षेपण संबंध (Dispersion Relations): एकल-प्रजाति (केवल इलेक्ट्रॉन) और दोहरी-प्रजाति (इलेक्ट्रॉन-आयन) प्लाज्मा के सिमुलेशन ने पृष्ठभूमि चुंबकीय क्षेत्र के समानांतर और लंबवत तरंगों के सैद्धांतिक विक्षेपण संबंधों को सफलतापूर्वक पुनरुत्पादित किया।
  • समानांतर प्रसार (Parallel Propagation): योजना ने इलेक्ट्रॉन साइक्लोट्रॉन वेव (ECW) और आयन साइक्लोट्रॉन वेव (ICW), साथ ही डैम्प्ड उच्च-क्रम अल्फवेन-साइक्लोट्रॉन मोड की शंक्वाकार संरचनाओं को सही ढंग से पकड़ा। जैसा कि अपेक्षित था, लैंगमुइर तरंगें और विद्युत चुम्बकीय L/R मोड अनुपस्थित थे।
  • लंबवत प्रसार (Perpendicular Propagation): सिमुलेशन ने बर्नस्टीन तरंगों (Bernstein waves) और साइक्लोट्रॉन हार्मोनिक अनुनाद को पुनरुत्पादित किया। दोहरी-प्रजाति मामलों में, आयन बर्नस्टीन तरंगों और हाइब्रिड वेव को देखा गया। अर्ध-तटस्थ धारणा के कारण O-मोड और X-मोड सही ढंग से अनुपस्थित थे।
• तरंग डैम्पिंग (Wave Damping): एक एकल साइनुसोइडल विक्षोभ (sinusoidal perturbation) वाले परीक्षण ने साइक्लोट्रॉन वेव की सैद्धांतिक रूप से अनुमानित डैम्पिंग दर को पकड़ने की कोड की क्षमता को प्रदर्शित किया, जो उच्च-क्रम गतिज प्रक्रियाओं को संभालने में योजना की सटीकता को प्रमाणित करता है।
• कंप्यूटेशनल दक्षता: विद्युत क्षेत्र के लिए पूर्ण कर्ल-कर्ल सिस्टम (Equation 5.12) को हल करने और एक सस्ते विकर्ण सन्निकटन (diagonal approximation) (Equation 5.13) के बीच तुलना ने दिखाया कि सन्निकटन ने परीक्षण किए गए मामलों के लिए नगण्य सटीकता हानि के साथ कंप्यूटेशनल लागत को लगभग 25% कम कर दिया।
• स्केलेबिलिटी: वितरित-मेमोरी आर्किटेक्चर पर स्ट्रॉन्ग और वीक स्केलिंग अध्ययनों ने 32 MPI कार्यों तक लगभग आदर्श स्पीडअप दिखाया और 128 कार्यों तक उच्च दक्षता बनाए रखी, जो दर्शाता है कि कोड बड़े पैमाने के सिमुलेशन के लिए उपयुक्त है।

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

लेखक सीमाओं और भविष्य के कार्य के क्षेत्रों को विनम्रतापूर्वक नोट करते हैं:

• सीमा स्थितियाँ (Boundary Conditions): वर्तमान कार्यान्वयन कूलम्ब गेज को आसानी से लागू करने और लैपलेस ऑपरेटर को उलटने (invert) के लिए आवधिक सीमा स्थितियों पर निर्भर करता है। गैर-आवधिक सीमाओं (जैसे फ्यूजन उपकरणों के लिए) तक इसका विस्तार करने के लिए अतिरिक्त बाधाओं और अनुकूलता विश्लेषण की आवश्यकता है।
• सटीकता का क्रम (Order of Accuracy): यह कार्य द्वितीय-क्रम Hodge ऑपरेटरों का उपयोग करता है। हालांकि मिमेटिक ढांचे के भीतर उच्च-क्रम ऑपरेटर संभव हैं, लेकिन अर्ध-तटस्थ सिमुलेशन (जिसमें उच्च-आवृत्ति मोड की कमी होती है) के लिए उनके विशिष्ट लाभ अभी भी एक खुला प्रश्न है।
• शॉक हैंडलिंग (Shock Handling): यह विधि सुचारू क्षेत्रों (smooth fields) को मानती है। हालांकि PIC विधियाँ आम तौर पर कण वितरण के माध्यम से शॉक को संभालती हैं, इस विशिष्ट ज्यामितीय विविक्तिकरण का मजबूत ग्रेडिएंट या शॉक वाले शासन (regimes) में व्यवहार की आगे जांच की आवश्यकता है।
• टाइम डिस्क्रीटाइजेशन (Time Discretization): ऊर्जा संरक्षण केवल अर्ध-असतत स्तर पर सटीक है; समय विविक्तिकरण के बाद, संरक्षण केवल समय इंटीग्रेटर के क्रम तक सटीक है। सटीक असतत ऊर्जा संरक्षण प्राप्त करने के लिए एक पूर्ण वैरियेशनल टाइम डिस्क्रीटाइजेशन का सुझाव दिया गया है।