A $C^0$ weak Galerkin method with preconditioning for constrained optimal control problems with general tracking

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "طريقة 'ويك جاليركن' من الفئة $C^0$ مع المسبق الشرطي لمسائل التحكم الأمثل المقيدة ذات التتبع العام."

1. بيان المسألة

تتناول هذه الورقة مسائل التحكم الأمثل (OCPs) المحكومة بمعادلات تفاضلية جزئية (PDEs) إهليلجية. وتحديداً، تركز على المسائل التي تتميز بـ:

• دوال تكلفة تتبع عامة: الهدف هو تقليل الفرق بين الحالة $y$ والحالة المرغوبة $y_d$، حيث يمكن تعريف $y_d$ فوق نقاط، أو منحنيات، أو مناطق فرعية (وليس فقط فوق النطاق بأكمله).
• قيود الحالة عند النقاط: يجب أن تحقق متغير الحالة $y(x)$ حدوداً دنيا وعليا ($\psi_-(x) \leq y(x) \leq \psi_+(x)$) في كل مكان تقريباً في النطاق.
• الصياغة الرياضية: من خلال استبعاد متغير التحكم عبر قيد المعادلة التفاضلية الجزئية، يتم اختزال المسألة إلى متباينة تباينية من الدرجة الرابعة تتضمن متغير الحالة.
• التحديات:
  • انخفاض الانتظام: وجود مصطلحات التتبع العامة (المقاييس) والقيود النقطية يؤدي إلى تدهور انتظام الحل الأمثل ومضاعفات لاغرانج، مما يؤدي غالباً إلى مضاعفات على شكل مقاييس (measures).
  • التعقيد العددي: تتطلب حلول المسائل من الدرجة الرابعة عادةً عناصر ذات استمرارية $C^1$، وهي معقدة التنفيذ. وتتطلب طرق "عقوبة الجزئي الداخلي" الحالية ($C^0\text{-IP}$) ضبطاً دقيقاً لمعاملات العقوبة وتتضمن حدود أثر (trace terms) تعقد عملية التجميع لكل عنصر.
  • سوء التكيّف (Ill-Conditioning): تكون الأنظمة الخطية الناتجة سيئة التكيّف بشكل حاد، خاصة عند وجود مصطلحات تتبع عامة، مما يعيق أداء الحلول التكرارية.

2. المنهجية

أ. طريقة $C^0$ ويك جاليركن ($C^0\text{-WG}$)

يقترح المؤلفون طريقة $C^0$ ويك جاليركن باستخدام عناصر لاغرانج التربيعية ($P_2$) المستمرة عالمياً.

• لابلاس الضعيف (Weak Laplacian): بدلاً من استخدام المشتقات القوية، تستخدم الطريقة لابلاس ضعيف متقطع ($\Delta_w$). بالنسبة للدالة $v_h$، يتم حساب لابلاس الضعيف على العنصر $T$ كدالة ثابتة مجزأة ($P_0$).
  • يتم اشتقاقه صراحةً باستخدام قيم الدوال عند مراكز الثقل للحواف/الأوجه المشتركة والبيانات الهندسية الأساسية (أطوال الحواف، مساحات العناصر).
  • الصيغة: $\Delta_w v_h|_T = \frac{1}{|T|} \sum_{e \subset \partial T} \frac{|e|}{2} (\nabla v_h^+ + \nabla v_h^-) \cdot n_T$.
• الصيغة الثنائية الخطية: تتكون الصيغة الثنائية الخطية المتقطعة $a_w(\cdot, \cdot)$ من:
  1. الضرب الداخلي لـ $L^2$ للابلاسات الضعيفة: $(\Delta_w u, \Delta_w v)_{\mathcal{T}_h}$.
  2. حد استقرار خالٍ من المعاملات (Parameter-free): $\frac{1}{4} \langle h_e^{-1} [\nabla u] \cdot n_e, [\nabla v] \cdot n_e \rangle_{\mathcal{E}_h^o}$.
• المزايا مقارنة بـ $C^0\text{-IP}$:
  • لا يوجد معامل عقوبة: على عكس طريقة $C^0\text{-IP}$، التي تتطلب معامل عقوبة $\rho$ كبيراً ومضبوطاً لضمان الاستقرار، فإن طريقة $C^0\text{-WG}$ خالية من المعاملات.
  • تجميع مبسط: تلغي هذه الصياغة حدود الأثر المعقدة (التي تتضمن متوسطات التدرجات) الموجودة في $C^0\text{-IP}$، مما يسمح بتجميع فعال لمصفوفة الصلابة لكل عنصر على حدة.

ب. المسبق الشرطي الجمعي لـ "Schwarz" (ASP)

لمعالجة سوء التكيّف في الأنظمة الخطية الناشئة عن صياغة الدرجة الرابعة:

