Distributed optimal control problems governed by poroelasticity equations

ملخص تقني: مسائل التحكم الأمثل الموزعة المحكومة بمعادلات المسامية المرنة

بيان المشكلة
تتناول هذه الورقة البحثية التقريب العددي لمسائل التحكم الأمثل الموزعة المحكومة بنموذج "بيو" (Biot) للتماسك في المسامية المرنة (poroelasticity). يصف النظام الفيزيائي التفاعل بين تدفق السائل الخلالي وتشوه الوسط المسامي، وهو نموذج بالغ الأهمية في تطبيقات الجيوفيزياء، وميكانيكا الأنسجة البيولوجية، وعزل الكربون. تهدف مسألة التحكم المحددة إلى تقليل دالة التكلفة التي تقيس الانحراف في إزاحة الجسم الصلب ($u$) وضغط السائل ($p$) عن القيم المستهدفة المطلوبة ($u_C, p_C$)، مع فرض جزاء على مقدار متغير التحكم. تم تحديد التحكم كعنصر المصدر/المصب للسائل ($m$) الموزع داخل النطاق. تركز الدراسة على حالة التماسك شبه الساكن (quasi-static) حيث يتم إهمال تأثيرات التسارع، وتتناول تحديداً التحديات الرياضية والعددية المرتبطة بـ "انغلاق المسامية المرنة" (poroelastic locking)، والذي يحدث عندما يكون معامل التخزين النوعي ($s_0$) قريباً من الصفر، مما يجعل السائل والوسط المرن غير قابلين للانضغاط تقريباً.

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

1. الصياغة المستمرة وصلاحية الحل (Well-Posedness):

   • أعاد المؤلفون صياغة المعادلات الحاكمة عن طريق أخذ المشتق الزمني لتوازن الزخم، وتحويل النظام إلى مجموعة من المعادلات التفاضلية.
   • حددوا فضاءات الدوال المناسبة ($V$ للإزاحة، و $Q$ للضغط) وأثبتوا صلاحية الحل الضعيف باستخدام نظرية النقاط السرجية (saddle-point theory) وشرط لادزينكايانا-بابوشكا-برزي (LBB).
   • تم إثبات وجود ووحدانية التحكم الأمثل باستخدام الطريقة المباشرة في حساب التغيرات، بالاستفادة من تقعر دالة التكلفة وخطية مؤثر الحالة.
   • تم اشتقاق شروط الأمثلية الضرورية والكافية من الدرجة الأولى عبر نهج الحالة المرافقة (adjoint state)، مما أدى إلى متباينة تباينية للتحكم.

2. المخطط المنفصل (Discrete Scheme):

   • التقسيم المكاني (Spatial Discretization): تم تقسيم النطاق باستخدام عناصر محدودة ($V_h \times Q_h$) تستوفي شرط (inf-sup) (على سبيل المثال، عناصر MINI أو عناصر تايلور-هود).
   • التقسيم الزمني (Temporal Discretization): تم استخدام مخطط "باك وارد أويلر" (وتحديداً طريقة dG(0) في الزمن)، مما نتج عنه نظام منفصل تماماً.
   • تقسيم التحكم (Control Discretization): تستخدم الورقة نهج التقسيم المتغير (الذي اقترحه هينزي) لمتغير التحكم. وهذا يعني أن فضاء التحكم لا يتم تقسيمه صراحة؛ بل يتم الحصول على التحكم المنفصل عن طريق إسقاط الحالة المرافقة المنفصلة على مجموعة التحكم المسموح بها. هذا يتجنب أخطاء التقريب التي تنتج عادةً عن تقسيم فضاء التحكم.

3. تحليل الخطأ:

   • اشتق المؤلفون تقديرات خطأ مسبقة (a priori error estimates) للمخطط المنفصل تماماً.
   • يفصل التحليل إجمالي الخطأ إلى أخطاء التقسيم (للتحكم الثابت) وأخطاء الاضطراب الناتجة عن الفرق بين التحكم الأمثل المستمر والمنفصل.
   • تم وضع تقديرات الاستقرار لكل من الأنظمة الأمامية والحالة المرافقة الخلفية لتحديد هذه الأخطاء.

المساهمات الرئيسية
تقدم الورقة ثلاث مساهمات رئيسية في الأدبيات:

• صياغة متناظرة جديدة: تقدم صياغة متناظرة ثنائية الحقل (الإزاحة والضغط) مصممة خصيصاً لمسائل التحكم الأمثل ضمن نموذج "بيو"، وهي تختلف عن الصياغات الرباعية أو المختلطة السابقة.
• الوجود والوحدانية: أثبتت الورقة بصرامة صلاحية الصياغة الجديدة وأثبتت وجود ووحدانية التحكم الأمثل، جنباً إلى جنب مع اشتقاق شروط الأمثلية من الدرجة الأولى.
• تقديرات الخطأ المسبقة: توفر الورقة تحليلاً كاملاً للخطأ المسبق للمخطط المنفصل تماماً. ومن الجدير بالذكر أنها اشتققت تقديرات الخطأ للإزاحة والضغط في معايير $L^\infty(0, T; H^1(\Omega))$ و $L^2(0, T; L^2(\Omega))$، وللمتغير التحكم في معيار $L^2(0, T; L^2(\Omega))$، محققة معدلات تقارب مثالية تحت افتراضات الانتظام القياسية.

النتائج

• نظرياً: تثبت الورقة أن المخطط العددي يتقارب برتبة $O(h^s + \Delta t)$، حيث $h$ هو حجم الشبكة المكانية، و $\Delta t$ هي الخطوة الزمنية، و $s$ تتعلق بانتظام الحل. وتحديداً، إذا كان الحل منتظماً بما يكفي ($s=1$)، فإن التقارب يكون خطياً في كل من الفضاء والزمن ($O(h + \Delta t)$).
• عددياً: أكدت تجربة عددية على نطاق مربع الوحدة النتائج النظرية. باستخدام مكتبة FEniCS مفتوحة المصدر، حل المؤلفون نظام الأمثلية باستخدام خوارزمية أمامية-خلفية (forward-backward algorithm). وتظهر النتائج أن المخطط المقترح يحقق معدلات التقارب المثالية المتوقعة للإزاحة والضغط ومتغيرات التحكم عند استخدام عناصر MINI و Taylor-Hood من الدرجة الدنيا.

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