Spectral interpolation in semi-implicit semi-Lagrangian methods for shallow water equations on the sphere

ملخص تقني: الاستيفاء الطيفي في طرق "شبه الضمنية شبه لاغرانج" لمعادلات المياه الضحلة على الكرة

بيان المشكلة
تُعد طرق "شبه الضمنية شبه لاغرانج" (SISL) المعيار المستخدم في نماذج الطقس والمناخ العالمية التشغيلية، لأنها تسمح بخطوات زمنية أكبر مما تسمح به شروط "كورانت-فريدريكس-ليوي" (CFL). ومع ذلك، يوجد تضارب عددي جوهري في المخططات التشغيلية الحالية؛ فبينما تعتمد عملية التقدير العددي لمشتقات المتغيرات التنبؤية مكانياً على طرق دقيقة طيفياً (مثل تحويلات الهارمونيك الكروي [SHT] أو طريقة "دبل فوريه سفير" [DFS])، فإن خطوة الاستيفاء المطلوبة لتقييم المتغيرات عند نقاط المغادرة البعيدة عن الشبكة تعتمد على تقنيات منخفضة الرتبة، مثل كثيرات حدود "لاغرانج" ذات الضرب التنسوري من الرتبة التكعيبية أو الخماسية.

يخلق هذا التباين عدم اتساقاً عددياً. ورغم أن الاستيفاء منخفض الرتبة غالباً ما يُبرر بمخاوف التكلفة الحسابية (بتعقيد قدره $O(n^2)$ مقابل $O(n^4)$ للاستيفاء الطيفي)، إلا أن هناك أدلة عددية محدودة في الأدبيات تبرر هذه المقايضة للأنظمة المعقدة مثل معادلات المياه الضحلة (SWE). علاوة على ذلك، يُعرف عن الاستيفاء منخفض الرتبة أنه يتسبب في انتشار عددي مفرط وأخطاء في الطور، مما قد يؤدي إلى تدهور الدقة وخصائص الحفظ عبر فترات التكامل الطويلة.

المنهجية
يقدم هذا العمل أول تحقيق عددي للاستيفاء الدقيق طيفياً ضمن مخططات SISL لمعادلات المياه الضحلة على الكرة. يبني المؤلفون عملهم على نموذج SISL القائم على طريقة DFS الذي طوره يوشيمورا [56]، والذي يستخدم أساس فوريه مستقراً ومخطط SETTLS ذو مستويين زمنيين.

الابتكار المنهجي الجوهري هو استبدال استيفاء "لاغرانج" القياسي بمخطط استيفاء دقيق طيفياً يعتمد على نفس أساس فوريه (DFS) المستخدم في المشتقات المكانية. ولمعالجة حاجز التكلفة الحسابية، استخدم المؤلفون "تحويل فوريه السريع غير المنتظم" (NUFFT).

• آلية الاستيفاء: يتم توسيع المتغيرات التنبؤية (تحديداً المصطلحات $R_v$ و $R_h$ في النظام شبه الضمني) باستخدام أساس DFS ذي الانتظام الجزئي. ثم يتم تقييم هذه التوسعات عند نقاط المغادرة غير الشبكية.
• تقليل التعقيد: إن التقييم المباشر لسلاسل فوريه ثنائية المتغير عند $O(n^2)$ من نقاط المغادرة سيتطلب عمليات بقدر $O(n^4)$. ومن خلال إعادة صياغة السلسلة في شكل أسي مركب واستخدام تقنية NUFFT (تحديداً مكتبة FINUFFT)، تم تقليل التعقيد الحسابي إلى $O(n^2 \log^2(1/\epsilon) + n^2 \log n)$، حيث $\epsilon$ هو حد دقة الخطأ. وهذا يجعل خطوة الاستيفاء متسقة حسابياً مع بقية مخطط SISL الطيفي.
• التنفيذ: قام المؤلفون بتنفيذ هذا الاستيفاء الطيفي في سياقين:
  1. نموذج SISL قائم على DFS (باستخدام أساس فوريه المستقر ليوشيمورا).
  2. نموذج SISL قائم على SHT (باستخدام شبكات غاوس-ليجندر)، حيث تم تكييف طريقة DFS لتوليد مستنتجات فوريه ثنائية المتغير من خطوط العرض غير المتساوية المسافات.

