Systematic comparison of VMEC and HINT equilibrium calculations for finite-beta LHD plasmas

तकनीकी सारांश: परिमित-$\beta$ LHD प्लाज्मा के लिए VMEC और HINT साम्यावस्था (Equilibrium) गणनाओं का व्यवस्थित तुलनात्मक अध्ययन

समस्या विवरण
स्टेलरैटर भौतिकी में, तीन-आयामी मैग्नेटोहाइड्रोडायनामिक (MHD) साम्यावस्था की गणना अक्सर VMEC सॉल्वर का उपयोग करके की जाती है। VMEC रूपांतरणात्मक सिद्धांत (variational principle) पर निर्भर करता है और कड़ाई से पूर्णतः संलग्नी फ्लक्स सतहों (perfectly nested flux surfaces) के अस्तित्व की धारणा रखता है। हालांकि, परिमित-$\beta$ प्लाज्मा में, बढ़ता हुआ प्लाज्मा दबाव चुंबकीय क्षेत्र की टोपोलॉजी को बदल सकता है, जिससे चुंबकीय द्वीप (magnetic islands) और स्टोकेस्टिक (stochastic) क्षेत्र उत्पन्न होते हैं। ये घटनाएँ, विशेष रूप से शून्य शुद्ध-धारा वाले उपकरणों में जहाँ वे पिरश-श्ल्यूटर (Pfirsch-Schlüter - PS) धाराओं द्वारा संचालित होती हैं, संलग्नी फ्लक्स सतह की धारणा की वैधता को बाधित कर सकती हैं। जब चुंबकीय द्वीप आपस में मिलते हैं, तो अंतिम बंद फ्लक्स सतह (LCFS) टूट जाती है, जिससे घिरा हुआ प्लाज्मा आयतन कम हो जाता है—यह एक ऐसी सीमा है जो शाफ्रानोव शिफ्ट (Shafranov shift) की तुलना में अधिक गंभीर देखी गई है। जबकि HINT कोड चुंबकीय क्षेत्रों और दबाव के गतिशील समीकरणों को हल करता है और इसमें संलग्नी सतहों की धारणा नहीं होती है, LHD में विभिन्न विन्यासों और बीटा मानों में VMEC और HINT के बीच एक व्यवस्थित तुलना करना इन धारणाओं के प्रभाव को स्पष्ट करने के लिए आवश्यक था।

कार्यप्रणाली
इस अध्ययन ने LHD प्लाज्मा के लिए VMEC और HINT के बीच साम्यावस्था गणनाओं का एक व्यवस्थित तुलनात्मक अध्ययन किया। जांच में तीन विशिष्ट निर्वात चुंबकीय-अक्ष विन्यासों पर ध्यान केंद्रित किया गया: $R_{ax}^V = 3.53$ मी (अंदर की ओर स्थानांतरित), $3.60$मी (मानक), और$3.85$मी (बाहर की ओर स्थानांतरित)। विश्लेषण में ऑन-एक्सिस बीटा मान ($\beta_0$) का दायरा$0.0%$से$5.0%$तक रखा गया, जिसे$0.5%$ के अंतराल में मापा गया।

दोनों कोडों ने एक प्रारंभिक दबाव प्रोफाइल का उपयोग किया जिसे $p = p_0(1-s)$ के रूप में परिभाषित किया गया था, जहाँ $s$ सामान्यीकृत टोरोइडल फ्लक्स है। अध्ययन ने तीन प्राथमिक मेट्रिक्स की तुलना की:

1. चुंबकीय अक्ष की स्थिति ($\langle R_{ax} \rangle$)।
2. अक्ष पर रोटेशनल ट्रांसफॉर्म ($\iota/2\pi$)।
3. स्पष्ट LCFS द्वारा घिरा आयतन।

HINT कोड, जो चुंबकीय द्वीपों और स्टोकेस्टिक क्षेत्रों को हल करने के लिए रिलैक्सेशन विधि का उपयोग करता है, की तुलना VMEC से की गई, जो संलग्नी फ्लक्स सतह की बाधा को लागू करता है।

मुख्य परिणाम

• निम्न-$\beta$ शासन (Low-$\beta$ Regime): निम्न $\beta_0$ मानों के लिए, VMEC और HINT ने सुसंगत साम्यावस्था समाधान प्रस्तुत किए। यह इंगित करता है कि इस शासन में संलग्नी फ्लक्स सतह की धारणा काफी हद तक वैध बनी रहती है, और दोनों कोडों के बीच 3D MHD प्रतिक्रियाएं तुलनीय हैं।
• क्रिटिकल $\beta_0$ और विचलन: एक विन्यास-निर्भर क्रिटिकल $\beta_0$ की पहचान की गई। इस थ्रेसहोल्ड के आगे, दोनों कोडों के समाधान अलग होने लगे। HINT गणनाओं में, चुंबकीय अक्ष ने VMEC की तुलना में बड़ा शाफ्रानोव शिफ्ट और उच्च रोटेशनल ट्रांसफॉर्म प्रदर्शित किया।
• फ्लक्स सतह का टूटना और आयतन: HINT में, उच्च बीटा मानों पर स्पष्ट LCFS द्वारा घिरा आयतन घट गया। इस कमी का कारण प्लाज्मा किनारे के पास स्टोकेस्टिक चुंबकीय क्षेत्रों का विकास है, जो फ्लक्स सतहों को तोड़ देता है। इसके विपरीत, VMEC, जो अपनी संलग्नी सतह की धारणा द्वारा बाधित है, $\beta_0$ बढ़ने के साथ बढ़े हुए आयतन में निरंतर वृद्धि दिखाता है, जिससे वह फ्लक्स सतह के टूटने को दर्शाने में विफल रहता है।
• विन्यास निर्भरता: विन्यास के अंदर की ओर से बाहर की ओर स्थानांतरित होने के साथ क्रिटिकल $\beta_0$ का मान कम होता गया। $R_{ax}^V = 3.85$ मी (बाहर की ओर स्थानांतरित) वाले विन्यास को 3D साम्यावस्था प्रतिक्रियाओं के प्रति सबसे अधिक संवेदनशील पाया गया। यह संवेदनशीलता इसलिए उत्पन्न होती है क्योंकि $R_{ax}^V$ के साथ हेलिकल रिपल (helical ripple) बढ़ता है, जो मजबूत PS धाराओं को प्रेरित करता है, जबकि चुंबकीय शियर (magnetic shear) कम हो जाता है। मजबूत PS धाराओं और कम शियर का संयोजन बड़े चुंबकीय द्वीपों और कम बीटा मानों पर शीघ्र स्टोकेस्टिसिटी की ओर ले जाता है।

महत्व और दावे
यह शोध स्थापित करता है कि LHD प्लाज्मा में 3D साम्यावस्था प्रतिक्रियाएं तब अधिक महत्वपूर्ण हो जाती हैं जब विन्यास 'इनवर्ड' से 'आउटवर्ड' की ओर स्थानांतरित होता है। अध्ययन यह प्रदर्शित करता है कि परिमित बीटा पर, PS करंट-संचालित विक्षोभों और परिणामी किनारे की स्टोकेस्टिसिटी के कारण VMEC में केंद्रीय संलग्नी फ्लक्स सतह की धारणा समझौतापूर्ण हो जाती है।

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