Non-uniform $\alpha$-Robust Alikhanov Mixed FEM with Optimal Convergence for the Time-Fractional Allen--Cahn Equation

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "طريقة العناصر المحدودة المختلطة (Mixed FEM) غير المنتظمة من نوع Alikhanov ذات المتانة $\alpha$ والتقارب الأمثل لمعادلة Allen–Cahn ذات الاشتقاق الزمني الكسري."

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

تتناول الورقة الحل العددي لـ معادلة Allen–Cahn ذات الاشتقاق الزمني الكسري، وهي نموذج لمجال الطور (phase-field model) يُستخدم لوصف انفصال الطور وتكون الحبيبات في المواد. تُعرف المعادلة على نطاق متعدد الأوجه محدب $\Omega \subset \mathbb{R}^d$ عبر فاصل زمني $J=(0, T]$:

$$\begin{cases} \partial^\alpha_t u - \kappa^2 \Delta u = -F'(u) := f(u) & \text{في } \Omega \times J, \\ u = 0 & \text{على } \partial\Omega \times J, \\ u(x, 0) = u_0(x) & \text{في } \Omega. \end{cases}$$

الخصائص الرئيسية:

• الاشتقاق الكسري: $\partial^\alpha_t$ هو مشتق كابوتو (Caputo) من الرتبة $0 < \alpha < 1$.
• اللاخطية: الحد $F(u) = \frac{1}{4}(1-u^2)^2$ يمثل طاقة الواجهة.
• تحدي الانتظام (Regularity): عادة ما تظهر حلول المسائل ذات الاشتقاق الزمني الكسري نقاط تفرد ضعيفة بالقرب من $t=0$ (أي أن المشتقات تتباعد)، حتى لو كانت البيانات الأولية سلسة. غالبًا ما تفشل مخططات الخطوات الزمنية المنتظمة القياسية في تحقيق معدلات تقارب مثالية في ظل هذه الظروف.
• الهدف: تطوير مخطط عددي كامل التمييز يوفر معدلات تقارب مثالية ويكون متيناً بالنسبة لـ $\alpha$ (أي أن ثوابت الخطأ تظل محدودة عندما $\alpha \to 1^-$)، حتى مع فرض قيود منخفضة نسبيًا على انتظام البيانات الأولية ($u_0 \in H^1_0(\Omega) \cap H^{3+\epsilon}(\Omega)$).

2. المنهجية

يقترح المؤلفون طريقة العناصر المحدودة المختلطة (Mixed FEM) مقترنة بـ مخطط Alikhanov للخطوات الزمنية غير المنتظمة.

أ. التمييز المكاني (Mixed FEM)

• الصياغة: تمت إعادة كتابة المؤثر المكاني من الدرجة الثانية كنظام من الدرجة الأولى عن طريق إدخال متغير التدفق $\sigma = \nabla u$.
• الفضاءات: يتم البحث عن زوج الحل $(u, \sigma)$ في $V \times W$، حيث $V = L^2(\Omega)$ و $W = H(\text{div}; \Omega)$.
• العناصر المحدودة: يستخدم المؤلفون فضاءًا جزئيًا من العناصر المحدودة $V_h \times W_h$ يستوفي خصائص محددة (مثل إسقاط Fortin، و $\nabla \cdot W_h \subseteq V_h$). ومن الأمثلة على ذلك عناصر Raviart-Thomas ($P1_{dc}$ و $RT1$).
• الميزة: تسمح طريقة FEM المختلطة بتقريب متزامن لمجال الطور $u$ وتدفقه $\sigma$ بدقة مثالية، وهو أمر بالغ الأهمية للكميات الفيزيائية مثل تدفق الكتلة.

ب. التمييز الزمني (Alikhanov Scheme)

• الشبكة المتدرجة (Graded Mesh): للتعامل مع نقطة التفرد الأولية، تم استخدام شبكة زمنية متدرجة: $t_n = (n/N)^\gamma T$، حيث $\gamma \ge 1$. تعمل هذه الشبكة على تكثيف الخطوات الزمنية بالقرب من $t=0$.
• صيغة Alikhanov: تم استخدام تقريب من الدرجة الثانية لدقة الاشتقاق الكسري عند مستويات زمنية مزاحة $t_{n-\nu}$ (حيث $\nu = \alpha/2$). وهذا يعد تعميماً لطريقة Crank-Nicolson للمشتقات الكسرية.
• الخطيّة (Linearization): تم خطيّة الحد غير الخطي $u^3$ باستخدام طريقة نيوتن (تحديداً، الخطية حول الخطوة الزمنية السابقة $u^{n-1}$)، مما يتجنب حل أنظمة غير خطية بالكامل في كل خطوة.

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

• تحليل الانتظام: تم وضع نتائج انتظام جديدة للحل $u$ وتدفقه $\sigma$ تحت فرضية $u_0 \in H^1_0(\Omega) \cap H^{3+\epsilon}(\Omega)$. توضح هذه النتائج سلوك التفرد بالقرب من $t=0$.
• متباينة Grönwall المتقطعة الكسرية المعدلة: المساهمة النظرية الرئيسية هي اشتقاق متباينة Grönwall المتقطعة الكسرية المنقحة. تتيح هذه الأداة للمؤلفين استنتاج تقديرات الاستقرار والتقارب تحت قيود أخف على الخطوة الزمنية مقارنة بالأعمال السابقة، مع ضمان عدم تضخم الثوابت عندما $\alpha \to 1$.

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

