مقارنة نماذج استجابة البلازما لتأثيرات RMP على المشتت وطبقة الخدش في KSTAR Comparison of plasma response models for RMP effects on the divertor and scrape-off layer in KSTAR

شارك:
المجلة: Nuclear Fusion، المجلد: 66، العدد: 7
DOI: https://doi.org/10.1088/1741-4326/ae6b9e
تاريخ النشر: 2026-05-11
المؤلف: H. Frerichs وآخرون
الموضوع الرئيسي: أبحاث الاندماج المغناطيسي المحتجز

نظرة عامة

تناقش هذه القسم تأثير الاضطرابات المغناطيسية الرنانة (RMPs) على الأنماط المحلية الحادة (ELMs) في التوكاماك، مع التركيز بشكل خاص على آثارها على أحمال الجسيمات والحرارة على أهداف المشتت. تكشف الدراسة، التي أجريت على بلازما H-mode المثبطة لـ ELM في KSTAR، عن تباين كبير في بصمات المغناطيس المحسوبة—تتراوح من 2 سم إلى 14 سم—استنادًا إلى نماذج استجابة بلازما مختلفة (GPEC، MARS-F، M3D-C1، وJOREK). يتم عكس هذا التباين في محاكاة EMC3-EIRENE، التي تشير إلى أن القيم القصوى ومدى الخطوط الحلزونية غالبًا ما يتم تقديرها بشكل مبالغ فيه مقارنةً بقياسات التصوير الحراري بالأشعة تحت الحمراء (IRTV). تشير النتائج إلى أن إدخال الطاقة الأولية إلى طبقة الخدش قد يتم تقديره بشكل مبالغ فيه، مع ملاحظة توافق أفضل تحت ظروف انخفاض إدخال الطاقة والنقل العرضي أو كثافة أعلى في مجرى المياه وفقدان إشعاعي أكبر.

تؤكد الاستنتاجات على الحاجة إلى تحسين جودة البيانات في التجارب المستقبلية لتدقيق المدخلات المتعلقة بكثافة الفاصل والنقل العرضي الشاذ. بالإضافة إلى ذلك، يُوصى بربط نموذج حدود البلازما مع نموذج القلب لتعزيز دقة إدخال الطاقة الصافية كشرط حدودي، خاصة بعد الترقية إلى مشتت من التنجستن. يبرز القسم التباينات بين نماذج استجابة البلازما، ويعزوها إلى الاختلافات في كيفية معالجة التوازن داخل وخارج الفاصل. بينما يوفر M3D-C1 أفضل تطابق لأحمال حرارة المشتت استنادًا إلى حسابات خطية، فإن المزيد من المعايير والتحليلات الحساسية لنماذج استجابة البلازما ضرورية. كما يشير النص إلى إمكانية استخدام JOREK-STARWALL لنهج حدودي متسق ذاتيًا، والذي لا يزال موضوعًا للبحث المستقبلي.

مقدمة

ت outlines المقدمة أهمية نظام الاحتواء العالي (H-mode) لتشغيل البلازما المحترقة في ITER، مع تسليط الضوء على التحدي الذي تطرحه الأنماط المحلية الحادة (ELMs) التي تتطلب خطط تحكم لمنع الأضرار الناتجة عن فقدان الطاقة السريع. تتضمن طريقة التحكم البارزة استخدام الاضطرابات المغناطيسية الرنانة (RMPs)، التي، على الرغم من فعاليتها في التوكاماك الحالية، تؤدي إلى أنماط خطوط حلزونية في أحمال الجسيمات والحرارة على أهداف المشتت. يمكن أن تؤثر هذه الأنماط، التي لوحظت في تجارب مثل تلك التي في DIII-D، على كفاءة تبديد الطاقة، مما يتطلب توقعات دقيقة لأحمال حرارة المشتت تأخذ في الاعتبار استجابات البلازما للاضطرابات الخارجية.

