Primitive variable regularization to derive novel Hyperbolic Shallow Water Moment Equations

إليك ملخص تقني مفصل لورقة البحث بعنوان "تنظيم المتغيرات الأولية لاشتقاق معادلات عزم المياه الضحلة الهيبيربولية الجديدة" للباحث جوليان كوليرماير.

1. بيان المشكلة

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

لمعالجة ذلك، تم تطوير معادلات عزم المياه الضحلة (SWME). تقوم هذه النماذج بتوسيع ملف تعريف السرعة الرأسية باستخدام كثيرات حدود ليجاندر المتعامدة، مما ينتج عنه نظام من $N+2$ من المعادلات التفاضلية الجزئية (PDEs) لعمق الماء ($h$)، ومتوسط سرعة العمق ($u_m$)، ومعاملات التوسع ($\alpha_i$).

التحدي الجوهري: تعاني نماذج SWME الحالية من قصور تحليلي حرج مزدوج:

1. الافتقار إلى الهيبيربولية العالمية (Global Hyperbolicity): بالنسبة للرتب $N \geq 2$، يكون نظام SWME الأصلي هيبيربوليًا محليًا فقط. فهو يفقد خاصية الهيبيربولية (مما يؤدي إلى عدم استقرار عددي) عند حدوث انحرافات صغيرة عن حالة التوازن أو عند وجود عزم رتب عليا كبير.
2. الافتقار إلى حالات مستقرة تحليلية: النماذج الحالية لا تسمح بوجود حالات مستقرة تحليلية يمكن اشتقاقها من شروط رانكين-هوجونيات (Rankine–Hugoniot)، مما يعيق تطوير مخططات عددية متوازنة جيدًا (well-balanced).

تضمنت المحاولات السابقة لمعالجة هذه المشكلات ما يلي:

• HSWME: تنظيم المتغيرات في الصيغة الحملية (convective variables) عبر ضبط المعاملات من الرتب العليا لتساوي صفرًا. يضمن هذا الهيبيربولية العالمية، ولكنه يفقد الدقة الفيزيائية ولا يزال يفتقر إلى حالات مستقرة تحليلية.
• SWLME: خطية العزوم من الرتب العليا للحصول على حالات مستقرة وهيبيربولية تحليلية، ولكن على حساب الدقة الفيزيائية الكبيرة (فقدان تأثيرات الزخم غير الخطية).

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

1. الحفاظ على توازن الكتلة والزخم.
2. تحقيق الهيبيربولية العالمية.
3. السماح بوجود حالات مستقرة تحليلية.
4. الحفاظ على دقة عالية مقارنة بنموذج SWME غير المنظم.

2. المنهجية

يقترح المؤلفون استراتيجية جديدة لـ "تنظيم المتغيرات الأولية" (Primitive Variable Regularization)، وذلك على خلاف الطرق السابقة التي قامت بتنظيم المتغيرات الحملية.

الخطوات الرئيسية:

1. تحويل المتغيرات: يتم تحويل النظام بين المتغيرات الأولية ($U_p = [h, u_m, \alpha_1, \dots, \alpha_N]$) والمتغيرات الحملية ($U_c = [h, hu_m, h\alpha_1, \dots, h\alpha_N]$) باستخدام مصفوفة تحويل قابلة للعكس $T$.
2. التنظيم في الفضاء الأولي: بدلاً من ضبط المعاملات من الرتب العليا لتساوي صفرًا ($\alpha_i = 0$ لـ $i \geq 2$) مباشرة في مصفوفة النظام الحملي (كما هو متبع في HSWME)، يقوم المؤلفون بالتنظيم في مصفوفة النظام الأولي ($A_p$).
   • يقومون بتقييم مصفوفة النظام الأولي حول ملفات تعريف سرعة خطية (بضبط $\alpha_i = 0$ لـ $i \geq 2$ في معادلات العزم).
3. التحويل العكسي: يتم تحويل مصفوفة النظام الأولي المنظمة مرة أخرى إلى المتغيرات الحملية. ومن الأهمية بمكان أنه بسبب الطبيعة غير الخطية للتحويل، تظهر معاملات الرتب العليا $\alpha_i$ مجددًا في مصفوفة النظام الحملي، مما يحافظ على حدود الاقتران (coupling terms) التي فُقدت في الطرق السابقة.
4. النهج الهجين (PMHSWME): يجمع المؤلفون بين التنظيم الأولي واستراتيجية الإبقاء على معادلة توازن الزخم دون تغيير (وهي ميزة في نماذج SWLME/MHSWME). ينتج عن ذلك نموذج يتم فيه تنظيم معادلات العزم ذات الرتب العليا فقط، بينما تظل معادلات الكتلة والزخم دقيقة تمامًا.

