Interpretation of stochastic primitive equations with relaxed hydrostatic assumption as a higher order approximation of 3D stochastic Navier-Stokes

ملخص تقني: تفسير معادلات العناصر الأولية العشوائية مع فرضية الهيدروستاتيكية المرنة

بيان المشكلة
تبحث هذه الورقة في تقارب الحلول لمعادلات نافير-ستوكس العشوائية ثلاثية الأبعاد (تحديداً ضمن إطار "عدم اليقين الموقعي" أو LU) نحو نظائرها من المعادلات الأولية. تتناول الدراسة أوجه القصور في "فرضية الهيدروستاتيكية القوية" الكلاسيكية، التي تفترض أن التسارع الرأسي مهمل. وبينما تعد هذه الفرضية صالحة لديناميكيات المحيطات واسعة النطاق، إلا أنها تفشل في رصد ظواهر مثل الحمل الحراري العميق والتدفق الصاعد/الهابط القوي. يهدف المؤلفون إلى تحديد ما إذا كان "الفرض الهيدروستاتيكي المرن" (أو الضعيف) — الذي يحتفظ بحدود مارتينجال محددة تمثل الانحرافات عن التوازن الهيدروستاتيكي — يوفر تقريباً أعلى رتبة من المعادلات الأولية لـمعادلات نافير-ستوكس العشوائية ثلاثية الأبعاد مقارنة بالنموذج الهيدروستاتيكي القوي القياسي.

المنهجية
يستخدم المؤلفون إطار "عدم اليقين الموقعي" (LU)، الذي ينمذج تدفق السوائل عبر تفكيك الإزاحة اللاغرانجية إلى سرعة أويلرية واسعة النطاق وحد من الضجيج عالي التذبذب صغير النطاق. تتضمن المنهجية الخطوات التالية:

1. التحجيم والتحليل التقاربي: ينظر المؤلفون في نطاق رقيق $S_\epsilon = S_H \times (-\epsilon, \epsilon)$ حيث $\epsilon$ هو نسبة الجانب. ويقومون باشتقاق نسخة مُحجمة من معادلات نافير-ستوكس (SNS) ذات عدم اليقين الموقعي، ويقدمون معامل تحجيم ضجيج رأسي $\alpha_\sigma$ للتمييز بين مساهمات الضجيج الأفقي والرأسي.
2. اشتقاق النموذج:
   • النموذج الهيدروستاتيكي القوي (PE): مشتق عبر الاستبعاد الصوري لحدود الرتبة $\epsilon^2$ في معادلة الزخم الرأسي.
   • النموذج الهيدروستاتيكي الضعيف ($PE^\epsilon_{\alpha_\sigma}$): مشتق عبر الاحتفاظ بالحدود من الرتبة $\alpha_\sigma \epsilon^2$ و $\alpha_\sigma^2 \epsilon^2$ في معادلة الزخم الرأسي. يتضمن هذا النموذج حدود مارتينجال تمثل انحرافات التسارع الرأسي. ولضمان جودة التعيين (well-posedness)، يقدم المؤلفون تنظيماً عبر نواة تلافيف (convolution kernel) $K$ لتصفية مساهمات الضجيج عالية التردد في معادلة الزخم الرأسي.
3. الأدوات الرياضية:
   • المساقط المعدلة: للتعامل مع حدود الضغط في السياق العشوائي، يعرّف المؤلفون "مساقط ليري المعدلة" ($P$ و $P_\epsilon$). تقوم هذه المساقط بإلغاء حدود التدرج المُحجمة الناشئة عن الضغط العشوائي، مما يسمح بتطبيق أدوات نظرية ليري الكلاسيكية رغم وجود حدود التغاير في حساب إيتو (Itô calculus).
   • أنظمة التقارب: يغطي التحليل نظامين من الشروط الحدية:
     • الغطاء الصلب (Rigid-lid): يُستخدم لدراسة تقارب $L^2-H^1$ للحلول الضعيفة.
     • الدورية الكاملة (Fully Periodic): تُستخدم لدراسة تقارب $H^1-H^2$ للحلول القوية.
4. تقدير الخطأ: يضع المؤلفون تقديرات طاقة للفرق بين حلول معادلات نافير-ستوكس المُحجمة والمنتظمة ($rSNS^\epsilon_{\alpha_\sigma}$) والمعادلات الأولية المُحجمة والمنتظمة ($PE^\epsilon_{\alpha_\sigma}$). ويستخدمون متباينات "جرونوال" العشوائية ومتباينات "بوركهولدر-ديفيدي-باغندي" لتحديد حدود العزوات الثانية للخطأ.

المساهمات والنتائج الرئيسية

1. التقريب من رتبة أعلى: تثبت الورقة أن المعادلات الأولية الهيدروستاتيكية الضعيفة توفر تقريباً من رتبة أعلى لمعادلات نافير-ستوكس (LU) ثلاثية الأبعاد مقارنة بالمعادلات الهيدروستاتيكية القوية، بشرط أن يحقق معامل الضجيج الرأسي $\alpha_\sigma$ شروط تحجيم محددة.
2. معدلات التقارب:
   • بالنسبة لـ النموذج الهيدروستاتيكي الضعيف، يتم حصر الخطأ التقاربي بحد من الرتبة $O(\alpha_\sigma^2 \epsilon^2)$. ويتحقق التقارب في الاحتمال عندما يكون $\alpha_\sigma = o(\epsilon^{-1})$.
   • بالنسبة لـ النموذج الهيدروستاتيكي القوي (أو نافير-ستوكس غير المنتظمة مقارنة بالنموذج الضعيف)، يتضمن حد الخطأ حداً من الرتبة $O(\alpha_\sigma^4 \epsilon^2)$. وبناءً على ذلك، يُظهر النموذج الهيدروستاتيكي الضعيف أنه تقريب أفضل بكثير في الأنظمة التي تكون فيها $\alpha_\sigma$ كبيرة ولكن $\alpha_\sigma \epsilon$ تظل صغيرة.
3. تحديد "المنطقة الرمادية": يحدد المؤلفون نظاماً معيناً حيث $\alpha_\sigma \sim \epsilon^{-\gamma}$ مع $\gamma \in [1/2, 1)$. في هذه المنطقة، تشكل المعادلات الأولية الهيدروستاتيكية الضعيفة تقريباً مناسباً لنظام نافير-ستوكس المنتظم، بينما لا تحقق المعادلات الهيدروستاتيكية القوية ذلك.
4. جودة التعيين ووقت الانفجار: تعمم الدراسة نتائج جودة التعيين الحتمية للنطاقات الرقيقة على السياق العشوائي. وقد ثبت أن زمن انفجار (blow-up time) معادلات نافير-ستوكس المُحجمة يؤول إلى اللانهاية في الاحتمال مع اقتراب نسبة الجانب $\epsilon \to 0$.
5. الابتكار التقني: يُبرز استخدام "مسقط ليري المعدل" كإسهام تقني محوري؛ حيث يسمح بإلغاء حدود تدرج الضغط العشوائية التي تنشأ من "ليما إيتو"، مما يسهل استخدام أدوات التحليل الدالي القياسية في سياق عشوائي قد تؤدي فيه حدود التغاير إلى تعقيد التقديرات.

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