Equivalence of mixed and nonconforming methods on general polytopal partitions. Part I: Multiscale and projection methods

إليك ملخص تقني مفصل للورقة البحثية بعنوان "تكافؤ الطرق المختلطة وغير المتوافقة على التقسيمات متعددة الأوجه: الجزء الأول: الطرق متعددة المقاييس وطرق الإسقاط" للباحث سيمون ليمير (Simon Lemaire).

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

تتقصى هذه الورقة التكافؤ الرياضي بين طرق العناصر المحدودة المختلطة (القائمة على التدفق - flux-based) والطرق الأولية غير المتوافقة (القائمة على الجهد - potential-based) لمسائل الانتشار المتغير على تقسيمات متعددة الأوجه العامة (general polytopal partitions).

المسألة المستمرة الأساسية (المسألة 0) هي معادلة الانتشار المتغير:
$$\mathbf{b}^{-1}\boldsymbol{\sigma} - \nabla u = 0, \quad -\text{div}\,\boldsymbol{\sigma} = f, \quad u|_{\partial\Omega} = 0$$
حيث $\boldsymbol{\sigma}$ هو التدفق، و $u$ هو الجهد، و $\mathbf{b}$ هو حقل موتر الحركية (mobility tensor field) متماثل وموجب التعريف.

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

• هل يمكننا بناء فضاءات ذات أبعاد محدودة تكون فيها الصيغ الأولية والمختلطة متكافئة؟
• كيف يظل هذا التكافؤ قائماً عندما يتم الحفاظ على العلاقة التكوينية ($\boldsymbol{\sigma} = \mathbf{b}\nabla u$) بدقة مقابل الحالات التي يتم فيها إسقاطها؟
• كيف ترتبط هذه البناءات النظرية بالطرق العملية متعددة المقاييس وطرق الإسقاط؟

2. المنهجية

يستخدم المؤلف نهج بناء الفضاء الافتراضي (virtual space construction)، منتقلاً من التحليل المستمر ذي الأبعاد اللانهائية إلى التحليل المنفصل ذي الأبعاد المحدودة.

أ. المستوى المستمر: التهجين (Hybridization)

تؤسس الورقة أولاً "حجر رشيد" (Rosetta stone) على المستوى المستمر عبر إدخال صيغ مهجنة:

1. الأولي المهجن (المسألة 3): يخفف من شرط التوافق في $H^1_0(\Omega)$، ويقدم أثر تدفق (flux trace) على الهيكل (skeleton).
2. المختلط المهجن (المسألة 4): يخفف من شرط التوافق في $\mathbf{H}(\text{div}; \Omega)$، ويقدم أثر جهد (potential trace) على الهيكل.
   باستخدام نظرية مضاعفات لاجرانج، يثبت المؤلف أن:

• الأولي $\iff$ الأولي المهجن
• المختلط $\iff$ المختلط المهجن
• بشكل حاسم: الأولي المهجن $\implies$ المختلط، والمختلط المهجن $\implies$ الأولي.
  هذا يثبت مخطط تكافؤ كاملاً (الشكل 1) يربط جميع الصيغ الأربع عبر العلاقة التكوينية $\boldsymbol{\sigma} = \mathbf{b}\nabla u$.

ب. المستوى المنفصل: البناء الافتراضي

يعرّف المؤلف نظائر ذات أبعاد محدودة باستخدام فضاءات متعددة الحدود مكسورة (broken polynomial spaces) وفضاءات العناصر الافتراضية:

• فضاء التدفق ($\boldsymbol{\Sigma}_{\ell,k}$): مُعرَّف محلياً بحيث يكون $\text{div}\,\boldsymbol{\tau} \in \mathbb{P}_{\ell-1}$ وآثار عمودية (normal traces) تنتمي إلى $\mathbb{P}_{k-1}$.
• فضاء الجهد ($U_{\ell,k}$): مُعرَّف ضمنياً عبر معادلة الـ PDE: حيث $v \in H^1(K)$ بحيث يكون $\text{div}(\mathbf{b}\nabla v) \in \mathbb{P}_{\ell-1}$ وآثار التدفق العمودي تنتمي إلى $\mathbb{P}_{k-1}$.
• الهوية الرئيسية: يضمن هذا البناء أن $\boldsymbol{\Sigma}_{\ell,k}(K) = \mathbf{b}\nabla U_{\ell,k}(K)$. وهذا يعني أن العلاقة التكوينية المنفصلة محفوظة بدقة.

يثبت المؤلف أن النسخ المنفصلة من المسائل المهجنة متكافئة، مما يؤسس لـ التكافؤ الأولي/المختلط المنفصل لهذه الفضاءات الافتراضية.

ج. التنفيذ العملي: طرق الإسقاط (Projection Methods)

بما أن الفضاءات $U_{\ell,k}$ هي عموماً غير متعددة حدود (افتراضية) ولا يمكن حسابها مباشرة، تبحث الورقة في سيناريوهين:

1. الحفظ الدقيق (الطرق متعددة المقاييس): في حالات محددة (مثل طريقة Raviart-Thomas من الدرجة الأدنى على الـ simplices مع $\mathbf{b}$ ثابت داخل الخلية)، يتقلص الفضاء الافتراضي إلى فضاء متعدد حدود معروف معزز بـ "فقاعات" (bubbles) (مثل Crouzeix-Raviart المعزز بفقاعات تربيعية). وهذا يستعيد التكافؤات المعروفة (مثل التكافؤ بين MsHHO و MHM).
2. الإسقاط (الحالة العامة): عندما يكون فضاء التدفق متعدد حدود ولكن العلاقة التكوينية لا يمكن حفظها بدقة (أي $\mathbf{b}\nabla U \not\subset \boldsymbol{\Sigma}$)، يقدم المؤلف طرق الإسقاط.
   • يُعاد صياغة المسألة الأولية باستخدام إسقاط متعامد موزون للتدرج على فضاء التدفق.
   • المسألة الناتجة (المسألة 45) هي "طريقة عناصر محدودة إسقاطية" حيث قد يظل فضاء الجهد افتراضياً بحتاً (معرفاً فقط بدرجات الحرية)، بينما يكون فضاء التدفق قابلاً للحساب.

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