1. نتائج انتظام جديدة: تثبت الورقة تقديرات انتظام محسنة للحل وتدفقه تحت فرضيات بيانات أولية أضعف ($H^{3+\epsilon}$) من تلك المطلوبة عادة في الأدبيات (التي تتطلب غالباً $H^6$).
2. تقديرات خطأ متينة لـ $\alpha$: يستنتج المؤلفون تقديرات خطأ مثالية في معيار $L^2$ لكل من الحل $u$ والتدفق $\sigma$. والأهم من ذلك، أن حدود الخطأ متينة بالنسبة لـ $\alpha$، مما يعني أن الثوابت في التقديرات تظل محدودة عندما $\alpha \to 1^-$.
3. رتبة تقارب مثالية: تحقق الطريقة معدل تقارب أمثل قدره $O(\log(N)(h^2 + N^{-\min\{\epsilon\gamma\alpha, 2\}}))$ في معيار $L^2$.
   • $h$: حجم الشبكة المكانية.
   • $N$: عدد الخطوات الزمنية.
   • العامل اللوغاريتمي $\log(N)$ هو عقوبة طفيفة نموذجية في مثل هذه التحليلات.
4. شرط خطوة زمنية مرن: باستخدام متباينة Grönwall المعدلة، تتجنب الطريقة القيود الصارمة للخطوة الزمنية ($\Delta t \le C \alpha$) التي توجد غالباً في تحليلات Alikhanov السابقة، مما يجعل المخطط أكثر عملية.

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

• الحد النظري للخطأ:
  $$\max_{1 \le n \le N} \|u_h^n - u(t_n)\| + \max_{1 \le n \le N} t_n^{\alpha/2} \|\sigma_h^n - \sigma(t_n)\| \le C \log(N) (h^2 + N^{-\min\{\epsilon\gamma\alpha, 2\}})$$
  يؤكد هذا تقارب الدرجة الثانية في الفضاء ($h^2$) وتقارب قريب من الدرجة الثانية في الزمن، بشرما تم اختيار معامل التدرج $\gamma$ بشكل مناسب (على سبيل المثال، $\gamma \approx 2/\alpha$).
• المتانة: ثوابت الخطأ $C$ مستقلة عن $\alpha$ عندما تقترب من 1، مما يضمن أن المخطط يعمل بشكل صحيح في الحد الذي يمثل معادلة Allen–Cahn الكلاسيكية.
• التحقق العددي:
  • أُجريت التجارب في بُعدين باستخدام FreeFem++ مع عناصر Raviart-Thomas.
  • المثال 5.1 (بيانات سلسة): أكد معدلات التقارب المثالية ($\approx 2$ في الفضاء والزمن).
  • المثال 5.2 و 5.3 (انتظام محدود): أظهر أن المخطط يحافظ على التقارب الأمثل حتى عندما تكون $u_0 \in H^3$ أو $H^{3.5}$ (أقل من $H^4$ أو $H^6$ المطلوبة في بعض الطرق السابقة).
  • المثال 5.4 (انتظام منخفض): أظهر الاستقرار والتقارب لـ $u_0 \in H^2$، مما يؤكد متانة الطريقة للظروف الأولية غير السلسة.
  • ظلت النتائج صالحة لمختلف الرتب الكسرية ($\alpha = 0.4, 0.6, 0.8, 0.99$).

5. الأهمية

• سد الفجوة: يسد هذا العمل فجوة في الأدبيات المتعلقة بتحليل الخطأ المتين لـ $\alpha$ لطرق العناصر المحدودة المختلطة المطبقة على معادلات Allen–Cahn ذات الاشتقاق الزمني الكسري. ركزت الأعمال السابقة غالبًا على FEM القياسي أو تطلبت فرضيات انتظام عالية جدًا.
• القابلية للتطبيق العملي: إن استخدام شبكة متدرجة ومخطط Alikhanov يجعل الطريقة فعالة للغاية للمسائل التي تعاني من تفردات أولية، وهي شائعة في نماذج الانتشار الكسري.
• الدقة النظرية: يوفر اشتقاق متباينة Grönwall المتقطعة المعدلة أداة قوية لتحليل المسائل الزمنية الكسرية غير الخطية الأخرى، مما قد يؤدي إلى تخفيف قيود الاستقرار في الأبحاث المستقبلية.
• الأهمية الفيزيائية: من خلال توفير تقريبات دقيقة للتدفق ($\sigma = \nabla u$) جنباً إلى جنب مع مجال الطور، فإن الطريقة أفضل ملاءمة للتطبيقات التي تتطلب حسابات دقيقة لديناميكيات الواجهة وديناميكيات نقل الكتلة.

في الختام، تقدم الورقة إطاراً رياضياً صارماً وفعالاً عددياً لحل معادلات Allen–Cahn ذات الاشتقاق الزمني الكسري، مما يوفر تقارباً أمثلاً مع متطلبات انتظام دنيا ومتانة عبر كامل نطاق الرتب الكسرية.