DOI: https://doi.org/10.1088/1741-4326/ae44ad
تاريخ النشر: 2026-02-11
المؤلف: Stuart Royce Sands Benjamin وآخرون
الموضوع الرئيسي: أبحاث الاندماج المغناطيسي المحتجز
نظرة عامة
في هذه الدراسة، تم تقييم مقياس الاستقرار لأنماط التمزق النيوكلاسيكية (NTMs)، وبشكل خاص الحد الأدنى لعرض الجزيرة المستقر هامشياً \( w^*_m \)، عبر مجموعة بيانات تضم 14,651 توازنات لمحطات توكاماك تجريبية عالية الحقل. أظهرت التحليلات أن التوكاماك الأكبر مع تقليل الاستطالة و/أو زيادة نصف القطر الصغير أظهر زيادة كبيرة في \( w^*_m \)، والتي تعزى أساساً إلى انخفاض في دفع البوتستراب. تم إنشاء التوازنات باستخدام طرق مونت كارلو، وتم تطبيعها لتحقيق أكثر من 400 ميغاوات من طاقة الاندماج، وخضعت لتحليل التمزق من خلال معادلة رذرفورد المعدلة لأنماط بوليدية وتوريدية رنانة مختلفة حتى \( n = 4 \).
أشارت النتائج إلى أن بعض تكوينات التصميم، وخاصة التوكاماك الكبيرة ومنخفضة الاستطالة، يمكن أن تحقق تحسينات كبيرة في \( w^*_m \)، وهو أمر حاسم لمنع نمو NTM. تم استخدام تقنيات التعلم الآلي لكشف الاعتماديات الأساسية في \( w^*_m \)، مما يبرز زيادة الاستقرار المرتبطة بزيادة الثلاثية بسبب انخفاض قيم \( \Delta’ \). تشير هذه النتائج إلى أنه قد يكون من الممكن تحقيق استقرار NTM السلبي في عمليات المحطة التجريبية من خلال تصميم سيناريوهات بعناية. ستسعى الأبحاث المستقبلية إلى إنشاء سير عمل للتحقق تجريبياً من \( w^*_m \) كمقياس موثوق لاستقرار NTM في تصميم التوكاماك.
مقدمة
تناقش المقدمة نمط التمزق النيوكلاسيكي (NTM)، وهو عدم استقرار بلازما مقاوم يشكل مخاطر كبيرة على تشغيل التوكاماك، خاصة في سياق محطات الاندماج التجريبية من الجيل التالي (FPPs) التي تهدف إلى سيناريوهات ذات عائد مرتفع. يتم بدء NTMs بواسطة تدرج ضغط البلازما وتتطلب “جزيرة” مغناطيسية “بذور” لنموها. طرق الاستقرار الحالية، مثل دفع تيار الإلكترون الدائري وحقن الشعاع المحايد خارج المحور، ليست ممكنة في التصاميم الأولية لمحطات الاندماج فائقة التوصيل عالية الحرارة بسبب القيود التكنولوجية والاقتصادية. يركز المؤلفون على تحدي استقرار NTMs في توكاماك H-mode عالية الحقل بدون هذه الطرق الشائعة للتحكم.
لمعالجة ذلك، تستخدم الدراسة معادلة رذرفورد المعدلة (MRE) لتحديد الظروف التي يبقى فيها حجم الجزر البذور صغيراً جداً لتحفيز نمو NTM. من خلال تحليل قاعدة بيانات توازن اصطناعية لمختلف تصميمات التوكاماك، يجد المؤلفون أن الحد الأدنى لعرض الجزيرة المستقر هامشياً، المشار إليه بـ $w^*_m$، يتأثر بشكل كبير بالمعلمات الهندسية مثل نصف القطر الرئيسي ($R_0$)، والاستطالة ($\kappa$)، ونصف القطر الصغير ($a$). تشير النتائج إلى أن تحسين هذه المعلمات يمكن أن يعزز $w^*_m$، مما يحسن استقرار NTM السلبي. يوضح البحث هيكله، مفصلاً الأقسام التالية التي ستتناول نموذج فيزياء التمزق، وطرق توليد قاعدة البيانات، ونتائج التحليل، ومناقشات حول العمل المستقبلي المتعلق بتوقع احتمالات بدء NTM.
النتائج
في هذا القسم، يتم تقديم نتائج تحليل التمزق الذي تم إجراؤه على قاعدة بيانات H-mode لمونت كارلو. شمل التحليل توليد 14,651 توازنات مستقرة مثالية-MHD، والتي تطابقت مع 7,418 هندسة توكاماك فريدة. بعد تطبيق معايير التبادل المقاوم، تم استبعاد 977 توازناً، مما ترك 13,719 سيناريو لمزيد من التحليل. من بين هذه، أظهرت 1,375 توازناً قيم إيجابية لمعامل الاستقرار $\Delta’$، ومع ذلك، وُجد أن 114 منها فقط غير مستقرة خطياً مع $\Delta’ > \Delta’_c$. بالإضافة إلى ذلك، تم اعتبار 844 توازناً غير مناسبة للتحليل بسبب كون عرض الجزيرة المستقر هامشياً أقل من عرض الطبقة الخطية ($w^*_m < \delta_r$). تشير النتائج إلى أن 2,136 توازناً، أو 16.9% من الحالات المتبقية، حققت استقراراً كاملاً، يتميز بشرط حيث $\frac{dw}{dt} < 0$ لجميع الأنماط داخل قاعدة البيانات. في هذه السيناريوهات المستقرة بالكامل، تم تعريف الحد الأدنى لعرض الجزيرة المستقر هامشياً بحجم الجهاز، مع $w^*_m = 1$. تسلط النتائج الضوء على الاتجاهات الكبيرة في استقرار التوازنات عبر مساحة الإدخال المحددة بواسطة معلمات مثل $R_0$، $\epsilon$، $\delta$، $\kappa$، والتيار المنظم $I/I^*$. ستتناول المناقشات الإضافية في الأقسام التالية الفيزياء الكامنة وراء هذه الاتجاهات واعتماديات $\Delta'$ ذات الهيكل الفردي عبر قاعدة البيانات.
مناقشة
في هذا القسم، يناقش المؤلفون الإطار لتحليل استقرار توازنات البلازما المحورية، مع التركيز على الديناميات الخطية وغير الخطية لأنماط التمزق. يعتمد التحليل على شرط التوازن \( J \times B = \nabla p \)، مع تمثيل المجال المغناطيسي من حيث التدفق المغناطيسي البوليدي \( \Psi \) وزاوية التوريد \( \phi \). يتم تقييم استقرار التمزق باستخدام كود DCON المقاوم، الذي يحسب الاضطرابات في البلازما ويحدد الظروف التي تصبح فيها هذه الاضطرابات غير مستقرة. يستنتج المؤلفون مجموعة من المعادلات التي تحكم ديناميات أنماط التمزق ذات الهيكل الفردي، مع التأكيد على أهمية الحد الأدنى لعرض الجزيرة المستقر هامشياً \( w^*_m \) كمقياس حاسم لتقييم الاستقرار ضد أنماط التمزق النيوكلاسيكية (NTMs).
تشير النتائج إلى أن الحد الأدنى لعرض الجزيرة المستقر هامشياً يتفاوت بشكل كبير عبر قاعدة بيانات H-mode الاصطناعية، مع قيم متوسطة، قصوى، ودنيا تبلغ \( 6.8 \times 10^{-3} \)، \( 5.8 \times 10^{-2} \)، و\( 9.4 \times 10^{-5} \) على التوالي. يكشف التحليل أن زيادة نصف القطر الرئيسي ونصف القطر الصغير أو تعديل الاستطالة يمكن أن يعزز استقرار NTM، مما يشير إلى أن التوكاماك الأكبر أو الأقل أداءً قد تكون أكثر ملاءمة لاستقرار العمليات. كما يشير المؤلفون إلى أن بعض التأثيرات الفيزيائية، مثل تيارات استقطاب الأيونات والدوران التفاضلي، مستبعدة من تحليلهم، مؤكدين أن استبعادها مبرر بالنظر إلى التركيز على تحديد أنظمة التصميم المواتية لمحطات الاندماج المستقبلية. بشكل عام، تسلط الدراسة الضوء على التفاعل المعقد بين التأثيرات المثبتة والمزعزعة المعتمدة على التوازن على NTMs، مما يوفر رؤى لتحسين تصميمات التوكاماك.
DOI: https://doi.org/10.1088/1741-4326/ae44ad
Publication Date: 2026-02-11
Author(s): Stuart Royce Sands Benjamin et al.
Primary Topic: Magnetic confinement fusion research
Overview
In this study, the stability metric for neoclassical tearing modes (NTMs), specifically the minimum marginally stable island width \( w^*_m \), was evaluated across a dataset of 14,651 inductive high-field tokamak pilot plant equilibria. The analysis revealed that larger tokamaks with reduced elongation and/or increased minor radius exhibited a significant increase in \( w^*_m \), attributed primarily to a decrease in bootstrap drive. The equilibria were generated using Monte Carlo methods, normalized to achieve over 400 MW of fusion power, and subjected to tearing analysis through the modified Rutherford equation for various resonant poloidal and toroidal modes up to \( n = 4 \).
The findings indicated that certain design configurations, particularly large and low-elongation tokamaks, could yield substantial improvements in \( w^*_m \), which is critical for preventing NTM growth. Machine learning techniques were employed to uncover underlying dependencies in \( w^*_m \), highlighting increased stability associated with higher triangularity due to lower \( \Delta’ \) values. These results suggest that it may be possible to achieve passive NTM stability in pilot plant operations through careful scenario design. Future research will aim to establish a workflow for experimentally validating \( w^*_m \) as a reliable NTM stability metric in tokamak design.
Introduction
The introduction discusses the neoclassical tearing mode (NTM), a resistive plasma instability that poses significant risks to the operation of tokamaks, particularly in the context of next-generation fusion pilot plants (FPPs) aiming for high-gain scenarios. NTMs are initiated by a plasma pressure gradient and require a ‘seed’ magnetic island for growth. Current stabilization methods, such as electron cyclotron current drive and off-axis neutral beam injection, are not feasible in the preliminary designs of high-temperature superconducting FPPs due to technological and economic constraints. The authors focus on the challenge of stabilizing NTMs in inductive, high-field, H-mode tokamaks without these common control methods.
To address this, the study employs the modified Rutherford equation (MRE) to identify conditions under which the size of seed islands remains too small to trigger NTM growth. By analyzing a synthetic equilibrium database of various tokamak designs, the authors find that the minimum marginally stable island width, denoted as $w^*_m$, is significantly influenced by geometric parameters such as major radius ($R_0$), elongation ($\kappa$), and minor radius ($a$). The results indicate that optimizing these parameters can enhance $w^*_m$, thereby improving passive NTM stability. The paper outlines its structure, detailing subsequent sections that will cover the tearing physics model, database generation methods, analysis results, and discussions on future work related to predicting NTM onset probabilities.
Results
In this section, the results of the tearing analysis conducted on the Monte Carlo H-mode database are presented. The analysis involved generating 14,651 ideal-MHD stable equilibria, which corresponded to 7,418 unique tokamak geometries. After applying resistive interchange criteria, 977 equilibria were excluded, leaving 13,719 scenarios for further analysis. Among these, 1,375 equilibria exhibited positive values of the stability parameter $\Delta’$, yet only 114 were found to be linearly unstable with $\Delta’ > \Delta’_c$. Additionally, 844 equilibria were deemed unsuitable for analysis due to their minimum marginally stable island widths being smaller than the linear layer width ($w^*_m < \delta_r$). The findings indicate that 2,136 equilibria, or 16.9% of the remaining cases, achieved full stabilization, characterized by a condition where $\frac{dw}{dt} < 0$ for all modes within the pedestal. In these fully stabilized scenarios, the minimum marginally stable island width was defined as the device size, with $w^*_m = 1$. The results highlight significant trends in the stability of the equilibria across the input space defined by parameters such as $R_0$, $\epsilon$, $\delta$, $\kappa$, and the normalized current $I/I^*$. Further discussions in the subsequent sections will elaborate on the physics underlying these trends and the dependencies of single-helicity toroidal $\Delta'$ across the database.
Discussion
In this section, the authors discuss the framework for analyzing the stability of axisymmetric plasma equilibria, focusing on the linear and nonlinear dynamics of tearing modes. The analysis is grounded in the equilibrium condition \( J \times B = \nabla p \), with the magnetic field represented in terms of poloidal magnetic flux \( \Psi \) and toroidal angle \( \phi \). The tearing stability is assessed using the resistive DCON code, which calculates perturbations in the plasma and identifies conditions under which these perturbations become unstable. The authors derive a set of equations governing the dynamics of single-helicity tearing modes, emphasizing the significance of the minimum marginally stable island width \( w^*_m \) as a critical metric for evaluating stability against neoclassical tearing modes (NTMs).
The findings indicate that the minimum marginally stable island width varies significantly across the synthetic H-mode database, with mean, maximum, and minimum values of \( 6.8 \times 10^{-3} \), \( 5.8 \times 10^{-2} \), and \( 9.4 \times 10^{-5} \), respectively. The analysis reveals that increasing the major radius and minor radius or adjusting elongation can enhance NTM stability, suggesting that larger or less performant tokamaks may be more favorable for operational stability. The authors also note that certain physics effects, such as ion polarization currents and differential rotation, are excluded from their analysis, arguing that their omission is justified given the focus on identifying favorable design regimes for future fusion pilot plants. Overall, the study highlights the complex interplay of equilibrium-dependent stabilizing and destabilizing effects on NTMs, providing insights for optimizing tokamak designs.