المساهمات الرئيسية والنتائج
قام المؤلفون بتقييم مخططات الاستيفاء الطيفي الجديدة مقابل استيفاء "لاغرانج" (LAG) منخفض الرتبة باستخدام مجموعة من اختبارات معادلات المياه الضحلة القياسية (Williamson TC1, TC2, TC5, TC6، واختبار Galewsky). وتظهر النتوات عدة مزايا متميزة:

1. الدقة والتقارب:

   • TC1 (نقل جرس جيب التمام): حافظ الاستيفاء الطيفي لـ DFS على نفس معدل التقارب المكاني لـ LAG، ولكنه حقق أخطاءً أقل بنحو رتبتين من المقدار. وعلى مدى فترات تكامل طويلة (10 دورات)، قلل منهج DFS الأخطاء بمقدار ثلاث رتب من المقدار مقارنة بـ LAG.
   • TC2 (الجريان الجيوستروفي): بينما فشل استيفاء LAG في إظهار التقارب الزمني المتوقع من الدرجة الثانية (ربما بسبب سيطرة الأخطاء المكانية)، حققت مخططات استيفاء DFS دقة الدرجة الثانية النظرية.
   • TC5 و TC6 (تدفق الجبال وموجات روسبي-هورنز): أنتج منهج DFS باستمرار أخطاءً أقل في حقول الارتفاع مقارنة بـ LAG، خاصة عند الدقة العالية والخطوات الزمنية الصغيرة. وفي TC6، حافظ منهج DFS بشكل أفضل على شكل قيعان الموجة، مما جعله يطابق الحلول المرجعية عالية الرتبة لـ Euler بشكل أوثق.

2. خصائص الحفظ:

   • حسن الاستيفاء الطيفي بشكل كبير من حفظ الكتلة والطاقة. في TC1، كان خطأ الكتلة المطلق بعد 120 يوماً هو $1.21 \times 10^{-11}$ لمنهج DFS مقابل $4.34 \times 10^{-4}$ لمنهج LAG.
   • في TC2 و TC5، حافظ منهج DFS على استقرار أخطاء حفظ الكتلة والطاقة بمرور الوقت، بينما كانت أخطاء LAG تميل إلى الزيادة أو التذبذب بشكل أكبر.

3. الانتشار العددي والديناميكيات صغيرة النطاق:

   • اختبار Galewsky: يتضمن هذا الاختبار تطور نفاث غير خطي إلى اضطراب. كشف التحليل الطيفي لطيف الطاقة الحركية أن استيفاء LAG تسبب في تبديد سريع للطاقة عند الأعداد الموجية العالية. في المقابل، حافظت مخططات استيفاء DFS على تسلسل الطاقة، مما سمح لطيف الأعداد الموجية العالية بالاقتراب من مقياس $k^{-5/3}$ النظري الملاحظ في التدفقات المضطربة الدوارة. فشل استيفاء LAG في التقاط هذا السلوك، مما يشير إلى وجود تخميد عددي مفرط.

4. الكفاءة الحسابية:

   • أظهرت مقارنات وقت التنفيذ الفعلي (Wall-clock time) أنه بالنسبة لمتطلبات الدقة المنخفضة، يكون استيفاء LAG أسرع. ومع ذلك، كلما زادت الدقة المطلوبة (مع تشديد حدود الخطأ)، يصبح منهج DFS أكثر كفاءة، حيث يصل إلى الدقة المستهدفة في وقت أقل من منهج LAG. وهذا يؤكد أن الاستيفاء الطيفي المعزز بـ NUFFT هو خيار قابل للتطبيق للتطبيقات عالية الدقة.

الأهمية
يزعم البحث أن دمج الاستيفاء الدقيق طيفياً في مخططات SISL يحل عدم التطابق العددي بين التقدير المكاني واستيفاء الإزاحة (Advection). وتثبت الدراسة أن هذا النهج ليس مجرد تحسين نظري، بل يقدم فوائد ملموسة:

• يقضي على "فجوة الدقة" بين حسابات المشتقات الطيفية وخطوات الاستيفاء.
• يقلل من الانتشار العددي، مما يحافظ على الهياكل دقيقة النطاق والسمات الموجية عبر فترات التكامل الطويلة.
• يعزز حفظ الكتلة والطاقة، وهو أمر بالغ الأهمية لمحاكاة المناخ طويلة الأمد.
• يحقق هذه التحسينات دون تكبد تكاليف حسابية باهظة، بفضل تسريع NUFFT.

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