محاكاة الميكروسكوبي-الماكروسكوبي لمعادلات العزوم للمياه الضحلة

المعلومات الأساسية

• معرّف الورقة: 2511.15737
• العنوان: محاكاة الميكروسكوبي-الماكروسكوبي لمعادلات العزوم للمياه الضحلة
• المؤلف: فيليم روژك (المشرفون: جوليان كويلرماير، ريك فيربيست)
• التصنيف: physics.flu-dyn, cs.NA, math.NA
• تاريخ النشر: 18 نوفمبر 2025
• رابط الورقة: https://arxiv.org/abs/2511.15737

الملخص

تُحكم تدفقات المياه الضحلة عادة بمعادلات نافييه-ستوكس، لكن النمذجة العملية غالباً ما تستخدم معادلات المياه الضحلة المبسطة (SWE)، مما يؤدي إلى نتائج غير دقيقة. لهذا السبب، طور الباحثون معادلات العزوم للمياه الضحلة (SWME)، وهي نموذج يستخدم معادلات ومتغيرات أكثر، مما يحسّن الدقة بشكل كبير، لكن التكلفة الحسابية تزداد أيضاً. تطبق هذه الورقة طريقة الميكروسكوبي-الماكروسكوبي (micro-macro) على SWME، مما يسمح بالتبديل بين نموذجين بمستويات دقة مختلفة، مما يتيح خطوات زمنية استقرار أكبر. تقدم الورقة تحليلاً نظرياً لوقت التشغيل وتعرض نتائج اختبارات انهيار السد ونقل الموجات. تحقق طريقة الميكروسكوبي-الماكروسكوبي تسريعاً حسابياً ملحوظاً مع الحفاظ على دقة كافية.

الخلفية البحثية والدافع

1. المشكلة المراد حلها

تواجه محاكاة تدفقات المياه الضحلة تناقضاً بين الدقة والكفاءة الحسابية:

• معادلات المياه الضحلة التقليدية (SWE): تفترض ملف السرعة الرأسي ثابتاً، سريعة الحساب لكن بأخطاء كبيرة عند تغير ملف السرعة
• معادلات العزوم للمياه الضحلة (SWME): تستخدم تقريباً متعدد الحدود لملف السرعة الرأسي، دقيقة لكن بتكلفة حسابية كبيرة

2. أهمية المشكلة

لمحاكاة تدفقات المياه الضحلة تطبيقات واسعة في عدة مجالات حرجة:

• التنبؤ بالطقس
• الهيدروليكا السطحية الحرة للأنهار والقنوات
• تقييم الانهيارات الجليدية والانزلاقات الأرضية
• حسابات نقل الجزيئات
• نمذجة عمليات الطلاء

3. قيود الطرق الموجودة

• SWE: الافتراض المتوسط العمقي يؤدي إلى عدم دقة عند تغير السرعة الرأسية
• SWME: زيادة المعاملات متعددة الحدود تزيد من عدد المتغيرات (N+2 معادلة)، مما يطيل وقت الحساب ويقيد خطوة الزمن
• طرق التقسيم الزمني الصريحة (مثل Forward Euler): رغم أن كل خطوة رخيصة حسابياً، لكن شروط الاستقرار صارمة وتتطلب خطوات زمنية كثيرة

4. الدافع البحثي

تم تطبيق طريقة الميكروسكوبي-الماكروسكوبي سابقاً فقط على نماذج العزوم للغازات النادرة، ولم تُطبق على تدفقات السطح الحر. تهدف هذه الورقة إلى توسيع هذه الطريقة إلى SWME لتسريع محاكاة تدفقات السطح الحر مع الحفاظ على دقة كافية.

المساهمات الأساسية

