Inverse problem for a multi-term time-fractional diffusion equation with the Caputo derivatives

إليك ملخص تقني مفصل للورقة البحثية بعنوان "المسألة العكسية لمعادلة الانتشار ذات الكسور الزمنية متعددة الحدود باستخدام مشتقات كابوتو" (Inverse Problem for a Multi-Term Time-Fractional Diffusion Equation with the Caputo Derivatives) من إعداد راوشان أشوروف ودامير شاموراتوف.

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

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

$$\sum_{j=1}^{M} q_j \partial_t^{\rho_j} u(x, t) + A(x, D)u(x, t) = f(x)g(t), \quad x \in \Omega, \, 0 < t \leq T$$

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

• المؤثرات: $\partial_t^{\rho_j}$ تمثل مشتقة كابوتو الكسرية من الرتبة $0 < \rho_M < \dots < \rho_1 < 1$. و $A(x, D)$ هو مؤثر تفاضلي متجانس، موجب، وإهليلجي صوري من الرتبة الزوجية $m \geq 2$.
• حد المصدر: المصدر قابل للفصل، $f(x)g(t)$، حيث الجزء المعتمد على الزمن $g(t)$ معلوم، والمكون المكاني $f(x)$ هو المجهول المراد إعادة بنائه.
• الشروط:
  • الشرط الابتدائي: $u(x, 0) = \phi(x)$.
  • الشروط الحدية: شروط متجانسة $B_j u(x, t) = 0$ على $\partial \Omega$.
  • التحديد الزائد: يُعطى شرط إضافي عند زمن داخلي ثابت $t_0 \in (0, T]$: $u(x, t_0) = \Psi(x)$.

الهدف هو تحديد الزوج $(u(x, t), f(x))$ الذي يحقق هذه الشروط. المسألة هي بطبيعتها سيئة التحديد (ill-posed) بسبب احتمال تلاشي المقام في صيغة إعادة البناء.

2. المنهجية

استخدم المؤلفون طريقة فورييه (التحلل الطيفي) مقترنة بالتحليل التقاربي المتقدم للدوال الخاصة.

1. التحلل الطيفي:

   • تم اختزال المسألة إلى سلسلة من المعادلات التفاضلية العادية عن طريق توسيع الحل $u(x, t)$ والمصدر $f(x)$ بدلالة الدوال الذاتية $\{v_k(x)\}$ والقيم الذاتية $\{\lambda_k\}$ للمؤثر $A(x, D)$.
   • حول هذا الإجراء المعادلة التفاضلية الجزئية إلى نظام من المعادلات التفاضلية العادية الكسرية متعددة الحدود للمعاملات الفورية $T_k(t)$ و $f_k$.

