Fully discrete error analysis of finite element discretizations of time-dependent Stokes equations in a stream-function formulation

بيان المشكلة
تتناول هذه الورقة تحليل الخطأ للتقريب المنفصل كلياً باستخدام العناصر المحدودة لمعادلات ستوكس المعتمدة على الزمن في بُعدين، وتحديداً باستخدام صياغة دالة الانسياب (stream-function). تتضمن مشكلة ستوكس الكلاسيكية إيجاد السرعة $\mathbf{u}$ والضغط $p$ اللذين يحققان معادلات الزخم والاستمرارية. في نهج دالة الانسياب، يتم تمثيل السرعة كدوران متجهي لدالة انسياب قياسية $\psi$ ($\mathbf{u} = \text{curl}(\psi)$)، مما يحقق شرط خلو التباعد تلقائياً. يؤدي هذا التحويل إلى اختزال النظام إلى معادلة معلمية (parabolic) من الرتبة الرابعة تتعلق بمؤثر بلهامي (biharmonic operator $\Delta^2$) والمشتق الزمني لـ لابلاس ($\partial_t \Delta$).

يسلط المؤلفون الضوء على قيد حرج في طرق العناصر المحدودة القياسية (مثل عناصر Taylor-Hood أو MINI) لمشكلات ستوكس: تقديرات الخطأ للسرعة غالباً ما تعتمد على تقريب الضغط. وبناءً على ذلك، فإن الاضطرابات في مصدر القوة الناتجة عن مجال تدرجي قد لا يتم استيعابها بواسطة الضغط المنفصل، مما يؤدي إلى أخطاء غير مثالية في السرعة. توفر صياغة دالة الانسياب مساراً نحو طرق "متينة تجاه الضغط" (pressure-robust)، ولكن كان هناك نقص في النظرية الصارمة لتقديرات الخطأ المنفصلة كلياً (التي تجمع بين التقطيع المكاني والزمني) لهذه المشكلة المحددة من الرتبة الرابعة المعتمدة على الزمن.

المنهجية
تستخدم الورقة إطار عمل "غاليركن" (Galerkin) منفصل كلياً يتكون من العناصر التالية:

1. التقطيع الزمني: يستخدم المؤلفون طريقة "غاليركن المنفصلة" (dG) من درجة $r$ اختيارية. يتضمن ذلك تقسيم الفاصل الزمني إلى فترات فرعية واستخدام كثيرات حدود مجزأة، مما يسمح بوجود قفزات عند عقد الزمن.
2. التقطيع المكاني: يتم تقديم التقطيع المكاني في إطار عام بدلاً من حصره في نوع عنصر محدد. يعتمد التحليل على صيغة ثنائية خطية إهليلجية متماثلة $a_h(\cdot, \cdot)$ ومؤثر منفصل $A_h$. تم تصميم الإطار لاستيعاب طرق متنوعة، بشرماً أن تحقق شروط تعامد غاليركن المحددة. تغطي الورقة صراحةً:
   • العناصر المتوافقة من الفئة $C^1$ (مثل عناصر Argyris).
   • طرق الجزاء الداخلي من الفئة $C^0$ (باستخدام عناصر لاغرانج القياسية من الدرجة $\ge 2$).
3. تحليل الانتظام (Regularity Analysis): يُخصص جزء كبير من العمل لإثبات انتظام الحل المستمر $\psi$. يثبت المؤلفون وجود، ووحدانية، ونتائج استقرار للمشكلة المستمرة تحت افتراضات دنيا للبيانات ($f \in L^2(I; H^{-1}(\Omega))$ و $\psi_0 \in H^2_0(\Omega)$). ومن الأهمية بمكان أنهم يوضحون أنه بينما يمتلك الحل خصائص استقرار معينة، فإن الانتظام الأعلى للمشتق الزمني $\partial_t \psi$ (تحديداً في $L^2$ أو $H^{-1}$) لا يمكن ضمانه تحت هذه الافتراضات الدنيا، وهي حقيقة أُثبتت عبر مثال مضاد في الملحق (أ).

المساهمات والنتائج الرئيسية
المساهمة الأساسية هي وضع تقديرات خطأ من نوع أفضل تقريب (best-approximation-type) للحل المنفصل كلياً $\psi_{kh}$. تنص النتيجة الرئيسية (النظرية 4.5) على أن الخطأ في تدرج دالة الانسياب محكوم بمجموع ثلاثة أخطاء تقريب:
$$\|\nabla(\psi - \psi_{kh})\|_{L^2(I \times \Omega)} \le C \left( \|\nabla(\psi - \chi_{kh})\|_{L^2(I \times \Omega)} + \|\nabla(\psi - R_h\psi)\|_{L^2(I \times \Omega)} + \|\nabla(\psi - \pi_k\psi)\|_{L^2(I \times \Omega)} \right)$$
حيث أن:

• $\chi_{kh}$ هو عنصر عشوائي من الفضاء المنفصل.
• $R_h$ هو الإسقاط المرتبط بمشكلة بلهام (biharmonic) الساكنة.
• $\pi_k$ هو إسقاط زمني محدد.

من خلال تطبيق خصائص التقريب القياسية لهذه الإسقاطات تحت افتراضات الانتظام المثبتة، يستنتج المؤلفون معدلات تقارب مثالية في النطاقات المضلعة المحدبة:
$$\|\nabla(\psi - \psi_{kh})\|_{L^2(I \times \Omega)} \le C(k + h^2) \left( \|f\|_{L^2(I; H^{-1}(\Omega))} + \|\psi_0\|_{H^2_0(\Omega)} \right)$$
وهذا يعني تقارباً من الرتبة $O(k)$ في الزمن و $O(h^2)$ في الفضاء لخطأ السرعة (بما أن $\|\mathbf{u} - \mathbf{u}_{kh}\| = \|\nabla(\psi - \psi_{kh})\|$).

الأهمية والادعاءات
تضع هذه الورقة نفسها كخطوة تأسيسية في التحليل العددي لصياغات دالة الانسياب للمشكلات المعتمدة على الزمن.

• المتانة تجاه الضغط (Pressure Robustness): لا تحتوي التقديرات المستمدة على حدود الضغط، مما يعالج المشكلة المعروفة حيث تعاني الطرق القياسية من أخطاء في السرعة تعتمد على تقريبات الضغط. وتؤكد التجارب العددية أن طريقة دالة الانسياب تظل متينة حتى عندما يتم تعديل مصدر القوة بواسطة مجال تدرجي، بينما تظهر الطرق المختلطة القياسية (عناصر MINI) أداءً متدهوراً.
• افتراضات الحد الأدنى من الانتظام: لا يتطلب التحليل افتراضات انتظام إضافية تتجاوز ما توفره البيانات والنطاق بشكل طبيعي. يجادل المؤلفون صراحةً بأن الأعمال السابقة غالباً ما فرضت افتراضات صارمة للغاية (مثل اشتراط كون $\partial_{tt}\psi$ محدوداً) لتحقيق تقديرات مماثلة.
• العمومية: الإطار قابل للتطبيق على مجموعة واسعة من طرق التقطيع المكاني، بشرط أن تحقق شروط التعامد اللازمة، مما يجعل النتائج متعددة الاستخدامات لطرق العناصر المحدودة المتوافقة $C^1$ وطرق الجزاء الداخلي غير المتوافقة $C^0$.
• التحكم الأمثل: يشير المؤلفون إلى أن افتراضات الانتظام المنخفضة تجعل هذه النتائج مناسبة للتطبيق في مشكلات التحكم الأمثل.

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