Gauss Principle in Incompressible Flow: Unified Variational Perspective on Pressure and Projection

إليك ملخص تقني مفصل للورقة البحثية بعنوان "مبدأ غاوس في التدفق غير القابل للانضغاط: منظور متغير موحد للضغط والإسقاط" للباحث كارتيك دورايسامي.

1. بيان المشكلة

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

• ما الذي يحدده المبدأ بالضبط عند لحظة زمنية ثابتة (اختيار حالة التسارع أم السرعة).
• كيف يرتبط المبدأ بطرق الإسقاط العددية الكلاسيكية (مثل طريقة شورين).
• التفسير الفيزيائي لتفكيك الضغط إلى مكونات "مفروضة" (خارجية) و"رد فعل" (قيد).
• ما إذا كان المبدأ يحل بطبيعته حالات عدم التحديد في التدفق الكامن المستقر (مثل اختيار الدوران/الدوامية).

يهدف المؤلف إلى توضيح أن مبدأ غاوس هو بيان لحظي متغير، يفرض قيوداً حركية (كينماتيكية) على التسارع، وليس آلية لاختيار حقول السرعة في الحالة المستقرة أو الدوران.

2. المنهجية

تستخدم الورقة نهج حساب التغيرات (Variational Calculus) الصارم المطبق على ميكانيكا السوائل عند زمن ثابت $t$:

• تجميد الحالة: يتم التعامل مع حقل السرعة $u(\cdot, t)$ كلقطة ثابتة ومعلومة. المتغير المتغير هنا هو التسارع المادي $a = \partial_t u + (u \cdot \nabla)u$، وتحديداً المشتق الزمني $\partial_t u$.
• دالة أبيل (The Appellian Functional): يعرّف المؤلف دالة تربيعية (دالة أبيل) تمثل معيار $L^2$ الموزون لانحراف التسارع:
  $$\mathcal{A}[u_t; u] = \frac{1}{2} \int_{\Omega} \rho |u_t + C|^2 dV$$
  حيث $C = (u \cdot \nabla)u + \frac{1}{\rho}\nabla P_F$ يمثل التسارع الحملي "المجمد" وأي قوى مفروضة ($P_F$).
• القيود: يتم الحل تحت قيود مستوى التسارع المستمدة من عدم الانضغاط ($\nabla \cdot u = 0$) وعدم نفاذية الجدران ($u \cdot n = 0$):
  $$\nabla \cdot u_t = 0 \quad \text{في } \Omega, \quad u_t \cdot n = 0 \quad \text{على } \partial\Omega$$
• التحسين (Optimization): يتم حل المشكلة باستخدام مضاعفات لاغرانج (شروط كاروش-كون-تاكر)، حيث يتم تحديد المضاعف المرتبط بقيد التباعد (divergence) على أنه ضغط رد الفعل ($P_R$).

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

أ. التوحيد مع إسقاط ليري-هودج (Leray-Hodge Projection)

تثبت الورقة أن شروط الأمثلية من الدرجة الأولى لتقليل دالة "أبيل" متطابقة رياضياً مع إسقاط ليري-هودج.

• يؤدي تقليل دالة "أبيل" إلى مسألة بواسون-نيومان لضغط رد الفعل $P_R$:
  $$\nabla^2 P_R = -\rho \nabla \cdot C, \quad \partial_n P_R = -\rho C \cdot n$$
• يعمل الحل $P_R$ كمضاعف لاغرانج يقوم بإسقاط الباقي المؤقت (الجزء غير التباعدي من التسارع) على فضاء الحقول التدرجية، مما يزيل التباعد والتدفق العمودي على الجدران.
• يثبت هذا أن خوارزمية "التنبؤ-بواسون-التصحيح" القياسية في الهيدروديناميكا العددية هي تجسيد منفصل مباشر لمبدأ غاوس.

ب. توضيح تفكيك الضغط

يعرف المؤلف بدقة الضغط الكلي $p$ كجموع لـ ضغط مفروض ($P_F$) وضغط رد فعل ($P_R$):

• $P_F$: يمثل قوى الجسم المحافظة (مثل الجاذبية) المفروضة بشكل مستقل عن قيود التدفق.
• $P_R$: هو مضاعف لاغرانج الفريد الذي يفرض عدم الانضغاط والتماس مع الجدران. وهو لا يبذل شغلاً على الحركات خالية التباعد والمتماسّة مع الجدران.
• الثبات: تثبت الورقة أن إزاحة جهد قياسي بين $P_F$ و$P_R$ هي خاصية جبرية ثابتة لشروط الأمثلية. ومع ذلك، فيزيائياً، $P_R$ هو قوة القيد الحقيقية فقط. التفكيك يكون فريداً بمجرد تعريف $P_F$ فيزيائياً.

