What is the Strouhal number of turbulence driven by supernovae?

ملخص تقني: عدد ستروهال للتدفق الاضطرابي المدفوع بالمستعرات العظمى

بيان المشكلة
في نماذج الاضطراب القياسية المدفوعة بالزخم، ولا سيما محاكاة صندوق الاضطراب (TB) المثالية للموائع الفيزيائية الفلكية، يتم بناء تسارع الدفع العشوائي عادةً كعملية "أورنستين-أولين" (Ornstein–Uhlenbeck). ويعد "عدد ستروهال" ($St = t_{cor}/t_{out}$) بارامتراً حرجاً في هذا البناء، حيث يمثل النسبة بين زمن الارتباط لعملية الدفع ($t_{cor}$) وزمن دوران الدوامة عند المقياس الخارجي ($t_{out}$). إن الممارسة القياسية في نماذج صندوق الاضطراب هي ضبط $St = 1$، مما يعني أن عملية الدفع تفقد ارتباطها في نفس المقياس الزمني لأكبر الدوامات. ومع ذلك، فإن هذا الخيار هو إلى حد كبير درجة حرية غير مقيدة بحسابات تحليلية أو عددية للوسط بين النجمي (ISM) البلازمي القابل للانضغاط والواقعي. وقد أظهرت أعمال حديثة (Grete et al. 2025; Scanneria et al. 2025) أن إحصائيات الاضطراب، وتحديداً دالة توزيع الاحتمال (PDF) لكثافة الكتلة في الاضطراب المدفوع انضغاطياً، حساسة للغاية لهذا الخيار. ولا يزال من غير المعروف ما إذا كان $St=1$ يتحقق فيزيائياً في اضطراب الوسط بين النجمي المدفوع بالمستعرات العظمى (SN)، وإذا لم يكن كذلك، فما هي المقاييس التي يعمل عندها زمن ارتباط الدفع فعلياً.

المنهجية
يقوم المؤلفون بحساب $St$ مباشرة من محاكاة عالية الدقة، متدرجة الطبقات ومتعددة الأطوار للوسط بين النجمي في مجرات تشبه درب التبانة (MW) ومجرات انفجارية (SB)، والتي تم تحليلها سابقاً في (Connor et al. 2026). تستخدم هذه المحاكاة كود ramses لنمذجة اضطراب جاذبية-هيدروديناميكي مثالي مدفوع بالمستعرات العظمى، مع دمج شبكة تبريد تعتمد على الزمن لالتقاط الأطوار الدافئة، والمتأينة، والساخنة للوسط بين النجمي. تستبعد هذه المحاكاة الجاذبية الذاتية، والمجالات المغناطيسية، والأشعة الكونية لعزل تأثير دفع المستعرات العظمى.

تتضمن المنهجية ما يلي:

1. تحديد المقاييس: يتم اشتقاق المقياس الخارجي ($\ell_{out}$) من تكامل طيف طاقة السرعة، ويُحسب زمن دوران المقياس الخارجي ($t_{out}$) كـ $\ell_{out}/\langle u^2 \rangle^{1/2}_V$.
2. قياس الارتباط: بدلاً من فرض ارتباط الدفع، يستنتج المؤلفون زمن إلغاء ارتباط الدفع الفعال ($t_{cor}$) من ممتد ارتباط السرعة الأويلري (Eulerian velocity correlation tensor) المقاس زمنياً، $R_{ij}(x, t) = \langle \delta u_i(x, \tau) \delta u_j(x, \tau + t) \rangle_\tau$.
3. تحليل الممتد (Tensor Analysis): يقومون بحساب ممتد الارتباط الموحد ودمجه حتى أول عبور للصفر لتحديد ممتد زمن الارتباط المحلي $t_{cor}^{ij}(x)$. يسمح هذا بحساب مكونات عدد ستروهال المتناحية، والمُسقطة (موازية/عمودية على جهد الجاذبية)، والقطرية.
4. إعادة البناء المعتمد على المقياس: لاختبار النماذج الفيزيائية، يقوم المؤلفون ببناء عدد ستروهال يعتمد على المقياس $St(\ell)$، بافتراض أن زمن الارتباط يتحدد بواسطة آلية الدفع (ثابت)، بينما يتغير زمن دوران الدوامة مع المقياس وفقاً لطيف السرعة المقاس ($u_\ell \propto \ell^{1/4}$).

المساهمات والنتائج الرئيسية
تقدم الورقة أول قياس مباشر لـ $St$ في وسط بين نجمي متعدد الأطوار ومدفوع بالمستعرات العظمى. وتتلخص النتائج الرئيسية فيما يلي:

