A monotone iterative reconstruction method for an inverse drift problem in a two-dimensional parabolic equation

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

تتناول الورقة البحثية مشكلة الانجراف العكسية (inverse drift problem) لمعادلة مكافئة (parabolic equation) ثنائية الأبعاد معرفة على المربع الوحدة $\Omega = (0,1) \times (0,1)$. المعادلة الحاكمة هي:
$$(\partial_t - \Delta + q(x)\partial_{x_1} + C_p)u(x, t) = f(x)$$
تخضع لشروط حدودية مختلطة (ديريكليه في الأعلى/الأسفل، ونيومان في اليسار/اليمين) وشرط ابتدائي $u(x,0) = u_0(x)$.

• المجهول: معامل الانجراف المعتمد على المكان $q(x)$ (تحديداً معامل الحد $\partial_{x_1}$).
• البيانات: الملاحظة النهائية (بيانات الزمن النهائي) $g(x) = u(x, T)$.
• الهدف: استعادة معامل الانجراف المجهول $q(x)$ باستخدام الحالة النهائية $g(x)$ فقط، مع العلم بأن المصدر $f$ والبيانات الابتدية $u_0$ والثابت المحتمل $C_p$ معروفة.

هذه المشكلة غير مستقرة (ill-posed) وغير خطية بطبيعتها، حيث يعتمد الحل بشكل غير خطي على المعامل، ويمكن أن تؤدي الاضطرابات الصغيرة في البيانات إلى أخطاء كبيرة في إعادة البناء.

2. المنهجية

يقترح المؤلفون طريقة إعادة بناء تكرارية رتيبة (monotone iterative reconstruction method) تعتمد على نظرية النقطة الثابتة والمبادئ القصوى. تتضمن المنهجية الجوهرية ثلاث خطوات رئيسية:

أ. بناء المؤثر $K$

يعرّف المؤلفون مؤثراً غير خطي $K$ يعمل على فضاء دوال الانجراف المقبولة. بالنسبة لانجراف تجريبي $\psi$، يُعرّف المؤثر كالتالي:
$$K\psi = \frac{f(x) - \partial_t u(x, T; \psi) + \Delta g(x) - C_p g(x)}{\partial_{x_1} g(x)}$$
حيث $u(x, T; \psi)$ هو حل المسألة الأمامية مع الانجراف $\psi$.

• خاصية النقطة الثابتة: تكون الدالة $q$ هي الانجراف الحقيقي إذا وفقط إذا كان $Kq = q$.
• صلاحية التعريف: تم إثبات أن المقام $\partial_{x_1} g(x)$ موجب تماماً بموجب افتراضات محددة على البيانات (الافتراض 2.2)، مما يضمن أن المؤثر معرف بشكل جيد.

ب. تحليل الرتابة

المساهمة النظرية المركزية هي إثبات أن المؤثر $K$ هو رتيب (monotone). وتحديداً، إذا كان $\psi_1 \leq \psi_2$، فإن $K\psi_1 \leq K\psi_2$.

• تم اشتقاق هذه الخاصية باستخدام المبادئ القصوى للمعادلات المكافئة ولمّة هوف (Hopf's lemma).
• يعتمد الإثبات على إثبات موجبية المشتقات المكانية ($\partial_{x_1} u$) والمشتقات المختلطة ($\partial_{x_1} \partial_t u$) للحل في ظل الشروط الحدودية والابتدية المعطاة.

ج. الخوارزمية التكرارية

بناءً على الرتابة، صاغ المؤلفون مخططاً تكرارياً متقارباً:

1. التهيئة: البدء بتقدير حد علوي $q_0 = \frac{f + \Delta g - C_p g}{\partial_{x_1} g}$.
2. التكرار: حساب $q_{n+1} = K q_n$.
3. التقارب: ثبت أن المتتالية $\{q_n\}$ تتناقص رتيباً وتتقارب نقطياً إلى الحل الفريد $q$.

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

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

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

أجرى المؤلفون تجارب عددية باستخدام طرق الفروق المحدودة (طريقة أويلر الخلفية في الزمن، والفروق المركزية في المكان) للتحقق من النظرية.

• حالات الاختبار:
  • معامل سلس: $q(x,y) = 1 + \sin(\pi x)\sin(\pi y)$.
  • ثابت مجزأ: تضمين مربع ذو قفزة حادة.
  • هندسة معقدة: منطقة على شكل "رمز صيني" تحتوي على زوايا وثقوب.
• الأداء:
  • التقارب: يتناقص خطأ $L^2$ النسبي بسرعة في التكرارات القليلة الأولى للبيانات الخالية من الضجيج.
  • الرتابة: تُظهر التكرارات تحسناً رتيباً، حيث يصبح الملف المستعاد أكثر حدة وأقرب إلى الحل الدقيق خطوة بخطوة.
  • المتانة تجاه الضجيج: تظل الطريقة فعالة تحت ظروف بيانات نهائية مشوبة بالضجيج (حتى 3%). تم استخدام استراتيجية إزالة الضجيج (تنعيم البيانات النهائية) للتعامل مع التفاضل المطلوب في المؤثر $K$.
  • الانقطاعات: نجحت الطريقة في استعادة المعاملات المنقطعة (الثابتة مجزأة والأشكال المعقدة) دون تذبذب كبير، مما حافظ على الواجهات الحادة.

5. الأهمية

تقدم هذه الورقة مساهمة هامة في مجال المسائل العكسية للمعادلات المكافئة من خلال:

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

باختوصر، نجحت الورقة في جسر الفجوة بين الوحدانية النظرية وإعادة البناء العملية لمسائل الانجراف العكسية في بعدين، مقدمة خوارزمية تكرارية رتيبة ومتينة تعمل بشكل جيد حتى مع البيانات المشوبة بالضجيج.