Screening effect of plasma flow on the resonant magnetic perturbation penetration in tokamak based on two-fluid model

1. بيان المشكلة

يتناول هذا البحث قضية حرجة وهي أنماط التمزق النيوكلاسيكية (NTMs) في بلازما التوكماك، والتي تعد عائقاً رئيسياً أمام درجة حرارة قلب البلازما وسبباً في حدوث اضطرابات كبرى. تُدفع أنماط NTMs بفعل فقدان تيار "بوتستراب" (bootstrap current) داخل الجزر المغناطيسية. ويتمثل أحد التحديات الرئيسية في التحكم في هذه الأنماط في فهم التفاعل بين الاضطرابات المغناطيسية الرنينية (RMPs) والجزر المغناطيسية.

تحديداً، يبحث البحث في ظاهرة اختراق النمط (mode penetration)، حيث يتغلب اضطراب مغناطيسي خارجي (RMP) على حجب تدفق البلازما ليؤدي إلى إعادة اتصال مغناطيسي وتشكيل جزيرة مغناطيسية. وبينما استُخدمت نماذج أحادية السيولة لدراسة ذلك، إلا أنها غالباً ما تفشل في رصد الفيزياء المعقدة المتعلقة بدوران البلازما وتأثيرات نصف قطر لارمور المحدود. يهدف المؤلفون إلى توضيح:

• التمييز بين "إعادة الاتصال المدفوع" والاختراق الحقيقي للنمط في وجود تيار بوتستراب.
• الدور المحدد لكل من تدفق انجراف $E \times B$ مقابل تدفق الانجراف الدايامغناطيسي (diamagnetic drift) في عملية الحجب.
• الآليات الفيزيائية الكامنة وراء التذبذبات المرصودة في عرض الجزيرة تحت ظروف التدفق العالي.

2. المنهجية

تستخدم الدراسة كوداً مطوراً للقيم الأولية، وهو MDC (MHD@Dalian Code)، القائم على نموذج MHD ثنائي السيولة رباعي الحقول.

• المعادلات الحاكمة: يحل النموذج أربع معادلات غير خطية مترابطة لكل من:
  1. الدوامة ($U$)
  2. التدفق المغناطيسي المستعرض ($\psi$)
  3. ضغط البلازما ($p$)
  4. السرعة الأيونية الموازية ($v$)
• الفيزياء الأساسية المتضمنة:
  • تيار بوتستراب ($j_b$): وهو أمر بالغ الأهمية لزعزعة استقرار NTM.
  • تأثيرات النقل: كل من النقل الحراري الموازي ($\chi_{\parallel}$) والمتعامد ($\chi_{\perp}$).
  • تأثيرات السيولة المزدوجة: يحتفظ بالانجراف الدايامغناطيسي للإلكترونات وتدفق $E \times B$، مع مراعاة تأثيرات نصف قطر لارمور المحدود عبر المعامل $\delta$.
  • الهندسة: هندسة أسطوانية $(r, \theta, z)$ مع شروط حدودية دورية.
• الإعداد العددي:
  • تركز المحاكاة على سيناريو تفريغ يتم تسخينه أومياً بكثافة منخفضة ($n_e \approx 2 \times 10^{19} \text{ m}^{-3}$، $B_0 = 2 \text{ T}$).
  • تم ضبط السطح الرنيني عند $q = 3/2$ ($r \approx 0.4$).
  • تُطبق الاضطرابات المغناطيسية الرنينية (RMPs) كشروط حدودية لمحاكاة الاضطراب الخارجي.
  • تدرس الدراسة بشكل منهجي تغير سعة RMP، تردد دوران البلازما، كسر تيار بوتستراب ($f_b$)، المقاومة ($\eta$)، ومعاملات النقل.

3. المساهمات والنتائج الرئيسية

أ. إعادة تقييم عتبات اختراق النمط

• نمط التمزق الكلاسيكي ($f_b = 0$): في غياب تيار بوتستراب، لا توجد عتبة لاختراق النمط. يتناسب عرض الجزيرة المغناطيسية خطياً مع سعة RMP (إعادة اتصال مدفوعة).
• نمط التمزق النيوكلاسيكي ($f_b \neq 0$): عند تضمين تيار بوتستراب، لوحظت ظاهرة "تشبه اختراق النمط".
  • عند سعات RMP المنخفضة، تظل الجزيرة صغيرة.
  • بمجرد تجاوز RMP لعتبة حرجة، تنمو الجزيرة بسرعة لتصل إلى حجم كبير.
  • التمييز: يوضح المؤلفون أن هذه ليست عتبة واحدة، بل هي انتقال من نظام إعادة الاتصال المدفوع إلى نظام نمو NTM. وهذا يفسر سبب اكتشاف التجارب غالباً لعتبة اختراق محددة حتى عند دوران صفري، وهو ما قد تخطئ النماذج أحادية السيولة في تفسيره.

