An IMEX-DG solver with non-conforming mesh refinement for atmospheric dynamics with rotation

ملخص تقني: حلّ IMEX-DG مع تحسين الشبكة غير المتوافقة للديناميكا الجوية مع الدوران

بيان المشكلة
تتميز التدفقات الجوية بالتفاعل بين مقاييس مكانية وزمنية متعددة، تتراوح من الدورات المتوازنة واسعة النطاق إلى الموجات الصوتية والجاذبية السريعة. في نماذج الغلاف الجوي القابلة للانضغاط بالكامل، يؤدي التعايش بين الأنماط الجيوفيزيائية البطيئة والأنماط الصوتية السريعة إلى صلابة شديدة (stiffness)، مما يفرض قيوداً صارمة على حجم الخطوة الزمنية في طرق التكامل الصريحة بالكامل. وبينما تم دمج التكامل الزمني الضمني-الصريح (IMEX) وتجزئة الفضاء باستخدام طريقة "ديسكونتنيوس غاليركن" (DG) عالية الرتبة بنجاح لمعالجة الصلابة الصوتية في التكوينات الكارتيزية غير الدوارة، فإن توسيع هذه الأطر لتشمل الدوران الكوكبي يعد أمراً ضرورياً للتنبؤ بالطقس العددي (NWP) والمعايير المرجعية القياسية. ويضيف إدراج قوة كوريوليس تعقيداً إضافياً، لا سيما فيما يتعلق بالحفاظ على التوازن الجيوفستيري والتنفيذ الفعال للمصطلحات الدورانية ضمن الحلول عالية الرتبة التي لا تعتمد على المصفوفات (matrix-free).

المنهجية
يقدم المؤلفون توسيعاً لمحلل IMEX-DG موجود بالفعل عالي الرتبة، تم تنفيذه ضمن مكتبة العناصر المحدودة deal.II للتعامل مع معادلات أويلر القابلة للانضغاط مع التأثيرات الدورانية.

• المعادلات الحاكمة: يعالج الحل معادلات أويلر القابلة للانضغاط وغير اللزجة على نطاق محدود، مع دمج تسارع كوريوليس ($-f \mathbf{k} \times \mathbf{u}$) وتسارع الجاذبية. يتم إغلاق النظام باستخدام معادلة حالة الغاز المثالي المتوافق حرارياً.
• التجزئة الزمنية: تم استخدام مخطط (IMEX) من نوع "أدديتيف رانج-كورتا" (Additive Runge–Kutta) من الدرجة الثانية (تحديداً طريقة TR-BDF2). تُعامل المساهمات الصلبة (تدرجات الضغط، الجاذبية، والدوران) ضمنياً، بينما تُعامل مصطلحات الحمل الحراري غير الصلبة صراحةً. وهذا يسمح بتحديد الخطوة الزمنية بناءً على مقاييس التدفق ذات الصلة بالأرصاد الجوية بدلاً من سرعة الصوت.
• التجزئة المكانية: تُستخدم طريقة DG عالية الرتبة على شبكات تكيفية غير متوافقة (رباعية الأضلاع في 2D، وسداسية الأوجه في 3D). يستخدم الصيغة المعتمدة على قاعدة عقدية (nodal basis) على عناصر حاصل الضرب التنسوري لتمكين التقييم الفعال للمؤثرات التي لا تعتمد على المصفوفات.
• معالجة المصطلح الدوراني: تم تقديم ومقارنة استراتيجيتين بديلتين للتعامل مع المصطلح الدوراني ضمن تكرار "بيكارد" غير الخطي (nonlinear fixed-point iteration):
  1. الاستراتيجية R1: يتم تأخير (lagged) المصطلحات الدورانية والجاذبية إلى الطرف الأيمن، حيث تُقيم باستخدام تكرار بيكارد السابق.
  2. الاستراتيجية R2: تُعامل المصطلحات الدورانية والجاذبية ضمنياً في الطرف الأيسف من نظام الزخم-الطاقة المقترن.
• التنفيذ الذي لا يعتمد على المصفوفات (Matrix-Free): تتمثل المساهمة التقنية الرئيسية في اشتقاق صيغة لمعكوس المؤثر $(A + R)$، حيث $A$ هي مصفوفة الكتلة و $R$ تمثل مصطلح كوريوليس. من خلال استغلال البنية الجبرية لمؤثر كوريوليس (تحديداً أن $R = \beta A J$، حيث $J$ هو دوران مستوٍ)، اشتق المؤلفون تحليلاً (factorization) هو $(A+R)^{-1} = (I + \beta J)^{-1}A^{-1}$. وهذا يسمح بتطبيق المعكوس عبر عمليات معيارية لا تعتمد على المصفوفات متبوعة بتحويل جبري محلي، مما يتجنب الحاجة إلى حل أنظمة خطية مع المؤثر المقترن الكامل.

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

