Explicit Discrete Solution for Some Optimization Problems and Estimations with Respect to the Exact Solution

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

تتناول الورقة البحثية مسائل التوصيل الحراري في الحالة المستقرة التي تحكمها معادلة بواسون في مجال متعدد الأبعاد $\Omega \subset \mathbb{R}^n$ (تم تحليلها تحديداً في مجال مستطيل). تركز الدراسة على نظامين متميزين:

• النظام (S): يتميز بشروط حدودية مختلطة: شرط ديريكليه (درجة حرارة ثابتة $b$) على $\Gamma_1$، وشرط نيومان (تدفق حراري ثابت $q$) على $\Gamma_2$، وشرط أديباتي (تدفق صفري) على $\Gamma_3$.
• النظام ($S_\alpha$): مشابه للنظام (S)، ولكن تم استبدال شرط ديريكليه على $\Gamma_1$ بـ شرط روبن الحدودي (تدفق حراري بالحمل الحراري) يتضمن معامل $\alpha > 0$.

مرتبطة بهذه الأنظمة ثلاثة مسائل تحكم أمثل لكل نظام ($P_i$ و $P_{i\alpha}$ لـ $i=1,2,3$)، حيث تكون متغيرات التحكم هي:

1. $g$: مصدر طاقة داخلي موزع.
2. $q$: التدفق الحراري الحدودي على $\Gamma_2$.
3. $b$: درجة الحرارة البيئية على $\Gamma_1$.

الهدف هو تقليل دالة التكلفة التربيعية $J$ التي توازن بين انحراف الحالة (درجة الحرارة) عن هدف مرغوب $z_d$ وبين مقدار متغير التحكم (التنظيم).

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

2. المنهجية

يستخدم المؤلفون طريقة الفروق المحدودة (FDM) لتقطيع الأنظمة المستمرة ومسائل التحسين.

• التقطيع (Discretization):

  • تم تقطيع المجال بحجم شبكة منتظم $h = x_0/n$.
  • العقد الداخلية: تم استخدام تقريب مركزي من الدرجة الثانية كلاسيكي للمشتقة الثانية ($\partial^2 u / \partial x^2$).
  • الشروط الحدودية (القياسية):
    • تم تقريب شروط نيومان باستخدام فرق خلفي من الدرجة الأولى.
    • تم تقريب شروط روبن باستخدام فرق أمامي من الدرجة الأولى.
  • يؤدي هذا إلى أنظمة خطية ($S_h$ و $S_{h\alpha}$) ذات مصفوفات معاملات ثلاثية الأقطار.

• التحسين (Optimization):

  • تم تقطيع دالات التكلفة المستمرة لتكوين دالات تكلفة منفصلة ($J^h_i$ و $J^h_{i\alpha}$).
  • بما أن متغيرات التحكم يُفترض أنها ثابتة، فإن مسائل التحسين المنفصلة تختزل إلى تقليل دالات تربيعية لمتغير واحد، مما يسمح باشتقاق صيغ صريحة للتحكمات المنفصلة المثلى ($g^h_{op}, q^h_{op}, b^h_{op}$).

• تحليل التقارب:

  • تشتق الورقة تقديرات خطأ صريحة تقارن بين الحلول المنفصلة ($u^h, g^h_{op}$, إلخ) والحلول المستمرة الدقيقة ($u, g_{op}$, إلخ).
  • كما تحلل سلوك التقارب المزدوج: عندما يقترب حجم الشبكة $h \to 0$ ويقترب المعامل الحملي $\alpha \to \infty$.

• التقريب المحسن (القسم 7):

  • لمعالجة انخفاض رتبة التقارب ($O(h)$) الناتج عن تقريبات الحدود من الدرجة الأولى، يقترح المؤلفون مخطط فروق محدودة ثلاثي النقاط لشروط نيومان وروبن الحدودية.
  • يضيف هذا التعديل صياغة "النقطة الشبح" (ghost point)، مما يزيد فعلياً من دقة معالجة الحدود.

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

1. الحلول المنفصلة الصريحة: توفر الورقة تعبيرات تحليلية مغلقة لمتغيرات الحالة المنفصلة ومتغيرات التحكم المثلى لجميع أنواع مسائل التحسين الثلاثة (المصدر الموزع، التدفق، ودرجة الحرارة). وتعد هذه مساهمة كبيرة لأنها تتجنب الحلول العددية التكرارية لخطوة التحسين.
2. تقدير الخطأ الصارم: يثبت المؤلفون أن الحلول المنفصلة تتقارب مع الحلول المستمرة بدقة من الدرجة الأولى ($O(h)$) في معيار $L^2$ عند استخدام التقريبات الحدودية القياسية. كما يوفرون ثوابت صريحة لهذه الحدود الخطئية.
3. تحليل التقارب المزدوج: تؤسس الدراسة لتبادلية النهايات: الحل المنفصل لمسألة روبن يتقارب مع الحل المنفصل لمسألة ديريكليه مع $\alpha \to \infty$، وكلاهما يتقارب مع الحل المستمر مع $h \to 0$.
4. تحسين الرتبة: من خلال استخدام التقريب ثلاثي النقاط للشروط الحدودية، يوضح المؤلفون أن رتبة التقارب العالمية لمتغير الحالة تتحسن من $O(h)$ إلى $O(h^2)$، بينما يظل خطأ المشتقة عند $O(h)$.

4. النتائج

• معدلات التقارب:
  • طريقة FDM القياسية: الخطأ $\|u - u^h\|_H$ يتصرف كـ $C_1 h$.
  • طريقة FDM المحسنة (ثلاثية النقاط): الخطأ $\|u - \tilde{u}^h\|_H$ يتصرف كـ $D_1 h^2$.
• التحكمات المثلى: تظهر الصيغ الصريحة للتحكمات المثلى أن التحكمات المنفصلة المثلى ($g^h_{op}$, إلخ) تختلف عن التحكمات المستمرة بمقادๆ تتناسب مع $h$.
• المحاكاة العددية:
  • تؤكد المحاكاة على مجال مربع الوحدة ($x_0=y_0=1$) معدلات التقارب النظرية.
  • توضح الجداول والأشكال أن تقليص حجم الشبكة $h$ إلى النصف يقلل خطأ $L^2$ إلى النصف تقريباً، مما يؤكد سلوك $O(h)$.
  • تتحقق المحاكاة من أنه مع زيادة $\alpha$، تتقارب حلول نظام روبن ($S_\alpha$) مع حلول نظام ديريكليه ($S$).

5. الأهمية

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

باختصار، تجسر هذه الورقة الفجوة بين نظرية التحكم الأمثل المستمرة والتطبيق العددي المنفصل من خلال توفير صيغ منفصلة دقيقة وإثبات خصائص تقاربها، مع تقديم تقنية عملية لتعزيز الدقة العددية من خلال تحسين تقريبات الحدود.