High order tracer variance stable transport with low order energy conserving dynamics for the thermal shallow water equations

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

بيان المشكلة
غالباً ما تتطلب النمذجة العددية لديناميكيات السوائل الجيوفيزيائية حلاً وسطاً بين الكفاءة الحسابية والدقة الفيزيائية. تُستخدم المخططات العددية منخفضة الرتبة (الفرق المحدود، الحجم المحدود، والعناصر المحدودة) بشكل متكرر في الحلول الجوية الضمنية لأن أرقام شرطها (condition numbers) المنخفضة تسهل سرعة التقارب وتجميع المصفوفات. ومع ذلك، تعتمد هذه المخططات عادةً على "الرفع العلوي" (upwinding) لكبح الآثار عالية التردد، مما يؤدي إلى تبديد عددي مفرط يقلل من دقة الميزات الاضطرابية. ومن ناحية أخرى، تعد الطرق عالية الرتبة (مثل طريقة غاليركين المتقطعة - Discontinuous Galerine) فعالة في حل الحركات العطالية وبناء التدفقات العددية الدقيقة، لكنها مكلفة حسابياً وغالباً ما تفتقر إلى خصائص الاستقرار القوية التي تتميز بها الحلول الضمنية منخفضة الرتبة. علاوة على ذلك، غالباً ما تواجه مخططات النقل عالية الرتبة القياسية صعوبة في الحفاظ على تباين المتتبع (وهو إنتروبيا عددية للنظام) دون إدخال عدم استقرار، بينما قد تفشل المخططات منخفضة الرتبة في الحفاظ على بنية حفظ الطاقة المطلوبة للاستقرار طويل الأمد في الأنظمة الهاملتونية غير القانونية (non-canonical Hamiltonian systems).

المنهجية
يقترح المؤلفون إطار عمل عددي هجين يربط بين حل العناصر المحدودة المختلطة منخفض الرتبة للديناميكيات، وطريقة "غاليركين المتقطعة" (DG) عالية الرتبة لنقل المتتبعات الديناميكية الحرارية، وتحديداً في سياق معادلات المياه الضحلة الحرارية.

• الديناميكيات منخفضة الرتبة: يستخدم حل الديناميكيات صيغة عناصر محدودة مختلطة (أدنى رتبة متعددة الحدود) تحافظ على البنية المحافظة للطاقة. يعمل هذا الحل على شبكة خشنة ويحسب المتغيرات المستقبلية للسرعة ($u$)، والعمق ($h$)، والطفو الموزون بالعمق ($S$).
• النقل عالي الرتبة: يتم التعامل مع نقل الطفو باستخدام مخطط DG عالي الرتبة يعمل على شبكة متداخلة ومكثفة. يُعامل الطفو كمتغير مادي ($s$) بدلاً من متغير تدفق ($S=sh$). يستخدم المخطط توسعاً متماثلاً مائلاً (skew-symmetric expansion) لمشغل الحمل في المستوى المستمر لضمان أن الصيغة المتقطعة يمكنها حفظ تباين المتتبع بشكل مثبت.
• استراتيجية الربط: يتم ربط الملحقين عبر تسلسل هرمي متداخل من الشبكات. تُدمج شبكة الرتبة المنخفضة داخل شبكة النقل عالية الرتبة بحيث تتطابق مراكز خلايا العناصر منخفضة الرتبة مع نقاط تربيع غاوس-ليجندر (Gauss-Legendre) للعناصر عالية الرتبة. وهذا يخلق تطابقاً بنسبة 1:1 بين درجات الحرية.
  • يتم تشخيص حقل الطفو عالي الرتبة من المتغيرات المستقبلية منخفضة الرتبة في بداية كل خطوة زمنية.
  • يتم دفع معادلة نقل الطفو باستخدام تدفق الكتلة وحقول العمق منخفضة الرتبة (المتوسطة عبر الخطوة الزمنية).
  • يتم إسقاط الطفو عالي الرتبة الناتج مرة أخرى إلى فضاء الرتبة المنخفضة لتحديث حل الديناميكيات.
• خصائص الحفظ:
  • الطاقة: يحافظ الربط على التماثل المضاد لتدفق نقل الطفو وتدرج الضغط، مما يضمن حفظ الطاقة المتقطعة للديناميكيات منخفضة الرتبة.
  • تباين المتتبع: تم تصميم مخطط النقل عالي الرتبة للحفاظ على تباين المتتبع في الفضاء. ولضمان الاستقرار ومنع ضجيج الشبكة من التزييف غير الخطي (nonlinear aliasing)، يتم تضمين حد رفع علوي (upwinding)، مما يسمح بتبديد محكوم لتباين المتتبع (تخميد) بدلاً من النمو غير المحدود.
• التكامل الزمني: يستخدم النظام مكامل (SSP-RK3) من الرتبة الثالثة مستقر للأنظمة الصلبة (stiffly stable) لخطوة النقل، ونهج "شبه نيوتن" (يُحل عبر GMRES) للديناميكيات المقترنة.

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

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

النتائج
تم التحقق من صحة المنهج من خلال ثلاث حالات اختبار قياسية:

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

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

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