1. التوسيع ليشمل التدفقات الدوارة: التكامل الناجح لقوة كوريوليس في إطار عمل IMEX-DG عالي الرتبة وقابل للانضغاط ضمن مكتبة deal.II.
2. محلل دوراني فعال لا يعتمد على المصفوفات: تطوير تحليل جبري محدد لمؤثر كوريوليس يحافظ على كفاءة التنفيذات التي لا تعتمد على المصفوفات، متجنباً التكلفة الحسابية لحل الأنظمة الخطية الكبيرة للمصطلح الدوراني.
3. مقارنة الاستراتيجيات التكرارية: مقارنة منهجية بين معالجتي "بيكارد" غير الخطيتين للمصطلحات الدورانية (R1 و R2)، وتقييمهما من حيث الدقة، والمتانة، والكفاءة الحسابية.
4. التحقق على الشبكات التكيفية: أول تطبيق لإطار عمل IMEX-DG الدوار على شبكات غير متوافقة تكيفية، مما يثبت القدرة على إجراء تحسين محلي بالقرب من التضاريس مع الحفاظ على الاستقرار والدقة.

النتائج العددية
تم التحقق من الصيغة المقترحة من خلال سلسلة من الاختبارات المرجعية ثنائية وثلاثية الأبعاد:

• موجات الجاذبية-القصور الذاتي في 2D: أظهرت دراسات التقارب في اختبار مرجعي لموجات الجاذبية-القصور الذاتي أن كلتا الاستراتيجيتين R1 و R2 تحققان معدلات تقارب مثالية تقريباً تتوافق مع الرتبة النظرية للمخطط. وأكدت النتائج أن الحل يلتقط انتشار الموجات والحفاظ على الديناميكا المتوازنة بدقة.
• الموجات واسعة النطاق: أظهر اختبار واسع النطاق (بامتداد أفقي 60,000 كم) متانة الحل عبر فترات زمنية طويلة (96 ساعة). ظل خطأ التوازن الجيوسفيري المنفصل محدوداً وأقل من $10^{-7}$، مما يشير إلى أن عدم وجود تجزئة متوازنة تماماً لم يؤدِ إلى تذبذبات زائفة أو عدم استقرار.
• التدفق ثلاثي الأبعاد فوق التضاريس:
  • تلة متوسطة الانحدار: أظهرت عمليات المحاكاة على نطاق 60×40×16 كم أن إدراج الدوران يكسر بشكل صحيح التماثل المحوري لأنماط الموجات والاضطرابات، بما يتوافق مع التوقعات الفيزيائية. كانت النتائج متسقة بين الشبكات المنتظمة وغير المتوافقة، حيث وفرت الأخيرة دقة أفضل للتدرجات بالقرب من الجبل بتكلفة حسابية أقل.
  • تلة واسعة النطاق: سلط اختبار جديد أكبر (240×240×30 كم) الضوء على عدم التماثل الواضح الناتج عن الدوران في حقول السرعة، والذي تم التقاطه بوضوح دون تذبذبات زائفة.
• القابلية للتوسع المتوازي: أظهرت اختبارات التوسع القوي (Strong-scaling) على ما يصل إلى 2048 نواة (باستخدام سوبركمبيوتر MeluXina) أن إدراج المصطلحات الدورانية لا يقلل بشكل كبير من الكفاءة المتوازية للمحلل الأصلي غير الدوار. حافظ الحل على سلوك قريب من المثالية في التوسع، وأظهرت تهيئة الشبكة غير المتوافقة قابلية توسع مماثلة لحالة الشبكة المنتظمة.

الأهمية والادعاءات
يدعي البحث أن إطار عمل IMEX-DG الدوار المقترح يوسع بنجاح قدرات حلول الغلاف الجوي القابلة للانضغاض الحالية لتشمل الدوران الكوكبي دون المساس بدقة أو متانة أو قابلية توسع الطريقة الأصلية. ويؤكد المؤلفون أن المعالجة التي لا تعتمد على المصفوفات لمؤثر الدوران هي مُمكّن حاسم للكفاءة، مما يسمح للمحلل بالتعامل مع تعقيد قوة كوريوليس بأقل قدر من العبء الإضافي.

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