Multispecies Bhatnagar-Gross-Krook models and the Onsager reciprocal relations

ملخص تقني: نماذج باتكنار-غروس-كروك متعددة الأنواع وعلاقات أونساجر المتبادلة

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

المنهجية
يصيغ البحث النموذج الأكثر عمومية لنموذج BGK متعدد الأنواع لخليط ثنائي من الغازات أحادية الذرة. تتضمن هذه الصيغة ترددات التصادم ($\nu_{ij}$) ومعاملات ($\alpha_i, \beta_i, \gamma_i$) التي تحدد الاسترخاء نحو توزيعات ماكسويل المختلطة ($M_{ij}$) ذات السرعات ودرجات الحرارة المعتمدة على النوع. يضمن المؤلف حفظ الكتلة والزخم والطاقة من خلال فرض قيود محددة على هذه المعاملات.

لتحليل خصائص النقل، يستخدم البحث تقريب الانتشار بدلاً من التقريب الهيدروديناميكي القياسي. يتضمن ذلك:

1. التحجيم: إدخال معامل صغير $\epsilon$ لتمديد الإحداثيات المكانية ($\nabla \to \epsilon \nabla$) والزمن ($\partial_t \to \epsilon^2 \partial_t$)، وهو ما يناسب تدفقات الانتشار البطيئة والضعيفة.
2. التوسع: توسيع دالات التوزيع، والسرعات الجهرية، ودرجات الحرارة، والكثافات في قوى $\epsilon$.
3. الاشتقاق: تعويض هذه التوسعات في معادلات BGK المعاد تحجيمها لإيجاد التصحيح من الدرجة الأولى لدالة التوزيع ($f^{(1)}_i$).
4. حساب التدفق: دمج $f^{(1)}_i$ لاشتقاق تعبيرات صريحة لتدفقات الكتلة ($J_i$) والتدفق الحراري ($Q$) بدلالة تدرجات الكثافة ودرجة الحرارة.
5. المقارنة: مقارنة تدفقات BGK المشتقة مع التعبيرات الهيدروديناميكية القياسية (نظرية إنسكوج-تشابمان) التي تحقق علاقات أونساجر.

المساهمات والنتائج الرئيسية
النتيجة المركزية للبحث هي إثبات أن النسخة الأكثر شيوعًا من نموذج BGK متعدد الأنواع — حيث تكون ترددات التصادم ومعاملات النموذج مستقلة عن السرعة الجزيئية — لا تحقق علاقات أونساجر المتبادلة.

على وجه التحديد، يكشف التحليل عن عدم اتساق جوهري:

• وفقًا لعلاقات أونساجر، يجب أن يكون المعامل الذي يربط تدرج درجة الحرارة بتدفق الكتلة (الانتشار الحراري، $B$) غير صفري إذا كان المعامل الذي يربط تدرج الكثافة بالتدفق الحراري (تأثير دوفور، $C$) غير صفري، ويتم ربطهما بالعلاقة $C = pB$.
• اشتقاق نموذج BGK يعطي تدفقًا كتليًا يكون فيه حد الانتشار الحراري صفرًا تمامًا ($B=0$) ما لم تُتخذ حدود غير فيزيائية محددة.
• ومع ذلك، يحتوي التدفق الحراري الناتج على حد يتناسب مع تدرج الكثافة (تأثير دوفور) وهو غير صفري ويتناسب مع فرق الكتلة بين الأنواع.
• بناءً على ذلك، يتنبأ نموذج BGK بتأثير دوفور غير صفري إلى جانب تأثير سوريت صفري، مما ينتهك علاقة أونساجر $C = pB$.

يوضح البحث أيضًا أن هذا التناقض لا يمكن حله عن طريق معايرة معاملات النموذج (مثل المعامل $\beta$ الذي يحكم استرخاء السرعة النسبية). فلا توجد قيمة لـ $\beta$ تحقق في آن واحد معادلات النقل القياسية وعلاقات أونساجر للخلائط العامة. وحتى في حالة السوائل ذات الانتشار الحراري الضئيل (حيث $B \approx 0$)، يفشل نموذج BGK في إعادة إنتاج سلوك التدفق الحراري الصحيح دون فرض تلاشي التدفق المنتشر بأكمله.

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

يشير المؤلف إلى أنه بينما يفشل النموذج في تحقيق علاقات أونساجر بدقة، فإنه يحققها تقاربيًا في حالتين محددتين:

1. عندما تكون نسبة الكتلة $m_1/m_2 \to 0$ (وهو ما ينطبق على البلازما المتأينة حيث كتلة الإلكترون مهملة مقارنة بالأيونات).
2. عندما تكون نسبة الكتلة $m_1/m_2 \to 1$ (حيث يتلاشى الحد الأول في تعبير التدفق الحراري).

يختتم البحث بسرد متطلبات نموذج حركي "مثالي" متعدد الأنواع، والتي تشمل قوانين الحفظ، ونظرية H، ومبدأ عدم التمايز، وإيجابية درجة الحرارة/التركيز، والقدرة على تمثيل أرقام براندتل/شميدت التعسفية، والامتثال لعلاقات أونساجر. ويذكر المؤلف أنه حتى الآن، لم يثبت لأي نموذج من نوع BGK أنه يحقق جميع هذه المتطلبات في آن واحد. وتشير النتائج إلى أنه بينما تعد نماذج BGK الحالية تقريبات مفيدة، فإن عدم قدرتها على تحقيق علاقات أونساجر يحد من دقتها الفيزيائية وفائدتها في المعايرة الدقيقة لتطبيقات ديناميكا السوائل العامة.