Structure-Preserving High-Order Methods for the Compressible Euler Equations in Potential Temperature Formulation for Atmospheric Flows
تُطوّر هذه الورقة وتتحقق من صحة طرق "غاليركن" المتقطعة من الرتب العالية والمحافظة على البنية لمعادلات "أويلر" القابلة للانضغاط بصيغة درجة الحرارة الكامنة، والتي تتميز بتدفقات عددية وتجزئات لـحدود المصدر تضمن حفظ الإنتروبيا، والطاقة الكلية، والطاقة الحركية على شبكات منحنية عامة مع الحفاظ على خصائص التوازن الجيد للتدفقات الجوية.
المؤلفون الأصليون:Marco Artiano, Oswald Knoth, Peter Spichtinger, Hendrik Ranocha
تخيل أنك تحاول التنبؤ بالطقس من خلال مراقبة محيط هائل وغير مرئي من الهواء يلتف حول الأرض. ليس هذا مجرد نسيم لطيف؛ بل هو رقصة فوضوية وعالية السرعة من الغاز الذي يخضع لقواعد فيزيائية صارمة. يستخدم العلماء حواسيب قوية لحل هذه القواعد، المعروفة باسم معادلات "أويلر"، للتنبؤ بالعواصف، والتيارات النفاثة، وأنماط المناخ. ومع ذلك، فإن هذه المعادلات معقدة؛ فهي تتضمن توازناً دقيقاً بين سرعة حركة الهواء (الطاقة الحركية)، ودرجة حرارته (درجة الحرارة)، وكثافته (الكثافة). إذا اختل توازن المحاكاة الحاسوبية ولو قليلاً — مثل راقص يفقد إيقاعه — فقد يؤدي ذلك إلى خلق عواصف وهمية أو جعل الهواء يختفي، مما يؤدي إلى إفساد التوقعات. ولإصلاح ذلك، طور الباحثون طرقاً خاصة "محافظة على البنية". فكر في هذه الطرق كأنها مجموعة من القواعد الذهبية التي تجبر الحاسوب دائماً على الحفاظ على إجمالي الطاقة و"النظام" في الهواء (الإنتروبيا) بشكل دقيق تماماً، بغض النظر عن مدى جموح المحاكاة. يتعمق هذا البحث في طريقة شائعة لوصف هذا الهواء باستخدام مفهوم يُسمى "درجة الحرارة الكامنة"، وهي بمثابة رمز سري يستخدمه خبراء الأرصاد الجوية لتتبع كيفية تسخين أو تبريد كتل الهواء أثناء صعودها وهبوطها.
لقد قام مؤلفو هذه الورقة، ماركو أرتيانو وفريقه، ببناء مجموعة جديدة من الأدوات الرياضية لجعل محاكاة الطقس هذه أكثر دقة واستقراراً، وتحديداً لمنهج "درجة الحرارة الكامنة". لقد صمموا ثلاثة "تدفقات عددية" خاصة — وهي بمثابة القواعد التي يستخدمها الحاسوب لتبادل المعلومات بين الخلايا الشبكية المتجاورة. تم هندسة هذه القواعد لضمان الحفاظ على إجمالي الطاقة والإنتروبيا الديناميكية الحرارية (مقياس الفوضى) بشكل مثالي، حتى على الخرائط المعقدة والمنحنية للأرض. وكانت إحدى العقبات الرئيسية التي واجهوها هي الجاذبية. في العالم الحقيقي، تسحب الجاذبية الهواء للأسفل، مما يخلق توازناً دقيقاً يُعرف بالتوازن الهيدروستاتيكي. وإذا لم يستطع النموذج الحاسوبي محاكاة هذا التوازن بدقة، فإنه سيخلق رياحاً وهمية تطغى على الطقس الحقيقي. وقد أثبت الفريق أنه من خلال التعامل مع حد الجاذبية بطريقة معينة وذكية (كـ "حاصل ضرب غير محافظ")، يمكن لطرقهم الجديدة الحفاظ على هذا التوازن على الشبكات المنحنية، سواء كان الهواء في حالة درجة حرارة ثابتة أو حالة درجة حرارة كامنة ثابتة.
ولم يتوقفوا عند الجانب النظري فحسب؛ بل اختبروا أفكارهم من خلال سلسلة من المحاكاة الصعبة. فقد أجروا طرقهم الجديدة مقابل الطرق القديمة باستخدام سيناريوهات جوية كلاسيكية، مثل الأمواج التي تضرب الجبال وتكون عدم الاستقرار الباروكليني (المحركات وراء عواصف خطوط العرض المتوسطة). وفي هذه الاختبارات، كان لأداء منهج "درجة الحرارة الكامنة" الجديد أداء يضاهي طريقة "الطاقة الكلية" التقليدية، ولكن مع الميزة الإضافية المتمثلة في سهولة تنفيذها في نماذج الطقس الحالية. على سبيل المثال، في اختبار يحاكي "موجة كثافة" حيث يكون الضغط والسرعة ثابتين، حافظت تدفقاتهم الجديدة على حفظ الطاقة والإنتروبيا بشكل مثالي، بينما عانت الطرق الأقدم ما لم تُطبق تعديلات رياضية محددة. وفي اختبارات أكثر تعقيداً، مثل "دوامة تايلور-غرين" (وهي معيار لتدفق السوائل الدوامية)، حافظت طرقهم على الحفظ لفترات طويلة دون أن تنفجر المحاكاة.
ومن الأهمية بمكان أن الورقة تسلط الضوء على وجود مقايضة. فمن أجل تحقيق حفظ مثالي لكل من الطاقة والإنتروبيا مع الحفاظ على استقرار ضغط الهواء، تتطلب الطرق الجديدة نوعاً معيناً من المتوسط الرياضي (باستخدام "المتوسط اللوغاريتمي") لكيفية حساب الكثافة. فإذا استخدموا متوسطاً أبسط، فقد يصبح الضغط سالباً نظرياً، وهو أمر مستحيل فيزيائياً. ويظهر المؤلفون أنه من خلال الاختيار الدقيق لهذه المتوسطات، يمكنهم الحفاظ على "توازن الضغط"، مما يعني أن الهواء لن يبدأ في التحرك تلقائياً عندما ينبغي أن يكون ساكناً. كما أظهروا أن طرقهم تعمل على "الشبكات المنحنية" — وهي الشبكات التي يتم تشويهها لتتبع انحناء الأرض أو أشكال الجبال — دون فقدان استقرارها. وفي عمليات محاكاة لتدفق الهواء فوق الجبال، أنتج نهجهم "غير المحافظ" للجاذبية نتائج سلسة وواقعية، بينما خلقت الطرق القديمة "النقطية" تموجات مزعجة ووهمية في البيانات.
في النهاية، يشير هذا العمل إلى أن صياغة درجة الحرارة الكامنة هي بديل قوي وقابل للتطبيق لنمذجة الغلاف الجوي عالية الدقة. ومن خلال إثبات أن هذه التدفقات الجديدة يمكنها الحفاظ على الطاقة، والإنتروبيا، والتوازن الهيدروستاتيكي على الشبكات المنحنية، يوفر المؤلفون أساساً متيناً لبناء نماذج طقس ومناخ أفضل. وتظهر عمليات المحاكاة التي أجروها أن هذه الطرق ليست سليمة من الناحية النظرية فحسب، بل هي فعالة عملياً أيضاً، مما يوفر طريقة لتشغيل محاكاة جوية طويلة الأمد وعالية الدقة تظل مستقرة ودقيقة دون الحاجة إلى "ضبط" إضافي أو إصلاحات اصطناعية. وبينما تركز الورقة على البناء الرياضي والاختبار العددي لهذه المخططات، فإن النتائج تشير إلى مسار واعد نحو جعل توقعاتنا الرقمية للطقس أكثر موثوقية واتساقاً من الناحية الفيزيائية.
ملخص تقني: طرق عالية الرتبة حافظة للبنية لمعادلات أويلر القابلة للانضغاط بصيغة درجة الحرارة الكامنة
بيان المشكلة تتناول الورقة البحثية تطوير طرق عددية حافظة للبنية لمعادلات أويلر القابلة للانضغاط، والمصممة خصيصاً لتدفقات الغلاف الجوي. وبينما تُستخدم طرق "ديسكونتنيوس غاليركين" (Discontinuous Galerkin - DG) على نطاق واسع نظراً لدقتها العالية ورتبتها المرتفعة وقابليتها للتوسع في الحوسبة المتوازية، إلا أنها غالباً ما تتطلب تثبيتاً إضافياً (مثل المرشحات أو اللزوجة الاصطناعية) للحفاظ على الاستقرار في المحاكاة طويلة الأمد، مما قد يؤدي إلى تدهور الدقة. ويتمثل تحدٍ محدد في نمذجة الغلاف الجوي في اختيار المتغيرات التنبؤية؛ فبينما تُعد صيغة الطاقة الكلية معياراً في ديناميكا السوائل الحسابية (CFD) العامة، تُفضل صيغة درجة الحرارة الكامنة (θ) في نماذج الغلاف الجوي نظراً لدورها النظري في الأرصاد الجوية وقدرتها على اختزال المعادلة الديناميكية الحرارية إلى معادلة نقل بسيطة. ومع ذلك، فإن بناء مخططات حافظة للبنية (تحفظ الإنتروبيا، والطاقة الكلية، والطاقة الحركية) باستخدام صيغة درجة الحرارة الكامنة، لا سيما على الشبكات المنحنية العامة مع حدود المصدر الجاذبية، يظل مهمة معقدة. علاوة على ذلك، تعاني المخططات الموجودة "جيدة التوازن" (well-balanced) من صعوبة في الحفاظ على التوازن الهيدروستاتيكي لكل من حالات الخلفية متساوية الحرارة (isothermal) وحالات درجة الحرارة الكامنة الثابتة على الشبكات المنحنية.
المنهجية يطور المؤلفون إطار عمل يعتمد على طريقة "ديسكونتنيوس غاليركين الطيفية العناصر" (DGSEM) باستخدام نهج فرق التدفق (flux-differencing). ويتضمن جوهر المنهجية ما يلي:
الصياغة: تركز الدراسة على مجموعتين من المعادلات: صيغة الطاقة الكلية التقليدية وصيغة درجة الحرارة الكامنة. وفي الأخيرة، المتغيرات التنبؤية هي الكثافة (ϱ)، والزخم (ϱV)، ودرجة الحرارة الكامنة (ϱθ).
حفظ الإنتروبيا والطاقة: اشتق المؤلفون تدفقات عددية جديدة ثنائية النقاط مصممة لتلبية خصائص حفظ محددة:
حفظ الإنتروبيا (EC): ضمان حفظ الإنتروبيا الديناميكية الحرارية.
حفظ الطاقة الكلية (TEC): ضمان حفظ الطاقة الكلية (الطاقة الداخلية + الحركية).
حفظ الطاقة الحركية (KEP): الحفاظ على الطاقة الحركية في غياب تدرجات الضغط.
حفظ توازن الضغط (PEP): الحفاظ على حالات الضغط والسرعة الثابتة. وتوضح الورقة أنه بالنسبة لصيغة درجة الحرارة الكامنة، فإن تحقيق EC وTEC وKEP وPEP في آن واحد يتطلب أن يعتمد تدفق الكثافة على الضغط، وهو مقايضة ضرورية تم تحديدها من خلال تحليل شروط القفزة ومتغيرات الإنتروبيا.
التعامل مع حد المصدر: للتعامل مع حد الجهد الجاذبي (ϕ)، يعالج المؤلفون حد المصدر كحاصل ضرب غير محافظ. وقد اعتمدوا نهج فرق التدفق حيث يتم تمثيل حد المصدر عددياً باستخدام قيم متوسطة محددة (المتوسط اللوغاريتمي للحالات متساوية الحرارة، ومتوسط ستولارسكي لحالات درجة الحرارة الكامنة الثابتة) لضمان أن يكون المخطط جيد التوازن.
التمديد للمنحنيات: تم تمديد الطرق إلى شبكات منحنية ثلاثية الأبعاد. كما أثبت المؤلفون أن التدفقات المشتقة وتمثيلات حد المصدر تحافظ على خصائصها الحافظة للبنية (TEC وKPEP والتوازن) تحت أي خرائط إحداثيات، بشرب توفر هويات المقياس المنفصلة.
التنفيذ: تم تنفيذ المخططات باستخدام لغة جوليا (Julia) ضمن إطار العمل Trixi.jl، مع استخدام التكامل الزمني SSPRK43.
المساهمات الرئيسية
تدفقات جديدة لدرجة الحرارة الكامنة: اشتقت الورقة ووصفَت تدفقات عددية جديدة (TEC وEC وتدفق ETEC مدمج) مخصصة لصيغة درجة الحرارة الكامنة. وأثبتت أن تدفق ETEC هو التدفق الفريد الذي يحقق EC وTEC وKEP وPEP عند السرعة الثابتة، رغم أنه يستلزم تدفق كثافة يعتمد على الضغط.
تعميم KPEP وTEC مع الجاذبية: عمم المؤلفون خاصية حفظ الطاقة الحركية والجهد (KPEP) ووضعوا شروطاً لحفظ الطاقة الكلية في وجود حد الجهد الجاذبي لكل من صيغتي الطاقة الكلية ودرجة الحرارة الكامنة.
مخططات عامة جيدة التوازن: تتمثل إحدى المساهمات الهامة في تطوير مخططات جيدة التوازن قادرة على حفظ التوازن الهيدروستاتيكي لكل من حالات الخلفية متساوية الحرارة (T= ثابت) ودرجة الحرارة الكامنة الثابتة (θ= ثابت) على الشبكات المنحنية. وغالباً ما كانت الأعمال السابقة تقتصر هذه القدرة على الشبكات الكارتيزية أو حالات خلفية محددة.
معالجة حد المصدر غير المحافظ: طبقت الورقة بصرامة تمثيل حاصل الضرب غير المحافظ لتمثيل حد المصدر الجاذبي ضمن إطار DGSEM، موضحين كيفية إعادة كتابة المعادلات شبه المنفصلة إلى صيغة فرق التدفق التي تحفظ الخصائص الهيكلية المطلوبة.
النتائج تحقق المؤلفون من النتائج النظرية من خلال سلسلة من التجارب العددية:
خصائص الحفظ: في غياب الجاذبية، تُظهر التدفقات الجديدة حفظاً دقيقاً للإنتروبيا والطاقة الكلية حتى دقة الآلة للتدفقات السلسة (موجة الكثافة) والتدفقات المضطربة (دوامة تايلور-غرين).
جودة التوازن (Well-Balancedness): على الشبكات المنحنية المشوهة، تحافظ المخططات على التوازن الهيدروستاتيكي لكل من حالات الخلفية متساوية الحرارة (isothermal) والآديباتية (ثابتة θ)، حيث تظل نسبة الخطأ عند دقة الآلة لفترات زمنية طويلة (500,000 خطوة زمنية)، مما يؤكد قوة تمثيل حد المصدر غير المحافظ.
التقارب: يُظهر تحليل التقارب لموجات القصور الذاتي والجاذبية أن التمثيل غير المحافظ لحد المصدر يعطي أخطاءً أقل بكثير (تصل إلى أربعة أوامر مقدارية في السرعة الرأسية) مقارنة بالتمثيل النقطي، خاصة في الشبكات المنحنية. وتحقق المخططات معدلات التقارب العالية المتوقعة.
سيناريوهات الغلاف الجوي: تم اختبار الطرق في اختبارات مرجعية كلاسيكية لموجات الجبل (الموجات الهيدروستاتيكية الخطية، وغير الهيدروستاتيكية الخطية، وجبل شآر). وتُظهر النتائج أن الصيغ غير المحافظة أكثر قوة، حيث تظهر تذبذبات زائفة أقل مقارنة بتمثيلات حد المصدر النقطية. ويؤدي أداء صيغة درجة الحرارة الكامنة أداءً مقارباً لصيغة الطاقة الكلية في هذه السيناريوهات.
عدم الاستقرار الباروكليني: بدأت الورقة باختبارات على معيار عدم الاستقرار الباروكليني (DCMIP 2016) على شبكة مكعبة (cubed-sphere)، مما يعزز التحقق من قابلية تطبيق هذه الطرق على تدفقات الغلاف الجوي العالمية المعقدة.
الأهمية تزعم الورقة أن هذه التطورات تحظى باهتمام كبير من مجتمع نمذجة الغلاف الجوي. فمن خلال توفير طرق حافظة للبنية مخصصة لصيغة درجة الحرارة الكامنة، يسهل هذا العمل دمج مخططات عددية عالية الرتبة، ومستقرة، ودقيقة في نماذج الغلاف الجوي الحالية دون الحاجة إلى الانتقال إلى متغيرات الطاقة الكلية. إن القدرة على حفظ التوازن الهيدروستاتيكي على الشبكات المنحبية العامة لحالات خلفية متعددة تُعد تقدماً حاسماً للمحاكاة الجوية الواقعية، حيث تعد التضاريس المعقدة والتكامل طويل الأمد أمراً معتاداً. يجسد هذا العمل الفجوة بين المزايا النظرية لصيغة درجة الحرارة الكامنة والمتطلبات العملية للتمثيل العددي الحافظ للبنية.