تهدف الدراسة إلى مقارنة نماذج استجابة البلازما المختلفة—خطية وغير خطية، مثالية ومقاومة MHD، سائل واحد وسوائل مزدوجة—تحديدًا فيما يتعلق بآثارها على بصمات المشتت في بلازما KSTAR H-mode مع تثبيط ELM عبر RMPs. تشمل النماذج المستخدمة GPEC، MARS-F، M3D-C1، وJOREK، جميعها تستخدم هندسة حلقية كاملة. ستقوم الأبحاث بتحليل تأثير هذه النماذج على بصمات المغناطيس وأحمال الحرارة الناتجة على أهداف المشتت، مع مناقشة إعدادات المحاكاة التفصيلية والنتائج في الأقسام اللاحقة.

مناقشة

تركز قسم المناقشة من ورقة البحث على الهندسة المغناطيسية وآثارها على سلوك البلازما في تفريغ KSTAR 30306، تحديدًا عند 7850 مللي ثانية. تكشف التحليل أن تيار البلازما هو $I_p = 0.51$ MA، مع حقل مغناطيسي قدره 1.79 T وعامل أمان قدره $q_{95} = 5.24$. يتميز التكوين بانحياز نقطة X أقل وتطبيق الاضطرابات المغناطيسية الرنانة (RMPs) مع رقم وضع حلقي $n = 1$. يؤدي وجود RMPs إلى تشكيل سلاسل جزر مغناطيسية وسلوك خطي فوضوي، بينما يحدد أيضًا فاصلًا مغناطيسيًا يحدد خطوط المجال داخل وخارج حدود البلازما. تقارن الدراسة بين نماذج استجابة البلازما المختلفة (GPEC، MARS-F، JOREK، وM3D-C1) وتجد أن GPEC وMARS-F تنتجان بصمات مغناطيسية أكبر، بأحجام 14.4 سم و13.9 سم، على التوالي، بينما تنتج JOREK وM3D-C1 بصمات أصغر قدرها 4.9 سم و1.7 سم.

تشير محاكاة أحمال الحرارة على هدف المشتت الخارجي إلى أن أنماط الخطوط وأحمال الحرارة القصوى تختلف بشكل كبير عبر النماذج، حيث تظهر MARS-F أكبر الأنماط. ومع ذلك، جميع النماذج تقدر أحمال الحرارة القصوى بشكل مبالغ فيه مقارنةً بقياسات التصوير الحراري بالأشعة تحت الحمراء (IRTV)، مما يشير إلى تقديرات محتملة مبالغ فيها في إدخال الطاقة إلى طبقة الخدش (SOL). تشير النتائج إلى أن أنماط أحمال الحرارة تتأثر بالنقل العرضي وكثافة البلازما، مع زيادة الكثافات العليا في مجرى المياه تؤدي إلى زيادة الفقدان الإشعاعي. يؤكد المؤلفون على الحاجة إلى تحسين جودة البيانات في التجارب المستقبلية لتدقيق معلمات الإدخال وتعزيز دقة المحاكاة. كما يبرزون أهمية ربط نماذج حدود البلازما مع نماذج القلب لتحقيق تمثيل أكثر دقة لإدخال الطاقة وآثار RMPs على استجابة البلازما. بشكل عام، تؤكد النتائج على تعقيد نمذجة سلوك البلازما وضرورة المزيد من المعايير لنماذج استجابة البلازما.

Journal: Nuclear Fusion, Volume: 66, Issue: 7
DOI: https://doi.org/10.1088/1741-4326/ae6b9e
Publication Date: 2026-05-11
Author(s): H. Frerichs et al.
Primary Topic: Magnetic confinement fusion research

Overview

