Solution of variable order fractional differential equations using Homotopy Analysis Method

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

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

تتناول الورقة تحدي حل المعادلات التفاضلية الكسرية ذات الرتب المتغيرة (VOFDEs)، مع التركيز بشكل خاص على معادلات الانتشار ذات الرتب المتغيرة.

• السياق: بينما يُستخدم الحساب الكسري ذو الرتب الثابتة على نطاق واسع لنمذجة الانتشار في الأوساط المتجانسة، فإنه يفشل في وصف العمليات في الأوساط غير المتجانسة (مثل الأوساط المسامية) حيث تتغير خصائص الانتشار مع المكان، أو الزمان، أو كلاهما.
• الفجوة: في مثل هذه الأنظمة المعقدة، لا تكون رتبة المشتقة الكسرية، والتي يُرمز لها بـ $\alpha$، ثابتة بل تتغير كدالة للمكان ($x$)، أو الزمان ($t$)، أو كليهما ($\alpha(x,t)$).
• الهدف: يهدف المؤلفان إلى حل نوعين محددين من معادلات الانتشار ذات الرتب المتغيرة:
  1. معادلة انتشار زمنية كسرية خطية ذات رتبة متغيرة (بدون حد مصدر).
  2. معادلة انتشار غير خطية ذات رتبة متغيرة تتضمن حداً تفاعلياً.
• الحداثة: يزعم المؤلفان أنه بينما استُخدمت الطرق العددية (مثل الفروق المحدودة) سابقاً، لم يحاول أي باحث حل هذه المعادلات المحددة ذات الرتب المتغيرة باستخدام طريقة تحليل الهوموتوبي (HAM) قبل هذا العمل.

2. المنهجية: طريقة تحليل الهומوتوبي (HAM)

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

الخطوات الرئيسية للمنهجية:

1. التعريفات: تستخدم الدراسة مشتق كابوتو الكسري ذو الرتبة المتغيرة وتكامل ريمان-ليوفيل ذو الرتبة المتغيرة.
2. بناء المؤثر:
   • يتم اختيار مؤثر خطي ($L$)، يعتمد عادةً على أعلى حد للمشتقة الكسرية.
   • يتم بناء مؤثر غير خطي ($N$) لتمثيل المعادلة التفاضلية بأكملها، بما في ذلك حدود المصدر واللاخطية.
3. معادلة تشوه الرتبة صفر: يتم إنشاء هوموتوبي لتحويل التخمين الأولي إلى الحل الدقيق.
4. معادلات التشوه من الرتبة $m$: يتم تفكيك المشكلة إلى سلسلة من المشكلات الخطية الفرعية. تُشتق معادلة التشوه من الرتبة $m$ كالتالي:
   $$L[u_m(x,t) - \chi_m u_{m-1}(x,t)] = \hbar R_m(\vec{u}_{m-1})$$
   حيث:
   • $\chi_m$ هو معامل (يساوي 0 إذا كان $m \leq 1$، و1 إذا كان $m > 1$).
   • $\hbar$ هو معامل التحكم في التقارب.
   • $R_m$ هو خطأ البواقي من الرتبة $(m-1)$.
5. الحل المتسلسلة: يُعبر عن الحل كسلسلة لانهائية $u(x,t) = \sum_{n=0}^{\infty} u_n(x,t)$.
6. التحكم في التقارب: الخطوة الحاسمة تتضمن تحديد القيمة المثلى لـ معامل التحكم في التقارب ($\hbar$). يقوم المؤلفون بتقليل متوسط مربع خطأ البواقي ($E_m$) لضمان تقارب سلسلة الحل بسرعة ودقة.

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

• أول تطبيق لـ HAM على VOFDEs: تؤسس الورقة لأول تطبيق ناجح لـ HAM لحل معادلات الانتشار الكسرية ذات الرتب المتغيرة، مما يوسع نطاق تطبيق HAM إلى ما وراء المشكلات ذات الرتب الثابتة.
• التعامل مع الرتب المتغيرة: تتعامل الطريقة بفعالية مع الحالات التي تعتمد فيها الرتبة $\alpha$ على الزمان ($t$)، أو المكان ($x$)، أو مزيج منهما ($\alpha = 0.8 + 0.2xt$).
• الديناميكيات غير الخطية: نجح المؤلفان في حل معادلة انتشار غير خطية ذات رتبة متغيرة تتضمن حداً تفاعلياً، مع تقديم تحليل الاستقرار لحالات مختلفة، وهو ما وُصف بأنه الأول من نوعه في سياق الأنظمة الكسرية ذات الرتب المتغيرة.
• الحل التحليلي المتسلسل: على عكس الطرق العددية البحتة التي توفر نقاطاً منفصلة، توفر HAM حلاً متسلسلاً يقدم رؤية أعمق للسلوك الوظيفي للنظام.

4. النتائج والمناقشة

تحقق المؤلفان من صحة منهجيتهما من خلال دراستي حالة:

الحالة 1: الانتشار الخطي ذو الرتبة المتغيرة

• المعادلة: $D_t^{\alpha(x,t)} u = K \frac{\partial^2 u}{\partial x^2}$ مع $\alpha(x,t) = 0.8 + 0.2xt$.
• التحقق: تمت مقارنة النتائج مع الحلول العددية الموجودة من قبل Sun et al. [25].
• التقارب: من خلال تقليل خطأ البواقي، وُجد أن القيمة المثلى لـ $\hbar$ هي -0.975296.
• الدقة: مع وجود خمسة حدود فقط في السلسلة، انخفض خطأ البواقي إلى $2.18 \times 10^{-16}$، مما أظهر دقة عالية للغاية. تطابقت منحنيات التطور تماماً مع الأدبيات الموجودة.

الحالة 2: الانتشار غير الخطي ذو الرتبة المتغيرة

• المعادلة: معادلة غير خطية تتضمن الانتشار وحداً تفاعلياً ($u(1-u)$).
• التقارب: تم تحديد القيمة المثلى لـ $\hbar$ لتكون -0.134256.
• الأداء: انخفض خطأ البواقي بشكل ملحوظ مع زيادة عدد الحدود (من 0.289 عند حدين إلى 0.0056 عند 4 حدود).
• التصور: قدمت الدراسة رسوماً بيانية ثلاثية الأبعاد توضح تطور $u(x,t)$ عبر المكان والزمان، مما أكد الاستقرار والأهمية الفيزيائية للحل.

5. الأهمية والخاتمة

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

باختصار، نجحت الورقة في جسر الفجوة بين نظرية الحساب الكسري المتقدمة وتقنيات الحل التحليلية العملية، مثبتة أن HAM طريقة قوية لمعالجة تعقيدات الانتشار الكسري ذي الرتب المتغيرة.