Numerical solution of elliptic distributed optimal control problems with boundary value tracking

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "الحل العددي لمسائل التحكم الأمثل التوزيعي الإهليلجي مع تتبع القيم الحدية".

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

تتناول الورقة مسألة التحكم الأمثل لتتبع القيم الحدية (OCP) المقيدة بمعادلة تفاضلية جزئية إهليلجية ذات شروط حدودية من نوع نيومان.

• الهدف: تقليل دالة التكلفة $J(y_\varrho, u_\varrho)$ التي توازن بين خطأ التتبع للحالة $y_\varrho$ تجاه هدف $y$ على الحدود $\Gamma$، وحداً للتنظيم (regularization) للتحكم $u_\varrho$:
  $$J(y_\varrho, u_\varrho) = \frac{1}{2} \|y_\varrho - y\|^2_{L^2(\Gamma)} + \frac{1}{2} \varrho \|u_\varrho\|^2_{U}$$
• القيود: تخضع الحالة لمعضلة القيمة الحدية لنيومان:
  $$-\Delta y_\varrho + y_\varrho = u_\varrho \quad \text{في } \Omega, \quad \partial_n y_\varrho = 0 \quad \text{على } \Gamma$$
• التمييز الجوهري: على عكس الأعمال السابقة التي استخدمت تنظيم $L^2(\Omega)$، تبحث هذه الورقة في تنظيم الطاقة (energy regularization) حيث يتم اختيار فضاء التحكم ليكون $U = \tilde{H}^{-1}(\Omega) := [H^1(\Omega)]^*$. يضمن هذا الاختيار أن تكون الخريطة من الحالة إلى التحكم (state-to-control map) هي "تشاكل" (isomorphism) ضمن $H^1(\Omega)$، مما يتجنب الحاجة إلى دوال أساس ذات رتب أعلى كما هو مطلوب في تنظيم $L^2$.

2. المنهجية

أ. إعادة الصياغة التباينية القائمة على الحالة

يعيد المؤلفون صياغة مسألة التحكم الأمثل (OCP) إلى مسألة تباينية قائمة على الحالة. ومن خلال استخدام حجج الثنائية (duality arguments)، يثبتون أن معيار التحكم في الفضاء المزدوج يكافئ معيار $H^1(\Omega)$ للحالة. يسمح هذا بحذف متغير التحكم $u_\varrho$، مما يقلل المسألة إلى تقليل دالة تعتمد فقط على الحالة $y_\varrho$:
$$\tilde{J}(y_\varrho) = \frac{1}{2} \|y_\varrho - y\|^2_{L^2(\Gamma)} + \frac{1}{2} \varrho \|y_\varrho\|^2_{H^1(\Omega)}$$
ويحقق المقلل (minimizer) معادلة تدرج (5) تتضمن الضرب الداخلي لـ $L^2(\Gamma)$ و $H^1(\Omega)$.

ب. التجزئة بالعناصر المحدودة (Finite Element Discretization)

تركز الورقة على تجزئة العناصر المحدودة المتوافقة ذات الضرب التنسوري (conforming tensor-product FE) على نطاق مكعب الوحدة $\Omega = (0,1)^d$.

• دوال الأساس: يقوم المؤلفون ببناء فضاء عناصر محدودة محدد $Y_h \subset Y$ باستخدام دوال أساس خطية مجزأة معدلة. ومن الأهمية بمكان أن هذه الدوال مصممة بحيث يكون لها مشتق طبيعي متلاشٍ ($\partial_n \phi = 0$) على الحدود، مما يحقق شرط نيومان بشكل طبيعي دون الحاجة إلى مضاعفات لاغرانج.
• النظام المنفصل: تؤدي التجزئة إلى نظام خطي على الصورة:
  $$(\hat{M}_h + \varrho(\tilde{K}_h + \mathring{K}_h))y_h = y_h$$
  حيث تمثل المصفوفات مساهمات كتلة الحدود، وصلابة الحدود، وصلابة الداخل.

ج. الحلول السريعة

لحل النظام الخطي الناتج بكفاءة، يقترح المؤلفون نهج متمم شات (Schur Complement - SC):

