Superelastic Heating in Treanor-Gordiets Plasmas: A Unified Analytic Closure

ملخص تقني: التسخين المرن الفائق في بلازما تريانور-غورديتس

بيان المشكلة
في البلازما غير المتوازنة حرارياً، تخطئ النماذج التوافقية التقليدية بشكل كبير في تقدير معدلات تسخين الطاقة من الإلكترونات إلى الاهتزازات ($Q_{v-e}$) فيما يُعرف بالتسخين فائق المرونة (superelastic heating). تفشل هذه النماذج عندما تتجاوز درجة الحرارة الاهتزازية ($T_v$) درجة حرارة الغاز ($T_g$)، وهي حالة شائعة في مراحل ما بين النبضات في الاحتراق بمساعدة البلازما (PAC)، وفي ذيول إعادة دخول الغلاف الجوي للمركبات فائقة السرعة، وفي تفريغات التصنيع الكيميائي غير المتوازنة. تحت هذه الظروف، يتبع التوزيع الاهتزازي توزيع "تريانور" الذي يتميز بفرط التعداد في المستويات عالية الارتفاع. وتغفل النماذج التوافقية التقليدية، التي تفترض تباعداً متساوياً في الطاقة بين المستويات، مساهمة التسخين الأسية من هذه الحالات غير المتوافقة (anharmonic) المفرطة التعداد، مما يؤدي إلى تقديرات أقل من حقيقة معدلات التسخين بعوامل تصل إلى خمسة أضعاف تقريباً. وعلى العكس من ذلك، عندما تكون $T_g > T_v$، تبالغ النماذج التوافقية في تقدير التسخين بنسبة 15% تقريباً. علاوة على ذلك، اقتصرت النماذج الموجودة المتسقة ديناميكياً حرارياً على أنظمة درجات الحرارة المنخفضة أو الانتقالات أحادية الكم، وفشلت في رصد انتقالات التوافقات العليا متعددة الكم والتشوهات غير المتوافقة الموجودة في التفريغات عالية الطاقة.

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

1. صياغة الطاقة غير المتوافقة: يتم وصف مستويات الطاقة الاهتزازية باستخدام توسيع "دونهام" المبتور من الدرجة الثانية. ينتج عن ذلك تعبير تحليلي لفجوة الطاقة ($\Delta E_{n,m}$) بين المستويات الاهتزازية $n$ و $n+m$، والذي يتضمن دالة عيب عدم التوافق $\delta(n,m)$ التي ترصد الانحرافات عن التوافق حتى الدرجة الثانية.
2. النمذجة الحركية للتعدادات: يفترض النموذج توزيع "تريانور" للحالات الاهتزازية لمراعاة فرط التعداد عندما تكون $T_v > T_g$. ولمعالجة التباعد الفيزيائي لتوزيع "تريانور" الصرف عند المستويات العالية، تدمج الصيغة نسبة "غورديتس" المصححة بالتدفق ($G(n,m)$). وتعمل هذه النسبة على نمذجة منطقة الهضبة وراء حد أدنى لـ "تريانور" ($n^*$)، حيث يتفوق استرخاء الاهتزاز-الترجمة (V-T) على عمليات الضخ العلوي للاهتزاز-الاهتزاز (V-V).
3. نسبة التعداد الموحدة: يتم بناء نسبة تعداد مستمرة، خالية من القطع (piecewise-free)، عبر قطع نظام "تريانور" ديناميكياً عند نقطة العبور الحرجة $n^*$. ويتضمن ذلك قفزة كمية فعالة ($m^*$) تنقل النظام من منحنى "تريانور" إلى هضبة "غورديتس".
4. اشتقاق عامل التصحيح: من خلال تطبيق التوازن التفصيلي على فجوات الطاقة غير المتوافقة ونسب التعداد الموحدة، يشتق المؤلف عامل تصحيح عدم توافق معمّم ومغلق الشكل ($\Phi_{anh}^{(n,m)}$). يقوم هذا العامل بتعديل حد التسخين التوافقي القياسي لمراعاة المنافسة الحركية بين الضخ العلوي (V-V up-pumping) والاسترخاء (V-T relaxation).
5. الفصل العياني: لتمكين التنفيذ الفعال في الحلول المائعية، يتم اشتقاق صيغة مفككة (decoupled formulation). تفترض هذه الصيغة أن معامل معدل التصادم الإلكتروني المباشر مستقل عن المستوى الاهتزازي الأولي، مما يسمح بفصل مقاطع التصادم العرضية عن الحالة الديناميكية الحرارية للمجموعة الاهتزازية. وهذا يقلل من العبء الحسابي من مصفوفة كاملة من حالة إلى حالة (StS) إلى مجموعة من المعدلات الخاصة بكل قناة.