ج. التمييز بين الإسقاط اللحظي واختيار الحالة

تتمثل إحدى المساهمات الحاسمة في توضيح نطاق المبدأ:

• الطبيعة اللحظية: يعمل مبدأ غاوس عند زمن ثابت على حقل سرعة مجمد. إنه يحدد تصحيح التسارع المطلوب لاستيفاء القيود.
• عدم اختيار الدوران: في التدفق المستقر، يعطي المبدำ $u_t^* = 0$ لأي حالة مستقرة مقبولة حركياً (بغض النظر عن الدوران $\Gamma$). لذلك، المبدأ لا يختار بحد ذاته الدوران الصحيح أو نقاط الركود.
• التحسين الثانوي: يتطلب اختيار الدوران (مثل استيفاء شرط كوتا) خطوة تحسين ثانوية عبر فضاء الحالات المستقرة، باستخدام دالة "أبيل" (أو تكاليف مشابهة) كدالة، كما اقترح غونزاليس وتيها (2022). خطوة غاوس نفسها توفر فقط الضغط المتوافق مع حالة معطاة.

د. التعامد من نوع ديريكليه كضبط للمقياس (Gauge Fix)

تتناول الورقة حرية التمثيل في تقسيم الضغط. وتوضح أن فرض تعامد ديريكليه ($\int \nabla P_R \cdot \nabla P_F dV = 0$) هو عرف رياضي (ضبط للمقياس) لاختيار ممثل فريد للتقسيم، بما يشبه تفكيك هيلمولتز، ولكنه ليس قيداً فيزيائياً على القوى.

4. النتائج

• اشتقاق معادلات أويلر: يستعيد الاشتقاق المتغير معادلات أويلر اللحظية ($\rho a = -\nabla p$) كشرط للاستقرار.
• المقياس التشخيصي: يعمل الحد الأدنى لقيمة دالة "أبيل" $\mathcal{A}^* = \frac{1}{2\rho} \|\nabla P_R\|_{L^2}^2$ كتشخيص كمي. وهو يقيس "جهد الإسقاط" المطلوب لفرض القيود.
  • إذا كان $\mathcal{A}^* = 0$، فإن التحديث المؤقت متوافق بالفعل مع القيود.
  • تشير القيم الكبيرة لـ $\mathcal{A}^*$ إلى وجود تباعد كبير أو عدم توافق في التدفق العمودي على الجدران، وهو أمر مفيد للكشف عن الأخطاء العددية أو عدم اتساق الحدود.
• التفرد الفيزيائي: ثبت أن ضغط رد الفعل $P_R$ هو قوة القيد الوحيدة (حتى ثابت معين) بمجرد تحديد القوى المفروضة، مما يحل النقاشات حول عدم تفرد الضغط في التدفقات غير القابلة للانضغاط.

5. الأهمية

• الوضوح النظري: تحل الورقة الارتباك في الأدبيات الحديثة من خلال الفصل الصارم بين فرض القيود اللحظي (غاوس/الإسقاط) وبين اختيار الحالة المستقرة (التحسين الثانوي). وهي توضح أن مبدأ غاوس يشرح كيف يحافظ الضغط على عدم الانضغاط لحظة بلحظة، وليس أي حالة تدفق مستقرة سيشغلها النظام.
• الرؤية العددية: توفر أساساً متغيراً صارماً لطرق الإسقاط المستخدمة في ديناميكا السوائل الحسابية (CFD)، وتربطها مباشرة بمبادئ الميكانيكا الكلاسيكية.
• الفائدة الحسابية: يوفر تحديد الحد الأدنى لدالة "أبيل" كمعيار $L^2$ لتدرج الضغط أداة تشخيصية عددية عملية لمراقبة جودة الحلول العددية واتساق الشروط الحدودية.
• إطار موحد: توحد التعامل مع القوى المفروضة وقوى رد الفعل، مما يوفر إطاراً واضحاً للنمذجة منخفضة الرتبة والمحاكاة المدعومة بالبيانات حيث يمكن "فرض" فيزياء معينة بينما تظل القيود دقيقة.

باخت-صار، يعيد عمل دورايسامي ترسيخ مبدأ غاوس كالتفسير المتغير الأساسي لمعادلة بواسون للضغط في التدفق غير القابل للانضغاط، مؤكداً على دوره كمصحح حركي لحظي وليس كمنتقٍ لطوبولوجيا التدفق العالمية.