A backward problem for the time-fractional pseudo-parabolic equation with a variable coefficient

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

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

تتناول الورقة مسألة عكسية (تحديداً مسألة خلفية في الزمن) لـ معادلة شبه بارابولية كسرية زمنية ذات معامل يعتمد على الزمن. المعادلة الحاكمة هي:

$$\partial^\alpha_t u(x, t) - u_{xx}(x, t) - \mu(t) u_{xxt}(x, t) = f(x, t), \quad (x, t) \in (0, l) \times (0, T)$$

المكونات الرئيسية:

• المؤثر: $\partial^\alpha_t$ هو مشتق كابوتو (Caputo) الكسري من الرتبة $\alpha \in (0, 1)$.
• المعامل: $\mu(t)$ هو معامل شبه بارابولي متغير زمنياً.
• الشروط الحدية: شروط ديريكليه المتجانسة $u(0, t) = u(l, t) = 0$.
  ها.
• البيانات: حد المصدر $f(x, t)$ وقياس الزمن النهائي $u(x, T) = \psi(x)$ معلومان.
• الهدف: إعادة بناء الحالة الأولية المجهولة $u_0(x) = u(x, 0)$ والمسار الكامل $u(x, t)$ على الفترة $[0, T]$.

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

2. المنهجية

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

أ. التحليل النظري (الوجود والوحدانية)

• التفكك الطيفي: يتم تحليل المسألة باستخدام الدوال الذاتية للمؤثر المكاني $-\partial_{xx}$ (توسع متسلسلة الجيب).
• اختزال فولتير: من خلال إسقاط المعادلة التفاضلية الجزئية على القاعدة الذاتية، تتحول المسألة إلى نظام من المعادلات التفاضلية العادية (ODEs) الكسرية. يتم تحويل هذه المعادلات إلى معادلات فولتير التكاملية من النوع الثاني ذات نوى ضعيفة التفرد.
• إعادة البناء الصريح: يستنتج المؤلفون تمثيلاً متسلسلياً صريحاً للمحلول $u(x, t)$ والبيانات الأولية $u_0(x)$.
• لمة الموجبية: تتضمن خطوة حاسمة إثبات أن دالة المعامل $A_k(T)$ في صيغة إعادة البناء موجبة تماماً، مما يضمن أن الحل محدد جيداً.
• الاستقرار: تم إثبات الاعتماد المستمر للحالة الأولية على البيانات النهائية، مما يوفر تقدير استقرار في معايير $L^2$.

ب. التقطيع العددي (المحلل الأمامي)

لحل المسألة المباشرة (الأمامية) عددياً، تم تطوير مخطط فرق متناهي (finite-difference scheme):

• التقطيع المكاني: تُستخدم الفروق المركزية القياسية من الدرجة الثانية لـ Laplacian ($u_{xx}$) وحد المشتق المختلط ($u_{xxt}$).
• التقديد الزمني: يتم تقريب مشتق كابوتو باستخدام مخطط L1 على شبكة زمنية متدرجة ($t_k = T(k/M)^r$). الشبكات المتدرجة ضرورية للتعامل مع التفرد الأولي النموذجي للمشتقات الكسرية ($\partial_t u \sim t^{\alpha-1}$).
• النظام الخطي: ينتج عن التقطيع نظام خطي ثلاثي الأقطار عند كل خطوة زمنية، يتم حله بكفاءة عبر خوارزمية توماس (Thomas algorithm).
• التحليل: ثبت أن المخطط مستقر غير مشروط ومتقارب. معدل التقارب مثالي للشبكة المتدرجة، وهو محدود أساساً بالدقة من الدرجة الأولى لتقطيع حد شبه البارابولي.

ج. خوارزمية إعادة البناء العكسية

تُحل المسألة العكسية باستخدام نهج بناء المؤثر + تنظيم تيكونو (Tikhonov regularization):

1. بناء الخريطة الأمامية: يتم بناء مؤثر الحل الأمامي المنفصل $F_h$ عن طريق تشغيل المحلل الأمامي $N-1$ من المرات مع شروط أولية نبضية (متجهات القاعدة الأساسية) وحد مصدر صفري.
2. التقسيم الأفيني: يتم تفكيك الحالة النهائية إلى جزء متجانس (يعتمد على $u_0$) وجزء مدفوع (يعتمد على $f$).
3. التنظيم: تُصاغ إعادة البناء كمسألة تقليل باستخدام تنظيم تيكونو من الدرجة صفر:
   $$J_\lambda(u_{0,h}) = \frac{1}{2} \| F_h u_{0,h} - d_h \|_2^2 + \frac{\lambda}{2} \| u_{0,h} \|_2^2$$
   حيث $d_h$ هو الفرق بين البيانات النهائية المقاسة والمساهمة المدفوعة. يعمل المعامل $\lambda$ على تثبيت النظام غير جيد التشكيل.

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

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

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

• التحقق الخالي من الضجيج: نجحت الطريقة في استعادة الحالة الأولية بدقة عالية. يظهر تحليل الخطأ معدلات التقارب المتوقعة مع تحسين الشبكات المكانية ($N$) والزمنية ($M$). تتطابق الحالة الأولية المعاد بناؤها والحالة النهائية المنتشرة للأمام مع الحلول الدقيقة بشكل وثيق.
• القوة تجاه البيانات المشوبة بالضجيج:
  • أُجريت التجارب مع ضجيج غاوسي مضاف بمستويات 1%، 3%، و5%.
  • النتائج: مع زيادة مستويات الضجيج، زاد الخطأ في الحالة الأولية المعاد بناؤها، لكن النمو كان محكوماً ومتناسباً (اعتماد مستقر)، بدلاً من الانفجار.
  • التمثيل البصري: حتى مع ضجيج بنسبة 5%، حافظ تنظيم تيكونو على الشكل الرئيسي لملف الحالة الأولية، مما منع حدوث تذبذبات زائفة.
  • الحساسية للرتبة الكسرية: وُجد أن الدقة غير حساسة نسبياً للرتبة الكسرية $\alpha$، حيث كان مستوى الضجيج هو العامل المهيمن في التأثير على مقدار الخطأ.

5. الأهمية

يسد هذا العمل فجوة كبيرة في الأدبيات المتعلقة بـ المسائل العكسية للمعادلات شبه الباربولية الكسرية الزمنية.

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

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