2. حل المعادلات التفاضلية العادية الكسرية:

   • يتم التعبير عن حل المعادلات التفاضلية العادية الكسرية باستخدام دالة ميتغ-ليفلر متعددة الحدود (multinomial Mittag-Leffler function)، ويُرمز لها بـ $E_{\rho', \beta}(z_1, \dots, z_M)$.
   • تتضمن صيغة إعادة بناء المعامل المجهول $f_k$ حداً في المقام $b_{k, \rho_1}(t_0)$، وهو تكامل يتضمن $g(t)$ ودالة ميتغ-ليفلر.

3. التحليل التقاربي (المساهمة الجوهرية):

   • اشتق المؤلفون توسعات تقاربية دقيقة وجديدة لدالة ميتغ-ليفلر متعددة الحدود للقيم الكبيرة.
   • تُستخدم هذه التوسعات لإنشاء حدود دنيا منتظمة للحد في المقام $b_{k, \rho_1}(t_0)$. هذه الخطوة حاسمة لإثبات تقارب السلسلة ووجود الحل.

4. براهن التقارب:

   • باستخدام الحدود الدنيا المشتقة وخصائص دالة ميتغ-ليفلر، أثبت المؤلفون التقارب المطلق والمنتظم لسلسلة فورييه لكل من الحل $u(x, t)$ والمصدر $f(x)$ في فضاءات سوبوليف والفضاءات المتصلة المناسبة.

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

• التقديرات التقاربية الجديدة: تقدم الورقة أول توسع تقاربي دقيق لهذا الشكل المحدد من دالة ميتغ-ليفلر متعددة الحدود المطلوبة لمعادلات الانتشار الكسرية متعددة الحدود. وهذا يسمح بالتحكم الصارم في المقامات الطيفية.
• إطار عمل عام للمؤثرات: على عكس الأعمال السابقة التي اقتصرت على مؤثر لابلاس أو مؤثر لابلاس الكسري، تتعامل هذه الدراسة مع مؤثر إهليلجي ذاتي متجانس عام $A(x, D)$ من الرتبة $m \geq 2$ مع شروط حدية عامة.
• التعامل مع سوء التحديد: عالج المؤلفون بدقة الطبيعة سيئة التحديد للمسألة العكسية. فقد حددوا الشروط التي تكون فيها المسألة جيدة التحديد (حل فريد) وشخصوا حالة عدم الوحدانية (عندما يتلاشى المقام).
• نظريات الوجود والوحدانية:
  • الوجود: مثبت تحت فرضيات النعومة على البيانات الأولية ($\phi, \Psi$) والدالة الزمنية $g(t)$.
  • الوحدانية: مثبتة إذا كانت مجموعة المؤشرات التي يتلاشى عندها المقام فارغة. إذا كانت المجموعة غير فارغة، فقد استُخلصت شروط التوافق الضرورية والكافية، ووُصف هيكل عدم الوحدانية بالكامل (توجد حلول لانهائية).

4. النتائج الرئيسية

أ. السلوك التقاربي (النظرية 3.1)

لقيم $|z_1|$ الكبيرة، تسلك دالة ميتغ-ليفلر متعددة الحدود كما يلي:
$$E_{\rho', \beta}(z_1, \dots, z_M) = -\frac{1}{z_1 \Gamma(\beta - \rho_1)} - \frac{1}{z_1^2} \left( \frac{1}{\Gamma(\beta - 2\rho_1)} - \sum \dots \right) + O(|z_1|^{-p-1})$$
هذا التوسع هو الأساس لتقديرات الحد الأدنى.

ب. الحدود الدنيا للمقام (القسم 4)

الحد $b_{k, \rho_1}(t_0)$ هو الذي يحدد قابلية الحل.

• الحالة 1 (دالة $g(t)$ حافظة للإشارة): إذا لم يتغير إشارة $g(t)$، فإن $|b_{k, \rho_1}(t_0)| \geq \frac{C}{\lambda_k}$. في هذه الحالة، لا يتلاشى المقام أبداً، مما يضمن حلاً فريداً.
• الحالة 2 (دالة $g(t)$ متغيرة الإشارة): إذا كان $g(0) \neq 0$ وتغيرت إشارة $g(t)$، فقد يتلاشى المقام لعدد محدود من المؤشرات $k$.
  • إذا كانت مجموعة المؤشرات التي يتلاشى عندها المقام $K_{0, \rho_1}$ فارغة، فإن الحل فريد.
  • إذا كانت $K_{0, \rho_1}$ غير فارغة، يوجد حل إذا وفقط إذا حققت البيانات شروط التعامد: $\Psi_k = \phi_k [1 - \lambda_k t_0^{\rho_1} E_{\dots}(-\lambda_k t_0^{\rho_1})]$ لجميع $k \in K_{0, \rho_1}$. في هذه الحالة، الحل ليس فريداً (يمكن إضافة معاملات اختيارية لـ $k \in K_{0, \rho_1}$).

ج. الوجود والوحدانية (النظريات 5.1، 5.2، 6.1)

• الوجود: يوجد حل كلاسيكي $(u, f)$ إذا كانت البيانات $\phi, \Psi$ تنتمي إلى المجال $D(\hat{A}^{\tau+1})$ مع $\tau > N/2m$، وكانت $g(t)$ تحقق شروط الإشارة.
• الوحدانية: الحل فريد إذا كان $b_{k, \rho_1}(t_0) \neq 0$ لجميع $k$. إذا وجدت أصفار، تفشل الوحدانية ما لم تتحقق شروط التوافق، وحتى حينها، لا يكون الحل فريداً.

د. حالة خاصة: الطارة (القسم 7)

تم تطبيق النتائج على طارة ذات $N$ أبعاد مع مؤثرات ذات معاملات ثابتة. تم تحديد مجال القوة الكسرية للمؤثر صراحةً مع فضاءات سوبوليف $L^{\tau m}_2(T^N)$، مما يوفر إطاراً تحليلياً وظيفياً ملموساً.

5. الأهمية

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

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