The (Limited) Effect of Viscosity in Multiphase Turbulent Mixing

بيان المشكلة
تنتشر الغازات متعددة الأطوار في كل مكان في البيئات الفيزيائية الفلكية مثل الوسط بين النجمي (ISM)، والوسط المحيط بالمجرات (CGM)، والوسط داخل العناقيد المجريّة (ICM). وتُحكم ديناميكيات هذه الأنظمة من خلال التفاعل بين الاضطراب، والتبريد الإشعاعي، وعمليات النقل. وبينما تمت دراسة دور المجالات المغناطيسية في كبح عدم استقرار كيلفن-هلمهولتز (KHI) وتغيير الخلط في طبقات الخلط الإشعاعي المضطربة (TRMLs)، لا يزال التأثير المحدد للزوجة محل نقاش. ومن المتوقع أن تكون البلازما المتأينة بالكامل عالية اللزوجة (Spitzer 1962)، ومع ذلك فإن صافي اللزوجة الفعالة في الأنظمة الفيزيائية الفلكية لا يزال غير مؤكد. وهناك سؤال جوهري حول ما إذا كانت اللزوجة تكبح الاضطراب والخلط الضروريين لعلاقات القياس "العالمية" بين السرعة المضطربة ومعدل التبريد (Fielding et al. 2020; Tan et al. 2021)، لا سيما في الأنظمة التي يشير فيها رقم كارلوفيتز ($Ka < 1$) إلى تدفق صفيحي (laminar).

المنهجية
استخدم المؤلفون محاكاة عددية مثالية باستخدام كود Athena لعزل تأثيرات اللزوجة. وقد سارت الدراسة عبر مرحلتين رئيسيتين:

1. محاكاة ثنائية الأبعاد ثابتة الحرارة (Adiabatic): لتحديد اللزوجة الحرجة ($\eta_{Crit}$) المطلوبة لكبح عدم استقرار كيلفن-هلمهولتز (KHI)، قام المؤلفون بمحاكاة تدفقات قص مع تغيير كثافة التجاوز ($\chi = \rho_{cold}/\rho_{hot}$) وأرقام ماخ ($M$). وقد استخدموا إعداد صفيحة مستوية لاختبار اعتماد كثافة التجاوز، وإعداد ورقة مستوية لاختبار اعتماد رقم ماخ، وذلك لتجنب آثار الحدود. كما اشتقوا تقديراً تحليلياً لـ $\eta_{Crit}$ بناءً على مقارنة الزمن الدوري للزوجة مع زمن نمو عدم استقرار KHI، باستخدام المتوسط التوافقي للزوجة الكينماتيكية لمراعاة اختلافات الزخم في السوائل متعددة الكثافة.
2. طبقات الخلط الإشعاعي المضطربة ثلاثية الأبعاد (TRMLs): باستخدام اللزوجة الحرجة المشتقة لكثافة تجاوز قدرها $\chi=100$، أجرى المؤلفون محاكاة ثلاثية الأبعاد لطبقات TRMLs مع لزوجة صريحة. وقاموا بتغيير كسر اللزوجة بالنسبة لـ $\eta_{Crit}$ (من 0% إلى 50%) ورقم دامكولر ($Da = t_{turb}/t_{cool}$) عن طريق تعديل طول موجة الاضطراب. سمح هذا باستكشاف كل من نظام التبريد الضعيف ($Da < 1$) ونظام التبريد القوي ($Da > 1$). تضمنت المحاكاة التبريد الإشعاعي (خوارزمية Townsend) ولكن دون توصيل حراري صريح، مما ضمن بقاء $Ka < 1$.

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

• تدرج اللزوجة الحرجة: تؤكد الدراسة أن اللزوجة الحرجة المطلوبة لكبح KHI تتدرج خطياً مع كثافة التجاوز ($\eta_{Crit} \propto \chi$) في حد كثافة التجاوز العالية. ويعود ذلك إلى أن الوسط البارد والكثيف يحمل معظم الزخم، مما يتطلب احتكاكاً أقوى لكبح عدم الاستقرار. كما تحقق المؤلفون من أنه عند أرقام ماخ العالية، تعمل تأثيرات الانضغاطية على استقرار KHI حتى بدون لزوجة، مما يحدد رقم ماخ حرج ($M_{Crit}$) يتجاوزه يتم فيه كبح عدم الاستقرار.
• نظام التبريد الضعيف ($Da < 1$): في هذا النظام، تكون أزمنة التبريد أطول من أزمنة الخلط المضطرب. وجد المؤلفون أن اللزوجة تكبح الاضطراب بقوة، مما ينقل التدفق من الحالة المضطربة إلى الحالة الصفيحية. وبناءً على ذلك، تنهار العلاقة العالمية بين السرعة المضطربة ($u'$) والسطوع السطحي ($Q$) المستمدة للسوائل عديمة اللزوجة. ومع ذلك، ومن ناحية غير متوقعة، يظل إجمالي اللمعان البولومتري غير متأثر إلى حد كبير. ويعزو المؤلفون ذلك إلى "تأثير عنق الزجاجة" حيث يحد زمن التبريد الطويل من فقدان الطاقة الإجمالي بغض النظر عن آلية الخلط (الانتشار مقابل الحمل الحراري).
• نظام التبريد القوي ($Da > 1$): في هذا النظام، تكون أزمنة التبريد أقصر من أزمنة اللزوجة. واكتشف المؤلفون أن للزوجة تأثيراً محدوداً على ديناميكيات النظام. فعلى الرغم من وجود اللزوجة، فإن عملية التبريد تعمل على شحذ تدرجات درجة الحرارة والسرعة عند الواجهة، مما يعادل فعلياً تأثير التنعيم الناتج عن اللزوجة. وهذا يسمح لـ KHI بالتطور والحفاظ على مستويات اضطراب مماثلة للحالات عديمة اللزوجة. ونتيجة لذلك، تظل علاقات القياس بين $u'$ و $Q$ قائمة، ويتصرف النظام فعلياً كأنه غير لزج.
• التنافس بين اللزوجة والتبريد: اشتق المؤلفون شرطاً تحليلياً ($Da > \zeta$، حيث $\zeta$ هو كسر اللزوجة الحرجة) لتحديد متى يهيمن التبريد على اللزوجة. وتتوافق نتائجهم العددية مع هذا التنبؤ. في نظام التبريد القوي، تؤدي زيادة كمية الغاز متوسط درجة الحرارة في عمليات التشغيل اللزجة (بسبب جبهة أعرض قليلاً) إلى زيادة طفيفة في السطوع السطحي (تصل إلى عامل $\sim 2$) مقارنة بالتشغيلات عديمة اللزوجة، لكن قوانين القياس الأساسية تظل قائمة.

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

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

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