1. تطبيق أول لطريقة الميكروسكوبي-الماكروسكوبي على معادلات العزوم للمياه الضحلة، ملء الفراغ في تطبيق هذه الطريقة في مجال تدفقات السطح الحر
2. اشتقاق كامل للخطوات الأربع لطريقة الميكروسكوبي-الماكروسكوبي لـ SWME: الخطوة الميكروسكوبية، خطوة التقييد، الخطوة الماكروسكوبية وخطوة المطابقة
3. الاستفادة من التعامد متعدد الحدود Legendre، لتحقيق خطوة مطابقة بسيطة وفعالة (من خلال تقليل معيار L2)
4. تحليل نظري للتعقيد الحسابي: إثبات أنه عندما M>>L، يكون تعقيد طريقة الميكروسكوبي-الماكروسكوبي O((M²+L²+M)/(2Δx))، أفضل من النموذج الميكروسكوبي البحت O(M²/Δx)
5. التحقق العددي: تحقيق تسريع حسابي بمقدار 2 مرة أو أكثر في حالتي اختبار انهيار السد ونقل الموجات، مع الحفاظ على دقة جيدة
6. التنفيذ مفتوح المصدر: توفير مكتبة كود Python، باستخدام تقسيم الحجم المحدود من الدرجة الأولى والتنسيق المقسم

شرح الطريقة

تعريف المهمة

المدخلات: الشروط الابتدائية لتدفق المياه الضحلة (الارتفاع h، متوسط السرعة um، المعاملات متعددة الحدود αj) والشروط الحدية
المخرجات: حل حقل التدفق المتطور زمنياً
القيود: تحقيق حفظ الكتلة والزخم وتطور ملف السرعة، مع تحقيق شرط الاستقرار CFL

معمارية النموذج

1. أساس معادلات العزوم للمياه الضحلة (SWME)

يتم توسيع السرعة الرأسية باستخدام متعدد حدود: $u(t, x, \zeta) = u_m(t, x) + \sum_{j=1}^{N} \alpha_j(t, x)\phi_j(\zeta)$

حيث:

• $\zeta = \frac{z - h_b(t,x)}{h(t,x)} \in [0,1]$ هي الإحداثي الرأسي المعياري
• $\phi_j$ هي متعددات Legendre المزاحة على 0,1 (معايرة بحيث $\phi_j(0)=1$)
• N هي درجة متعدد الحدود، والنظام يحتوي على N+2 متغير

مجموعة المعادلات الحاكمة:

• حفظ الكتلة: $\partial_t h + \partial_x(hu_m) = 0$
• حفظ الزخم: $\partial_t(hu_m) + \partial_x[h(u_m^2 + \sum_{j=1}^{N}\frac{\alpha_j^2}{2j+1}) + \frac{g}{2}e_z h^2] = -\frac{\nu}{\lambda}(u_m + \sum_{j=1}^{N}\alpha_j) + hg(e_x - e_z\partial_x h_b)$
• تطور المعاملات (i=1,...,N): معادلات مقترنة معقدة تتضمن الثوابت $A_{ijk}, B_{ijk}, C_{ij}$

2. الخطوات الأربع لطريقة الميكروسكوبي-الماكروسكوبي

الخطوة 1: الخطوة الميكروسكوبية (Microscopic Step)

• استخدام النموذج الميكروسكوبي بـ M متغير
• تقسيم زمني Forward Euler بخطوة δt
• حساب من $w^n$ إلى الحل الوسيط $w^{n,*}$

$w^{n,*} = w^n + \delta t \cdot f(w^n)$

الخطوة 2: خطوة التقييد (Restriction)

• تقييد الحل الميكروسكوبي M-البعد إلى حل ماكروسكوبي L-البعد (M>L)
• الاستفادة من خصائص الأساس المتعامد، القطع المباشر للمعاملات ذات الرتبة العالية
• $W^{n,*} = (h, hu_m, h\alpha_1, ..., h\alpha_{L-2})^{n,*}$

الخطوة 3: الخطوة الماكروسكوبية (Macroscopic Step)

• استخدام النموذج الماكروسكوبي بـ L متغير
• تقسيم زمني Forward Euler بخطوة Δt (عادة Δt>δt)
• حساب من $W^{n,*}$ إلى $W^{n+1}$

$W^{n+1} = W^{n,*} + \Delta t \cdot F(W^{n,*})$

الخطوة 4: خطوة المطابقة (Matching)

• إعادة بناء الحل الميكروسكوبي M-البعد من الحل الماكروسكوبي L-البعد
• الابتكار الرئيسي: الاستفادة من تعامد متعددات Legendre، تقليل مسافة L2

$u^{n+1} = \arg\min_{\bar{u} \in V(W^{n+1})} d(\bar{u}, u^{n,*})$