The section discusses the impact of resonant magnetic perturbations (RMPs) on edge localized modes (ELMs) in tokamaks, particularly focusing on their effects on particle and heat loads on divertor targets. The study, conducted on an ELM suppressed H-mode plasma at KSTAR, reveals significant variability in the computed magnetic footprints—ranging from 2 cm to 14 cm—based on different plasma response models (GPEC, MARS-F, M3D-C1, and JOREK). This variability is mirrored in EMC3-EIRENE simulations, which indicate that the peak values and extent of helical striations are often overestimated compared to infrared thermography (IRTV) measurements. The findings suggest that the initial power input into the scrape-off layer may be overestimated, with better alignment observed under conditions of lower input power and cross-field transport or higher upstream density and radiative losses.

The conclusions emphasize the need for improved data quality in future experiments to refine inputs related to separatrix density and anomalous cross-field transport. Additionally, coupling the plasma boundary model with a core model is recommended to enhance the accuracy of net power input as a boundary condition, especially following the upgrade to a tungsten divertor. The section highlights discrepancies among plasma response models, attributing them to differences in how equilibrium is treated within and beyond the separatrix. While M3D-C1 provides the best match for divertor heat loads based on linear, single-fluid calculations, further benchmarking and sensitivity analyses of plasma response models are necessary. The text also notes the potential for using JOREK-STARWALL for a self-consistent boundary approach, which remains a topic for future research.

Introduction

The introduction outlines the significance of the high confinement (H-mode) regime for burning plasma operation in ITER, highlighting the challenge posed by edge localized modes (ELMs) that necessitate control schemes to prevent damage from rapid energy loss. A prominent control method involves the use of resonant magnetic perturbations (RMPs), which, while effective in current tokamaks, induce helical striation patterns in particle and heat loads on divertor targets. These patterns, observed in experiments such as those at DIII-D, can influence power dissipation efficiency, necessitating accurate predictions of divertor heat loads that account for plasma responses to external perturbations.

The study aims to compare various plasma response models—linear and non-linear, ideal and resistive MHD, single fluid and two fluid—specifically regarding their effects on divertor footprints in a KSTAR H-mode plasma with ELM suppression via RMPs. The models employed include GPEC, MARS-F, M3D-C1, and JOREK, all utilizing full toroidal geometry. The research will analyze the impact of these models on the magnetic footprints and resulting heat loads on divertor targets, with detailed simulation setups and results discussed in subsequent sections.

Discussion

The discussion section of the research paper focuses on the magnetic geometry and its implications for plasma behavior in KSTAR discharge 30306, specifically at 7850 ms. The analysis reveals that the plasma current is $I_p = 0.51$ MA, with a magnetic field of 1.79 T and a safety factor of $q_{95} = 5.24$. The configuration is characterized by a lower X-point bias and the application of resonant magnetic perturbations (RMPs) with toroidal mode number $n = 1$. The presence of RMPs leads to the formation of magnetic island chains and chaotic field line behavior, while also establishing a magnetic separatrix that delineates field lines inside and outside the plasma boundary. The study compares different plasma response models (GPEC, MARS-F, JOREK, and M3D-C1) and finds that GPEC and MARS-F produce larger magnetic footprints, with sizes of 14.4 cm and 13.9 cm, respectively, while JOREK and M3D-C1 yield smaller footprints of 4.9 cm and 1.7 cm.

The simulations of heat loads on the outer divertor target indicate that the striation patterns and peak heat loads vary significantly across the models, with MARS-F showing the largest patterns. However, all models overestimate peak heat loads compared to infrared thermography (IRTV) measurements, suggesting potential overestimations in the power input to the scrape-off layer (SOL). The results indicate that the heat load patterns are influenced by cross-field transport and the density of the plasma, with higher upstream densities leading to increased radiative losses. The authors emphasize the need for improved data quality in future experiments to refine input parameters and enhance the accuracy of the simulations. They also highlight the importance of coupling plasma boundary models with core models to achieve a more accurate representation of power input and the effects of RMPs on plasma response. Overall, the findings underscore the complexity of modeling plasma behavior and the necessity for further benchmarking of plasma response models.