1. إطار توحيد شامل: توفر الورقة برهاناً صارماً للتكافؤ بين الطرق الأولية والمختلطة على شبكات متعددة الأوجه عامة مع معاملات انتشار متغيرة. وهي توحد مختلف الطرق متعددة المقاييس (MsFE, MsHHO, MHM) تحت مظلة نظرية واحدة.
2. بناء الفضاء الافتراضي: تقدم الوربة بناءً افتراضياً محدداً ($\boldsymbol{\Sigma}_{\ell,k}$ و $U_{\ell,k}$) يحفظ العلاقة التكوينية $\boldsymbol{\sigma} = \mathbf{b}\nabla u$ بدقة. وهذا يسد الفجوة بين الطرق متعددة المقاييس والعناصر المحدودة الكلاسيكية.
3. تعميم نتائج "الرصاصة الفضية" (Silver Bullet): توضح الورقة أن التكافؤ بين طريقة Raviart-Thomas من الدرجة الأدنى (المختلطة) وطريقة Crouzeix-Raviert المعززة بالفقاعات (الأولية) ينطبق على أي simplex (ليس فقط المنتظمة) ويمتد ليشمل متوازي المستطيلات مع موترات حركية قطرية.
4. إعادة تفسير الطرق ذات الاستقرار الذاتي: توضح الورقة أن طرق Hybrid High-Order (HHO) و Weak Galerkin (WG) ذات الاستقرار الذاتي هي في الأساس طرق إسقاط. ومن خلال رؤيتها عبر عدسة التكافؤ المستمد، تقدم الورقة معادلاتها المختلطة ذات العناصر المحدودة، والتي كانت غير معروفة أو غير واضحة سابقاً.
5. معيار عملي لصلاحية الحل (Well-posedness): تستنتج الورقة معياراً بسيطاً وقابلاً للفحص لصلاحية الحل والتكافؤ لطرق الإسقاط على الشبكات متعددة الأوجه العامة. يعتمد هذا المعيار على تطابق (surjectivity) خريطة التدفق (التي تربط التدفقات بتدرجها وآثارها العمودية)، وهو ما يعد إعادة صياغة لشرط LBB.

4. النتائج الرئيسية

• مبرهنة (التكافؤ المنفصل): بالنسبة للفضاءات الافتراضية المعرفة، المسألة المختلطة (المسألة 6) والمسألة الأولية (المسألة 5) متكافئتان. حل المسألة المختلطة يعطي تدفقاً $\boldsymbol{\sigma} = \mathbf{b}\nabla u_K$ وجهداً $u = \pi_{K}(u_K)$ يحل المسألة الأولية.
• تكافؤ الطرق متعددة المقاييس: تؤكد الورقة التكافؤ من المستوى الأول بين طريقة MsHHO (Multiscale Hybrid High-Order) وطريقة MHM (Multiscale Hybrid Mixed)، مما يوسع النتائج السابقة.
• تعزيز الفقاعات (Bubble Enrichment): ثبت أن طريقة Raviart-Thomas من الدرجة الأدنى على أي simplex مكافئة لطريقة Crouzeix-Raviart المعززة بفقاعة "ضعيفة" تربيعية محددة، دون اشتراط الانتظام الهندسي للشبكة.
• تكافؤ الإسقاط: بالنسبة لفضاءات BDM أو Raviart-Thomas من درجات أعلى حيث لا يمكن حفظ العلاقة التكوينية بدقة، فإن الطريقة المختلطة تكافئ طريقة إسقاط أولية (المعادلة 45). يوفر هذا أساساً نظرياً لطرق HHO و WG ذات الاستقرار الذاتي.

5. الأهمية

• التوحيد النظري: يوحد هذا العمل المجالات التي تبدو متباعدة مثل الطرق متعددة المقاييس، والطرق الهجينة، وطرق الإسقاط، مبيناً أنها أوجه مختلفة لنفس البنية التباينية الأساسية.
• البصيرة الخوارزمية: من خلال إرساء التكافؤ، تسمح الورقة للممارسين باستخدام الصيغة الأكثر ملاءمة. على سبيل المثال، يمكن حل مسألة أولية (غالباً ما تكون أسهل في التكييف/Preconditioning) ثم إعادة بناء تدفق متوافق، أو العكس.
• التحقق من صحة الطرق ذات الاستقرار الذاتي: توفر الورقة تبريراً صارماً للطرق الهجينة "الخالية من الاستقرار" (HHO/WG) عبر ربطها بنظريات العناصر المحددة المختلطة المفهومة جيداً. وهذا يثبت استقرارها وخصائص تقاربها.
• تمهيد للجزء الثاني: تضع هذه الورقة المسرح للجزء الثاني، والذي سيتناول طرق العناصر متعددة الأوجه (Polytopal Element Methods) حيث لا يكون فضاء التدفق العالمي متوافقاً في $\mathbf{H}(\text{div})$، مما يخفف من الافتراضات الموضوعة هنا.