من خلال التوسع والاستفادة من التعامد: $d(\bar{u}, u^{n,*}) = \int_0^1 (\bar{u}_m - u_m)^2 d\zeta + \sum_{j=1}^{N}(\bar{\alpha}_j - \alpha_j)^2 \int_0^1 \phi_j^2(\zeta) d\zeta$

الحل الأمثل هو: $\bar{\alpha}_i = \alpha_i$ (i=L-1,...,M-2)، أي وراثة مباشرة للمعاملات ذات الرتبة العالية من الحل الميكروسكوبي الوسيط

نقاط الابتكار التقني

1. الاستفادة الماهرة من الأساس المتعامد: تعامد متعددات Legendre يجعل خطوات التقييد والمطابقة بسيطة وفعالة، بدون الحاجة إلى خوارزميات تحسين معقدة
2. ميزة خطوة الزمن:
   • الخطوة الميكروسكوبية: $\delta t = \frac{\Delta x \cdot C}{|\lambda|_{max,M}}$
   • الخطوة الماكروسكوبية: $\Delta t = \frac{\Delta x \cdot C}{|\lambda|_{max,L}}$
   • بما أن $|\lambda|_{max,L} < |\lambda|_{max,M}$، يمكن استخدام خطوة زمنية أكبر في الخطوة الماكروسكوبية
3. تقريب القيم الذاتية: استخدام صيغة تقريبية لتجنب حساب القيم الذاتية المكلف $|\lambda|_{max} = |u_m| + \sqrt{gh + \sum_{j=1}^{N}\alpha_j^2}$ دقيقة لـ N=0,1، محافظة قليلاً لـ N>1 لضمان الاستقرار
4. صيغة التقسيم: تقسيم PDE إلى حدود الحمل والحدود المصدرية وحلها بشكل منفصل، مما يزيد المرونة

إعداد التجارب

مجموعات البيانات (حالات الاختبار)

الحالة 1: اختبار انهيار السد (Dam Break Test)

• الشروط الابتدائية: $w(0,x) = \begin{cases}(2,0,...,0)^T & x \leq 0 \\ (1,0,...,0)^T & x > 0\end{cases}$
• المجال المكاني: -4, 4، 400 خلية شبكية
• المعاملات: λ=ν=0.1، وقت الإنهاء t=2

الحالة 2: اختبار نقل الموجات (Wave Transport Test)

• الشروط الابتدائية: $w(0,x) = (3+e^{-1.5x^2}, 0,...,0)^T$
• المجال المكاني: -4, 4، 400 خلية شبكية
• المعاملات: λ=ν=0.1، وقت الإنهاء t=1

مؤشرات التقييم

1. مؤشرات الدقة:
   • دقة ملف السرعة u(ζ) في مواقع محددة (x=2)
   • التوزيع المكاني لارتفاع الماء h(x)
   • التقارب مع الحل المرجعي (نموذج ميكروسكوبي M=7)
2. مؤشرات الكفاءة:
   • نسبة التسريع (Speed-up) بالنسبة للنموذج الميكروسكوبي المرجعي
   • مقارنة وقت الحساب

طرق المقارنة

1. الحل المرجعي (Micro, M=7): نموذج SWME ميكروسكوبي بحت بـ 7 متغيرات
2. نموذج ماكروسكوبي أساسي (Macro, L=2 أو L=3): نموذج ماكروسكوبي بحت بـ 2-3 متغيرات
3. طريقة الميكروسكوبي-الماكروسكوبي: مجموعات معاملات مختلفة
   • تثبيت L، تغيير M: L=2 أو 3، M=4,5,6,7
   • تثبيت M، تغيير L: M=6 أو 7، L=2,3,4,5

تفاصيل التنفيذ

• التقسيم المكاني: طريقة الحجم المحدود من الدرجة الأولى
• التقسيم الزمني: صيغة Forward Euler الصريحة
• ثابت CFL: C≤1
• لغة البرمجة: Python
• الكود مفتوح المصدر: https://github.com/verbiestrik/PhD-RUG

نتائج التجارب

النتائج الرئيسية

اختبار انهيار السد

تحليل الدقة:

1. تثبيت L وتغيير M (الشكل 3.1-3.2):
   • عند L=2: زيادة M تحسّن ملف السرعة، لكن لا يتقارب مع الحل المرجعي
   • عند L=3: النتائج أقرب بكثير من الحل المرجعي، لكن ملف السرعة لا يزال لم يتقارب تماماً