ب. تأثيرات الحجب: $E \times B$ مقابل الانجراف الدايامغناطيسي

تقارن الدراسة بين قدرات الحجب لتدفق $E \times B$ وتدفق الانجراف الدايامغناطيسي:

• عتبة الاختراق: يوفر كلا التدفقين تأثيرات حجب متشابهة فيما يتعلق بالعتبة المطلية لاختراق البلازما. تعتمد العتبة على فرق الدوران بين السطح الرنيني والـ RMP، بغض النظر عن نوع التدفق.
• الاستقرار: يظهر اختلاف جوهري في عرض الجزيرة المشبع.
  • تدفق $E \times B$: يؤدي إلى عروض جزيرة مشبعة أكبر.
  • تدفق الانجراف الدايامغناطيسي: يمارس تأثيراً استقرارياً قوياً، مما يؤدي إلى عروض جزيرة مشبعة أصغر بكثير مقارنة بتدفق $E \times B$ عند نفس تردد الدوران. ويُعزى ذلك إلى حد $\delta \nabla_{\parallel} p$ في قانون أوم المعمم.

ج. اكتشاف تذبذب عرض الجزيرة

في الأنظمة ذات أرقام لوندكويست ($S$) العالية ونسب النقل الموازي إلى المتعامد ($\chi_{\parallel}/\chi_{\perp}$) العالية، اكتشف المؤلفون تذبذباً دورياً في عرض الجزيرة المغناطيسية بعد الاختراق.

• الآلية: يُدفع التذبذب بواسطة حلقة تغذية راجعة سلبية بين الجزيرة المغناطيسية وضغط البلازما.
  1. مع نمو الجزيرة، يتسطح تدرج الضغط داخل الجزيرة، ولكن التدرج عند السطح الرنيني يتغير.
  2. هذا يغير تردد الانجراف الدايامغناطيسي ($\omega_{dia}$).
  3. يؤثر التغيير في التدفق على عرض الجزيرة عبر حدود السيولة المزدوجة.
  4. يتذبذب النظام حيث ينظم تدرج الضغط وعرض الجزيرة بعضهما البعض.
• التحقق: يختفي التذبذب عندما:
  • تزداد المقاومة ($\eta$) (الانتشار يخمد التأثير).
  • ينخفض كسر تيار بوتستراب ($f_b$).
  • ينخفض النقل الموازي ($\chi_{\parallel}$) أو يزدزد النقل المتعامد ($\chi_{\perp}$).

4. الأهمية والاستنتاج

يقدم هذا العمل رؤى بالغة الأهمية لتشغيل أجهزة التوكماك المتقدمة المستقبلية (مثل ITER وCFETR)، والتي تعتمد بشكل كبير على كسور عالية من تيار بوتستراب.

1. استراتيجيات التحكم: تشير النتائج إلى أنه بينما يمكن استخدام RMPs لقفل وتحديد مواقع NTMs من أجل كبتها بواسطة درع التيار الإلكتروني الحلقي (ECCD)، فإن الانجراف الدايامغناطيسي يلعب دوراً استقرارياً حيوياً يجب الاستفادة منه.
2. النمذجة الفيزيائية: تسلط الدراسة الضوء على ضرورة استخدام نماذج السيولة المزدوجة بدلاً من النماذج أحادية السيولة للتنبؤ بدقة بعتبات NTM وديناميكيات الجزيرة، لا سيما فيما يتعلق بالتأثير الاستقراري للتدفقات الدايامغناطيسية.
3. النوافذ التشغيلية: يساعد تحديد أنظمة التذبذب في تحديد الحدود التشغيلية حيث قد تتصرف الجزر المغناطيسية بشكل غير متوقع بسبب الاقتران بين الضغط والتدفق، مما يساعد في تصميم أنظمة قوية لتخفيف الاضطرابات.

باختصار، تثبت الورقة أن حجب تدفق البلازما ليس موحداً؛ فالانجراف الدايامغناطيسي أكثر فعالية في استقرار الجزر المغناطيسية من تدفق $E \times B$، وأن التفاعل بين تدرجات الضغط والتدفق يمكن أن يؤدي إلى سلوكيات تذبذبية معقدة في سيناريوهات التوكماك عالية الأداء.