Steady Incremental Viscosity Splitting Method for solving the stationary Navier-Stokes equation

ملخص تقني: طريقة تقسيم اللزوجة المتزايدة المستقرة لمعادلات نافييه-ستوكس المستقرة

بيان المشكلة
تتناول الورقة البحثية الحل العددي لمعادلات نافييه-ستوكس المستقرة وغير القابلة للانضغاط مع تثبيت (grad-div). المعادلات الحاكمة هي:
$$-\nu\Delta \mathbf{u} + \mathbf{u} \cdot \nabla \mathbf{u} + \nabla p - \gamma\nabla(\nabla \cdot \mathbf{u}) = \mathbf{f}, \quad \nabla \cdot \mathbf{u} = 0$$
في نطاق $\Omega$، مع شروط حدودية ديريكليه متجانسة. يشير المؤلفون إلى أن حل الأنظمة الخطية الناتجة ذات نقطة السرج غير المتماثلة (الناشئة عن خطية بيكار القياسية) مكلف حاسوبياً، لا سيما في المسائل ثلاثية الأبعاد حيث تفشل الحلول المباشرة بسبب قيود الذاكرة وتواجه الحلول التكرارية صعوبة في التعامل مع بنية نقطة السرج. وبينما تُعد طرق التقسيم (splitting methods) مستخدمة جيداً في التدفقات غير المستقرة، فإن تطبيقها على المسائل المستقرة تماماً أقل تطوراً.

المنهجية
يقترح المؤلفون طريقة تقسيم اللزوجة المتزايدة المستقرة (SIVS). تعتمد هذه الطريقة على تكييف مخطط تقسيم اللزوجة المتزايدة (IVS)، المستخدم سابقاً في التدفقات غير المستقرة، لتناسب سياق الحالة المستقرة. تعمل الطريقة كـمخطط نقطة ثابت تكراري:

1. الخطوة 1 (تنبؤ السرعة): بالاعتماد على التكرار السابق $(\mathbf{u}^{k-1}, p^{k-1})$، يتم حل معادلة تفاضلية جزئية إهليلجية للحصول على سرعة وسيطة $\tilde{\mathbf{u}}^k$:
   $$a(\tilde{\mathbf{u}}^k, \mathbf{v}) + c^*(\mathbf{u}^{k-1}, \tilde{\mathbf{u}}^k, \mathbf{v}) + b(p^{k-1}, \mathbf{v}) = \langle \mathbf{f}, \mathbf{v} \rangle$$
   هنا، يتضمن الشكل الثنائي الخطية $a(\cdot, \cdot)$ الحد اللزج وحن حد التثبيت (grad-div) ($\gamma$).

2. الخطوة 2 (تصحيح الضغط والإسقاط): حل السرعة المحدثة $\mathbf{u}^k$ والضغط $p^k$ من خلال تصحيح الحل الوسيط:
   $$a(\mathbf{u}^k - \tilde{\mathbf{u}}^k, \mathbf{v}) + b(p^k - p^{k-1}, \mathbf{v}) = 0$$
   $$b(q, \mathbf{u}^k) = 0$$
   تضمن هذه الخطوة أن تكون السرعة النهائية خالية من التباعد (divergence-free) وتلبي الشروط الحدودية.

التفسير الجبري
يفسر البحث طريقة SIVS كطريقة تقسيم جبري. من خلال إجراء تحليل LU غير كامل (ILU) لمصفوفة النظام، يوضح المؤلفون أن SIVS تقابل تحليل كتلة (block factorization) محدد حيث تكون مصفوفة متمم شور $S = -B\tilde{A}^{-1}B^T$ (حيث $\tilde{A}$ تمثل مؤثر اللزوجة الخطي) مكافئة طيفياً لمصفوفة كتلة الضغط. والأهم من ذلك، أنه بخلاف تكرارات بيكار القياسية التي تتطلب عكس متمم شور مختلف في كل خطوة، تتطلب طريقة SIVS حل نظام باستخدام مصفوفة (SPD) ثابتة لمتغير الضغط عبر جميع التكرارات غير الخطية. وهذا يسمح بالحساب المسبق وإعادة استخدام الممهدات (preconditioners)، مما يقلل التكلفة الحسابية بشكل كبير.

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

• مخطط مبتكر: تقديم طريقة SIVS لمعادلات نافييه-ستوكس المستقرة مع تثبيت grad-div، مما يضمن أن مجال السرعة المنفصل خالٍ من التباعد ويلبي شروط ديريكليه الكاملة.
• تحليل التقارب: أثبت المؤلفون الحد المكتوم الموحد (uniform boundedness) للتكرارات والتقارب الهندسي للمخطط تحت شرط صغر البيانات ($\Lambda < 1/\sqrt{2}$).
• الكفاءة الجبرية: أُثبت أن طريقة SIVS هي نهج تقسيم جبري يعمل على فك ارتباط حلول السرعة والضغط إلى أنظمة خطية أبسط مقارنة بنظام نقطة السرج المقترن القياسي.

النتائج العددية
تحقق المؤلفون من صحة الطريقة باستخدام FreeFem++ وعناصر Taylor-Hood ($P_2, P_1$) عبر ثلاث حالات اختبار:

1. اختبار التقارب: باستخدام حل مُصنّع، حققت الطريقة معدلات التقارب المثلى المتوقعة للسرعة والضغط وأخطاء التباعد. ظل عدد التكرارات غير الخطية ثابتاً عبر تحسينات الشبكة (mesh refinements).
2. تجربة تجويف المحرك ثنائي الأبعاد (2D Lid-Driven Cavity): تم اختبار المخطط عند أرقام رينولدز تصل إلى $Re=5000$. طابقت نتائج ملفات السرعة والخطوط الانسيابية البيانات المرجعية (Ghia et al.) بدقة. كانت أعداد التكرارات قابلة للمقارنة مع طريقة (Incremental Picard-Yosida - IPY)، لكن المؤلفين يتوقعون وقتاً أقل في المعالجة المركزية (CPU time) بسبب الأنظمة الخطية الأبسط.
3. تدفق ثلاثي الأبعاد حول أسطوانة: تم حل مسألة مرجعية عند $Re \approx 450$ (بناءً على $\nu=10^{-3}$). تقارب المخطط في 11 تكراراً غير خطياً، منتجاً معاملات السحب والرفع وانخفاض الضغط ضمن الفترات المرجعية المذكورة في الأدبيات.
4. تأثير المعاملات: أشارت الاختبارات على معامل grad-div $\gamma$ إلى أن القيم الأكبر لـ $\gamma$ تسرع التقارب إلى الحالة المستقرة، رغم أنها قد تؤدي إلى تدهور حالة التكييف (conditioning) للنظام الخطي في الخطوة 1.

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