2. تثبيت M وتغيير L (الشكل 3.3-3.4):
   • عند M=6: مع زيادة L، يتقارب ملف السرعة وارتفاع الماء مع الحل المرجعي
   • عند M=7: النتائج أفضل قليلاً، سلوك التقارب واضح
   • الاكتشاف الرئيسي: عدد المتغيرات الماكروسكوبية L هو العامل الرئيسي المؤثر على الدقة

تأثير التسريع (الجدول 3.1-3.2):

تثبيت M=7، تغيير L:

اختبار نقل الموجات

تحليل الدقة:

1. ارتفاع الماء h(x): يظهر تقارباً جيداً مع زيادة M أو L
2. ملف السرعة u(ζ): تقارب أضعف، خاصة عند تثبيت L وتغيير M
3. الاتجاه العام متسق مع اختبار انهيار السد

تأثير التسريع (الجدول 3.3-3.4):

• متطابق جداً مع اختبار انهيار السد (الفرق <5%)
• اكتشاف مهم: نسبة التسريع لا تعتمد بشكل أساسي على الشروط الابتدائية، لها عمومية
• أفضل تكوين (M=4, L=2): تسريع بمقدار 2.482×

تجارب الاستئصال

التحقق من خطوة المطابقة (الشكل 2.1):

• إعداد الاختبار: M=7، الحل الدقيق هو 1.2 مرة من الحل الوسيط
• مع زيادة L من 2 إلى 6، يتقارب ملف السرعة المطابق تدريجياً مع الحل الدقيق
• حتى عند L=2 (مطابقة مباشرة لـ 2 متغير فقط)، نتيجة المطابقة دقيقة جداً
• الخلاصة: خطوة المطابقة تعمل بشكل جيد، لا تدخل أخطاء كبيرة

تحليل الحالات

الظواهر الفيزيائية في اختبار انهيار السد:

• عدم الاستمرارية الأولية في الارتفاع تؤدي إلى انتشار الموجات
• ملف السرعة عند x=2 يظهر توزيعاً رأسياً غير منتظم
• طريقة الميكروسكوبي-الماكروسكوبي (M=7, L=5) تستطيع التقاط هذه الميزات بدقة

الملاحظات في اختبار نقل الموجات:

• انتشار وتوسع حزمة غاوسية
• ملف السرعة بسيط نسبياً، النموذج الماكروسكوبي يؤدي أداء أفضل
• تقارب ارتفاع الماء أفضل من ملف السرعة

اكتشافات التجارب

1. الدقة تحددها بشكل أساسي النموذج الماكروسكوبي: زيادة L لها تأثير أكبر على تحسين الدقة من زيادة M
2. قانونية نسبة التسريع:
   • كلما كان M أصغر، كان التسريع أوضح (M=4 يحقق أقصى 2.5×)
   • كلما كان L أصغر، كان التسريع أوضح (لكن الدقة تنخفض)
   • التوازن الأمثل: M=6-7, L=3-4، تسريع حوالي 1.5-2×
3. استقلالية المشكلة: نسبة التسريع في حالتي الاختبار المختلفتين متطابقة تقريباً، مما يشير إلى عمومية الطريقة
4. اختلاف سلوك التقارب: ارتفاع الماء يتقارب أسهل من ملف السرعة

الأعمال ذات الصلة

نمذجة تدفقات المياه الضحلة

1. معادلات المياه الضحلة الكلاسيكية (SWE):
   • مشتقة من معادلات نافييه-ستوكس من خلال متوسط العمق
   • تفترض السرعة الرأسية ثابتة
   • فعالة حسابياً لكن دقة محدودة1,5
2. معادلات العزوم للمياه الضحلة (SWME):
   • طورها كوالسكي وتوريليون (2019)1
   • تستخدم توسيع متعدد الحدود لملف السرعة الرأسي
   • تحسّن الدقة بشكل كبير لكن التكلفة الحسابية تزداد
3. مجالات التطبيق:
   • التنبؤ بالطقس2
   • الهيدروليكا3
   • محاكاة الانهيارات الجليدية4

طريقة الميكروسكوبي-الماكروسكوبي

1. التطبيقات في ديناميكا الغازات:
   • كويلرماير وفاندكاستيل (2023) لنماذج عزوم الغازات النادرة6
   • استراتيجية التسريع الهرمي