• عدد ستروهال العالمي: قيم عدد ستروهال المتوسطة المقاسة أقل بكثير من الوحدة. بالنسبة لنموذج درب التبانة (MW)، فإن الوسيط المتناحي هو $St = 0.26^{+0.30}_{-0.16}$، وبالنسبة لنموذج المجرة الانفجارية (SB)، هو $St = 0.25^{+0.11}_{-0.12}$. القيم المُسقطة (الموازية والعمودية على القرص) هي أيضاً دون الوحدة ($\sim 0.22–0.25$). يشير هذا إلى أن الدفع الفعال في اضطراب المستعرات العظمى يفقد ارتباطه بشكل أكبر بكثير من زمن دوران الدوامة عند المقياس الخارجي.
• التفسير الفيزيائي (نموذج نصف قطر التبريد): يفسر المؤلفون قيم $St$ المنخفضة عند المقياس الخارجي من خلال نموذج "نصف قطر التبريد/طبقة التلامس". في هذا الإطار، يتم حقن الاضطراب محلياً بالقرب من نصف قطر التبريد ($R_{cool}$) لبقايا المستعر الأعظم، حيث تولد طبقة التلامس غير المستقرة اضطراباً دورانيّاً (solenoidal). عند مقياس الحقن هذا ($\ell_{inj} \sim R_{cool}$)، يكون عدد ستروهال المحلي طبيعياً $St(\ell_{inj}) \approx 1$ لأن زمن نمو عدم الاستقرار يحدد كلاً من ارتباط الدفع وزمن دوران الدوامة.
• الاعتماد على المقياس: من خلال نشر شرط $St(\ell_{inj}) \approx 1$ عبر الشلال الخالي من المقياس (حيث يتبع طيف السرعة $P_u(k) \propto k^{-3/2}$)، يُتوقع الحصول على تدرج $St(\ell) \propto \ell^{-3/4}$. ويتحقق المؤلفون من هذا التنبؤ عبر إعادة بناء $St(\ell)$ من طيف السرعة. تصل المنحنيات المعاد بناؤها إلى $St=1$ عند مقياس لابعوي شبه عالمي $\ell^*/\ell_{out} \approx 0.12–0.13$.
• مقاييس الحقن الفعالة: تم استرداد المقاييس الفيزيائية حيث $St=1$ لتكون $\ell^*_{MW} \approx 25$ فرسخاً فلكياً و $\ell^*_{SB} \approx 32$ فرسخاً فلكياً. وتتوافق هذه القيم مع النطاق المعياري لنصف قطر تبريد بقايا المستعر الأعظم ($R_{cool} \sim 10–30$ فرسخاً فلكياً)، مما يدعم النموذج القائل بأن اضطراب المستعرات العظمى يتم دفعه محلياً عند طبقة التلامس الإشعاعية بدلاً من المقياس الخارجي العالمي.

الأهمية والادعاءات
تدعي الورقة أن المعيار القياسي لـ $St=1$ في نماذج صندوق الاضطراب لا يمثل فيزياء المقياس الخارجي لاضطراب الوسط بين النجمي المدفوع بالمستعرات العظمى. بدلاً من ذلك، يعد $St=1$ تقريباً للمقياس المحلي المرتبط بحقن الاضطراب بالقرب من نصف قطر التبريد للمستعرات العظمى.

تكمن أهمية هذه النتيجة في آثارها على نماذج الوسط بين النجمي (ISM)، والوسط المحيط بالمجرة (CGM)، والوسط بين العناقيد (ICM). وبما أن دالة توزيع الاحتمال لكثافة الكتلة ونماذج معدل تكوين النجوم حساسة للارتباطات الزمنية لقوة الدفع (كما أظهر Scannier et al. 2025)، فإن اكتشاف أن اضطراب المستعرات العظمى يعمل في نظام "الارتباط القصير" ($St \ll 1$) عند المقياس الخارجي يشير إلى أن النماذج الحالية التي تمت معايرتها باستخدام $St=1$ قد تتطلب إعادة تقييم. وتحديداً، تضع القيم المقاسة اضطراب المستعرات العظمى في جانب "الارتباط القصير" من الفضاء البارامتري، مما ينتج دوال توزيع كثافة أضيق وفراغات أقل كثافة مقارنة بالدفع ذي الارتباط الطويل. ويخلص المؤلفون إلى أن نموذج نصف قطر التبريد/طبقة التلامس يوفر تفسيراً مبرراً فيزيائياً لأعداد ستروهال المرصودة، وهو قابل للتطبيق على أي اضطراب مدفوع بموجات صدمية تبريد إشعاعي.