Update-Magnitude State Redistribution (UM-SRD): A Shut-off Extension of Weighted SRD for Cut-Cell Methods

ملخص تقني: إعادة توزيع الحالة بناءً على مقدار التحديث (UM-SRD)

بيان المشكلة
تسمح طرق الخلايا المقطوعة (Cut-cell methods) لخطط الحجم المحدود بالتعامل مع الهندسات المعقدة دون الحاجة إلى شبكات متوافقة مع الأجسام، لكنها تسبب تحديات استقرار شديدة بسبب الخلايا المقطوعة الصغيرة بشكل عشوائي. تعالج طرق إعادة توزيع الحالة (SRD)، وتحديداً إطار عمل الـ SRD الموزون الذي طوره بيرجر وجولياني (2024)، مشكلات الاستقرار هذه عن طريق دمج الخلايا الصغيرة في جيرانها، وتطوير المتوسطات على الشبكة المدمجة، ثم إعادة توزيع النتائج مرة أخرى على الشبكة الأصلية.

ومع ذلك، فإن هناك قصوراً حرجاً في صيغة الـ SRD الموزون الحالية يتمثل في غياب آلية "الإيقاف" (shut-off). فخلافاً لطرق إعادة توزيع التدفق الكلاسيكية، حيث تتلاشى عملية ما بعد المعالجة طبيعياً عندما تكون اختلالات التدفق صفراً، يطبق الـ SRD الموزون مؤثر إعادة توزيع غير بديهي عند كل خطوة زمنية، بغض النظر عن مقدار تحديث الحجم المحدود. يؤدي هذا التفعيل المستمر إلى قصورين عمليين:

1. الفشل في الحفاظ على الحالات المستقرة: حتى عندما يكون توازن التدفق صفراً تماماً، يقوم الـ SRD باستبدال القيم النقطية بمتوسطات موزونة حجمياً، مما يمنع التقارب إلى دقة الآلة ويؤدي إلى أخطاء من رتبة $O(1)$ في المناطق الساكنة.
2. التبدد غير الضروري: في المناطق الناعمة ذات التحديثات الصغيرة ولكن غير الصفرية، يؤدي التطبيق المستمر لإعادة التوزيع إلى إدخال انتشار عددي يتناسب مع تدرج الجوار، وبشكل مستقل عن الديناميكيات الفيزيائية الفعلية.

المنهجية
يقدم هذا البحث إعادة توزيع الحالة بناءً على مقدار التحديث (UM-SRD)، وهو امتداد للإيقاف لتقنية الـ SRD الموزون، مصمم لدمج مؤثر الـ SRD مع مؤثر الوحدة بناءً على مقدار التحديث المحلي للحجم المحدود.

تعمل الطريقة ضمن إطار عمل الخطوات الزمنية القياسي المكون من ثلاث مراحل:

1. تحديث الحجم المحدود: حساب الحالة ما قبل إعادة التوزيع $U^*$ باستخدام المخطط الأساسي.
2. تجميع الجوار: حساب متوسطات الجوار الموزونة $\hat{Q}^*_j$ على الجوارات المدمجة $M_j$.
3. إعادة التوزيع المدمج: بدلاً من تطبيق مؤثر الـ SRD $S_j$ مباشرة، يطبق UM-SRD مزيجاً محدباً:
   $$R_j(s_j) = (1 - s_j)I + s_j S_j$$
   حيث $I$ هو مؤثر الوحدة و $s_j \in [0, 1]$ هو معامل الدمج.

يتم تحديد معامل الدمج $s_j$ بواسطة دالة إيقاف ناعمة محددة محلياً $h(\eta_j)$:
$$s_j = \frac{\eta_j^p}{\eta_j^p + \tau^p}$$
حيث $\eta_j$ هو مؤشر معياري لأقصى مقدار مطلق للتحديث داخل الجوار $M_j$ (تحديداً $\eta_j = \frac{\max \|\Delta U_i\|}{\epsilon + \max \|\Delta U_k\|}$).