2. المعادلات التفاضلية العشوائية:
   • ديبرابانت وآخرون (2017) لمحاكاة مونت كارلو7
3. ابتكار هذه الورقة: أول تطبيق للطريقة على تدفقات السطح الحر

الطرق العددية

1. التقسيم الزمني: طريقة Forward Euler الصريحة8
2. التقسيم المكاني: طريقة الحجم المحدود
3. تحليل الاستقرار: شرط CFL والاستقرار المتوازن10

الاستنتاجات والمناقشة

الاستنتاجات الرئيسية

1. جدوى الطريقة: تطبيق ناجح لطريقة الميكروسكوبي-الماكروسكوبي على SWME، خطوات التقييد والمطابقة بسيطة وفعالة بناءً على أساس Legendre المتعامد
2. تحسن الأداء: تحقيق تسريع حسابي بمقدار 2 مرة أو أكثر في حالتي الاختبار مع الحفاظ على دقة جيدة
3. التحكم في الدقة: اختيار النموذج الماكروسكوبي (قيمة L) هو العامل الرئيسي المحدد للدقة
4. العمومية: نسبة التسريع مستقلة عن الشروط الابتدائية، لها قابلية تطبيق واسعة

القيود

1. مشاكل الاستقرار:
   • استخدام شرط CFL فقط، لم يتم النظر في استقرار الحدود المصدرية
   • قد تكون غير مستقرة في الحالات التالية:
     • معامل الاحتكاك ν/λ كبير جداً
     • حجم الشبكة Δx كبير جداً
     • h, um, α1,...,αN قريبة جداً من الصفر
2. قيود التقارب:
   • ملف السرعة لا يتقارب تماماً في بعض التكوينات
   • تقارب أضعف عند تثبيت L وتغيير M
3. قيود الطريقة:
   • اختبار تقسيم مكاني من الدرجة الأولى فقط
   • لم يتم النظر في خطوات زمنية تكيفية
   • افتراضات مبسطة لقاع مستو والجاذبية نحو الأسفل
4. التعقيد الحسابي:
   • التحليل النظري يحذف عوامل ثابتة
   • قد تتأثر نسبة التسريع الفعلية بتفاصيل التنفيذ

الاتجاهات المستقبلية

1. استراتيجيات تبديل النموذج:
   • تقسيم مجال تكيفي: استخدام النموذج الماكروسكوبي في المناطق ذات ملف السرعة البسيط، والنموذج الميكروسكوبي في المناطق المعقدة
   • تكيف زمني: تعديل ديناميكي لـ M و L بناءً على خصائص حقل التدفق
2. تحسينات الاستقرار:
   • النظر في شروط استقرار الحدود المصدرية10
   • صيغ زمنية ضمنية أو شبه ضمنية للتعامل مع الحدود المصدرية الصلبة
3. طرق ذات رتبة أعلى:
   • تقسيم مكاني ذو رتبة عالية (مثل WENO)
   • تكامل زمني ذو رتبة عالية (مثل Runge-Kutta)
4. متغيرات دالة المسافة:
   • استكشاف معايير مطابقة بخلاف معيار L2
   • مطابقة بقيود فيزيائية (مثل حفظ الزخم)
5. سيناريوهات أكثر تعقيداً:
   • تضاريس قاع غير مستوية
   • مشاكل ثنائية الأبعاد
   • نماذج الاضطراب

التقييم المتعمق

المميزات

1. ابتكار الطريقة:
   • أول توسيع لطريقة الميكروسكوبي-الماكروسكوبي إلى تدفقات السطح الحر
   • استخدام ماهر لتعامد Legendre لتحقيق مطابقة فعالة
   • ربط وثيق بين النظرية والممارسة
2. المساهمات النظرية:
   • تحليل شامل للتعقيد الحسابي
   • اشتقاق رياضي واضح
   • حل صريح لخطوة المطابقة
3. كفاية التجارب:
   • اختبار في سيناريوهات فيزيائية مختلفة
   • دراسة نظامية للمعاملات (تغيير M و L)
   • تقييم مزدوج للدقة والكفاءة
   • التحقق المستقل من خطوة المطابقة
4. القيمة العملية:
   • تأثير تسريع كبير (2× أو أكثر)
   • كود مفتوح المصدر يسهل الاستنساخ
   • الطريقة لها عمومية (مستقلة عن المشكلة)