المساهمات الرئيسية

• إغلاق تحليلي موحد: تقدم الورقة معادلة حاكمة واحدة لتبادل الحرارة بين الإلكترونات والاهتزازات تظل صالحة عبر كل من الأنظمة منخفضة الحرارة (الشبه بولتزمانية) والأنظمة عالية الحرارة (تريانور-غورديتس).
• الاتساق الديناميكي الحراري: تفرض الصيغة بصرامة القانون الثاني للديناميكا الحرارية. وقد ثبت أنه عندما تكون $T_e = T_v = T_g$، فإن صافي انتقال الطاقة يتلاشى، ويساوي معدل التسخين تماماً معدل التبريد ($Q_{v-e} = Q_{e-v}$)، وهو شرط غالباً ما تنتهكه النماذج السابقة.
• عامل تصحيح عدم التوافق: يجسد إدخال $\Phi_{anh}$ فيزيائية التشوه غير المتوافق بشكل صريح، مما يصحح نقص تقدير التسخين في الأنظمة التي تكون فيها $T_v > T_g$، والمبالغة في تقدير التسخين في الأنظمة التي تكون فيها $T_g > T_v$.
• الكفاءة الحسابية: يحقق النموذج العياني المفكك انخفاضاً في تقييم المعدلات يتراوح بين 40 إلى 70 ضعفاً مقارنة بمرجعيات (StS) الكاملة مع الحفاظ على دقة عالية (بانحراف لا يتجاوز 15% عن حلول StS الدقيقة).

النتائج
تم التحقق من صحة النموذج مقابل مراجع الجمع المباشر من حالة إلى حالة (StS) لبلازما النيتروجين ($N_2$) وأول أكسيد الكربون ($CO$) عبر نظامين متميزين:

1. نظام الغاز البارد ($T_g = 300$ كلفن، $T_e = 2.0$ إلكترون فولت): في هذا النظام، الذي يميز مراحل ما بين النبضات في الاحتراق بمساعدة البلازما، يمكن أن تتجاوز $T_v$ درجة حرارة $T_g$ بشكل كبير. لقد قلل النموذج التوافقي من معدل التسخين بخمسة أضعاف. وقد نجح النموذج غير المتوافق المقترح، مع $\Phi_{anh} > 1$، في استعادة مرجع StS بدقة من خلال مراعاة الفائض فائق المرونة الناتج عن المستويات عالية الارتفاع المفرطة التعداد.
2. نظام الغاز الساخن ($T_g = 20,000$ كلفن، $T_e = 0.4$ إلكترون فولت): في هذا النظام، النموذجي لتسخين الصدمات في السرعات الفائقة، تكون $T_g > T_v$. بالغ النموذج التوافقي في تقدير التسخين بنسبة 15% تقريباً. وقد قام النموذج المعمم بتقليص تدفق التسخين بشكل صحيح ($\Phi_{anh} < 1$) ليتطابق مع مرجع StS.
3. أداء النموذج المفكك: طابقت نتائج النموذج العياني المفكك، الذي يستخدم مستويات تمثيلية موزونة ديناميكياً، نتائج الجمع الدقيقة، مما أثبت إمكانية التقاط فيزياء عدم التوافق دون الحاجة إلى حل كل انتقال فردي من حالة إلى حالة.

الأهمية
تزعم الورقة أن هذا الإطار الموحد يوفر معادلة حاكمة لتبادل الحرارة بين الإلكترونات والحالات المثارة في البيئات غير المتوازنة، بما في ذلك الاحتراق بمساعدة البلازما والتدفقات فائقة السرعة. ومن خلال معالجة قيود النماذج التوافقية التقليدية وضمان الاتساق الديناميكي الحراري، يتيح هذا العمل تطوير نماذج مختزلة دقيقة ومحدودة المعدل للحلول المائعية العيانية. وهذا أمر مهم بشكل خاص للتطبيقات التي تكون فيها بيانات الحركة (StS) الكاملة غير مكتملة أو مكلفة حسابياً، مما يوفر مساراً قوياً نحو محاكاة التفاعلات المعقدة بين البلازما والغاز بدقة عالية وتكلفة حسابية منخفضة.