3. المساهمات الرئيسية

أ. اشتقاق نماذج جديدة

تشتق الورقة تسلسلين جديدين من النماذج:

1. PHSWME (معادلات عزم المياه الضحلة الهيبيربولية الأولية): نموذج منظم بالكامل في المتغيرات الأولية.
2. PMHSWME (معادلات عزم المياه الضحلية الهيبيربولية الأولية الهجينة): نموذج هجين ينظم معادلات العزم الأولية ولكنه يُبقي على معادلة الزخم الأصلية دون تغيير.

ب. البراهين التحليلية

• الهيبيربولية العالمية: يثبت المؤلفون أن كلا النموذجين (PHSWWME و PMHSWME) هما هيبيربوليان عالميًا لأي رتبة $N$. وقد أظهرت القيم الذاتية أنها حقيقية ومتمايزة، وهي مشتقة من كثيرات حدود مميزة تتضمن مشتقات كثيرات حدود ليجاندر ($P'_{N+1}$).
• الحالات المستقرة التحليلية: على عكس النماذج السابقة، تسمح هذه الأنظمة الجديدة باشتقاق حالات مستقرة تحليلية. تعتمد هذه الحالات المستقرة على رقم فرود (Froude number) وأرقام العزم غير البعدية، مما يعمم شروط القفزة الهيدروليكية القياسية في معادلات SWE.

ج. الرؤية النظرية

توضح الورقة أن اختيار فضاء المتغيرات للتنظيم أمر بالغ الأهمية. فالتنظيم في المتغيرات الحملية (HSWME/MHSWME) يؤدي إلى هيبيربولية محلية أو حالات مستقرة غير قابلة للتطبيق. أما التنظيم في المتغيرات الأولية فيسمح بإجراء تحويل تشابه يحافظ على البنية الهيبيربولية مع استعادة حدود الاقتران الضرورية في الصيغة الحملية.

4. النتائج

التحقق العددي (حالة اختبار كسر السد)

اختبر المؤلفون النماذج مقابل مشكلة كسر سد قياسية مع ملف تعريف سرعة أولي خطي، وقارنوها بنموذج مرجعي عالي الدقة.

• الدقة:
  • SWLME: كان الأداء الأسوأ، حيث أظهر أخطاءً كبيرة بسبب الخطية المفرطة.
  • HSWME: أدى أداء ضعيف لارتفاع الماء ($h$) ومتوسط السرعة ($u_m$) لأنه غيّر معادلة الزخم.
  • PHSWME: أظهر تحسنًا عن HSWWME ولكنه لا يزال يُظهر أخطاء في $h$ و $u_m$ بسبب التغييرات في اقتران الزخم.
  • PMHSWME: حقق أعلى دقة لجميع المتغيرات ($h, u_m, \alpha_1, \alpha_2$). من خلال الحفاظ على معادلة الزخم الدقيقة مع تنظيم العزوم من الرتب العليا فقط، قلل هذا النموذج من الانحراف عن SWME الأصلي.
• الاستقرار: ظلت النماذج الجديدة (PHSWME و PMHSWME) مستقرة، بينما أصبح نموذج SWME الأصلي غير مستقر للرتب العالية ($N \geq 5$) في حالة الاختبار.

التحليل المقارن

يسلط جدول ملخص (الجدول 1 في الورقة) الضوء على أن PMHSWME هو النموذج الوحيد الذي يلبي جميع الخصائص المطلوبة:

• الهيبيربولية العالمية: نعم
• الحفاظ على معادلة الزخم: نعم
• حالات مستقرة تحليلية: نعم
• معادلات العزم غير الخطية: نعم

5. الأهمية

يعالج هذا العمل مقايضة طويلة الأمد في النمذجة منخفضة الرتبة لتدفقات السطح الحر.

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

في الختام، نجحت الورقة في اشتقاق فئة جديدة من معادلات عزم المياه الضحلة تتغلب على قيود النماذج السابقة من خلال نقل نموذج التنظيم من المتغيرات الحملية إلى المتغيرات الأولية، مما نتج عنه نموذج PMHSWME الذي يتميز بالهيبيربولية العالمية، والقابلية للحل تحليليًا للحالات المستقرة، والدقة العالية.