تأثير الانكماش في البلازما: كيف يقوم التوكاماك بسحب الجسيمات نحو المركز ذاتيًا
السياق
تقع هذه الدراسة في موقع متقدم بين أبحاث الاضطراب والنقل في بلازما حافة التوكاماك. معظم الدراسات المجاورة تبحث في النقل المضطرب للحرارة والجسيمات نحو الخارج، مع تثبيت الكثافة عند الحدود الداخلية أو استخدام نماذج مبسطة. هنا، ولأول مرة في محاكاة سائلة شاملة مع ملف كثافة متوافق ذاتيًا، يُظهر أن الانكماش ينشأ طبيعيًا دون افتراضات مصطنعة. وهذا يربط العمل بأسئلة قديمة في فيزياء الحبس، مثل تشكيل قاعدة الكثافة في نمط H.
المنهجية
استُخدم كود GDB، وهو نموذج كهرومغناطيسي ثلاثي الأبعاد شامل مبني على معادلات براغينسكي المختزلة بالانجراف. حُوكت بلازما بمعلمات قريبة من تلك الخاصة بـ Alcator C-Mod: نصف قطر كبير 68.5 سم، ونصف قطر صغير 22 سم، ومجال 5.3 تسلا. غطت منطقة الحساب الطبقة الحلقية الطورية بأكملها من 5.3 سم داخل آخر سطح مغناطيسي مغلق إلى 3.3 سم في طبقة الكشط. كان ملف الكثافة الابتدائي مسطحًا، ومصدر الجسيمات موضعًا عند الفاصل، بينما مصادر الحرارة عند الحدود الداخلية. استمرت المحاكاة 8.4 مللي ثانية من الزمن النقلي، مما سمح لملف الكثافة بالوصول إلى حالة شبه مستقرة. لتحليل الآليات، أُجريت حسابات إضافية مع تعطيل الموصلية الحرارية الإلكترونية والتدفقات الحرارية العرضية.
النتائج
أظهرت المحاكاة أنه من ملف ابتدائي مسطح يتشكل خلال بضعة مللي ثانية ملف بذروة كثافة في مركز منطقة الحساب. يوجد تدفق جسيمات نحو الداخل طوال التطور، لكن طبيعته تتغير. في المرحلة المبكرة، عندما يكون تدرج الكثافة موجبًا (الكثافة تزداد نحو الخارج)، يهيمن التدفق المضطرب المرتبط بالانتشار الحراري الإلكتروني. أظهر التحليل الخطي لنموذج مبسط أنه عندما تكون نسبة تدرجي الكثافة ودرجة الحرارة سالبة (η<0)، تُحدث الموجات الانجرافية إزاحة طورية بين تذبذبات الكثافة والجهد، مما يعطي تدفقًا سالبًا (نحو الداخل). في المرحلة المتأخرة، عندما يصبح ملف الكثافة مسطحًا، يضعف المساهمة المضطربة، لكن ينشأ تدفق E×B توازني ناتج عن عدم التماثل القطبي في الكثافة والجهد. هذا عدم التماثل يولده تدفق حرارة بفيرش-شلوتر، لكن النظرية الكلاسيكية الجديدة (neoclassical) تتنبأ بتدفق صافٍ صفري بسبب استجابة بولتزمان للإلكترونات. في نموذج GDB، تُؤخذ في الاعتبار تأثيرات غير مثالية (المقاومة، قصور الإلكترونات، الحث الكهرومغناطيسي) والتي تُحدث إزاحة طورية قطبية صغيرة δs ≈ −0.08π بين الكثافة والجهد. تؤدي هذه الإزاحة إلى تدفق داخلي محدود يفسر تراكم الكثافة كميًا. أظهرت حسابات مراقبة مع تعطيل التدفقات الحرارية العرضية أن ذروة الكثافة تتبدد، مما يؤكد دور عدم التماثل القطبي.
الأهمية
تختلف هذه الدراسة عن الدراسات المجاورة في أنها لا تفرض ملف الكثافة من الخارج، بل تسمح له بالتطور بشكل متوافق ذاتيًا. هذا هو أول عرض في إطار نموذج سائل شامل للحافة أن انكماش الجسيمات ينشأ طبيعيًا ويؤدي إلى تركيز الكثافة. لقد تبين أن الصورة الكلاسيكية الجديدة غير مكتملة: إن أخذ اللاتماثلية في قانون أوم بعين الاعتبار يفتح قناة نقل إضافية. هذا يغير فهم توازن الجسيمات في بلازما الحافة ويوفر أساسًا لتنبؤات أكثر دقة لملفات الكثافة.
التطور المستقبلي
يمكن أن يتطور الموضوع نحو أخذ هندسة أكثر واقعية في الاعتبار (مثل الاستطالة، والمحوّل) وتأثيرات حركية للبلازما الساخنة. من المهم التحقق مما إذا كان انكماش E×B غير المثالي المكتشف سيستمر في النماذج الجيروكينتيكية وعند معدلات تصادم أقل. كما أن دراسة دور التأين القطبي غير المتماثل للنيوترونات واعدة، والذي قد يعزز أو يضعف الانكماش في توكاماك حقيقي.
التأثير
ستؤثر النتائج على نمذجة حبس البلازما في التوكاماك والستيلاراتور، خاصة على التنبؤ بملفات الكثافة والانتقال إلى أنماط الحبس المحسّن. كما أنها مفيدة لتفسير تجارب حقن الغاز وحقن الكريات.
الخطوات التالية
تشمل الخطوات التالية مقارنة ملفات الكثافة التجريبية على Alcator C-Mod ومنشآت أخرى، بالإضافة إلى توسيع النموذج ليشمل التصحيحات الحركية والهندسة الحقيقية.
المسائل المفتوحة الرئيسية
ترتبط هذه الدراسة ارتباطًا مباشرًا بالمشكلة غير المحلولة المتعلقة بتشكيل قاعدة الكثافة في نمط H وبسؤال حول الآليات التي تحدد تركيز الكثافة في أنماط الحبس المختلفة. كما تلقي الضوء على لغز قديم: لماذا يفتقر نمط I إلى حاجز للجسيمات بينما يوجد حاجز للحرارة.
في المحاكاة، تسبح الجسيمات "ضد التيار": فهي تتحرك من الحافة نحو المركز، رغم أن تدرج الكثافة في البداية موجه نحو الخارج. هذا النقل غير البديهي ممكن بفضل الارتباط بين تذبذبات الكثافة والمجال الكهربائي، وليس الانتشار العادي.
فكرة الانضغاط الذاتي للبلازما نحو المركز تتردد مع مفهوم "الانكماش" في الخيال العلمي، كما في رواية لاري نيفن "عالم-الحلقة"، حيث تحصر المجالات المغناطيسية البلازما في مفاعل.
أرقام رئيسية
- مدة المحاكاة: 8.4 مللي ثانية
- إزاحة الطور δs: ≈ −0.08π
- ذروة الكثافة: 1.6×1020 م^-3
- معدل نمو الكثافة: خلال 4 مللي ثانية يصبح الملف مسطحًا، وخلال 8 مللي ثانية يصبح ذرويًا
- نسبة التدرجات η: من −2 إلى +2