Cavity shape reconstruction with a homogeneous Robin condition via a constrained coupled complex boundary method with ADMM

ملخص تقني: إعادة بناء شكل التجويف باستخدام شرط روبن المتجانس عبر طريقة الحدود المعقدة المترابطة (CCBM) وخوارزمية (ADMM)

بيان المشكلة
تتناول هذه الدراسة مشكلة تحديد الحدود العكسية لإعادة بناء حدود داخلية مجهولة $\Gamma$ لتجويف داخل جسم موصل $\Omega$. النطاق متصل مزدوجاً، يحده حد خارجي يمكن الوصول إليه $\Sigma$ وحد داخلي لا يمكن الوصول إليه $\Gamma$. المتغير الحالة $u$ (الذي يمثل الجهد الكهربائي أو درجة الحرارة) يحقق معادلة لابلاس $-\Delta u = 0$ في $\Omega$. على الحد الخارجي المتاح $\Sigma$، يتم تحديد بيانات كوشي $(f, g)$، حيث $u=f$ و $\partial_n u = g$. وعلى الحد المجهول $\Gamma$، تحقق الدالة شرط روبن المتجانس: $\partial_n u + \alpha u = 0$، حيث $\alpha$ ثابت موجب تماماً.

هذه المشكلة غير محددة جيداً (ill-posed) بطبيعتها. وبينما قد يتحقق التفرد المحلي لنطاقات معينة، فإن زوجاً واحداً من بيانات كوشي يؤدي عموماً إلى عدد لا نهائي من النطاقات الممكنة، مما يجعل إعادة بناء $\Gamma$ غير محددة دون قيود إضافية أو تنظيم (regularization). علاوة على ذلك، غالباً ما تواجه طرق تحسين الشكل الكلاسيكية صعوبة في استعادة الميزات المقعرة أو غير المحدبة للحدود، لا سيما في حالات القياس الواحد.

المنهجية
يقترح المؤلفون تطبيقاً مبتكرًا لـ طريقة الحدود المعقدة المترابطة (CCBM) بالتزامن مع إطار عمل خوارزمية اتجاه المتغيرات المضافة (ADMM) لمعالجة هذه التحديات.

1. صياغة CCBM:
   تمت إعادة صياغة مشكلة القيم الحدية الحقيقية المفرطة في التحديد كمسألة قيم حدية معقدة. يتم ربط بيانات كوشي على $\Sigma$ عبر شرط روبن معقد: $\partial_n u + i\rho u = g + i\rho f$ على $\Sigma$، بينما يظل شرط روبن الفيزيائي $\partial_n u + \alpha u = 0$ على $\Gamma$. هنا، $u = u_1 + iu_2$ هو حل ذو قيم معقدة.
   تتحول المشكلة العكسية إلى مشكلة تحسين شكل: البحث عن النطاق $\Omega$ بحيث يتلاشى الجزء التخيلي من الحل، $u_2$، في $\Omega$. دالة التكلفة المراد تقليلها هي $J(\Omega) = \frac{1}{2} \int_\Omega |u_2|^2 dx$.

2. تحليل حساسية الشكل:
   يشتق البحث مشتق الشكل لدالة التكلفة باستخدام حساب الشكل. يتم التعبير عن التدرج كتكامل حدودي على $\Gamma$ يتضمن حل الحالة $u$ وحالة مُلحقة (adjoint state) $p$. النظام الملحق هو معادلة تفاضلية جزئية معقدة مع شروط حدية محددة مستمدة من لاغرانج للمسألة التحسينية.

3. التحسين المقيد عبر ADMM:
   لتحسين الاستقرار ضد الضوضاء والبدايات السيئة، خاصة بالنسبة للهندسات المقعرة، أدخل المؤلفون قيد عدم تساوٍ على الجزء الحقيقي من متغير الحالة: $a \leq \Re\{u\} \leq b$ في كل مكان في $\Omega$.
   يؤدي هذا إلى مسألة تحسين شكل مقيدة (المسألة 3). ولحل هذه المسألة، استخدم المؤلفون استراتيجية ADMM:

• فصل المتغيرات: يتم إدخال متغير مساعد $v$ لفك الارتباط بين قيد عدم التساوي ومعادلة الحالة.
• لاغرانج المضاف: يتم حل المسألة عن طريق تقليل دالة لاغرانج المضافة.
• المخطط التكراري: تقوم الخوارزمية بالتبادل بين:
  • تحديث الشكل (SP1): تقليل لاغرانج المضاف بالنسبة للشكل $\omega$ (باستخدام تدرج الشكل وطريقة تدرج Sobolev لتشويه الشبكة).
  • تحديث المتغير (SP2): تحديث المتغير المساعد $v$ عبر الإسقاط على المجموعة المحدبة المسموح بها $K = \{v \in L^2(\Omega) : a \leq v \leq b\}$.
  • تحديث المضاعف: تحديث مضاعف لاغرانج $\lambda$.

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

• تطبيق مبتكر لـ CCBM: يوسع هذا العمل إطار عمل CCBM، الذي استُخدم سابقاً لمشكلات المصدر والموصلية العكسية، ليشمل تحديد حدود روبن.
• الصياغة المقيدة: يعد دمج قيود عدم التساوي في إطار عمل CCBM عبر ADMM مساهمة أساسية. هذا النهج يفرض حدوداً فيزيائية على متغير الحالة دون إدخال مصطلحات تنظيم هندسي إضافية (مثل عقوبة المحيط)، مما يحافظ على بنية CCBM الأصلية مع تعزيز الاستقرار.
• اشتقاق مشتق الشكل: يقدم البحث اشتقاقاً دقيقاً لتدرج الشكل لنظام الحالة المعقد المرتبط بفك ADMM، بما في ذلك تعبيرات بديلة متعددة للتدرج بناءً على صياغات مختلفة للملحق (adjoint).
• التنفيذ العددي: تم تنفيذ الطريقة باستخدام طريقة العناصر المحدودة (FEM) في برنامج FreeFem++، باستخدام تمثيل $H^1$-Riesz لتدرج الشكل لضمان تشوهات سلسة للشبكة.

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

• كشف التقعر: أظهرت طريقة CCBM غير المقيدة أداءً أفضل مقارنة بالمناهج الكلاسيكية في اكتشاف المناطق المقعرة، لكنها لا تزال تظهر قيوداً في تحديد عمق التقعرات، خاصة مع معاملات روبن الصغيرة ($\alpha$) أو البيانات المشوبة بالضوضاء.
• تأثير القيود: أظهر نهج ADMM المقيد استقراراً فائقاً. فقد نجح في كبح التذبذبات الزائفة الناتجة عن الضوضاء ووفر عمليات إعادة بناء أكثر استقراراً عبر مختلف البدايات الأولية.
• حساسية المعاملات: وُجد أن معامل الوزن $\rho$ في الربط المعقد يؤثر على الحساسية للميزات الهندسية؛ حيث تعمل القيم الأكبر على تحسين دقة التفاصيل الدقيقة ولكنها تتطلب ضبطاً دقيقاً لتجنب تضخيم الضوضاء.
• تأثير البيانات: اعتمدت جودة إعادة البناء بشكل كبير على بيانات ديريكليه $f$ المحددة. وقد أعطت الإثارات المتغيرة مكانياً (مثل الدوال المثلثية) استعادة أفضل بكثير للميزات المقعرة مقارنة بالبيانات الثابتة.
• المتانة ضد الضوضاء: في مستويات الضوضاء المنخفضة ($\delta = 0.005, 0.01$)، حافظت طريقة ADMM المقيدة على الاتساق الهندسي، بينما أظهرت الطريقة غير المقيدة حساسية أكبر للبداية الأولية والضوضاء.

الأهمية والادعاءات
يزعم البحث أن إطار عمل ADMM-CCBM المقترح يوفر نهجاً أكثر استقراراً ومتانة لمشكلات إعادة بناء الشكل العكسية التي تتضمن شروط روبن المتجانسة مقارنة بصياغة CCBM القياسية. وتكمن الأهمية في القدرة على التعامل مع عدم التحديد الشديد لمشكلات القياس الواحد، لا سيما في استعادة الهندسة غير المحدبة.

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