• المنهج: تم تطوير مسبق شرطي جمعي لـ "Schwarz" من مستوى واحد.
• الآلية: يتم تقسيم النطاق العالمي إلى نطاقات فرعية متداخلة. يتم تفكيك المسألة العالمية إلى مسائل فرعية محلية تُحل بشكل مستقل، ثم تُجمع حلولها جمعياً.
• الحد النظري: ثبت أن رقم الحالة (condition number) $\kappa(B_h A_h)$ للنظام المسبق الشرطي محدود بـ $O(\delta^{-4})$، حيث $\delta$ هو حجم التداخل. ومع وجود نطاقات فرعية منتظمة الشكل، يمكن تحسين ذلك إلى $O(H^{-1}\delta^{-3})$.

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

1. التقطيع المبتكر ($C^0\text{-WG}$):

   • تطوير طريقة $C^0\text{-WG}$ مخصصة لمسائل التحكم الأمثل المقيدة ذات التتبع العام.
   • استخدام العناصر التربيعية المستمرة عالمياً لتبسيط حساب لابلاس الضعيف.
   • إلغاء معاملات العقوبة وحدود الأثر، مما يؤدي إلى تنفيذ خالٍ من المعاملات، وقوي، وفعال حاسوبياً.
   • تحليل خطأ صارم يثبت جودة المسألة (well-posedness) ومعدلات التقارب.

2. التمهيد الشرطي المتخصص:

   • تصميم مسبق شرطي جمعي لـ "Schwarz" مصمم خصيصاً لصياغة $C^0\text{-WG}$.
   • إثبات أن المسبق الشرطي يقلل رقم الحالة بشكل كبير، مما يسمح بحلول تكرارية قابلة للتوسع للأنظمة سيئة التكيّف النموذجية لمسائل التحكم الأمثل من الدرجة الرابعة.

3. تحليل التقارب الصارم:

   • إثبات رتب التقارب المثلى في المعايير المعتمدة على الشبكة.
   • إثبات التكافؤ بين معيار $C^0\text{-WG}$ ومعايير "سوبولف المكسورة" (broken Sobolev norms) القياسية المستخدمة في تحليل $C^0\text{-IP}$.

4. النتائج العددية

أجرى المؤلفون تجارب عددية في بيئات ثنائية الأبعاد لكل من مسائل "بيهارمونيك" (biharmonic) ومسائل التحكم الأمثل المقيدة ذات التتبع العام.

• مسألة بيهارمونيك:
  • حققت طريقة $C^0\text{-WG}$ تقارباً من الدرجة الأولى (ولوحظ وجود تقارب فائق/superconvergence) في المعيار المعتمد على الشبكة، وهو ما يتوافق مع التوقعات النظرية ويماثل طريقة $C^0\text{-IP}$.
• مسألة التحكم الأمثل المقيدة ذات التتبع العام:
  • الدقة: أنتجت طريقة $C^0\text{-WG}$ باستمرار تقريبات أكثر دقة من طريقة $C^0\text{-IP}$ عبر معايير $L^\infty$ و $L^2$ و $H^1$.
  • معدلات التقارب: أظهرت $C^0\text{-WG}$ رتب تقارب أعلى في معايير $L^\infty$ و $L^2$ و $H^1$ مقارنة بـ $C^0\text{-IP}$.
• أداء المسبق الشرطي:
  • أرقام الحالة: بدون المسبق الشرطي، نما رقم الحالة كـ $O(h^{-4})$.
  • مع استخدام ASP: تم تخفيف نمو رقم الحالة بشكل كبير (بمقدار أكثر من رتبتين عشريتين)، حيث نما تقريباً كـ $O(h^{-3})$ أو أفضل اعتماداً على التداخل $\delta$.
  • تأثير التداخل: أدى زيادة حجم التداخل ($\delta = 2h$ مقابل $\delta = h$) إلى تقليل رقم الحالة بشكل أكبر، مما يؤكد فعالية استراتيجية تفكيك النطاق.

5. الأهمية

تمثل هذه الورقة تقدماً كبيراً في الحل العددي لمسائل التحكم الأمثل عالية الدرجة:

• الكفاءة: من خلال إزالة الحاجة إلى ضبط معامل العقوبة وتبسيط تجميع المصفوفات، تقلل طريقة $C^0\text{-WG}$ من تعقيد التنفيذ والأعباء الحسابية.
• القوة (Robustness): تتعامل الطريقة بفعالية مع انخفاض انتظام الحلول (الناتج عن القيود النقطية والتتبع العام) دون التضحية بالاستقرار.
• القابلية للتوسع: يعمل دمج المسبق الشرطي الجمعي لـ "Schwarz" على حل مشكلة سوء التكيّف الحرجة، مما يجعل الطريقة مناسبة للمحاكاة واسعة النطاق والحوسبة المتوازية.
• الإطار العملي: يضع هذا العمل إطاراً حاسوبياً قوياً يمكن توسيعه ليشمل سيناريوهات أكثر تعقيداً، بما في ذلك المسائل المعتمدة على الزمن والمسائل غير الخطية.