1. الحذف: يتم حذف درجات الحرية الداخلية، مما يقلص النظام إلى نظام متمم شات حدودي $S_{BB}y_B = y_B$.
2. التكافؤ الطيفي: يثبت المؤلفون أنه بالنسبة لـ $\varrho \leq h$، فإن متمم شات $S_{BB}$ يكافئ طيفياً مصفوفة كتلة الحدود $M_{BB}$.
3. الحلول التكرارية:
   • يتم حل نظام SC باستخدام طريقة الاشتقاق المترافق (Conjugate Gradient - CG).
   • التحسين المسبق (Preconditioning): تُستخدم مصفوفة كتلة مجمعة (lumped mass matrix) كمحسن مسبق قطري لطريقة PCG، مما يضمن قياساً (scaling) أمثلاً.
   • المحلل الداخلي: يتم تقريب تأثير معكوس مصفوفة الصلابة الداخلية (المطلوب لحساب متمم شات) باستخدام طريقة الاشتقاق المترافق المسبقة بمتعدد الشبكات الجبري (Algebraic MultiGrid - AMG).

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

1. تقديرات الخطأ النظرية: تشتق الورقة تقديرات خطأ صارمة للحالة المنفصلة $y_{\varrho h}$ بالنسبة للهدف $y$. تعتمد معدلات التقارب على انتظام (regularity) الهدف:
   • إذا كان $y \in H^{1/2}(\Gamma)$، فإن الخطأ هو $O(h^{1/2})$ (مع $\varrho = h$).
   • إذا كان $y \in H^1(\Gamma)$، فإن الخطأ هو $O(h)$ (مع $\varrho = h$).
   • إذا كان $y \in H^2(\Gamma)$ (هدف سلس)، فإن الخطأ هو $O(h^2)$ (مع $\varrho = h^2$).
2. بناء محلل أمثل: اشتقاق محلل سريع يعتمد على متمم شات مع محسن مسبق للكتلة المجمعة. يظهر التحليل النظري أن عدد الشرط (condition number) محدود بشكل مستقل عن حجم الشبكة $h$ ومعلمة التنظيم $\varrho$ (بشرط $\varrho \leq h$).
3. التعامل مع شروط نيومان: بناء فضاء عناصر محدودة متوافق مع دوال أساس تحقق صراحةً شروط نيومان المتجانسة، مما يبسط التنفيذ مقارنة بالطرق التي تتطلب مضاعفات لاغرانج.

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

يتحقق المؤلفون من نظريتهم من خلال تجارب عددية في 3D ($\Omega = (0,1)^3$) باستخدام ثلاث دوال أهداف مختلفة:

• الحالة 1: هدف سلس ($y \in H^2(\Gamma)$)

  • الهدف: $y = \cos(\pi x_1)\cos(\pi x_2)\cos(\pi x_3)$.
  • اختيار المعلمة: $\varrho = h^2$.
  • النتيجة: تحقيق تقارب من الدرجة الثانية ($O(h^2)$) في خطأ $L^2(\Gamma)$.
  • أداء المحلل: ظل عدد تكرارات PCG لنظام متمم شات مستقلاً عن مستوى الشبكة (حوالي 4–29 تكراراً اعتماداً على المحسن المسبق)، مما يؤكد قوة المحلل.

• الحالة 2: هدف أقل انتظاماً ($y \in H^{3/2-\epsilon}(\Gamma)$)

  • الهدف: كثير حدود من الدرجة الثانية لا يحقق شروط نيومان المتجانسة.
  • اختيار المعلمة: $\varrho = h^{3/2}$.
  • النتيجة: لوحظ معدل تقارب قدره $O(h^{1.5})$، وهو ما يطابق التوقعات النظرية.

• الحالة 3: هدف غير مستمر ($y \in L^2(\Gamma)$)

  • الهدف: دالة ثابتة قطعة قطعة (1 داخل منطقة فرعية، و0 خارجها).
  • اختيار المعلمة: $\varrho = h$.
  • النتيجة: لوحظ معدل تقارب قدره $O(h^{0.5})$، وهو ما يتوافق مع حالة الانتظام الأدنى.

في جميع الحالات، أظهرت محللات AMG-PCG و Schur Complement PCG أعداد تكرار مستقلة عن الشبكة، مما يثبت كفاءة الحلول السريعة المقترحة.

5. الأهمية والآفاق المستقبلية

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