5. جودة الكتابة:
   • هيكل واضح، منطق صارم
   • تعبير رياضي دقيق
   • رسوم بيانية غنية بالمعلومات

أوجه القصور

1. تحليل الاستقرار غير مكتمل:
   • عدم معالجة كافية لاستقرار الحدود المصدرية
   • يحد من نطاق تطبيق الطريقة
   • نقص مناقشة حالات الفشل
2. تحليل التقارب:
   • السبب وراء عدم تقارب ملف السرعة تماماً لم يتم استكشافه بعمق
   • نقص تحليل الخطأ الكمي (مثل خطأ L2)
   • عدم توفير تقديرات نظرية لرتبة التقارب
3. قيود إعداد التجارب:
   • اختبار طريقة من الدرجة الأولى فقط
   • دقة الشبكة ثابتة (400 خلية)
   • عدم مناقشة تمثيلية اختيار المعاملات (λ=ν=0.1)
4. عدم كفاية المقارنة:
   • عدم المقارنة مع طرق تسريع أخرى (مثل الطرق الضمنية، الشبكات المتعددة)
   • نقص مقارنة الدقة مع SWE البحت
   • عدم تقييم الحمل الإضافي للذاكرة
5. عيوب التحليل النظري:
   • التعقيد الحسابي يحذف عوامل ثابتة، القدرة التنبؤية محدودة
   • عدم تحليل انتشار الخطأ
   • خصائص تقريب شرط CFL قد تكون محافظة جداً

التأثير

1. المساهمة الأكاديمية:
   • توفير أداة جديدة لمحاكاة تدفقات المياه الضحلة
   • توسيع نطاق تطبيق طريقة الميكروسكوبي-الماكروسكوبي
   • وضع أساس للأبحاث اللاحقة
2. القيمة العملية:
   • يمكن تطبيقها مباشرة على التنبؤ بالطقس والهندسة المائية وغيرها
   • التسريع بمقدار 2 مرة له أهمية كبيرة في المحاكاة الكبيرة الحجم
   • الكود مفتوح المصدر يعزز التطبيق العملي
3. قابلية الاستنساخ:
   • توفير مستودع GitHub كامل
   • وصف الطريقة مفصل
   • إعدادات المعاملات واضحة
4. القيود:
   • تحتاج إلى هندسة إضافية قبل الاستخدام في بيئة الإنتاج
   • مشاكل الاستقرار تحد من سيناريوهات التطبيق
   • الدقة من الدرجة الأولى قد لا تفي باحتياجات بعض التطبيقات

السيناريوهات المناسبة

موصى به الاستخدام:

1. تدفقات المياه الضحلة بمعامل احتكاك متوسط (ν/λ معتدل)
2. المشاكل التي تتطلب التقاط ملف السرعة الرأسي
3. السيناريوهات ذات الموارد الحسابية المحدودة لكن تحتاج دقة معينة
4. التطوير الأولي والتقييم السريع

غير موصى به الاستخدام:

1. سيناريوهات معامل الاحتكاك العالي (خطر الاستقرار)
2. الشبكات الدقيقة جداً أو المتغيرات قريبة من الصفر
3. التطبيقات التي تتطلب دقة عالية الرتبة
4. الحساب في الوقت الفعلي (التسريع بمقدار 2 مرة قد لا يكون كافياً)

توصيات التكوين الأمثل:

• التوازن بين الدقة والسرعة: M=6-7, L=3-4
• الأولوية للسرعة: M=4-5, L=2
• الأولوية للدقة: M=7, L=5 (لكن التسريع محدود)

المراجع

1 J. Kowalski and M. Torrilhon, "Moment approximations and model cascades for shallow flow," Commun. Comput. Phys., 2019.

6 J. Koellermeier and H. Vandecasteele, "Hierarchical micro-macro acceleration for moment models of kinetic equations," Journal of Computational Physics, 2023.

7 K. Debrabant, G. Samaey, and P. Zielinski, "A micro-macro acceleration method for the monte carlo simulation of stochastic differential equations," SIAM Journal on Numerical Analysis, 2017.

10 J. Koellermeier and Q. Huang, "Equilibrium Stability Analysis of Hyperbolic Shallow Water Moment Equations," Math. Method. Appl. Sci., 2022.

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