• نظام التحديث الصفري: إذا كانت جميع التحديثات في الجوار عند دقة الآلة ($\Delta U_i \approx 0$)، فإن $\eta_j = 0$، وبالتالي $s_j = 0$، وتتحول الطريقة تماماً إلى مؤثر الوحدة ($U^{n+1} = U^*$).
• النظام النشط: إذا حدثت تحديثات كبيرة (على سبيل المثال، بالقرب من الانقطاعات أو عندما يكون تيبس الخلية الصغيرة نشطاً)، فإن $\eta_j \approx 1$، وبالتالي $s_j \approx 1$، وتتصرف الطريقة كـ SRD موزون كامل.

المساهمات التحليلية الرئيسية
يقدم البحث تحليلاً دقيقاً لـ UM-SRD لنموذج حمل تيار خطي أحادي الأبعاد مع خلية مقطوعة واحدة:

• الحفظ: ثبت أن المؤثر المدمج محافظ على الكتلة/الزخم بالبناء لأي $s_j \in [0, 1]$.
• تقليل التباين الكلي (TVD): تحت نفس شرط CFL للمخطط التصاعدي من الدرجة الأولى، ثبت أن UM-SRD هو TVD على شبكة الخلية المقطوعة الأصلية. ويعتمد هذا على المزيج المحدب بين مؤثر الوحدة ومؤثر الـ SRD الذي يحقق خاصية TVD.
• الحفاظ على الحالة المستقرة: ثبت أن UM-SRD يحافظ بدقة على الحالات المستقرة للمخطط الأساسي. إذا كان تحديث الحجم المحدود صفراً، فإن $s_j=0$، ولا تحدث أي إعادة توزيع، مما يتجنب الانجراف الزائف المتأصل في الـ SRD القياسي.
• الدقة: في المناطق الناعمة، يُظهر UM-SRD أنه اضطراب للمخطط الأساسي. إذا تناسب معامل الدمج مع حجم الشبكة (أو عبر المؤشر المعياري المستخدم في الممارسة العملية)، فإن المساهمة الإضافية في خطأ القطع المحلي تكون من رتبة أعلى من مخطط الدرجة الأولى الأساسي. وبالتالي، يحتفظ UM-SRD بمعدل تقارب الدرجة الأولى للمنهج التصاعدي الأساسي.

النتائج العددية
تؤكد التجارب العددية على مسائل نموذجية أحادية وثنائية الأبعاد التنبؤات النظرية:

• التقارب: يظهر UM-SRD تقارباً من الدرجة الأولى في معايير $L_1$ و $L_\infty$ للحمل الناعم، مع أخطاء لا يمكن تمييزها عن المخطط الأساسي والـ SRD الموزون القياسي (الفروقات $< 0.01\%$).
• سلوك الإيقاف: في اختبارات النبضات المدمجة، ينخفض معامل الدمج $s_j$ إلى الصفر تماماً عند غياب النبضة، مستعيداً المخطط الأساسي، ويرتفع إلى $\approx 0.99$ عندما تتفاعل النبضة مع الخلية الصغيرة.
• الاستقرار: يعمل UM-SRD على تثبيت المخطط الأساسي تحت ظروف خطوات زمنية للشبكة الكاملة (حيث يكون شرط CFL المحلي عند الخلية الصغيرة 10 أضعاف حد الاستقرار)، مما يمنع التباعد الأسي الملاحظ في المنهج التصاعدي الأساسي.
• الخصائص غير التذبذبية: تظهر الاختبارات التي تحتوي على بيانات منقطعة أن UM-SRD يظل غير تذبذبي، مطابِقاً لسلوك الـ SRD القياسي بالقرب من الخلية الصغيرة.
• الحفاظ على الحالة المستقرة: تؤكد التجارب أنه بينما يسبب الـ SRD القياسي أخطاء من رتبة $O(1)$ في الملفات التعريفية الساكنة، يحافظ UM-SRD على الملف التعريفي الأولي بدقة الآلة.
• التوسع لثنائي الأبعاد: تم تطبيق الخوارزمية بنجاح على شبكة ديكارتية ثنائية الأبعاد مع حدود مدمجة، مما أظهر تقارباً من الدرجة الأولى دون تعديل في منطق الدمج.

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